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.

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.

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.

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.

Genomic context sensitizes regulatory elements to genetic disruption.

Genomic context critically modulates regulatory function but is difficult to manipulate systematically. The murine insulin-like growth factor 2 (Igf2)/H19 locus is a paradigmatic model of enhancer selectivity, whereby CTCF occupancy at an imprinting control region directs downstream enhancers to activate either H19 or Igf2. We used synthetic regulatory genomics to repeatedly replace the native locus with 157-kb payloads, and we systematically dissected its architecture. Enhancer deletion and ectopic delivery revealed previously uncharacterized long-range regulatory dependencies at the native locus. Exchanging the H19 enhancer cluster with the Sox2 locus control region (LCR) showed that the H19 enhancers relied on their native surroundings while the Sox2 LCR functioned autonomously. Analysis of regulatory DNA actuation across cell types revealed that these enhancer clusters typify broader classes of context sensitivity genome wide. These results show that unexpected dependencies influence even well-studied loci, and our approach permits large-scale manipulation of complete loci to investigate the relationship between regulatory architecture and function.

Transcriptional reprogramming at the intersection of the heat shock response and proteostasis.

Cellular homeostasis is constantly challenged by a myriad of extrinsic and intrinsic stressors. To mitigate the stress-induced damage, cells activate transient survival programs. The heat shock response (HSR) is an evolutionarily well-conserved survival program that is activated in response to proteotoxic stress. The HSR encompasses a dual regulation of transcription, characterized by rapid activation of genes encoding molecular chaperones and concomitant global attenuation of non-chaperone genes. Recent genome-wide approaches have delineated the molecular depth of stress-induced transcriptional reprogramming. The dramatic rewiring of gene and enhancer networks is driven by key transcription factors, including heat shock factors (HSFs), that together with chromatin-modifying enzymes remodel the 3D chromatin architecture, determining the selection of either gene activation or repression. Here, we highlight the current advancements of molecular mechanisms driving transcriptional reprogramming during acute heat stress. We also discuss the emerging implications of HSF-mediated stress signaling in the context of physiological and pathological conditions.

Dissection of a CTCF topological boundary uncovers principles of enhancer-oncogene regulation.

Enhancer-gene communication is dependent on topologically associating domains (TADs) and boundaries enforced by the CCCTC-binding factor (CTCF) insulator, but the underlying structures and mechanisms remain controversial. Here, we investigate a boundary that typically insulates fibroblast growth factor (FGF) oncogenes but is disrupted by DNA hypermethylation in gastrointestinal stromal tumors (GISTs). The boundary contains an array of CTCF sites that enforce adjacent TADs, one containing FGF genes and the other containing ANO1 and its putative enhancers, which are specifically active in GIST and its likely cell of origin. We show that coordinate disruption of four CTCF motifs in the boundary fuses the adjacent TADs, allows the ANO1 enhancer to contact FGF3, and causes its robust induction. High-resolution micro-C maps reveal specific contact between transcription initiation sites in the ANO1 enhancer and FGF3 promoter that quantitatively scales with FGF3 induction such that modest changes in contact frequency result in strong changes in expression, consistent with a causal relationship.

Enhancer switching in cell lineage priming is linked to eRNA, Brg1's AT-hook, and SWI/SNF recruitment.

RNA transcribed from enhancers, i.e., eRNA, has been suggested to directly activate transcription by recruiting transcription factors and co-activators. Although there have been specific examples of eRNA functioning in this way, it is not clear how general this may be. We find that the AT-hook of SWI/SNF preferentially binds RNA and, as part of the esBAF complex, associates with eRNA transcribed from intronic and intergenic regions. Our data suggest that SWI/SNF is globally recruited in cis by eRNA to cell-type-specific enhancers, representative of two distinct stages that mimic early mammalian development, and not at enhancers that are shared between the two stages. In this manner, SWI/SNF facilitates recruitment and/or activation of MLL3/4, p300/CBP, and Mediator to stage-specific enhancers and super-enhancers that regulate the transcription of metabolic and cell lineage priming-related genes. These findings highlight a connection between ATP-dependent chromatin remodeling and eRNA in cell identity and typical- and super-enhancer activation.

