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.

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.

Rules of engagement for condensins and cohesins guide mitotic chromosome formation.

During mitosis, interphase chromatin is rapidly converted into rod-shaped mitotic chromosomes. Using Hi-C, imaging, proteomics and polymer modeling, we determine how the activity and interplay between loop-extruding SMC motors accomplishes this dramatic transition. Our work reveals rules of engagement for SMC complexes that are critical for allowing cells to refold interphase chromatin into mitotic chromosomes. We find that condensin disassembles interphase chromatin loop organization by evicting or displacing extrusive cohesin. In contrast, condensin bypasses cohesive cohesins, thereby maintaining sister chromatid cohesion while separating the sisters. Studies of mitotic chromosomes formed by cohesin, condensin II and condensin I alone or in combination allow us to develop new models of mitotic chromosome conformation. In these models, loops are consecutive and not overlapping, implying that condensins do not freely pass one another but stall upon encountering each other. The dynamics of Hi-C interactions and chromosome morphology reveal that during prophase loops are extruded in vivo at ~1-3 kb/sec by condensins as they form a disordered discontinuous helical scaffold within individual chromatids.

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.

Loop extrusion rules: the next generation.

The interphase genome of vertebrates contains roughly 100 000 dynamic loops formed by cohesins. These loops are thought to play important roles in many functions, but their exact contribution in each case remains hotly disputed. The key challenge in studying these loops is the lack of a single experimental technique that could reliably and comprehensively visualize their locations and dynamics. Yet, we can infer them using theoretical models that integrate complementary experimental observations. Modeling proved instrumental in showing that cohesins form loops via extrusion. The loop extrusion model made numerous successful qualitative and quantitative predictions and inspired many experiments. However, it also demonstrated limited accuracy in predicting contact maps. Recent research suggests that the original model did not fully account for the intricate details of the mechanism of loop extrusion and its complex regulation. Here, we review the progress in visualizing extrusion and characterizing the cohesin cofactors. These discoveries can be summarized as 'rules' of cohesin movement along chromosomes and incorporated into the next generation of models. Such improved models will enable more accurate inferences of positions and dynamics of cohesin loops and generate better predictions for designing experiments.

Igh and Igk loci use different folding principles for V gene recombination due to distinct chromosomal architectures of pro-B and pre-B cells.

Extended loop extrusion across the immunoglobulin heavy-chain (Igh) locus facilitates VH-DJH recombination following downregulation of the cohesin-release factor Wapl by Pax5, resulting in global changes in the chromosomal architecture of pro-B cells. Here, we demonstrate that chromatin looping and VK-JK recombination at the Igk locus were insensitive to Wapl upregulation in pre-B cells. Notably, the Wapl protein was expressed at a 2.2-fold higher level in pre-B cells compared with pro-B cells, which resulted in a distinct chromosomal architecture with normal loop sizes in pre-B cells. High-resolution chromosomal contact analysis of the Igk locus identified multiple internal loops, which likely juxtapose VK and JK elements to facilitate VK-JK recombination. The higher Wapl expression in Igµ-transgenic pre-B cells prevented extended loop extrusion at the Igh locus, leading to recombination of only the 6 most 3' proximal VH genes and likely to allelic exclusion of all other VH genes in pre-B cells. These results suggest that pro-B and pre-B cells with their distinct chromosomal architectures use different chromatin folding principles for V gene recombination, thereby enabling allelic exclusion at the Igh locus, when the Igk locus is recombined.

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.

A de novo transcription-dependent TAD boundary underpins critical multiway interactions during antibody class switch recombination.

