Insulator function and topological domain border strength scale with architectural protein occupancy

Characterization of Tau95 led to the identification of a four-subunit TFIIIC complex in trypanosomatid parasites

RNA polymerase III (RNAP III) synthetizes small essential non-coding RNA molecules such as tRNAs and 5S rRNA. In yeast and vertebrates, RNAP III needs general transcription factors TFIIIA, TFIIIB, and TFIIIC to initiate transcription. TFIIIC, composed of six subunits, binds to internal promoter elements in RNAP III-dependent genes. Limited information is available about RNAP III transcription in the trypanosomatid protozoa Trypanosoma brucei and Leishmania major, which diverged early from the eukaryotic lineage. Analyses of the first published draft of the trypanosomatid genome sequences failed to recognize orthologs of any of the TFIIIC subunits, suggesting that this transcription factor is absent in these parasites. However, a putative TFIIIC subunit was recently annotated in the databases. Here we characterize this subunit in T. brucei and L. major and demonstrate that it corresponds to Tau95. In silico analyses showed that both proteins possess the typical Tau95 sequences: the DNA binding region and the dimerization domain. As anticipated for a transcription factor, Tau95 localized to the nucleus in insect forms of both parasites. Chromatin immunoprecipitation (ChIP) assays demonstrated that Tau95 binds to tRNA and U2 snRNA genes in T. brucei. Remarkably, by performing tandem affinity purifications we identified orthologs of TFIIIC subunits Tau55, Tau131, and Tau138 in T. brucei and L. major. Thus, contrary to what was assumed, trypanosomatid parasites do possess a TFIIIC complex. Other putative interacting partners of Tau95 were identified in T. brucei and L. major. KEY POINTS: • A four-subunit TFIIIC complex is present in T. brucei and L. major • TbTau95 associates with tRNA and U2 snRNA genes • Putative interacting partners of Tau95 might include some RNAP II regulators.

Chromatin insulator mechanisms ensure accurate gene expression by controlling overall 3D genome organization

Chromatin insulators are DNA-protein complexes that promote specificity of enhancer-promoter interactions and maintain distinct transcriptional states through control of 3D genome organization. In this review, we highlight recent work visualizing how mammalian CCCTC-binding factor acts as a boundary to dynamic DNA loop extrusion mediated by cohesin. We also discuss new studies in both mammals and Drosophila that elucidate biological redundancy of chromatin insulator function and interplay with transcription with respect to topologically associating domain formation. Finally, we present novel concepts in spatiotemporal regulation of chromatin insulator function during differentiation and development and possible consequences of disrupted insulator activity on cellular proliferation.

DNA choreography: correlating mobility and organization of DNA across different resolutions from loops to chromosomes

The dynamics of DNA in the cell nucleus plays a role in cellular processes and fates but the interplay of DNA mobility with the hierarchical levels of DNA organization is still underexplored. Here, we made use of DNA replication to directly label genomic DNA in an unbiased genome-wide manner. This was followed by live-cell time-lapse microscopy of the labeled DNA combining imaging at different resolutions levels simultaneously and allowing one to trace DNA motion across organization levels within the same cells. Quantification of the labeled DNA segments at different microscopic resolution levels revealed sizes comparable to the ones reported for DNA loops using 3D super-resolution microscopy, topologically associated domains (TAD) using 3D widefield microscopy, and also entire chromosomes. By employing advanced chromatin tracking and image registration, we discovered that DNA exhibited higher mobility at the individual loop level compared to the TAD level and even less at the chromosome level. Additionally, our findings indicate that chromatin movement, regardless of the resolution, slowed down during the S phase of the cell cycle compared to the G1/G2 phases. Furthermore, we found that a fraction of DNA loops and TADs exhibited directed movement with the majority depicting constrained movement. Our data also indicated spatial mobility differences with DNA loops and TADs at the nuclear periphery and the nuclear interior exhibiting lower velocity and radius of gyration than the intermediate locations. On the basis of these insights, we propose that there is a link between DNA mobility and its organizational structure including spatial distribution, which impacts cellular processes.

Genome organization regulates nuclear pore complex formation and promotes differentiation during Drosophila oogenesis

Genome organization can regulate gene expression and promote cell fate transitions. The differentiation of germline stem cells (GSCs) to oocytes in Drosophila involves changes in genome organization mediated by heterochromatin and the nuclear pore complex (NPC). Heterochromatin represses germ cell genes during differentiation, and NPCs anchor these silenced genes to the nuclear periphery, maintaining silencing to allow for oocyte development. Surprisingly, we found that genome organization also contributes to NPC formation, mediated by the transcription factor Stonewall (Stwl). As GSCs differentiate, Stwl accumulates at boundaries between silenced and active gene compartments. Stwl at these boundaries plays a pivotal role in transitioning germ cell genes into a silenced state and activating a group of oocyte genes and nucleoporins (Nups). The upregulation of these Nups during differentiation is crucial for NPC formation and further genome organization. Thus, cross-talk between genome architecture and NPCs is essential for successful cell fate transitions.

The zinc-finger protein Z4 cooperates with condensin II to regulate somatic chromosome pairing and 3D chromatin organization

Chromosome pairing constitutes an important level of genome organization, yet the mechanisms that regulate pairing in somatic cells and the impact on 3D chromatin organization are still poorly understood. Here, we address these questions in Drosophila, an organism with robust somatic pairing. In Drosophila, pairing preferentially occurs at loci consisting of numerous architectural protein binding sites (APBSs), suggesting a role of architectural proteins (APs) in pairing regulation. Amongst these, the anti-pairing function of the condensin II subunit CAP-H2 is well established. However, the factors that regulate CAP-H2 localization and action at APBSs remain largely unknown. Here, we identify two factors that control CAP-H2 occupancy at APBSs and, therefore, regulate pairing. We show that Z4, interacts with CAP-H2 and is required for its localization at APBSs. We also show that hyperosmotic cellular stress induces fast and reversible unpairing in a Z4/CAP-H2 dependent manner. Moreover, by combining the opposite effects of Z4 depletion and osmostress, we show that pairing correlates with the strength of intrachromosomal 3D interactions, such as active (A) compartment interactions, intragenic gene-loops, and polycomb (Pc)-mediated chromatin loops. Altogether, our results reveal new players in CAP-H2-mediated pairing regulation and the intimate interplay between inter-chromosomal and intra-chromosomal 3D interactions.

Members of an array of zinc finger proteins specify distinct Hox chromatin boundaries

The partitioning of repressive from actively transcribed chromatin domains in mammalian cells fosters cell-type specific gene expression patterns. During differentiation, this partitioning is reconstructed, reflecting gene expression profiles appropriate to new cellular identities. Yet, the chromatin occupancy of the key chromatin insulator, CTCF, at the developmentally important Hox clusters remains unchanged during differentiation. Thus, dynamic changes in chromatin boundaries must entail other activities, such as the previously identified MAZ insulator. Given its requirement for chromatin loop formation, we examined cohesin-based chromatin occupancy without the known insulators, CTCF and MAZ, and identified a novel family of zinc finger proteins (ZNFs), some of which exhibit tissue-specific expression. Thus far, two of these novel ZNFs facilitate the formation of chromatin boundaries at the Hox clusters that are distinct from each other and from that of MAZ. PATZ1 was critical to the thoracolumbar boundary of Hox clusters in differentiating motor neurons in vitro , skeletal development in vivo, and to looping interactions within the genome. On the other hand, ZNF263 contributed to cervicothoracic boundaries in motor neurons. We propose that these novel insulating activities act in concert with cohesin, alone or combinatorially, with or without CTCF, to implement precise positional identity and cell fate during development.

