Computer vision for pattern detection in chromosome contact maps

scGrapHiC: deep learning-based graph deconvolution for Hi-C using single cell gene expression

SUMMARY

Single-cell Hi-C (scHi-C) protocol helps identify cell-type-specific chromatin interactions and sheds light on cell differentiation and disease progression. Despite providing crucial insights, scHi-C data is often underutilized due to the high cost and the complexity of the experimental protocol. We present a deep learning framework, scGrapHiC, that predicts pseudo-bulk scHi-C contact maps using pseudo-bulk scRNA-seq data. Specifically, scGrapHiC performs graph deconvolution to extract genome-wide single-cell interactions from a bulk Hi-C contact map using scRNA-seq as a guiding signal. Our evaluations show that scGrapHiC, trained on seven cell-type co-assay datasets, outperforms typical sequence encoder approaches. For example, scGrapHiC achieves a substantial improvement of 23.2% in recovering cell-type-specific Topologically Associating Domains over the baselines. It also generalizes to unseen embryo and brain tissue samples. scGrapHiC is a novel method to generate cell-type-specific scHi-C contact maps using widely available genomic signals that enables the study of cell-type-specific chromatin interactions.

AVAILABILITY AND IMPLEMENTATION

The GitHub link: https://github.com/rsinghlab/scGrapHiC contains the source code of scGrapHiC and associated scripts to preprocess publicly available datasets to produce the results and visualizations we have discuss in this manuscript.

Enhancing Hi-C contact matrices for loop detection with Capricorn: a multiview diffusion model

MOTIVATION

High-resolution Hi-C contact matrices reveal the detailed three-dimensional architecture of the genome, but high-coverage experimental Hi-C data are expensive to generate. Simultaneously, chromatin structure analyses struggle with extremely sparse contact matrices. To address this problem, computational methods to enhance low-coverage contact matrices have been developed, but existing methods are largely based on resolution enhancement methods for natural images and hence often employ models that do not distinguish between biologically meaningful contacts, such as loops and other stochastic contacts.

RESULTS

We present Capricorn, a machine learning model for Hi-C resolution enhancement that incorporates small-scale chromatin features as additional views of the input Hi-C contact matrix and leverages a diffusion probability model backbone to generate a high-coverage matrix. We show that Capricorn outperforms the state of the art in a cross-cell-line setting, improving on existing methods by 17% in mean squared error and 26% in F1 score for chromatin loop identification from the generated high-coverage data. We also demonstrate that Capricorn performs well in the cross-chromosome setting and cross-chromosome, cross-cell-line setting, improving the downstream loop F1 score by 14% relative to existing methods. We further show that our multiview idea can also be used to improve several existing methods, HiCARN and HiCNN, indicating the wide applicability of this approach. Finally, we use DNA sequence to validate discovered loops and find that the fraction of CTCF-supported loops from Capricorn is similar to those identified from the high-coverage data. Capricorn is a powerful Hi-C resolution enhancement method that enables scientists to find chromatin features that cannot be identified in the low-coverage contact matrix.

AVAILABILITY AND IMPLEMENTATION

Implementation of Capricorn and source code for reproducing all figures in this paper are available at https://github.com/CHNFTQ/Capricorn.

Inhibiting AGTR1 reduces AML burden and protects the heart from cardiotoxicity in mouse models

Clinical treatment of acute myeloid leukemia (AML) largely relies on intensive chemotherapy. However, the application of chemotherapy is often hindered by cardiotoxicity. Patient sequence data revealed that angiotensin II receptor type 1 (AGTR1) is a shared target between AML and cardiovascular disease (CVD). We found that inhibiting AGTR1 sensitized AML to chemotherapy and protected the heart against chemotherapy-induced cardiotoxicity in a human AML cell-transplanted mouse model. These effects were regulated by the AGTR1-Notch1 axis in AML cells and cardiomyocytes from mice. In mouse cardiomyocytes, AGTR1 was hyperactivated by AML and chemotherapy. AML leukemogenesis increased the expression of the angiotensin-converting enzyme and led to increased production of angiotensin II, the ligand of AGTR1, in an MLL-AF9-driven AML mouse model. In this model, the AGTR1-Notch1 axis regulated a variety of genes involved with cell stemness and chemotherapy resistance. AML cell stemness was reduced after Agtr1a deletion in the mouse AML cell transplant model. Mechanistically, Agtr1a deletion decreased γ-secretase formation, which is required for transmembrane Notch1 cleavage and release of the Notch1 intracellular domain into the nucleus. Using multiomics, we identified AGTR1-Notch1 signaling downstream genes and found decreased binding between these gene sequences with Notch1 and chromatin enhancers, as well as increased binding with silencers. These findings describe an AML/CVD association that may be used to improve AML treatment.

Light control of three-dimensional chromatin organization in soybean

Higher-order chromatin structure is critical for regulation of gene expression. In plants, light profoundly affects the morphogenesis of emerging seedlings as well as global gene expression to ensure optimal adaptation to environmental conditions. However, the changes and functional significance of chromatin organization in response to light during seedling development are not well documented. We constructed Hi-C contact maps for the cotyledon, apical hook and hypocotyl of soybean subjected to dark and light conditions. The resulting high-resolution Hi-C contact maps identified chromosome territories, A/B compartments, A/B sub-compartments, TADs (Topologically Associated Domains) and chromatin loops in each organ. We observed increased chromatin compaction under light and we found that domains that switched from B sub-compartments in darkness to A sub-compartments under light contained genes that were activated during photomorphogenesis. At the local scale, we identified a group of TADs constructed by gene clusters consisting of different numbers of Small Auxin-Upregulated RNAs (SAURs), which exhibited strict co-expression in the hook and hypocotyl in response to light stimulation. In the hypocotyl, RNA polymerase II (RNAPII) regulated the transcription of a SAURs cluster under light via TAD condensation. Our results suggest that the 3D genome is involved in the regulation of light-related gene expression in a tissue-specific manner.

A continuum of zinc finger transcription factor retention on native chromatin underlies dynamic genome organization

