Transcriptional coupling of distant regulatory genes in living embryos

Bridging spatial and temporal scales of developmental gene regulation

The development of multicellular organisms relies on the precise coordination of molecular events across multiple spatial and temporal scales. Understanding how information flows from molecular interactions to cellular processes and tissue organization during development is crucial for explaining the remarkable reproducibility of complex organisms. This review explores how chromatin-encoded information is transduced from localized transcriptional events to global gene expression patterns, highlighting the challenge of bridging these scales. We discuss recent experimental findings and theoretical frameworks, emphasizing polymer physics as a tool for describing the relationship between chromatin structure and dynamics across scales. By integrating these perspectives, we aim to clarify how gene regulation is coordinated across levels of biological organization and suggest strategies for future experimental approaches.

Massively parallel reporter assays and mouse transgenic assays provide correlated and complementary information about neuronal enhancer activity

High-throughput massively parallel reporter assays (MPRAs) and phenotype-rich in vivo transgenic mouse assays are two potentially complementary ways to study the impact of noncoding variants associated with psychiatric diseases. Here, we investigate the utility of combining these assays. Specifically, we carry out an MPRA in induced human neurons on over 50,000 sequences derived from fetal neuronal ATAC-seq datasets and enhancers validated in mouse assays. We also test the impact of over 20,000 variants, including synthetic mutations and 167 common variants associated with psychiatric disorders. We find a strong and specific correlation between MPRA and mouse neuronal enhancer activity. Four out of five tested variants with significant MPRA effects affected neuronal enhancer activity in mouse embryos. Mouse assays also reveal pleiotropic variant effects that could not be observed in MPRA. Our work provides a catalog of functional neuronal enhancers and variant effects and highlights the effectiveness of combining MPRAs and mouse transgenic assays.

Solid phase transitions as a solution to the genome folding paradox

Ultra-long-range genomic contacts, which are key components of neuronal genome architecture1-3, constitute a biochemical enigma. This is because regulatory DNA elements make selective and stable contacts with DNA sequences located hundreds of kilobases away, instead of interacting with proximal sequences occupied by the exact same transcription factors1,4. This is exemplified in olfactory sensory neurons (OSNs), in which only a fraction of LHX2-, EBF1- and LDB1-bound sites interact with each other, converging into highly selective multi-chromosomal enhancer hubs5. To obtain biochemical insight into this process, here we assembled olfactory receptor (OR) enhancer hubs in vitro with recombinant proteins and enhancer DNA. Cell-free reconstitution of enhancer hubs revealed that OR enhancers form nucleoprotein condensates with unusual, solid-like characteristics. Assembly of these solid condensates is orchestrated by specific DNA motifs enriched in OR enhancers, which are likely to confer distinct homotypic properties on their resident LHX2-EBF1-LDB1 complexes. Single-molecule tracking and pulse-chase experiments in vivo confirmed that LHX2 and EBF1 assemble OR-transcription-competent condensates with solid properties in OSN nuclei, under physiological concentrations of protein. Thus, homophilic nucleoprotein interactions that are influenced by DNA sequence generate new types of biomolecular condensate, which might provide a generalizable explanation for the stability and specificity of long-range genomic contacts across cell types.

Coarse-grained chromatin dynamics by tracking multiple similarly labeled gene loci

The "holy grail" of chromatin research would be to follow the chromatin configuration in individual live cells over time. One way to achieve this goal would be to track the positions of multiple loci arranged along the chromatin polymer with fluorescent labels. Using distinguishable labels would define each locus uniquely in a microscopic image but would restrict the number of loci that could be observed simultaneously due to experimental limits to the number of distinguishable labels. Using the same label for all loci circumvents this limitation but requires a (currently lacking) framework for how to establish each observed locus identity, i.e., to which genomic position it corresponds. Here, we analyze theoretically, using simulations of Rouse model polymers, how single-particle tracking of multiple identically labeled loci enables the determination of loci identity. We show that the probability of correctly assigning observed loci to genomic positions converges exponentially to unity as the number of observed loci configurations increases. The convergence rate depends only weakly on the number of labeled loci, so that even large numbers of loci can be identified with high fidelity by tracking them across about eight independent chromatin configurations. In the case of two distinct labels that alternate along the chromatin polymer, we find that the probability of the correct assignment converges faster than for same-labeled loci, requiring observation of fewer independent chromatin configurations to establish loci identities. Finally, for a modified Rouse model polymer, which realizes a population of dynamic loops, we find that the success probability also converges to unity exponentially as the number of observed loci configurations increases, albeit slightly more slowly than for a classical Rouse model polymer. Altogether, these results establish particle tracking of multiple identically or alternately labeled loci over time as a feasible way to infer temporal dynamics of the coarse-grained configuration of the chromatin polymer in individual living cells.

Chromatin loops are an ancestral hallmark of the animal regulatory genome

In bilaterian animals, gene regulation is shaped by a combination of linear and spatial regulatory information. Regulatory elements along the genome are integrated into gene regulatory landscapes through chromatin compartmentalization1,2, insulation of neighbouring genomic regions3,4 and chromatin looping that brings together distal cis-regulatory sequences5. However, the evolution of these regulatory features is unknown because the three-dimensional genome architecture of most animal lineages remains unexplored6,7. To trace the evolutionary origins of animal genome regulation, here we characterized the physical organization of the genome in non-bilaterian animals (sponges, ctenophores, placozoans and cnidarians)8,9 and their closest unicellular relatives (ichthyosporeans, filastereans and choanoflagellates)10 by combining high-resolution chromosome conformation capture11,12 with epigenomic marks and gene expression data. Our comparative analysis showed that chromatin looping is a conserved feature of genome architecture in ctenophores, placozoans and cnidarians. These sequence-determined distal contacts involve both promoter-enhancer and promoter-promoter interactions. By contrast, chromatin loops are absent in the unicellular relatives of animals. Our findings indicate that spatial genome regulation emerged early in animal evolution. This evolutionary innovation introduced regulatory complexity, ultimately facilitating the diversification of animal developmental programmes and cell type repertoires.

Three-dimensional regulatory hubs support oncogenic programs in glioblastoma

