Sub-kb Hi-C in D. melanogaster reveals conserved characteristics of TADs between insect and mammalian cells

Heterochromatic 3D genome organization is directed by HP1a- and H3K9-dependent and independent mechanisms

Whether and how histone post-translational modifications and the proteins that bind them drive 3D genome organization remains unanswered. Here, we evaluate the contribution of H3K9-methylated constitutive heterochromatin to 3D genome organization in Drosophila tissues. We find that the predominant organizational feature of wild-type tissues is the segregation of euchromatic chromosome arms from heterochromatic pericentromeres. Reciprocal perturbation of HP1a⋅H3K9me binding, using a point mutation in the HP1a chromodomain or replacement of the replication-dependent histone H3 with H3K9R mutant histones, revealed that HP1a binding to methylated H3K9 in constitutive heterochromatin is required to limit contact frequency between pericentromeres and chromosome arms and regulate the distance between arm and pericentromeric regions. Surprisingly, the self-association of pericentromeric regions is largely preserved despite the loss of H3K9 methylation and HP1a occupancy. Thus, the HP1a⋅H3K9 interaction contributes to but does not solely drive the segregation of euchromatin and heterochromatin inside the nucleus.

Development of a New Model System to Study Long-Distance Interactions Supported by Architectural Proteins

Chromatin architecture is critical for the temporal and tissue-specific activation of genes that determine eukaryotic development. The functional interaction between enhancers and promoters is controlled by insulators and tethering elements that support specific long-distance interactions. However, the mechanisms of the formation and maintenance of long-range interactions between genome regulatory elements remain poorly understood, primarily due to the lack of convenient model systems. Drosophila became the first model organism in which architectural proteins that determine the activity of insulators were described. In Drosophila, one of the best-studied DNA-binding architectural proteins, Su(Hw), forms a complex with Mod(mdg4)-67.2 and CP190 proteins. Using a combination of CRISPR/Cas9 genome editing and attP-dependent integration technologies, we created a model system in which the promoters and enhancers of two reporter genes are separated by 28 kb. In this case, enhancers effectively stimulate reporter gene promoters in cis and trans only in the presence of artificial Su(Hw) binding sites (SBS), in both constructs. The expression of the mutant Su(Hw) protein, which cannot interact with CP190, and the mutation inactivating Mod(mdg4)-67.2, lead to the complete loss or significant weakening of enhancer-promoter interactions, respectively. The results indicate that the new model system effectively identifies the role of individual subunits of architectural protein complexes in forming and maintaining specific long-distance interactions in the D. melanogaster model.

The N-terminal dimerization domains of human and Drosophila CTCF have similar functionality

BACKGROUND

CTCF is highly likely to be the ancestor of proteins that contain large clusters of C2H2 zinc finger domains, and its conservation is observed across most bilaterian organisms. In mammals, CTCF is the primary architectural protein involved in organizing chromosome topology and mediating enhancer-promoter interactions over long distances. In Drosophila, CTCF (dCTCF) cooperates with other architectural proteins to establish long-range interactions and chromatin boundaries. CTCFs of various organisms contain an unstructured N-terminal dimerization domain (DD) and clusters comprising eleven zinc-finger domains of the C2H2 type. The Drosophila (dCTCF) and human (hCTCF) CTCFs share sequence homology in only five C2H2 domains that specifically bind to a conserved 15 bp motif.

RESULTS

Previously, we demonstrated that CTCFs from different organisms carry unstructured N-terminal dimerization domains (DDs) that lack sequence homology. Here we used the CTCFattP(mCh) platform to introduce desired changes in the Drosophila CTCF gene and generated a series of transgenic lines expressing dCTCF with different variants of the N-terminal domain. Our findings revealed that the functionality of dCTCF is significantly affected by the deletion of the N-terminal DD. Additionally, we observed a strong impact on the binding of the dCTCF mutant to chromatin upon deletion of the DD. However, chromatin binding was restored in transgenic flies expressing a chimeric CTCF protein with the DD of hCTCF. Although the chimeric protein exhibited lower expression levels than those of the dCTCF variants, it efficiently bound to chromatin similarly to the wild type (wt) protein.

CONCLUSIONS

Our findings suggest that one of the evolutionarily conserved functions of the unstructured N-terminal dimerization domain is to recruit dCTCF to its genomic sites in vivo.

Smoother: on-the-fly processing of interactome data using prefix sums

Nucleic acid interactome data, such as chromosome conformation capture data and RNA-DNA interactome data, are currently analyzed via pipelines that must be rerun for each new parameter set. A more dynamic approach is desirable since the optimal parameter set is commonly unknown ahead of time and rerunning pipelines is a time-consuming process. We have developed an approach fast enough to process interactome data on-the-fly using a sparse prefix sum index. With this index, we created Smoother, a flexible, multifeatured visualization and analysis tool that allows interactive filtering, e.g. by mapping quality, almost instant comparisons between different normalization approaches, e.g. iterative correction, and ploidy correction. Further, Smoother can overlay other sequencing data or genomic annotations, compare different samples, and perform virtual 4C analysis. Smoother permits a novel way to interact with and explore interactome data, fostering comprehensive, high-quality data analysis. Smoother is available at https://github.com/Siegel-Lab/BioSmoother under the MIT license.