3D epigenomics and 3D epigenopathies

Mammalian genomes are intricately compacted to form sophisticated 3-dimensional structures within the tiny nucleus, so called 3D genome folding. Despite their shapes reminiscent of an entangled yarn, the rapid development of molecular and next-generation sequencing technologies (NGS) has revealed that mammalian genomes are highly organized in a hierarchical order that delicately affects transcription activities. An increasing amount of evidence suggests that 3D genome folding is implicated in diseases, giving us a clue on how to identify novel therapeutic approaches. In this review, we will study what 3D genome folding means in epigenetics, what types of 3D genome structures there are, how they are formed, and how the technologies have developed to explore them. We will also discuss the pathological implications of 3D genome folding. Finally, we will discuss how to leverage 3D genome folding and engineering for future studies. [BMB Reports 2024; 57(5): 216-231].

The scales, mechanisms, and dynamics of the genome architecture

Even when split into several chromosomes, DNA molecules that make up our genome are too long to fit into the cell nuclei unless massively folded. Such folding must accommodate the need for timely access to selected parts of the genome by transcription factors, RNA polymerases, and DNA replication machinery. Here, we review our current understanding of the genome folding inside the interphase nuclei. We consider the resulting genome architecture at three scales with a particular focus on the intermediate (meso) scale and summarize the insights gained from recent experimental observations and diverse computational models.

Enhancer-promoter interactions become more instructive in the transition from cell-fate specification to tissue differentiation

To regulate expression, enhancers must come in proximity to their target gene. However, the relationship between the timing of enhancer-promoter (E-P) proximity and activity remains unclear, with examples of uncoupled, anticorrelated and correlated interactions. To assess this, we selected 600 characterized enhancers or promoters with tissue-specific activity in Drosophila embryos and performed Capture-C in FACS-purified myogenic or neurogenic cells during specification and tissue differentiation. This enabled direct comparison between E-P proximity and activity transitioning from OFF-to-ON and ON-to-OFF states across developmental conditions. This showed remarkably similar E-P topologies between specified muscle and neuronal cells, which are uncoupled from activity. During tissue differentiation, many new distal interactions emerge where changes in E-P proximity reflect changes in activity. The mode of E-P regulation therefore appears to change as embryogenesis proceeds, from largely permissive topologies during cell-fate specification to more instructive regulation during terminal tissue differentiation, when E-P proximity is coupled to activation.

Kdm1a safeguards the topological boundaries of PRC2-repressed genes and prevents aging-related euchromatinization in neurons

Kdm1a is a histone demethylase linked to intellectual disability with essential roles during gastrulation and the terminal differentiation of specialized cell types, including neurons, that remains highly expressed in the adult brain. To explore Kdm1a's function in adult neurons, we develop inducible and forebrain-restricted Kdm1a knockouts. By applying multi-omic transcriptome, epigenome and chromatin conformation data, combined with super-resolution microscopy, we find that Kdm1a elimination causes the neuronal activation of nonneuronal genes that are silenced by the polycomb repressor complex and interspersed with active genes. Functional assays demonstrate that the N-terminus of Kdm1a contains an intrinsically disordered region that is essential to segregate Kdm1a-repressed genes from the neighboring active chromatin environment. Finally, we show that the segregation of Kdm1a-target genes is weakened in neurons during natural aging, underscoring the role of Kdm1a safeguarding neuronal genome organization and gene silencing throughout life.

Machine and deep learning methods for predicting 3D genome organization

Three-Dimensional (3D) chromatin interactions, such as enhancer-promoter interactions (EPIs), loops, Topologically Associating Domains (TADs), and A/B compartments play critical roles in a wide range of cellular processes by regulating gene expression. Recent development of chromatin conformation capture technologies has enabled genome-wide profiling of various 3D structures, even with single cells. However, current catalogs of 3D structures remain incomplete and unreliable due to differences in technology, tools, and low data resolution. Machine learning methods have emerged as an alternative to obtain missing 3D interactions and/or improve resolution. Such methods frequently use genome annotation data (ChIP-seq, DNAse-seq, etc.), DNA sequencing information (k-mers, Transcription Factor Binding Site (TFBS) motifs), and other genomic properties to learn the associations between genomic features and chromatin interactions. In this review, we discuss computational tools for predicting three types of 3D interactions (EPIs, chromatin interactions, TAD boundaries) and analyze their pros and cons. We also point out obstacles of computational prediction of 3D interactions and suggest future research directions.

Enhancer-promoter interactions can form independently of genomic distance and be functional across TAD boundaries

Topologically Associating Domains (TADs) have been suggested to facilitate and constrain enhancer-promoter interactions. However, the role of TAD boundaries in effectively restricting these interactions remains unclear. Here, we show that a significant proportion of enhancer-promoter interactions are established across TAD boundaries in Drosophila embryos, but that developmental genes are strikingly enriched in intra- but not inter-TAD interactions. We pursued this observation using the twist locus, a master regulator of mesoderm development, and systematically relocated one of its enhancers to various genomic locations. While this developmental gene can establish inter-TAD interactions with its enhancer, the functionality of these interactions remains limited, highlighting the existence of topological constraints. Furthermore, contrary to intra-TAD interactions, the formation of inter-TAD enhancer-promoter interactions is not solely driven by genomic distance, with distal interactions sometimes favored over proximal ones. These observations suggest that other general mechanisms must exist to establish and maintain specific enhancer-promoter interactions across large distances.

Effective modeling of the chromatin structure by coarse-grained methods

