Transcription-driven genome organization: a model for chromosome structure and the regulation of gene expression tested through simulations

Genome-wide chromosome architecture prediction reveals biophysical principles underlying gene structure

Classical observations suggest a connection between 3D gene structure and function, but testing this hypothesis has been challenging due to technical limitations. To explore this, we developed epigenetic highly predictive heteromorphic polymer (e-HiP-HoP), a model based on genome organization principles to predict the 3D structure of human chromatin. We defined a new 3D structural unit, a "topos," which represents the regulatory landscape around gene promoters. Using GM12878 cells, we predicted the 3D structure of over 10,000 active gene topoi and stored them in the 3DGene database. Data mining revealed folding motifs and their link to Gene Ontology features. We computed a structural diversity score and identified influential nodes-chromatin sites that frequently interact with gene promoters, acting as key regulators. These nodes drive structural diversity and are tied to gene function. e-HiP-HoP provides a framework for modeling high-resolution chromatin structure and a mechanistic basis for chromatin contact networks that link 3D gene structure with function.

Diverse epistatic effects in barley-powdery mildew interactions localize to host chromosome hotspots

Barley Mildew locus a (Mla) encodes a multi-allelic series of nucleotide-binding leucine-rich repeat (NLR) receptors that specify recognition to diverse cereal diseases. We exploited time-course transcriptome dynamics of barley and derived immune mutants infected with the powdery mildew fungus, Blumeria hordei (Bh), to infer gene effects governed by Mla6 and two other loci significant to disease development, Blufensin1 (Bln1), and Required for Mla6 resistance3 (rar3 = Sgt1 ΔKL308-309 ). Interactions of Mla6 and Bln1 resulted in diverse epistatic effects on the Bh-induced barley transcriptome, differential immunity to Pseudomonas syringae expressing the effector protease AvrPphB, and reaction to Bh. From a total of 468 barley NLRs, 115 were grouped under different gene effect models; genes classified under these models localized to host chromosome hotspots. The corresponding Bh infection transcriptome was classified into nine co-expressed modules, linking differential expression with pathogen structures, signifying that disease is regulated by an inter-organismal network that diversifies the response.

Enhancers on the edge - how the nuclear envelope controls gene regulatory elements

Precise temporal and sequential control of gene expression during development and in response to environmental stimuli requires tight regulation of the physical contact between gene regulatory elements and promoters. Current models describing how the genome folds in 3D space to establish these interactions often ignore the role of the most stable structural nuclear feature - the nuclear envelope. While contributions of 3D folding within/between topologically associated domains (TADs) have been extensively described, mechanical contributions from the nuclear envelope can impact enhancer-promoter interactions both directly and indirectly through influencing intra/inter-TAD interactions. Importantly, these nuclear envelope contributions clearly link this mechanism to development and, when defective, to human disease. Here, we discuss evidence for nuclear envelope regulation of tissue-specific enhancer-promoter pairings, potential mechanisms for this regulation, exciting recent findings that other regulatory elements such as microRNAs and long noncoding RNAs are under nuclear envelope regulation, the possible involvement of condensates, and how disruption of this regulation can lead to disease.

Transcription regulates the spatio-temporal dynamics of genes through micro-compartmentalization

Although our understanding of the involvement of heterochromatin architectural factors in shaping nuclear organization is improving, there is still ongoing debate regarding the role of active genes in this process. In this study, we utilize publicly-available Micro-C data from mouse embryonic stem cells to investigate the relationship between gene transcription and 3D gene folding. Our analysis uncovers a nonmonotonic - globally positive - correlation between intragenic contact density and Pol II occupancy, independent of cohesin-based loop extrusion. Through the development of a biophysical model integrating the role of transcription dynamics within a polymer model of chromosome organization, we demonstrate that Pol II-mediated attractive interactions with limited valency between transcribed regions yield quantitative predictions consistent with chromosome-conformation-capture and live-imaging experiments. Our work provides compelling evidence that transcriptional activity shapes the 4D genome through Pol II-mediated micro-compartmentalization.

Topological Spectra and Entropy of Chromatin Loop Networks

The 3D folding of a mammalian gene can be studied by a polymer model, where the chromatin fiber is represented by a semiflexible polymer which interacts with multivalent proteins, representing complexes of DNA-binding transcription factors and RNA polymerases. This physical model leads to the natural emergence of clusters of proteins and binding sites, accompanied by the folding of chromatin into a set of topologies, each associated with a different network of loops. Here, we combine numerics and analytics to first classify these networks and then find their relative importance or statistical weight, when the properties of the underlying polymer are those relevant to chromatin. Unlike polymer networks previously studied, our chromatin networks have finite average distances between successive binding sites, and this leads to giant differences between the weights of topologies with the same number of edges and nodes but different wiring. These weights strongly favor rosettelike structures with a local cloud of loops with respect to more complicated nonlocal topologies. Our results suggest that genes should overwhelmingly fold into a small fraction of all possible 3D topologies, which can be robustly characterized by the framework we propose here.

Combinatorics and topological weights of chromatin loop networks

Polymer physics models suggest that chromatin spontaneously folds into loop networks with transcription units (TUs), such as enhancers and promoters, as anchors. Here we use combinatoric arguments to enumerate the emergent chromatin loop networks, both in the case where TUs are labeled and where they are unlabeled. We then combine these mathematical results with those of computer simulations aimed at finding the inter-TU energy required to form a target loop network. We show that different topologies are vastly different in terms of both their combinatorial weight and energy of formation. We explain the latter result qualitatively by computing the topological weight of a given network-i.e., its partition function in statistical mechanics language-in the approximation where excluded volume interactions are neglected. Our results show that networks featuring local loops are statistically more likely with respect to networks including more nonlocal contacts. We suggest our classification of loop networks, together with our estimate of the combinatorial and topological weight of each network, will be relevant to catalog three-dimensional structures of chromatin fibers around eukaryotic genes, and to estimate their relative frequency in both simulations and experiments.

The chromatin - triple helix connection

