MCM complexes are barriers that restrict cohesin-mediated loop extrusion.

Eukaryotic genomes are compacted into loops and topologically associating domains (TADs)1-3, which contribute to transcription, recombination and genomic stability4,5. Cohesin extrudes DNA into loops that are thought to lengthen until CTCF boundaries are encountered6-12. Little is known about whether loop extrusion is impeded by DNA-bound machines. Here we show that the minichromosome maintenance (MCM) complex is a barrier that restricts loop extrusion in G1 phase. Single-nucleus Hi-C (high-resolution chromosome conformation capture) of mouse zygotes reveals that MCM loading reduces CTCF-anchored loops and decreases TAD boundary insulation, which suggests that loop extrusion is impeded before reaching CTCF. This effect extends to HCT116 cells, in which MCMs affect the number of CTCF-anchored loops and gene expression. Simulations suggest that MCMs are abundant, randomly positioned and partially permeable barriers. Single-molecule imaging shows that MCMs are physical barriers that frequently constrain cohesin translocation in vitro. Notably, chimeric yeast MCMs that contain a cohesin-interaction motif from human MCM3 induce cohesin pausing, indicating that MCMs are 'active' barriers with binding sites. These findings raise the possibility that cohesin can arrive by loop extrusion at MCMs, which determine the genomic sites at which sister chromatid cohesion is established. On the basis of in vivo, in silico and in vitro data, we conclude that distinct loop extrusion barriers shape the three-dimensional genome.

A central role for canonical PRC1 in shaping the 3D nuclear landscape.

Polycomb group (PcG) proteins silence gene expression by chemically and physically modifying chromatin. A subset of PcG target loci are compacted and cluster in the nucleus; a conformation that is thought to contribute to gene silencing. However, how these interactions influence gross nuclear organization and their relationship with transcription remains poorly understood. Here we examine the role of Polycomb-repressive complex 1 (PRC1) in shaping 3D genome organization in mouse embryonic stem cells (mESCs). Using a combination of imaging and Hi-C analyses, we show that PRC1-mediated long-range interactions are independent of CTCF and can bridge sites at a megabase scale. Impairment of PRC1 enzymatic activity does not directly disrupt these interactions. We demonstrate that PcG targets coalesce in vivo, and that developmentally induced expression of one of the target loci disrupts this spatial arrangement. Finally, we show that transcriptional activation and the loss of PRC1-mediated interactions are separable events. These findings provide important insights into the function of PRC1, while highlighting the complexity of this regulatory system.

Ultra-long-range interactions between active regulatory elements.

Contacts between enhancers and promoters are thought to relate to their ability to activate transcription. Investigating factors that contribute to such chromatin interactions is therefore important for understanding gene regulation. Here, we have determined contact frequencies between millions of pairs of cis-regulatory elements from chromosome conformation capture data sets and analyzed a collection of hundreds of DNA-binding factors for binding at regions of enriched contacts. This analysis revealed enriched contacts at sites bound by many factors associated with active transcription. We show that active regulatory elements, independent of cohesin and polycomb, interact with each other across distances of tens of megabases in vertebrate and invertebrate genomes and that interactions correlate and change with activity. However, these ultra-long-range interactions are not dependent on RNA polymerase II transcription or individual transcription cofactors. Using simulations, we show that a model of chromatin and multivalent binding factors can give rise to long-range interactions via bridging-induced clustering. We propose that long-range interactions between cis-regulatory elements are driven by at least three distinct processes: cohesin-mediated loop extrusion, polycomb contacts, and clustering of active regions.

Bioframe: operations on genomic intervals in Pandas dataframes.

