Multi-enhancer transcriptional hubs confer phenotypic robustness

Transcription Factor Condensates Mediate Clustering of MET Regulon and Enhancement in Gene Expression

Some transcription factors (TFs) can form liquid-liquid phase separated (LLPS) condensates. However, the functions of these TF condensates in 3D genome organization and gene regulation remain elusive. In response to methionine (met) starvation, budding yeast TF Met4 and a few co-activators, including Met32, induce a set of genes involved in met biosynthesis. Here, we show that the endogenous Met4 and Met32 form co-localized puncta-like structures in yeast nuclei upon met depletion. Recombinant Met4 and Met32 form mixed droplets with LLPS properties in vitro. In relation to chromatin, Met4 puncta co-localize with target genes, and at least a subset of these target genes is clustered in 3D in a Met4-dependent manner. A MET3pr-GFP reporter inserted near several native Met4 binding sites becomes co-localized with Met4 puncta and displays enhanced transcriptional activity. A Met4 variant with a partial truncation of an intrinsically disordered region (IDR) shows less puncta formation, and this mutant selectively reduces the reporter activity near Met4 binding sites to the basal level. Overall, these results support a model where Met4 and co-activators form condensates to bring multiple target genes into a vicinity with higher local TF concentrations, which facilitates a strong response to methionine depletion.

Enhancer selectivity across cell types delineates three functionally distinct enhancer-promoter regulation patterns

BACKGROUND

Multiple enhancers co-regulating the same gene is prevalent and plays a crucial role during development and disease. However, how multiple enhancers coordinate the same gene expression across various cell types remains largely unexplored at genome scale.

RESULTS

We develop a computational approach that enables the quantitative assessment of enhancer specificity and selectivity across diverse cell types, leveraging enhancer-promoter (E-P) interactions data. We observe two well-known gene regulation patterns controlled by enhancer clusters, which regulate the same gene either in a limited number of cell types (Specific pattern, Spe) or in the majority of cell types (Conserved pattern, Con), both of which are enriched for super-enhancers (SEs). We identify a previously overlooked pattern (Variable pattern, Var) that multiple enhancers link to the same gene, but rarely coexist in the same cell type. These three patterns control the genes associating with distinct biological function and exhibit unique epigenetic features. Specifically, we discover a subset of Var patterns contains Shared enhancers with stable enhancer-promoter interactions in the majority of cell types, which might contribute to maintaining gene expression by recruiting abundant CTCF.

CONCLUSIONS

Together, our findings reveal three distinct E-P regulation patterns across different cell types, providing insights into deciphering the complexity of gene transcriptional regulation.

Dynamic modes of Notch transcription hubs conferring memory and stochastic activation revealed by live imaging the co-activator Mastermind

Developmental programming involves the accurate conversion of signalling levels and dynamics to transcriptional outputs. The transcriptional relay in the Notch pathway relies on nuclear complexes containing the co-activator Mastermind (Mam). By tracking these complexes in real time, we reveal that they promote the formation of a dynamic transcription hub in Notch ON nuclei which concentrates key factors including the Mediator CDK module. The composition of the hub is labile and persists after Notch withdrawal conferring a memory that enables rapid reformation. Surprisingly, only a third of Notch ON hubs progress to a state with nascent transcription, which correlates with polymerase II and core Mediator recruitment. This probability is increased by a second signal. The discovery that target-gene transcription is probabilistic has far-reaching implications because it implies that stochastic differences in Notch pathway output can arise downstream of receptor activation.

A fine kinetic balance of interactions directs transcription factor hubs to genes

Eukaryotic gene regulation relies on the binding of sequence-specific transcription factors (TFs). TFs bind chromatin transiently yet occupy their target sites by forming high-local concentration microenvironments (hubs and condensates) that increase the frequency of binding events. Despite their ubiquity, such microenvironments have been difficult to study in endogenous contexts due to technical limitations. Here, we overcome these limitations and investigate how hubs drive TF occupancy at their targets. Using a DNA binding perturbation to a hub-forming TF, Zelda, in Drosophila embryos, we find that hub properties, including the stability and frequencies of associations to targets, are key determinants of TF occupancy. Our data suggest that the targeting of these hubs is driven not just by specific DNA motif recognition, but also by a fine-tuned kinetic balance of interactions between TFs and their co-binding partners.

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.

Cell-type-directed design of synthetic enhancers

Transcriptional enhancers act as docking stations for combinations of transcription factors and thereby regulate spatiotemporal activation of their target genes1. It has been a long-standing goal in the field to decode the regulatory logic of an enhancer and to understand the details of how spatiotemporal gene expression is encoded in an enhancer sequence. Here we show that deep learning models2-6, can be used to efficiently design synthetic, cell-type-specific enhancers, starting from random sequences, and that this optimization process allows detailed tracing of enhancer features at single-nucleotide resolution. We evaluate the function of fully synthetic enhancers to specifically target Kenyon cells or glial cells in the fruit fly brain using transgenic animals. We further exploit enhancer design to create 'dual-code' enhancers that target two cell types and minimal enhancers smaller than 50 base pairs that are fully functional. By examining the state space searches towards local optima, we characterize enhancer codes through the strength, combination and arrangement of transcription factor activator and transcription factor repressor motifs. Finally, we apply the same strategies to successfully design human enhancers, which adhere to enhancer rules similar to those of Drosophila enhancers. Enhancer design guided by deep learning leads to better understanding of how enhancers work and shows that their code can be exploited to manipulate cell states.

Affinity-optimizing enhancer variants disrupt development