Dysregulation of enhancer-promoter communication in the three-dimensional (3D) nucleus is increasingly recognized as a potential driver of oncogenic programs. Here, we profiled the 3D enhancer-promoter networks of patient-derived glioblastoma stem cells to identify central regulatory nodes. We focused on hyperconnected 3D hubs and demonstrated that hub-interacting genes exhibit high and coordinated expression at the single-cell level and are associated with oncogenic programs that distinguish glioblastoma from low-grade glioma. Epigenetic silencing of a recurrent hub-with an uncharacterized role in glioblastoma-was sufficient to cause downregulation of hub-connected genes, shifts in transcriptional states, and reduced clonogenicity. Integration of datasets across 16 cancers identified "universal" and cancer-type-specific 3D hubs that enrich for oncogenic programs and factors associated with worse prognosis. Genetic alterations could explain only a small fraction of hub hyperconnectivity and increased activity. Overall, our study provides strong support for the potential central role of 3D regulatory hubs in controlling oncogenic programs and properties.

Principles of long-range gene regulation

Transcription from gene promoters occurs in specific spatiotemporal patterns in multicellular organisms, controlled by genomic regulatory elements. The communication between a regulatory element and a promoter requires a certain degree of physical proximity between them; hence, most gene regulation occurs locally in the genome. However, recent discoveries have revealed long-range gene regulation strategies that enhance interactions between regulatory elements and promoters by overcoming the distances between them in the linear genome. These new findings challenge the traditional view of how gene expression patterns are controlled. This review examines long-range gene regulation strategies recently reported in Drosophila and mammals, offering insights into their mechanisms and evolution.

Regulation of 3D genome organization during T cell activation

Within the three-dimensional (3D) nuclear space, the genome organizes into a series of orderly structures that impose important influences on gene regulation. T lymphocytes, crucial players in adaptive immune responses, undergo intricate transcriptional remodeling upon activation, leading to differentiation into specific effector and memory T cell subsets. Recent evidence suggests that T cell activation is accompanied by dynamic changes in genome architecture at multiple levels, providing a unique biological context to explore the functional relevance and molecular mechanisms of 3D genome organization. Here, we summarize recent advances that link the reorganization of genome architecture to the remodeling of transcriptional programs and conversion of cell fates during T cell activation and differentiation. We further discuss how various chromatin architecture regulators, including CCCTC-binding factor and several transcription factors, collectively modulate the genome architecture during this process.

Loop Catalog: a comprehensive HiChIP database of human and mouse samples