Stage-specific dynamic reorganization of genome topology shapes transcriptional neighborhoods in developing human retinal organoids

We have generated a high-resolution Hi-C map of developing human retinal organoids to elucidate spatiotemporal dynamics of genomic architecture and its relationship with gene expression patterns. We demonstrate progressive stage-specific alterations in DNA topology and correlate these changes with transcription of cell-type-restricted gene markers during retinal differentiation. Temporal Hi-C reveals a shift toward A compartment for protein-coding genes and B compartment for non-coding RNAs, displaying high and low expression, respectively. Notably, retina-enriched genes are clustered near lost boundaries of topologically associated domains (TADs), and higher-order assemblages (i.e., TAD cliques) localize in active chromatin regions with binding sites for eye-field transcription factors. These genes gain chromatin contacts at their transcription start site as organoid differentiation proceeds. Our study provides a global view of chromatin architecture dynamics associated with diversification of cell types during retinal development and serves as a foundational resource for in-depth functional investigations of retinal developmental traits.

RNA polymerases reshape chromatin and coordinate transcription on individual fibers

During eukaryotic transcription, RNA polymerases must initiate and pause within a crowded, complex environment, surrounded by nucleosomes and other transcriptional activity. This environment creates a spatial arrangement along individual chromatin fibers ripe for both competition and coordination, yet these relationships remain largely unknown owing to the inherent limitations of traditional structural and sequencing methodologies. To address these limitations, we employed long-read chromatin fiber sequencing (Fiber-seq) to visualize RNA polymerases within their native chromatin context at single-molecule and near single-nucleotide resolution along up to 30 kb fibers. We demonstrate that Fiber-seq enables the identification of single-molecule RNA Polymerase (Pol) II and III transcription associated footprints, which, in aggregate, mirror bulk short-read sequencing-based measurements of transcription. We show that Pol II pausing destabilizes downstream nucleosomes, with frequently paused genes maintaining a short-term memory of these destabilized nucleosomes. Furthermore, we demonstrate pervasive direct coordination and anti-coordination between nearby Pol II genes, Pol III genes, transcribed enhancers, and insulator elements. This coordination is largely limited to spatially organized elements within 5 kb of each other, implicating short-range chromatin environments as a predominant determinant of coordinated polymerase initiation. Overall, transcription initiation reshapes surrounding nucleosome architecture and coordinates nearby transcriptional machinery along individual chromatin fibers.

ZFH-2 is required for Drosophila ovarian follicle development and is expressed at the band/interband boundaries of polytene chromosomes

The transcription factor ZFH-2 has well-documented roles in Drosophila neurogenesis and other developmental processes. Here we provide the first evidence that ZFH-2 has a role in oogenesis. We demonstrate that ZFH-2 is expressed in the wild-type ovary and that a loss of zfh-2 function produces a mutant ovary phenotype where egg chambers are reduced in number and fused. We also show that a loss of zfh-2 function can suppress a daughterless loss-of-function ovary phenotype suggesting a possible genetic relationship between these two genes in the ovary. We also show that ZFH-2 is located at the boundary between bands and interbands on polytene chromosomes and that at a subset of these sites ZFH-2 colocalizes with the insulator/promoter cofactor CP190.

The N-Terminal Part of Drosophila CP190 Is a Platform for Interaction with Multiple Architectural Proteins

CP190 is a co-factor in many Drosophila architectural proteins, being involved in the formation of active promoters and insulators. CP190 contains the N-terminal BTB/POZ (Broad-Complex, Tramtrack and Bric a brac/POxvirus and Zinc finger) domain and adjacent conserved regions involved in protein interactions. Here, we examined the functional roles of these domains of CP190 in vivo. The best-characterized architectural proteins with insulator functions, Pita, Su(Hw), and dCTCF, interacted predominantly with the BTB domain of CP190. Due to the difficulty of mutating the BTB domain, we obtained a transgenic line expressing a chimeric CP190 with the BTB domain of the human protein Kaiso. Another group of architectural proteins, M1BP, Opbp, and ZIPIC, interacted with one or both of the highly conserved regions in the N-terminal part of CP190. Transgenic lines of D. melanogaster expressing CP190 mutants with a deletion of each of these domains were obtained. The results showed that these mutant proteins only partially compensated for the functions of CP190, weakly binding to selective chromatin sites. Further analysis confirmed the essential role of these domains in recruitment to regulatory regions associated with architectural proteins. We also found that the N-terminal of CP190 was sufficient for recruiting Z4 and Chromator proteins and successfully achieving chromatin opening. Taken together, our results and the results of previous studies showed that the N-terminal region of CP190 is a platform for simultaneous interaction with various DNA-binding architectural proteins and transcription complexes.

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.

CTCF and R-loops are boundaries of cohesin-mediated DNA looping