The interphase chromatin structure is extremely complex, precise and dynamic. Experimental methods can only show the frequency of interaction of the various parts of the chromatin. Therefore, it is extremely important to develop theoretical methods to predict the chromatin structure. In this publication, we implemented an extended version of the SBS model described by Barbieri et al. and created the ChroMC program that is easy to use and freely available (https://github.com/regulomics/chroMC) to other users. We also describe the necessary factors for the effective modeling of the chromatin structure in Drosophila melanogaster. We compared results of chromatin structure predictions using two methods: Monte Carlo and Molecular Dynamic. Our simulations suggest that incorporating black, non-reactive chromatin is necessary for successful prediction of chromatin structure, while the loop extrusion model with a long range attraction potential or Lennard-Jones (with local attraction force) as well as using Hi-C data as input are not essential for the basic structure reconstruction. We also proposed a new way to calculate the similarity of the properties of contact maps including the calculation of local similarity.Communicated by Ramaswamy H. Sarma.

Genome organization regulates nuclear pore complex formation and promotes differentiation during Drosophila oogenesis

Genome organization can regulate gene expression and promote cell fate transitions. The differentiation of germline stem cells (GSCs) to oocytes in Drosophila involves changes in genome organization mediated by heterochromatin and the nuclear pore complex (NPC). Heterochromatin represses germ-cell genes during differentiation and NPCs anchor these silenced genes to the nuclear periphery, maintaining silencing to allow for oocyte development. Surprisingly, we find that genome organization also contributes to NPC formation, mediated by the transcription factor Stonewall (Stwl). As GSCs differentiate, Stwl accumulates at boundaries between silenced and active gene compartments. Stwl at these boundaries plays a pivotal role in transitioning germ-cell genes into a silenced state and activating a group of oocyte genes and Nucleoporins (Nups). The upregulation of these Nups during differentiation is crucial for NPC formation and further genome organization. Thus, crosstalk between genome architecture and NPCs is essential for successful cell fate transitions.

Mining the Utricularia gibba genome for insulator-like elements for genetic engineering

Introduction

Gene expression is often controlled via cis-regulatory elements (CREs) that modulate the production of transcripts. For multi-gene genetic engineering and synthetic biology, precise control of transcription is crucial, both to insulate the transgenes from unwanted native regulation and to prevent readthrough or cross-regulation of transgenes within a multi-gene cassette. To prevent this activity, insulator-like elements, more properly referred to as transcriptional blockers, could be inserted to separate the transgenes so that they are independently regulated. However, only a few validated insulator-like elements are available for plants, and they tend to be larger than ideal.

Methods

To identify additional potential insulator-like sequences, we conducted a genome-wide analysis of Utricularia gibba (humped bladderwort), one of the smallest known plant genomes, with genes that are naturally close together. The 10 best insulator-like candidates were evaluated in vivo for insulator-like activity.

Results

We identified a total of 4,656 intergenic regions with expression profiles suggesting insulator-like activity. Comparisons of these regions across 45 other plant species (representing Monocots, Asterids, and Rosids) show low levels of syntenic conservation of these regions. Genome-wide analysis of unmethylated regions (UMRs) indicates ~87% of the targeted regions are unmethylated; however, interpretation of this is complicated because U. gibba has remarkably low levels of methylation across the genome, so that large UMRs frequently extend over multiple genes and intergenic spaces. We also could not identify any conserved motifs among our selected intergenic regions or shared with existing insulator-like elements for plants. Despite this lack of conservation, however, testing of 10 selected intergenic regions for insulator-like activity found two elements on par with a previously published element (EXOB) while being significantly smaller.

Discussion

Given the small number of insulator-like elements currently available for plants, our results make a significant addition to available tools. The high hit rate (2 out of 10) also implies that more useful sequences are likely present in our selected intergenic regions; additional validation work will be required to identify which will be most useful for plant genetic engineering.

3D chromatin interactions involving Drosophila insulators are infrequent but preferential and arise before TADs and transcription

In mammals, insulators contribute to the regulation of loop extrusion to organize chromatin into topologically associating domains. In Drosophila the role of insulators in 3D genome organization is, however, under current debate. Here, we addressed this question by combining bioinformatics analysis and multiplexed chromatin imaging. We describe a class of Drosophila insulators enriched at regions forming preferential chromatin interactions genome-wide. Notably, most of these 3D interactions do not involve TAD borders. Multiplexed imaging shows that these interactions occur infrequently, and only rarely involve multiple genomic regions coalescing together in space in single cells. Finally, we show that non-border preferential 3D interactions enriched in this class of insulators are present before TADs and transcription during Drosophila development. Our results are inconsistent with insulators forming stable hubs in single cells, and instead suggest that they fine-tune existing 3D chromatin interactions, providing an additional regulatory layer for transcriptional regulation.

In silico discovery of repetitive elements as key sequence determinants of 3D genome folding

Natural and experimental genetic variants can modify DNA loops and insulating boundaries to tune transcription, but it is unknown how sequence perturbations affect chromatin organization genome wide. We developed a deep-learning strategy to quantify the effect of any insertion, deletion, or substitution on chromatin contacts and systematically scored millions of synthetic variants. While most genetic manipulations have little impact, regions with CTCF motifs and active transcription are highly sensitive, as expected. Our unbiased screen and subsequent targeted experiments also point to noncoding RNA genes and several families of repetitive elements as CTCF-motif-free DNA sequences with particularly large effects on nearby chromatin interactions, sometimes exceeding the effects of CTCF sites and explaining interactions that lack CTCF. We anticipate that our disruption tracks may be of broad interest and utility as a measure of 3D genome sensitivity, and our computational strategies may serve as a template for biological inquiry with deep learning.

The Conserved Chromatin Remodeler SMARCAD1 Interacts with TFIIIC and Architectural Proteins in Human and Mouse

In vertebrates, SMARCAD1 participates in transcriptional regulation, heterochromatin maintenance, DNA repair, and replication. The molecular basis underlying its involvement in these processes is not well understood. We identified the RNA polymerase III general transcription factor TFIIIC as an interaction partner of native SMARCAD1 in mouse and human models using endogenous co-immunoprecipitations. TFIIIC has dual functionality, acting as a general transcription factor and as a genome organizer separating chromatin domains. We found that its partnership with SMARCAD1 is conserved across different mammalian cell types, from somatic to pluripotent cells. Using purified proteins, we confirmed that their interaction is direct. A gene expression analysis suggested that SMARCAD1 is dispensable for TFIIIC function as an RNA polymerase III transcription factor in mouse ESCs. The distribution of TFIIIC and SMARCAD1 in the ESC genome is distinct, and unlike in yeast, SMARCAD1 is not enriched at active tRNA genes. Further analysis of SMARCAD1-binding partners in pluripotent and differentiated mammalian cells reveals that SMARCAD1 associates with several factors that have key regulatory roles in chromatin organization, such as cohesin, laminB, and DDX5. Together, our work suggests for the first time that the SMARCAD1 enzyme participates in genome organization in mammalian nuclei through interactions with architectural proteins.

The crucial role of CTCF in mitotic progression during early development of sea urchin

CCCTC-binding factor (CTCF), an insulator protein with 11 zinc fingers, is enriched at the boundaries of topologically associated domains (TADs) in eukaryotic genomes. In this study, we isolated and analyzed the cDNAs encoding HpCTCF, the CTCF homolog in the sea urchin Hemicentrotus pulcherrimus, to investigate its expression patterns and functions during the early development of sea urchin. HpCTCF contains nine zinc fingers corresponding to fingers 2-10 of the vertebrate CTCF. Expression pattern analysis revealed that HpCTCF mRNA was detected at all developmental stages and in the entire embryo. Upon expressing the HpCTCF-GFP fusion protein in early embryos, we observed its uniform distribution within interphase nuclei. However, during mitosis, it disappeared from the chromosomes and subsequently reassembled on the chromosome during telophase. Moreover, the morpholino-mediated knockdown of HpCTCF resulted in mitotic arrest during the morula to blastula stage. Most of the arrested chromosomes were not phospholylated at serine 10 of histone H3, indicating that mitosis was arrested at the telophase by HpCTCF depletion. Furthermore, impaired sister chromatid segregation was observed using time-lapse imaging of HpCTCF-knockdown embryos. Thus, HpCTCF is essential for mitotic progression during the early development of sea urchins, especially during the telophase-to-interphase transition. However, the normal development of pluteus larvae in CRISPR-mediated HpCTCF-knockout embryos suggests that disruption of zygotic HpCTCF expression has little effect on embryonic and larval development.

Multiomic analysis of cohesin reveals that ZBTB transcription factors contribute to chromatin interactions

One bottleneck in understanding the principles of 3D chromatin structures is caused by the paucity of known regulators. Cohesin is essential for 3D chromatin organization, and its interacting partners are candidate regulators. Here, we performed proteomic profiling of the cohesin in chromatin and identified transcription factors, RNA-binding proteins and chromatin regulators associated with cohesin. Acute protein degradation followed by time-series genomic binding quantitation and BAT Hi-C analysis were conducted, and the results showed that the transcription factor ZBTB21 contributes to cohesin chromatin binding, 3D chromatin interactions and transcriptional repression. Strikingly, multiomic analyses revealed that the other four ZBTB factors interacted with cohesin, and double degradation of ZBTB21 and ZBTB7B led to a further decrease in cohesin chromatin occupancy. We propose that multiple ZBTB transcription factors orchestrate the chromatin binding of cohesin to regulate chromatin interactions, and we provide a catalog of many additional proteins associated with cohesin that warrant further investigation.

Structural insights into human TFIIIC promoter recognition

Transcription factor (TF) IIIC recruits RNA polymerase (Pol) III to most of its target genes. Recognition of intragenic A- and B-box motifs in transfer RNA (tRNA) genes by TFIIIC modules τA and τB is the first critical step for tRNA synthesis but is mechanistically poorly understood. Here, we report cryo-electron microscopy structures of the six-subunit human TFIIIC complex unbound and bound to a tRNA gene. The τB module recognizes the B-box via DNA shape and sequence readout through the assembly of multiple winged-helix domains. TFIIIC220 forms an integral part of both τA and τB connecting the two subcomplexes via a ~550-amino acid residue flexible linker. Our data provide a structural mechanism by which high-affinity B-box recognition anchors TFIIIC to promoter DNA and permits scanning for low-affinity A-boxes and TFIIIB for Pol III activation.

LPAD: using network construction and label propagation to detect topologically associating domains from Hi-C data

With the development of chromosome conformation capture technique, the study of spatial conformation of a genome based on Hi-C technique has made a quantum leap. Previous studies reveal that genomes are folded into hierarchy of three-dimensional (3D) structures associated with topologically associating domains (TADs), and detecting TAD boundaries is of great significance in the chromosome-level analysis of 3D genome architecture. In this paper, we propose a novel TAD identification method, LPAD, which first extracts node correlations from global interactions of chromosomes based on the random walk with restart and then builds an undirected graph from Hi-C contact matrix. Next, LPAD designs a label propagation-based approach to discover communities and generates TADs. Experimental results verify the effectiveness and quality of TAD detections compared with existing methods. Furthermore, experimental evaluation of chromatin immunoprecipitation sequencing data shows that LPAD performs high enrichment of histone modifications remarkably nearby the TAD boundaries, and these results demonstrate LPAD's advantages on TAD identification accuracy.

Topologically associating domain boundaries are required for normal genome function

Topologically associating domain (TAD) boundaries partition the genome into distinct regulatory territories. Anecdotal evidence suggests that their disruption may interfere with normal gene expression and cause disease phenotypes^1-3, but the overall extent to which this occurs remains unknown. Here we demonstrate that targeted deletions of TAD boundaries cause a range of disruptions to normal in vivo genome function and organismal development. We used CRISPR genome editing in mice to individually delete eight TAD boundaries (11-80 kb in size) from the genome. All deletions examined resulted in detectable molecular or organismal phenotypes, which included altered chromatin interactions or gene expression, reduced viability, and anatomical phenotypes. We observed changes in local 3D chromatin architecture in 7 of 8 (88%) cases, including the merging of TADs and altered contact frequencies within TADs adjacent to the deleted boundary. For 5 of 8 (63%) loci examined, boundary deletions were associated with increased embryonic lethality or other developmental phenotypes. For example, a TAD boundary deletion near Smad3/Smad6 caused complete embryonic lethality, while a deletion near Tbx5/Lhx5 resulted in a severe lung malformation. Our findings demonstrate the importance of TAD boundary sequences for in vivo genome function and reinforce the critical need to carefully consider the potential pathogenicity of noncoding deletions affecting TAD boundaries in clinical genetics screening.

Somatic chromosome pairing has a determinant impact on 3D chromatin organization

In the nucleus, chromatin is intricately structured into multiple layers of 3D organization important for genome activity. How distinct layers influence each other is not well understood. In particular, the contribution of chromosome pairing to 3D chromatin organization has been largely neglected. Here, we address this question in Drosophila, an organism that shows robust chromosome pairing in interphasic somatic cells. The extent of chromosome pairing depends on the balance between pairing and anti-pairing factors, with the anti-pairing activity of the CAP-H2 condensin II subunit being the best documented. Here, we identify the zinc-finger protein Z4 as a strong anti-pairer that interacts with and mediates the chromatin binding of CAP-H2. We also report that hyperosmotic cellular stress induces fast and reversible chromosome unpairing that depends on Z4/CAP-H2. And, most important, by combining Z4 depletion and osmostress, we show that chromosome pairing reinforces intrachromosomal 3D interactions. On the one hand, pairing facilitates RNAPII occupancy that correlates with enhanced intragenic gene-loop interactions. In addition, acting at a distance, pairing reinforces chromatin-loop interactions mediated by Polycomb (Pc). In contrast, chromosome pairing does not affect which genomic intervals segregate to active (A) and inactive (B) compartments, with only minimal effects on the strength of A-A compartmental interactions. Altogether, our results unveil the intimate interplay between inter-chromosomal and intra-chromosomal 3D interactions, unraveling the interwoven relationship between different layers of chromatin organization and the essential contribution of chromosome pairing.

3D genomics and its applications in precision medicine

Three-dimensional (3D) genomics is an emerging discipline that studies the three-dimensional structure of chromatin and the three-dimensional and functions of genomes. It mainly studies the three-dimensional conformation and functional regulation of intranuclear genomes, such as DNA replication, DNA recombination, genome folding, gene expression regulation, transcription factor regulation mechanism, and the maintenance of three-dimensional conformation of genomes. Self-chromosomal conformation capture (3C) technology has been developed, and 3D genomics and related fields have developed rapidly. In addition, chromatin interaction analysis techniques developed by 3C technologies, such as paired-end tag sequencing (ChIA-PET) and whole-genome chromosome conformation capture (Hi-C), enable scientists to further study the relationship between chromatin conformation and gene regulation in different species. Thus, the spatial conformation of plant, animal, and microbial genomes, transcriptional regulation mechanisms, interaction patterns of chromosomes, and the formation mechanism of spatiotemporal specificity of genomes are revealed. With the help of new experimental technologies, the identification of key genes and signal pathways related to life activities and diseases is sustaining the rapid development of life science, agriculture, and medicine. In this paper, the concept and development of 3D genomics and its application in agricultural science, life science, and medicine are introduced, which provides a theoretical basis for the study of biological life processes.

CTCF, BEAF-32, and CP190 are not required for the establishment of TADs in early Drosophila embryos but have locus-specific roles

The boundaries of topologically associating domains (TADs) are delimited by insulators and/or active promoters; however, how they are initially established during embryogenesis remains unclear. Here, we examined this during the first hours of Drosophila embryogenesis. DNA-FISH confirms that intra-TAD pairwise proximity is established during zygotic genome activation (ZGA) but with extensive cell-to-cell heterogeneity. Most newly formed boundaries are occupied by combinations of CTCF, BEAF-32, and/or CP190. Depleting each insulator individually from chromatin revealed that TADs can still establish, although with lower insulation, with a subset of boundaries (~10%) being more dependent on specific insulators. Some weakened boundaries have aberrant gene expression due to unconstrained enhancer activity. However, the majority of misexpressed genes have no obvious direct relationship to changes in domain-boundary insulation. Deletion of an active promoter (thereby blocking transcription) at one boundary had a greater impact than deleting the insulator-bound region itself. This suggests that cross-talk between insulators and active promoters and/or transcription might reinforce domain boundary insulation during embryogenesis.

Nuclear architecture and the structural basis of mitotic memory

The nucleus is a complex organelle that hosts the genome and is essential for vital processes like DNA replication, DNA repair, transcription, and splicing. The genome is non-randomly organized in the three-dimensional space of the nucleus. This functional sub-compartmentalization was thought to be organized on the framework of nuclear matrix (NuMat), a non-chromatin scaffold that functions as a substratum for various molecular processes of the nucleus. More recently, nuclear bodies or membrane-less subcompartments of the nucleus are thought to arise due to phase separation of chromatin, RNA, and proteins. The nuclear architecture is an amalgamation of the relative organization of chromatin, epigenetic landscape, the nuclear bodies, and the nucleoskeleton in the three-dimensional space of the nucleus. During mitosis, the nucleus undergoes drastic changes in morphology to the degree that it ceases to exist as such; various nuclear components, including the envelope that defines the nucleus, disintegrate, and the chromatin acquires mitosis-specific epigenetic marks and condenses to form chromosome. Upon mitotic exit, chromosomes are decondensed, re-establish hierarchical genome organization, and regain epigenetic and transcriptional status similar to that of the mother cell. How this mitotic memory is inherited during cell division remains a puzzle. NuMat components that are a part of the mitotic chromosome in the form of mitotic chromosome scaffold (MiCS) could potentially be the seeds that guide the relative re-establishment of the epigenome, chromosome territories, and the nuclear bodies. Here, we synthesize the advances towards understanding cellular memory of nuclear architecture across mitosis and propose a hypothesis that a subset of NuMat proteome essential for nucleation of various nuclear bodies are retained in MiCS to serve as seeds of mitotic memory, thus ensuring the daughter cells re-establish the complex status of nuclear architecture similar to that of the mother cells, thereby maintaining the pre-mitotic transcriptional status.

Drosophila SUMM4 complex couples insulator function and DNA replication control

Asynchronous replication of chromosome domains during S phase is essential for eukaryotic genome function, but the mechanisms establishing which domains replicate early versus late in different cell types remain incompletely understood. Intercalary heterochromatin domains replicate very late in both diploid chromosomes of dividing cells and in endoreplicating polytene chromosomes where they are also underreplicated. Drosophila SNF2-related factor SUUR imparts locus-specific underreplication of polytene chromosomes. SUUR negatively regulates DNA replication fork progression; however, its mechanism of action remains obscure. Here, we developed a novel method termed MS-Enabled Rapid protein Complex Identification (MERCI) to isolate a stable stoichiometric native complex SUMM4 that comprises SUUR and a chromatin boundary protein Mod(Mdg4)-67.2. Mod(Mdg4) stimulates SUUR ATPase activity and is required for a normal spatiotemporal distribution of SUUR in vivo. SUUR and Mod(Mdg4)-67.2 together mediate the activities of gypsy insulator that prevent certain enhancer-promoter interactions and establish euchromatin-heterochromatin barriers in the genome. Furthermore, SuUR or mod(mdg4) mutations reverse underreplication of intercalary heterochromatin. Thus, SUMM4 can impart late replication of intercalary heterochromatin by attenuating the progression of replication forks through euchromatin/heterochromatin boundaries. Our findings implicate a SNF2 family ATP-dependent motor protein SUUR in the insulator function, reveal that DNA replication can be delayed by a chromatin barrier, and uncover a critical role for architectural proteins in replication control. They suggest a mechanism for the establishment of late replication that does not depend on an asynchronous firing of late replication origins.

Drosophila insulator proteins exhibit in vivo liquid-liquid phase separation properties

Drosophila insulator proteins and the cohesin subunit Rad21 coalesce in vivo to form liquid-droplet condensates, suggesting that liquid–liquid phase separation mediates their function in 3D genome organization.

Mounting evidence implicates liquid–liquid phase separation (LLPS), the condensation of biomolecules into liquid-like droplets in the formation and dissolution of membraneless intracellular organelles (MLOs). Cells use MLOs or condensates for various biological processes, including emergency signaling and spatiotemporal control over steady-state biochemical reactions and heterochromatin formation. Insulator proteins are architectural elements involved in establishing independent domains of transcriptional activity within eukaryotic genomes. In Drosophila, insulator proteins form nuclear foci known as insulator bodies in response to osmotic stress. However, the mechanism through which insulator proteins assemble into bodies is yet to be investigated. Here, we identify signatures of LLPS by insulator bodies, including high disorder tendency in insulator proteins, scaffold–client–dependent assembly, extensive fusion behavior, sphericity, and sensitivity to 1,6-hexanediol. We also show that the cohesin subunit Rad21 is a component of insulator bodies, adding to the known insulator protein constituents and γH2Av. Our data suggest a concerted role of cohesin and insulator proteins in insulator body formation and under physiological conditions. We propose a mechanism whereby these architectural proteins modulate 3D genome organization through LLPS.

Phosphorylated histone variant γH2Av is associated with chromatin insulators in Drosophila

Chromatin insulators are responsible for orchestrating long-range interactions between enhancers and promoters throughout the genome and align with the boundaries of Topologically Associating Domains (TADs). Here, we demonstrate an association between gypsy insulator proteins and the phosphorylated histone variant H2Av (γH2Av), normally a marker of DNA double strand breaks. Gypsy insulator components colocalize with γH2Av throughout the genome, in polytene chromosomes and in diploid cells in which Chromatin IP data shows it is enriched at TAD boundaries. Mutation of insulator components su(Hw) and Cp190 results in a significant reduction in γH2Av levels in chromatin and phosphatase inhibition strengthens the association between insulator components and γH2Av and rescues γH2Av localization in insulator mutants. We also show that γH2Av, but not H2Av, is a component of insulator bodies, which are protein condensates that form during osmotic stress. Phosphatase activity is required for insulator body dissolution after stress recovery. Together, our results implicate the H2A variant with a novel mechanism of insulator function and boundary formation.

CohesinDB: a comprehensive database for decoding cohesin-related epigenomes, 3D genomes and transcriptomes in human cells

Cohesin is a multifunctional protein responsible for transcriptional regulation and chromatin organization. Cohesin binds to chromatin at tens of thousands of distinct sites in a conserved or tissue-specific manner, whereas the function of cohesin varies greatly depending on the epigenetic properties of specific chromatin loci. Cohesin also extensively mediates cis-regulatory modules (CRMs) and chromatin loops. Even though next-generation sequencing technologies have provided a wealth of information on different aspects of cohesin, the integration and exploration of the resultant massive cohesin datasets are not straightforward. Here, we present CohesinDB (https://cohesindb.iqb.u-tokyo.ac.jp), a comprehensive multiomics cohesin database in human cells. CohesinDB includes 2043 epigenomics, transcriptomics and 3D genomics datasets from 530 studies involving 176 cell types. By integrating these large-scale data, CohesinDB summarizes three types of 'cohesin objects': 751 590 cohesin binding sites, 957 868 cohesin-related chromatin loops and 2 229 500 cohesin-related CRMs. Each cohesin object is annotated with locus, cell type, classification, function, 3D genomics and cis-regulatory information. CohesinDB features a user-friendly interface for browsing, searching, analyzing, visualizing and downloading the desired information. CohesinDB contributes a valuable resource for all researchers studying cohesin, epigenomics, transcriptional regulation and chromatin organization.

NURF301 contributes to gypsy chromatin insulator-mediated nuclear organization

Chromatin insulators are DNA-protein complexes that can prevent the spread of repressive chromatin and block communication between enhancers and promoters to regulate gene expression. In Drosophila, the gypsy chromatin insulator complex consists of three core proteins: CP190, Su(Hw), and Mod(mdg4)67.2. These factors concentrate at nuclear foci termed insulator bodies, and changes in insulator body localization have been observed in mutants defective for insulator function. Here, we identified NURF301/E(bx), a nucleosome remodeling factor, as a novel regulator of gypsy insulator body localization through a high-throughput RNAi imaging screen. NURF301 promotes gypsy-dependent insulator barrier activity and physically interacts with gypsy insulator proteins. Using ChIP-seq, we found that NURF301 co-localizes with insulator proteins genome-wide, and NURF301 promotes chromatin association of Su(Hw) and CP190 at gypsy insulator binding sites. These effects correlate with NURF301-dependent nucleosome repositioning. At the same time, CP190 and Su(Hw) both facilitate recruitment of NURF301 to chromatin. Finally, Oligopaint FISH combined with immunofluorescence revealed that NURF301 promotes 3D contact between insulator bodies and gypsy insulator DNA binding sites, and NURF301 is required for proper nuclear positioning of gypsy binding sites. Our data provide new insights into how a nucleosome remodeling factor and insulator proteins cooperatively contribute to nuclear organization.

Breaking boundaries: Pan BETi disrupt 3D chromatin structure, BD2-selective BETi are strictly epigenetic transcriptional regulators

BACKGROUND

Bromodomain and extraterminal proteins (BETs) are more than just epigenetic regulators of transcription. Here we highlight a new role for the BET protein BRD4 in the maintenance of higher order chromatin structure at Topologically Associating Domain Boundaries (TADBs). BD2-selective and pan (non-selective) BET inhibitors (BETi) differentially support chromatin structure, selectively affecting transcription and cell viability.

METHODS

Using RNA-seq and BRD4 ChIP-seq, the differential effect of BETi treatment on the transcriptome and BRD4 chromatin occupancy of human aortic endothelial cells from diabetic patients (dHAECs) stimulated with TNFα was evaluated. Chromatin decondensation and DNA fragmentation was assessed by immunofluorescence imaging and quantification. Key dHAEC findings were verified in proliferating monocyte-like THP-1 cells using real time-PCR, BRD4 co-immunoprecipitation studies, western blots, proliferation and apoptosis assays.

FINDINGS

We discovered that 1) BRD4 co-localizes with Ying-Yang 1 (YY1) at TADBs, critical chromatin structure complexes proximal to many DNA repair genes. 2) BD2-selective BETi enrich BRD4/YY1 associations, while pan-BETi do not. 3) Failure to support chromatin structures through BRD4/YY1 enrichment inhibits DNA repair gene transcription, which induces DNA damage responses, and causes widespread chromatin decondensation, DNA fragmentation, and apoptosis. 4) BD2-selective BETi maintain high order chromatin structure and cell viability, while reducing deleterious pro-inflammatory transcription.