HiChIP enables cost-effective and high-resolution profiling of chromatin loops. To leverage the increasing number of HiChIP datasets, we developed Loop Catalog (https://loopcatalog.lji.org), a web-based database featuring loop calls from 1000+ distinct human and mouse HiChIP samples from 152 studies plus 44 high-resolution Hi-C samples. We demonstrate its utility for interpreting GWAS and eQTL variants through SNP-to-gene linking, identifying enriched sequence motifs and motif pairs, and generating regulatory networks and 2D representations of chromatin structure. Our catalog spans over 4.19M unique loops, and with embedded analysis modules, constitutes an important resource for the field.

Two unrelated distal genes activated by a shared enhancer benefit from localizing inside the same small topological domain

Enhancers are tissue-specific regulatory DNA elements that can activate transcription of genes over distance. Their target genes most often are located in the same contact domain-chromosomal entities formed by cohesin DNA loop extrusion and typically flanked by CTCF-bound boundaries. Enhancers shared by multiple unrelated genes are underexplored but may be more common than anticipated. Here, we analyzed the interplay between an enhancer and two distal functionally unrelated genes residing at opposite domain boundaries. The enhancer strongly activated their expression and supported their frequent interactions. Cohesin structured the domain and supported their transcription, but the genes did not rely on each other's transcription or show gene competition. Deleting either domain boundary not only extended the contact domain but led to reduced contacts within the original domain and reduction in the expression of both genes. Conversely, by isolating either gene with the enhancer in shorter domains, through insertion of new CTCF boundaries, intradomain contact frequencies increased, and the gene isolated with the enhancer was upregulated. Collectively, this shows that an enhancer can independently activate unrelated distal genes and that long-range gene regulation benefits from operating in small contact domains.

Coarse-grained chromatin dynamics by tracking multiple similarly labeled gene loci

The "holy grail" of chromatin research would be to follow the chromatin configuration in individual live cells over time. One way to achieve this goal would be to track the positions of multiple loci arranged along the chromatin polymer with fluorescent labels. Use of distinguishable labels would define each locus uniquely in a microscopic image but would restrict the number of loci that could be observed simultaneously, because of experimental limits to the number of distinguishable labels. Use of the same label for all loci circumvents this limitation but requires a (currently lacking) framework for how to establish each observed locus identity, i.e. to which genomic position it corresponds. Here we analyze theoretically, using simulations of Rouse-model polymers, how single-particle-tracking of multiple identically-labeled loci enables determination of loci identity. We show that the probability of correctly assigning observed loci to genomic positions converges exponentially to unity as the number of observed loci configurations increases. The convergence rate depends only weakly on the number of labeled loci, so that even large numbers of loci can be identified with high fidelity by tracking them across about 8 independent chromatin configurations. In the case of two distinct labels that alternate along the chromatin polymer, we find that the probability of the correct assignment converges faster than for same-labeled loci, requiring observation of fewer independent chromatin configurations to establish loci identities. Finally, for a modified Rouse-model polymer, that realizes a population of dynamic loops, we find that the success probability also converges to unity exponentially as the number of observed loci configurations increases, albeit slightly more slowly than for a classical Rouse model polymer. Altogether, these results establish particle tracking of multiple identically- or alternately-labeled loci over time as a feasible way to infer temporal dynamics of the coarse-grained configuration of the chromatin polymer in individual living cells.

The spatial choreography of mRNA biosynthesis

Properly timed gene expression is essential for all aspects of organismal physiology. Despite significant progress, our understanding of the complex mechanisms governing the dynamics of gene regulation in response to internal and external signals remains incomplete. Over the past decade, advances in technologies like light and cryo-electron microscopy (Cryo-EM), cryo-electron tomography (Cryo-ET) and high-throughput sequencing have spurred new insights into traditional paradigms of gene expression. In this Review, we delve into recent concepts addressing 'where' and 'when' gene transcription and RNA splicing occur within cells, emphasizing the dynamic spatial and temporal organization of the cell nucleus.

The role of transcription bodies in gene expression: what embryos teach us

Transcription does not occur diffusely throughout the nucleus but is concentrated in specific areas. Areas of accumulated transcriptional machinery have been called clusters, hubs, or condensates, while transcriptionally active areas have been referred to as transcription factories or transcription bodies. Despite the widespread occurrence of transcription bodies, it has been difficult to study their assembly, function, and effect on gene expression. This review highlights the advantages of developmental model systems such as zebrafish and fruit fly embryos, in addressing these questions. We focus on three important discoveries that were made in embryos. (i) It had previously been suggested that, in transcription bodies, the different steps of the transcription process are organized in space. We explore how work in embryos has revealed that they can also be organized in time. In this case, transcription bodies mature from transcription factor clusters to elongating transcription bodies. This type of organization has important implications for transcription body function. (ii) The relevance of clustering for in vivo gene regulation has benefited greatly from studies in embryos. We discuss examples in which transcription bodies regulate developmental gene expression by compensating for low transcription factor concentrations and low-affinity enhancers. Finally, (iii) while accumulations of transcriptional machinery can facilitate transcription locally, work in embryos showed that transcription bodies can also sequester the transcriptional machinery, modulating the availability for activity at other sites. In brief, the reviewed literature highlights the properties of developmental model organisms that make them powerful systems for uncovering the form and function of transcription bodies.

Exploring the interplay between enhancer-promoter interactions and transcription

Enhancers in metazoan genomes are known to activate their target genes across both short and long genomic distances. Recent advances in chromosome conformation capture assays and single-cell imaging have shed light on the underlying chromatin contacts and dynamics. Yet the relationship between 3D physical enhancer-promoter (E-P) interactions and transcriptional activation remains unresolved. In this brief review, we discuss recent studies exploring this relationship across scales: from developmental stages to the minutes surrounding transcriptional activation and from the tissue level to single-allele subcellular dynamics. We discuss how seemingly contradictory observations might be reconciled and contribute to a refined causal relationship between E-P interactions and transcription, with mutual influences.

Bridging spatial and temporal scales of developmental gene regulation

The development of multicellular organisms relies on the precise coordination of molecular events across multiple spatial and temporal scales. Understanding how information flows from molecular interactions to cellular processes and tissue organization during development is crucial for explaining the remarkable reproducibility of complex organisms. This review explores how chromatin-encoded information is transduced from localized transcriptional events to global gene expression patterns, highlighting the challenge of bridging these scales. We discuss recent experimental findings and theoretical frameworks, emphasizing polymer physics as a tool for describing the relationship between chromatin structure and dynamics across scales. By integrating these perspectives, we aim to clarify how gene regulation is coordinated across levels of biological organization and suggest strategies for future experimental approaches.

Transcription factor clusters as information transfer agents

Deciphering how genes interpret information from transcription factor (TF) concentrations within the cell nucleus remains a fundamental question in gene regulation. Recent advancements have revealed the heterogeneous distribution of TF molecules, posing challenges to precisely decoding concentration signals. Using high-resolution single-cell imaging of the fluorescently tagged TF Bicoid in living Drosophila embryos, we show that Bicoid accumulation in submicrometer clusters preserves the spatial information of the maternal Bicoid gradient. These clusters provide precise spatial cues through intensity, size, and frequency. We further discover that Bicoid target genes colocalize with these clusters in an enhancer-binding affinity-dependent manner. Our modeling suggests that clustering offers a faster sensing mechanism for global nuclear concentrations than freely diffusing TF molecules detected by simple enhancers.

Binding profiles for 961 Drosophila and C. elegans transcription factors reveal tissue-specific regulatory relationships

A catalog of transcription factor (TF) binding sites in the genome is critical for deciphering regulatory relationships. Here, we present the culmination of the efforts of the modENCODE (model organism Encyclopedia of DNA Elements) and modERN (model organism Encyclopedia of Regulatory Networks) consortia to systematically assay TF binding events in vivo in two major model organisms, Drosophila melanogaster (fly) and Caenorhabditis elegans (worm). These data sets comprise 605 TFs identifying 3.6 M sites in the fly and 356 TFs identifying 0.9 M sites in the worm, and represent the majority of the regulatory space in each genome. We demonstrate that TFs associate with chromatin in clusters termed "metapeaks," that larger metapeaks have characteristics of high-occupancy target (HOT) regions, and that the importance of consensus sequence motifs bound by TFs depends on metapeak size and complexity. Combining ChIP-seq data with single-cell RNA-seq data in a machine-learning model identifies TFs with a prominent role in promoting target gene expression in specific cell types, even differentiating between parent-daughter cells during embryogenesis. These data are a rich resource for the community that should fuel and guide future investigations into TF function. To facilitate data accessibility and utility, all strains expressing green fluorescent protein (GFP)-tagged TFs are available at the stock centers for each organism. The chromatin immunoprecipitation sequencing data are available through the ENCODE Data Coordinating Center, GEO, and through a direct interface that provides rapid access to processed data sets and summary analyses, as well as widgets to probe the cell-type-specific TF-target relationships.

Three-dimensional regulatory hubs support oncogenic programs in glioblastoma

Dysregulation of enhancer-promoter communication in the context of the three-dimensional (3D) nucleus is increasingly recognized as a potential driver of oncogenic programs. Here, we profiled the 3D enhancer-promoter networks of primary patient-derived glioblastoma stem cells (GSCs) in comparison with neuronal stem cells (NSCs) to identify potential central nodes and vulnerabilities in the regulatory logic of this devastating cancer. Specifically, we focused on hyperconnected 3D regulatory hubs and demonstrated that hub-interacting genes exhibit high and coordinated expression at the single-cell level and strong association with oncogenic programs that distinguish IDH-wt glioblastoma patients from low-grade glioma. Epigenetic silencing of a recurrent 3D enhancer hub-with an uncharacterized role in glioblastoma-was sufficient to cause concordant downregulation of multiple hub-connected genes along with significant shifts in transcriptional states and reduced clonogenicity. By integrating published datasets from other cancer types, we also identified both universal and cancer type-specific 3D regulatory hubs which enrich for varying oncogenic programs and nominate specific factors associated with worse outcomes. Genetic alterations, such as focal duplications, could explain only a small fraction of the detected hyperconnected hubs and their increased activity. Overall, our study provides computational and experimental support for the potential central role of 3D regulatory hubs in controlling oncogenic programs and properties.

Gene mobility elements mediate cell type specific genome organization and radial gene movement in vivo

Understanding the level of genome organization that governs gene regulation remains a challenge despite advancements in chromatin profiling techniques. Cell type specific chromatin architectures may be obscured by averaging heterogeneous cell populations. Here we took a reductionist perspective, starting with the relocation of the hunchback gene to the nuclear lamina in Drosophila neuroblasts. We previously found that this event terminates competence to produce early-born neurons and is mediated by an intronic 250 base-pair element, which we term gene mobility element (GME). Here we found over 800 putative GMEs globally that are chromatin accessible and are Polycomb (PcG) target sites. GMEs appear to be distinct from PcG response elements, however, which are largely chromatin inaccessible in neuroblasts. Performing in situ Hi-C of purified neuroblasts, we found that GMEs form megabase-scale chromatin interactions, spanning multiple topologically associated domain borders, preferentially contacting other GMEs. These interactions are cell type and stage-specific. Notably, GMEs undergo developmentally- timed mobilization to/from the neuroblast nuclear lamina, and domain swapping a GFP reporter transgene intron with a GME relocates the transgene to the nuclear lamina in embryos. We propose that GMEs constitute a genome organizational framework and mediate gene-to-lamina mobilization during progenitor competence state transitions in vivo .

The cells are all-right: Regulation of the Lefty genes by separate enhancers in mouse embryonic stem cells

Enhancers play a critical role in regulating precise gene expression patterns essential for development and cellular identity; however, how gene-enhancer specificity is encoded within the genome is not clearly defined. To investigate how this specificity arises within topologically associated domains (TAD), we performed allele-specific genome editing of sequences surrounding the Lefty1 and Lefty2 paralogs in mouse embryonic stem cells. The Lefty genes arose from a tandem duplication event and these genes interact with each other in chromosome conformation capture assays which place these genes within the same TAD. Despite their physical proximity, we demonstrate that these genes are primarily regulated by separate enhancer elements. Through CRISPR-Cas9 mediated deletions to remove the intervening chromatin between the Lefty genes, we reveal a distance-dependent dosage effect of the Lefty2 enhancer on Lefty1 expression. These findings indicate a role for chromatin distance in insulating gene expression domains in the Lefty locus in the absence of architectural insulation.

Cohesin distribution alone predicts chromatin organization in yeast via conserved-current loop extrusion

BACKGROUND

Inhomogeneous patterns of chromatin-chromatin contacts within 10-100-kb-sized regions of the genome are a generic feature of chromatin spatial organization. These features, termed topologically associating domains (TADs), have led to the loop extrusion factor (LEF) model. Currently, our ability to model TADs relies on the observation that in vertebrates TAD boundaries are correlated with DNA sequences that bind CTCF, which therefore is inferred to block loop extrusion. However, although TADs feature prominently in their Hi-C maps, non-vertebrate eukaryotes either do not express CTCF or show few TAD boundaries that correlate with CTCF sites. In all of these organisms, the counterparts of CTCF remain unknown, frustrating comparisons between Hi-C data and simulations.

RESULTS

To extend the LEF model across the tree of life, here, we propose the conserved-current loop extrusion (CCLE) model that interprets loop-extruding cohesin as a nearly conserved probability current. From cohesin ChIP-seq data alone, we derive a position-dependent loop extrusion rate, allowing for a modified paradigm for loop extrusion, that goes beyond solely localized barriers to also include loop extrusion rates that vary continuously. We show that CCLE accurately predicts the TAD-scale Hi-C maps of interphase Schizosaccharomyces pombe, as well as those of meiotic and mitotic Saccharomyces cerevisiae, demonstrating its utility in organisms lacking CTCF.

CONCLUSIONS

The success of CCLE in yeasts suggests that loop extrusion by cohesin is indeed the primary mechanism underlying TADs in these systems. CCLE allows us to obtain loop extrusion parameters such as the LEF density and processivity, which compare well to independent estimates.

Transcription factor clusters as information transfer agents

Deciphering how genes interpret information from the concentration of transcription factors (TFs) within the cell nucleus remains a fundamental question in gene regulation. Recent advancements have unveiled the heterogeneous distribution of TF molecules in the nucleus, posing challenges to the precise decoding of concentration signals. To explore this phenomenon, we employ high-resolution single-cell imaging of a fluorescently tagged TF protein, Bicoid, in living fly embryos. We show that accumulation of Bicoid in submicron clusters preserves the spatial information of the maternal Bicoid gradient, and that cluster intensity, size, and frequency offer remarkably precise spatial cues. We further discover that various known gene targets of Bicoid activation colocalize with clusters and that for the target gene Hunchback, this colocalization is dependent on its enhancer binding affinity. Modeling information transfer through these clusters suggests that clustering offers a more rapid sensing mechanism for global nuclear concentrations than freely diffusing TF molecules detected by simple enhancers.

Enhancers: A Focus on Synthetic Biology and Correlated Gene Expression

Enhancers are central for the regulation of metazoan transcription but have proven difficult to study, primarily due to a myriad of interdependent variables shaping their activity. Consequently, synthetic biology has emerged as the main approach for dissecting mechanisms of enhancer function. We start by reviewing simple but highly parallel reporter assays, which have been successful in quantifying the complexity of the activator/coactivator mechanisms at enhancers. We then describe studies that examine how enhancers function in the genomic context and in combination with other enhancers, revealing that they activate genes through a variety of different mechanisms, working together as a system. Here, we primarily focus on synthetic reporter genes that can quantify the dynamics of enhancer biology through time. We end by considering the consequences of having many genes and enhancers within a 'local environment', which we believe leads to correlated gene expression and likely reports on the general principles of enhancer biology.

Single-cell new RNA sequencing reveals principles of transcription at the resolution of individual bursts

Analyses of transcriptional bursting from single-cell RNA-sequencing data have revealed patterns of variation and regulation in the kinetic parameters that could be inferred. Here we profiled newly transcribed (4-thiouridine-labelled) RNA across 10,000 individual primary mouse fibroblasts to more broadly infer bursting kinetics and coordination. We demonstrate that inference from new RNA profiles could separate the kinetic parameters that together specify the burst size, and that the synthesis rate (and not the transcriptional off rate) controls the burst size. Importantly, transcriptome-wide inference of transcriptional on and off rates provided conclusive evidence that RNA polymerase II transcribes genes in bursts. Recent reports identified examples of transcriptional co-bursting, yet no global analyses have been performed. The deep new RNA profiles we generated with allelic resolution demonstrated that co-bursting rarely appears more frequently than expected by chance, except for certain gene pairs, notably paralogues located in close genomic proximity. Altogether, new RNA single-cell profiling critically improves the inference of transcriptional bursting and provides strong evidence for independent transcriptional bursting of mammalian genes.

Regulatory landscape of enhancer-mediated transcriptional activation

Enhancers are noncoding regulatory elements that instruct spatial and temporal specificity of gene transcription in response to a variety of intrinsic and extrinsic signals during development. Although it has long been postulated that enhancers physically interact with target promoters through the formation of stable loops, recent studies have changed this static view: sequence-specific transcription factors (TFs) and coactivators are dynamically recruited to enhancers and assemble so-called transcription hubs. Dynamic assembly of transcription hubs appears to serve as a key scaffold to integrate regulatory information encoded by surrounding genome and biophysical properties of transcription machineries. In this review, we outline emerging new models of transcriptional regulation by enhancers and discuss future perspectives.

RNA polymerases reshape chromatin architecture and couple transcription on individual fibers

RNA polymerases must initiate and pause within a complex chromatin environment, surrounded by nucleosomes and other transcriptional machinery. This environment creates a spatial arrangement along individual chromatin fibers ripe for both competition and coordination, yet these relationships remain largely unknown owing to the inherent limitations of traditional structural and sequencing methodologies. To address this, we employed long-read chromatin fiber sequencing (Fiber-seq) in Drosophila to visualize RNA polymerase (Pol) within its native chromatin context with single-molecule precision along up to 30 kb fibers. We demonstrate that Fiber-seq enables the identification of individual Pol II, nucleosome, and transcription factor footprints, revealing Pol II pausing-driven destabilization of downstream nucleosomes. Furthermore, we demonstrate pervasive direct distance-dependent transcriptional coupling between nearby Pol II genes, Pol III genes, and transcribed enhancers, modulated by local chromatin architecture. Overall, transcription initiation reshapes surrounding nucleosome architecture and couples nearby transcriptional machinery along individual chromatin fibers.

Chromosome structure in Drosophila is determined by boundary pairing not loop extrusion

Two different models have been proposed to explain how the endpoints of chromatin looped domains ('TADs') in eukaryotic chromosomes are determined. In the first, a cohesin complex extrudes a loop until it encounters a boundary element roadblock, generating a stem-loop. In this model, boundaries are functionally autonomous: they have an intrinsic ability to halt the movement of incoming cohesin complexes that is independent of the properties of neighboring boundaries. In the second, loops are generated by boundary:boundary pairing. In this model, boundaries are functionally non-autonomous, and their ability to form a loop depends upon how well they match with their neighbors. Moreover, unlike the loop-extrusion model, pairing interactions can generate both stem-loops and circle-loops. We have used a combination of MicroC to analyze how TADs are organized, and experimental manipulations of the even skipped TAD boundary, homie, to test the predictions of the 'loop-extrusion' and the 'boundary-pairing' models. Our findings are incompatible with the loop-extrusion model, and instead suggest that the endpoints of TADs in flies are determined by a mechanism in which boundary elements physically pair with their partners, either head-to-head or head-to-tail, with varying degrees of specificity. Although our experiments do not address how partners find each other, the mechanism is unlikely to require loop extrusion.

Stem-loop and circle-loop TADs generated by directional pairing of boundary elements have distinct physical and regulatory properties

The chromosomes in multicellular eukaryotes are organized into a series of topologically independent loops called TADs. In flies, TADs are formed by physical interactions between neighboring boundaries. Fly boundaries exhibit distinct partner preferences, and pairing interactions between boundaries are typically orientation-dependent. Pairing can be head-to-tail or head-to-head. The former generates a stem-loop TAD, while the latter gives a circle-loop TAD. The TAD that encompasses the Drosophila even skipped (eve) gene is formed by the head-to-tail pairing of the nhomie and homie boundaries. To explore the relationship between loop topology and the physical and regulatory landscape, we flanked the nhomie boundary region with two attP sites. The attP sites were then used to generate four boundary replacements: λ DNA, nhomie forward (WT orientation), nhomie reverse (opposite of WT orientation), and homie forward (same orientation as WT homie). The nhomie forward replacement restores the WT physical and regulatory landscape: in MicroC experiments, the eve TAD is a 'volcano' triangle topped by a plume, and the eve gene and its regulatory elements are sequestered from interactions with neighbors. The λ DNA replacement lacks boundary function: the endpoint of the 'new' eve TAD on the nhomie side is ill-defined, and eve stripe enhancers activate a nearby gene, eIF3j. While nhomie reverse and homie forward restore the eve TAD, the topology is a circle-loop, and this changes the local physical and regulatory landscape. In MicroC experiments, the eve TAD interacts with its neighbors, and the plume at the top of the eve triangle peak is converted to a pair of 'clouds' of contacts with the next-door TADs. Consistent with the loss of isolation afforded by the stem-loop topology, the eve enhancers weakly activate genes in the neighboring TADs. Conversely, eve function is partially disrupted.

Molecular models of bidirectional promoter regulation

Approximately 11% of human genes are transcribed by a bidirectional promoter (BDP), defined as two genes with <1 kb between their transcription start sites. Despite their evolutionary conservation and enrichment for housekeeping genes and oncogenes, the regulatory role of BDPs remains unclear. BDPs have been suggested to facilitate gene coregulation and/or decrease expression noise. This review discusses these potential regulatory functions through the context of six prospective underlying mechanistic models: a single nucleosome free region, shared transcription factor/regulator binding, cooperative negative supercoiling, bimodal histone marks, joint activation by enhancer(s), and RNA-mediated recruitment of regulators. These molecular mechanisms may act independently and/or cooperatively to facilitate the coregulation and/or decreased expression noise predicted of BDPs.

How Transcription Factor Clusters Shape the Transcriptional Landscape

In eukaryotic cells, gene transcription typically occurs in discrete periods of promoter activity, interspersed with intervals of inactivity. This pattern deviates from simple stochastic events and warrants a closer examination of the molecular interactions that activate the promoter. Recent studies have identified transcription factor (TF) clusters as key precursors to transcriptional bursting. Often, these TF clusters form at chromatin segments that are physically distant from the promoter, making changes in chromatin conformation crucial for promoter-TF cluster interactions. In this review, I explore the formation and constituents of TF clusters, examining how the dynamic interplay between chromatin architecture and TF clustering influences transcriptional bursting. Additionally, I discuss techniques for visualizing TF clusters and provide an outlook on understanding the remaining gaps in this field.

New Drosophila promoter-associated architectural protein Mzfp1 interacts with CP190 and is required for housekeeping gene expression and insulator activity

In Drosophila, a group of zinc finger architectural proteins recruits the CP190 protein to the chromatin, an interaction that is essential for the functional activity of promoters and insulators. In this study, we describe a new architectural C2H2 protein called Madf and Zinc-Finger Protein 1 (Mzfp1) that interacts with CP190. Mzfp1 has an unusual structure that includes six C2H2 domains organized in a C-terminal cluster and two tandem MADF domains. Mzfp1 predominantly binds to housekeeping gene promoters located in both euchromatin and heterochromatin genome regions. In vivo mutagenesis studies showed that Mzfp1 is an essential protein, and both MADF domains and the CP190 interaction region are required for its functional activity. The C2H2 cluster is sufficient for the specific binding of Mzfp1 to regulatory elements, while the second MADF domain is required for Mzfp1 recruitment to heterochromatin. Mzfp1 binds to the proximal part of the Fub boundary that separates regulatory domains of the Ubx and abd-A genes in the Bithorax complex. Mzfp1 participates in Fub functions in cooperation with the architectural proteins Pita and Su(Hw). Thus, Mzfp1 is a new architectural C2H2 protein involved in the organization of active promoters and insulators in Drosophila.

Single-cell nascent RNA sequencing unveils coordinated global transcription

Transcription is the primary regulatory step in gene expression. Divergent transcription initiation from promoters and enhancers produces stable RNAs from genes and unstable RNAs from enhancers1,2. Nascent RNA capture and sequencing assays simultaneously measure gene and enhancer activity in cell populations3. However, fundamental questions about the temporal regulation of transcription and enhancer-gene coordination remain unanswered, primarily because of the absence of a single-cell perspective on active transcription. In this study, we present scGRO-seq-a new single-cell nascent RNA sequencing assay that uses click chemistry-and unveil coordinated transcription throughout the genome. We demonstrate the episodic nature of transcription and the co-transcription of functionally related genes. scGRO-seq can estimate burst size and frequency by directly quantifying transcribing RNA polymerases in individual cells and can leverage replication-dependent non-polyadenylated histone gene transcription to elucidate cell cycle dynamics. The single-nucleotide spatial and temporal resolution of scGRO-seq enables the identification of networks of enhancers and genes. Our results suggest that the bursting of transcription at super-enhancers precedes bursting from associated genes. By imparting insights into the dynamic nature of global transcription and the origin and propagation of transcription signals, we demonstrate the ability of scGRO-seq to investigate the mechanisms of transcription regulation and the role of enhancers in gene expression.

Exploring the reciprocity between pioneer factors and development

Development is regulated by coordinated changes in gene expression. Control of these changes in expression is largely governed by the binding of transcription factors to specific regulatory elements. However, the packaging of DNA into chromatin prevents the binding of many transcription factors. Pioneer factors overcome this barrier owing to unique properties that enable them to bind closed chromatin, promote accessibility and, in so doing, mediate binding of additional factors that activate gene expression. Because of these properties, pioneer factors act at the top of gene-regulatory networks and drive developmental transitions. Despite the ability to bind target motifs in closed chromatin, pioneer factors have cell type-specific chromatin occupancy and activity. Thus, developmental context clearly shapes pioneer-factor function. Here, we discuss this reciprocal interplay between pioneer factors and development: how pioneer factors control changes in cell fate and how cellular environment influences pioneer-factor binding and activity.

Gene activity fully predicts transcriptional bursting dynamics

Transcription commonly occurs in bursts, with alternating productive (ON) and quiescent (OFF) periods, governing mRNA production rates. Yet, how transcription is regulated through bursting dynamics remains unresolved. Here, we conduct real-time measurements of endogenous transcriptional bursting with single-mRNA sensitivity. Leveraging the diverse transcriptional activities in early fly embryos, we uncover stringent relationships between bursting parameters. Specifically, we find that the durations of ON and OFF periods are linked. Regardless of the developmental stage or body-axis position, gene activity levels predict individual alleles' average ON and OFF periods. Lowly transcribing alleles predominantly modulate OFF periods (burst frequency), while highly transcribing alleles primarily tune ON periods (burst size). These relationships persist even under perturbations of cis-regulatory elements or trans-factors and account for bursting dynamics measured in other species. Our results suggest a novel mechanistic constraint governing bursting dynamics rather than a modular control of distinct parameters by distinct regulatory processes.

Dynamic microenvironments shape nuclear organization and gene expression

Live imaging has revealed that the regulation of gene expression is largely driven by transient interactions. For example, many regulatory proteins bind chromatin for just seconds, and loop-like genomic contacts are rare and last only minutes. These discoveries have been difficult to reconcile with our canonical models that are predicated on stable and hierarchical interactions. Proteomic microenvironments that concentrate nuclear factors may explain how brief interactions can still mediate gene regulation by creating conditions where reactions occur more frequently. Here, we summarize new imaging technologies and recent discoveries implicating microenvironments as a potential driver of nuclear function. Finally, we propose that key properties of proteomic microenvironments, such as their size, enrichment, and lifetimes, are directly linked to regulatory function.

Conserved Cis-Acting Range Extender Element Mediates Extreme Long-Range Enhancer Activity in Mammals

While most mammalian enhancers regulate their cognate promoters over moderate distances of tens of kilobases (kb), some enhancers act over distances in the megabase range. The sequence features enabling such extreme-distance enhancer-promoter interactions remain elusive. Here, we used in vivo enhancer replacement experiments in mice to show that short- and medium-range enhancers cannot initiate gene expression at extreme-distance range. We uncover a novel conserved cis-acting element, Range EXtender (REX), that confers extreme-distance regulatory activity and is located next to a long-range enhancer of Sall1. The REX element itself has no endogenous enhancer activity. However, addition of the REX to other short- and mid-range enhancers substantially increases their genomic interaction range. In the most extreme example observed, addition of the REX increased the range of an enhancer by an order of magnitude, from its native 71kb to 840kb. The REX element contains highly conserved [C/T]AATTA homeodomain motifs. These motifs are enriched around long-range limb enhancers genome-wide, including the ZRS, a benchmark long-range limb enhancer of Shh. Mutating the [C/T]AATTA motifs within the ZRS does not affect its limb-specific enhancer activity at short range, but selectively abolishes its long-range activity, resulting in severe limb reduction in knock-in mice. In summary, we identify a sequence signature globally associated with long-range enhancer-promoter interactions and describe a prototypical REX element that is necessary and sufficient to confer extreme-distance gene activation by remote enhancers.

Polymer Collapse & Liquid-Liquid Phase-Separation are Coupled in a Generalized Prewetting Transition

The three-dimensional organization of chromatin is thought to play an important role in controlling gene expression. Specificity in expression is achieved through the interaction of transcription factors and other nuclear proteins with particular sequences of DNA. At unphysiological concentrations many of these nuclear proteins can phase-separate in the absence of DNA, and it has been hypothesized that, in vivo, the thermodynamic forces driving these phases help determine chromosomal organization. However it is unclear how DNA, itself a long polymer subject to configurational transitions, interacts with three-dimensional protein phases. Here we show that a long compressible polymer can be coupled to interacting protein mixtures, leading to a generalized prewetting transition where polymer collapse is coincident with a locally stabilized liquid droplet. We use lattice Monte-Carlo simulations and a mean-field theory to show that these phases can be stable even in regimes where both polymer collapse and coexisting liquid phases are unstable in isolation, and that these new transitions can be either abrupt or continuous. For polymers with internal linear structure we further show that changes in the concentration of bulk components can lead to changes in three-dimensional polymer structure. In the nucleus there are many distinct proteins that interact with many different regions of chromatin, potentially giving rise to many different Prewet phases. The simple systems we consider here highlight chromatin's role as a lower-dimensional surface whose interactions with proteins are required for these novel phases.

Polymer Collapse & Liquid-Liquid Phase-Separation are Coupled in a Generalized Prewetting Transition

The three-dimensional organization of chromatin is thought to play an important role in controlling gene expression. Specificity in expression is achieved through the interaction of transcription factors and other nuclear proteins with particular sequences of DNA. At unphysiological concentrations many of these nuclear proteins can phase-separate in the absence of DNA, and it has been hypothesized that, in vivo, the thermodynamic forces driving these phases help determine chromosomal organization. However it is unclear how DNA, itself a long polymer subject to configurational transitions, interacts with three-dimensional protein phases. Here we show that a long compressible polymer can be coupled to interacting protein mixtures, leading to a generalized prewetting transition where polymer collapse is coincident with a locally stabilized liquid droplet. We use lattice Monte-Carlo simulations and a mean-field theory to show that these phases can be stable even in regimes where both polymer collapse and coexisting liquid phases are unstable in isolation, and that these new transitions can be either abrupt or continuous. For polymers with internal linear structure we further show that changes in the concentration of bulk components can lead to changes in three-dimensional polymer structure. In the nucleus there are many distinct proteins that interact with many different regions of chromatin, potentially giving rise to many different Prewet phases. The simple systems we consider here highlight chromatin's role as a lower-dimensional surface whose interactions with proteins are required for these novel phases.

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.

Real-time single-molecule imaging of transcriptional regulatory networks in living cells

Gene regulatory networks drive the specific transcriptional programmes responsible for the diversification of cell types during the development of multicellular organisms. Although our knowledge of the genes involved in these dynamic networks has expanded rapidly, our understanding of how transcription is spatiotemporally regulated at the molecular level over a wide range of timescales in the small volume of the nucleus remains limited. Over the past few decades, advances in the field of single-molecule fluorescence imaging have enabled real-time behaviours of individual transcriptional components to be measured in living cells and organisms. These efforts are now shedding light on the dynamic mechanisms of transcription, revealing not only the temporal rules but also the spatial coordination of underlying molecular interactions during various biological events.

Enhancer-promoter specificity in gene transcription: molecular mechanisms and disease associations

Although often located at a distance from their target gene promoters, enhancers are the primary genomic determinants of temporal and spatial transcriptional specificity in metazoans. Since the discovery of the first enhancer element in simian virus 40, there has been substantial interest in unraveling the mechanism(s) by which enhancers communicate with their partner promoters to ensure proper gene expression. These research efforts have benefited considerably from the application of increasingly sophisticated sequencing- and imaging-based approaches in conjunction with innovative (epi)genome-editing technologies; however, despite various proposed models, the principles of enhancer-promoter interaction have still not been fully elucidated. In this review, we provide an overview of recent progress in the eukaryotic gene transcription field pertaining to enhancer-promoter specificity. A better understanding of the mechanistic basis of lineage- and context-dependent enhancer-promoter engagement, along with the continued identification of functional enhancers, will provide key insights into the spatiotemporal control of gene expression that can reveal therapeutic opportunities for a range of enhancer-related diseases.

Chromatin gene-gene loops support the cross-regulation of genes with related function

Chromatin loops between gene pairs have been observed in diverse contexts in both flies and vertebrates. Combining high-resolution Capture-C, DNA fluorescence in situ hybridization, and genetic perturbations, we dissect the functional role of three loops between genes with related function during Drosophila embryogenesis. By mutating the loop anchor (but not the gene) or the gene (but not loop anchor), we disentangle loop formation and gene expression and show that the 3D proximity of paralogous gene loci supports their co-regulation. Breaking the loop leads to either an attenuation or enhancement of expression and perturbs their relative levels of expression and cross-regulation. Although many loops appear constitutive across embryogenesis, their function can change in different developmental contexts. Taken together, our results indicate that chromatin gene-gene loops act as architectural scaffolds that can be used in different ways in different contexts to fine-tune the coordinated expression of genes with related functions and sustain their cross-regulation.

A Hierarchical Mechanotransduction System: From Macro to Micro

Mechanotransduction is a strictly regulated process whereby mechanical stimuli, including mechanical forces and properties, are sensed and translated into biochemical signals. Increasing data demonstrate that mechanotransduction is crucial for regulating macroscopic and microscopic dynamics and functionalities. However, the actions and mechanisms of mechanotransduction across multiple hierarchies, from molecules, subcellular structures, cells, tissues/organs, to the whole-body level, have not been yet comprehensively documented. Herein, the biological roles and operational mechanisms of mechanotransduction from macro to micro are revisited, with a focus on the orchestrations across diverse hierarchies. The implications, applications, and challenges of mechanotransduction in human diseases are also summarized and discussed. Together, this knowledge from a hierarchical perspective has the potential to refresh insights into mechanotransduction regulation and disease pathogenesis and therapy, and ultimately revolutionize the prevention, diagnosis, and treatment of human diseases.

Opticool: Cutting-edge transgenic optical tools

Only a few short decades have passed since the sequencing of GFP, yet the modern repertoire of transgenically encoded optical tools implies an exponential proliferation of ever improving constructions to interrogate the subcellular environment. A myriad of tags for labeling proteins, RNA, or DNA have arisen in the last few decades, facilitating unprecedented visualization of subcellular components and processes. Development of a broad array of modern genetically encoded sensors allows real-time, in vivo detection of molecule levels, pH, forces, enzyme activity, and other subcellular and extracellular phenomena in ever expanding contexts. Optogenetic, genetically encoded optically controlled manipulation systems have gained traction in the biological research community and facilitate single-cell, real-time modulation of protein function in vivo in ever broadening, novel applications. While this field continues to explosively expand, references are needed to assist scientists seeking to use and improve these transgenic devices in new and exciting ways to interrogate development and disease. In this review, we endeavor to highlight the state and trajectory of the field of in vivo transgenic optical tools.

The landscape of super-enhancer regulates remote target gene transcription through loop domains in adipose tissue of pig

BACKGROUND

A super-enhancer (SE) is a huge cluster of multiple enhancers that control the key genes for cell identity and function. The rise of advanced chromatin immunoprecipitation sequencing (ChIP-seq) technology such as Cleavage Under Targets and Tagmentation (CUT&Tag) allows more SEs to be discovered. However, SE studies in Luchuan and Duroc pigs are very rare in animal husbandry.

RESULTS

We used the CUT&Tag technique to identify 145 and 378 SEs from the adipose tissues of Luchuan and Duroc pigs, respectively. There were significant differences in the peak coverage ratio of SE peaks in the gene promoter region between the two breeds. Not only that, peak signals at the start and end point of the SE peak profile showed obvious spikes. The proximal target genes of SE were highly expressed compared with the background genes and the typical enhancer target genes. Subsequently, in conjoint analysis with high-throughput chromosome conformation capture sequencing (Hi-C seq) data, we predicted the remote regulatory genes of SE and found that their expression level was related to the distance of SE extended to the loop's anchor, but not the length of loops. According to our prediction model, SEs can maintain promoter accessibility of partial remote target genes through loop domains. Finally, a batch of SEs closely related to fat metabolism traits were obtained by performing a coalition analysis of quantitative trait loci and SE data.

CONCLUSIONS

This work enabled us to obtain hundreds of SEs from Luchuan and Duroc pigs. Our model provides a new method for predicting the SE remote target genes based on loop domains, and to further explore the potential role of super-enhancer in the regulation of fat metabolism.

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.

What are tethering elements?

High-resolution Micro-C maps identified a specialized class of regulatory DNAs termed 'tethering elements' (TEs) in Drosophila. These 300-500-bp elements facilitate specific long-range genomic associations or loops. The POZ-containing transcription factor GAF (GAGA-associated factor) contributes to loop formation. Tether-tether interactions accelerate Hox gene activation by distal enhancers, and coordinate transcription of duplicated genes (paralogs) through promoter-promoter associations. Some TEs engage in ultra-long-range enhancer-promoter and promoter-promoter interactions (meta-loops) in the Drosophila brain. We discuss the basis for tether-tether specificity and speculate on the occurrence of similar elements in vertebrate genomes.

Binding profiles for 954 Drosophila and C. elegans transcription factors reveal tissue specific regulatory relationships

A catalog of transcription factor (TF) binding sites in the genome is critical for deciphering regulatory relationships. Here we present the culmination of the modERN (model organism Encyclopedia of Regulatory Networks) consortium that systematically assayed TF binding events in vivo in two major model organisms, Drosophila melanogaster (fly) and Caenorhabditis elegans (worm). We describe key features of these datasets, comprising 604 TFs identifying 3.6M sites in the fly and 350 TFs identifying 0.9 M sites in the worm. Applying a machine learning model to these data identifies sets of TFs with a prominent role in promoting target gene expression in specific cell types. TF binding data are available through the ENCODE Data Coordinating Center and at https://epic.gs.washington.edu/modERNresource, which provides access to processed and summary data, as well as widgets to probe cell type-specific TF-target relationships. These data are a rich resource that should fuel investigations into TF function during 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.

Functional analysis of the Drosophila eve locus in response to non-canonical combinations of gap gene expression levels

Transcription factor combinations play a key role in shaping cellular identity. However, the precise relationship between specific combinations and downstream effects remains elusive. Here, we investigate this relationship within the context of the Drosophila eve locus, which is controlled by gap genes. We measure spatiotemporal levels of four gap genes in heterozygous and homozygous gap mutant embryos and correlate them with the striped eve activity pattern. Although changes in gap gene expression extend beyond the manipulated gene, the spatial patterns of Eve expression closely mirror canonical activation levels in wild type. Interestingly, some combinations deviate from the wild-type repertoire but still drive eve activation. Although in homozygous mutants some Eve stripes exhibit partial penetrance, stripes consistently emerge at reproducible positions, even with varying gap gene levels. Our findings suggest a robust molecular canalization of cell fates in gap mutants and provide insights into the regulatory constraints governing multi-enhancer gene loci.