Structural insights into the cooperative nucleosome recognition and chromatin opening by FOXA1 and GATA4.

Mouse FOXA1 and GATA4 are prototypes of pioneer factors, initiating liver cell development by binding to the N1 nucleosome in the enhancer of the ALB1 gene. Using cryoelectron microscopy (cryo-EM), we determined the structures of the free N1 nucleosome and its complexes with FOXA1 and GATA4, both individually and in combination. We found that the DNA-binding domains of FOXA1 and GATA4 mainly recognize the linker DNA and an internal site in the nucleosome, respectively, whereas their intrinsically disordered regions interact with the acidic patch on histone H2A-H2B. FOXA1 efficiently enhances GATA4 binding by repositioning the N1 nucleosome. In vivo DNA editing and bioinformatics analyses suggest that the co-binding mode of FOXA1 and GATA4 plays important roles in regulating genes involved in liver cell functions. Our results reveal the mechanism whereby FOXA1 and GATA4 cooperatively bind to the nucleosome through nucleosome repositioning, opening chromatin by bending linker DNA and obstructing nucleosome packing.

H3K4me1 facilitates promoter-enhancer interactions and gene activation during embryonic stem cell differentiation.

Histone H3 lysine 4 mono-methylation (H3K4me1) marks poised or active enhancers. KMT2C (MLL3) and KMT2D (MLL4) catalyze H3K4me1, but their histone methyltransferase activities are largely dispensable for transcription during early embryogenesis in mammals. To better understand the role of H3K4me1 in enhancer function, we analyze dynamic enhancer-promoter (E-P) interactions and gene expression during neural differentiation of the mouse embryonic stem cells. We found that KMT2C/D catalytic activities were only required for H3K4me1 and E-P contacts at a subset of candidate enhancers, induced upon neural differentiation. By contrast, a majority of enhancers retained H3K4me1 in KMT2C/D catalytic mutant cells. Surprisingly, H3K4me1 signals at these KMT2C/D-independent sites were reduced after acute depletion of KMT2B, resulting in aggravated transcriptional defects. Our observations therefore implicate KMT2B in the catalysis of H3K4me1 at enhancers and provide additional support for an active role of H3K4me1 in enhancer-promoter interactions and transcription in mammalian cells.

High-throughput PRIME-editing screens identify functional DNA variants in the human genome.

Despite tremendous progress in detecting DNA variants associated with human disease, interpreting their functional impact in a high-throughput and single-base resolution manner remains challenging. Here, we develop a pooled prime-editing screen method, PRIME, that can be applied to characterize thousands of coding and non-coding variants in a single experiment with high reproducibility. To showcase its applications, we first identified essential nucleotides for a 716 bp MYC enhancer via PRIME-mediated single-base resolution analysis. Next, we applied PRIME to functionally characterize 1,304 genome-wide association study (GWAS)-identified non-coding variants associated with breast cancer and 3,699 variants from ClinVar. We discovered that 103 non-coding variants and 156 variants of uncertain significance are functional via affecting cell fitness. Collectively, we demonstrate that PRIME is capable of characterizing genetic variants at single-base resolution and scale, advancing accurate genome annotation for disease risk prediction, diagnosis, and therapeutic target identification.

Time will tell: comparing timescales to gain insight into transcriptional bursting.

Recent imaging studies have captured the dynamics of regulatory events of transcription inside living cells. These events include transcription factor (TF) DNA binding, chromatin remodeling and modification, enhancer-promoter (E-P) proximity, cluster formation, and preinitiation complex (PIC) assembly. Together, these molecular events culminate in stochastic bursts of RNA synthesis, but their kinetic relationship remains largely unclear. In this review, we compare the timescales of upstream regulatory steps (input) with the kinetics of transcriptional bursting (output) to generate mechanistic models of transcription dynamics in single cells. We highlight open questions and potential technical advances to guide future endeavors toward a quantitative and kinetic understanding of transcription regulation.

