Detecting chromatin interactions between and along sister chromatids with SisterC

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.

Cohesin-Dependent Loop Extrusion: Molecular Mechanics and Role in Cell Physiology

The most prominent representatives of multisubunit SMC complexes, cohesin and condensin, are best known as structural components of mitotic chromosomes. It turned out that these complexes, as well as their bacterial homologues, are molecular motors, the ATP-dependent movement of these complexes along DNA threads leads to the formation of DNA loops. In recent years, we have witnessed an avalanche-like accumulation of data on the process of SMC dependent DNA looping, also known as loop extrusion. This review briefly summarizes the current understanding of the place and role of cohesin-dependent extrusion in cell physiology and presents a number of models describing the potential molecular mechanism of extrusion in a most compelling way. We conclude the review with a discussion of how the capacity of cohesin to extrude DNA loops may be mechanistically linked to its involvement in sister chromatid cohesion.

Acute multi-level response to defective de novo chromatin assembly in S-phase

Long-term perturbation of de novo chromatin assembly during DNA replication has profound effects on epigenome maintenance and cell fate. The early mechanistic origin of these defects is unknown. Here, we combine acute degradation of Chromatin Assembly Factor 1 (CAF-1), a key player in de novo chromatin assembly, with single-cell genomics, quantitative proteomics, and live-microscopy to uncover these initiating mechanisms in human cells. CAF-1 loss immediately slows down DNA replication speed and renders nascent DNA hyperaccessible. A rapid cellular response, distinct from canonical DNA damage signaling, is triggered and lowers histone mRNAs. As a result, histone variants usage and their modifications are altered, limiting transcriptional fidelity and delaying chromatin maturation within a single S-phase. This multi-level response induces a cell-cycle arrest after mitosis. Our work reveals the immediate consequences of defective de novo chromatin assembly during DNA replication, explaining how at later times the epigenome and cell fate can be altered.

Chromosome structure modeling tools and their evaluation in bacteria

The three-dimensional (3D) structure of bacterial chromosomes is crucial for understanding chromosome function. With the growing availability of high-throughput chromosome conformation capture (3C/Hi-C) data, the 3D structure reconstruction algorithms have become powerful tools to study bacterial chromosome structure and function. It is highly desired to have a recommendation on the chromosome structure reconstruction tools to facilitate the prokaryotic 3D genomics. In this work, we review existing chromosome 3D structure reconstruction algorithms and classify them based on their underlying computational models into two categories: constraint-based modeling and thermodynamics-based modeling. We briefly compare these algorithms utilizing 3C/Hi-C datasets and fluorescence microscopy data obtained from Escherichia coli and Caulobacter crescentus, as well as simulated datasets. We discuss current challenges in the 3D reconstruction algorithms for bacterial chromosomes, primarily focusing on software usability. Finally, we briefly prospect future research directions for bacterial chromosome structure reconstruction algorithms.

Replicating Chromosomes in Whole-Cell Models of Bacteria

Computational models of cells cannot be considered complete unless they include the most fundamental process of life, the replication of genetic material. In a recent study, we presented a computational framework to model systems of replicating bacterial chromosomes as polymers at 10 bp resolution with Brownian dynamics. This approach was used to investigate changes in chromosome organization during replication and extend the applicability of an existing whole-cell model (WCM) for a genetically minimal bacterium, JCVI-syn3A, to the entire cell cycle. To achieve cell-scale chromosome structures that are realistic, we modeled the chromosome as a self-avoiding homopolymer with bending and torsional stiffnesses that capture the essential mechanical properties of dsDNA in Syn3A. Additionally, the polymer interacts with ribosomes distributed according to cryo-electron tomograms of Syn3A. The polymer model was further augmented by computational models of loop extrusion by structural maintenance of chromosomes (SMC) protein complexes and topoisomerase action, and the modeling and analysis of multi-fork replication states.

3D models of fungal chromosomes to enhance visual integration of omics data