Cohesin and CCCTC-binding factor (CTCF) are key regulatory proteins of three-dimensional (3D) genome organization. Cohesin extrudes DNA loops that are anchored by CTCF in a polar orientation. Here, we present direct evidence that CTCF binding polarity controls cohesin-mediated DNA looping. Using single-molecule imaging, we demonstrate that a critical N-terminal motif of CTCF blocks cohesin translocation and DNA looping. The cryo-EM structure of the cohesin-CTCF complex reveals that this CTCF motif ahead of zinc fingers can only reach its binding site on the STAG1 cohesin subunit when the N terminus of CTCF faces cohesin. Remarkably, a C-terminally oriented CTCF accelerates DNA compaction by cohesin. DNA-bound Cas9 and Cas12a ribonucleoproteins are also polar cohesin barriers, indicating that stalling may be intrinsic to cohesin itself. Finally, we show that RNA-DNA hybrids (R-loops) block cohesin-mediated DNA compaction in vitro and are enriched with cohesin subunits in vivo, likely forming TAD boundaries.

CapsNetYY1: identifying YY1-mediated chromatin loops based on a capsule network architecture

BACKGROUND

Previous studies have identified that chromosome structure plays a very important role in gene control. The transcription factor Yin Yang 1 (YY1), a multifunctional DNA binding protein, could form a dimer to mediate chromatin loops and active enhancer-promoter interactions. The deletion of YY1 or point mutations at the YY1 binding sites significantly inhibit the enhancer-promoter interactions and affect gene expression. To date, only a few computational methods are available for identifying YY1-mediated chromatin loops.

RESULTS

We proposed a novel model named CapsNetYY1, which was based on capsule network architecture to identify whether a pair of YY1 motifs can form a chromatin loop. Firstly, we encode the DNA sequence using one-hot encoding method. Secondly, multi-scale convolution layer is used to extract local features of the sequence, and bidirectional gated recurrent unit is used to learn the features across time steps. Finally, capsule networks (convolution capsule layer and digital capsule layer) used to extract higher level features and recognize YY1-mediated chromatin loops. Compared with DeepYY1, the only prediction for YY1-mediated chromatin loops, our model CapsNetYY1 achieved the better performance on the independent datasets (AUC [Formula: see text]).

CONCLUSION

The results indicate that CapsNetYY1 is an excellent method for identifying YY1-mediated chromatin loops. We believe that the CapsNetYY1 method will be used for predictive classification of other DNA sequences.

Boundary bypass activity in the abdominal-B region of the Drosophila bithorax complex is position dependent and regulated

Expression of Abdominal-B (Abd-B) in abdominal segments A5-A8 is controlled by four regulatory domains, iab-5-iab-8. Each domain has an initiator element (which sets the activity state), elements that maintain this state and tissue-specific enhancers. To ensure their functional autonomy, each domain is bracketed by boundary elements (Mcp, Fab-7, Fab-7 and Fab-8). In addition to blocking crosstalk between adjacent regulatory domains, the Fab boundaries must also have bypass activity so the relevant regulatory domains can 'jump over' intervening boundaries and activate the Abd-B promoter. In the studies reported here we have investigated the parameters governing bypass activity. We find that the bypass elements in the Fab-7 and Fab-8 boundaries must be located in the regulatory domain that is responsible for driving Abd-B expression. We suggest that bypass activity may also be subject to regulation.

Super-enhancers complexes zoom in transcription in cancer

Super-enhancers (SEs) consist of multiple typical enhancers enriched at high density with transcription factors, histone-modifying enzymes and cofactors. Oncogenic SEs promote tumorigenesis and malignancy by altering protein-coding gene expression and noncoding regulatory element function. Therefore, they play central roles in the treatment of cancer. Here, we review the structural characteristics, organization, identification, and functions of SEs and the underlying molecular mechanism by which SEs drive oncogenic transcription in tumor cells. We then summarize abnormal SE complexes, SE-driven coding genes, and noncoding RNAs involved in tumor development. In summary, we believe that SEs show great potential as biomarkers and therapeutic targets.

Binding by the Polycomb complex component BMI1 and H2A monoubiquitination shape local and long-range interactions in the Arabidopsis genome

Three-dimensional (3D) chromatin organization is highly dynamic during development and seems to play a crucial role in regulating gene expression. Self-interacting domains, commonly called topologically associating domains (TADs) or compartment domains (CDs), have been proposed as the basic structural units of chromatin organization. Surprisingly, although these units have been found in several plant species, they escaped detection in Arabidopsis (Arabidopsis thaliana). Here, we show that the Arabidopsis genome is partitioned into contiguous CDs with different epigenetic features, which are required to maintain appropriate intra-CD and long-range interactions. Consistent with this notion, the histone-modifying Polycomb group machinery is involved in 3D chromatin organization. Yet, while it is clear that Polycomb repressive complex 2 (PRC2)-mediated trimethylation of histone H3 on lysine 27 (H3K27me3) helps establish local and long-range chromatin interactions in plants, the implications of PRC1-mediated histone H2A monoubiquitination on lysine 121 (H2AK121ub) are unclear. We found that PRC1, together with PRC2, maintains intra-CD interactions, but it also hinders the formation of H3K4me3-enriched local chromatin loops when acting independently of PRC2. Moreover, the loss of PRC1 or PRC2 activity differentially affects long-range chromatin interactions, and these 3D changes differentially affect gene expression. Our results suggest that H2AK121ub helps prevent the formation of transposable element/H3K27me1-rich long loops and serves as a docking point for H3K27me3 incorporation.

Boundary Bypass Activity in the Abdominal-B Region of the Drosophila Bithorax Complex is Position Dependent and Regulated

Expression of Abdominal-B ( Abd-B ) in abdominal segments A5 - A8 is controlled by four regulatory domains, iab-5 - iab-8 . Each domain has an initiator element (which sets the activity state), elements that maintain this state and tissue-specific enhancers. To ensure their functional autonomy, each domain is bracketed by boundary elements ( Mcp , Fab-7 , Fab-7 and Fab-8 ). In addition to blocking crosstalk between adjacent regulatory domains, the Fab boundaries must also have bypass activity so the relevant regulatory domains can "jump over" intervening boundaries and activate the Abd-B promoter. In the studies reported here we have investigated the parameters governing bypass activity. We find that the bypass elements in the Fab-7 and Fab-8 boundaries must be located in the regulatory domain that is responsible for driving Abd-B expression. We suggest that bypass activity may also be subject to regulation.

Summary Statement

Boundaries separating Abd-B regulatory domains block crosstalk between domains and mediate their interactions with Abd-B . The latter function is location but not orientation dependent.

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.

Mechanisms of Interaction between Enhancers and Promoters in Three Drosophila Model Systems

In higher eukaryotes, the regulation of developmental gene expression is determined by enhancers, which are often located at a large distance from the promoters they regulate. Therefore, the architecture of chromosomes and the mechanisms that determine the functional interaction between enhancers and promoters are of decisive importance in the development of organisms. Mammals and the model animal Drosophila have homologous key architectural proteins and similar mechanisms in the organization of chromosome architecture. This review describes the current progress in understanding the mechanisms of the formation and regulation of long-range interactions between enhancers and promoters at three well-studied key regulatory loci in Drosophila.

Mechanisms of enhancer-promoter communication and chromosomal architecture in mammals and Drosophila

The spatial organization of chromosomes is involved in regulating the majority of intranuclear processes in higher eukaryotes, including gene expression. Drosophila was used as a model to discover many transcription factors whose homologs play a key role in regulation of gene expression in mammals. According to modern views, a cohesin complex mostly determines the architecture of mammalian chromosomes by forming chromatin loops on anchors created by the CTCF DNA-binding architectural protein. The role of the cohesin complex in chromosome architecture is poorly understood in Drosophila, and CTCF is merely one of many Drosophila architectural proteins with a proven potential to organize specific long-range interactions between regulatory elements in the genome. The review compares the mechanisms responsible for long-range interactions and chromosome architecture between mammals and Drosophila.

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.

CW198 acts as a genetic insulator to block enhancer-promoter interaction in plants

Insulators in vertebrates play a role in genome architecture and orchestrate temporo-spatial enhancer-promoter interactions. In plants, insulators and their associated binding factors have not been documented as of yet, largely as a result of a lack of characterized insulators. In this study, we took a comprehensive strategy to identify and validate the enhancer-blocking insulator CW198. We show that a 1.08-kb CW198 fragment from Arabidopsis can, when interposed between an enhancer and a promoter, efficiently abrogate the activation function of both constitutive and floral organ-specific enhancers in transgenic Arabidopsis and tobacco plants. In plants, both transcriptional crosstalk and spreading of histone modifications were rarely detectable across CW198, which resembles the insulation property observed across the CTCF insulator in the mammalian genome. Taken together, our findings support that CW198 acts as an enhancer-blocking insulator in both Arabidopsis and tobacco. The significance of the present findings and their relevance to the mitigation of mutual interference between enhancers and promoters, as well as multiple promoters in transgenes, is discussed.

A global high-density chromatin interaction network reveals functional long-range and trans-chromosomal relationships

BACKGROUND

Chromatin contacts are essential for gene-expression regulation; however, obtaining a high-resolution genome-wide chromatin contact map is still prohibitively expensive owing to large genome sizes and the quadratic scale of pairwise data. Chromosome conformation capture (3C)-based methods such as Hi-C have been extensively used to obtain chromatin contacts. However, since the sparsity of these maps increases with an increase in genomic distance between contacts, long-range or trans-chromatin contacts are especially challenging to sample.

RESULTS

Here, we create a high-density reference genome-wide chromatin contact map using a meta-analytic approach. We integrate 3600 human, 6700 mouse, and 500 fly Hi-C experiments to create species-specific meta-Hi-C chromatin contact maps with 304 billion, 193 billion, and 19 billion contacts in respective species. We validate that meta-Hi-C contact maps are uniquely powered to capture functional chromatin contacts in both cis and trans. We find that while individual dataset Hi-C networks are largely unable to predict any long-range coexpression (median 0.54 AUC), meta-Hi-C networks perform comparably in both cis and trans (0.65 AUC vs 0.64 AUC). Similarly, for long-range expression quantitative trait loci (eQTL), meta-Hi-C contacts outperform all individual Hi-C experiments, providing an improvement over the conventionally used linear genomic distance-based association. Assessing between species, we find patterns of chromatin contact conservation in both cis and trans and strong associations with coexpression even in species for which Hi-C data is lacking.

CONCLUSIONS

We have generated an integrated chromatin interaction network which complements a large number of methodological and analytic approaches focused on improved specificity or interpretation. This high-depth "super-experiment" is surprisingly powerful in capturing long-range functional relationships of chromatin interactions, which are now able to predict coexpression, eQTLs, and cross-species relationships. The meta-Hi-C networks are available at https://labshare.cshl.edu/shares/gillislab/resource/HiC/ .

The 3D architecture of the pepper genome and its relationship to function and evolution

The organization of chromatin into self-interacting domains is universal among eukaryotic genomes, though how and why they form varies considerably. Here we report a chromosome-scale reference genome assembly of pepper (Capsicum annuum) and explore its 3D organization through integrating high-resolution Hi-C maps with epigenomic, transcriptomic, and genetic variation data. Chromatin folding domains in pepper are as prominent as TADs in mammals but exhibit unique characteristics. They tend to coincide with heterochromatic regions enriched with retrotransposons and are frequently embedded in loops, which may correlate with transcription factories. Their boundaries are hotspots for chromosome rearrangements but are otherwise depleted for genetic variation. While chromatin conformation broadly affects transcription variance, it does not predict differential gene expression between tissues. Our results suggest that pepper genome organization is explained by a model of heterochromatin-driven folding promoted by transcription factories and that such spatial architecture is under structural and functional constraints.

The organization of chromatin into self-interacting domains is universal among eukaryotic genomes. Here, the authors report a reference-grade pepper genome assembly and use this reference to help describe the relationship among 3D chromatin conformation, chromatin function, and gene expression.

Decoding the PITX2-controlled genetic network in atrial fibrillation

Atrial fibrillation (AF), the most common sustained cardiac arrhythmia and a major risk factor for stroke, often arises through ectopic electrical impulses derived from the pulmonary veins (PVs). Sequence variants in enhancers controlling expression of the transcription factor PITX2, which is expressed in the cardiomyocytes (CMs) of the PV and left atrium (LA), have been implicated in AF predisposition. Single nuclei multiomic profiling of RNA and analysis of chromatin accessibility combined with spectral clustering uncovered distinct PV- and LA-enriched CM cell states. Pitx2-mutant PV and LA CMs exhibited gene expression changes consistent with cardiac dysfunction through cell type-distinct, PITX2-directed, cis-regulatory grammars controlling target gene expression. The perturbed network targets in each CM were enriched in distinct human AF predisposition genes, suggesting combinatorial risk for AF genesis. Our data further reveal that PV and LA Pitx2-mutant CMs signal to endothelial and endocardial cells through BMP10 signaling with pathogenic potential. This work provides a multiomic framework for interrogating the basis of AF predisposition in the PVs of humans.

3D or Not 3D: Shaping the Genome during Development

One of the most fundamental questions in developmental biology is how one fertilized cell can give rise to a fully mature organism and how gene regulation governs this process. Precise spatiotemporal gene expression is required for development and is believed to be achieved through a complex interplay of sequence-specific information, epigenetic modifications, trans-acting factors, and chromatin folding. Here we review the role of chromatin folding during development, the mechanisms governing 3D genome organization, and how it is established in the embryo. Furthermore, we discuss recent advances and debated questions regarding the contribution of the 3D genome to gene regulation during organogenesis. Finally, we describe the mechanisms that can reshape the 3D genome, including disease-causing structural variations and the emerging view that transposable elements contribute to chromatin organization.

Single-Molecule Micromanipulation and Super-Resolution Imaging Resolve Nanodomains Underlying Chromatin Folding in Mitotic Chromosomes

The folding of interphase chromatin into highly compact mitotic chromosomes is one of the most recognizable changes during the cell cycle. However, the structural organization underlying this drastic compaction remains elusive. Here, we combine several super resolution methods, including structured illumination microscopy (SIM), binding-activated localization microscopy (BALM), and atomic force microscopy (AFM), to examine the structural details of the DNA within the mitotic chromosome, both in the native state and after up to 30-fold extension using single-molecule micromanipulation. Images of native chromosomes reveal widespread ∼125 nm compact granules (CGs) throughout the metaphase chromosome. However, at maximal extensions, we find exclusively ∼90 nm domains (mitotic nanodomains, MNDs) that are unexpectedly resistant to extensive forces of tens of nanonewtons. The DNA content of the MNDs is estimated to be predominantly ∼80 kb, which is comparable to the size of the inner loops predicted by a recent nested loop model of the mitotic chromosome. With this DNA content, the total volume expected of the human genome assuming closely packed MNDs is nearly identical to what is observed. Thus, altogether, these results suggest that these mechanically stable MNDs, and their higher-order assembly into CGs, are the dominant higher-level structures that underlie the compaction of chromatin from interphase to metaphase.

Anopheles mosquitoes reveal new principles of 3D genome organization in insects