INTERPRETATION

BRD4 plays a previously unrecognized role at TADBs. BETi differentially impact TADB stability. Our results provide translational insight for the development of BETi as therapeutics for a range of diseases including CVD, chronic kidney disease, cancer, and COVID-19.

TFIIIC-based chromatin insulators through eukaryotic evolution

Eukaryotic chromosomes are divided into domains with distinct structural and functional properties, such as differing levels of chromatin compaction and gene transcription. Domains of relatively compact chromatin and minimal transcription are termed heterochromatic, whereas euchromatin is more open and actively transcribed. Insulators separate these domains and maintain their distinct features. Disruption of insulators can cause diseases such as cancer. Many insulators contain tRNA genes (tDNAs), examples of which have been shown to block the spread of activating or silencing activities. This characteristic of specific tDNAs is conserved through evolution, such that human tDNAs can serve as barriers to the spread of silencing in fission yeast. Here we demonstrate that tDNAs from the methylotrophic fungus Pichia pastoris can function effectively as insulators in distantly-related budding yeast. Key to the function of tDNAs as insulators is TFIIIC, a transcription factor that is also required for their expression. TFIIIC binds additional loci besides tDNAs, some of which have insulator activity. Although the mechanistic basis of TFIIIC-based insulation has been studied extensively in yeast, it is largely uncharacterized in metazoa. Utilising publicly-available genome-wide ChIP-seq data, we consider the extent to which mechanisms conserved from yeast to man may suffice to allow efficient insulation by TFIIIC in the more challenging chromatin environments of metazoa and suggest features that may have been acquired during evolution to cope with new challenges. We demonstrate the widespread presence at human tDNAs of USF1, a transcription factor with well-established barrier activity in vertebrates. We predict that tDNA-based insulators in higher organisms have evolved through incorporation of modules, such as binding sites for factors like USF1 and CTCF that are absent from yeasts, thereby strengthening function and providing opportunities for regulation between cell types.