Interactions between transcription and cohesin-mediated loop extrusion can influence 3D chromatin architecture. However, their relevance in biology is unclear. Here, we report a direct role for such interactions in the mechanism of antibody class switch recombination (CSR) at the murine immunoglobulin heavy chain locus (Igh). Using Tri-C to measure higher-order multiway interactions on single alleles, we find that the juxtaposition (synapsis) of transcriptionally active donor and acceptor Igh switch (S) sequences, an essential step in CSR, occurs via the interaction of loop extrusion complexes with a de novo topologically associating domain (TAD) boundary formed via transcriptional activity across S regions. Surprisingly, synapsis occurs predominantly in proximity to the 3' CTCF-binding element (3'CBE) rather than the Igh super-enhancer, suggesting a two-step mechanism whereby transcription of S regions is not topologically coupled to synapsis, as has been previously proposed. Altogether, these insights advance our understanding of how 3D chromatin architecture regulates CSR.

Conformation of sister chromatids in the replicated human genome.

The three-dimensional organization of the genome supports regulated gene expression, recombination, DNA repair, and chromosome segregation during mitosis. Chromosome conformation capture (Hi-C)1,2 analysis has revealed a complex genomic landscape of internal chromosomal structures in vertebrate cells3-7, but the identical sequence of sister chromatids has made it difficult to determine how they topologically interact in replicated chromosomes. Here we describe sister-chromatid-sensitive Hi-C (scsHi-C), which is based on labelling of nascent DNA with 4-thio-thymidine and nucleoside conversion chemistry. Genome-wide conformation maps of human chromosomes reveal that sister-chromatid pairs interact most frequently at the boundaries of topologically associating domains (TADs). Continuous loading of a dynamic cohesin pool separates sister-chromatid pairs inside TADs and is required to focus sister-chromatid contacts at TAD boundaries. We identified a subset of TADs that are overall highly paired and are characterized by facultative heterochromatin and insulated topological domains that form separately within individual sister chromatids. The rich pattern of sister-chromatid topologies and our scsHi-C technology will make it possible to investigate how physical interactions between identical DNA molecules contribute to DNA repair, gene expression, chromosome segregation, and potentially other biological processes.

Condensin: The Little Motor Protein that Could.

Kong et al. (2020) present the low-resolution structure of the ATPÉ£S-bound human condensin I and II complexes and demonstrate that human condensins can extrude DNA loops in a symmetric and asymmetric fashion and compact nucleosome-bound DNA.

The genome-wide multi-layered architecture of chromosome pairing in early Drosophila embryos.

Genome organization involves cis and trans chromosomal interactions, both implicated in gene regulation, development, and disease. Here, we focus on trans interactions in Drosophila, where homologous chromosomes are paired in somatic cells from embryogenesis through adulthood. We first address long-standing questions regarding the structure of embryonic homolog pairing and, to this end, develop a haplotype-resolved Hi-C approach to minimize homolog misassignment and thus robustly distinguish trans-homolog from cis contacts. This computational approach, which we call Ohm, reveals pairing to be surprisingly structured genome-wide, with trans-homolog domains, compartments, and interaction peaks, many coinciding with analogous cis features. We also find a significant genome-wide correlation between pairing, transcription during zygotic genome activation, and binding of the pioneer factor Zelda. Our findings reveal a complex, highly structured organization underlying homolog pairing, first discovered a century ago in Drosophila. Finally, we demonstrate the versatility of our haplotype-resolved approach by applying it to mammalian embryos.

Highly structured homolog pairing reflects functional organization of the Drosophila genome.

Trans-homolog interactions have been studied extensively in Drosophila, where homologs are paired in somatic cells and transvection is prevalent. Nevertheless, the detailed structure of pairing and its functional impact have not been thoroughly investigated. Accordingly, we generated a diploid cell line from divergent parents and applied haplotype-resolved Hi-C, showing that homologs pair with varying precision genome-wide, in addition to establishing trans-homolog domains and compartments. We also elucidate the structure of pairing with unprecedented detail, observing significant variation across the genome and revealing at least two forms of pairing: tight pairing, spanning contiguous small domains, and loose pairing, consisting of single larger domains. Strikingly, active genomic regions (A-type compartments, active chromatin, expressed genes) correlated with tight pairing, suggesting that pairing has a functional implication genome-wide. Finally, using RNAi and haplotype-resolved Hi-C, we show that disruption of pairing-promoting factors results in global changes in pairing, including the disruption of some interaction peaks.