Uncovering dark matter in cancer by identifying epigenetic drivers.

The complex relationship between chromatin accessibility, transcriptional regulation, and cancer transitions presents a daunting puzzle. Terekhanova et al. created a pan-cancer epigenetic and transcriptomic atlas at single-cell resolution, yielding important insights into the underlying chromatin architecture of cancer transitions and novel discoveries with the potential to advance precision medicine.

Neuronal MAPT expression is mediated by long-range interactions with cis-regulatory elements.

Tauopathies are a group of neurodegenerative diseases defined by abnormal aggregates of tau, a microtubule-associated protein encoded by MAPT. MAPT expression is near absent in neural progenitor cells (NPCs) and increases during differentiation. This temporally dynamic expression pattern suggests that MAPT expression could be controlled by transcription factors and cis-regulatory elements specific to differentiated cell types. Given the relevance of MAPT expression to neurodegeneration pathogenesis, identification of such elements is relevant to understanding disease risk and pathogenesis. Here, we performed chromatin conformation assays (HiC & Capture-C), single-nucleus multiomics (RNA-seq+ATAC-seq), bulk ATAC-seq, and ChIP-seq for H3K27ac and CTCF in NPCs and differentiated neurons to nominate candidate cis-regulatory elements (cCREs). We assayed these cCREs using luciferase assays and CRISPR interference (CRISPRi) experiments to measure their effects on MAPT expression. Finally, we integrated cCRE annotations into an analysis of genetic variation in neurodegeneration-affected individuals and control subjects. We identified both proximal and distal regulatory elements for MAPT and confirmed the regulatory function for several regions, including three regions centromeric to MAPT beyond the H1/H2 haplotype inversion breakpoint. We also found that rare and predicted damaging genetic variation in nominated CREs was nominally depleted in dementia-affected individuals relative to control subjects, consistent with the hypothesis that variants that disrupt MAPT enhancer activity, and thereby reduced MAPT expression, may be protective against neurodegenerative disease. Overall, this study provides compelling evidence for pursuing detailed knowledge of CREs for genes of interest to permit better understanding of disease risk.

Transcriptome- and proteome-wide association studies identify genes associated with renal cell carcinoma.

We performed a series of integrative analyses including transcriptome-wide association studies (TWASs) and proteome-wide association studies (PWASs) of renal cell carcinoma (RCC) to nominate and prioritize molecular targets for laboratory investigation. On the basis of a genome-wide association study (GWAS) of 29,020 affected individuals and 835,670 control individuals and prediction models trained in transcriptomic reference models, our TWAS across four kidney transcriptomes (GTEx kidney cortex, kidney tubules, TCGA-KIRC [The Cancer Genome Atlas kidney renal clear-cell carcinoma], and TCGA-KIRP [TCGA kidney renal papillary cell carcinoma]) identified 38 gene associations (false-discovery rate <5%) in at least two of four transcriptomic panels and identified 12 genes that were independent of GWAS susceptibility regions. Analyses combining TWAS associations across 48 tissues from GTEx identified associations that were replicable in tumor transcriptomes for 23 additional genes. Analyses by the two major histologic types (clear-cell RCC and papillary RCC) revealed subtype-specific associations, although at least three gene associations were common to both subtypes. PWAS identified 13 associated proteins, all mapping to GWAS-significant loci. TWAS-identified genes were enriched for active enhancer or promoter regions in RCC tumors and hypoxia-inducible factor binding sites in relevant cell lines. Using gene expression correlation, common cancers (breast and prostate) and RCC risk factors (e.g., hypertension and BMI) display genetic contributions shared with RCC. Our work identifies potential molecular targets for RCC susceptibility for downstream functional investigation.