Chromosomes are hierarchically folded within cell nuclei into territories, domains and subdomains, but the functional importance and evolutionary dynamics of these hierarchies are poorly defined. Here, we comprehensively profile genome organizations of five Anopheles mosquito species and show how different levels of chromatin architecture influence each other. Patterns observed on Hi-C maps are associated with known cytological structures, epigenetic profiles, and gene expression levels. Evolutionary analysis reveals conservation of chromatin architecture within synteny blocks for tens of millions of years and enrichment of synteny breakpoints in regions with increased genomic insulation. However, in-depth analysis shows a confounding effect of gene density on both insulation and distribution of synteny breakpoints, suggesting limited causal relationship between breakpoints and regions with increased genomic insulation. At the level of individual loci, we identify specific, extremely long-ranged looping interactions, conserved for ~100 million years. We demonstrate that the mechanisms underlying these looping contacts differ from previously described Polycomb-dependent interactions and clustering of active chromatin.

Anopheles mosquitoes are vectors of human malaria, and better understanding of them has implications for public health. Here, the authors apply Hi-C, FISH, RNA-seq, and ChIP-seq techniques to comprehensively characterize chromatin architecture and its evolutionary dynamics in five Anopheles species.

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.

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.

Understanding 3D Genome Organization and Its Effect on Transcriptional Gene Regulation Under Environmental Stress in Plant: A Chromatin Perspective

The genome of a eukaryotic organism is comprised of a supra-molecular complex of chromatin fibers and intricately folded three-dimensional (3D) structures. Chromosomal interactions and topological changes in response to the developmental and/or environmental stimuli affect gene expression. Chromatin architecture plays important roles in DNA replication, gene expression, and genome integrity. Higher-order chromatin organizations like chromosome territories (CTs), A/B compartments, topologically associating domains (TADs), and chromatin loops vary among cells, tissues, and species depending on the developmental stage and/or environmental conditions (4D genomics). Every chromosome occupies a separate territory in the interphase nucleus and forms the top layer of hierarchical structure (CTs) in most of the eukaryotes. While the A and B compartments are associated with active (euchromatic) and inactive (heterochromatic) chromatin, respectively, having well-defined genomic/epigenomic features, TADs are the structural units of chromatin. Chromatin architecture like TADs as well as the local interactions between promoter and regulatory elements correlates with the chromatin activity, which alters during environmental stresses due to relocalization of the architectural proteins. Moreover, chromatin looping brings the gene and regulatory elements in close proximity for interactions. The intricate relationship between nucleotide sequence and chromatin architecture requires a more comprehensive understanding to unravel the genome organization and genetic plasticity. During the last decade, advances in chromatin conformation capture techniques for unravelling 3D genome organizations have improved our understanding of genome biology. However, the recent advances, such as Hi-C and ChIA-PET, have substantially increased the resolution, throughput as well our interest in analysing genome organizations. The present review provides an overview of the historical and contemporary perspectives of chromosome conformation capture technologies, their applications in functional genomics, and the constraints in predicting 3D genome organization. We also discuss the future perspectives of understanding high-order chromatin organizations in deciphering transcriptional regulation of gene expression under environmental stress (4D genomics). These might help design the climate-smart crop to meet the ever-growing demands of food, feed, and fodder.

HiCRes: a computational method to estimate and predict the genomic resolution of Hi-C libraries

Three-dimensional (3D) conformation of the chromatin is crucial to stringently regulate gene expression patterns and DNA replication in a cell-type specific manner. Hi-C is a key technique for measuring 3D chromatin interactions genome wide. Estimating and predicting the resolution of a library is an essential step in any Hi-C experimental design. Here, we present the mathematical concepts to estimate the resolution of a dataset and predict whether deeper sequencing would enhance the resolution. We have developed HiCRes, a docker pipeline, by applying these concepts to several Hi-C libraries.

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.

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.

Mechanisms of CP190 Interaction with Architectural Proteins in Drosophila Melanogaster

Most of the known Drosophila architectural proteins interact with an important cofactor, CP190, that contains three domains (BTB, M, and D) that are involved in protein–protein interactions. The highly conserved N-terminal CP190 BTB domain forms a stable homodimer that interacts with unstructured regions in the three best-characterized architectural proteins: dCTCF, Su(Hw), and Pita. Here, we identified two new CP190 partners, CG4730 and CG31365, that interact with the BTB domain. The CP190 BTB resembles the previously characterized human BCL6 BTB domain, which uses its hydrophobic groove to specifically associate with unstructured regions of several transcriptional repressors. Using GST pull-down and yeast two-hybrid assays, we demonstrated that mutations in the hydrophobic groove strongly affect the affinity of CP190 BTB for the architectural proteins. In the yeast two-hybrid assay, we found that architectural proteins use various mechanisms to improve the efficiency of interaction with CP190. Pita and Su(Hw) have two unstructured regions that appear to simultaneously interact with hydrophobic grooves in the BTB dimer. In dCTCF and CG31365, two adjacent regions interact simultaneously with the hydrophobic groove of the BTB and the M domain of CP190. Finally, CG4730 interacts with the BTB, M, and D domains of CP190 simultaneously. These results suggest that architectural proteins use different mechanisms to increase the efficiency of interaction with CP190.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Monocytic THP-1 cells diverge significantly from their primary counterparts: a comparative examination of the chromosomal conformations and transcriptomes