Large-scale chromatin reorganization reactivates placenta-specific genes that drive cellular aging

Nuclear deformation, a hallmark frequently observed in senescent cells, is presumed to be associated with the erosion of chromatin organization at the nuclear periphery. However, how such gradual changes in higher-order genome organization impinge on local epigenetic modifications to drive cellular mechanisms of aging has remained enigmatic. Here, through large-scale epigenomic analyses of isogenic young, senescent, and progeroid human mesenchymal progenitor cells (hMPCs), we delineate a hierarchy of integrated structural state changes that manifest as heterochromatin loss in repressive compartments, euchromatin weakening in active compartments, switching in interfacing topological compartments, and increasing epigenetic entropy. We found that the epigenetic de-repression unlocks the expression of pregnancy-specific beta-1 glycoprotein (PSG) genes that exacerbate hMPC aging and serve as potential aging biomarkers. Our analyses provide a rich resource for uncovering the principles of epigenomic landscape organization and its changes in cellular aging and for identifying aging drivers and intervention targets with a genome-topology-based mechanism.

Essential role of Cp190 in physical and regulatory boundary formation

Boundaries in animal genomes delimit contact domains with enhanced internal contact frequencies and have debated functions in limiting regulatory cross-talk between domains and guiding enhancers to target promoters. Most mammalian boundaries form by stalling of chromosomal loop-extruding cohesin by CTCF, but most Drosophila boundaries form CTCF independently. However, how CTCF-independent boundaries form and function remains largely unexplored. Here, we assess genome folding and developmental gene expression in fly embryos lacking the ubiquitous boundary-associated factor Cp190. We find that sequence-specific DNA binding proteins such as CTCF and Su(Hw) directly interact with and recruit Cp190 to form most promoter-distal boundaries. Cp190 is essential for early development and prevents regulatory cross-talk between specific gene loci that pattern the embryo. Cp190 was, in contrast, dispensable for long-range enhancer-promoter communication at tested loci. Cp190 is thus currently the major player in fly boundary formation and function, revealing that diverse mechanisms evolved to partition genomes into independent regulatory domains.

