Progestogen-driven B7-H4 contributes to onco-fetal immune tolerance.

Immune tolerance mechanisms are shared in cancer and pregnancy. Through cross-analyzing single-cell RNA-sequencing data from multiple human cancer types and the maternal-fetal interface, we found B7-H4 (VTCN1) is an onco-fetal immune tolerance checkpoint. We showed that genetic deficiency of B7-H4 resulted in immune activation and fetal resorption in allogeneic pregnancy models. Analogously, B7-H4 contributed to MPA/DMBA-induced breast cancer progression, accompanied by CD8+ T cell exhaustion. Female hormone screening revealed that progesterone stimulated B7-H4 expression in placental and breast cancer cells. Mechanistically, progesterone receptor (PR) bound to a newly identified -58 kb enhancer, thereby mediating B7-H4 transcription via the PR-P300-BRD4 axis. PR antagonist or BRD4 degrader potentiated immunotherapy in a murine B7-H4+ breast cancer model. Thus, our work unravels a mechanistic and biological connection of a female sex hormone (progesterone) to onco-fetal immune tolerance via B7-H4 and suggests that the PR-P300-BRD4 axis is targetable for treating B7-H4+ cancer.

Direct observation of a condensate effect on super-enhancer controlled gene bursting.

Enhancers are distal DNA elements believed to loop and contact promoters to control gene expression. Recently, we found diffraction-sized transcriptional condensates at genes controlled by clusters of enhancers (super-enhancers). However, a direct function of endogenous condensates in controlling gene expression remains elusive. Here, we develop live-cell super-resolution and multi-color 3D-imaging approaches to investigate putative roles of endogenous condensates in the regulation of super-enhancer controlled gene Sox2. In contrast to enhancer distance, we find instead that the condensate's positional dynamics are a better predictor of gene expression. A basal gene bursting occurs when the condensate is far (>1 µm), but burst size and frequency are enhanced when the condensate moves in proximity (<1 µm). Perturbations of cohesin and local DNA elements do not prevent basal bursting but affect the condensate and its burst enhancement. We propose a three-way kissing model whereby the condensate interacts transiently with gene locus and regulatory DNA elements to control gene bursting.

The primitive endoderm supports lineage plasticity to enable regulative development.

Mammalian blastocyst formation involves the specification of the trophectoderm followed by the differentiation of the inner cell mass into embryonic epiblast and extra-embryonic primitive endoderm (PrE). During this time, the embryo maintains a window of plasticity and can redirect its cellular fate when challenged experimentally. In this context, we found that the PrE alone was sufficient to regenerate a complete blastocyst and continue post-implantation development. We identify an in vitro population similar to the early PrE in vivo that exhibits the same embryonic and extra-embryonic potency and can form complete stem cell-based embryo models, termed blastoids. Commitment in the PrE is suppressed by JAK/STAT signaling, collaborating with OCT4 and the sustained expression of a subset of pluripotency-related transcription factors that safeguard an enhancer landscape permissive for multi-lineage differentiation. Our observations support the notion that transcription factor persistence underlies plasticity in regulative development and highlight the importance of the PrE in perturbed development.

Super-enhancers include classical enhancers and facilitators to fully activate gene expression.

Super-enhancers are compound regulatory elements that control expression of key cell identity genes. They recruit high levels of tissue-specific transcription factors and co-activators such as the Mediator complex and contact target gene promoters with high frequency. Most super-enhancers contain multiple constituent regulatory elements, but it is unclear whether these elements have distinct roles in activating target gene expression. Here, by rebuilding the endogenous multipartite α-globin super-enhancer, we show that it contains bioinformatically equivalent but functionally distinct element types: classical enhancers and facilitator elements. Facilitators have no intrinsic enhancer activity, yet in their absence, classical enhancers are unable to fully upregulate their target genes. Without facilitators, classical enhancers exhibit reduced Mediator recruitment, enhancer RNA transcription, and enhancer-promoter interactions. Facilitators are interchangeable but display functional hierarchy based on their position within a multipartite enhancer. Facilitators thus play an important role in potentiating the activity of classical enhancers and ensuring robust activation of target genes.