Systematic identification of genotype-dependent enhancer variants in eosinophilic esophagitis.

Eosinophilic esophagitis (EoE) is a rare atopic disorder associated with esophageal dysfunction, including difficulty swallowing, food impaction, and inflammation, that develops in a small subset of people with food allergies. Genome-wide association studies (GWASs) have identified 9 independent EoE risk loci reaching genome-wide significance (p < 5 × 10-8) and 27 additional loci of suggestive significance (5 × 10-8 < p < 1 × 10-5). In the current study, we perform linkage disequilibrium (LD) expansion of these loci to nominate a set of 531 variants that are potentially causal. To systematically interrogate the gene regulatory activity of these variants, we designed a massively parallel reporter assay (MPRA) containing the alleles of each variant within their genomic sequence context cloned into a GFP reporter library. Analysis of reporter gene expression in TE-7, HaCaT, and Jurkat cells revealed cell-type-specific gene regulation. We identify 32 allelic enhancer variants, representing 6 genome-wide significant EoE loci and 7 suggestive EoE loci, that regulate reporter gene expression in a genotype-dependent manner in at least one cellular context. By annotating these variants with expression quantitative trait loci (eQTL) and chromatin looping data in related tissues and cell types, we identify putative target genes affected by genetic variation in individuals with EoE. Transcription factor enrichment analyses reveal possible roles for cell-type-specific regulators, including GATA3. Our approach reduces the large set of EoE-associated variants to a set of 32 with allelic regulatory activity, providing functional insights into the effects of genetic variation in this disease.

Evolutionary conservation of VSX2 super-enhancer modules in retinal development.

Super-enhancers (SEs) are expansive regions of genomic DNA that regulate the expression of genes involved in cell identity and cell fate. We recently identified developmental stage- and cell type-specific modules within the murine Vsx2 SE. Here, we show that the human VSX2 SE modules have similar developmental stage- and cell type-specific activity in reporter gene assays. By inserting the human sequence of one VSX2 SE module into a mouse with microphthalmia, eye size was rescued. To understand the function of these SE modules during human retinal development, we deleted individual modules in human embryonic stem cells and generated retinal organoids. Deleting one module results in small organoids, recapitulating the small-eyed phenotype of mice with microphthalmia, while deletion of the other module led to disruptions in bipolar neuron development. This prototypical SE serves as a model for understanding developmental stage- and cell type-specific effects of neurogenic transcription factors with complex expression patterns. Moreover, by elucidating the gene regulatory mechanisms, we can begin to examine how dysregulation of these mechanisms contributes to phenotypic diversity and disease.

ETS transcription factors regulate precise matrix metalloproteinase expression and follicle rupture in Drosophila.

Drosophila matrix metalloproteinase 2 (MMP2) is specifically expressed in posterior follicle cells of stage-14 egg chambers (mature follicles) and is crucial for the breakdown of the follicular wall during ovulation, a process that is highly conserved from flies to mammals. The factors that regulate spatiotemporal expression of MMP2 in follicle cells remain unknown. Here, we demonstrate crucial roles for the ETS-family transcriptional activator Pointed (Pnt) and its endogenous repressor Yan in the regulation of MMP2 expression. We found that Pnt is expressed in posterior follicle cells and overlaps with MMP2 expression in mature follicles. Genetic analysis demonstrated that pnt is both required and sufficient for MMP2 expression in follicle cells. In addition, Yan was temporally upregulated in stage-13 follicle cells to fine-tune Pnt activity and MMP2 expression. Furthermore, we identified a 1.1 kb core enhancer that is responsible for the spatiotemporal expression of MMP2 and contains multiple pnt/yan binding motifs. Mutation of pnt/yan binding sites significantly impaired the Mmp2 enhancer activity. Our data reveal a mechanism of transcriptional regulation of Mmp2 expression in Drosophila ovulation, which could be conserved in other biological systems.