REVIEW article

Front. Immunol., 04 June 2021

Sec. B Cell Biology

Volume 12 - 2021 | https://doi.org/10.3389/fimmu.2021.682397

The Causes and Consequences of Spatial Organization of the Genome in Regulation of Gene Expression

  • Biomedical Research Foundation, Academy of Athens, Athens, Greece

Abstract

Regulation of gene expression in time, space and quantity is orchestrated by the functional interplay of cis-acting elements and trans-acting factors. Our current view postulates that transcription factors recognize enhancer DNA and read the transcriptional regulatory code by cooperative DNA binding to specific DNA motifs, thus instructing the recruitment of transcriptional regulatory complexes forming a plethora of higher-ordered multi-protein-DNA and protein-protein complexes. Here, we reviewed the formation of multi-dimensional chromatin assemblies implicated in gene expression with emphasis on the regulatory role of enhancer hubs as coordinators of stochastic gene expression. Enhancer hubs contain many interacting regulatory elements and represent a remarkably dynamic and heterogeneous network of multivalent interactions. A functional consequence of such complex interaction networks could be that individual enhancers function synergistically to ensure coordination, tight control and robustness in regulation of expression of spatially connected genes. In this review, we discuss fundamental paradigms of such inter- and intra- chromosomal associations both in the context of immune-related genes and beyond.

What Are the Enhancers

Regulation of gene expression in time, space and quantity is orchestrated by the functional interplay of cis-acting regulatory elements and trans-acting factors. The establishment of cell identity during development and the response to environmental cues such as pathogens requires transcriptional reprogramming (1, 2). Cells respond to signals by engaging various receptors that recognize and interpret the environment by initiating signal transduction cascades, culminating in the activation of transcription factors (TFs), which bind to specific regulatory DNA sequences to control the expression of genes, thus orchestrating a dynamic interplay between genome form and function (36). Inappropriate regulation of gene expression during development or adult life can cause abnormal embryogenesis and disease (7, 8). Enhancers and promoters represent the major classes of DNA regulatory elements responsible for warranting the execution of precise gene expression programs by functioning as information processing units interpreting the extra- and intra-cellular signals in the form of transcription factor binding events, which are followed by the recruitment of various transcriptional regulatory proteins such as cofactors, chromatin modifiers and the transcriptional apparatus at the core promoter (915). Core promoters are DNA regions, surrounding the start site of transcription (TSS), directing the accurate initiation of transcription by interacting with the various components of the general transcription machinery, including RNA Pol II (16). In addition, the DNA region (100-200 bp) located immediately upstream of the TSS often bears critical transcription factor binding sites and it is required to increase the rate of transcription (17, 18). In some cases, upstream promoter elements can function in a cell-type specific manner to mediate the interactions with distal regulatory elements like the enhancers (19, 20).

Enhancers activate transcription of target genes independent of their distance, position and orientation relative to the promoter of the affected gene (10, 21, 22). In certain cases, enhancers can act synergistically to ensure for proper gene activation (23). Our current model envisions that TFs recognize enhancer DNA in a combinatorial manner by reading the transcriptional regulatory code through the interplay of nucleotide sequence recognition, sequence context, 3D DNA structure (shape), DNA binding cooperativity, interacting cofactors and nucleosome occupancy (2429). Establishment and maintenance of cell identity and function arise from integrated actions of transcription factors with chromatin–associated proteins, coactivators, and with the spatially organized genome (30). Combinatorial TF DNA binding leads to a synchronized targeting of enhancers by multifunctional transcriptional regulatory proteins such as coactivators, corepressors, chromatin modifiers and remodelers through direct and indirect physical associations, to assemble a plethora of high-order multi-protein-DNA and protein-protein complexes required for regulation of transcription (4, 9, 14). TFs belong to a diverse class of DNA binding protein families with distinct structural characteristics and individual members from the same family can bind similar DNA motifs, yet carry out distinct functions, in vivo (3133). In addition, as the number of DNA binding sites in the genome for most TFs is greater than the number of the corresponding protein molecules per cell, it appears that only some TF molecules occupy most of their functional relevant sites in physiological conditions (3438). TFs typically recognize and bind to short 6-12bp-long degenerate DNA motifs. Consistently, TFs bind to thousands of genomic sites across chromatin landscapes, but the number of genes changing their expression is frequently limited to a few hundred (3942). Cells, at least in part, resolve this apparent paradox through combinatorial strategies where TFs must work together when bound to enhancers and promoters to achieve the specificity required for the selection of the appropriate genes to become activated. Cooperative DNA binding of TFs to low-affinity sites provides the means to target the correct genes for expression in time, space and in response to signals and is mediated by direct or indirect DNA-and/or chromatin-facilitating contacts between the proteins (33). Interestingly, high throughput ChIP-seq studies identified that binding of TFs to certain genomic regions is highly clustered, each cluster is composed of both functionally related and unrelated TFs, representing strong binding of one or a few TFs, and weaker binding to more degenerate motifs by many other TFs, implying a substantial degree of cooperation between individual motifs (43, 44). Large genomic regions with an unusually high degree of enrichment for lineage-specific TFs, cofactors (e.g. Med1) and histone modifications (H3K27ac) have been recently identified and characterized as super-enhancers (SEs), different from classical enhancers. These clusters of closely spaced enhancer elements often regulate the expression of genes determining cell identity and fate and driving oncogenic transcription (4548). In general, super-enhancers are involved in strong intradomain interactions that is, between its individual enhancer constituents, which are stronger compared to typical enhancer-promoter interactions (15, 49). As super-enhancers consist of multiple enhancer elements, we assume that there will be complex functional relationships among the different constituents of any given super-enhancer working, depending on the biological context, either in a synergistic or additive manner or a combination thereof, thus resembling the arrangement and functional relationships of TF motifs building a typical enhancer.

One of the best-characterized examples of assembled super-enhancers is provided after stimulation of endothelial cells with TNFα, which leads to the formation of super-enhancers specific for inflammatory genes (50). Major players for the inducible assembly of these super-enhancers are the transcription factor NF-κB and the coactivator BRD4, a member of the bromodomain and extraterminal domain (BET) family of factors that is typically recruited to enhancers and super-enhancers. In brief, stimulation of endothelial cells with TNFα promotes the rapid redistribution of NF-κB across diverse genomic loci and subsequently the assembly of inflammatory SE chromatin configurations. This is accompanied by striking colocalization of the NF-κB p65 subunit and BRD4 in regulatory sequences marked by histone acetylation (H3K27ac). Thus, NF-κB efficiently targets its binding sites and promotes a pronounced recruitment of BRD4 across inflammatory SEs, which coordinate the development of a pro-inflammatory gene expression program (50).

The blood cell-specific SE cluster, known as BENC, is located 1.7 MB downstream of the Myc locus, and consists of multiple enhancer modules working in a combinatorial manner. BENC is bound by the transcription factors GFI1b, RUNX1 and MYB (51) and is essential for the regulation of Myc expression in hematopoietic and leukemic stem cells critical for hematopoietic malignancies. BENC deletion in hematopoietic stem cells closely mimics the conditional deletion of the Myc gene. Thus, overlapping, adjacent and/or partially redundant regulatory units arranged in clusters and bound by a variety of transcription factors functioning in a combinatorial manner, are essential for the generation of cell-type specific transcriptional regulatory outputs.

Enhancer-Promoter Networks