Transcription factor (TF) residence on chromatin translates into quantitative transcriptional or structural outcomes on genome. Commonly used formaldehyde crosslinking fixes TF-DNA interactions cumulatively and compromises the measured occupancy level. Here we mapped the occupancy level of global or individual zinc finger TFs like CTCF and MAZ, in the form of highly resolved footprints, on native chromatin. By incorporating reinforcing perturbation conditions, we established S-score, a quantitative metric to proxy the continuum of CTCF or MAZ retention across different motifs on native chromatin. The native chromatin-retained CTCF sites harbor sequence features within CTCF motifs better explained by S-score than the metrics obtained from other crosslinking or native assays. CTCF retention on native chromatin correlates with local SUMOylation level, and anti-correlates with transcriptional activity. The S-score successfully delineates the otherwise-masked differential stability of chromatin structures mediated by CTCF, or by MAZ independent of CTCF. Overall, our study established a paradigm continuum of TF retention across binding sites on native chromatin, explaining the dynamic genome organization.

Genome assembly of the southern pine beetle (Dendroctonus frontalis Zimmerman) reveals the origins of gene content reduction in Dendroctonus

Dendroctonus frontalis, also known as southern pine beetle (SPB), represents the most damaging forest pest in the southeastern United States. Strategies to predict, monitor and suppress SPB outbreaks have had limited success. Genomic data are critical to inform on pest biology and to identify molecular targets to develop improved management approaches. Here, we produced a chromosome-level genome assembly of SPB using long-read sequencing data. Synteny analyses confirmed the conservation of the core coleopteran Stevens elements and validated the bona fide SPB X chromosome. Transcriptomic data were used to obtain 39,588 transcripts corresponding to 13,354 putative protein-coding loci. Comparative analyses of gene content across 14 beetle and 3 other insects revealed several losses of conserved genes in the Dendroctonus clade and gene gains in SPB and Dendroctonus that were enriched for loci encoding membrane proteins and extracellular matrix proteins. While lineage-specific gene losses contributed to the gene content reduction observed in Dendroctonus, we also showed that widespread misannotation of transposable elements represents a major cause of the apparent gene expansion in several non-Dendroctonus species. Our findings uncovered distinctive features of the SPB gene complement and disentangled the role of biological and annotation-related factors contributing to gene content variation across beetles.

CD-Loop: a chromatin loop detection method based on the diffusion model

Motivation

In recent years, there have been significant advances in various chromatin conformation capture techniques, and annotating the topological structure from Hi-C contact maps has become crucial for studying the three-dimensional structure of chromosomes. However, the structure and function of chromatin loops are highly dynamic and diverse, influenced by multiple factors. Therefore, obtaining the three-dimensional structure of the genome remains a challenging task. Among many chromatin loop prediction methods, it is difficult to fully extract features from the contact map and make accurate predictions at low sequencing depths.

Results

In this study, we put forward a deep learning framework based on the diffusion model called CD-Loop for predicting accurate chromatin loops. First, by pre-training the input data, we obtain prior probabilities for predicting the classification of the Hi-C contact map. Then, by combining the denoising process based on the diffusion model and the prior probability obtained by pre-training, candidate loops were predicted from the input Hi-C contact map. Finally, CD-Loop uses a density-based clustering algorithm to cluster the candidate chromatin loops and predict the final chromatin loops. We compared CD-Loop with the currently popular methods, such as Peakachu, Chromosight, and Mustache, and found that in different cell types, species, and sequencing depths, CD-Loop outperforms other methods in loop annotation. We conclude that CD-Loop can accurately predict chromatin loops and reveal cell-type specificity. The code is available at https://github.com/wangyang199897/CD-Loop.

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.

High resolution maps of chromatin reorganization through mouse meiosis reveal novel features of the 3D meiotic structure

When germ cells transition from the mitotic cycle into meiotic prophase I (MPI), chromosomes condense into an array of chromatin loops that are required to promote homolog pairing and genetic recombination. To identify the changes in chromosomal conformation, we isolated nuclei on a trajectory from spermatogonia to the end of MPI. At each stage along this trajectory, we built genomic interaction maps with the highest temporal and spatial resolution to date. The changes in chromatin folding coincided with a concurrent decline in mitotic cohesion and a rise in meiotic cohesin complexes. We found that the stereotypical large-scale A and B compartmentalization was lost during meiotic prophase I alongside the loss of topological associating domains (TADs). Still, local subcompartments were detected and maintained throughout meiosis. The enhanced Micro-C resolution revealed that, despite the loss of TADs, higher frequency contact sites between two loci were detectable during meiotic prophase I coinciding with CTCF bound sites. The pattern of interactions around these CTCF sites with their neighboring loci showed that CTCF sites were often anchoring the meiotic loops. Additionally, the localization of CTCF to the meiotic axes indicated that these anchors were at the base of loops. Strikingly, even in the face of the dramatic reconfiguration of interphase chromatin into a condensed loop-array, the interactions between regulatory elements remained well preserved. This establishes a potential mechanism for how the meiotic chromatin maintains active transcription within a highly structured genome. In summary, the high temporal and spatial resolution of these data revealed previously unappreciated aspects of mammalian meiotic chromatin organization.

Sister chromatid cohesion halts DNA loop expansion

Eukaryotic genomes are folded into DNA loops mediated by structural maintenance of chromosomes (SMC) complexes such as cohesin, condensin, and Smc5/6. This organization regulates different DNA-related processes along the cell cycle, such as transcription, recombination, segregation, and DNA repair. During the G2 stage, SMC-mediated DNA loops coexist with cohesin complexes involved in sister chromatid cohesion (SCC). However, the articulation between the establishment of SCC and the formation of SMC-mediated DNA loops along the chromatin remains unknown. Here, we show that SCC is indeed a barrier to cohesin-mediated DNA loop expansion along G2/M Saccharomyces cerevisiae chromosomes.

Comparative study on chromatin loop callers using Hi-C data reveals their effectiveness

BACKGROUND