Linking chromatin acylation mark-defined proteome and genome in living cells.

A generalizable strategy with programmable site specificity for in situ profiling of histone modifications on unperturbed chromatin remains highly desirable but challenging. We herein developed a single-site-resolved multi-omics (SiTomics) strategy for systematic mapping of dynamic modifications and subsequent profiling of chromatinized proteome and genome defined by specific chromatin acylations in living cells. By leveraging the genetic code expansion strategy, our SiTomics toolkit revealed distinct crotonylation (e.g., H3K56cr) and ß-hydroxybutyrylation (e.g., H3K56bhb) upon short chain fatty acids stimulation and established linkages for chromatin acylation mark-defined proteome, genome, and functions. This led to the identification of GLYR1 as a distinct interacting protein in modulating H3K56cr's gene body localization as well as the discovery of an elevated super-enhancer repertoire underlying bhb-mediated chromatin modulations. SiTomics offers a platform technology for elucidating the "metabolites-modification-regulation" axis, which is widely applicable for multi-omics profiling and functional dissection of modifications beyond acylations and proteins beyond histones.

Lineage-specific 3D genome organization is assembled at multiple scales by IKAROS.

A generic level of chromatin organization generated by the interplay between cohesin and CTCF suffices to limit promiscuous interactions between regulatory elements, but a lineage-specific chromatin assembly that supersedes these constraints is required to configure the genome to guide gene expression changes that drive faithful lineage progression. Loss-of-function approaches in B cell precursors show that IKAROS assembles interactions across megabase distances in preparation for lymphoid development. Interactions emanating from IKAROS-bound enhancers override CTCF-imposed boundaries to assemble lineage-specific regulatory units built on a backbone of smaller invariant topological domains. Gain of function in epithelial cells confirms IKAROS' ability to reconfigure chromatin architecture at multiple scales. Although the compaction of the Igκ locus required for genome editing represents a function of IKAROS unique to lymphocytes, the more general function to preconfigure the genome to support lineage-specific gene expression and suppress activation of extra-lineage genes provides a paradigm for lineage restriction.

Transcriptional regulation of innate and adaptive lymphocyte lineages.

The lymphocyte family has expanded significantly in recent years to include not only the adaptive lymphocytes (T cells, B cells) and NK cells, but also several additional innate lymphoid cell (ILC) types. ILCs lack clonally distributed antigen receptors characteristic of adaptive lymphocytes and instead respond exclusively to signaling via germline-encoded receptors. ILCs resemble T cells more closely than any other leukocyte lineage at the transcriptome level and express many elements of the core T cell transcriptional program, including Notch, Gata3, Tcf7, and Bcl11b. We present our current understanding of the shared and distinct transcriptional regulatory mechanisms involved in the development of adaptive T lymphocytes and closely related ILCs. We discuss the possibility that a core set of transcriptional regulators common to ILCs and T cells establish enhancers that enable implementation of closely aligned effector pathways. Studies of the transcriptional regulation of lymphopoiesis will support the development of novel therapeutic approaches to correct early lymphoid developmental defects and aberrant lymphocyte function.

Spatially resolved epigenomic profiling of single cells in complex tissues.

The recent development of spatial omics methods has enabled single-cell profiling of the transcriptome and 3D genome organization with high spatial resolution. Expanding the repertoire of spatial omics tools, a spatially resolved single-cell epigenomics method will accelerate understanding of the spatial regulation of cell and tissue functions. Here, we report a method for spatially resolved epigenomic profiling of single cells using in situ tagmentation and transcription followed by multiplexed imaging. We demonstrated the ability to profile histone modifications marking active promoters, putative enhancers, and silent promoters in individual cells, and generated high-resolution spatial atlas of hundreds of active promoters and putative enhancers in embryonic and adult mouse brains. Our results suggested putative promoter-enhancer pairs and enhancer hubs regulating developmentally important genes. We envision this approach will be generally applicable to spatial profiling of epigenetic modifications and DNA-binding proteins, advancing our understanding of how gene expression is spatiotemporally regulated by the epigenome.

Repression and 3D-restructuring resolves regulatory conflicts in evolutionarily rearranged genomes.