Enhancers control the location and timing of gene expression and contain the majority of variants associated with disease1-3. The ZRS is arguably the most well-studied vertebrate enhancer and mediates the expression of Shh in the developing limb4. Thirty-one human single-nucleotide variants (SNVs) within the ZRS are associated with polydactyly4-6. However, how this enhancer encodes tissue-specific activity, and the mechanisms by which SNVs alter the number of digits, are poorly understood. Here we show that the ETS sites within the ZRS are low affinity, and identify a functional ETS site, ETS-A, with extremely low affinity. Two human SNVs and a synthetic variant optimize the binding affinity of ETS-A subtly from 15% to around 25% relative to the strongest ETS binding sequence, and cause polydactyly with the same penetrance and severity. A greater increase in affinity results in phenotypes that are more penetrant and more severe. Affinity-optimizing SNVs in other ETS sites in the ZRS, as well as in ETS, interferon regulatory factor (IRF), HOX and activator protein 1 (AP-1) sites within a wide variety of enhancers, cause gain-of-function gene expression. The prevalence of binding sites with suboptimal affinity in enhancers creates a vulnerability in genomes whereby SNVs that optimize affinity, even slightly, can be pathogenic. Searching for affinity-optimizing SNVs in genomes could provide a mechanistic approach to identify causal variants that underlie enhanceropathies.

Homeodomain complex formation and biomolecular condensates in Hox gene regulation

Hox genes are a family of homeodomain transcription factors that regulate specialized morphological structures along the anterior-posterior axis of metazoans. Over the past few decades, researchers have focused on defining how Hox factors with similar in vitro DNA binding activities achieve sufficient target specificity to regulate distinct cell fates in vivo. In this review, we highlight how protein interactions with other transcription factors, many of which are also homeodomain proteins, result in the formation of transcription factor complexes with enhanced DNA binding specificity. These findings suggest that Hox-regulated enhancers utilize distinct combinations of homeodomain binding sites, many of which are low-affinity, to recruit specific Hox complexes. However, low-affinity sites can only yield reproducible responses with high transcription factor concentrations. To overcome this limitation, recent studies revealed how transcription factors, including Hox factors, use intrinsically disordered domains (IDRs) to form biomolecular condensates that increase protein concentrations. Moreover, Hox factors with altered IDRs have been associated with altered transcriptional activity and human disease states, demonstrating the importance of IDRs in mediating essential Hox output. Collectively, these studies highlight how Hox factors use their DNA binding domains, protein-protein interaction domains, and IDRs to form specific transcription factor complexes that yield accurate gene expression.

Transcription factor condensates and signaling driven transcription

Transcription Factor (TF) condensates are a heterogenous mix of RNA, DNA, and multiple co-factor proteins capable of modulating the transcriptional response of the cell. The dynamic nature and the spatial location of TF-condensates in the 3D nuclear space is believed to provide a fast response, which is on the same pace as the signaling cascade and yet ever-so-specific in the crowded environment of the nucleus. However, the current understanding of how TF-condensates can achieve these feet so quickly and efficiently is still unclear. In this review, we draw parallels with other protein condensates and share our speculations on how the nucleus uses these TF-condensates to achieve high transcriptional specificity and fidelity. We discuss the various constituents of TF-condensates, their properties, and the known and unknown functions of TF-condensates with a particular focus on steroid signaling-induced transcriptional programs.

Spatial promoter-enhancer hubs in cancer: organization, regulation, and function

Transcriptional dysregulation is a hallmark of cancer and can be driven by altered enhancer landscapes. Recent studies in genome organization have revealed that multiple enhancers and promoters can spatially coalesce to form dynamic topological assemblies, known as promoter-enhancer hubs, which strongly correlate with elevated gene expression. In this review, we discuss the structure and complexity of promoter-enhancer hubs recently identified in multiple cancer types. We further discuss underlying mechanisms driving dysregulation of promoter-enhancer hubs and speculate on their functional role in pathogenesis. Understanding the role of promoter-enhancer hubs in transcriptional dysregulation can provide insight into new therapeutic approaches to target these complex features of genome organization.

Integrative approaches based on genomic techniques in the functional studies on enhancers

With the development of sequencing technology and the dramatic drop in sequencing cost, the functions of noncoding genes are being characterized in a wide variety of fields (e.g. biomedicine). Enhancers are noncoding DNA elements with vital transcription regulation functions. Tens of thousands of enhancers have been identified in the human genome; however, the location, function, target genes and regulatory mechanisms of most enhancers have not been elucidated thus far. As high-throughput sequencing techniques have leapt forwards, omics approaches have been extensively employed in enhancer research. Multidimensional genomic data integration enables the full exploration of the data and provides novel perspectives for screening, identification and characterization of the function and regulatory mechanisms of unknown enhancers. However, multidimensional genomic data are still difficult to integrate genome wide due to complex varieties, massive amounts, high rarity, etc. To facilitate the appropriate methods for studying enhancers with high efficacy, we delineate the principles, data processing modes and progress of various omics approaches to study enhancers and summarize the applications of traditional machine learning and deep learning in multi-omics integration in the enhancer field. In addition, the challenges encountered during the integration of multiple omics data are addressed. Overall, this review provides a comprehensive foundation for enhancer analysis.

Transcript accumulation rates in the early Caenorhabditis elegans embryo

Dynamic transcriptional changes are widespread in rapidly dividing developing embryos when cell fate decisions are made quickly. The Caenorhabditis elegans embryo overcomes these constraints partly through the rapid production of high levels of transcription factor mRNAs. Transcript accumulation rates for some developmental genes are known at single-cell resolution, but genome-scale measurements are lacking. We estimate zygotic mRNA accumulation rates from single-cell RNA sequencing data calibrated with single-molecule transcript imaging. Rapid transcription is common in the early C. elegans embryo with rates highest soon after zygotic transcription begins. High-rate genes are enriched for recently duplicated cell-fate regulators and share common genomic features. We identify core promoter elements associated with high rate and measure their contributions for two early endomesodermal genes, ceh-51 and sdz-31. Individual motifs modestly affect accumulation rates, suggesting multifactorial control. These results are a step toward estimating absolute transcription kinetics and understanding how transcript dosage drives developmental decisions.

Quantitative control of Ets1 dosage by a multi-enhancer hub promotes Th1 cell differentiation and protects from allergic inflammation