Transcriptional coupling of distant regulatory genes in living embryos

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

The importance of considering regulatory domains in genome-wide analyses - the nearest gene is often wrong!

The expression of a large number of genes is regulated by regulatory elements that are located far away from their promoters. Identifying which gene is the target of a specific regulatory element or is affected by a non-coding mutation is often accomplished by assigning these regions to the nearest gene in the genome. However, this heuristic ignores key features of genome organisation and gene regulation; in that the genome is partitioned into regulatory domains, which at some loci directly coincide with the span of topologically associated domains (TADs), and that genes are regulated by enhancers located throughout these regions, even across intervening genes. In this review, we examine the results from genome-wide studies using chromosome conformation capture technologies and from those dissecting individual gene regulatory domains, to highlight that the phenomenon of enhancer skipping is pervasive and affects multiple types of genes. We discuss how simply assigning a genomic region of interest to its nearest gene is problematic and often leads to incorrect predictions and highlight that where possible information on both the conservation and topological organisation of the genome should be used to generate better hypotheses.

The article has an associated Future Leader to Watch interview.

Summary: Identifying which gene is the target of an enhancer is often accomplished by assigning it to the nearest gene, here we discuss how this heuristic can lead to incorrect predictions.

The role of insulators and transcription in 3D chromatin organization of flies