Regulatory landscapes drive complex developmental gene expression, but it remains unclear how their integrity is maintained when incorporating novel genes and functions during evolution. Here, we investigated how a placental mammal-specific gene, Zfp42, emerged in an ancient vertebrate topologically associated domain (TAD) without adopting or disrupting the conserved expression of its gene, Fat1. In ESCs, physical TAD partitioning separates Zfp42 and Fat1 with distinct local enhancers that drive their independent expression. This separation is driven by chromatin activity and not CTCF/cohesin. In contrast, in embryonic limbs, inactive Zfp42 shares Fat1's intact TAD without responding to active Fat1 enhancers. However, neither Fat1 enhancer-incompatibility nor nuclear envelope-attachment account for Zfp42's unresponsiveness. Rather, Zfp42's promoter is rendered inert to enhancers by context-dependent DNA methylation. Thus, diverse mechanisms enabled the integration of independent Zfp42 regulation in the Fat1 locus. Critically, such regulatory complexity appears common in evolution as, genome wide, most TADs contain multiple independently expressed genes.

RNA-Mediated Feedback Control of Transcriptional Condensates.

Regulation of biological processes typically incorporates mechanisms that initiate and terminate the process and, where understood, these mechanisms often involve feedback control. Regulation of transcription is a fundamental cellular process where the mechanisms involved in initiation have been studied extensively, but those involved in arresting the process are poorly understood. Modeling of the potential roles of RNA in transcriptional control suggested a non-equilibrium feedback control mechanism where low levels of RNA promote condensates formed by electrostatic interactions whereas relatively high levels promote dissolution of these condensates. Evidence from in vitro and in vivo experiments support a model where RNAs produced during early steps in transcription initiation stimulate condensate formation, whereas the burst of RNAs produced during elongation stimulate condensate dissolution. We propose that transcriptional regulation incorporates a feedback mechanism whereby transcribed RNAs initially stimulate but then ultimately arrest the process.

Comprehensive In Vivo Interrogation Reveals Phenotypic Impact of Human Enhancer Variants.

Establishing causal links between non-coding variants and human phenotypes is an increasing challenge. Here, we introduce a high-throughput mouse reporter assay for assessing the pathogenic potential of human enhancer variants in vivo and examine nearly a thousand variants in an enhancer repeatedly linked to polydactyly. We show that 71% of all rare non-coding variants previously proposed as causal lead to reporter gene expression in a pattern consistent with their pathogenic role. Variants observed to alter enhancer activity were further confirmed to cause polydactyly in knockin mice. We also used combinatorial and single-nucleotide mutagenesis to evaluate the in vivo impact of mutations affecting all positions of the enhancer and identified additional functional substitutions, including potentially pathogenic variants hitherto not observed in humans. Our results uncover the functional consequences of hundreds of mutations in a phenotype-associated enhancer and establish a widely applicable strategy for systematic in vivo evaluation of human enhancer variants.

The cis-Regulatory Atlas of the Mouse Immune System.

A complete chart of cis-regulatory elements and their dynamic activity is necessary to understand the transcriptional basis of differentiation and function of an organ system. We generated matched epigenome and transcriptome measurements in 86 primary cell types that span the mouse immune system and its differentiation cascades. This breadth of data enable variance components analysis that suggests that genes fall into two distinct classes, controlled by either enhancer- or promoter-driven logic, and multiple regression that connects genes to the enhancers that regulate them. Relating transcription factor (TF) expression to the genome-wide accessibility of their binding motifs classifies them as predominantly openers or closers of local chromatin accessibility, pinpointing specific cis-regulatory elements where binding of given TFs is likely functionally relevant, validated by chromatin immunoprecipitation sequencing (ChIP-seq). Overall, this cis-regulatory atlas provides a trove of information on transcriptional regulation through immune differentiation and a foundational scaffold to define key regulatory events throughout the immunological genome.

Functional Enhancers Shape Extrachromosomal Oncogene Amplifications.