Multi-enhancer hubs are spatial clusters of enhancers present across numerous developmental programs. Here, we studied the functional relevance of these three-dimensional structures in T cell biology. Mathematical modeling identified a highly connected multi-enhancer hub at the Ets1 locus, comprising a noncoding regulatory element that was a hotspot for sequence variation associated with allergic disease in humans. Deletion of this regulatory element in mice revealed that the multi-enhancer connectivity was dispensable for T cell development but required for CD4+ T helper 1 (Th1) differentiation. These mice were protected from Th1-mediated colitis but exhibited overt allergic responses. Mechanistically, the multi-enhancer hub controlled the dosage of Ets1 that was required for CTCF recruitment and assembly of Th1-specific genome topology. Our findings establish a paradigm wherein multi-enhancer hubs control cellular competence to respond to an inductive cue through quantitative control of gene dosage and provide insight into how sequence variation within noncoding elements at the Ets1 locus predisposes individuals to allergic responses.

Mutational scans reveal differential evolvability of Drosophila promoters and enhancers

Rapid enhancer and slow promoter evolution have been demonstrated through comparative genomics. However, it is not clear how this information is encoded genetically and if this can be used to place evolution in a predictive context. Part of the challenge is that our understanding of the potential for regulatory evolution is biased primarily toward natural variation or limited experimental perturbations. Here, to explore the evolutionary capacity of promoter variation, we surveyed an unbiased mutation library for three promoters in Drosophila melanogaster. We found that mutations in promoters had limited to no effect on spatial patterns of gene expression. Compared to developmental enhancers, promoters are more robust to mutations and have more access to mutations that can increase gene expression, suggesting that their low activity might be a result of selection. Consistent with these observations, increasing the promoter activity at the endogenous locus of shavenbaby led to increased transcription yet limited phenotypic changes. Taken together, developmental promoters may encode robust transcriptional outputs allowing evolvability through the integration of diverse developmental enhancers.

This article is part of the theme issue ‘Interdisciplinary approaches to predicting evolutionary biology’.

Shared retinoic acid responsive enhancers coordinately regulate nascent transcription of Hoxb coding and non-coding RNAs in the developing mouse neural tube

Signaling pathways regulate the patterns of Hox gene expression that underlie their functions in the specification of axial identity. Little is known about the properties of cis-regulatory elements and underlying transcriptional mechanisms that integrate graded signaling inputs to coordinately control Hox expression. Here, we optimized a single molecule fluorescent in situ hybridization (smFISH) technique with probes spanning introns to evaluate how three shared retinoic acid response element (RARE)-dependent enhancers in the Hoxb cluster regulate patterns of nascent transcription in vivo at the level of single cells in wild-type and mutant embryos. We predominately detect nascent transcription of only a single Hoxb gene in each cell, with no evidence for simultaneous co-transcriptional coupling of all or specific subsets of genes. Single and/or compound RARE mutations indicate that each enhancer differentially impacts global and local patterns of nascent transcription, suggesting that selectivity and competitive interactions between these enhancers is important to robustly maintain the proper levels and patterns of nascent Hoxb transcription. This implies that rapid and dynamic regulatory interactions potentiate transcription of genes through combined inputs from these enhancers in coordinating the retinoic acid response.

Decoding enhancer complexity with machine learning and high-throughput discovery

Enhancers are genomic DNA elements controlling spatiotemporal gene expression. Their flexible organization and functional redundancies make deciphering their sequence-function relationships challenging. This article provides an overview of the current understanding of enhancer organization and evolution, with an emphasis on factors that influence these relationships. Technological advancements, particularly in machine learning and synthetic biology, are discussed in light of how they provide new ways to understand this complexity. Exciting opportunities lie ahead as we continue to unravel the intricacies of enhancer function.

Connectome and regulatory hubs of CAGE highly active enhancers

Evidence indicates that enhancers are transcriptionally active. Herein, we investigated transcriptionally active enhancers by using cap analysis of gene expression (CAGE) combined with epigenetic marks and chromatin interactions. We identified CAGE-tag highly active (CHA) enhancers as distant regulatory elements with CAGE-tag ≥ 90th percentile and overlapping with H3K27ac peaks (4.5% of enhancers). CHA enhancers were conserved between mouse and man and were independent from super-enhancers in predicting cell identity with lower P-values. CHA enhancers had increased open chromatin and a higher recruitment of cell-specific transcription factors as well as molecules involved in 3D genome interactions. HiChIP analysis of enhancer-promoter looping indicated that CHA enhancers had a higher density of anchor loops when compared to regular enhancers. A subset of CHA enhancers and promoters characterized by a high density of chromatin loops and forming hub regulatory units were connected to the promoter of immediate early response genes, genes involved in cancer and encoding for transcription factors. Promoter of genes within hub CHA regulatory units were less likely to be paused. CHA enhancers were enriched in gene variants associated with autoimmune disorders and had looping with causal candidate genes as revealed by Mendelian randomization. Hence, CHA enhancers form a dense hierarchical network of chromatin interactions between regulatory elements and genes involved in cell identity and disorders.

Differential modularity of the mammalian Engrailed 1 enhancer network directs sweat gland development

Enhancers are context-specific regulators of expression that drive biological complexity and variation through the redeployment of conserved genes. An example of this is the enhancer-mediated control of Engrailed 1 (EN1), a pleiotropic gene whose expression is required for the formation of mammalian eccrine sweat glands. We previously identified the En1 candidate enhancer (ECE) 18 cis-regulatory element that has been highly and repeatedly derived on the human lineage to potentiate ectodermal EN1 and induce our species' uniquely high eccrine gland density. Intriguingly, ECE18 quantitative activity is negligible outside of primates and ECE18 is not required for En1 regulation and eccrine gland formation in mice, raising the possibility that distinct enhancers have evolved to modulate the same trait. Here we report the identification of the ECE20 enhancer and show it has conserved functionality in mouse and human developing skin ectoderm. Unlike ECE18, knock-out of ECE20 in mice reduces ectodermal En1 and eccrine gland number. Notably, we find ECE20, but not ECE18, is also required for En1 expression in the embryonic mouse brain, demonstrating that ECE20 is a pleiotropic En1 enhancer. Finally, that ECE18 deletion does not potentiate the eccrine phenotype of ECE20 knock-out mice supports the secondary incorporation of ECE18 into the regulation of this trait in primates. Our findings reveal that the mammalian En1 regulatory machinery diversified to incorporate both shared and lineage-restricted enhancers to regulate the same phenotype, and also have implications for understanding the forces that shape the robustness and evolvability of developmental traits.