The application of high-throughput sequencing-based chromatin profiling (e.g. ChIP-seq, DNAseI-seq), ideal for measuring the potential regulatory activity, indicates that the human genome contains about a million candidate different enhancers (52, 53) interspersed at regions located proximally, distantly, or within genes, and that in each cell type a fraction of these putative enhancers is active. Enhancers regulate gene expression in quite diverse manners. For example, enhancers do not necessarily regulate their closest gene (10), each gene may be regulated by more than one enhancer (54, 55), and a single enhancer can regulate multiple genes (56). The above mechanistic principles suggest the existence of complex and highly dynamic enhancer-promoter networks and diverse modes of enhancer action. The multiplicity of the enhancer-promoter functional interactions is highly regulated by transcription factors and cofactors associated with them to direct the re-configuration and the spatial folding of the 3D chromatin interaction, thus bridging remotely located transcriptional regulatory elements to a functional assembly. Although there is strong evidence that enhancer-promoter interactions are usually established concurrently with gene activation, and thus comprising an integral part of global genomic regulation, it is not yet proven whether enhancer-promoter contacts are the cause or the consequence of regulated gene expression. Some enhancer-promoter paired configurations are formed in the absence of any transcriptional events, thus suggesting that, simply, their identification does not necessarily have functional implications (57). However, forced chromatin looping between an enhancer and a promoter led to induction of transcription at high levels (58), thus strongly suggesting the requirement for establishing functional 3D structures connecting enhancers and promoters that could be integrated into a higher order condensed multi-genic 3D transcriptionally-active hubs at specific nuclear positions to coordinate gene expression of functionally related genes (59, 60).

Furthermore, although the identification of specific enhancer-promoter contacts provides important information regarding the assignment of enhancers to the correct target genes, these interactions are not necessarily predictive of functional regulation. For example, recent studies have shown that sonic hedgehog (Shh) expression in the developing embryo is regulated by multiple enhancer elements. Using a combination of 3D-DNA FISH and chromosome conformation capture approaches, Benabdallah et al., demonstrated that Shh expression is incompatible with the classic enhancer-promoter looping model. The authors found that, practically, there is no spatial proximity between the enhancer and the promoter during the differentiation of embryonic stem cells to progenitor neural cells (61). In contrast, however, Shh expression in the zone of polarized activity in the developing limb in mice, which is controlled by a different enhancer called ZRS, is characterized by the highest enhancer-promoter proximity as compared to other non-expressing tissues, a finding consistent with the formation of a tight chromosomal loop between the enhancer and the promoter (62). Similar to the findings of Shh expression in ES cells, live-cell imaging experiments revealed that the expression of the key pluripotency transcription factor Sox2 in ES cells is not controlled by the spatial proximity between its control region SCR and the promoter, despite the fact that 3C assays have identified enriched contacts between the enhancer and the promoter (63). Although it is not yet fully understood how a spatially remote enhancer could affect the activity of its target promoter, we speculate that the assembly of condensates composed of transcriptional regulatory proteins could trigger the formation of large macromolecular bridges or hubs, which may alter the configuration of the intervening chromatin structure, thus leading to increased enhancer-promoter spatial distances, despite the fact that these regulatory sequences are functionally linked.

Essential questions related to the formation of dense enhancer-promoter or enhancer-enhancer networks are how they are established, maintained, or dissolved. We speculate that it would be difficult or nearly impossible for any distal enhancer to locate its target promoter if the genome was largely homogeneously structured, thus lacking an internal higher order organization. We now know that, in general, the hierarchical organization of chromatin involves the A (euchromatin) and B (heterochromatin) compartments, with compartment B marked by a more compact chromatin packaging as compared to compartment A (64), which is characterized by increased chromatin accessibility and higher gene expression levels (65). Chromatin contacts within each compartment are specifically enriched, but contacts between genomic regions present in the A and B compartments are not favored. Each chromatin compartment is further substructured to additional hierarchical layers termed Topological Associated Domains (TADs) (66, 67), where regions within the same TAD tend to interact more frequently with each other than with regions belonging to different TADs. Thus, within TADs, numerous genomic interactions (loops) are formed connecting regulatory elements such as enhancers and promoters, with a significantly smaller number of interactions spanning across the boundaries of adjacent TADs. TAD borders are frequently demarcated by the binding of CTCF, a TF functioning as a chromatin organizer (68). TAD anchors are frequently bound by the CTCF-cohesin complex usually in a convergent motif orientation (66) through a loop extrusion process mediated by the cohesin ring until the process is blocked at CTCF sites arranged in a convergent orientation (69). Mutations at the TAD boundaries can lead to chromatin reconfiguration at the TAD borders causing enhancers present in one TAD to activate transcription of genes present in the neighboring TAD (70). More specifically, rearrangement of TAD boundaries at the WNT6/IHH/EPHA4/PAX3 locus caused an extensive rewiring of the interactions between enhancers and promoters leading to limb malformations (70). The functional insulatory role of CTCF in gene regulation was further investigated in a study where the disruption of CTCF binding at the TAD boundaries via CRISPR led to the aberrant interaction of a constitutive enhancer with the PDGFRA oncogene in glioma (71). Interestingly, removal of CTCF binding sites both at the boundary and within TADs at the Sox9-Kcnj2 locus led to fusion of neighboring TADs, which, however, did not result in major effects on developmental gene regulation (72). In contrast, however, inversion of CTCF binding sites within the protocadherin enhancer rewires the enhancer-promoter interactions and alters gene expression patterns (73). Taken together, the data presented above indicate that the precise role of chromatin topology in enhancer/promoter function is determined by the specific genomic context.

An important issue regarding the multiplicity of enhancer-promoter network of interactions relates to the mechanisms by which an enhancer selects the correct promoter for activation. While within TADs, the number of candidate promoter targets for a given enhancer is relatively limited, in many cases enhancers and their target promoters are not adjacent in the primary DNA sequence, but interrupted by one or more genes that are not co-expressed or co-regulated with the enhancer’s target gene. For example, expression of the Sex comps reduced (Scr) gene during Drosophila embryogenesis is controlled by an enhancer which is located far upstream. Importantly, the continuity of the Scr gene and its enhancer is interrupted by the localization of the fushi tarazu (ftz) gene, which is regulated by an enhancer positioned at close distance from the Scr gene, whereas the Scr enhancer is present at the 3’ of the ftz gene. In other words, the Scr enhancer is closer to ftz than to Scr promoter, and the ftz enhancer is flanked by the Scr and ftz promoters. However, both genes are expressed at different times during embryogenesis, and their expression is controlled by their individual enhancers. The selectivity of the Scr enhancer for activating transcription from the Scr promoter and not from the closest ftzp promoter depends on promoter proximal DNA elements called tethering elements, which when bound by specific transcription factors might mediate highly specific enhancer-promoter interactions between the Scr enhancer and promoter (19, 74). These data are compatible with the looping and relocation models for enhancer-promoter communication, because they can explain how an enhancer can skip more proximal promoters (75).

The structural chromatin loops responsible for forming TADs are dependent on ubiquitously expressed CTCF and cohesin proteins and therefore are distinct from those assembled by the interactions of regulatory sequences. The latter are generally mediated by enhancer- and promoter-bound transcription factors and cofactors and therefore their formation is dynamic and developmentally regulated (30). Our current understanding is that TADs define regions of the genome where enhancer-promoter interactions are allowed to occur to ensure robust and reliable transcriptional activation, thus precluding inappropriate regulatory interactions between enhancers and promoters located in different TADs (76). Despite the established roles of the architectural protein CTCF in supporting regulatory enhancer-promoter interactions (66, 73), its depletion although leads to a significant reduction of TAD structures results in moderate rather than dramatic changes in global gene expression (77, 78). An emerging idea derived from single-cell 3C and advanced imaging experiments (79, 80) to explain this unexpected finding is that the cell-to-cell variability in the overall TAD structure indicates that the existence of a highly heterogeneous structure of the chromatin could provide a reasonable explanation for the relatively modest effects in transcription upon TAD disruption. Thus, it is not yet fully understood how TADs and TAD boundaries affect the regulation of gene expression, and the regulatory logic that rules their involvement.

Understanding the Regulatory Genome in 3D