MOTIVATION: Genomic intervals are one of the most prevalent data structures in computational genome biology, and used to represent features ranging from genes, to DNA binding sites, to disease variants. Operations on genomic intervals provide a language for asking questions about relationships between features. While there are excellent interval arithmetic tools for the command line, they are not smoothly integrated into Python, one of the most popular general-purpose computational and visualization environments. RESULTS: Bioframe is a library to enable flexible and performant operations on genomic interval dataframes in Python. Bioframe extends the Python data science stack to use cases for computational genome biology by building directly on top of two of the most commonly-used Python libraries, NumPy and Pandas. The bioframe API enables flexible name and column orders, and decouples operations from data formats to avoid unnecessary conversions, a common scourge for bioinformaticians. Bioframe achieves these goals while maintaining high performance and a rich set of features. AVAILABILITY AND IMPLEMENTATION: Bioframe is open-source under MIT license, cross-platform, and can be installed from the Python Package Index. The source code is maintained by Open2C on GitHub at https://github.com/open2c/bioframe.

Pairtools: From sequencing data to chromosome contacts.

The field of 3D genome organization produces large amounts of sequencing data from Hi-C and a rapidly-expanding set of other chromosome conformation protocols (3C+). Massive and heterogeneous 3C+ data require high-performance and flexible processing of sequenced reads into contact pairs. To meet these challenges, we present pairtools-a flexible suite of tools for contact extraction from sequencing data. Pairtools provides modular command-line interface (CLI) tools that can be flexibly chained into data processing pipelines. The core operations provided by pairtools are parsing of.sam alignments into Hi-C pairs, sorting and removal of PCR duplicates. In addition, pairtools provides auxiliary tools for building feature-rich 3C+ pipelines, including contact pair manipulation, filtration, and quality control. Benchmarking pairtools against popular 3C+ data pipelines shows advantages of pairtools for high-performance and flexible 3C+ analysis. Finally, pairtools provides protocol-specific tools for restriction-based protocols, haplotype-resolved contacts, and single-cell Hi-C. The combination of CLI tools and tight integration with Python data analysis libraries makes pairtools a versatile foundation for a broad range of 3C+ pipelines.

Cooltools: Enabling high-resolution Hi-C analysis in Python.

Chromosome conformation capture (3C) technologies reveal the incredible complexity of genome organization. Maps of increasing size, depth, and resolution are now used to probe genome architecture across cell states, types, and organisms. Larger datasets add challenges at each step of computational analysis, from storage and memory constraints to researchers' time; however, analysis tools that meet these increased resource demands have not kept pace. Furthermore, existing tools offer limited support for customizing analysis for specific use cases or new biology. Here we introduce cooltools (https://github.com/open2c/cooltools), a suite of computational tools that enables flexible, scalable, and reproducible analysis of high-resolution contact frequency data. Cooltools leverages the widely-adopted cooler format which handles storage and access for high-resolution datasets. Cooltools provides a paired command line interface (CLI) and Python application programming interface (API), which respectively facilitate workflows on high-performance computing clusters and in interactive analysis environments. In short, cooltools enables the effective use of the latest and largest genome folding datasets.

Single-nucleus Hi-C reveals unique chromatin reorganization at oocyte-to-zygote transition.

Chromatin is reprogrammed after fertilization to produce a totipotent zygote with the potential to generate a new organism. The maternal genome inherited from the oocyte and the paternal genome provided by sperm coexist as separate haploid nuclei in the zygote. How these two epigenetically distinct genomes are spatially organized is poorly understood. Existing chromosome conformation capture-based methods are not applicable to oocytes and zygotes owing to a paucity of material. To study three-dimensional chromatin organization in rare cell types, we developed a single-nucleus Hi-C (high-resolution chromosome conformation capture) protocol that provides greater than tenfold more contacts per cell than the previous method. Here we show that chromatin architecture is uniquely reorganized during the oocyte-to-zygote transition in mice and is distinct in paternal and maternal nuclei within single-cell zygotes. Features of genomic organization including compartments, topologically associating domains (TADs) and loops are present in individual oocytes when averaged over the genome, but the presence of each feature at a locus varies between cells. At the sub-megabase level, we observed stochastic clusters of contacts that can occur across TAD boundaries but average into TADs. Notably, we found that TADs and loops, but not compartments, are present in zygotic maternal chromatin, suggesting that these are generated by different mechanisms. Our results demonstrate that the global chromatin organization of zygote nuclei is fundamentally different from that of other interphase cells. An understanding of this zygotic chromatin 'ground state' could potentially provide insights into reprogramming cells to a state of totipotency.