A pathway for mitotic chromosome formation.

Mitotic chromosomes fold as compact arrays of chromatin loops. To identify the pathway of mitotic chromosome formation, we combined imaging and Hi-C analysis of synchronous DT40 cell cultures with polymer simulations. Here we show that in prophase, the interphase organization is rapidly lost in a condensin-dependent manner, and arrays of consecutive 60-kilobase (kb) loops are formed. During prometaphase, ~80-kb inner loops are nested within ~400-kb outer loops. The loop array acquires a helical arrangement with consecutive loops emanating from a central "spiral staircase" condensin scaffold. The size of helical turns progressively increases to ~12 megabases during prometaphase. Acute depletion of condensin I or II shows that nested loops form by differential action of the two condensins, whereas condensin II is required for helical winding.

Two independent modes of chromatin organization revealed by cohesin removal.

Imaging and chromosome conformation capture studies have revealed several layers of chromosome organization, including segregation into megabase-sized active and inactive compartments, and partitioning into sub-megabase domains (TADs). It remains unclear, however, how these layers of organization form, interact with one another and influence genome function. Here we show that deletion of the cohesin-loading factor Nipbl in mouse liver leads to a marked reorganization of chromosomal folding. TADs and associated Hi-C peaks vanish globally, even in the absence of transcriptional changes. By contrast, compartmental segregation is preserved and even reinforced. Strikingly, the disappearance of TADs unmasks a finer compartment structure that accurately reflects the underlying epigenetic landscape. These observations demonstrate that the three-dimensional organization of the genome results from the interplay of two independent mechanisms: cohesin-independent segregation of the genome into fine-scale compartments, defined by chromatin state; and cohesin-dependent formation of TADs, possibly by loop extrusion, which helps to guide distant enhancers to their target genes.

SMC complexes differentially compact mitotic chromosomes according to genomic context.

Structural maintenance of chromosomes (SMC) protein complexes are key determinants of chromosome conformation. Using Hi-C and polymer modelling, we study how cohesin and condensin, two deeply conserved SMC complexes, organize chromosomes in the budding yeast Saccharomyces cerevisiae. The canonical role of cohesin is to co-align sister chromatids, while condensin generally compacts mitotic chromosomes. We find strikingly different roles for the two complexes in budding yeast mitosis. First, cohesin is responsible for compacting mitotic chromosome arms, independently of sister chromatid cohesion. Polymer simulations demonstrate that this role can be fully accounted for through cis-looping of chromatin. Second, condensin is generally dispensable for compaction along chromosome arms. Instead, it plays a targeted role compacting the rDNA proximal regions and promoting resolution of peri-centromeric regions. Our results argue that the conserved mechanism of SMC complexes is to form chromatin loops and that distinct SMC-dependent looping activities are selectively deployed to appropriately compact chromosomes.

Targeted Degradation of CTCF Decouples Local Insulation of Chromosome Domains from Genomic Compartmentalization.

The molecular mechanisms underlying folding of mammalian chromosomes remain poorly understood. The transcription factor CTCF is a candidate regulator of chromosomal structure. Using the auxin-inducible degron system in mouse embryonic stem cells, we show that CTCF is absolutely and dose-dependently required for looping between CTCF target sites and insulation of topologically associating domains (TADs). Restoring CTCF reinstates proper architecture on altered chromosomes, indicating a powerful instructive function for CTCF in chromatin folding. CTCF remains essential for TAD organization in non-dividing cells. Surprisingly, active and inactive genome compartments remain properly segregated upon CTCF depletion, revealing that compartmentalization of mammalian chromosomes emerges independently of proper insulation of TADs. Furthermore, our data support that CTCF mediates transcriptional insulator function through enhancer blocking but not as a direct barrier to heterochromatin spreading. Beyond defining the functions of CTCF in chromosome folding, these results provide new fundamental insights into the rules governing mammalian genome organization.