The DNA in many organisms, including humans, is shown to be organized in topologically associating domains (TADs). In Drosophila, several architectural proteins are enriched at TAD borders, but it is still unclear whether these proteins play a functional role in the formation and maintenance of TADs. Here, we show that depletion of BEAF-32, Cp190, Chro, and Dref leads to changes in TAD organization and chromatin loops. Their depletion predominantly affects TAD borders located in regions moderately enriched in repressive modifications and depleted in active ones, whereas TAD borders located in euchromatin are resilient to these knockdowns. Furthermore, transcriptomic data has revealed hundreds of genes displaying differential expression in these knockdowns and showed that the majority of differentially expressed genes are located within reorganized TADs. Our work identifies a novel and functional role for architectural proteins at TAD borders in Drosophila and a link between TAD reorganization and subsequent changes in gene expression.

Functional Similarities of Protein-Coding Genes in Topologically Associating Domains and Spatially-Proximate Genomic Regions

Topologically associating domains (TADs) are the structural and functional units of the genome. However, the functions of protein-coding genes existing in the same or different TADs have not been fully investigated. We compared the functional similarities of protein-coding genes existing in the same TAD and between different TADs, and also in the same gap region (the region between two consecutive TADs) and between different gap regions. We found that the protein-coding genes from the same TAD or gap region are more likely to share similar protein functions, and this trend is more obvious with TADs than the gap regions. We further created two types of gene–gene spatial interaction networks: the first type is based on Hi-C contacts, whereas the second type is based on both Hi-C contacts and the relationship of being in the same TAD. A graph auto-encoder was applied to learn the network topology, reconstruct the two types of networks, and predict the functions of the central genes/nodes based on the functions of the neighboring genes/nodes. It was found that better performance was achieved with the second type of network. Furthermore, we detected long-range spatially-interactive regions based on Hi-C contacts and calculated the functional similarities of the gene pairs from these regions.

CRISPR and biochemical screens identify MAZ as a cofactor in CTCF-mediated insulation at Hox clusters

CCCTC-binding factor (CTCF) is critical to three-dimensional genome organization. Upon differentiation, CTCF insulates active and repressed genes within Hox gene clusters. We conducted a genome-wide CRISPR knockout (KO) screen to identify genes required for CTCF-boundary activity at the HoxA cluster, complemented by biochemical approaches. Among the candidates, we identified Myc-associated zinc-finger protein (MAZ) as a cofactor in CTCF insulation. MAZ colocalizes with CTCF at chromatin borders and, similar to CTCF, interacts with the cohesin subunit RAD21. MAZ KO disrupts gene expression and local contacts within topologically associating domains. Similar to CTCF motif deletions, MAZ motif deletions lead to derepression of posterior Hox genes immediately after CTCF boundaries upon differentiation, giving rise to homeotic transformations in mouse. Thus, MAZ is a factor contributing to appropriate insulation, gene expression and genomic architecture during development.

Rice RS2-9, which is bound by transcription factor OSH1, blocks enhancer-promoter interactions in plants

Insulators characterized in Drosophila and mammals have been shown to play a key role in the restriction of promiscuous enhancer-promoter interactions, as well as reshaping the topological landscape of chromosomes. Yet the role of insulators in plants remains poorly understood, in large part because of a lack of well-characterized insulators and binding factor(s). In this study, we isolated a 1.2-kb RS2-9 insulator from the Oryza sativa (rice) genome that can, when interposed between an enhancer and promoter, efficiently block the activation function of both constitutive and floral organ-specific enhancers in transgenic Arabidopsis and Nicotiana tabacum (tobacco). In the rice genome, the genes flanking RS2-9 exhibit an absence of mutual transcriptional interactions, as well as a lack of histone modification spread. We further determined that O. sativa Homeobox 1 (OSH1) bound two regions of RS2-9, as well as over 50 000 additional sites in the rice genome, the majority of which resided in intergenic regions. Mutation of one of the two OSH1-binding sites in RS2-9 impaired insulation activity by up to 60%, whereas the mutation of both binding sites virtually abolished insulator function. We also demonstrated that OSH1 binding sites were associated with 72% of the boundaries of topologically associated domains (TADs) identified in the rice genome, which is comparable to the 77% of TAD boundaries bound by the insulator CCCTC-binding factor (CTCF) in mammals. Taken together, our findings indicate that OSH1-RS2-9 acts as a true insulator in plants, and highlight a potential role for OSH1 in gene insulation and topological organization in plant genomes.

OUP accepted manuscript

Eukaryotic chromosomes are spatially segregated into topologically associating domains (TADs). Some TADs are attached to the nuclear lamina (NL) through lamina-associated domains (LADs). Here, we identified LADs and TADs at two stages of Drosophila spermatogenesis – in bam^Δ86 mutant testes which is the commonly used model of spermatogonia (SpG) and in larval testes mainly filled with spermatocytes (SpCs). We found that initiation of SpC-specific transcription correlates with promoters’ detachment from the NL and with local spatial insulation of adjacent regions. However, this insulation does not result in the partitioning of inactive TADs into sub-TADs. We also revealed an increased contact frequency between SpC-specific genes in SpCs implying their de novo gathering into transcription factories. In addition, we uncovered the specific X chromosome organization in the male germline. In SpG and SpCs, a single X chromosome is stronger associated with the NL than autosomes. Nevertheless, active chromatin regions in the X chromosome interact with each other more frequently than in autosomes. Moreover, despite the absence of dosage compensation complex in the male germline, randomly inserted SpG-specific reporter is expressed higher in the X chromosome than in autosomes, thus evidencing that non-canonical dosage compensation operates in SpG.