Mammalian genomes are extensively transcribed, producing a large number of coding and non-coding transcripts. A large fraction of the nuclear RNAs is physically associated with chromatin, functioning in gene activation and silencing, shaping higher-order genome organisation, such as involvement in long-range enhancer-promoter interactions, transcription hubs, heterochromatin, nuclear bodies and phase transitions. Different mechanisms allow the tethering of these chromatin-associated RNAs (caRNA) to chromosomes, including RNA binding proteins, the RNA polymerases and R-loops. In this review, we focus on the sequence-specific targeting of RNA to DNA by forming triple helical structures and describe its interplay with chromatin. It turns out that nucleosome positioning at triple helix target sites and the nucleosome itself are essential factors in determining the formation and stability of triple helices. The histone H3-tail plays a critical role in triple helix stabilisation, and the role of its epigenetic modifications in this process is discussed.

Shifting landscapes: the role of 3D genomic organizations in gene regulatory strategies

3D genome folding enables the physical storage of chromosomes into the compact volume of a cell's nucleus, allows for the accurate segregation of chromatin to daughter cells, and has been shown to be tightly coupled to the way in which genetic information is converted into transcriptional programs [1-3]. Importantly, this link between chromatin architecture and gene regulation is a selectable feature in which modifications to chromatin organization accompany, or perhaps even drive the establishment of new regulatory strategies with enduring impacts on animal body plan complexity. Here, we discuss the nature of different 3D genome folding systems found across the tree of life, with particular emphasis on metazoans, and the relative influence of these systems on gene regulation. We suggest how the properties of these folding systems have influenced regulatory strategies employed by different lineages and may have catalyzed the partitioning and specialization of genetic programs that enabled multicellularity and organ-grade body plan complexity.

G-Quadruplexes in Nuclear Biomolecular Condensates

G-quadruplexes (G4s) have long been implicated in the regulation of chromatin packaging and gene expression. These processes require or are accelerated by the separation of related proteins into liquid condensates on DNA/RNA matrices. While cytoplasmic G4s are acknowledged scaffolds of potentially pathogenic condensates, the possible contribution of G4s to phase transitions in the nucleus has only recently come to light. In this review, we summarize the growing evidence for the G4-dependent assembly of biomolecular condensates at telomeres and transcription initiation sites, as well as nucleoli, speckles, and paraspeckles. The limitations of the underlying assays and the remaining open questions are outlined. We also discuss the molecular basis for the apparent permissive role of G4s in the in vitro condensate assembly based on the interactome data. To highlight the prospects and risks of G4-targeting therapies with respect to the phase transitions, we also touch upon the reported effects of G4-stabilizing small molecules on nuclear biomolecular condensates.

Enhancer-promoter contact formation requires RNAPII and antagonizes loop extrusion

Homotypic chromatin interactions and loop extrusion are thought to be the two main drivers of mammalian chromosome folding. Here we tested the role of RNA polymerase II (RNAPII) across different scales of interphase chromatin organization in a cellular system allowing for its rapid, auxin-mediated degradation. We combined Micro-C and computational modeling to characterize subsets of loops differentially gained or lost upon RNAPII depletion. Gained loops, extrusion of which was antagonized by RNAPII, almost invariably formed by engaging new or rewired CTCF anchors. Lost loops selectively affected contacts between enhancers and promoters anchored by RNAPII, explaining the repression of most genes. Surprisingly, promoter-promoter interactions remained essentially unaffected by polymerase depletion, and cohesin occupancy was sustained. Together, our findings reconcile the role of RNAPII in transcription with its direct involvement in setting-up regulatory three-dimensional chromatin contacts genome wide, while also revealing an impact on cohesin loop extrusion.

Dynamic chromatin architectures provide insights into the genetics of cattle myogenesis

BACKGROUND

Sharply increased beef consumption is propelling the genetic improvement projects of beef cattle in China. Three-dimensional genome structure is confirmed to be an important layer of transcription regulation. Although genome-wide interaction data of several livestock species have already been produced, the genome structure states and its regulatory rules in cattle muscle are still limited.

RESULTS

Here we present the first 3D genome data in Longissimus dorsi muscle of fetal and adult cattle (Bos taurus). We showed that compartments, topologically associating domains (TADs), and loop undergo re-organization and the structure dynamics were consistent with transcriptomic divergence during muscle development. Furthermore, we annotated cis-regulatory elements in cattle genome during myogenesis and demonstrated the enrichments of promoter and enhancer in selection sweeps. We further validated the regulatory function of one HMGA2 intronic enhancer near a strong sweep region on primary bovine myoblast proliferation.

CONCLUSIONS

Our data provide key insights of the regulatory function of high order chromatin structure and cattle myogenic biology, which will benefit the progress of genetic improvement of beef cattle.

MINE is a method for detecting spatial density of regulatory chromatin interactions based on a multi-modal network

Chromatin interactions play essential roles in chromatin conformation and gene expression. However, few tools exist to analyze the spatial density of regulatory chromatin interactions (SD-RCI). Here, we present the multi-modal network (MINE) toolkit, including MINE-Loop, MINE-Density, and MINE-Viewer. The MINE-Loop network aims to enhance the detection of RCIs, MINE-Density quantifies the SD--RCI, and MINE-Viewer facilitates 3D visualization of the density of chromatin interactions and participating regulatory factors (e.g., transcription factors). We applied MINE to investigate the relationship between the SD-RCI and chromatin volume change in HeLa cells before and after liquid-liquid phase separation. Changes in SD-RCI before and after treating the HeLa cells with 1,6-hexanediol suggest that changes in chromatin organization was related to the degree of activation or repression of genes. Together, the MINE toolkit enables quantitative studies on different aspects of chromatin conformation and regulatory activity.

The miR-430 locus with extreme promoter density forms a transcription body during the minor wave of zygotic genome activation

In anamniote embryos, the major wave of zygotic genome activation starts during the mid-blastula transition. However, some genes escape global genome repression, are activated substantially earlier, and contribute to the minor wave of genome activation. The mechanisms underlying the minor wave of genome activation are little understood. We explored the genomic organization and cis-regulatory mechanisms of a transcription body, in which the minor wave of genome activation is first detected in zebrafish. We identified the miR-430 cluster as having excessive copy number and the highest density of Pol-II-transcribed promoters in the genome, and this is required for forming the transcription body. However, this transcription body is not essential for, nor does it encompasse, minor wave transcription globally. Instead, distinct minor-wave-specific promoter architecture suggests that promoter-autonomous mechanisms regulate the minor wave of genome activation. The minor-wave-specific features also suggest distinct transcription initiation mechanisms between the minor and major waves of genome activation.

