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.

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.

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.

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.

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.

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.

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.

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.

Nuclear organisation and replication timing are coupled through RIF1-PP1 interaction.

Three-dimensional genome organisation and replication timing are known to be correlated, however, it remains unknown whether nuclear architecture overall plays an instructive role in the replication-timing programme and, if so, how. Here we demonstrate that RIF1 is a molecular hub that co-regulates both processes. Both nuclear organisation and replication timing depend upon the interaction between RIF1 and PP1. However, whereas nuclear architecture requires the full complement of RIF1 and its interaction with PP1, replication timing is not sensitive to RIF1 dosage. The role of RIF1 in replication timing also extends beyond its interaction with PP1. Availing of this separation-of-function approach, we have therefore identified in RIF1 dual function the molecular bases of the co-dependency of the replication-timing programme and nuclear architecture.

From bedside to bench: regulation of host factors in SARS-CoV-2 infection.

The zoonotic coronavirus SARS-CoV-2 (severe acute respiratory syndrome coronavirus-2), which causes COVID-19 (coronavirus disease-2019), has resulted in a pandemic. This has led to an urgent need to understand the molecular determinants of SARS-CoV-2 infection, factors associated with COVID-19 heterogeneity and severity, and therapeutic options for these patients. In this review, we discuss the role of host factors in SARS-CoV-2 infection and describe variations in host factor expression as mechanisms underlying the symptoms and severity of COVID-19. We focus on two host factors, angiotensin-converting enzyme 2 (ACE2) and transmembrane serine protease 2 (TMPRSS2), implicated in SARS-CoV-2 infection. We also discuss genetic variants associated with COVID-19 severity revealed in selected patients and based on genome-wide association studies (GWASs). Furthermore, we highlight important advances in cell and chromatin biology, such as single-cell RNA and chromatin sequencing and chromosomal conformation assays, as methods that may aid in the discovery of viral-host interactions in COVID-19. Understanding how regulation of host factor genes varies in physiological and pathological states might explain the heterogeneity observed in SARS-CoV-2 infection, help identify pathways for therapeutic development, and identify patients most likely to progress to severe COVID-19.

Order and stochasticity in the folding of individual Drosophila genomes.

Mammalian and Drosophila genomes are partitioned into topologically associating domains (TADs). Although this partitioning has been reported to be functionally relevant, it is unclear whether TADs represent true physical units located at the same genomic positions in each cell nucleus or emerge as an average of numerous alternative chromatin folding patterns in a cell population. Here, we use a single-nucleus Hi-C technique to construct high-resolution Hi-C maps in individual Drosophila genomes. These maps demonstrate chromatin compartmentalization at the megabase scale and partitioning of the genome into non-hierarchical TADs at the scale of 100 kb, which closely resembles the TAD profile in the bulk in situ Hi-C data. Over 40% of TAD boundaries are conserved between individual nuclei and possess a high level of active epigenetic marks. Polymer simulations demonstrate that chromatin folding is best described by the random walk model within TADs and is most suitably approximated by a crumpled globule build of Gaussian blobs at longer distances. We observe prominent cell-to-cell variability in the long-range contacts between either active genome loci or between Polycomb-bound regions, suggesting an important contribution of stochastic processes to the formation of the Drosophila 3D 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.

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.

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.

DNA Methylation Directs Polycomb-Dependent 3D Genome Re-organization in Naive Pluripotency.

The DNA hypomethylation that occurs when embryonic stem cells (ESCs) are directed to the ground state of naive pluripotency by culturing in two small molecule inhibitors (2i) results in redistribution of polycomb (H3K27me3) away from its target loci. Here, we demonstrate that 3D genome organization is also altered in 2i, with chromatin decompaction at polycomb target loci and a loss of long-range polycomb interactions. By preventing DNA hypomethylation during the transition to the ground state, we are able to restore to ESC in 2i the H3K27me3 distribution, as well as polycomb-mediated 3D genome organization that is characteristic of primed ESCs grown in serum. However, these cells retain the functional characteristics of 2i ground-state ESCs. Our findings demonstrate the central role of DNA methylation in shaping major aspects of 3D genome organization but caution against assuming causal roles for the epigenome and 3D genome in gene regulation and function in ESCs.

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.

Nuclear lamina integrity is required for proper spatial organization of chromatin in Drosophila.

How the nuclear lamina (NL) impacts on global chromatin architecture is poorly understood. Here, we show that NL disruption in Drosophila S2 cells leads to chromatin compaction and repositioning from the nuclear envelope. This increases the chromatin density in a fraction of topologically-associating domains (TADs) enriched in active chromatin and enhances interactions between active and inactive chromatin. Importantly, upon NL disruption the NL-associated TADs become more acetylated at histone H3 and less compact, while background transcription is derepressed. Two-colour FISH confirms that a TAD becomes less compact following its release from the NL. Finally, polymer simulations show that chromatin binding to the NL can per se compact attached TADs. Collectively, our findings demonstrate a dual function of the NL in shaping the 3D genome. Attachment of TADs to the NL makes them more condensed but decreases the overall chromatin density in the nucleus by stretching interphase chromosomes.