Complexity of enhancer networks predicts cell identity and disease genes revealed by single-cell multi-omics analysis

Many enhancers exist as clusters in the genome and control cell identity and disease genes; however, the underlying mechanism remains largely unknown. Here, we introduce an algorithm, eNet, to build enhancer networks by integrating single-cell chromatin accessibility and gene expression profiles. The complexity of enhancer networks is assessed by two metrics: the number of enhancers and the frequency of predicted enhancer interactions (PEIs) based on chromatin co-accessibility. We apply eNet algorithm to a human blood dataset and find cell identity and disease genes tend to be regulated by complex enhancer networks. The network hub enhancers (enhancers with frequent PEIs) are the most functionally important. Compared with super-enhancers, enhancer networks show better performance in predicting cell identity and disease genes. eNet is robust and widely applicable in various human or mouse tissues datasets. Thus, we propose a model of enhancer networks containing three modes: Simple, Multiple and Complex, which are distinguished by their complexity in regulating gene expression. Taken together, our work provides an unsupervised approach to simultaneously identify key cell identity and disease genes and explore the underlying regulatory relationships among enhancers in single cells.

Developmental phenomics suggests that H3K4 monomethylation confers multi-level phenotypic robustness

How histone modifications affect animal development remains difficult to ascertain. Despite the prevalence of histone 3 lysine 4 monomethylation (H3K4me1) on enhancers, hypomethylation appears to have minor effects on phenotype and viability. Here, we genetically reduce H3K4me1 deposition in Drosophila melanogaster and find that hypomethylation reduces transcription factor enrichment in nuclear microenvironments, disrupts gene expression, and reduces phenotypic robustness. Using a developmental phenomics approach, we find changes in morphology, metabolism, behavior, and offspring production. However, many phenotypic changes are only detected when hypomethylated flies develop outside of standard laboratory environments or with specific genetic backgrounds. Therefore, quantitative phenomics measurements can unravel how pleiotropic modulators of gene expression affect developmental robustness under conditions resembling the natural environments of a species.

Unveiling dynamic enhancer–promoter interactions in Drosophila melanogaster

Proper enhancer–promoter interactions are essential to maintaining specific transcriptional patterns and preventing ectopic gene expression. Drosophila is an ideal model organism to study transcriptional regulation due to extensively characterized regulatory regions and the ease of implementing new genetic and molecular techniques for quantitative analysis. The mechanisms of enhancer–promoter interactions have been investigated over a range of length scales. At a DNA level, compositions of both enhancer and promoter sequences affect transcriptional dynamics, including duration, amplitude, and frequency of transcriptional bursting. 3D chromatin topology is also important for proper enhancer–promoter contacts. By working competitively or cooperatively with one another, multiple, simultaneous enhancer–enhancer, enhancer–promoter, and promoter–promoter interactions often occur to maintain appropriate levels of mRNAs. For some long-range enhancer–promoter interactions, extra regulatory elements like insulators and tethering elements are required to promote proper interactions while blocking aberrant ones. This review provides an overview of our current understanding of the mechanism of enhancer–promoter interactions and how perturbations of such interactions affect transcription and subsequent physiological outcomes.

Shadow enhancers modulate distinct transcriptional parameters that differentially effect downstream patterning events

Transcription in the early Drosophila blastoderm is coordinated by the collective action of hundreds of enhancers. Many genes are controlled by so-called 'shadow enhancers', which provide resilience to environment or genetic insult, allowing the embryo to robustly generate a precise transcriptional pattern. Emerging evidence suggests that many shadow enhancer pairs do not drive identical expression patterns, but the biological significance of this remains unclear. In this study, we characterize the shadow enhancer pair controlling the gene short gastrulation (sog). We removed either the intronic proximal enhancer or the upstream distal enhancer and monitored sog transcriptional kinetics. Notably, each enhancer differs in sog spatial expression, timing of activation and RNA Polymerase II loading rates. In addition, modeling of individual enhancer activities demonstrates that these enhancers integrate activation and repression signals differently. Whereas activation is due to the sum of the two enhancer activities, repression appears to depend on synergistic effects between enhancers. Finally, we examined the downstream signaling consequences resulting from the loss of either enhancer, and found changes in tissue patterning that can be explained by the differences in transcriptional kinetics measured.

Deep cis-regulatory homology of the butterfly wing pattern ground plan

Butterfly wing patterns derive from a deeply conserved developmental ground plan yet are diverse and evolve rapidly. It is poorly understood how gene regulatory architectures can accommodate both deep homology and adaptive change. To address this, we characterized the cis-regulatory evolution of the color pattern gene WntA in nymphalid butterflies. Comparative assay for transposase-accessible chromatin using sequencing (ATAC-seq) and in vivo deletions spanning 46 cis-regulatory elements across five species revealed deep homology of ground plan-determining sequences, except in monarch butterflies. Furthermore, noncoding deletions displayed both positive and negative regulatory effects that were often broad in nature. Our results provide little support for models predicting rapid enhancer turnover and suggest that deeply ancestral, multifunctional noncoding elements can underlie rapidly evolving trait systems.

Shared Transcriptional Machinery at Homologous Alleles Leads to Reduced Transcription in Early Drosophila Embryos