Non-coding regions amplified beyond oncogene borders have largely been ignored. Using a computational approach, we find signatures of significant co-amplification of non-coding DNA beyond the boundaries of amplified oncogenes across five cancer types. In glioblastoma, EGFR is preferentially co-amplified with its two endogenous enhancer elements active in the cell type of origin. These regulatory elements, their contacts, and their contribution to cell fitness are preserved on high-level circular extrachromosomal DNA amplifications. Interrogating the locus with a CRISPR interference screening approach reveals a diversity of additional elements that impact cell fitness. The pattern of fitness dependencies mirrors the rearrangement of regulatory elements and accompanying rewiring of the chromatin topology on the extrachromosomal amplicon. Our studies indicate that oncogene amplifications are shaped by regulatory dependencies in the non-coding genome.

A Genome-wide Framework for Mapping Gene Regulation via Cellular Genetic Screens.

Over one million candidate regulatory elements have been identified across the human genome, but nearly all are unvalidated and their target genes uncertain. Approaches based on human genetics are limited in scope to common variants and in resolution by linkage disequilibrium. We present a multiplex, expression quantitative trait locus (eQTL)-inspired framework for mapping enhancer-gene pairs by introducing random combinations of CRISPR/Cas9-mediated perturbations to each of many cells, followed by single-cell RNA sequencing (RNA-seq). Across two experiments, we used dCas9-KRAB to perturb 5,920 candidate enhancers with no strong a priori hypothesis as to their target gene(s), measuring effects by profiling 254,974 single-cell transcriptomes. We identified 664 (470 high-confidence) cis enhancer-gene pairs, which were enriched for specific transcription factors, non-housekeeping status, and genomic and 3D conformational proximity to their target genes. This framework will facilitate the large-scale mapping of enhancer-gene regulatory interactions, a critical yet largely uncharted component of the cis-regulatory landscape of the human genome.

Identifying cis Elements for Spatiotemporal Control of Mammalian DNA Replication.

The temporal order of DNA replication (replication timing [RT]) is highly coupled with genome architecture, but cis-elements regulating either remain elusive. We created a series of CRISPR-mediated deletions and inversions of a pluripotency-associated topologically associating domain (TAD) in mouse ESCs. CTCF-associated domain boundaries were dispensable for RT. CTCF protein depletion weakened most TAD boundaries but had no effect on RT or A/B compartmentalization genome-wide. By contrast, deletion of three intra-TAD CTCF-independent 3D contact sites caused a domain-wide early-to-late RT shift, an A-to-B compartment switch, weakening of TAD architecture, and loss of transcription. The dispensability of TAD boundaries and the necessity of these "early replication control elements" (ERCEs) was validated by deletions and inversions at additional domains. Our results demonstrate that discrete cis-regulatory elements orchestrate domain-wide RT, A/B compartmentalization, TAD architecture, and transcription, revealing fundamental principles linking genome structure and function.

Regulation of T helper cell differentiation by the interplay between histone modification and chromatin interaction.

The development and function of the immune system are controlled by temporospatial gene expression programs, which are regulated by cis-regulatory elements, chromatin structure, and trans-acting factors. In this study, we cataloged the dynamic histone modifications and chromatin interactions at regulatory regions during T helper (Th) cell differentiation. Our data revealed that the H3K4me1 landscape established by MLL4 in naive CD4+ T cells is critical for restructuring the regulatory interaction network and orchestrating gene expression during the early phase of Th differentiation. GATA3 plays a crucial role in further configuring H3K4me1 modification and the chromatin interaction network during Th2 differentiation. Furthermore, we demonstrated that HSS3-anchored chromatin loops function to restrict the activity of the Th2 locus control region (LCR), thus coordinating the expression of Th2 cytokines. Our results provide insights into the mechanisms of how the interplay between histone modifications, chromatin looping, and trans-acting factors contributes to the differentiation of Th cells.

Collaboration between distinct SWI/SNF chromatin remodeling complexes directs enhancer selection and activation of macrophage inflammatory genes.