The balance between gyrase and topoisomerase I activities determines levels of supercoiling, nucleoid compaction, and viability in bacteria

Two enzymes are responsible for maintaining supercoiling in the human pathogen Streptococcus pneumoniae, gyrase (GyrA₂GyrB₂) and topoisomerase I. To attain diverse levels of topoisomerase I (TopoI, encoded by topA), two isogenic strains derived from wild-type strain R6 were constructed: P_Zn topA, carrying an ectopic topA copy under the control of the ZnSO₄-regulated P_Zn promoter and its derivative ΔtopAP_Zn topA, which carries a topA deletion at its native chromosomal location. We estimated the number of TopoI and GyrA molecules per cell by using Western-blot and CFUs counting, and correlated these values with supercoiling levels. Supercoiling was estimated in two ways. We used classical 2D-agarose gel electrophoresis of plasmid topoisomers to determine supercoiling density (σ) and we measured compaction of nucleoids using for the first time super-resolution confocal microscopy. Notably, we observed a good correlation between both supercoiling calculations. In R6, with σ = -0.057, the average number of GyrA molecules per cell (2,184) was higher than that of TopoI (1,432), being the GyrA:TopoI proportion of 1:0.65. In ΔtopAP_Zn topA, the number of TopoI molecules depended, as expected, on ZnSO₄ concentration in the culture media, being the proportions of GyrA:TopoI molecules in 75, 150, and 300 μM ZnSO₄ of 1:0.43, 1:0.47, and 1:0.63, respectively, which allowed normal supercoiling and growth. However, in the absence of ZnSO₄, a higher GyrA:TopoI ratio (1:0.09) caused hyper-supercoiling (σ = -0.086) and lethality. Likewise, growth of ΔtopAP_Zn topA in the absence of ZnSO₄ was restored when gyrase was inhibited with novobiocin, coincidentally with the resolution of hyper-supercoiling (σ change from -0.080 to -0.068). Given that TopoI is a monomer and two molecules of GyrA are present in the gyrase heterotetramer, the gyrase:TopoI enzymes proportion would be 1:1.30 (wild type R6) or of 1:1.26-0.86 (ΔtopAP_Zn topA under viable conditions). Higher proportions, such as 1:0.18 observed in ΔtopAP_Zn topA in the absence of ZnSO₄ yielded to hyper-supercoiling and lethality. These results support a role of the equilibrium between gyrase and TopoI activities in supercoiling maintenance, nucleoid compaction, and viability. Our results shed new light on the mechanism of action of topoisomerase-targeting antibiotics, paving the way for the use of combination therapies.

Spatial genome organization, TGFβ, and biomolecular condensates: Do they talk during development?

Cis-regulatory elements govern gene expression programs to determine cell identity during development. Recently, the possibility that multiple enhancers are orchestrated in clusters of enhancers has been suggested. How these elements are arranged in the 3D space to control the activation of a specific promoter remains unclear. Our recent work revealed that the TGFβ pathway drives the assembly of enhancer clusters and precise gene activation during neurogenesis. We discovered that the TGFβ pathway coactivator JMJD3 was essential in maintaining these structures in the 3D space. To do that, JMJD3 required an intrinsically disordered region involved in forming phase-separated biomolecular condensates found in the enhancer clusters. Our data support the existence of a relationship between 3D-conformation of the chromatin, biomolecular condensates, and TGFβ-driven response during mammalian neurogenesis. In this review, we discuss how signaling (TGFβ), epigenetics (JMJD3), and biochemical properties (biomolecular condensates nucleation) are coordinated to modulate the genome structure to guarantee proper neural development. Moreover, we comment on the potential underlying mechanisms and implications of the enhancer-mediated regulation. Finally, we point out the knowledge gaps that still need to be addressed.

Surfactancy in a tadpole model of proteins

We model the environment of eukaryotic nuclei by representing macromolecules by only their entropic properties, with globular molecules represented by spherical colloids and flexible molecules by polymers. We put particular focus on proteins with both globular and intrinsically disordered regions, which we represent with 'tadpole' constructed by grafting single polymers and colloids together. In Monte Carlo simulations, we find these tadpoles support phase separation via depletion flocculation, and demonstrate several surfactant behaviours, including being found preferentially at interfaces and forming micelles in single phase solution. Furthermore, the model parameters can be tuned to give a tadpole a preference for either bulk phase. However, we find entropy too weak to drive these behaviours by itself at likely biological concentrations.

The transcriptional coactivator RUVBL2 regulates Pol II clustering with diverse transcription factors

RNA polymerase II (Pol II) apparatuses are compartmentalized into transcriptional clusters. Whether protein factors control these clusters remains unknown. In this study, we find that the ATPase-associated with diverse cellular activities (AAA + ) ATPase RUVBL2 co-occupies promoters with Pol II and various transcription factors. RUVBL2 interacts with unphosphorylated Pol II in chromatin to promote RPB1 carboxy-terminal domain (CTD) clustering and transcription initiation. Rapid depletion of RUVBL2 leads to a decrease in the number of Pol II clusters and inhibits nascent RNA synthesis, and tethering RUVBL2 to an active promoter enhances Pol II clustering at the promoter. We also identify target genes that are directly linked to the RUVBL2-Pol II axis. Many of these genes are hallmarks of cancers and encode proteins with diverse cellular functions. Our results demonstrate an emerging activity for RUVBL2 in regulating Pol II cluster formation in the nucleus.

Deciphering multi-way interactions in the human genome