Trisomies Reorganize Human 3D Genome.

Trisomy is the presence of one extra copy of an entire chromosome or its part in a cell nucleus. In humans, autosomal trisomies are associated with severe developmental abnormalities leading to embryonic lethality, miscarriage or pronounced deviations of various organs and systems at birth. Trisomies are characterized by alterations in gene expression level, not exclusively on the trisomic chromosome, but throughout the genome. Here, we applied the high-throughput chromosome conformation capture technique (Hi-C) to study chromatin 3D structure in human chorion cells carrying either additional chromosome 13 (Patau syndrome) or chromosome 16 and in cultured fibroblasts with extra chromosome 18 (Edwards syndrome). The presence of extra chromosomes results in systematic changes of contact frequencies between small and large chromosomes. Analyzing the behavior of individual chromosomes, we found that a limited number of chromosomes change their contact patterns stochastically in trisomic cells and that it could be associated with lamina-associated domains (LAD) and gene content. For trisomy 13 and 18, but not for trisomy 16, the proportion of compacted loci on a chromosome is correlated with LAD content. We also found that regions of the genome that become more compact in trisomic cells are enriched in housekeeping genes, indicating a possible decrease in chromatin accessibility and transcription level of these genes. These results provide a framework for understanding the mechanisms of pan-genome transcription dysregulation in trisomies in the context of chromatin spatial organization.

Developmentally regulated Shh expression is robust to TAD perturbations.

Topologically associating domains (TADs) have been proposed to both guide and constrain enhancer activity. Shh is located within a TAD known to contain all its enhancers. To investigate the importance of chromatin conformation and TAD integrity on developmental gene regulation, we have manipulated the Shh TAD - creating internal deletions, deleting CTCF sites, and deleting and inverting sequences at TAD boundaries. Chromosome conformation capture and fluorescence in situ hybridisation assays were used to investigate the changes in chromatin conformation that result from these manipulations. Our data suggest that these substantial alterations in TAD structure have no readily detectable effect on Shh expression patterns or levels of Shh expression during development - except where enhancers are deleted - and result in no detectable phenotypes. Only in the case of a larger deletion at one TAD boundary could ectopic influence of the Shh limb enhancer be detected on a gene (Mnx1) in the neighbouring TAD. Our data suggests that, contrary to expectations, the developmental regulation of Shh expression is remarkably robust to TAD perturbations.

A mechanism of cohesin-dependent loop extrusion organizes zygotic genome architecture.

Fertilization triggers assembly of higher-order chromatin structure from a condensed maternal and a naïve paternal genome to generate a totipotent embryo. Chromatin loops and domains have been detected in mouse zygotes by single-nucleus Hi-C (snHi-C), but not bulk Hi-C. It is therefore unclear when and how embryonic chromatin conformations are assembled. Here, we investigated whether a mechanism of cohesin-dependent loop extrusion generates higher-order chromatin structures within the one-cell embryo. Using snHi-C of mouse knockout embryos, we demonstrate that the zygotic genome folds into loops and domains that critically depend on Scc1-cohesin and that are regulated in size and linear density by Wapl. Remarkably, we discovered distinct effects on maternal and paternal chromatin loop sizes, likely reflecting differences in loop extrusion dynamics and epigenetic reprogramming. Dynamic polymer models of chromosomes reproduce changes in snHi-C, suggesting a mechanism where cohesin locally compacts chromatin by active loop extrusion, whose processivity is controlled by Wapl. Our simulations and experimental data provide evidence that cohesin-dependent loop extrusion organizes mammalian genomes over multiple scales from the one-cell embryo onward.

Coolpup.py: versatile pile-up analysis of Hi-C data.

MOTIVATION: Hi-C is currently the method of choice to investigate the global 3D organization of the genome. A major limitation of Hi-C is the sequencing depth required to robustly detect loops in the data. A popular approach used to mitigate this issue, even in single-cell Hi-C data, is genome-wide averaging (piling-up) of peaks, or other features, annotated in high-resolution datasets, to measure their prominence in less deeply sequenced data. However, current tools do not provide a computationally efficient and versatile implementation of this approach. RESULTS: Here, we describe coolpup.py-a versatile tool to perform pile-up analysis on Hi-C data. We demonstrate its utility by replicating previously published findings regarding the role of cohesin and CTCF in 3D genome organization, as well as discovering novel details of Polycomb-driven interactions. We also present a novel variation of the pile-up approach that can aid the statistical analysis of looping interactions. We anticipate that coolpup.py will aid in Hi-C data analysis by allowing easy to use, versatile and efficient generation of pile-ups. AVAILABILITY AND IMPLEMENTATION: Coolpup.py is cross-platform, open-source and free (MIT licensed) software. Source code is available from https://github.com/Phlya/coolpuppy and it can be installed from the Python Packaging Index.

C-TALE, a new cost-effective method for targeted enrichment of Hi-C/3C-seq libraries.

Studies performed using Hi-C and other high-throughput whole-genome C-methods have demonstrated that 3D organization of eukaryotic genomes is functionally relevant. Unfortunately, ultra-deep sequencing of Hi-C libraries necessary to detect loop structures in large vertebrate genomes remains rather expensive. However, many studies are in fact aimed at determining the fine-scale 3D structure of comparatively small genomic regions up to several Mb in length. Such studies typically focus on the spatial structure of domains of coregulated genes, molecular mechanisms of loop formation, and interrogation of functional significance of GWAS-revealed polymorphisms. Therefore, a handful of molecular techniques based on Hi-C have been developed to address such issues. These techniques commonly rely on in-solution hybridization of Hi-C/3C-seq libraries with pools of biotinylated baits covering the region of interest, followed by deep sequencing of the enriched library. Here, we describe a new protocol of this kind, C-TALE (Chromatin TArget Ligation Enrichment). Preparation of hybridization probes from bacterial artificial chromosomes and an additional round of enrichment make C-TALE a cost-effective alternative to existing many-versus-all C-methods.

Bivalent promoter hypermethylation in cancer is linked to the H327me3/H3K4me3 ratio in embryonic stem cells.

BACKGROUND: Thousands of mammalian promoters are defined by co-enrichment of the histone tail modifications H3K27me3 (repressive) and H3K4me3 (activating) and are thus termed bivalent. It was previously observed that bivalent genes in human ES cells (hESC) are frequent targets for hypermethylation in human cancers, and depletion of DNA methylation in mouse embryonic stem cells has a marked impact on H3K27me3 distribution at bivalent promoters. However, only a fraction of bivalent genes in stem cells are targets of hypermethylation in cancer, and it is currently unclear whether all bivalent promoters are equally sensitive to DNA hypomethylation and whether H3K4me3 levels play a role in the interplay between DNA methylation and H3K27me3. RESULTS: We report the sub-classification of bivalent promoters into two groups-promoters with a high H3K27me3:H3K4me3 (hiBiv) ratio or promoters with a low H3K27me3:H3K4me3 ratio (loBiv). HiBiv are enriched in canonical Polycomb components, show a higher degree of local intrachromosomal contacts and are highly sensitive to DNA hypomethylation in terms of H3K27me3 depletion from broad Polycomb domains. In contrast, loBiv promoters are enriched in non-canonical Polycomb components, show lower intrachromosomal contacts and are less sensitive to DNA hypomethylation at the same genomic resolution. Multiple systems reveal that hiBiv promoters are more depleted of Polycomb complexes than loBiv promoters following a reduction in DNA methylation, and we demonstrate that H3K27me3 re-accumulates at promoters when DNA methylation is restored. In human cancer, we show that hiBiv promoters lose H3K27me3 and are more susceptible to DNA hypermethylation than loBiv promoters. CONCLUSION: We conclude that bivalency as a general term to describe mammalian promoters is an over-simplification and our sub-classification has revealed novel insights into the interplay between the largely antagonistic presence of DNA methylation and Polycomb systems at bivalent promoters. This approach redefines molecular pathologies underlying disease in which global DNA methylation is aberrant or where Polycomb mutations are present.

Low Level of Expression of C-Terminally Truncated Human FUS Causes Extensive Changes in the Spinal Cord Transcriptome of Asymptomatic Transgenic Mice.

A number of mutations in a gene encoding RNA-binding protein FUS have been linked to the development of a familial form of amyotrophic lateral sclerosis known as FUS-ALS. C-terminal truncations of FUS by either nonsense or frameshift mutations lead to the development of FUS-ALS with a particularly early onset and fast progression. However, even in patients bearing these highly pathogenic mutations the function of motor neurons is not noticeably compromised for at least a couple of decades, suggesting that until cytoplasmic levels of FUS lacking its C-terminal nuclear localisation signal reaches a critical threshold, motor neurons are able to tolerate its permanent production. In order to identify how the nervous system responds to low levels of pathogenic variants of FUS we produced and characterised a mouse line, L-FUS[1-359], with a low neuronal expression level of a highly aggregation-prone and pathogenic form of C-terminally truncated FUS. In contrast to mice that express substantially higher level of the same FUS variant and develop severe early onset motor neuron pathology, L-FUS[1-359] mice do not develop any clinical or histopathological signs of motor neuron deficiency even at old age. Nevertheless, we detected substantial changes in the spinal cord transcriptome of these mice compared to their wild type littermates. We suggest that at least some of these changes reflect activation of cellular mechanisms compensating for the potentially damaging effect of pathogenic FUS production. Further studies of these mechanism might reveal effective targets for therapy of FUS-ALS and possibly, other forms of ALS.

Quantitative differences in TAD border strength underly the TAD hierarchy in Drosophila chromosomes.

Chromosomes in many organisms, including Drosophila and mammals, are folded into topologically associating domains (TADs). Increasing evidence suggests that TAD folding is hierarchical, wherein subdomains combine to form larger superdomains, instead of a sequence of nonoverlapping domains. Here, we studied the hierarchical structure of TADs in Drosophila. We show that the boundaries of TADs of different hierarchical levels are characterized by the presence of different portions of active chromatin, but do not vary in the binding of architectural proteins, such as CCCTC binding factor or cohesin. The apparent hierarchy of TADs in Drosophila chromosomes is not likely to have functional importance but rather reflects various options of long-range chromatin folding directed by the distribution of active and inactive chromatin segments and may represent population average.

Active chromatin and transcription play a key role in chromosome partitioning into topologically associating domains.

Recent advances enabled by the Hi-C technique have unraveled many principles of chromosomal folding that were subsequently linked to disease and gene regulation. In particular, Hi-C revealed that chromosomes of animals are organized into topologically associating domains (TADs), evolutionary conserved compact chromatin domains that influence gene expression. Mechanisms that underlie partitioning of the genome into TADs remain poorly understood. To explore principles of TAD folding in Drosophila melanogaster, we performed Hi-C and poly(A)(+) RNA-seq in four cell lines of various origins (S2, Kc167, DmBG3-c2, and OSC). Contrary to previous studies, we find that regions between TADs (i.e., the inter-TADs and TAD boundaries) in Drosophila are only weakly enriched with the insulator protein dCTCF, while another insulator protein Su(Hw) is preferentially present within TADs. However, Drosophila inter-TADs harbor active chromatin and constitutively transcribed (housekeeping) genes. Accordingly, we find that binding of insulator proteins dCTCF and Su(Hw) predicts TAD boundaries much worse than active chromatin marks do. Interestingly, inter-TADs correspond to decompacted inter-bands of polytene chromosomes, whereas TADs mostly correspond to densely packed bands. Collectively, our results suggest that TADs are condensed chromatin domains depleted in active chromatin marks, separated by regions of active chromatin. We propose the mechanism of TAD self-assembly based on the ability of nucleosomes from inactive chromatin to aggregate, and lack of this ability in acetylated nucleosomal arrays. Finally, we test this hypothesis by polymer simulations and find that TAD partitioning may be explained by different modes of inter-nucleosomal interactions for active and inactive chromatin.

Evolution of the Genome 3D Organization: Comparison of Fused and Segregated Globin Gene Clusters.

The genomes are folded in a complex three-dimensional (3D) structure. Some features of this organization are common for all eukaryotes, but little is known about its evolution. Here, we have studied the 3D organization and regulation of zebrafish globin gene domain and compared its organization and regulation with those of other vertebrate species. In birds and mammals, the α- and ß-globin genes are segregated into separate clusters located on different chromosomes and organized into chromatin domains of different types, whereas in cold-blooded vertebrates, including Danio rerio, α- and ß-globin genes are organized into common clusters. The major globin gene locus of Danio rerio is of particular interest as it is located in a genomic area that is syntenic in vertebrates and is controlled by a conserved enhancer. We have found that the major globin gene locus of Danio rerio is structurally and functionally segregated into two spatially distinct subloci harboring either adult or embryo-larval globin genes. These subloci demonstrate different organization at the level of chromatin domains and different modes of spatial organization, which appears to be due to selective interaction of the upstream enhancer with the sublocus harboring globin genes of the adult type. These data are discussed in terms of evolution of linear and 3D organization of gene clusters in vertebrates.

Systematic assessment of ISWI subunits shows that NURF creates local accessibility for CTCF.

Catalytic activity of the imitation switch (ISWI) family of remodelers is critical for nucleosomal organization and DNA binding of certain transcription factors, including the insulator protein CTCF. Here we define the contribution of individual subcomplexes by deriving a panel of isogenic mouse stem cell lines, each lacking one of six ISWI accessory subunits. Individual deletions of subunits of either CERF, RSF, ACF, WICH or NoRC subcomplexes only moderately affect the chromatin landscape, while removal of the NURF-specific subunit BPTF leads to a strong reduction in chromatin accessibility and SNF2H ATPase localization around CTCF sites. This affects adjacent nucleosome occupancy and CTCF binding. At a group of sites with reduced chromatin accessibility, CTCF binding persists but cohesin occupancy is reduced, resulting in decreased insulation. These results suggest that CTCF binding can be separated from its function as an insulator in nuclear organization and identify a specific role for NURF in mediating SNF2H localization and chromatin opening at bound CTCF sites.

Pairtools: from sequencing data to chromosome contacts.

The field of 3D genome organization produces large amounts of sequencing data from Hi-C and a rapidly-expanding set of other chromosome conformation protocols (3C+). Massive and heterogeneous 3C+ data require high-performance and flexible processing of sequenced reads into contact pairs. To meet these challenges, we present pairtools - a flexible suite of tools for contact extraction from sequencing data. Pairtools provides modular command-line interface (CLI) tools that can be flexibly chained into data processing pipelines. Pairtools provides both crucial core tools as well as auxiliary tools for building feature-rich 3C+ pipelines, including contact pair manipulation, filtration, and quality control. Benchmarking pairtools against popular 3C+ data pipelines shows advantages of pairtools for high-performance and flexible 3C+ analysis. Finally, pairtools provides protocol-specific tools for multi-way contacts, haplotype-resolved contacts, and single-cell Hi-C. The combination of CLI tools and tight integration with Python data analysis libraries makes pairtools a versatile foundation for a broad range of 3C+ pipelines.

Cohesin is required for long-range enhancer action at the Shh locus.

The regulatory landscapes of developmental genes in mammals can be complex, with enhancers spread over many hundreds of kilobases. It has been suggested that three-dimensional genome organization, particularly topologically associating domains formed by cohesin-mediated loop extrusion, is important for enhancers to act over such large genomic distances. By coupling acute protein degradation with synthetic activation by targeted transcription factor recruitment, here we show that cohesin, but not CTCF, is required for activation of the target gene Shh by distant enhancers in mouse embryonic stem cells. Cohesin is not required for activation directly at the promoter or by an enhancer located closer to the Shh gene. Our findings support the hypothesis that chromatin compaction via cohesin-mediated loop extrusion allows for genes to be activated by enhancers that are located many hundreds of kilobases away in the linear genome and suggests that cohesin is dispensable for enhancers located more proximally.