Macrophages elicit immune responses to pathogens through induction of inflammatory genes. Here, we examined the role of three variants of the SWI/SNF nucleosome remodeling complex-cBAF, ncBAF, and PBAF-in the macrophage response to bacterial endotoxin (lipid A). All three SWI/SNF variants were prebound in macrophages and retargeted to genomic sites undergoing changes in chromatin accessibility following stimulation. Cooperative binding of all three variants associated with de novo chromatin opening and latent enhancer activation. Isolated binding of ncBAF and PBAF, in contrast, associated with activation and repression of active enhancers, respectively. Chemical and genetic perturbations of variant-specific subunits revealed pathway-specific regulation in the activation of lipid A response genes, corresponding to requirement for cBAF and ncBAF in inflammatory and interferon-stimulated gene (ISG) activation, respectively, consistent with differential engagement of SWI/SNF variants by signal-responsive transcription factors. Thus, functional diversity among SWI/SNF variants enables increased regulatory control of innate immune transcriptional programs, with potential for specific therapeutic targeting.

CBFß-SMMHC Inhibition Triggers Apoptosis by Disrupting MYC Chromatin Dynamics in Acute Myeloid Leukemia.

The fusion oncoprotein CBFß-SMMHC, expressed in leukemia cases with chromosome 16 inversion, drives leukemia development and maintenance by altering the activity of the transcription factor RUNX1. Here, we demonstrate that CBFß-SMMHC maintains cell viability by neutralizing RUNX1-mediated repression of MYC expression. Upon pharmacologic inhibition of the CBFß-SMMHC/RUNX1 interaction, RUNX1 shows increased binding at three MYC distal enhancers, where it represses MYC expression by mediating the replacement of the SWI/SNF complex component BRG1 with the polycomb-repressive complex component RING1B, leading to apoptosis. Combining the CBFß-SMMHC inhibitor with the BET inhibitor JQ1 eliminates inv(16) leukemia in human cells and a mouse model. Enhancer-interaction analysis indicated that the three enhancers are physically connected with the MYC promoter, and genome-editing analysis demonstrated that they are functionally implicated in deregulation of MYC expression. This study reveals a mechanism whereby CBFß-SMMHC drives leukemia maintenance and suggests that inhibitors targeting chromatin activity may prove effective in inv(16) leukemia therapy.

Analysis of Genetically Diverse Macrophages Reveals Local and Domain-wide Mechanisms that Control Transcription Factor Binding and Function.

Non-coding genetic variation is a major driver of phenotypic diversity and allows the investigation of mechanisms that control gene expression. Here, we systematically investigated the effects of >50 million variations from five strains of mice on mRNA, nascent transcription, transcription start sites, and transcription factor binding in resting and activated macrophages. We observed substantial differences associated with distinct molecular pathways. Evaluating genetic variation provided evidence for roles of ∼100 TFs in shaping lineage-determining factor binding. Unexpectedly, a substantial fraction of strain-specific factor binding could not be explained by local mutations. Integration of genomic features with chromatin interaction data provided evidence for hundreds of connected cis-regulatory domains associated with differences in transcription factor binding and gene expression. This system and the >250 datasets establish a substantial new resource for investigation of how genetic variation affects cellular phenotypes.

Time-Resolved Analysis Reveals Rapid Dynamics and Broad Scope of the CBP/p300 Acetylome.

The acetyltransferases CBP and p300 are multifunctional transcriptional co-activators. Here, we combined quantitative proteomics with CBP/p300-specific catalytic inhibitors, bromodomain inhibitor, and gene knockout to reveal a comprehensive map of regulated acetylation sites and their dynamic turnover rates. CBP/p300 acetylates thousands of sites, including signature histone sites and a multitude of sites on signaling effectors and enhancer-associated transcriptional regulators. Time-resolved acetylome analyses identified a subset of CBP/p300-regulated sites with very rapid (<30 min) acetylation turnover, revealing a dynamic balance between acetylation and deacetylation. Quantification of acetylation, mRNA, and protein abundance after CBP/p300 inhibition reveals a kinetically competent network of gene expression that strictly depends on CBP/p300-catalyzed rapid acetylation. Collectively, our in-depth acetylome analyses reveal systems attributes of CBP/p300 targets, and the resource dataset provides a framework for investigating CBP/p300 functions and for understanding the impact of small-molecule inhibitors targeting its catalytic and bromodomain activities.