Chromatin architecture, a key regulator of gene expression, can be inferred using chromatin contact data from chromosome conformation capture, or Hi-C. However, classical Hi-C does not preserve multi-way contacts. Here we use long sequencing reads to map genome-wide multi-way contacts and investigate higher order chromatin organization in the human genome. We use hypergraph theory for data representation and analysis, and quantify higher order structures in neonatal fibroblasts, biopsied adult fibroblasts, and B lymphocytes. By integrating multi-way contacts with chromatin accessibility, gene expression, and transcription factor binding, we introduce a data-driven method to identify cell type-specific transcription clusters. We provide transcription factor-mediated functional building blocks for cell identity that serve as a global signature for cell types.

Mapping higher order chromatin architecture is important. Here the authors use long sequencing reads to map genome-wide multi-way contacts and investigate higher order chromatin organisation; they use hypergraph theory for data representation and analysis, and apply this to different cell types.

Polymer architecture orchestrates the segregation and spatial organization of replicating E. coli chromosomes in slow growth

The mechanism of chromosome segregation and organization in the bacterial cell cycle of E. coli is one of the least understood aspects in its life cycle. The E. coli chromosome is often modelled as a bead spring ring polymer. We introduce cross-links in the DNA-ring polymer, resulting in the formation of loops within each replicating bacterial chromosome. We use simulations to show that the chosen polymer-topology ensures its self-organization along the cell long-axis, such that various chromosomal loci get spatially localized as seen in vivo. The localization of loci arises due to entropic repulsion between polymer loops within each daughter DNA confined in a cylinder. The cellular addresses of the loci in our model are in fair agreement with those seen in experiments as given in J. A. Cass et al., Biophys. J., 2016, 110, 2597-2609. We also show that the adoption of such modified polymer architectures by the daughter DNAs leads to an enhanced propensity of their spatial segregation. Secondly, we match other experimentally reported results, including observation of the cohesion time and the ter-transition. Additionally, the contact map generated from our simulations reproduces the macro-domain like organization as seen in the experimentally obtained Hi-C map. Lastly, we have also proposed a plausible reconciliation of the 'Train Track' and the 'Replication Factory' models which provide conflicting descriptions of the spatial organization of the replication forks. Thus, we reconcile observations from complementary experimental techniques probing bacterial chromosome organization.

Domain Model of Eukaryotic Genome Organization: From DNA Loops Fixed on the Nuclear Matrix to TADs

The article reviews the development of ideas on the domain organization of eukaryotic genome, with special attention on the studies of DNA loops anchored to the nuclear matrix and their role in the emergence of the modern model of eukaryotic genome spatial organization. Critical analysis of results demonstrating that topologically associated chromatin domains are structural-functional blocks of the genome supports the notion that these blocks are fundamentally different from domains whose existence was proposed by the domain hypothesis of eukaryotic genome organization formulated in the 1980s. Based on the discussed evidence, it is concluded that the model postulating that eukaryotic genome is built from uniformly organized structural-functional blocks has proven to be untenable.

Bacterial H-NS contacts DNA at the same irregularly spaced sites in both bridged and hemi-sequestered linear filaments

Gene silencing in bacteria is mediated by chromatin proteins, of which Escherichia coli H-NS is a paradigmatic example. H-NS forms nucleoprotein filaments with either one or two DNA duplexes. However, the structures, arrangements of DNA-binding domains (DBDs), and positions of DBD-DNA contacts in linear and bridged filaments are uncertain. To characterize the H-NS DBD contacts that silence transcription by RNA polymerase, we combined ·OH footprinting, molecular dynamics, statistical modeling, and DBD mapping using a chemical nuclease (Fe²⁺-EDTA) tethered to the DBDs (TEN-map). We find that H-NS DBDs contact DNA at indistinguishable locations in bridged or linear filaments and that the DBDs vary in orientation and position with ∼10-bp average spacing. Our results support a hemi-sequestration model of linear-to-bridged H-NS switching. Linear filaments able to inhibit only transcription initiation switch to bridged filaments able to inhibit both initiation and elongation using the same irregularly spaced DNA contacts.

How to tame your genes: mechanisms of inflammatory gene repression by glucocorticoids

Glucocorticoids (GCs) are widely used therapeutic agents to treat a broad range of inflammatory conditions. Their functional effects are elicited by binding to the glucocorticoid receptor (GR), which regulates transcription of distinct gene networks in response to ligand. However, the mechanisms governing various aspects of undesired side effects versus beneficial immunomodulation upon GR activation remain complex and incompletely understood. In this review, we discuss emerging models of inflammatory gene regulation by GR, highlighting GR's regulatory specificity conferred by context-dependent changes in chromatin architecture and transcription factor or co-regulator dynamics. GR controls both gene activation and repression, with the repression mechanism being central to favorable clinical outcomes. We describe current knowledge about 3D genome organization and its role in spatiotemporal transcriptional control by GR. Looking beyond, we summarize the evidence for dynamics in gene regulation by GR through cooperative convergence of epigenetic modifications, transcription factor crosstalk, molecular condensate formation and chromatin looping. Further characterizing these genomic events will reframe our understanding of mechanisms of transcriptional repression by GR.

Wheat Breeding, Transcription Factories, and Genetic Interactions: New Perspectives

Epistatic interactions and negative heterosis have been shown to be associated with interchromosomal interactions in wheat. Physical gene-gene interactions between co-regulated genes clustered in "transcription factories" have been documented, and a genome-wide atlas of functionally paired, interacting regulatory elements and genes of wheat recently produced. Integration of these new studies on gene and regulatory element interactions, co-regulation of gene expression in "transcription factories," and epigenetics generates new perspectives for wheat breeding and trait enhancement.

Functional organization of RNA polymerase II in nuclear subcompartments

Distinct clusters of RNA polymerase II are responsible for gene transcription inside eukaryotic cell nuclei. Despite the functional implications of such subnuclear organization, the attributes of these clusters and the mechanisms underlying their formation remain only partially understood. Recently, the concept of proteins and RNA phase-separating into liquid-like droplets was proposed to drive the formation of transcriptionally-active subcompartments. Here, we attempt to reconcile previous with more recent findings, and discuss how the different ways of assembling the active RNA polymerase II transcriptional machinery relate to nuclear compartmentalization.

Liquid-Liquid Phase Separation in Chromatin