The functions of eukaryotic chromosomes and their spatial architecture in the nucleus are reciprocally dependent. Hi-C experiments are routinely used to study chromosome 3D organization by probing chromatin interactions. Standard representation of the data has relied on contact maps that show the frequency of interactions between parts of the genome. In parallel, it has become easier to build 3D models of the entire genome based on the same Hi-C data, and thus benefit from the methodology and visualization tools developed for structural biology. 3D modeling of entire genomes leverages the understanding of their spatial organization. However, this opportunity for original and insightful modeling is underexploited. In this paper, we show how seeing the spatial organization of chromosomes can bring new perspectives to omics data integration. We assembled state-of-the-art tools into a workflow that goes from Hi-C raw data to fully annotated 3D models and we re-analysed public omics datasets available for three fungal species. Besides the well-described properties of the spatial organization of their chromosomes (Rabl conformation, hypercoiling and chromosome territories), our results highlighted (i) in Saccharomyces cerevisiae, the backbones of the cohesin anchor regions, which were aligned all along the chromosomes, (ii) in Schizosaccharomyces pombe, the oscillations of the coiling of chromosome arms throughout the cell cycle and (iii) in Neurospora crassa, the massive relocalization of histone marks in mutants of heterochromatin regulators. 3D modeling of the chromosomes brings new opportunities for visual integration of omics data. This holistic perspective supports intuition and lays the foundation for building new concepts.

Meiosis: Dances Between Homologs

The raison d'être of meiosis is shuffling of genetic information via Mendelian segregation and, within individual chromosomes, by DNA crossing-over. These outcomes are enabled by a complex cellular program in which interactions between homologous chromosomes play a central role. We first provide a background regarding the basic principles of this program. We then summarize the current understanding of the DNA events of recombination and of three processes that involve whole chromosomes: homolog pairing, crossover interference, and chiasma maturation. All of these processes are implemented by direct physical interaction of recombination complexes with underlying chromosome structures. Finally, we present convergent lines of evidence that the meiotic program may have evolved by coupling of this interaction to late-stage mitotic chromosome morphogenesis.

Dynamics of chromosome organization in a minimal bacterial cell

Computational models of cells cannot be considered complete unless they include the most fundamental process of life, the replication and inheritance of genetic material. By creating a computational framework to model systems of replicating bacterial chromosomes as polymers at 10 bp resolution with Brownian dynamics, we investigate changes in chromosome organization during replication and extend the applicability of an existing whole-cell model (WCM) for a genetically minimal bacterium, JCVI-syn3A, to the entire cell-cycle. To achieve cell-scale chromosome structures that are realistic, we model the chromosome as a self-avoiding homopolymer with bending and torsional stiffnesses that capture the essential mechanical properties of dsDNA in Syn3A. In addition, the conformations of the circular DNA must avoid overlapping with ribosomes identitied in cryo-electron tomograms. While Syn3A lacks the complex regulatory systems known to orchestrate chromosome segregation in other bacteria, its minimized genome retains essential loop-extruding structural maintenance of chromosomes (SMC) protein complexes (SMC-scpAB) and topoisomerases. Through implementing the effects of these proteins in our simulations of replicating chromosomes, we find that they alone are sufficient for simultaneous chromosome segregation across all generations within nested theta structures. This supports previous studies suggesting loop-extrusion serves as a near-universal mechanism for chromosome organization within bacterial and eukaryotic cells. Furthermore, we analyze ribosome diffusion under the influence of the chromosome and calculate in silico chromosome contact maps that capture inter-daughter interactions. Finally, we present a methodology to map the polymer model of the chromosome to a Martini coarse-grained representation to prepare molecular dynamics models of entire Syn3A cells, which serves as an ultimate means of validation for cell states predicted by the WCM.

A spatial genome aligner for resolving chromatin architectures from multiplexed DNA FISH

Multiplexed fluorescence in situ hybridization (FISH) is a widely used approach for analyzing three-dimensional genome organization, but it is challenging to derive chromosomal conformations from noisy fluorescence signals, and tracing chromatin is not straightforward. Here we report a spatial genome aligner that parses true chromatin signal from noise by aligning signals to a DNA polymer model. Using genomic distances separating imaged loci, our aligner estimates spatial distances expected to separate loci on a polymer in three-dimensional space. Our aligner then evaluates the physical probability observed signals belonging to these loci are connected, thereby tracing chromatin structures. We demonstrate that this spatial genome aligner can efficiently model chromosome architectures from DNA FISH data across multiple scales and be used to predict chromosome ploidies de novo in interphase cells. Reprocessing of previous whole-genome chromosome tracing data with this method indicates the spatial aggregation of sister chromatids in S/G2 phase cells in asynchronous mouse embryonic stem cells and provides evidence for extranumerary chromosomes that remain tightly paired in postmitotic neurons of the adult mouse cortex.

Transcription-Driven Translocation of Cohesive and Non-Cohesive Cohesin In Vivo

Cohesin is a central architectural element of chromosomes that regulates numerous DNA-based events. The complex holds sister chromatids together until anaphase onset and organizes individual chromosomal DNAs into loops and self-associating domains. Purified cohesin diffuses along DNA in an ATP-independent manner but can be propelled by transcribing RNA polymerase. In conjunction with a cofactor, the complex also extrudes DNA loops in an ATP-dependent manner. In this study we examine transcription-driven translocation of cohesin under various conditions in yeast. To this end, obstacles of increasing size were tethered to DNA to act as roadblocks to complexes mobilized by an inducible gene. The obstacles were built from a GFP-lacI core fused to one or more mCherries. A chimera with four mCherries blocked cohesin passage in late G1. During M phase, the threshold barrier depended on the state of cohesion: non-cohesive complexes were also blocked by four mCherries whereas cohesive complexes were blocked by as few as three mCherries. Furthermore cohesive complexes that were stalled at obstacles, in turn, blocked the passage of non-cohesive complexes. That synthetic barriers capture mobilized cohesin demonstrates that transcription-driven complexes translocate processively in vivo. Together, this study reveals unexplored limitations to cohesin movement on chromosomes.

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.

Direct observation of surface charge and stiffness of human metaphase chromosomes

Metaphase chromosomes in which both polynucleotides and proteins are condensed with hierarchies are closely related to life phenomena such as cell division, cancer development, and cellular senescence. Nevertheless, their nature is rarely revealed, owing to their structural complexity and technical limitations in analytical methods. In this study, we used surface potential and nanomechanics mapping technology based on atomic force microscopy to measure the surface charge and intrinsic stiffness of metaphase chromosomes. We found that extra materials covering the chromosomes after the extraction process were positively charged. With the covering materials, the chromosomes were positively charged (ca. 44.9 ± 16.48 mV) and showed uniform stiffness (ca. 6.23 ± 1.98 MPa). In contrast, after getting rid of the extra materials through treatment with RNase and protease, the chromosomes were strongly negatively charged (ca. -197.4 ± 77.87 mV) and showed relatively non-uniform and augmented stiffness (ca. 36.87 ± 17.56 MPa). The results suggested undulating but compact coordination of condensed chromosomes. Additionally, excessive treatment with RNase and protease could destroy the chromosomal structure, providing an exceptional opportunity for multiscale stiffness mapping of polynucleotides, nucleosomes, chromatin fibers, and chromosomes in a single image. Our approach offers a new horizon in terms of an analytical technique for studying chromosome-related diseases.

Hypothesis-driven probabilistic modelling enables a principled perspective of genomic compartments

The Hi-C method has revolutionized the study of genome organization, yet interpretation of Hi-C interaction frequency maps remains a major challenge. Genomic compartments are a checkered Hi-C interaction pattern suggested to represent the partitioning of the genome into two self-interacting states associated with active and inactive chromatin. Based on a few elementary mechanistic assumptions, we derive a generative probabilistic model of genomic compartments, called deGeco. Testing our model, we find it can explain observed Hi-C interaction maps in a highly robust manner, allowing accurate inference of interaction probability maps from extremely sparse data without any training of parameters. Taking advantage of the interpretability of the model parameters, we then test hypotheses regarding the nature of genomic compartments. We find clear evidence of multiple states, and that these states self-interact with different affinities. We also find that the interaction rules of chromatin states differ considerably within and between chromosomes. Inspecting the molecular underpinnings of a four-state model, we show that a simple classifier can use histone marks to predict the underlying states with 87% accuracy. Finally, we observe instances of mixed-state loci and analyze these loci in single-cell Hi-C maps, finding that mixing of states occurs mainly at the cell level.

Regulation of the mitotic chromosome folding machines

Over the last several years enormous progress has been made in identifying the molecular machines, including condensins and topoisomerases that fold mitotic chromosomes. The discovery that condensins generate chromatin loops through loop extrusion has revolutionized, and energized, the field of chromosome folding. To understand how these machines fold chromosomes with the appropriate dimensions, while disentangling sister chromatids, it needs to be determined how they are regulated and deployed. Here, we outline the current understanding of how these machines and factors are regulated through cell cycle dependent expression, chromatin localization, activation and inactivation through post-translational modifications, and through associations with each other, with other factors and with the chromatin template itself. There are still many open questions about how condensins and topoisomerases are regulated but given the pace of progress in the chromosome folding field, it seems likely that many of these will be answered in the years ahead.

The 3D genome landscape: Diverse chromosomal interactions and their functional implications

Genome organization includes contacts both within a single chromosome and between distinct chromosomes. Thus, regulatory organization in the nucleus may include interplay of these two types of chromosomal interactions with genome activity. Emerging advances in omics and single-cell imaging technologies have allowed new insights into chromosomal contacts, including those of homologs and sister chromatids, and their significance to genome function. In this review, we highlight recent studies in this field and discuss their impact on understanding the principles of chromosome organization and associated functional implications in diverse cellular processes. Specifically, we describe the contributions of intra-chromosomal, inter-homolog, and inter-sister chromatid contacts to genome organization and gene expression.

Chromosome organization through the cell cycle at a glance

Genome organization and the three-dimensional folding of chromosomes are now seen as major contributors to nearly all nuclear functions including gene regulation, replication and repair. Recent studies have shown that in addition to the dramatic metamorphoses in chromosome conformation associated with entry to, and exit from mitosis, chromosomes undergo continual conformational changes throughout interphase with differential dynamics in loop structure, topological domains, compartments and lamina-associated domains. Understanding and accounting for these cell-cycle-dependent conformational changes is essential for the interpretation of data from a growing array of powerful molecular techniques to investigate genome conformation function, and to identify the molecules and mechanisms that drive chromosome conformational changes. In this Cell Science at a Glance article and the accompanying poster, we review Hi-C and microscopy studies describing cell-cycle-dependent conformational changes in chromosome structure.

Sister chromatid-sensitive Hi-C to map the conformation of replicated genomes

Chromosome conformation capture (Hi-C) techniques map the 3D organization of entire genomes. How sister chromatids fold in replicated chromosomes, however, cannot be determined with conventional Hi-C because of the identical DNA sequences of sister chromatids. Here, we present a protocol for sister chromatid-sensitive Hi-C (scsHi-C) that enables the distinction of DNA contacts within individual sister chromatids (cis sister contacts) from those between sister chromatids (trans sister contacts), thereby allowing investigation of the organization of replicated genomes. scsHi-C is based on live-cell labeling of nascent DNA by the synthetic nucleoside 4-thio-thymidine (4sT), which incorporates into a distinct DNA strand on each sister chromatid because of semi-conservative DNA replication. After purification of genomic DNA and in situ Hi-C library preparation, 4sT is chemically converted into 5-methyl-cytosine in the presence of OsO₄/NH₄Cl to introduce T-to-C signature point mutations on 4sT-labeled DNA. The Hi-C library is then sequenced, and ligated fragments are assigned to sister chromatids on the basis of strand orientation and the presence of signature mutations. The ensemble of scsHi-C contacts thereby represents genome-wide contact probabilities within and across sister chromatids. scsHi-C can be completed in 2 weeks, has been successfully applied in HeLa cells and can potentially be established for any cell type that allows proper cell cycle synchronization and incorporation of sufficient amounts of 4sT. The genome-wide maps of replicated chromosomes detected by scsHi-C enable investigation of the molecular mechanisms shaping sister chromatid topologies and the relevance of sister chromatid conformation in crucial processes like DNA repair, mitotic chromosome formation and potentially other biological processes.

Perspectives for the reconstruction of 3D chromatin conformation using single cell Hi-C data