Immortalized cell lines have long been used as model systems to systematically investigate biological processes under controlled and reproducible conditions, providing insights that have greatly advanced cellular biology and medical sciences. Recently, the widely used monocytic leukemia cell line, THP-1, was comprehensively examined to understand mechanistic relationships between the 3D chromatin structure and transcription during the trans-differentiation of monocytes to macrophages. To corroborate these observations in primary cells, we analyze in situ Hi-C and RNA-seq data of human primary monocytes and their differentiated macrophages in comparison to that obtained from the monocytic/macrophagic THP-1 cells. Surprisingly, we find significant differences between the primary cells and the THP-1 cells at all levels of chromatin structure, from loops to topologically associated domains to compartments. Importantly, the compartment-level differences correlate significantly with transcription: those genes that are in A-compartments in the primary cells but are in B-compartments in the THP-1 cells exhibit a higher level of expression in the primary cells than in the THP-1 cells, and vice versa. Overall, the genes in these different compartments are enriched for a wide range of pathways, and, at least in the case of the monocytic cells, their altered expression in certain pathways in the THP-1 cells argues for a less immune cell-like phenotype, suggesting that immortalization or prolonged culturing of THP-1 caused a divergence of these cells from their primary counterparts. It is thus essential to reexamine phenotypic details observed in cell lines with their primary counterparts so as to ensure a proper understanding of functional cell states in vivo.

FISH Going Meso-Scale: A Microscopic Search for Chromatin Domains

The intimate relationships between genome structure and function direct efforts toward deciphering three-dimensional chromatin organization within the interphase nuclei at different genomic length scales. For decades, major insights into chromatin structure at the level of large-scale euchromatin and heterochromatin compartments, chromosome territories, and subchromosomal regions resulted from the evolution of light microscopy and fluorescence in situ hybridization. Studies of nanoscale nucleosomal chromatin organization benefited from a variety of electron microscopy techniques. Recent breakthroughs in the investigation of mesoscale chromatin structures have emerged from chromatin conformation capture methods (C-methods). Chromatin has been found to form hierarchical domains with high frequency of local interactions from loop domains to topologically associating domains and compartments. During the last decade, advances in super-resolution light microscopy made these levels of chromatin folding amenable for microscopic examination. Here we are reviewing recent developments in FISH-based approaches for detection, quantitative measurements, and validation of contact chromatin domains deduced from C-based data. We specifically focus on the design and application of Oligopaint probes, which marked the latest progress in the imaging of chromatin domains. Vivid examples of chromatin domain FISH-visualization by means of conventional, super-resolution light and electron microscopy in different model organisms are provided.

Genomic organization of the autonomous regulatory domain of eyeless locus in Drosophila melanogaster

In Drosophila, expression of eyeless (ey) gene is restricted to the developing eyes and central nervous system. However, the flanking genes, myoglianin (myo), and bent (bt) have different temporal and spatial expression patterns as compared to the ey. How distinct regulation of ey is maintained is mostly unknown. Earlier, we have identified a boundary element intervening myo and ey genes (ME boundary) that prevents the crosstalk between the cis-regulatory elements of myo and ey genes. In the present study, we further searched for the cis-elements that define the domain of ey and maintain its expression pattern. We identify another boundary element between ey and bt, the EB boundary. The EB boundary separates the regulatory landscapes of ey and bt genes. The two boundaries, ME and EB, show a long-range interaction as well as interact with the nuclear architecture. This suggests functional autonomy of the ey locus and its insulation from differentially regulated flanking regions. We also identify a new Polycomb Response Element, the ey-PRE, within the ey domain. The expression state of the ey gene, once established during early development is likely to be maintained with the help of ey-PRE. Our study proposes a general regulatory mechanism by which a gene can be maintained in a functionally independent chromatin domain in gene-rich euchromatin.

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.

Chromatin Organization and Function in Drosophila

Eukaryotic genomes are packaged into high-order chromatin structures organized in discrete territories inside the cell nucleus, which is surrounded by the nuclear envelope acting as a barrier. This chromatin organization is complex and dynamic and, thus, determining the spatial and temporal distribution and folding of chromosomes within the nucleus is critical for understanding the role of chromatin topology in genome function. Primarily focusing on the regulation of gene expression, we review here how the genome of Drosophila melanogaster is organized into the cell nucleus, from small scale histone–DNA interactions to chromosome and lamina interactions in the nuclear space.

Seeing the forest through the trees: prioritising potentially functional interactions from Hi-C

Eukaryotic genomes are highly organised within the nucleus of a cell, allowing widely dispersed regulatory elements such as enhancers to interact with gene promoters through physical contacts in three-dimensional space. Recent chromosome conformation capture methodologies such as Hi-C have enabled the analysis of interacting regions of the genome providing a valuable insight into the three-dimensional organisation of the chromatin in the nucleus, including chromosome compartmentalisation and gene expression. Complicating the analysis of Hi-C data, however, is the massive amount of identified interactions, many of which do not directly drive gene function, thus hindering the identification of potentially biologically functional 3D interactions. In this review, we collate and examine the downstream analysis of Hi-C data with particular focus on methods that prioritise potentially functional interactions. We classify three groups of approaches: structural-based discovery methods, e.g. A/B compartments and topologically associated domains, detection of statistically significant chromatin interactions, and the use of epigenomic data integration to narrow down useful interaction information. Careful use of these three approaches is crucial to successfully identifying potentially functional interactions within the genome.

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.