Chromosome is one of the most fundamental part of cell biology where DNA holds the hierarchical information. DNA compacts its size by forming loops, and these regions house various protein particles, including CTCF, SMC3, H3 histone. Numerous sequencing methods, such as Hi-C, ChIP-seq, and Micro-C, have been developed to investigate these properties. Utilizing these data, scientists have developed a variety of loop prediction techniques that have greatly improved their methods for characterizing loop prediction and related aspects.

RESULTS

In this study, we categorized 22 loop calling methods and conducted a comprehensive study of 11 of them. Additionally, we have provided detailed insights into the methodologies underlying these algorithms for loop detection, categorizing them into five distinct groups based on their fundamental approaches. Furthermore, we have included critical information such as resolution, input and output formats, and parameters. For this analysis, we utilized the GM12878 Hi-C datasets at 5 KB, 10 KB, 100 KB and 250 KB resolutions. Our evaluation criteria encompassed various factors, including memory usages, running time, sequencing depth, and recovery of protein-specific sites such as CTCF, H3K27ac, and RNAPII.

CONCLUSION

This analysis offers insights into the loop detection processes of each method, along with the strengths and weaknesses of each, enabling readers to effectively choose suitable methods for their datasets. We evaluate the capabilities of these tools and introduce a novel Biological, Consistency, and Computational robustness score ( B C C score ) to measure their overall robustness ensuring a comprehensive evaluation of their performance.

Absence of chromosome axis protein recruitment prevents meiotic recombination chromosome-wide in the budding yeast Lachancea kluyveri

Meiotic recombination shows broad variations across species and along chromosomes and is often suppressed at and around genomic regions determining sexual compatibility such as mating type loci in fungi. Here, we show that the absence of Spo11-DSBs and meiotic recombination on Lakl0C-left, the chromosome arm containing the sex locus of the Lachancea kluyveri budding yeast, results from the absence of recruitment of the two chromosome axis proteins Red1 and Hop1, essential for proper Spo11-DSBs formation. Furthermore, cytological observation of spread pachytene meiotic chromosomes reveals that Lakl0C-left does not undergo synapsis. However, we show that the behavior of Lakl0C-left is independent of its particularly early replication timing and is not accompanied by any peculiar chromosome structure as detectable by Hi-C in this yet poorly studied yeast. Finally, we observed an accumulation of heterozygous mutations on Lakl0C-left and a sexual dimorphism of the haploid meiotic offspring, supporting a direct effect of this absence of meiotic recombination on L. kluyveri genome evolution and fitness. Because suppression of meiotic recombination on sex chromosomes is widely observed across eukaryotes, the mechanism for recombination suppression described here may apply to other species, with the potential to impact sex chromosome evolution.

Transcription-induced domains form the elementary constraining building blocks of bacterial chromosomes

Transcription generates local topological and mechanical constraints on the DNA fiber, leading to the generation of supercoiled chromosome domains in bacteria. However, the global impact of transcription on chromosome organization remains elusive, as the scale of genes and operons in bacteria remains well below the resolution of chromosomal contact maps generated using Hi-C (~5-10 kb). Here we combined sub-kb Hi-C contact maps and chromosome engineering to visualize individual transcriptional units. We show that transcriptional units form discrete three-dimensional transcription-induced domains that impose mechanical and topological constraints on their neighboring sequences at larger scales, modifying their localization and dynamics. These results show that transcriptional domains constitute primary building blocks of bacterial chromosome folding and locally impose structural and dynamic constraints.

Orchestrating chromosome conformation capture analysis with Bioconductor