The mechanism by which transcriptional machinery is recruited to enhancers and promoters to regulate gene expression is one of the most challenging and extensively studied questions in modern biology. We explored the possibility that interallelic interactions between two homologous alleles might affect gene regulation. Using an MS2- and PP7-based, allele-specific live imaging assay, we visualized de novo transcripts of a reporter gene in hemizygous and homozygous Drosophila embryos. Surprisingly, each homozygous allele produced fewer RNAs than the corresponding hemizygous allele, suggesting the possibility of allelic competition in homozygotes. However, the competition was not observed when the enhancer-promoter interaction was weakened by placing the reporter construct in a different chromosome location or by moving the enhancer further away from the promoter. Moreover, the reporter gene showed reduced transcriptional activity when a partial transcription unit (either an enhancer or reporter gene only) was in the homologous position. We propose that the transcriptional machinery that binds both the enhancer and promoter regions, such as RNA Pol II or preinitiation complexes, may be responsible for the allelic competition. We showed that the degree of allelic interference increased over developmental time as more Pol II was needed to activate zygotic genes. Such allelic competition was observed for an endogenous gene as well. Our study provides new insights into the role of 3D interallelic interactions in gene regulation.

A synthetic synthesis to explore animal evolution and development

Identifying the general principles by which genotypes are converted into phenotypes remains a challenge in the post-genomic era. We still lack a predictive understanding of how genes shape interactions among cells and tissues in response to signalling and environmental cues, and hence how regulatory networks generate the phenotypic variation required for adaptive evolution. Here, we discuss how techniques borrowed from synthetic biology may facilitate a systematic exploration of evolvability across biological scales. Synthetic approaches permit controlled manipulation of both endogenous and fully engineered systems, providing a flexible platform for investigating causal mechanisms in vivo. Combining synthetic approaches with multi-level phenotyping (phenomics) will supply a detailed, quantitative characterization of how internal and external stimuli shape the morphology and behaviour of living organisms. We advocate integrating high-throughput experimental data with mathematical and computational techniques from a variety of disciplines in order to pursue a comprehensive theory of evolution.

This article is part of the theme issue ‘Genetic basis of adaptation and speciation: from loci to causative mutations’.

Transcriptional coupling of distant regulatory genes in living embryos

The prevailing view of metazoan gene regulation is that individual genes are independently regulated by their own dedicated sets of transcriptional enhancers. Past studies have reported long-range gene-gene associations^1-3, but their functional importance in regulating transcription remains unclear. Here we used quantitative single-cell live imaging methods to provide a demonstration of co-dependent transcriptional dynamics of genes separated by large genomic distances in living Drosophila embryos. We find extensive physical and functional associations of distant paralogous genes, including co-regulation by shared enhancers and co-transcriptional initiation over distances of nearly 250 kilobases. Regulatory interconnectivity depends on promoter-proximal tethering elements, and perturbations in these elements uncouple transcription and alter the bursting dynamics of distant genes, suggesting a role of genome topology in the formation and stability of co-transcriptional hubs. Transcriptional coupling is detected throughout the fly genome and encompasses a broad spectrum of conserved developmental processes, suggesting a general strategy for long-range integration of gene activity.

Cis-regulatory hubs: a new 3D model of complex disease genetics with an application to schizophrenia

Genes and their regulatory elements are organized in neurons within 3D networks which model functional structures and explain schizophrenia genetic etiology.

The 3D conformation of the chromatin creates complex networks of noncoding regulatory regions (distal elements) and promoters impacting gene regulation. Despite the importance of the role of noncoding regions in complex diseases, little is known about their interplay within regulatory hubs and implication in multigenic diseases such as schizophrenia. Here we show that cis-regulatory hubs (CRHs) in neurons highlight functional interactions between distal elements and promoters, providing a model to explain epigenetic mechanisms involved in complex diseases. CRHs represent a new 3D model, where distal elements interact to create a complex network of active genes. In a disease context, CRHs highlighted strong enrichments in schizophrenia-associated genes, schizophrenia-associated SNPs, and schizophrenia heritability compared with equivalent structures. Finally, CRHs exhibit larger proportions of genes differentially expressed in schizophrenia compared with promoter-distal element pairs or TADs. CRHs thus capture causal regulatory processes improving the understanding of complex disease etiology such as schizophrenia. These multiple lines of genetic and statistical evidence support CRHs as 3D models to study dysregulation of gene expression in complex diseases more generally.

Olfactory receptor choice: a case study for gene regulation in a multi-enhancer system

The mammalian genome possesses >2000 olfactory receptor (OR) alleles regulated by 63 known OR-Enhancer elements, yet each olfactory sensory neuron (OSN) expresses only a single OR allele. Choreographed changes to OSN nuclear architecture are evidently necessary for OR expression. Additionally, the insulated organization of OR-enhancers around an OR allele is a hallmark of the chosen OR. However, the biology guiding OR choice itself is unclear. Innovations in single-cell and biophysics-based analysis of nuclear architecture are revising previous models of the nucleus to include its dynamic and probabilistic nature. In this review, we ground current knowledge of OR gene regulation in these emerging theories to speculate on mechanisms that may give rise to diverse and singular OR expression.

Mouse strain-specific polymorphic provirus functions as cis-regulatory element leading to epigenomic and transcriptomic variations

Polymorphic integrations of endogenous retroviruses (ERVs) have been previously detected in mouse and human genomes. While most are inert, a subset can influence the activity of the host genes. However, the molecular mechanism underlying how such elements affect the epigenome and transcriptome and their roles in driving intra-specific variation remain unclear. Here, by utilizing wildtype murine embryonic stem cells (mESCs) derived from distinct genetic backgrounds, we discover a polymorphic MMERGLN (GLN) element capable of regulating H3K27ac enrichment and transcription of neighboring loci. We demonstrate that this polymorphic element can enhance the neighboring Klhdc4 gene expression in cis, which alters the activity of downstream stress response genes. These results suggest that the polymorphic ERV-derived cis-regulatory element contributes to differential phenotypes from stimuli between mouse strains. Moreover, we identify thousands of potential polymorphic ERVs in mESCs, a subset of which show an association between proviral activity and nearby chromatin states and transcription. Overall, our findings elucidate the mechanism of how polymorphic ERVs can shape the epigenome and transcriptional networks that give rise to phenotypic divergence between individuals.

An explainable artificial intelligence approach for decoding the enhancer histone modifications code and identification of novel enhancers in Drosophila

BACKGROUND