The spatial arrangement of the regulatory genome in 3D has been mainly studied by Chromatin Conformation Capture (3C) methodologies and its ChIP or probe-based derivatives and advanced high-resolution microscopy techniques (62, 66, 8185). Notably, some of the first examples of well-characterized mechanisms of gene regulation by DNA-DNA contacts in mammals, were derived from studies in hematopoietic and immune cells regarding the 3D organization of the β-globin gene locus (86) and the interchromosomal association between the IFN-gamma gene locus and regulatory regions of the T(H)2 cytokine locus in a dynamic and cell-type-specific fashion (87). Additional studies have provided evidence that transcription factors bound to enhancers and promoters mediate this type of DNA-DNA communication by acting as the components exposing complementary interacting protein surfaces required for the establishment and the specificity for the formation of looped structures (20). A well-characterized example of a looping factor essential for chromatin interactions at the β-globin locus is the cofactor Ldb1 (58). Ldb1 together with GATA1 acts as a molecular bridge, which facilitates the communication between the Locus Control Region (LCR) and the β-globin promoter.

In the case of the virus inducible human IFNβ stochastic gene expression, the IFNβ enhancer contributes to the formation of highly specialized 3D associations required for singular and stochastic expression of the IFNβ gene in a subset of infected cells (35, 88, 89). Previous work has indicated that virus infection induces the coordinated activation of three distinct sets of transcription factors, NF-κB, IRF3 and 7, and ATF2/cJun, which together with the architectural protein HMG I(Y) bind cooperatively to the nucleosome-free enhancer/promoter of the IFNβ gene and form a nucleoprotein structure known as the enhanceosome (24, 9095) (Figure 1). Following its assembly, the IFNβ enhanceosome exposes a continuous recognition surface used to recruit the transcriptional coactivator CBP (93). Enhanceosomes with similar properties have now been identified and characterized on several model enhancers in both Drosophila and mammalian systems. For example, cell type-specific and/or inducible expression of the MHC class I and MHC class II expression, which is crucial for the initiation and regulation of adaptive immunity, requires the assembly of enhanceosomes which activate transcription by recruiting the coactivators NLRC5- and CIITA, respectively (96). As it is the case for the IFNβ enhanceosome, the NLRC5 and CIITA coactivators are recruited to their target promoters by interacting simultaneously with a continuous surface of the MHC Class I and Class II enhanceosome components, respectively.

Figure 1

The human IFNβ promoter contains one well-positioned nucleosome physically obscuring the TATA box and the start site of transcription. Once assembled, the IFNβ enhanceosome dynamically recruits histone acetyltransferases and chromatin remodelers causing sliding of the nucleosome masking the TATA to a new position 36bp downstream, thus unmasking the TATA box and the start site of transcription (91, 97, 98). The gene is initially expressed from a single allele, but its expression is subsequently becoming biallelic (88, 99). The discovery of the early monoallelic IFNβ expression suggested that gene activation itself is successfully completed only with a certain probability in a given cell. A similar mode of monoallelic stochastic activation has also been described for the IFNα family of genes and it is a typical and general feature of cytokine expression during immune responses (100, 101).

Given the pivotal role of IFNs in antiviral host responses, it has come as a surprise that, at the single-cell level, IFNs are expressed only by a fraction of virus-infected cells. Indeed, approximately ~20% of the virus-infected cells are expressing IFNβ. Moreover, the induction of Interferon Stimulated Genes (ISGs) after acute viral infections has also been found to be highly heterogeneous at the single-cell level (99), with only a fraction of cells inducing an antiviral gene expression program at physiological concentrations of IFNs.

Recent findings have indicated that IFNβ expression heterogeneity is due, at least in part, to stochastic interchromosomal interactions driving the variable patterns of IFNβ gene expression in response to virus infection (35, 88, 89). The important finding was that stochasticity in expression of human IFNβ is due to low intra-cellular concentrations of the transcription factor NF-κB, which is captured onto specialized DNA sites called NRCs (NF-κB Reception Centers) belonging to Alu-like repetitive elements and then it is delivered to the IFNβ enhancer, via stochastic interchromosomal interactions. It was shown that the transcription factor ThPOK (89), a GAGA binding factor, binds cooperatively with NF-κB to NRCs and mediates their physical proximity with the IFNβ gene and other co-regulated genes via its ability to oligomerize when bound to DNA, as it was shown by in vitro experiments. Furthermore, ThPOK knockdown experiments significantly decreased the frequency of interchromosomal interactions, NF-κB DNA binding to the IFNβ enhancer, and virus-induced IFNβ gene activation. Cooperative DNA binding between ThPOK and NF-κB, in vitro and in vivo, on the same face of the double DNA helix, is required for interchromosomal interactions and this distinguishes NRCs from various other Alu elements bearing interspersed NF-κB and ThPOK sites (89) (Figure 1A). These studies showed how DNA binding cooperativity of stereospecifically aligned transcription factors provides the necessary ultrasensitivity for the all-or-none stochastic cell responses to virus infection. The above mechanistic insights strongly suggest that the primary DNA sequence of the interacting DNA elements acts as a blueprint, which contains a specific linear genomic code that is interpreted via the construction of highly specialized high-order nucleoprotein complexes. Additional studies have also shown that stochastic IFNβ expression in response to virus infection is due to cell-to-cell variability due to limiting quantities of components ranging from the recognition of viral RNA by host factors and the activation of signaling pathways, to the exact levels of activated transcription factors (102).

Additional examples of architectural proteins acting as organizers of 3D genome architecture include YY1 (103), ZNF143 (104), and SATB1 (105). In contrast to CTCF, YY1 preferentially binds to active enhancers and promoters mediating their 3D interaction via its ability to dimerize (103). ZNF143 binds to the anchors of chromatin interactions and facilitates the formation of secure connections between promoters and distal regulatory elements (104). Finally, SATB1 can act as a tissue-specific organizer of gene expression with well-established roles in the regulation of TH2 cytokine gene locus (105, 106) by cooperative binding to multiple closely spaced consensus sites in nucleosomes (107). In summary, the examples described above demonstrate that the formation of 3D DNA-DNA nuclear contacts is not a purely random phenomenon but it depends on sequences hosting chromatin architectural proteins, which in many cases are in direct interaction with TFs.

The molecular mechanisms controlling olfactory receptor (OR) gene transcription (108) inducing the monoallelic expression of one out of more than 1,000 genes in each cell, share many common features with the cellular responses leading to virus-induced stochastic monoallelic IFNβ expression. The singular and stochastic expression of OR genes is controlled by a dense network of 3D DNA interactions occurring exclusively in olfactory sensory neurons (109111). The OR genes, which are interspersed in various chromosomes, aggregate in one nuclear compartment, where they are repressed by acquiring chromatin marks characteristic of heterochromatin. However, this aggregation brings in close proximity many OR enhancers located adjacent to the OR genes via multiple interchromosomal interactions forming a multi-enhancer hub. Monoallelic OR transcription is activated by the multi-enhancer hub, which associates stochastically with the corresponding gene chosen for expression. The TFs Lhx2 and Ebf bind cooperatively in a stereospecific manner on OR enhancers followed by the recruitment of the LIM domain protein Ldb1, which facilitates these interchromosomal interactions. Thus, DNA binding cooperativity between transcription factors, as it is also the case for IFNβ expression, provides the necessary specificity to drive spatiotemporal alterations of the 3D chromatin for the cell-type specific expression of OR genes and cellular responses upon viral infections, respectively.

Many mechanistic similarities to the stochastic monoallelic OR and IFNβ expression are also shared by the NF-κB-dependent TNFα-induced hierarchical co-expression of SAMD4A, TNFAIP2, and SLC6A5 genes (112, 113), which also form a tight signal-dependent multigene complex through intra- and interchromosomal interactions. TNFα induction causes the hierarchical co-transcription of these genes in a small fraction of cells containing tightly associated genes in a multigene complex.