Heterochromatin Networks: Topology, Dynamics, and Function (a Working Hypothesis)

Open systems can only exist by self-organization as pulsing structures exchanging matter and energy with the outer world. This review is an attempt to reveal the organizational principles of the heterochromatin supra-intra-chromosomal network in terms of nonlinear thermodynamics. The accessibility of the linear information of the genetic code is regulated by constitutive heterochromatin (CHR) creating the positional information in a system of coordinates. These features include scale-free splitting-fusing of CHR with the boundary constraints of the nucleolus and nuclear envelope. The analysis of both the literature and our own data suggests a radial-concentric network as the main structural organization principle of CHR regulating transcriptional pulsing. The dynamic CHR network is likely created together with nucleolus-associated chromatin domains, while the alveoli of this network, including springy splicing speckles, are the pulsing transcription hubs. CHR contributes to this regulation due to the silencing position variegation effect, stickiness, and flexible rigidity determined by the positioning of nucleosomes. The whole system acts in concert with the elastic nuclear actomyosin network which also emerges by self-organization during the transcriptional pulsing process. We hypothesize that the the transcriptional pulsing, in turn, adjusts its frequency/amplitudes specified by topologically associating domains to the replication timing code that determines epigenetic differentiation memory.

Plant 3D genomics: the exploration and application of chromatin organization

Eukaryotic genomes are highly folded for packing into higher-order chromatin structures in the nucleus. With the emergence of state-of-the-art chromosome conformation capture methods and microscopic imaging techniques, the spatial organization of chromatin and its functional implications have been interrogated. Our knowledge of 3D chromatin organization in plants has improved dramatically in the past few years, building on the early advances in animal systems. Here, we review recent advances in 3D genome mapping approaches, our understanding of the sophisticated organization of spatial structures, and the application of 3D genomic principles in plants. We also discuss directions for future developments in 3D genomics in plants.

Chromatin architectural proteins regulate flowering time by precluding gene looping

Chromatin structure is critical for gene expression and many other cellular processes. In Arabidopsis thaliana, the floral repressor FLC adopts a self-loop chromatin structure via bridging of its flanking regions. This local gene loop is necessary for active FLC expression. However, the molecular mechanism underlying the formation of this class of gene loops is unknown. Here, we report the characterization of a group of linker histone-like proteins, named the GH1-HMGA family in Arabidopsis, which act as chromatin architecture modulators. We demonstrate that these family members redundantly promote the floral transition through the repression of FLC A genome-wide study revealed that this family preferentially binds to the 5' and 3' ends of gene bodies. The loss of this binding increases FLC expression by stabilizing the FLC 5' to 3' gene looping. Our study provides mechanistic insights into how a family of evolutionarily conserved proteins regulates the formation of local gene loops.

Multi-Scale Organization of the Drosophila melanogaster Genome

Interphase chromatin, despite its appearance, is a highly organized framework of loops and bends. Chromosomes are folded into topologically associating domains, or TADs, and each chromosome and its homolog occupy a distinct territory within the nucleus. In Drosophila, genome organization is exceptional because homologous chromosome pairing is in both germline and somatic tissues, which promote interhomolog interactions such as transvection that can affect gene expression in trans. In this review, we focus on what is known about genome organization in Drosophila and discuss it from TADs to territory. We start by examining intrachromosomal organization at the sub-chromosome level into TADs, followed by a comprehensive analysis of the known proteins that play a key role in TAD formation and boundary establishment. We then zoom out to examine interhomolog interactions such as pairing and transvection that are abundant in Drosophila but rare in other model systems. Finally, we discuss chromosome territories that form within the nucleus, resulting in a complete picture of the multi-scale organization of the Drosophila genome.

The two waves in single-cell 3D genomics

For decades, biochemical methods for the analysis of genome structure and function provided cell-population-averaged data that allowed general principles and tendencies to be disclosed. Microscopy-based studies, which immanently involve single-cell analysis, did not provide sufficient spatial resolution to investigate the particularly small details of 3D genome folding. Nevertheless, these studies demonstrated that mutual positions of chromosome territories within cell nuclei and individual genomic loci within chromosomal territories can vary significantly in individual cells. The development of new technologies in biochemistry and the advent of super-resolution microscopy in the last decade have made possible the full-scale study of 3D genome organization in individual cells. Maps of the 3D genome build based on C-data and super-resolution microscopy are highly consistent and, therefore, biologically relevant. The internal structures of individual chromosomes, loci, and topologically associating domains (TADs) are resolved as well as cell-cycle dynamics. 3D modeling allows one to investigate the physical mechanisms underlying genome folding. Finally, joint profiling of genome topology and epigenetic features will allow 3D genomics to handle complex cell-to-cell heterogeneity. In this review, we summarize the present state of studies into 3D genome organization in individual cells, analyze the technical problems of single-cell studies, and outline perspectives of 3D genomics.