Enhancers are non-coding regions of the genome that control the activity of target genes. Recent efforts to identify active enhancers experimentally and in silico have proven effective. While these tools can predict the locations of enhancers with a high degree of accuracy, the mechanisms underpinning the activity of enhancers are often unclear.

RESULTS

Using machine learning (ML) and a rule-based explainable artificial intelligence (XAI) model, we demonstrate that we can predict the location of known enhancers in Drosophila with a high degree of accuracy. Most importantly, we use the rules of the XAI model to provide insight into the underlying combinatorial histone modifications code of enhancers. In addition, we identified a large set of putative enhancers that display the same epigenetic signature as enhancers identified experimentally. These putative enhancers are enriched in nascent transcription, divergent transcription and have 3D contacts with promoters of transcribed genes. However, they display only intermediary enrichment of mediator and cohesin complexes compared to previously characterised active enhancers. We also found that 10-15% of the predicted enhancers display similar characteristics to super enhancers observed in other species.

CONCLUSIONS

Here, we applied an explainable AI model to predict enhancers with high accuracy. Most importantly, we identified that different combinations of epigenetic marks characterise different groups of enhancers. Finally, we discovered a large set of putative enhancers which display similar characteristics with previously characterised active enhancers.

Transcription Factor Dynamics

To predict transcription, one needs a mechanistic understanding of how the numerous required transcription factors (TFs) explore the nuclear space to find their target genes, assemble, cooperate, and compete with one another. Advances in fluorescence microscopy have made it possible to visualize real-time TF dynamics in living cells, leading to two intriguing observations: first, most TFs contact chromatin only transiently; and second, TFs can assemble into clusters through their intrinsically disordered regions. These findings suggest that highly dynamic events and spatially structured nuclear microenvironments might play key roles in transcription regulation that are not yet fully understood. The emerging model is that while some promoters directly convert TF-binding events into on/off cycles of transcription, many others apply complex regulatory layers that ultimately lead to diverse phenotypic outputs. Cracking this kinetic code is an ongoing and challenging task that is made possible by combining innovative imaging approaches with biophysical models.

Molecular competition can shape enhancer activity in the Drosophila embryo

Transgenic reporters allow the measurement of regulatory DNA activity in vivo and consequently have long been useful tools for studying enhancers. Despite their utility, few studies have investigated the effects these reporters may have on the expression of other genes. Understanding these effects is required to accurately interpret reporter data and characterize gene regulatory mechanisms. By measuring the expression of Kruppel (Kr) enhancer reporters in live Drosophila embryos, we find reporters inhibit one another's expression and that of a nearby endogenous gene. Using synthetic transcription factor (TF) binding site arrays, we present evidence that competition for TFs is partially responsible for the observed transcriptional inhibition. We develop a simple thermodynamic model that predicts competition of the measured magnitude specifically when TF binding is restricted to distinct nuclear subregions. Our findings underline an unexpected role of the non-homogenous nature of the nucleus in regulating gene expression.

Variation in Pleiotropic Hub Gene Expression Is Associated with Interspecific Differences in Head Shape and Eye Size in Drosophila

Revealing the mechanisms underlying the breathtaking morphological diversity observed in nature is a major challenge in Biology. It has been established that recurrent mutations in hotspot genes cause the repeated evolution of morphological traits, such as body pigmentation or the gain and loss of structures. To date, however, it remains elusive whether hotspot genes contribute to natural variation in the size and shape of organs. As natural variation in head morphology is pervasive in Drosophila, we studied the molecular and developmental basis of differences in compound eye size and head shape in two closely related Drosophila species. We show differences in the progression of retinal differentiation between species and we applied comparative transcriptomics and chromatin accessibility data to identify the GATA transcription factor Pannier (Pnr) as central factor associated with these differences. Although the genetic manipulation of Pnr affected multiple aspects of dorsal head development, the effect of natural variation is restricted to a subset of the phenotypic space. We present data suggesting that this developmental constraint is caused by the coevolution of expression of pnr and its cofactor u-shaped (ush). We propose that natural variation in expression or function of highly connected developmental regulators with pleiotropic functions is a major driver for morphological evolution and we discuss implications on gene regulatory network evolution. In comparison to previous findings, our data strongly suggest that evolutionary hotspots are not the only contributors to the repeated evolution of eye size and head shape in Drosophila.

Decoding the organization, dynamics, and function of the 4D genome

Understanding how complex cell-fate decisions emerge at the molecular level is a key challenge in developmental biology. Despite remarkable progress in decoding the contribution of the linear epigenome, how spatial genome architecture functionally informs changes in gene expression remains unclear. In this review, we discuss recent insights in elucidating the molecular landscape of genome folding, emphasizing the multilayered nature of the 3D genome, its importance for gene regulation, and its spatiotemporal dynamics. Finally, we discuss how these new concepts and emergent technologies will enable us to address some of the outstanding questions in development and disease.

Enhancer redundancy in development and disease

Shadow enhancers are seemingly redundant transcriptional cis-regulatory elements that regulate the same gene and drive overlapping expression patterns. Recent studies have shown that shadow enhancers are remarkably abundant and control most developmental gene expression in both invertebrates and vertebrates, including mammals. Shadow enhancers might provide an important mechanism for buffering gene expression against mutations in non-coding regulatory regions of genes implicated in human disease. Technological advances in genome editing and live imaging have shed light on how shadow enhancers establish precise gene expression patterns and confer phenotypic robustness. Shadow enhancers can interact in complex ways and may also help to drive the formation of transcriptional hubs within the nucleus. Despite their apparent redundancy, the prevalence and evolutionary conservation of shadow enhancers underscore their key role in emerging metazoan gene regulatory networks.

Protecting P-type plasma membrane H+-ATPases from ROS

P-type ATPase are ubiquitous transport proteins across all kingdoms of life. These proteins share a common mechanism involving phosphorylation of an invariant aspartate to facilitate movement of substrates from protons to phospholipids across cellular membranes. In this issue of the Biochemical Journal, Welle et al. identify a conserved cysteine near the functionally critical aspartate of P-type plasma membrane H+-ATPases that protects the protein from reactive oxygen species.