In eukaryotic cells, protein and RNA factors involved in genome activities like transcription, RNA processing, DNA replication, and repair accumulate in self-organizing membraneless chromatin subcompartments. These structures contribute to efficiently conduct chromatin-mediated reactions and to establish specific cellular programs. However, the underlying mechanisms for their formation are only partly understood. Recent studies invoke liquid-liquid phase separation (LLPS) of proteins and RNAs in the establishment of chromatin activity patterns. At the same time, the folding of chromatin in the nucleus can drive genome partitioning into spatially distinct domains. Here, the interplay between chromatin organization, chromatin binding, and LLPS is discussed by comparing and contrasting three prototypical chromatin subcompartments: the nucleolus, clusters of active RNA polymerase II, and pericentric heterochromatin domains. It is discussed how the different ways of chromatin compartmentalization are linked to transcription regulation, the targeting of soluble factors to certain parts of the genome, and to disease-causing genetic aberrations.

Merging Established Mechanisms with New Insights: Condensates, Hubs, and the Regulation of RNA Polymerase II Transcription

The regulation of RNA polymerase II (pol II) transcription requires a complex and context-specific array of proteins and protein complexes, as well as nucleic acids and metabolites. Every major physiological process requires coordinated transcription of specific sets of genes at the appropriate time, and a breakdown in this regulation is a hallmark of human disease. A proliferation of recent studies has revealed that many general transcription components, including sequence-specific, DNA-binding transcription factors, Mediator, and pol II itself, are capable of liquid-liquid phase separation, to form condensates that partition these factors away from the bulk aqueous phase. These findings hold great promise for next-level understanding of pol II transcription; however, many mechanistic aspects align with more conventional models, and whether phase separation per se regulates pol II activity in cells remains controversial. In this review, we describe the conventional and condensate-dependent models, and why their similarities and differences are important. We also compare and contrast these models in the context of genome organization and pol II transcription (initiation, elongation, and termination), and highlight the central role of RNA in these processes. Finally, we discuss mutations that disrupt normal partitioning of transcription factors, and how this may contribute to disease.

A tandem duplication in Drosophila melanogaster shows enhanced expression beyond the gene copy number

Tandem duplicated genes are common features of genomes, but the phenotypic consequences of their origins are not well understood. It is not known whether a simple doubling of gene expression should be expected, or else some other expression outcome. This study describes an experimental framework using engineered deletions to assess any contribution of locally acting cis- and globally acting trans-regulatory factors to expression interactions of particular tandem duplicated genes. Acsx1L (CG6300) and Acsx1R (CG11659) are tandem duplicates of a putative acyl-CoA synthetase gene found in Drosophila melanogaster. Experimental deletions of the duplicated segments were used to investigate whether the presence of 1 tandem duplicated block influences the expression of its neighbor. Acsx1L, the gene in the left block, shows much higher expression than either its duplicate Acsx1R or the single Acsx1 in Drosophila simulans. Acsx1L expression decreases drastically upon deleting the right-hand duplicated block. Crosses among wildtype and deletion strains show that high tandem expression is primarily due to cis-acting interactions between the duplicated blocks. No effect of these genes on cuticular hydrocarbons was detected. Sequence and phylogenetic analysis suggest that the duplication rose to fixation in D. melanogaster and has been subject to extensive gene conversion. Some strains actually carry 3 tandem copies, yet strains with 3 Acsx1s do not have higher expression levels than strains with 2. Surveys of tandem duplicate expression have typically not found the expected 2-fold increase in expression. This study suggests that cis-regulatory interactions between duplicated blocks could be responsible for this trend.

A Shift in Paradigms: Spatial Genomics Approaches to Reveal Single-Cell Principles of Genome Organization

The genome tridimensional (3D) organization and its role towards the regulation of key cell processes such as transcription is currently a main question in biology. Interphase chromosomes are spatially segregated into "territories," epigenetically-defined large domains of chromatin that interact to form "compartments" with common transcriptional status, and insulator-flanked domains called "topologically associating domains" (TADs). Moreover, chromatin organizes around nuclear structures such as lamina, speckles, or the nucleolus to acquire a higher-order genome organization. Due to recent technological advances, the different hierarchies are being solved. Particularly, advances in microscopy technologies are shedding light on the genome structure at multiple levels. Intriguingly, more and more reports point to high variability and stochasticity at the single-cell level. However, the functional consequences of such variability in genome conformation are still unsolved. Here, I will discuss the implication of the cell-to-cell heterogeneity at the different scales in the context of newly developed imaging approaches, particularly multiplexed Fluorescence in situ hybridization methods that enabled "chromatin tracing." Extensions of these methods are now combining spatial information of dozens to thousands of genomic loci with the localization of nuclear features such as the nucleolus, nuclear speckles, or even histone modifications, creating the fast-moving field of "spatial genomics." As our view of genome organization shifts the focus from ensemble to single-cell, new insights to fundamental questions begin to emerge.

4DNvestigator: time series genomic data analysis toolbox

Data on genome organization and output over time, or the 4D Nucleome (4DN), require synthesis for meaningful interpretation. Development of tools for the efficient integration of these data is needed, especially for the time dimension. We present the '4DNvestigator', a user-friendly network-based toolbox for the analysis of time series genome-wide genome structure (Hi-C) and gene expression (RNA-seq) data. Additionally, we provide methods to quantify network entropy, tensor entropy, and statistically significant changes in time series Hi-C data at different genomic scales.

Predicting genome organisation and function with mechanistic modelling

Fitting-free mechanistic models based on polymer simulations predict chromatin folding in 3D by focussing on the underlying biophysical mechanisms. This class of models has been increasingly used in conjunction with experiments to study the spatial organisation of eukaryotic chromosomes. Feedback from experiments to models leads to successive model refinement and has previously led to the discovery of new principles for genome organisation. Here, we review the basis of mechanistic polymer simulations, explain some of the more recent approaches and the contexts in which they have been useful to explain chromosome biology, and speculate on how they might be used in the future.

Spatial organization of transcribing loci during early genome activation in Drosophila