Photoreceptor specification in the Drosophila eye depends on the stochastic and cell-autonomous expression of the transcription factor Spineless (Ss) (114116). The stochastic Ss expression in individual cells is regulated by a 60 kb region harboring the R7/8 enhancer and the silencer elements 1 and 2 (115, 116). It was shown that the stochastic Ss expression pattern depends on interchromosomal communication and on inter-allelic crosstalk, which averages the frequency of expression of each allele (116). Thus, although each allele follows an independent-stochastic choice of expression, 3D interchromosomal communication coordinates the expression status between alleles, thus ensuring for the proper expression of both alleles within the same randomly selected subset of photoreceptor cells. Conclusively, both the 3D interchromosomal interactions and the action of enhancers from a distance determine stochastic photoreceptor specification.

In summary, the examples described above when taken together with many other similar studies, suggest a general model according to which cell activation by environmental signals or during development and differentiation leads to a cascade of multiprotein and multigene high order dynamic assemblies, ranging from signalsomes to enhanceosomes and chromatin-modifying machines, thus providing local regulatory input from different genomic elements. This mode of information flow ensures that stochastic molecular interactions are canalized to produce specific, robust and reliable outputs. Furthermore, these observations strongly support the view that transcription factors are capable of inducing extensive chromatin rewiring, in many cases stochastic, at specific genomic loci, thus modifying gene expression patterns in a cell-specific manner.

Enhancer Hubs-Transcriptional Condensates

Enhancer hubs containing highly interacting regulatory elements represent a remarkably dynamic and heterogeneous network of multivalent macromolecular interactions. Most of these interactions could be established in a stochastic manner and individual cells within a homogeneous population may contain unique combinations of two or more interacting partners. The higher the number of interacting partners the lower the probability of a cell in the population to contain multi-partner hubs. We speculate that enhancer hubs by containing multiple enhancers targeting a single promoter or many promoters targeted by a single enhancer or many interconnecting enhancers targeting many different promoters, could aggregate into a single nuclear position. A functional consequence of such complex interaction networks could be that different enhancers function synergistically to ensure coordination, tight control and robustness in regulation of expression of spatially connected genes. Thus, enhancer-promoter hubs ensure that all target genes are “on” and are appropriately regulated in a spatiotemporal manner, with stable relative levels of expression. Thus, each gene is subject to the same intrinsic and extrinsic noise e.g. fluctuations in the concentration of transcriptional regulators etc. The high density of transcription factors bound to enhancers and promoters in the context of enhancer hubs or super-enhancers leads to the cooperative and synergistic recruitment of high concentration of coactivators, components of the transcriptional machinery and other cofactors (117). Conclusively, an important implication derived from the assembly of enhancer hubs is the accumulation of transcriptional regulators at high local concentrations, which then are equally available for binding to the enhancer and promoter elements of the associated genes in a locus-specific manner.

Since many transcriptional regulators typically contain large low complexity and/or intrinsically disordered regions (IDRs), like the activation domains of many TFs, the probability to undergo hetero- or homo-typic multivalent interaction or oligomerization is high (118). IDRs have been described as functional protein segments that are not likely to form a defined 3D structure (119) and have been classified mainly by their amino acid profile or by their hypothetical shape (120). High densities of proteins bearing IDRs and nucleic acids have been implicated in the formation of membraneless compartments through the well-understood mechanism of liquid-liquid phase separation (LLPS) (121125). The term phase separation describes a process by which distinct liquid droplets bearing high concentration of molecules are formed. The formation of such nuclear compartments increases the likelihood for certain biochemical reactions to occur, such as transcription, thus bypassing partially or entirely the necessity for direct physical contacts between individual enhancers and promoters (Figure 1B). Theoretically, an additional advantage of regulating coordination of transcription via a phase separation model is the plasticity and reversibility by which LLPS formations are driven, and thus provide a reliable mechanism of transcriptional initiation and termination of many coregulated genes at once. For example, clusters of enhancer elements have been proposed to regulate biological processes implicated in development and disease (47), as well as regulating mammalian cell identity (45). We speculate that in the case of the SEs assembling upon TNFα induction of endothelial cells (see above), the pronounced recruitment of cofactors, like BRD4, etc. could lead to the additional recruitment of RNA Pol II and cofactors via an extensive network of intermolecular interactions between specific regions of intrinsically disordered hydrophilic activation domains (IDRs), thus leading to the formation of phase-separated nuclear liquid droplets (126) (Figure 1B).

The proposed model describes that the highly condensed concentration of biochemical interactions between TFs and coactivators, occurring on TFBSs-enriched extended 3D chromatin configurations, promote the colocalization/co-recruitment of diverse SEs at specific nuclear sites and subsequently the formation of 3D chromatin condensates (Figure 1B).

The ability of eukaryotic cells to assemble high-order transcription condensates ensures that enhancers bearing low-affinity binding sites can efficiently recruit RNA Pol II molecules to target genes to support robust patterns of gene expression (27, 125, 126). It is reasonable to assume that the individual functions of SEs could become modified upon their recruitment to the 3D chromatin condensates, where they could acquire novel regulatory characteristics hardwired in dense networks of macromolecular interactions.

Funding