Emerging Evidence of Chromosome Folding by Loop Extrusion.

Chromosome organization poses a remarkable physical problem with many biological consequences: How can molecular interactions between proteins at the nanometer scale organize micron-long chromatinized DNA molecules, insulating or facilitating interactions between specific genomic elements? The mechanism of active loop extrusion holds great promise for explaining interphase and mitotic chromosome folding, yet remains difficult to assay directly. We discuss predictions from our polymer models of loop extrusion with barrier elements and review recent experimental studies that provide strong support for loop extrusion, focusing on perturbations to CTCF and cohesin assayed via Hi-C in interphase. Finally, we discuss a likely molecular mechanism of loop extrusion by structural maintenance of chromosomes complexes.

Micro-C XL: assaying chromosome conformation from the nucleosome to the entire genome.

We present Micro-C XL, an improved method for analysis of chromosome folding at mononucleosome resolution. Using long crosslinkers and isolation of insoluble chromatin, Micro-C XL increases signal-to-noise ratio. Micro-C XL maps of budding and fission yeast genomes capture both short-range chromosome fiber features such as chromosomally interacting domains and higher order features such as centromere clustering. Micro-C XL provides a single assay to interrogate chromosome folding at length scales from the nucleosome to the full genome.

Hematopoietic Stem Cells Are the Major Source of Multilineage Hematopoiesis in Adult Animals.

Hematopoietic stem cells (HSCs) sustain long-term reconstitution of hematopoiesis in transplantation recipients, yet their role in the endogenous steady-state hematopoiesis remains unclear. In particular, recent studies suggested that HSCs provide a relatively minor contribution to immune cell development in adults. We directed transgene expression in a fraction of HSCs that maintained reconstituting activity during serial transplantations. Inducible genetic labeling showed that transgene-expressing HSCs gave rise to other phenotypic HSCs, confirming their top position in the differentiation hierarchy. The labeled HSCs rapidly contributed to committed progenitors of all lineages and to mature myeloid cells and lymphocytes, but not to B-1a cells or tissue macrophages. Importantly, labeled HSCs gave rise to more than two-thirds of all myeloid cells and platelets in adult mice, and this contribution could be accelerated by an induced interferon response. Thus, classically defined HSCs maintain immune cell development in the steady state and during systemic cytokine responses.

Chromosome Compaction by Active Loop Extrusion.

During cell division, chromosomes are compacted in length by more than a 100-fold. A wide range of experiments demonstrated that in their compacted state, mammalian chromosomes form arrays of closely stacked consecutive ∼100 kb loops. The mechanism underlying the active process of chromosome compaction into a stack of loops is unknown. Here we test the hypothesis that chromosomes are compacted by enzymatic machines that actively extrude chromatin loops. When such loop-extruding factors (LEF) bind to chromosomes, they progressively bridge sites that are further away along the chromosome, thus extruding a loop. We demonstrate that collective action of LEFs leads to formation of a dynamic array of consecutive loops. Simulations and an analytically solved model identify two distinct steady states: a sparse state, where loops are highly dynamic but provide little compaction; and a dense state, where there are more stable loops and dramatic chromosome compaction. We find that human chromosomes operate at the border of the dense steady state. Our analysis also shows how the macroscopic characteristics of the loop array are determined by the microscopic properties of LEFs and their abundance. When the number of LEFs are used that match experimentally based estimates, the model can quantitatively reproduce the average loop length, the degree of compaction, and the general loop-array morphology of compact human chromosomes. Our study demonstrates that efficient chromosome compaction can be achieved solely by an active loop-extrusion process.