Cell-type-specific chromatin occupancy by the pioneer factor Zelda drives key developmental transitions in Drosophila

During Drosophila embryogenesis, the essential pioneer factor Zelda defines hundreds of cis-regulatory regions and in doing so reprograms the zygotic transcriptome. While Zelda is essential later in development, it is unclear how the ability of Zelda to define cis-regulatory regions is shaped by cell-type-specific chromatin architecture. Asymmetric division of neural stem cells (neuroblasts) in the fly brain provide an excellent paradigm for investigating the cell-type-specific functions of this pioneer factor. We show that Zelda synergistically functions with Notch to maintain neuroblasts in an undifferentiated state. Zelda misexpression reprograms progenitor cells to neuroblasts, but this capacity is limited by transcriptional repressors critical for progenitor commitment. Zelda genomic occupancy in neuroblasts is reorganized as compared to the embryo, and this reorganization is correlated with differences in chromatin accessibility and cofactor availability. We propose that Zelda regulates essential transitions in the neuroblasts and embryo through a shared gene-regulatory network driven by cell-type-specific enhancers.

HiC1Dmetrics: framework to extract various one-dimensional features from chromosome structure data

Eukaryotic genomes are organized in a three-dimensional spatial structure. In this regard, the development of chromosome conformation capture methods has enabled studies of chromosome organization on a genomic scale. Hi-C, the high-throughput chromosome conformation capture method, can reveal a population-averaged, hierarchical chromatin structure. The typical Hi-C analysis uses a two-dimensional (2D) contact matrix that indicates contact frequencies between all possible genomic position pairs. Oftentimes, however, such a 2D matrix is not amenable to handling quantitative comparisons, visualizations and integrations across multiple datasets. Although several one-dimensional (1D) metrics have been proposed to depict structural information in Hi-C data, their effectiveness is still underappreciated. Here, we first review the currently available 1D metrics for individual Hi-C samples or two-sample comparisons and then discuss their validity and suitable analysis scenarios. We also propose several new 1D metrics to identify additional unique features of chromosome structures. We highlight that the 1D metrics are reproducible and robust for comparing and visualizing multiple Hi-C samples. Moreover, we show that 1D metrics can be easily combined with epigenome tracks to annotate chromatin states in greater details. We develop a new framework, called HiC1Dmetrics, to summarize all 1D metrics discussed in this study. HiC1Dmetrics is open-source (github.com/wangjk321/HiC1Dmetrics) and can be accessed from both command-line and web-based interfaces. Our tool constitutes a useful resource for the community of chromosome-organization researchers.

Insulators in Plants: Progress and Open Questions

The genomes of higher eukaryotes are partitioned into topologically associated domains or TADs, and insulators (also known as boundary elements) are the key elements responsible for their formation and maintenance. Insulators were first identified and extensively studied in Drosophila as well as mammalian genomes, and have also been described in yeast and plants. In addition, many insulator proteins are known in Drosophila, and some have been investigated in mammals. However, much less is known about this important class of non-coding DNA elements in plant genomes. In this review, we take a detailed look at known plant insulators across different species and provide an overview of potential determinants of plant insulator functions, including cis-elements and boundary proteins. We also discuss methods previously used in attempts to identify plant insulators, provide a perspective on their importance for research and biotechnology, and discuss areas of potential future research.

Drosophila architectural protein CTCF is not essential for fly survival and is able to function independently of CP190

CTCF is the most likely ancestor of proteins that contain large clusters of C2H2 zinc finger domains (C2H2) and is conserved among most bilateral organisms. In mammals, CTCF functions as the main architectural protein involved in the organization of topology-associated domains (TADs). In vertebrates and Drosophila, CTCF is involved in the regulation of homeotic genes. Previously, it was found that null mutations in the dCTCF gene died as pharate adults, which failed to eclose from their pupal case, or shortly after hatching of adults. Here, we obtained several new null dCTCF mutations and found that the complete inactivation of dCTCF appears is limited mainly to phenotypic manifestations of the Abd-B gene and fertility of adult flies. Many modifiers that are not associated with an independent phenotypic manifestation can significantly enhance the expressivity of the null dCTCF mutations, indicating that other architectural proteins are able to functionally compensate for dCTCF inactivation in Drosophila. We also mapped the 715-735 aa region of dCTCF as being essential for the interaction with the BTB (Broad-Complex, Tramtrack, and Bric a brac) and microtubule-targeting (M) domains of the CP190 protein, which binds to many architectural proteins. However, the mutational analysis showed that the interaction with CP190 was not important for the functional activity of dCTCF in vivo.

3D genome organisation in Drosophila

Ever since Thomas Hunt Morgan's discovery of the chromosomal basis of inheritance by using Drosophila melanogaster as a model organism, the fruit fly has remained an essential model system in studies of genome biology, including chromatin organisation. Very much as in vertebrates, in Drosophila, the genome is organised in territories, compartments and topologically associating domains (TADs). However, these domains might be formed through a slightly different mechanism than in vertebrates due to the presence of a large and potentially redundant set of insulator proteins and the minor role of dCTCF in TAD boundary formation. Here, we review the different levels of chromatin organisation in Drosophila and discuss mechanisms and factors that might be involved in TAD formation. The dynamics of TADs and enhancer-promoter interactions in the context of transcription are covered in the light of currently conflicting results. Finally, we illustrate the value of polymer modelling approaches to infer the principles governing the three-dimensional organisation of the Drosophila genome.

M1BP cooperates with CP190 to activate transcription at TAD borders and promote chromatin insulator activity

Genome organization is driven by forces affecting transcriptional state, but the relationship between transcription and genome architecture remains unclear. Here, we identified the Drosophila transcription factor Motif 1 Binding Protein (M1BP) in physical association with the gypsy chromatin insulator core complex, including the universal insulator protein CP190. M1BP is required for enhancer-blocking and barrier activities of the gypsy insulator as well as its proper nuclear localization. Genome-wide, M1BP specifically colocalizes with CP190 at Motif 1-containing promoters, which are enriched at topologically associating domain (TAD) borders. M1BP facilitates CP190 chromatin binding at many shared sites and vice versa. Both factors promote Motif 1-dependent gene expression and transcription near TAD borders genome-wide. Finally, loss of M1BP reduces chromatin accessibility and increases both inter- and intra-TAD local genome compaction. Our results reveal physical and functional interaction between CP190 and M1BP to activate transcription at TAD borders and mediate chromatin insulator-dependent genome organization.

Transcriptional state plays a role in genome organization, however factors that link these processes are not well known. Here, the authors show Drosophila transcription factor Motif 1-binding protein (M1BP) interacts with the insulator protein CP190 to promote insulator function and activate Motif 1-dependent transcription at topologically associating domain (TAD) borders.