Cis-regulatory chromatin loops arise before TADs and gene activation, and are independent of cell fate during early Drosophila development

Acquisition of cell fate is thought to rely on the specific interaction of remote cis-regulatory modules (CRMs), for example, enhancers and target promoters. However, the precise interplay between chromatin structure and gene expression is still unclear, particularly within multicellular developing organisms. In the present study, we employ Hi-M, a single-cell spatial genomics approach, to detect CRM-promoter looping interactions within topologically associating domains (TADs) during early Drosophila development. By comparing cis-regulatory loops in alternate cell types, we show that physical proximity does not necessarily instruct transcriptional states. Moreover, multi-way analyses reveal that multiple CRMs spatially coalesce to form hubs. Loops and CRM hubs are established early during development, before the emergence of TADs. Moreover, CRM hubs are formed, in part, via the action of the pioneer transcription factor Zelda and precede transcriptional activation. Our approach provides insight into the role of CRM-promoter interactions in defining transcriptional states, as well as distinct cell types.

Fra-1 regulates its target genes via binding to remote enhancers without exerting major control on chromatin architecture in triple negative breast cancers

The ubiquitous family of dimeric transcription factors AP-1 is made up of Fos and Jun family proteins. It has long been thought to operate principally at gene promoters and how it controls transcription is still ill-understood. The Fos family protein Fra-1 is overexpressed in triple negative breast cancers (TNBCs) where it contributes to tumor aggressiveness. To address its transcriptional actions in TNBCs, we combined transcriptomics, ChIP-seqs, machine learning and NG Capture-C. Additionally, we studied its Fos family kin Fra-2 also expressed in TNBCs, albeit much less. Consistently with their pleiotropic effects, Fra-1 and Fra-2 up- and downregulate individually, together or redundantly many genes associated with a wide range of biological processes. Target gene regulation is principally due to binding of Fra-1 and Fra-2 at regulatory elements located distantly from cognate promoters where Fra-1 modulates the recruitment of the transcriptional co-regulator p300/CBP and where differences in AP-1 variant motif recognition can underlie preferential Fra-1- or Fra-2 bindings. Our work also shows no major role for Fra-1 in chromatin architecture control at target gene loci, but suggests collaboration between Fra-1-bound and -unbound enhancers within chromatin hubs sometimes including promoters for other Fra-1-regulated genes. Our work impacts our view of AP-1.

Enhancer grammar in development, evolution, and disease: dependencies and interplay

Each language has standard books describing that language's grammatical rules. Biologists have searched for similar, albeit more complex, principles relating enhancer sequence to gene expression. Here, we review the literature on enhancer grammar. We introduce dependency grammar, a model where enhancers encode information based on dependencies between enhancer features shaped by mechanistic, evolutionary, and biological constraints. Classifying enhancers based on the types of dependencies may identify unifying principles relating enhancer sequence to gene expression. Such rules would allow us to read the instructions for development within genomes and pinpoint causal enhancer variants underlying disease and evolutionary changes.

Precision of tissue patterning is controlled by dynamical properties of gene regulatory networks

During development, gene regulatory networks allocate cell fates by partitioning tissues into spatially organised domains of gene expression. How the sharp boundaries that delineate these gene expression patterns arise, despite the stochasticity associated with gene regulation, is poorly understood. We show, in the vertebrate neural tube, using perturbations of coding and regulatory regions, that the structure of the regulatory network contributes to boundary precision. This is achieved, not by reducing noise in individual genes, but by the configuration of the network modulating the ability of stochastic fluctuations to initiate gene expression changes. We use a computational screen to identify network properties that influence boundary precision, revealing two dynamical mechanisms by which small gene circuits attenuate the effect of noise in order to increase patterning precision. These results highlight design principles of gene regulatory networks that produce precise patterns of gene expression.

To loop or not to loop: what is the role of TADs in enhancer function and gene regulation?

The past decade has seen a huge jump in the resolution and scale at which we can interrogate the three-dimensional properties of the genome. This revealed different types of chromatin structures including topologically associating domains, partitioning genes and their enhancers into interacting domains. While the visualisation of these topologies and their dynamics has dramatically improved, our understanding of their underlying mechanisms and functional roles in gene expression has lagged behind. A suite of recent studies have addressed this using genetic manipulations to perturb topological features and loops at different scales. Here we assess the new biological insights gained on the functional relationship between genome topology and gene expression, with a particular focus on enhancer function.

Single-molecule tracking of transcription protein dynamics in living cells: seeing is believing, but what are we seeing?

A universe of transcription factors (TFs), cofactors, as well as chromatin remodeling and modifying enzymes combine or compete on chromatin to control transcription. Measuring quantitatively how these proteins dynamically interact is required in order to formulate models with predictive ability to elucidate transcription control mechanisms. Single molecule tracking (SMT) provides a powerful tool towards this goal: it is a fluorescence microscopy approach that measures the location and mobility of individual TF molecules, as well as their rates of association with and dissociation from chromatin in the physiological context of the living cell. Here we review SMT principles, and discuss key TF properties uncovered by live-cell SMT, such as fast turnover (seconds), and formation of clusters that locally increase activity.

Saccharomyces cerevisiae rDNA as super-hub: the region where replication, transcription and recombination meet

Saccharomyces cerevisiae ribosomal DNA, the repeated region where rRNAs are synthesized by about 150 encoding units, hosts all the protein machineries responsible for the main DNA transactions such as replication, transcription and recombination. This and its repetitive nature make rDNA a unique and complex genetic locus compared to any other. All the different molecular machineries acting in this locus need to be accurately and finely controlled and coordinated and for this reason rDNA is one of the most impressive examples of highly complex molecular regulated loci. The region in which the large molecular complexes involved in rDNA activity and/or regulation are recruited is extremely small: that is, the 2.5 kb long intergenic spacer, interrupting each 35S RNA coding unit from the next. All S. cerevisiae RNA polymerases (I, II and III) transcribing the different genetic rDNA elements are recruited here; a sequence responsible for each rDNA unit replication, which needs its molecular apparatus, also localizes here; moreover, it is noteworthy that the rDNA replication proceeds almost unidirectionally because each replication fork is stopped in the so-called replication fork barrier. These localized fork blocking events induce, with a given frequency, the homologous recombination process by which cells maintain a high identity among the rDNA repeated units. Here, we describe the different processes involving the rDNA locus, how they influence each other and how these mutual interferences are highly regulated and coordinated. We propose that an rDNA conformation as a super-hub could help in optimizing the micro-environment for all basic DNA transactions.

Navigating the crowd: visualizing coordination between genome dynamics, structure, and transcription

The eukaryotic genome is hierarchically structured yet highly dynamic. Regulating transcription in this environment demands a high level of coordination to permit many proteins to interact with chromatin fiber at appropriate sites in a timely manner. We describe how recent advances in quantitative imaging techniques overcome caveats of sequencing-based methods (Hi-C and related) by enabling direct visualization of transcription factors and chromatin at high resolution, from single genes to the whole nucleus. We discuss the contribution of fluorescence imaging to deciphering the principles underlying this coordination within the crowded nuclear space in living cells and discuss challenges ahead.

Dense and pleiotropic regulatory information in a developmental enhancer

Changes in gene regulation underlie much of phenotypic evolution¹. However, our understanding of the potential for regulatory evolution is biased, because most evidence comes from either natural variation or limited experimental perturbations². Using an automated robotics pipeline, we surveyed an unbiased mutation library for a developmental enhancer in Drosophila melanogaster. We found that almost all mutations altered gene expression and that parameters of gene expression-levels, location, and state-were convolved. The widespread pleiotropic effects of most mutations may constrain the evolvability of developmental enhancers. Consistent with these observations, comparisons of diverse Drosophila larvae revealed apparent biases in the phenotypes influenced by the enhancer. Developmental enhancers may encode a higher density of regulatory information than has been appreciated previously, imposing constraints on regulatory evolution.

Robust and efficient gene regulation through localized nuclear microenvironments

Developmental enhancers drive gene expression in specific cell types during animal development. They integrate signals from many different sources mediated through the binding of transcription factors, producing specific responses in gene expression. Transcription factors often bind low-affinity sequences for only short durations. How brief, low-affinity interactions drive efficient transcription and robust gene expression is a central question in developmental biology. Localized high concentrations of transcription factors have been suggested as a possible mechanism by which to use these enhancer sites effectively. Here, we discuss the evidence for such transcriptional microenvironments, mechanisms for their formation and the biological consequences of such sub-nuclear compartmentalization for developmental decisions and evolution.

Seven myths of how transcription factors read the cis-regulatory code

Genomics data are now being generated at large quantities, of exquisite high resolution and from single cells. They offer a unique opportunity to develop powerful machine learning algorithms, including neural networks, to uncover the rules of the cis-regulatory code. However, current modeling assumptions are often not based on state-of-the-art knowledge of the cis-regulatory code from transcription, developmental genetics, imaging and structural studies. Here I aim to fill this gap by giving a brief historical overview of the field, describing common misconceptions and providing knowledge that might help to guide computational approaches. I will describe the principles and mechanisms involved in the combinatorial requirement of transcription factor binding motifs for enhancer activity, including the role of chromatin accessibility, repressors and low-affinity motifs in the cis-regulatory code. Deciphering the cis-regulatory code would unlock an enormous amount of regulatory information in the genome and would allow us to locate cis-regulatory genetic variants involved in development and disease.

Transcriptional Enhancers in Drosophila

Key discoveries in Drosophila have shaped our understanding of cellular "enhancers." With a special focus on the fly, this chapter surveys properties of these adaptable cis-regulatory elements, whose actions are critical for the complex spatial/temporal transcriptional regulation of gene expression in metazoa. The powerful combination of genetics, molecular biology, and genomics available in Drosophila has provided an arena in which the developmental role of enhancers can be explored. Enhancers are characterized by diverse low- or high-throughput assays, which are challenging to interpret, as not all of these methods of identifying enhancers produce concordant results. As a model metazoan, the fly offers important advantages to comprehensive analysis of the central functions that enhancers play in gene expression, and their critical role in mediating the production of phenotypes from genotype and environmental inputs. A major challenge moving forward will be obtaining a quantitative understanding of how these cis-regulatory elements operate in development and disease.

CDK-Mediator and FBXL19 prime developmental genes for activation by promoting atypical regulatory interactions

Appropriate developmental gene regulation relies on the capacity of gene promoters to integrate inputs from distal regulatory elements, yet how this is achieved remains poorly understood. In embryonic stem cells (ESCs), a subset of silent developmental gene promoters are primed for activation by FBXL19, a CpG island binding protein, through its capacity to recruit CDK-Mediator. How mechanistically these proteins function together to prime genes for activation during differentiation is unknown. Here we discover that in mouse ESCs FBXL19 and CDK-Mediator support long-range interactions between silent gene promoters that rely on FBXL19 for their induction during differentiation and gene regulatory elements. During gene induction, these distal regulatory elements behave in an atypical manner, in that the majority do not acquire histone H3 lysine 27 acetylation and no longer interact with their target gene promoter following gene activation. Despite these atypical features, we demonstrate by targeted deletions that these distal elements are required for appropriate gene induction during differentiation. Together these discoveries demonstrate that CpG-island associated gene promoters can prime genes for activation by communicating with atypical distal gene regulatory elements to achieve appropriate gene expression.

Imaging of DNA and RNA in Living Eukaryotic Cells to Reveal Spatiotemporal Dynamics of Gene Expression

This review focuses on imaging DNA and single RNA molecules in living cells to define eukaryotic functional organization and dynamic processes. The latest advances in technologies to visualize individual DNA loci and RNAs in real time are discussed. Single-molecule fluorescence microscopy provides the spatial and temporal resolution to reveal mechanisms regulating fundamental cell functions. Novel insights into the regulation of nuclear architecture, transcription, posttranscriptional RNA processing, and RNA localization provided by multicolor fluorescence microscopy are reviewed. A perspective on the future use of live imaging technologies and overcoming their current limitations is provided.