The early Drosophila embryo provides unique experimental advantages for addressing fundamental questions of gene regulation at multiple levels of organization, from individual gene loci to the entire genome. Using 1.5-h-old Drosophila embryos undergoing the first wave of genome activation,¹ we detected ∼110 discrete "speckles" of RNA polymerase II (RNA Pol II) per nucleus, two of which were larger and localized to the histone locus bodies (HLBs).^2,3 In the absence of the primary driver of Drosophila genome activation, the pioneer factor Zelda (Zld),^1,4,5 70% fewer speckles were present; however, the HLBs tended to be larger than wild-type (WT) HLBs, indicating that RNA Pol II accumulates at the HLBs in the absence of robust early-gene transcription. We observed a uniform distribution of distances between active genes in the nuclei of both WT and zld mutant embryos, indicating that early co-regulated genes do not cluster into nuclear sub-domains. However, in instances whereby transcribing genes did come into close 3D proximity (within 400 nm), they were found to have distinct RNA Pol II speckles. In contrast to the emerging model whereby active genes are clustered to facilitate co-regulation and sharing of transcriptional resources, our data support an "individualist" model of gene control at early genome activation in Drosophila. This model is in contrast to a "collectivist" model, where active genes are spatially clustered and share transcriptional resources, motivating rigorous tests of both models in other experimental systems.

Epromoters function as a hub to recruit key transcription factors required for the inflammatory response

Gene expression is controlled by the involvement of gene-proximal (promoters) and distal (enhancers) regulatory elements. Our previous results demonstrated that a subset of gene promoters, termed Epromoters, work as bona fide enhancers and regulate distal gene expression. Here, we hypothesized that Epromoters play a key role in the coordination of rapid gene induction during the inflammatory response. Using a high-throughput reporter assay we explored the function of Epromoters in response to type I interferon. We find that clusters of IFNa-induced genes are frequently associated with Epromoters and that these regulatory elements preferentially recruit the STAT1/2 and IRF transcription factors and distally regulate the activation of interferon-response genes. Consistently, we identified and validated the involvement of Epromoter-containing clusters in the regulation of LPS-stimulated macrophages. Our findings suggest that Epromoters function as a local hub recruiting the key TFs required for coordinated regulation of gene clusters during the inflammatory response.

Isolation of Nucleoid Fraction from Mycoplasma gallisepticum Cells with Synchronized Cell Division

It is assumed that unknown mechanisms can be involved in adaptation Mycoplasma gallisepticum to unfavorable factors, one of these can be local rearrangements of the structure and spatial organization of the chromosome. To study these mechanisms, we obtained a culture of M. gallisepticum with synchronized division and isolated the nucleoid fraction from this culture by the method of mild cell lysis and centrifugation in a sucrose gradient. Liquid chromatography-mass spectrometry analysis of the proteome showed that in comparison with the cell lysate, the nucleoid fraction was enriched with DNA-binding proteins. This analysis will help to find new nucleoid-associated proteins and to study their dynamics, distribution, and their role during infection and under stress conditions.

Suppression of liquid-liquid phase separation by 1,6-hexanediol partially compromises the 3D genome organization in living cells

Liquid-liquid phase separation (LLPS) contributes to the spatial and functional segregation of molecular processes within the cell nucleus. However, the role played by LLPS in chromatin folding in living cells remains unclear. Here, using stochastic optical reconstruction microscopy (STORM) and Hi-C techniques, we studied the effects of 1,6-hexanediol (1,6-HD)-mediated LLPS disruption/modulation on higher-order chromatin organization in living cells. We found that 1,6-HD treatment caused the enlargement of nucleosome clutches and their more uniform distribution in the nuclear space. At a megabase-scale, chromatin underwent moderate but irreversible perturbations that resulted in the partial mixing of A and B compartments. The removal of 1,6-HD from the culture medium did not allow chromatin to acquire initial configurations, and resulted in more compact repressed chromatin than in untreated cells. 1,6-HD treatment also weakened enhancer-promoter interactions and TAD insulation but did not considerably affect CTCF-dependent loops. Our results suggest that 1,6-HD-sensitive LLPS plays a limited role in chromatin spatial organization by constraining its folding patterns and facilitating compartmentalization at different levels.

Complex small-world regulatory networks emerge from the 3D organisation of the human genome

The discovery that overexpressing one or a few critical transcription factors can switch cell state suggests that gene regulatory networks are relatively simple. In contrast, genome-wide association studies (GWAS) point to complex phenotypes being determined by hundreds of loci that rarely encode transcription factors and which individually have small effects. Here, we use computer simulations and a simple fitting-free polymer model of chromosomes to show that spatial correlations arising from 3D genome organisation naturally lead to stochastic and bursty transcription as well as complex small-world regulatory networks (where the transcriptional activity of each genomic region subtly affects almost all others). These effects require factors to be present at sub-saturating levels; increasing levels dramatically simplifies networks as more transcription units are pressed into use. Consequently, results from GWAS can be reconciled with those involving overexpression. We apply this pan-genomic model to predict patterns of transcriptional activity in whole human chromosomes, and, as an example, the effects of the deletion causing the diGeorge syndrome.

The Role of Human Satellite III (1q12) Copy Number Variation in the Adaptive Response during Aging, Stress, and Pathology: A Pendulum Model

The pericentric satellite III (SatIII or Sat3) and II tandem repeats recently appeared to be transcribed under stress conditions, and the transcripts were shown to play an essential role in the universal stress response. In this paper, we review the role of human-specific SatIII copy number variation (CNV) in normal stress response, aging and pathology, with a focus on 1q12 loci. We postulate a close link between transcription of SatII/III repeats and their CNV. The accrued body of data suggests a hypothetical universal mechanism, which provides for SatIII copy gain during the stress response, alongside with another, more hypothetical reverse mechanism that might reduce the mean SatIII copy number, likely via the selection of cells with excessively large 1q12 loci. Both mechanisms, working alternatively like swings of the pendulum, may ensure the balance of SatIII copy numbers and optimum stress resistance. This model is verified on the most recent data on SatIII CNV in pathology and therapy, aging, senescence and response to genotoxic stress in vitro.

Coordination of transcription, processing, and export of highly expressed RNAs by distinct biomolecular condensates