DT was supported from the Greek Secretariat for Research and Technology (Cooperative Grants SYNERGASIA I #969, Excellence Award ARISTEIA I #1567, European Committee (FP7 projects INTEGER, NANOMA, PREDICTA and BIOFOS), from the European Economic Area (EL0084), BIO-IMAGING GR (MIS-5002755), and from the KMW offsets program. MA was supported from research grants from the Bodossaki Foundation, Fondation Santé, Gilead Hellas, and Antonios and Ioannis Angelicoussis Foundation. The funder was not involved in the study design, collection, analysis, interpretation of data, the writing of this article or the decision to submit it for publication.

Statements

Author contributions

All authors contributed to the article and approved the submitted version.

Acknowledgments

We thank Dimitris Valakos and Marianna Koutsi and members of our lab for helpful discussions. We apologize to the authors of many relevant and important studies for not citing their work due to space limitations.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  • 1

    SungMH. Transcriptional Reprogramming. In: Dubitzky W., Wolkenhauer O., Cho KH., Yokota H. (eds) Encyclopedia of Systems Biology. New York, NY: Springer (2013). pp. 2258–9. doi: 10.1007/978-1-4419-9863-7_756

  • 2

    PopeSDMedzhitovR. Emerging Principles of Gene Expression Programs and Their Regulation. Mol Cell (2018) 71(3):389–97. doi: 10.1016/j.molcel.2018.07.017

  • 3

    WeakeVMWorkmanJL. Inducible Gene Expression: Diverse Regulatory Mechanisms. Nat Rev Genet (2010) 11(6):426–37. doi: 10.1038/nrg2781

  • 4

    SmaleSTNatoliG. Transcriptional Control of Inflammatory Responses. Cold Spring Harb Perspect Biol (2014) 6(11):a016261. doi: 10.1101/cshperspect.a016261

  • 5

    FurlongEEMLevineM. Developmental Enhancers and Chromosome Topology. Science (2018) 361(6409):1341–5. doi: 10.1126/science.aau0320

  • 6

    VossTCHagerGL. Dynamic Regulation of Transcriptional States by Chromatin and Transcription Factors. Nat Rev Genet (2014) 15:6981. doi: 10.1038/nrg3623

  • 7

    LeeTIYoungRA. Transcriptional Regulation and its Misregulation in Disease. Cell (2013) 152(6):1237–51. doi: 10.1016/j.cell.2013.02.014

  • 8

    RichelsRShilatifardA. Enhancer Logic and Mechanics in Development and Disease. Trends Cell Bio (2018) 28(8):608–30. doi: 10.1016/j.tcb.2018.04.003

  • 9

    SpitzFFurlongEE. Transcription Factors: From Enhancer Binding to Developmental Control. Nat Rev Genet (2012) 13(9):613–26. doi: 10.1038/nrg3207

  • 10

    PennacchioLABickmoreWDeanANobregaMABejeranoG. Enhancers: Five Essential Questions. Nat Rev Genet (2013) 14(4):288–95. doi: 10.1038/nrg3458

  • 11

    CaloEWysockaJ. Modification of Enhancer Chromatin: What, How, and Why? Mol Cell (2013) 49(5):825–37. doi: 10.1016/j.molcel.2013.01.038

  • 12

    HeinzSRomanoskiCEBennerCGlassCK. The Selection and Function of Cell Type-Specific Enhancers. Nat Rev Mol Cell Biol (2015) 16(3):144–54. doi: 10.1038/nrm3949

  • 13

    ReiterFWienerroitherSStarkA. Combinatorial Function of Transcription Factors and Cofactors. Curr Opin Genet Dev (2017) 43:7381. doi: 10.1016/j.gde.2016.12.007

  • 14

    CatarinoRRStarkA. Assessing Sufficiency and Necessity of Enhancer Activities for Gene Expression and the Mechanisms of Transcription Activation. Genes Dev (2018) 32(3-4):202–23. doi: 10.1101/gad.310367.117

  • 15

    OudelaarAMDaviesJOJHanssenLLPTeleniusJMSchwessingerRLiuYet al. Single-Allele Chromatin Interactions Identify Regulatory Hubs in Dynamic Compartmentalized Domains. Nat Genet (2018) 50:1744–51. doi: 10.1038/s41588-018-0253-2

  • 16

    BhuiyanTTimmersHTM. Promoter Recognition: Putting TFIID on the Spot. Trends Cell Biol (2019) 29(9):752–63. doi: 10.1016/j.tcb.2019.06.004

  • 17

    McKnightSLKingsburyR. Transcriptional Control Signals of a Eukaryotic Protein-Coding Gene. Science (1982) 217(4557):316–24. doi: 10.1126/science.6283634

  • 18

    MyersRMTillyKManiatisT. Fine Structure Genetic Analysis of a Beta-Globin Promoter. Science (1986) 232(4750):613–8. doi: 10.1126/science.3457470

  • 19

    CalhounVCStathopoulosALevineM. Promoter-Proximal Tethering Elements Regulate Enhancer-Promoter Specificity in the Drosophila Antennapedia Complex. Proc Natl Acad Sci U S A (2002) 99(14):9243–7. doi: 10.1073/pnas.142291299

  • 20

    NolisIKMcKayDJMantouvalouELomvardasSMerikaMThanosD. Transcription Factors Mediate Long-Range Enhancer-Promoter Interactions. Proc Natl Acad Sci U S A (2009) 106(48):20222–7. doi: 10.1073/pnas.0902454106

  • 21

    ViselARubinEMPennacchioLA. Genomic Views of Distant-Acting Enhancers. Nature (2009) 461(7261):199205. doi: 10.1038/nature08451

  • 22

    SchoenfelderSFraserP. Long-Range Enhancer-Promoter Contacts in Gene Expression Control. Nat Rev Genet (2019) 20(8):437–55. doi: 10.1038/s41576-019-0128-0

  • 23

    GuerreroLMarco-FerreresRSerranoALArredondoJJCerveraM. Secondary Enhancers Synergise With Primary Enhancers to Guarantee Fine-Tuned Muscle Gene Expression. Dev Biol (2010) 337(1):1628. doi: 10.1016/j.ydbio.2009.10.006

  • 24

    ThanosDManiatisT. Virus Induction of Human IFN Beta Gene Expression Requires the Assembly of an Enhanceosome. Cell (1995) 83(7):1091–100. doi: 10.1016/0092-8674(95)90136-1

  • 25

    AgelopoulosMThanosD. Epigenetic Determination of a Cell-Specific Gene Expression Program by ATF-2 and the Histone Variant Macroh2a. EMBO J (2006) 25(20):4843–53. doi: 10.1038/sj.emboj.7601364

  • 26

    SlatteryMRileyTLiuPAbeNGomez-AlcalaPDrorIet al. Cofactor Binding Evokes Latent Differences in DNA Binding Specificity Between Hox Proteins. Cell (2011) 147(6):1270–82. doi: 10.1016/j.cell.2011.10.053

  • 27

    FarleyEKOlsonKMZhangWRokhsarDSLevineMS. Syntax Compensates for Poor Binding Sites to Encode Tissue Specificity of Developmental Enhancers. Proc Natl Acad Sci U S A (2016) 113(23):6508–13. doi: 10.1073/pnas.1605085113

  • 28

    LavigneMDVatsellasGPolyzosAMantouvalouESianidisGMaraziotisIet al. Composite Macroh2a/NRF-1 Nucleosomes Suppress Noise and Generate Robustness in Gene Expression. Cell Rep (2015) 11(7):1090–101. doi: 10.1016/j.celrep.2015.04.022

  • 29

    InukaiSKockKHBulykML. Transcription factor-DNA Binding: Beyond Binding Site Motifs. Curr Opin Genet Dev (2017) 43:110–9. doi: 10.1016/j.gde.2017.02.007

  • 30

    StadhoudersRFillionGJGrafT. Transcription Factors and 3D Genome Conformation in Cell-Fate Decisions. Nature (2019) 569:345–54. doi: 10.1038/s41586-019-1182-7

  • 31

    VaquerizasJMKummerfeldSKTeichmannSALuscombeNM. A Census of Human Transcription Factors: Function, Expression and Evolution. Nat Rev Genet (2009) 10(4):252–63. doi: 10.1038/nrg2538

  • 32

    LambertSAJolmaACampitelliLFDasPKYinYAlbuMet al. The Human Transcription Factors. Cell (2018) 172(4):650–65. doi: 10.1016/j.cell.2018.01.029

  • 33

    KribelbauerJFRastogiCBussemakerHJMannRS. Low-Affinity Binding Sites and the Transcription Factor Specificity Paradox in Eukaryotes. Annu Rev Cell Dev Biol (2019) 35:357–79. doi: 10.1146/annurev-cellbio-100617-062719

  • 34

    HottigerMOFelzienLKNabelGJ. Modulation of Cytokine-Induced HIV Gene Expression by Competitive Binding of Transcription Factors to the Coactivator P300. EMBO J (1998) 17(11):3124–34. doi: 10.1093/emboj/17.11.3124

  • 35

    AntonakiADemetriadesCPolyzosABanosAVatsellasGLavigneMDet al. Genomic Analysis Reveals a Novel Nuclear Factor-κb (Nf-κb)-Binding Site in Alu-repetitive Elements. J Biol Chem (2011) 286(44):38768–82. doi: 10.1074/jbc.M111.234161

  • 36

    WongDTeixeiraAOikonomopoulosSHumburgPLoneINSalibaDet al. Extensive Characterization of NF-κb Binding Uncovers non-Canonical Motifs and Advances the Interpretation of Genetic Functional Traits. Genome Biol (2011) 12(7):R70. doi: 10.1186/gb-2011-12-7-r70

  • 37

    SimicevicJSchmidAWGilardoniPAZollerBRaghavSKKrierIet al. Absolute Quantification of Transcription Factors During Cellular Differentiation Using Multiplexed Targeted Proteomics. Nat Methods (2013) 10(6):570–6. doi: 10.1038/nmeth.2441

  • 38

    GillespieMAPaliiCGSanchez-TaltavullDShannonPLongabaughWJRDownesDJet al. Absolute Quantification of Transcription Factors Reveals Principles of Gene Regulation in Erythropoiesis. Mol Cell (2020) 78(5):96074.e11. doi: 10.1016/j.molcel.2020.03.031

  • 39

    VockleyCMD’IppolitoAMMcDowellICMajorosWHSafiASongLet al. Direct GR Binding Sites Potentiate Clusters of TF Binding Across the Human Genome. Cell (2016) 166(5):126981.e19. doi: 10.1016/j.cell.2016.07.049

  • 40

    ReddyTEPauliFSprouseRONeffNFNewberryKMGarabedianMJet al. Genomic Determination of the Glucocorticoid Response Reveals Unexpected Mechanisms of Gene Regulation. Genome Res (2009) 19(12):2163–71. doi: 10.1101/gr.097022.109

  • 41

    FreaneyJEKimRMandhanaRHorvathCM. Extensive Cooperation of Immune Master Regulators IRF3 and Nfκb in RNA Pol II Recruitment and Pause Release in Human Innate Antiviral Transcription. Cell Rep (2013) 4(5):959–73. doi: 10.1016/j.celrep.2013.07.043

  • 42

    BarakatTSHalbritterFZhangMRendeiroAFPerenthalerEBockCet al. Functional Dissection of the Enhancer Repertoire in Human Embryonic Stem Cells. Cell Stem Cell (2018) 23(2):276288.e8. doi: 10.1016/j.stem.2018.06.014

  • 43

    YanJEngeMWhitingtonTDaveKLiuJSurIet al. Transcription Factor Binding in Human Cells Occurs in Dense Clusters Formed Around Cohesin Anchor Sites. Cell (2013) 154(4):801–13. doi: 10.1016/j.cell.2013.07.034

  • 44

    PartridgeECChhetriSBProkopJWRamakerRCJansenCSGohSTet al. Occupancy Maps of 208 Chromatin-Associated Proteins in One Human Cell Type. Nature (2020) 583(7818):720–8. doi: 10.1038/s41586-020-2023-4

  • 45

    WhyteWAOrlandoDAHniszDAbrahamBJLinCYKageyMHet al. Master Transcription Factors and Mediator Establish Super-Enhancers At Key Cell Identity Genes. Cell (2013) 153(2):307–19. doi: 10.1016/j.cell.2013.03.035

  • 46

    ShinHYWilliMHyunYooKZengXWangCMetserGet al. Hierarchy Within the Mammary STAT5-driven Wap Super-Enhancer. Nat Genet (2016) 48(8):904–11. doi: 10.1038/ng.3606

  • 47

    HayDHughesJRBabbsCDaviesJOJGrahamBJHanssenLet al. Genetic Dissection of the α-Globin Super-Enhancer In Vivo. Nat Genet (2016) 48(8):895903. doi: 10.1038/ng.3605

  • 48

    DuklerNGulkoBHuangYFSiepelA. Is a Super-Enhancer Greater Than the Sum of its Parts? Nat Genet (2016) 49(1):23. doi: 10.1038/ng.3759

  • 49

    Ing-SimmonsESeitanVCFaureAJFlicekPCarrollTDekkerJet al. Spatial Enhancer Clustering and Regulation of Enhancer-Proximal Genes by Cohesion. Genome Res (2015) 25:504–13. doi: 10.1101/gr.184986.114

  • 50

    BrownJDLinCYDuanQGriffinGFederationAParanalRMet al. Nf-κb Directs Dynamic Super Enhancer Formation in Inflammation and Atherogenesis. Mol Cell (2014) 56(2):219–31. doi: 10.1016/j.molcel.2014.08.024

  • 51

    BahrCvon PaleskeLUsluVVRemeseiroSTakayamaNNgSWet al. A Myc Enhancer Cluster Regulates Normal and Leukaemic Haematopoietic Stem Cell Hierarchies. Nature (2018) 553(7689):515–20. doi: 10.1038/nature25193

  • 52

    ENCODE Project Consortium. An Integrated Encyclopedia of DNA Elements in the Human Genome. Nature (2012) 489(7414):5774. doi: 10.1038/nature11247

  • 53

    ENCODE Project ConsortiumMooreJEPurcaroMJPrattHEEpsteinCBShoreshNet al. Expanded Encyclopaedias of DNA Elements in the Human and Mouse Genomes. Nature (2020) 583(7818):699710. doi: 10.1038/s41586-020-2493-4

  • 54

    AgelopoulosMMcKayDJMannRS. Developmental Regulation of Chromatin Conformation by Hox Proteins in Drosophila. Cell Rep (2012) 1(4):350–9. doi: 10.1016/j.celrep.2012.03.003

  • 55

    McKayDJEstellaCMannRS. The Origins of the Drosophila Leg Revealed by the Cis-Regulatory Architecture of the Distalless Gene. Development (2009) 136(1):6171. doi: 10.1242/dev.029975

  • 56

    MohrsMBlankespoorCMWangZELootsGGAfzalVHadeibaHet al. Deletion of a Coordinate Regulator of Type 2 Cytokine Expression in Mice. Nat Immunol (2001) 2(9):842–7. doi: 10.1038/ni0901-842

  • 57

    Ghavi-HelmYKleinFAPakozdiTCiglarLNoordermeerDHuberWet al. Enhancer Loops Appear Stable During Development and are Associated With Paused Polymerase. Nature (2014) 512(7512):96100. doi: 10.1038/nature13417

  • 58

    DengWLeeJWangHMillerJReikAGregoryPDet al. Controlling Long-Range Genomic Interactions At a Native Locus by Targeted Tethering of a Looping Factor. Cell (2012) 149(6):1233–44. doi: 10.1016/j.cell.2012.03.051

  • 59

    Dall’AgneseACaputoLNicolettiCdi IulioJSchmittAGattoSet al. Transcription Factor-Directed Re-Wiring of Chromatin Architecture for Somatic Cell Nuclear Reprogramming Toward Trans-Differentiation. Mol Cell (2019) 76(3):453472.e8. doi: 10.1016/j.molcel.2019.07.036

  • 60

    Di GiammartinoDCKloetgenAPolyzosALiuYKimDMurphyDet al. KLF4 is Involved in the Organization and Regulation of Pluripotency-Associated Three-Dimensional Enhancer Networks. Nat Cell Biol (2019) 21(10):1179–90. doi: 10.1038/s41556-019-0390-6

  • 61

    BenabdallahNSWilliamsonIIllingworthRSKaneLBoyleSSenguptaDet al. Decreased Enhancer-Promoter Proximity Accompanying Enhancer Activation. Mol Cell (2019) 76(3):47384.e7. doi: 10.1016/j.molcel.2019.07.038

  • 62

    WilliamsonILetticeLAHillREBickmoreWA. Shh and ZRS Enhancer Colocalisation is Specific to the Zone of Polarising Activity. Development (2016) 143(16):29943001. doi: 10.1242/dev.139188

  • 63

    AlexanderJMGuanJLiBMaliskovaLSongMShenYet al. Live-Cell Imaging Reveals Enhancer-Dependent Sox2 Transcription in the Absence of Enhancer Proximity. Elife (2019) 8:e41769. doi: 10.7554/eLife.41769

  • 64

    Lieberman-AidenEvan BerkumNLWilliamsLImakaevMRagoczyTTellingAet al. Comprehensive Mapping of Long-Range Interactions Reveals Folding Principles of the Human Genome. Science (2009) 326(5950):289–93. doi: 10.1126/science.1181369

  • 65

    RowleyMJCorcesVG. Organizational Principles of 3D Genome Architecture. Nat Rev Genet (2018) 19(12):789800. doi: 10.1038/s41576-018-0060-8

  • 66

    RaoSSHuntleyMHDurandNCStamenovaEKBochkovIDRobinsonJTet al. A 3D Map of the Human Genome At Kilobase Resolution Reveals Principles of Chromatin Looping. Cell (2014) 159(7):1665–80. doi: 10.1016/j.cell.2014.11.021

  • 67

    DixonJRSelvarajSYueFKimALiYShenYet al. Topological Domains in Mammalian Genomes Identified by Analysis of Chromatin Interactions. Nature (2012) 485(7398):376–80. doi: 10.1038/nature11082

  • 68

    GhirlandoRFelsenfeldG. CTCF: Making the Right Connections. Genes Dev (2016) 30(8):881–91. doi: 10.1101/gad.277863.116

  • 69

    FudenbergGImakaevMLuCGoloborodkoAAbdennurNMirnyLA. Formation of Chromosomal Domains by Loop Extrusion. Cell Rep (2016) 15(9):2038–49. doi: 10.1016/j.celrep.2016.04.085

  • 70

    LupiáñezDGKraftKHeinrichVKrawitzPBrancatiFKlopockiEet al. Disruptions of Topological Chromatin Domains Cause Pathogenic Rewiring of Gene-Enhancer Interactions. Cell (2015) 161(5):1012–25. doi: 10.1016/j.cell.2015.04.004

  • 71

    FlavahanWADrierYLiauBBGillespieSMVenteicherASStemmer-RachamimovAOet al. Insulator Dysfunction and Oncogene Activation in IDH Mutant Gliomas. Nature (2016) 529(7584):110–4. doi: 10.1038/nature16490

  • 72

    DespangASchöpflinRFrankeMAliSJerkovićIPaliouCet al. Functional Dissection of the Sox9-Kcnj2 Locus Identifies Nonessential and Instructive Roles of TAD Architecture. Nat Genet (2019) 51(8):1263–71. doi: 10.1038/s41588-019-0466-z

  • 73

    GuoYXuQCanzioDShouJLiJGorkinDUet al. Crispr Inversion of CTCF Sites Alters Genome Topology and Enhancer/Promoter Function. Cell (2015) 162(4):900–10. doi: 10.1016/j.cell.2015.07.038

  • 74

    ZhouJLevineM. A Novel Cis-Regulatory Element, the PTS, Mediates an Anti-Insulator Activity in the Drosophila Embryo. Cell (1999) 99(6):567–75. doi: 10.1016/s0092-8674(00)81546-9

  • 75

    GrosveldFvan StaalduinenJStadhoudersR. Transcriptional Regulation by (Super)Enhancers: From Discovery to Mechanisms. Annu Rev Genomics Hum Genet (2021) 22:13.1–20. doi: 10.1146/annurev-genom-122220-093818

  • 76

    SymmonsOPanLRemeseiroSAktasTKleinFHuberWet al. The Shh Topological Domain Facilitates the Action of Remote Enhancers by Reducing the Effects of Genomic Distances. Dev Cell (2016) 39(5):529–43. doi: 10.1016/j.devcel.2016.10.015

  • 77

    NoraEPGoloborodkoAValtonALGibcusJHUebersohnAAbdennurNet al. Targeted Degradation of CTCF Decouples Local Insulation of Chromosome Domains From Genomic Compartmentalization. Cell (2017) 169(5):93044.e22. doi: 10.1016/j.cell.2017.05.004

  • 78

    RaoSSPHuangSCGlenn St HilaireBEngreitzJMPerezEMKieffer-KwonKRet al. Cohesin Loss Eliminates All Loop Domains. Cell (2017) 171(2):305320.e24. doi: 10.1016/j.cell.2017.09.026

  • 79

    CattoniDICardozo GizziAMGeorgievaMDi StefanoMValeriAChamoussetDet al. Single-Cell Absolute Contact Probability Detection Reveals Chromosomes are Organized by Multiple Low-Frequency Yet Specific Interactions. Nat Commun (2017) 8(1):1753. doi: 10.1038/s41467-017-01962-x

  • 80

    FinnEHPegoraroGBrandãoHBValtonALOomenMEDekkerJet al. Extensive Heterogeneity and Intrinsic Variation in Spatial Genome Organization. Cell (2019) 176(6):15021515.e10. doi: 10.1016/j.cell.2019.01.020

  • 81

    DekkerJRippeKDekkerMKlecknerN. Capturing Chromosome Conformation. Science (2002) 295(5558):1306–11. doi: 10.1126/science.1067799

  • 82

    FullwoodMJLiuMHPanYFLiuJXuHMohamedYBet al. An Oestrogen-Receptor-Alpha-Bound Human Chromatin Interactome. Nature (2009) 462(7269):5864. doi: 10.1038/nature08497

  • 83

    HughesJRRobertsNMcGowanSHayDGiannoulatouELynchMet al. Analysis of Hundreds of Cis-Regulatory Landscapes At High Resolution in a Single, High-Throughput Experiment. Nat Genet (2014) 46(2):205–12. doi: 10.1038/ng.2871

  • 84

    MumbachMRRubinAJFlynnRADaiCKhavariPAGreenleafWJet al. HiChIP: Efficient and Sensitive Analysis of Protein-Directed Genome Architecture. Nat Methods (2016) 13(11):919–22. doi: 10.1038/nmeth.3999

  • 85

    ChenHLevoMBarinovLFujiokaMJaynesJBGregorT. Dynamic Interplay Between Enhancer-Promoter Topology and Gene Activity. Nat Genet (2018) 50(9):1296–303. doi: 10.1038/s41588-018-0175-z

  • 86

    TolhuisBPalstraRJSplinterEGrosveldFde LaatW. Looping and Interaction Between Hypersensitive Sites in the Active Beta-Globin Locus. Mol Cell (2002) 10(6):1453–65. doi: 10.1016/s1097-2765(02)00781-5

  • 87

    SpilianakisCGLaliotiMDTownTLeeGRFlavellRA. Interchromosomal Associations Between Alternatively Expressed Loci. Nature (2005) 435(7042):637–45. doi: 10.1038/nature03574

  • 88

    ApostolouEThanosD. Virus Infection Induces NF-kappaB-dependent Interchromosomal Associations Mediating Monoallelic IFN-beta Gene Expression. Cell (2008) 134(1):8596. doi: 10.1016/j.cell.2008.05.052

  • 89

    NikopoulouCPanagopoulosGSianidisGPsarraEFordEThanosD. The Transcription Factor ThPOK Orchestrates Stochastic Interchromosomal Interactions Required for IFNB1 Virus-Inducible Gene Expression. Mol Cell (2018) 71(2):35261.e5. doi: 10.1016/j.molcel.2018.06.019

  • 90

    ThanosDManiatisT. The High Mobility Group Protein HMG I(Y) is Required for NF-kappa B-Dependent Virus Induction of the Human IFN-beta Gene. Cell (1992) 71(5):777–89. doi: 10.1016/0092-8674(92)90554-p

  • 91

    AgaliotiTLomvardasSParekhBYieJManiatisTThanosD. Ordered Recruitment of Chromatin Modifying and General Transcription Factors to the IFN-beta Promoter. Cell (2000) 103(4):667–78. doi: 10.1016/s0092-8674(00)00169-0

  • 92

    MunshiNAgaliotiTLomvardasSMerikaMChenGThanosD. Coordination of a Transcriptional Switch by HMGI(Y) Acetylation. Science (2001) 293(5532):1133–6. doi: 10.1126/science.293.5532.1133

  • 93

    MerikaMWilliamsAJChenGCollinsTThanosD. Recruitment of CBP/p300 by the IFN Beta Enhanceosome is Required for Synergistic Activation of Transcription. Mol Cell (1998) 1(2):277–87. doi: 10.1016/s1097-2765(00)80028-3

  • 94

    MerikaMThanosD. Enhanceosomes. Curr Opin Genet Dev (2001) 11(2):205–8. doi: 10.1016/s0959-437x(00)00180-5

  • 95

    PanneDManiatisTHarrisonSC. An Atomic Model of the Interferon-Beta Enhanceosome. Cell (2007) 129(6):1111–23. doi: 10.1016/j.cell.2007.05.019

  • 96

    KobayashiKSvan den ElsenPJ. NLRC5: A Key Regulator of MHC Class I-dependent Immune Responses. Nat Rev Immunol (2012) 12(12):813–20. doi: 10.1038/nri3339

  • 97

    LomvardasSThanosD. Nucleosome Sliding Via TBP DNA Binding In Vivo. Cell (2001) 106(6):685–96. doi: 10.1016/s0092-8674(01)00490-1

  • 98

    LomvardasSThanosD. Modifying Gene Expression Programs by Altering Core Promoter Chromatin Architecture. Cell (2002) 110(2):261–71. doi: 10.1016/s0092-8674(02)00822-x

  • 99

    RandURinasMSchwerkJNöhrenGLinnesMKrögerAet al. Multi-Layered Stochasticity and Paracrine Signal Propagation Shape the Type-I Interferon Response. Mol Syst Biol (2012) 8:584. doi: 10.1038/msb.2012.17

  • 100

    WimmersFSubediNvan BuuringenNHeisterDViviéJBeeren-ReinierenIet al. Single-Cell Analysis Reveals That Stochasticity and Paracrine Signaling Control Interferon-Alpha Production by Plasmacytoid Dendritic Cells. Nat Commun (2018) 9(1):3317. doi: 10.1038/s41467-018-05784-3

  • 101

    FangMXieHDouganSKPloeghHvan OudenaardenA. Stochastic Cytokine Expression Induces Mixed T Helper Cell States. PloS Biol (2013) 11(7):e1001618. doi: 10.1371/journal.pbio.1001618

  • 102

    ZhaoMZhangJPhatnaniHScheuSManiatisT. Stochastic Expression of the Interferon-β Gene. PloS Biol (2012) 10(1):e1001249. doi: 10.1371/journal.pbio.1001249

  • 103

    WeintraubASLiCHZamudioAVSigovaAAHannettNMDayDSet al. Yy1 Is a Structural Regulator of Enhancer-Promoter Loops. Cell (2017) 171(7):157388.e28. doi: 10.1016/j.cell.2017.11.008

  • 104

    BaileySDZhangXDesaiKAidMCorradinOCowper-Sal LariRet al. ZNF143 Provides Sequence Specificity to Secure Chromatin Interactions At Gene Promoters. Nat Commun (2015) 2:6186. doi: 10.1038/ncomms7186

  • 105

    CaiSHanHJKohwi-ShigematsuT. Tissue-Specific Nuclear Architecture and Gene Expression Regulated by SATB1. Nat Genet (2003) 34(1):4251. doi: 10.1038/ng1146

  • 106

    CaiSLeeCCKohwi-ShigematsuT. SATB1 Packages Densely Looped, Transcriptionally Active Chromatin for Coordinated Expression of Cytokine Genes. Nat Genet (2006) 38(11):1278–88. doi: 10.1038/ng1913

  • 107

    GhoshRPShiQYangLReddickMPNikitinaTZhurkinVBet al. Satb1 Integrates DNA Binding Site Geometry and Torsional Stress to Differentially Target Nucleosome-Dense Regions. Nat Commun (2019) 10(1):3221. doi: 10.1038/s41467-019-11118-8

  • 108

    LomvardasSBarneaGPisapiaDJMendelsohnMKirklandJAxelR. Interchromosomal Interactions and Olfactory Receptor Choice. Cell (2006) 126(2):403–13. doi: 10.1016/j.cell.2006.06.035

  • 109

    Markenscoff-PapadimitriouEAllenWEColquittBMGohTMurphyKKMonahanKet al. Enhancer Interaction Networks as a Means for Singular Olfactory Receptor Expression. Cell (2014) 159(3):543–57. doi: 10.1016/j.cell.2014.09.033

  • 110

    MonahanKSchierenICheungJMumbey-WafulaAMonukiESLomvardasS. Cooperative Interactions Enable Singular Olfactory Receptor Expression in Mouse Olfactory Neurons. Elife (2017) 6:e28620. doi: 10.7554/eLife.28620

  • 111

    MonahanKHortaALomvardasS. LHX2- and LDB1-mediated Trans Interactions Regulate Olfactory Receptor Choice. Nature (2019) 565(7740):448–53. doi: 10.1038/s41586-018-0845-0

  • 112

    FanucchiSShibayamaYBurdSWeinbergMSMhlangaMM. Chromosomal Contact Permits Transcription Between Coregulated Genes. Cell (2013) 155(3):606–20. doi: 10.1016/j.cell.2013.09.051

  • 113

    BanosAAgelopoulosMThanosD. Stochastic Responses are Not Left to Pure “Chance”. Cell (2013) 155(3):499502. doi: 10.1016/j.cell.2013.10.002

  • 114

    WernetMFMazzoniEOÇelikADuncanDMDuncanIDesplanC. Stochastic Spineless Expression Creates the Retinal Mosaic for Colour Vision. Nature (2006) 440(7081):174–80. doi: 10.1038/nature04615

  • 115

    AndersonCReissIZhouCChoASiddiqiHMormannBet al. Natural Variation in Stochastic Photoreceptor Specification and Color Preference in Drosophila. Elife (2017) 6:e29593. doi: 10.7554/eLife.29593

  • 116

    JohnstonRJJrDesplanC. Interchromosomal Communication Coordinates Intrinsically Stochastic Expression Between Alleles. Science (2014) 343(6171):661–5. doi: 10.1126/science.1243039

  • 117

    ShrinivasKSabariBRCoffeyELKleinIABoijaAZamudioAVet al. Enhancer Features That Drive Formation of Transcriptional Condensates. Mol Cell (2019) 75(3):54961.e7. doi: 10.1016/j.molcel.2019.07.009

  • 118

    ChoW-KSpilleJ-HHechtMLeeCLiCGrubeVet al. Mediator and RNA Polymerase II Clusters Associate in Transcription-Dependent Condensates. Science (2018) 361, 6400:412–5. doi: 10.1126/science.aar4199

  • 119

    van der LeeRBuljanMLangBWeatherittRJDaughdrillGWDunkerAKet al. Classification of Intrinsically Disordered Regions and Proteins. Chem Rev (2014) 114(13):6589–631. doi: 10.1021/cr400525m

  • 120

    BoijaAKleinIASabariBRDall’AgneseACoffeyELZamudioAVet al. Transcription Factors Activate Genes Through the Phase-Separation Capacity of Their Activation Domains. Cell (2018) 175(7):1842–55. doi: 10.1016/j.cell.2018.10.042

  • 121

    BrangwynneCPEckmannCRCoursonDSRybarskaAHoegeCGharakhaniJet al. Germline P Granules are Liquid Droplets That Localize by Controlled Dissolution/Condensation. Science (2009) 324(5935):1729–32. doi: 10.1126/science.1172046

  • 122

    KloetgenAThandapaniPTsirigosAAifantisI. 3d Chromosomal Landscapes in Hematopoiesis and Immunity. Trends Immunol (2019) 40(9):809–24. doi: 10.1016/j.it.2019.07.003

  • 123

    HniszDShrinivasKYoungRAChakrabortyAKSharpPA. A Phase Separation Model for Transcriptional Control. Cell (2017) 169(1):1323. doi: 10.1016/j.cell.2017.02.007

  • 124

    NairSJYangLMeluzziDOhSYangFFriedmanMJet al. Phase Separation of Ligand-Activated Enhancers Licenses Cooperative Chromosomal Enhancer Assembly. Nat Struct Mol Biol (2019) 26(3):193203. doi: 10.1038/s41594-019-0190-5

  • 125

    ZamudioAVDall’AgneseAHenningerJEManteigaJCAfeyanLKHannettNMet al. Mediator Condensates Localize Signaling Factors to Key Cell Identity Genes. Mol Cell (2019) 76(5):75366.e6. doi: 10.1016/j.molcel.2019.08.016

  • 126

    SabariBRDall’AgneseABoijaAKleinIACoffeyELShrinivasKet al. Coactivator Condensation At Super-Enhancers Links Phase Separation and Gene Control. Science (2018) 361(6400):eaar3958. doi: 10.1126/science.aar3958

Summary

Keywords

regulation of transcription, enhancers, transcription factors, chromatin, enhancer hubs, stochastic expression

Citation

Agelopoulos M, Foutadakis S and Thanos D (2021) The Causes and Consequences of Spatial Organization of the Genome in Regulation of Gene Expression. Front. Immunol. 12:682397. doi: 10.3389/fimmu.2021.682397

Received

18 March 2021

Accepted

18 May 2021

Published

04 June 2021

Volume

12 - 2021

Edited by

Ageliki Tsagaratou, University of North Carolina at Chapel Hill, United States

Reviewed by

Chris Benner, University of California, San Diego, United States; Wing Fuk Chan, Walter and Eliza Hall Institute of Medical Research, Australia

Updates

Copyright

*Correspondence: Dimitris Thanos,

This article was submitted to B Cell Biology, a section of the journal Frontiers in Immunology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics