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.

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 V_H-DJ_H 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 V_K-J_K 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 V_K and J_K elements to facilitate V_K-J_K 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 V_H genes and likely to allelic exclusion of all other V_H 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.

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.

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.

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.

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.

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.

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.

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.

Anxiety-depressive disorders in children: Neurobiological and neurohormonal aspects

Objective

The need for diagnosis and correct classification of depression among children is dictated by its burdeness with age, high risk of recurrence at further stages of child development and propensity to suicidal behavior.

Materials and methods

One hundred and sixty adolescents with ADD were included in our study. The study design comprised: clinicopsychopathological, somatoneurological, psychological and neurohormonal methods.

Results

Clustering symptom of anxiety and depression have shown that for children in early puberty is typical more formation of somatic (35.2%), behavioral (21.6%), phobic (21.6%) variants of depression, less-anxiety (13.5%), asthenia (8.1%). For children in puberty – apathetic (30.1%), anxiety (28.9%), dismorfofobic (27.7%), behavioral (13.3%) variants. The proportion of suicidal behavior of depression increases in proportion to age, mainly due to suicidal thoughts, sayings, auto-aggressive behavior. Analysis of the formation conditions of anxiety and depression in children showed a significant correlation of genetic, biological and socio-environmental components.

Symptoms of the minimal brain dysfunction (MBD) at an early ontogeny (prognostic value = +4.8), loaded natal period (PV = +4.2), frequent colds in the medical history (PV = 3.7), signs of cerebro-organic failure (PV = +3.8) and obesity as an endocrine disorder (PV = +2.1), sex and age of the manifestation of the first depressive episode in the early stages ontogenesis in boys aged 7 years (PS = +5.3), in girls aged 9 years (PS = +3.9) have been registered among biological risk factors for the ADD formation. Prognostic significance of neurohormonal parameters as regards the risk for ADD formation in adolescents is based on a decreased serotonin level in patients with depressive anxiety (PV = +2.8) and a reduced melatonin level, irrespective of the variant of clinical depression (PV = +2.4).

Conclusion

Determination of neurobiological and neurohormonal risk markers for the development of the ADD makes it possible to carry out psychoprophylactic measures.

Disclosure of interest

The authors have not supplied their declaration of competing interest.

A technology for providing therapeutic training of children with the anxiety-depressive disorders, comorbid with obesity

Background and aims

The study considers the problem of the anxiety-depressive disorders (ADD) in children with concomitant obesity in the context of the search for effective methods of their correction.

Materials and methods

Included 64 children at puberty with the anxiety-depressive disorders, comorbid with obesity. Design comprised: clinicopsychopathological, somatoneurological, psychological monitoring (CDRS-R); the Spielberger, “Me and my illness” and “Man in the Rain” projective drawing tests.

Results

The symptoms of emotional and vegetative tension, depression registered in all studied children.

In the developed model of support, the interventions of primary level are aimed on family and closest encirclement of the child (family psychotherapy, psycho-educational programs for children and their parents). Strategies:

– stressful situation in the family (practical, psychological support and education of parents);

– for relief of depressive and vegetative-anxious symptoms (using art therapy, cognitive behavioral therapy);

– for education of patients in the wider context – quality of life, coping, motivation and modification of inappropriate behaviors (modeling of new cognitions).

The technology includes the diagnosis of the psychosomatic core of the disease, areas of neurotic fixation of the child with ADD and obesity, correction of emotional homeostasis and cognitive imbalance by activation of the personality individual resources as a result using cognitive-behavioral and art-dynamic therapy, as well as a “Control of eating behavior” training program, and solving certain situational problems by a child with the disease.

Conclusion

Our method makes the tactics of psychotherapeutic and psycho-educational intervention at an early stage of the disease.

Disclosure of interest

The authors have not supplied their declaration of competing interest.

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.

Formation of Chromosomal Domains by Loop Extrusion

Topologically associating domains (TADs) are fundamental structural and functional building blocks of human interphase chromosomes, yet the mechanisms of TAD formation remain unclear. Here, we propose that loop extrusion underlies TAD formation. In this process, cis-acting loop-extruding factors, likely cohesins, form progressively larger loops but stall at TAD boundaries due to interactions with boundary proteins, including CTCF. Using polymer simulations, we show that this model produces TADs and finer-scale features of Hi-C data. Each TAD emerges from multiple loops dynamically formed through extrusion, contrary to typical illustrations of single static loops. Loop extrusion both explains diverse experimental observations-including the preferential orientation of CTCF motifs, enrichments of architectural proteins at TAD boundaries, and boundary deletion experiments-and makes specific predictions for the depletion of CTCF versus cohesin. Finally, loop extrusion has potentially far-ranging consequences for processes such as enhancer-promoter interactions, orientation-specific chromosomal looping, and compaction of mitotic chromosomes.

Compaction and segregation of sister chromatids via active loop extrusion

The mechanism by which chromatids and chromosomes are segregated during mitosis and meiosis is a major puzzle of biology and biophysics. Using polymer simulations of chromosome dynamics, we show that a single mechanism of loop extrusion by condensins can robustly compact, segregate and disentangle chromosomes, arriving at individualized chromatids with morphology observed in vivo. Our model resolves the paradox of topological simplification concomitant with chromosome 'condensation', and explains how enzymes a few nanometers in size are able to control chromosome geometry and topology at micron length scales. We suggest that loop extrusion is a universal mechanism of genome folding that mediates functional interactions during interphase and compacts chromosomes during mitosis.

DOI:http://dx.doi.org/10.7554/eLife.14864.001