Genes under control of super-enhancers are expressed at extremely high levels and are frequently associated with nuclear speckles. Recent data suggest that the high concentration of unphosphorylated RNA polymerase II (Pol II) and Mediator recruited to super-enhancers create phase-separated condensates. Transcription initiates within or at the surface of these phase-separated droplets and the phosphorylation of Pol II, associated with transcription initiation and elongation, dissociates Pol II from these domains leading to engagement with nuclear speckles, which are enriched with RNA processing factors. The transitioning of Pol II from transcription initiation domains to RNA processing domains effectively co-ordinates transcription and processing of highly expressed RNAs which are then rapidly exported into the cytoplasm.

Three-dimensional chromosome organization in flowering plants

Research on plant three-dimensional (3D) genome architecture made rapid progress over the past 5 years. Numerous Hi-C interaction data sets were generated in a wide range of plant species, allowing for a comprehensive overview on 3D chromosome folding principles in the plant kingdom. Plants lack important genes reported to be vital for chromosome folding in animals. However, similar 3D structures such as topologically associating domains and chromatin loops were identified. Recent studies in Arabidopsis thaliana revealed how chromosomal regions are positioned within the nucleus by determining their association with both, the nuclear periphery and the nucleolus. Additionally, many plant species exhibit high-frequency interactions among KNOT entangled elements, which are associated with safeguarding the genome from invasive DNA elements. Many of the recently published Hi-C data sets were generated to aid de novo genome assembly and remain to date little explored. These data sets represent a valuable resource for future comparative studies, which may lead to a more profound understanding of the evolution of 3D chromosome organization in plants.

Collaboration through chromatin: motors of transcription and chromatin structure

Packaging of the eukaryotic genome into chromatin places fundamental physical constraints on transcription. Clarifying how transcription operates within these constraints is essential to understand how eukaryotic gene expression programs are established and maintained. Here we review what is known about the mechanisms of transcription on chromatin templates. Current models indicate that transcription through chromatin is accomplished by the combination of an inherent nucleosome disrupting activity of RNA polymerase and the action of ATP-dependent chromatin remodeling motors. Collaboration between these two types of molecular motors is proposed to occur at all stages of transcription through diverse mechanisms. Further investigation of how these two motors combine their basic activities is essential to clarify the interdependent relationship between genome structure and transcription.

Models behind the mystery of establishing enhancer-promoter interactions

Enhancers and promoters are transcriptional regulatory elements whose facilitated interactions increase gene expression. Enhancer DNA sequences can be located far away from the promoter sequences that they regulate. Currently, the mechanism facilitating the establishment of enhancer-promoter interactions remains unclear. However, mutations causing errors in these interactions have been linked to cancer and disease, further conveying the need to understand the full mechanism. This review discusses multiple models that have been proposed to describe how enhancers go the distance to interact with promoters. Evidence supporting loop formation models is reviewed in addition to more complex hypotheses involving aspects of 3D chromatin organization and phase separation.

Non-coding RNAs in chromatin folding and nuclear organization

One of the most intriguing questions facing modern biology concerns how the genome directs the construction of cells, tissues, and whole organisms. It is tempting to suggest that the part of the genome that does not encode proteins contains architectural plans. We are still far from understanding how these plans work at the level of building tissues and the body as a whole. However, the results of recent studies demonstrate that at the cellular level, special non-coding RNAs serve as scaffolds for the construction of various intracellular structures. The term "architectural RNAs" was proposed to designate this subset of non-coding RNAs. In this review, we discuss the role of architectural RNAs in the construction of the cell nucleus and maintenance of the three-dimensional organization of the genome.

New insights into promoter-enhancer communication mechanisms revealed by dynamic single-molecule imaging

Establishing cell-type-specific gene expression programs relies on the action of distal enhancers, cis-regulatory elements that can activate target genes over large genomic distances - up to Mega-bases away. How distal enhancers physically relay regulatory information to target promoters has remained a mystery. Here, we review the latest developments and insights into promoter-enhancer communication mechanisms revealed by live-cell, real-time single-molecule imaging approaches.

Computational methods for predicting 3D genomic organization from high-resolution chromosome conformation capture data

The advent of high-resolution chromosome conformation capture assays (such as 5C, Hi-C and Pore-C) has allowed for unprecedented sequence-level investigations into the structure-function relationship of the genome. In order to comprehensively understand this relationship, computational tools are required that utilize data generated from these assays to predict 3D genome organization (the 3D genome reconstruction problem). Many computational tools have been developed that answer this need, but a comprehensive comparison of their underlying algorithmic approaches has not been conducted. This manuscript provides a comprehensive review of the existing computational tools (from November 2006 to September 2019, inclusive) that can be used to predict 3D genome organizations from high-resolution chromosome conformation capture data. Overall, existing tools were found to use a relatively small set of algorithms from one or more of the following categories: dimensionality reduction, graph/network theory, maximum likelihood estimation (MLE) and statistical modeling. Solutions in each category are far from maturity, and the breadth and depth of various algorithmic categories have not been fully explored. While the tools for predicting 3D structure for a genomic region or single chromosome are diverse, there is a general lack of algorithmic diversity among computational tools for predicting the complete 3D genome organization from high-resolution chromosome conformation capture data.

Homotypic clustering of L1 and B1/Alu repeats compartmentalizes the 3D genome

Organization of the genome into euchromatin and heterochromatin appears to be evolutionarily conserved and relatively stable during lineage differentiation. In an effort to unravel the basic principle underlying genome folding, here we focus on the genome itself and report a fundamental role for L1 (LINE1 or LINE-1) and B1/Alu retrotransposons, the most abundant subclasses of repetitive sequences, in chromatin compartmentalization. We find that homotypic clustering of L1 and B1/Alu demarcates the genome into grossly exclusive domains, and characterizes and predicts Hi-C compartments. Spatial segregation of L1-rich sequences in the nuclear and nucleolar peripheries and B1/Alu-rich sequences in the nuclear interior is conserved in mouse and human cells and occurs dynamically during the cell cycle. In addition, de novo establishment of L1 and B1 nuclear segregation is coincident with the formation of higher-order chromatin structures during early embryogenesis and appears to be critically regulated by L1 and B1 transcripts. Importantly, depletion of L1 transcripts in embryonic stem cells drastically weakens homotypic repeat contacts and compartmental strength, and disrupts the nuclear segregation of L1- or B1-rich chromosomal sequences at genome-wide and individual sites. Mechanistically, nuclear co-localization and liquid droplet formation of L1 repeat DNA and RNA with heterochromatin protein HP1α suggest a phase-separation mechanism by which L1 promotes heterochromatin compartmentalization. Taken together, we propose a genetically encoded model in which L1 and B1/Alu repeats blueprint chromatin macrostructure. Our model explains the robustness of genome folding into a common conserved core, on which dynamic gene regulation is overlaid across cells.

Principles of 3D compartmentalization of the human genome

Chromatin is organized in the nucleus via CTCF loops and compartmental domains. Here, we compare different cell types to identify distinct paradigms of compartmental domain formation in human tissues. We identify and quantify compartmental forces correlated with histone modifications characteristic of transcriptional activity and previously underappreciated roles for distinct compartmental domains correlated with the presence of H3K27me3 and H3K9me3, respectively. We present a computer simulation model capable of predicting compartmental organization based on the biochemical characteristics of independent chromatin features. Using this model, we show that the underlying forces responsible for compartmental domain formation in human cells are conserved and that the diverse compartmentalization patterns seen across cell types are due to differences in chromatin features. We extend these findings to Drosophila to suggest that the same principles are at work beyond humans. These results offer mechanistic insights into the fundamental forces driving the 3D organization of the genome.

Using high-resolution Hi-C data and computer simulations, Nichols and Corces show that compartments arise as a consequence of interactions among proteins that correlate with the presence of H3K27ac, H3K27me3, and H3K9me3, suggesting that human cells contain at least three distinct compartments. The same principles apply to other organisms.

Random sub-diffusion and capture of genes by the nuclear pore reduces dynamics and coordinates inter-chromosomal movement

Hundreds of genes interact with the yeast nuclear pore complex (NPC), localizing at the nuclear periphery and clustering with co-regulated genes. Dynamic tracking of peripheral genes shows that they cycle on and off the NPC and that interaction with the NPC slows their sub-diffusive movement. Furthermore, NPC-dependent inter-chromosomal clustering leads to coordinated movement of pairs of loci separated by hundreds of nanometers. We developed fractional Brownian motion simulations for chromosomal loci in the nucleoplasm and interacting with NPCs. These simulations predict the rate and nature of random sub-diffusion during repositioning from nucleoplasm to periphery and match measurements from two different experimental models, arguing that recruitment to the nuclear periphery is due to random sub-diffusion and transient capture by NPCs. Finally, the simulations do not lead to inter-chromosomal clustering or coordinated movement, suggesting that interaction with the NPC is necessary, but not sufficient, to cause clustering.

Specific Features of RNA Polymerases I and III: Structure and Assembly

RNA polymerase I (RNAPI) and RNAPIII are multi-heterogenic protein complexes that specialize in the transcription of highly abundant non-coding RNAs, such as ribosomal RNA (rRNA) and transfer RNA (tRNA). In terms of subunit number and structure, RNAPI and RNAPIII are more complex than RNAPII that synthesizes thousands of different mRNAs. Specific subunits of the yeast RNAPI and RNAPIII form associated subcomplexes that are related to parts of the RNAPII initiation factors. Prior to their delivery to the nucleus where they function, RNAP complexes are assembled at least partially in the cytoplasm. Yeast RNAPI and RNAPIII share heterodimer Rpc40-Rpc19, a functional equivalent to the αα homodimer which initiates assembly of prokaryotic RNAP. In the process of yeast RNAPI and RNAPIII biogenesis, Rpc40 and Rpc19 form the assembly platform together with two small, bona fide eukaryotic subunits, Rpb10 and Rpb12. We propose that this assembly platform is co-translationally seeded while the Rpb10 subunit is synthesized by cytoplasmic ribosome machinery. The translation of Rpb10 is stimulated by Rbs1 protein, which binds to the 3′-untranslated region of RPB10 mRNA and hypothetically brings together Rpc19 and Rpc40 subunits to form the αα-like heterodimer. We suggest that such a co-translational mechanism is involved in the assembly of RNAPI and RNAPIII complexes.

Coordination of transcription and processing of tRNA

Coordination of transcription and processing of RNA is a basic principle in regulation of gene expression in eukaryotes. In the case of mRNA, coordination is primarily founded on a co-transcriptional processing mechanism by which a nascent precursor mRNA undergoes maturation via cleavage and modification by the transcription machinery. A similar mechanism controls the biosynthesis of rRNA. However, the coordination of transcription and processing of tRNA, a rather short transcript, remains unknown. Here, we present a model for high molecular weight initiation complexes of human RNA polymerase III that assemble on tRNA genes and process precursor transcripts to mature forms. These multifunctional initiation complexes may support co-transcriptional processing, such as the removal of the 5' leader of precursor tRNA by RNase P. Based on this model, maturation of tRNA is predetermined prior to transcription initiation.

Repression of transcription by the glucocorticoid receptor: A parsimonious model for the genomics era

Glucocorticoids are potent anti-inflammatory drugs that are used to treat an extraordinary range of human disease, including COVID-19, underscoring the ongoing importance of understanding their molecular mechanisms. Early studies of GR signaling led to broad acceptance of models in which glucocorticoid receptor (GR) monomers tether repressively to inflammatory transcription factors, thus abrogating inflammatory gene expression. However, newer data challenge this core concept and present an exciting opportunity to reframe our understanding of GR signaling. Here, we present an alternate, two-part model for transcriptional repression by glucocorticoids. First, widespread GR-mediated induction of transcription results in rapid, primary repression of inflammatory gene transcription and associated enhancers through competition-based mechanisms. Second, a subset of GR-induced genes, including targets that are regulated in coordination with inflammatory transcription factors such as NF-κB, exerts secondary repressive effects on inflammatory gene expression. Within this framework, emerging data indicate that the gene set regulated through the cooperative convergence of GR and NF-κB signaling is central to the broad clinical effectiveness of glucocorticoids in terminating inflammation and promoting tissue repair.