Genome-wide chromatin conformation capture assays provide formidable insights into the spatial organization of genomes. However, due to the complexity of the data structure, their integration in multi-omics workflows remains challenging. We present data structures, computational methods and visualization tools available in Bioconductor to investigate Hi-C, micro-C and other 3C-related data, in R. An online book ( https://bioconductor.org/books/OHCA/ ) further provides prospective end users with a number of workflows to process, import, analyze and visualize any type of chromosome conformation capture data.

ZNF143 deletion alters enhancer/promoter looping and CTCF/cohesin geometry

The transcription factor ZNF143 contains a central domain of seven zinc fingers in a tandem array and is involved in 3D genome construction. However, the mechanism by which ZNF143 functions in chromatin looping remains unclear. Here, we show that ZNF143 directionally recognizes a diverse range of genomic sites directly within enhancers and promoters and is required for chromatin looping between these sites. In addition, ZNF143 is located between CTCF and cohesin at numerous CTCF sites, and ZNF143 removal narrows the space between CTCF and cohesin. Moreover, genetic deletion of ZNF143, in conjunction with acute CTCF degradation, reveals that ZNF143 and CTCF collaborate to regulate higher-order topological chromatin organization. Finally, CTCF depletion enlarges direct ZNF143 chromatin looping. Thus, ZNF143 is recruited by CTCF to the CTCF sites to regulate CTCF/cohesin configuration and TAD (topologically associating domain) formation, whereas directional recognition of genomic DNA motifs directly by ZNF143 itself regulates promoter activity via chromatin looping.

Phage-inducible chromosomal minimalist islands (PICMIs), a novel family of small marine satellites of virulent phages

Phage satellites are bacterial genetic elements that co-opt phage machinery for their own dissemination. Here we identify a family of satellites, named Phage-Inducible Chromosomal Minimalist Islands (PICMIs), that are broadly distributed in marine bacteria of the family Vibrionaceae. A typical PICMI is characterized by reduced gene content, does not encode genes for capsid remodelling, and packages its DNA as a concatemer. PICMIs integrate in the bacterial host genome next to the fis regulator, and encode three core proteins necessary for excision and replication. PICMIs are dependent on virulent phage particles to spread to other bacteria, and protect their hosts from other competitive phages without interfering with their helper phage. Thus, our work broadens our understanding of phage satellites and narrows down the minimal number of functions necessary to hijack a tailed phage.

Chromosome architecture in an archaeal species naturally lacking structural maintenance of chromosomes proteins

Proteins in the structural maintenance of chromosomes (SMC) superfamily play key roles in chromosome organization and are ubiquitous across all domains of life. However, SMC proteins are notably absent in the Desulfurococcales of phylum Crenarchaeota. Intrigued by this observation, we performed chromosome conformation capture experiments in the model Desulfurococcales species Aeropyrum pernix. As in other archaea, we observe chromosomal interaction domains across the chromosome. The boundaries between chromosomal interaction domains show a dependence on transcription and translation for their definition. Importantly, however, we reveal an additional higher-order, bipartite organization of the chromosome-with a small high-gene-expression and self-interacting domain that is defined by transcriptional activity and loop structures. Viewing these data in the context of the distribution of SMC superfamily proteins in the Crenarchaeota, we suggest that the organization of the Aeropyrum genome represents an evolutionary antecedent of the compartmentalized architecture observed in the Sulfolobus lineage.

A Comprehensive Evaluation of Generalizability of Deep Learning-Based Hi-C Resolution Improvement Methods

Hi-C is a widely used technique to study the 3D organization of the genome. Due to its high sequencing cost, most of the generated datasets are of a coarse resolution, which makes it impractical to study finer chromatin features such as Topologically Associating Domains (TADs) and chromatin loops. Multiple deep learning-based methods have recently been proposed to increase the resolution of these datasets by imputing Hi-C reads (typically called upscaling). However, the existing works evaluate these methods on either synthetically downsampled datasets, or a small subset of experimentally generated sparse Hi-C datasets, making it hard to establish their generalizability in the real-world use case. We present our framework-Hi-CY-that compares existing Hi-C resolution upscaling methods on seven experimentally generated low-resolution Hi-C datasets belonging to various levels of read sparsities originating from three cell lines on a comprehensive set of evaluation metrics. Hi-CY also includes four downstream analysis tasks, such as TAD and chromatin loops recall, to provide a thorough report on the generalizability of these methods. We observe that existing deep learning methods fail to generalize to experimentally generated sparse Hi-C datasets, showing a performance reduction of up to 57%. As a potential solution, we find that retraining deep learning-based methods with experimentally generated Hi-C datasets improves performance by up to 31%. More importantly, Hi-CY shows that even with retraining, the existing deep learning-based methods struggle to recover biological features such as chromatin loops and TADs when provided with sparse Hi-C datasets. Our study, through the Hi-CY framework, highlights the need for rigorous evaluation in the future. We identify specific avenues for improvements in the current deep learning-based Hi-C upscaling methods, including but not limited to using experimentally generated datasets for training.

H3K27me3 spreading organizes canonical PRC1 chromatin architecture to regulate developmental programs

Polycomb Repressive Complex 2 (PRC2)-mediated histone H3K27 tri-methylation (H3K27me3) recruits canonical PRC1 (cPRC1) to maintain heterochromatin. In early development, polycomb-regulated genes are connected through long-range 3D interactions which resolve upon differentiation. Here, we report that polycomb looping is controlled by H3K27me3 spreading and regulates target gene silencing and cell fate specification. Using glioma-derived H3 Lys-27-Met (H3K27M) mutations as tools to restrict H3K27me3 deposition, we show that H3K27me3 confinement concentrates the chromatin pool of cPRC1, resulting in heightened 3D interactions mirroring chromatin architecture of pluripotency, and stringent gene repression that maintains cells in progenitor states to facilitate tumor development. Conversely, H3K27me3 spread in pluripotent stem cells, following neural differentiation or loss of the H3K36 methyltransferase NSD1, dilutes cPRC1 concentration and dissolves polycomb loops. These results identify the regulatory principles and disease implications of polycomb looping and nominate histone modification-guided distribution of reader complexes as an important mechanism for nuclear compartment organization.

Highlights

The confinement of H3K27me3 at PRC2 nucleation sites without its spreading correlates with increased 3D chromatin interactions.The H3K27M oncohistone concentrates canonical PRC1 that anchors chromatin loop interactions in gliomas, silencing developmental programs.Stem and progenitor cells require factors promoting H3K27me3 confinement, including H3K36me2, to maintain cPRC1 loop architecture.The cPRC1-H3K27me3 interaction is a targetable driver of aberrant self-renewal in tumor cells.

3D chromatin interactions involving Drosophila insulators are infrequent but preferential and arise before TADs and transcription

In mammals, insulators contribute to the regulation of loop extrusion to organize chromatin into topologically associating domains. In Drosophila the role of insulators in 3D genome organization is, however, under current debate. Here, we addressed this question by combining bioinformatics analysis and multiplexed chromatin imaging. We describe a class of Drosophila insulators enriched at regions forming preferential chromatin interactions genome-wide. Notably, most of these 3D interactions do not involve TAD borders. Multiplexed imaging shows that these interactions occur infrequently, and only rarely involve multiple genomic regions coalescing together in space in single cells. Finally, we show that non-border preferential 3D interactions enriched in this class of insulators are present before TADs and transcription during Drosophila development. Our results are inconsistent with insulators forming stable hubs in single cells, and instead suggest that they fine-tune existing 3D chromatin interactions, providing an additional regulatory layer for transcriptional regulation.

The shape of chromatin: insights from computational recognition of geometric patterns in Hi-C data

The three-dimensional organization of chromatin plays a crucial role in gene regulation and cellular processes like deoxyribonucleic acid (DNA) transcription, replication and repair. Hi-C and related techniques provide detailed views of spatial proximities within the nucleus. However, data analysis is challenging partially due to a lack of well-defined, underpinning mathematical frameworks. Recently, recognizing and analyzing geometric patterns in Hi-C data has emerged as a powerful approach. This review provides a summary of algorithms for automatic recognition and analysis of geometric patterns in Hi-C data and their correspondence with chromatin structure. We classify existing algorithms on the basis of the data representation and pattern recognition paradigm they make use of. Finally, we outline some of the challenges ahead and promising future directions.

Conserved H3K27me3-associated chromatin looping mediates physical interactions of gene clusters in plants

Higher-order chromatin organization is essential for transcriptional regulation, genome stability maintenance, and other genome functions. Increasing evidence has revealed significant differences in 3D chromatin organization between plants and animals. However, the extent, pattern, and rules of chromatin organization in plants are still unclear. In this study, we systematically identified and characterized long-range chromatin loops in the Arabidopsis 3D genome. We identified hundreds of long-range cis chromatin loops and found their anchor regions are closely associated with H3K27me3 epigenetic modifications. Furthermore, we demonstrated that these chromatin loops are dependent on Polycomb group (PcG) proteins, suggesting that the Polycomb repressive complex 2 (PRC2) complex is essential for establishing and maintaining these novel loops. Although most of these PcG-medicated chromatin loops are stable, many of these loops are tissue-specific or dynamically regulated by different treatments. Interestingly, tandemly arrayed gene clusters and metabolic gene clusters are enriched in anchor regions. Long-range H3K27me3-marked chromatin interactions are associated with the coregulation of specific gene clusters. Finally, we also identified H3K27me3-associated chromatin loops associated with gene clusters in Oryza sativa and Glycine max, indicating that these long-range chromatin loops are conserved in plants. Our results provide novel insights into genome evolution and transcriptional coregulation in plants.

Extrusion fountains are hallmarks of chromosome organization emerging upon zygotic genome activation

The first activation of gene expression during development (zygotic genome activation, ZGA) is accompanied by massive changes in chromosome organization. The connection between these two processes remains unknown. Using Hi-C for zebrafish embryos, we found that chromosome folding starts by establishing "fountains", novel elements of chromosome organization, emerging selectively at enhancers upon ZGA. Using polymer simulations, we demonstrate that fountains can emerge as sites of targeted cohesin loading and require two-sided, yet desynchronized, loop extrusion. Specific loss of fountains upon loss of pioneer transcription factors that drive ZGA reveals a causal connection between enhancer activity and fountain formation. Finally, we show that fountains emerge in early Medaka and Xenopus embryos; moreover, we found cohesin-dependent fountain pattern on enhancers of mouse embryonic stem cells. Taken together, fountains are the first enhancer-specific elements of chromosome organization; they constitute starting points of chromosome folding during early development, likely serving as sites of targeted cohesin loading.

The genome of the reef-building glass sponge Aphrocallistes vastus provides insights into silica biomineralization

Well-annotated and contiguous genomes are an indispensable resource for understanding the evolution, development, and metabolic capacities of organisms. Sponges, an ecologically important non-bilaterian group of primarily filter-feeding sessile aquatic organisms, are underrepresented with respect to available genomic resources. Here we provide a high-quality and well-annotated genome of Aphrocallistes vastus, a glass sponge (Porifera: Hexactinellida) that forms large reef structures off the coast of British Columbia (Canada). We show that its genome is approximately 80 Mb, small compared to most other metazoans, and contains nearly 2500 nested genes, more than other genomes. Hexactinellida is characterized by a unique skeletal architecture made of amorphous silicon dioxide (SiO₂), and we identified 419 differentially expressed genes between the osculum, i.e. the vertical growth zone of the sponge, and the main body. Among the upregulated ones, mineralization-related genes such as glassin, as well as collagens and actins, dominate the expression profile during growth. Silicateins, suggested being involved in silica mineralization, especially in demosponges, were not found at all in the A. vastus genome and suggests that the underlying mechanisms of SiO₂ deposition in the Silicea sensu stricto (Hexactinellida + Demospongiae) may not be homologous.

CTCF, BEAF-32, and CP190 are not required for the establishment of TADs in early Drosophila embryos but have locus-specific roles

The boundaries of topologically associating domains (TADs) are delimited by insulators and/or active promoters; however, how they are initially established during embryogenesis remains unclear. Here, we examined this during the first hours of Drosophila embryogenesis. DNA-FISH confirms that intra-TAD pairwise proximity is established during zygotic genome activation (ZGA) but with extensive cell-to-cell heterogeneity. Most newly formed boundaries are occupied by combinations of CTCF, BEAF-32, and/or CP190. Depleting each insulator individually from chromatin revealed that TADs can still establish, although with lower insulation, with a subset of boundaries (~10%) being more dependent on specific insulators. Some weakened boundaries have aberrant gene expression due to unconstrained enhancer activity. However, the majority of misexpressed genes have no obvious direct relationship to changes in domain-boundary insulation. Deletion of an active promoter (thereby blocking transcription) at one boundary had a greater impact than deleting the insulator-bound region itself. This suggests that cross-talk between insulators and active promoters and/or transcription might reinforce domain boundary insulation during embryogenesis.

Chromosome arm length, and a species-specific determinant, define chromosome arm width

Mitotic chromosomes in different organisms adopt various dimensions. What defines these dimensions is scarcely understood. Here, we compare mitotic chromosomes in budding and fission yeasts harboring similarly sized genomes distributed among 16 or 3 chromosomes, respectively. Hi-C analyses and superresolution microscopy reveal that budding yeast chromosomes are characterized by shorter-ranging mitotic chromatin contacts and are thinner compared with the thicker fission yeast chromosomes that contain longer-ranging mitotic contacts. These distinctions persist even after budding yeast chromosomes are fused to form three fission-yeast-length entities, revealing a species-specific organizing principle. Species-specific widths correlate with the known binding site intervals of the chromosomal condensin complex. Unexpectedly, within each species, we find that longer chromosome arms are always thicker and harbor longer-ranging contacts, a trend that we also observe with human chromosomes. Arm length as a chromosome width determinant informs mitotic chromosome formation models.

Reference panel guided topological structure annotation of Hi-C data

Accurately annotating topological structures (e.g., loops and topologically associating domains) from Hi-C data is critical for understanding the role of 3D genome organization in gene regulation. This is a challenging task, especially at high resolution, in part due to the limited sequencing coverage of Hi-C data. Current approaches focus on the analysis of individual Hi-C data sets of interest, without taking advantage of the facts that (i) several hundred Hi-C contact maps are publicly available, and (ii) the vast majority of topological structures are conserved across multiple cell types. Here, we present RefHiC, an attention-based deep learning framework that uses a reference panel of Hi-C datasets to facilitate topological structure annotation from a given study sample. We compare RefHiC against tools that do not use reference samples and find that RefHiC outperforms other programs at both topological associating domain and loop annotation across different cell types, species, and sequencing depths.

Chromosome-scale assemblies of Acanthamoeba castellanii genomes provide insights into Legionella pneumophila infection-related chromatin reorganization

The unicellular amoeba Acanthamoeba castellanii is ubiquitous in aquatic environments, where it preys on bacteria. The organism also hosts bacterial endosymbionts, some of which are parasitic, including human pathogens such as Chlamydia and Legionella spp. Here we report complete, high-quality genome sequences for two extensively studied A. castellanii strains, Neff and C3. Combining long- and short-read data with Hi-C, we generated near chromosome-level assemblies for both strains with 90% of the genome contained in 29 scaffolds for the Neff strain and 31 for the C3 strain. Comparative genomics revealed strain-specific functional enrichment, most notably genes related to signal transduction in the C3 strain and to viral replication in Neff. Furthermore, we characterized the spatial organization of the A. castellanii genome and showed that it is reorganized during infection by Legionella pneumophila Infection-dependent chromatin loops were found to be enriched in genes for signal transduction and phosphorylation processes. In genomic regions where chromatin organization changed during Legionella infection, we found functional enrichment for genes associated with metabolism, organelle assembly, and cytoskeleton organization. Given Legionella infection is known to alter its host's cell cycle, to exploit the host's organelles, and to modulate the host's metabolism in its favor, these changes in chromatin organization may partly be related to mechanisms of host control during Legionella infection.

Karyotype engineering reveals spatio-temporal control of replication firing and gene contacts

Eukaryotic genomes vary in terms of size, chromosome number, and genetic complexity. Their temporal organization is complex, reflecting coordination between DNA folding and function. Here, we used fused karyotypes of budding yeast to characterize the effects of chromosome length on nuclear architecture. We found that size-matched megachromosomes expand to occupy a larger fraction of the enlarged nucleus. Hi-C maps reveal changes in the three-dimensional structure corresponding to inactivated centromeres and telomeres. De-clustering of inactive centromeres results in their loss of early replication, highlighting a functional correlation between genome organization and replication timing. Repositioning of former telomere-proximal regions on chromosome arms exposed a subset of contacts between flocculin genes. Chromatin reorganization of megachromosomes during cell division remained unperturbed, and it revealed that centromere-rDNA contacts in anaphase, extending over 0.3 Mb on wild-type chromosome, cannot exceed ∼1.7 Mb. Our results highlight the relevance of engineered karyotypes to unveiling relationships between genome organization and function.

Loop detection using Hi-C data with HiCExplorer

BACKGROUND

Chromatin loops are an essential factor in the structural organization of the genome; however, their detection in Hi-C interaction matrices is a challenging and compute-intensive task. The approach presented here, integrated into the HiCExplorer software, shows a chromatin loop detection algorithm that applies a strict candidate selection based on continuous negative binomial distributions and performs a Wilcoxon rank-sum test to detect enriched Hi-C interactions.

RESULTS

HiCExplorer's loop detection has a high detection rate and accuracy. It is the fastest available CPU implementation and utilizes all threads offered by modern multicore platforms.

CONCLUSIONS

HiCExplorer's method to detect loops by using a continuous negative binomial function combined with the donut approach from HiCCUPS leads to reliable and fast computation of loops. All the loop-calling algorithms investigated provide differing results, which intersect by $\sim 50\%$ at most. The tested in situ Hi-C data contain a large amount of noise; achieving better agreement between loop calling algorithms will require cleaner Hi-C data and therefore future improvements to the experimental methods that generate the data.

Nucleome programming is required for the foundation of totipotency in mammalian germline development

Germ cells are unique in engendering totipotency, yet the mechanisms underlying this capacity remain elusive. Here, we perform comprehensive and in-depth nucleome analysis of mouse germ-cell development in vitro, encompassing pluripotent precursors, primordial germ cells (PGCs) before and after epigenetic reprogramming, and spermatogonia/spermatogonial stem cells (SSCs). Although epigenetic reprogramming, including genome-wide DNA de-methylation, creates broadly open chromatin with abundant enhancer-like signatures, the augmented chromatin insulation safeguards transcriptional fidelity. These insulatory constraints are then erased en masse for spermatogonial development. Notably, despite distinguishing epigenetic programming, including global DNA re-methylation, the PGCs-to-spermatogonia/SSCs development entails further euchromatization. This accompanies substantial erasure of lamina-associated domains, generating spermatogonia/SSCs with a minimal peripheral attachment of chromatin except for pericentromeres-an architecture conserved in primates. Accordingly, faulty nucleome maturation, including persistent insulation and improper euchromatization, leads to impaired spermatogenic potential. Given that PGCs after epigenetic reprogramming serve as oogenic progenitors as well, our findings elucidate a principle for the nucleome programming that creates gametogenic progenitors in both sexes, defining a basis for nuclear totipotency.

Smc3 acetylation, Pds5 and Scc2 control the translocase activity that establishes cohesin-dependent chromatin loops

Cohesin is a DNA translocase that is instrumental in the folding of the genome into chromatin loops, with functional consequences on DNA-related processes. Chromatin loop length and organization likely depend on cohesin processivity, translocation rate and stability on DNA. Here, we investigate and provide a comprehensive overview of the roles of various cohesin regulators in tuning chromatin loop expansion in budding yeast Saccharomyces cerevisiae. We demonstrate that Scc2, which stimulates cohesin ATPase activity, is also essential for cohesin translocation, driving loop expansion in vivo. Smc3 acetylation during the S phase counteracts this activity through the stabilization of Pds5, which finely tunes the size and stability of loops in G2.

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.

Stripenn detects architectural stripes from chromatin conformation data using computer vision

Architectural stripes tend to form at genomic regions harboring genes with salient roles in cell identity and function. Therefore, the accurate identification and quantification of these features are essential for understanding lineage-specific gene regulation. Here, we present Stripenn, an algorithm rooted in computer vision to systematically detect and quantitate architectural stripes from chromatin conformation measurements using various technologies. We demonstrate that Stripenn outperforms existing methods and highlight its biological applications in the context of B and T lymphocytes. By comparing stripes across distinct cell types and different species, we find that these chromatin features are highly conserved and form at genes with prominent roles in cell-type-specific processes. In summary, Stripenn is a computational method that borrows concepts from widely used image processing techniques to demarcate and quantify architectural stripes.

Chromosome conformation capture techniques have recently revealed features beyond chromatin loops such as architectural stripes. Here the authors present their stripe detection tool ‘Stripenn’ to detect and quantitate stripes from any type of chromatin conformation capture data. They show that architectural stripes are enriched at transcriptionally active and accessible genomic regions.

Stripenn detects architectural stripes from chromatin conformation data using computer vision

Architectural stripes tend to form at genomic regions harboring genes with salient roles in cell identity and function. Therefore, the accurate identification and quantification of these features are essential for understanding lineage-specific gene regulation. Here, we present Stripenn, an algorithm rooted in computer vision to systematically detect and quantitate architectural stripes from chromatin conformation measurements using various technologies. We demonstrate that Stripenn outperforms existing methods and highlight its biological applications in the context of B and T lymphocytes. By comparing stripes across distinct cell types and different species, we find that these chromatin features are highly conserved and form at genes with prominent roles in cell-type-specific processes. In summary, Stripenn is a computational method that borrows concepts from widely used image processing techniques to demarcate and quantify architectural stripes.

Interferons reshape the 3D conformation and accessibility of macrophage chromatin

Engagement of macrophages in innate immune responses is directed by type I and type II interferons (IFN-I and IFN-γ, respectively). IFN triggers drastic changes in cellular transcriptomes, executed by JAK-STAT signal transduction and the transcriptional control of interferon-stimulated genes (ISG) by STAT transcription factors. Here, we study the immediate-early nuclear response to IFN-I and IFN-γ in murine macrophages. We show that the mechanism of gene control by both cytokines includes a rapid increase of DNA accessibility and rearrangement of the 3D chromatin contacts particularly between open chromatin of ISG loci. IFN-stimulated gene factor 3 (ISGF3), the major transcriptional regulator of ISG, controlled homeostatic and, most notably, induced-state DNA accessibility at a subset of ISG. Increases in DNA accessibility correlated with the appearance of activating histone marks at surrounding nucleosomes. Collectively our data emphasize changes in the three-dimensional nuclear space and epigenome as an important facet of transcriptional control by the IFN-induced JAK-STAT pathway.

Comparative analysis of the Mercenaria mercenaria genome provides insights into the diversity of transposable elements and immune molecules in bivalve mollusks

BACKGROUND

The hard clam Mercenaria mercenaria is a major marine resource along the Atlantic coasts of North America and has been introduced to other continents for resource restoration or aquaculture activities. Significant mortality events have been reported in the species throughout its native range as a result of diseases (microbial infections, leukemia) and acute environmental stress. In this context, the characterization of the hard clam genome can provide highly needed resources to enable basic (e.g., oncogenesis and cancer transmission, adaptation biology) and applied (clam stock enhancement, genomic selection) sciences.

RESULTS

Using a combination of long and short-read sequencing technologies, a 1.86 Gb chromosome-level assembly of the clam genome was generated. The assembly was scaffolded into 19 chromosomes, with an N50 of 83 Mb. Genome annotation yielded 34,728 predicted protein-coding genes, markedly more than the few other members of the Venerida sequenced so far, with coding regions representing only 2% of the assembly. Indeed, more than half of the genome is composed of repeated elements, including transposable elements. Major chromosome rearrangements were detected between this assembly and another recent assembly derived from a genetically segregated clam stock. Comparative analysis of the clam genome allowed the identification of a marked diversification in immune-related proteins, particularly extensive tandem duplications and expansions in tumor necrosis factors (TNFs) and C1q domain-containing proteins, some of which were previously shown to play a role in clam interactions with infectious microbes. The study also generated a comparative repertoire highlighting the diversity and, in some instances, the specificity of LTR-retrotransposons elements, particularly Steamer elements in bivalves.

CONCLUSIONS

The diversity of immune molecules in M. mercenaria may allow this species to cope with varying and complex microbial and environmental landscapes. The repertoire of transposable elements identified in this study, particularly Steamer elements, should be a prime target for the investigation of cancer cell development and transmission among bivalve mollusks.

Enhancer-promoter interactions and transcription are largely maintained upon acute loss of CTCF, cohesin, WAPL or YY1

It remains unclear why acute depletion of CTCF (CCCTC-binding factor) and cohesin only marginally affects expression of most genes despite substantially perturbing three-dimensional (3D) genome folding at the level of domains and structural loops. To address this conundrum, we used high-resolution Micro-C and nascent transcript profiling in mouse embryonic stem cells. We find that enhancer-promoter (E-P) interactions are largely insensitive to acute (3-h) depletion of CTCF, cohesin or WAPL. YY1 has been proposed as a structural regulator of E-P loops, but acute YY1 depletion also had minimal effects on E-P loops, transcription and 3D genome folding. Strikingly, live-cell, single-molecule imaging revealed that cohesin depletion reduced transcription factor (TF) binding to chromatin. Thus, although CTCF, cohesin, WAPL or YY1 is not required for the short-term maintenance of most E-P interactions and gene expression, our results suggest that cohesin may facilitate TFs to search for and bind their targets more efficiently.

HiC1Dmetrics: framework to extract various one-dimensional features from chromosome structure data

Eukaryotic genomes are organized in a three-dimensional spatial structure. In this regard, the development of chromosome conformation capture methods has enabled studies of chromosome organization on a genomic scale. Hi-C, the high-throughput chromosome conformation capture method, can reveal a population-averaged, hierarchical chromatin structure. The typical Hi-C analysis uses a two-dimensional (2D) contact matrix that indicates contact frequencies between all possible genomic position pairs. Oftentimes, however, such a 2D matrix is not amenable to handling quantitative comparisons, visualizations and integrations across multiple datasets. Although several one-dimensional (1D) metrics have been proposed to depict structural information in Hi-C data, their effectiveness is still underappreciated. Here, we first review the currently available 1D metrics for individual Hi-C samples or two-sample comparisons and then discuss their validity and suitable analysis scenarios. We also propose several new 1D metrics to identify additional unique features of chromosome structures. We highlight that the 1D metrics are reproducible and robust for comparing and visualizing multiple Hi-C samples. Moreover, we show that 1D metrics can be easily combined with epigenome tracks to annotate chromatin states in greater details. We develop a new framework, called HiC1Dmetrics, to summarize all 1D metrics discussed in this study. HiC1Dmetrics is open-source (github.com/wangjk321/HiC1Dmetrics) and can be accessed from both command-line and web-based interfaces. Our tool constitutes a useful resource for the community of chromosome-organization researchers.

A compendium of chromatin contact maps reflecting regulation by chromatin remodelers in budding yeast

We herein employ in situ Hi-C with an auxin-inducible degron (AID) system to examine the effect of chromatin remodeling on 3D genome organization in yeast. Eight selected ATP-dependent chromatin remodelers representing various subfamilies contribute to 3D genome organization differently. Among the studied remodelers, the temporary depletions of Chd1p, Swr1p, and Sth1p (a catalytic subunit of the Remodeling the Structure of Chromatin [RSC] complex) cause the most significant defects in intra-chromosomal contacts, and the regulatory roles of these three remodelers in 3D genome organization differ depending on the chromosomal context and cell cycle stage. Furthermore, even though Chd1p and Isw1p are known to share functional similarities/redundancies, their depletions lead to distinct effects on 3D structures. The RSC and cohesin complexes also differentially modulate 3D genome organization within chromosome arm regions, whereas RSC appears to support the function of cohesin in centromeric clustering at G2 phase. Our work suggests that the ATP-dependent chromatin remodelers control the 3D genome organization of yeast through their chromatin-remodeling activities.

Cohesin regulates homology search during recombinational DNA repair

Homologous recombination repairs DNA double-strand breaks (DSB) using an intact dsDNA molecule as a template. It entails a homology search step, carried out along a conserved RecA/Rad51-ssDNA filament assembled on each DSB end. Whether, how and to what extent a DSB impacts chromatin folding, and how this (re)organization in turns influences the homology search process, remain ill-defined. Here we characterize two layers of spatial chromatin reorganization following DSB formation in Saccharomyces cerevisiae. Although cohesin folds chromosomes into cohesive arrays of ~20-kb-long chromatin loops as cells arrest in G2/M, the DSB-flanking regions interact locally in a resection- and 9-1-1 clamp-dependent manner, independently of cohesin, Mec1ATR, Rad52 and Rad51. This local structure blocks cohesin progression, constraining the DSB region at the base of a loop. Functionally, cohesin promotes DSB-dsDNA interactions and donor identification in cis, while inhibiting them in trans. This study identifies multiple direct and indirect ways by which cohesin regulates homology search during recombinational DNA repair.

Generation of gene-level resolution chromosome contact maps in bacteria and archaea

Chromosome conformation capture (Hi-C) has become a routine method for probing the 3D organization of genomes. However, when applied to bacteria and archaea, current protocols are expensive and limited in their resolution. By dissecting the different steps of published eukaryotic and prokaryotic Hi-C protocols, we have developed a cost- and time-effective approach to generate high-resolution (down to 500 bp – 1 kb) contact matrices of both bacteria and archaea genomes.

For complete details on the use and execution of this protocol, please refer to Cockram et al. (2020).

•

Optimized Hi-C protocol for archaeal and bacterial genomes

•

Generation of genome-wide contact maps up to 1 kb resolution

•

Detailed description of steps from cell fixation to sequencing library preparation

•

A cost- and time-effective approach offering gene-level resolution contact maps

Multi-scale architecture of archaeal chromosomes

Chromosome conformation capture (3C) technologies have identified topologically associating domains (TADs) and larger A/B compartments as two salient structural features of eukaryotic chromosomes. These structures are sculpted by the combined actions of transcription and structural maintenance of chromosomes (SMC) superfamily proteins. Bacterial chromosomes fold into TAD-like chromosomal interaction domains (CIDs) but do not display A/B compartment-type organization. We reveal that chromosomes of Sulfolobus archaea are organized into CID-like topological domains in addition to previously described larger A/B compartment-type structures. We uncover local rules governing the identity of the topological domains and their boundaries. We also identify long-range loop structures and provide evidence of a hub-like structure that colocalizes genes involved in ribosome biogenesis. In addition to providing high-resolution descriptions of archaeal chromosome architectures, our data provide evidence of multiple modes of organization in prokaryotic chromosomes and yield insights into the evolution of eukaryotic chromosome conformation.

Euryarchaeal genomes are folded into SMC-dependent loops and domains, but lack transcription-mediated compartmentalization

Hi-C has become a routine method for probing the 3D organization of genomes. However, when applied to prokaryotes and archaea, the current protocols are expensive and limited in their resolution. We develop a cost-effective Hi-C protocol to explore chromosome conformations of these two kingdoms at the gene or operon level. We first validate it on E. coli and V. cholera, generating sub-kilobase-resolution contact maps, and then apply it to the euryarchaeota H. volcanii, Hbt. salinarum, and T. kodakaraensis. With a resolution of up to 1 kb, we explore the diversity of chromosome folding in this phylum. In contrast to crenarchaeota, these euryarchaeota lack (active/inactive) compartment-like structures. Instead, their genomes are composed of self-interacting domains and chromatin loops. In H. volcanii, these structures are regulated by transcription and the archaeal structural maintenance of chromosomes (SMC) protein, further supporting the ubiquitous role of these processes in shaping the higher-order organization of genomes.