Construction of chromosomes 3D models based on single cell Hi-C data constitute an important challenge. We present a reconstruction approach, DPDchrom, that incorporates basic knowledge whether the reconstructed conformation should be coil-like or globular and spring relaxation at contact sites. In contrast to previously published protocols, DPDchrom can naturally form globular conformation due to the presence of explicit solvent. Benchmarking of this and several other methods on artificial polymer models reveals similar reconstruction accuracy at high contact density and DPDchrom advantage at low contact density. To compare 3D structures insensitively to spatial orientation and scale, we propose the Modified Jaccard Index. We analyzed two sources of the contact dropout: contact radius change and random contact sampling. We found that the reconstruction accuracy exponentially depends on the number of contacts per genomic bin allowing to estimate the reconstruction accuracy in advance. We applied DPDchrom to model chromosome configurations based on single-cell Hi-C data of mouse oocytes and found that these configurations differ significantly from a random one, that is consistent with other studies.

Conformational statistics of non-equilibrium polymer loops in Rouse model with active loop extrusion

Motivated by the recent experimental observations of the DNA loop extrusion by protein motors, in this paper, we investigate the statistical properties of the growing polymer loops within the ideal chain model. The loop conformation is characterized statistically by the mean gyration radius and the pairwise contact probabilities. It turns out that a single dimensionless parameter, which is given by the ratio of the loop relaxation time over the time elapsed since the start of extrusion, controls the crossover between near-equilibrium and highly non-equilibrium asymptotics in the statistics of the extruded loop, regardless of the specific time dependence of the extrusion velocity. In addition, we show that two-sided and one-sided loop extruding motors produce the loops with almost identical properties. Our predictions are based on two rigorous semi-analytical methods accompanied by asymptotic analysis of slow and fast extrusion limits.

Understanding 3D genome organization by multidisciplinary methods

Understanding how chromatin is folded in the nucleus is fundamental to understanding its function. Although 3D genome organization has been historically difficult to study owing to a lack of relevant methodologies, major technological breakthroughs in genome-wide mapping of chromatin contacts and advances in imaging technologies in the twenty-first century considerably improved our understanding of chromosome conformation and nuclear architecture. In this Review, we discuss methods of 3D genome organization analysis, including sequencing-based techniques, such as Hi-C and its derivatives, Micro-C, DamID and others; microscopy-based techniques, such as super-resolution imaging coupled with fluorescence in situ hybridization (FISH), multiplex FISH, in situ genome sequencing and live microscopy methods; and computational and modelling approaches. We describe the most commonly used techniques and their contribution to our current knowledge of nuclear architecture and, finally, we provide a perspective on up-and-coming methods that open possibilities for future major discoveries.

Mitotic chromosomes

Our understanding of the structure and function of mitotic chromosomes has come a long way since these iconic objects were first recognized more than 140 years ago, though many details remain to be elucidated. In this chapter, we start with the early history of chromosome studies and then describe the path that led to our current understanding of the formation and structure of mitotic chromosomes. We also discuss some of the remaining questions. It is now well established that each mitotic chromatid consists of a central organizing region containing a so-called "chromosome scaffold" from which loops of DNA project radially. Only a few key non-histone proteins and protein complexes are required to form the chromosome: topoisomerase IIα, cohesin, condensin I and condensin II, and the chromokinesin KIF4A. These proteins are concentrated along the axis of the chromatid. Condensins I and II are primarily responsible for shaping the chromosome and the scaffold, and they produce the loops of DNA by an ATP-dependent process known as loop extrusion. Modelling of Hi-C data suggests that condensin II adopts a spiral staircase arrangement with an extruded loop extending out from each step in a roughly helical pattern. Condensin I then forms loops nested within these larger condensin II loops, thereby giving rise to the final compaction of the mitotic chromosome in a process that requires Topo IIα.

The Spatial Organization of Sister Chromatids

Understanding how genomes are spatially organized is central to many aspects of cell biology. However, it has been difficult to study the relationships between sister chromatids because sequencing-based techniques such as Hi-C could not distinguish identical sister DNAs. Here, I discuss recent developments that provide insights into sister chromatid organization.