ScHiCEmbed: Bin-Specific Embeddings of Single-Cell Hi-C Data Using Graph Auto-Encoders

TADs and Their Borders: Free Movement or Building a Wall?

Recently a biochemistry experiment named methyl-3C was developed to simultaneously capture the chromosomal conformations and DNA methylation levels on individual single cells. However, the number of data sets generated from this experiment is still small in the scientific community compared with the greater amount of single-cell Hi-C data generated from separate single cells. Therefore, a computational tool to predict single-cell methylation levels based on single-cell Hi-C data on the same individual cells is needed. We developed a graph transformer named scHiMe to accurately predict the base-pair-specific (bp-specific) methylation levels based on both single-cell Hi-C data and DNA nucleotide sequences. We benchmarked scHiMe for predicting the bp-specific methylation levels on all of the promoters of the human genome, all of the promoter regions together with the corresponding first exon and intron regions, and random regions on the whole genome. Our evaluation showed a high consistency between the predicted and methyl-3C-detected methylation levels. Moreover, the predicted DNA methylation levels resulted in accurate classifications of cells into different cell types, which indicated that our algorithm successfully captured the cell-to-cell variability in the single-cell Hi-C data. scHiMe is freely available at http://dna.cs.miami.edu/scHiMe/.

Animal genomes are organized into topologically associated domains (TADs). TADs are thought to contribute to gene regulation by facilitating enhancer-promoter (E-P) contacts within a TAD and preventing these contacts across TAD borders. However, the absolute difference in contact frequency across TAD boundaries is usually less than 2-fold, even though disruptions of TAD borders can change gene expression by 10-fold. Existing models fail to explain this hypersensitive response. Here, we propose a futile cycle model of enhancer-mediated regulation that can exhibit hypersensitivity through bistability and hysteresis. Consistent with recent experiments, this regulation does not exhibit strong correlation between E-P contact and promoter activity, even though regulation occurs through contact. Through mathematical analysis and stochastic simulation, we show that this system can create an illusion of E-P biochemical specificity and explain the importance of weak TAD boundaries. It also offers a mechanism to reconcile apparently contradictory results from recent global TAD disruption with local TAD boundary deletion experiments. Together, these analyses advance our understanding of cis-regulatory contacts in controlling gene expression and suggest new experimental directions.

BackgroundHi-C data have been widely used to reconstruct chromosomal three-dimensional (3D) structures. One of the key limitations of Hi-C is the unclear relationship between spatial distance and the number of Hi-C contacts. Many methods used a fixed parameter when converting the number of Hi-C contacts to wish distances. However, a single parameter cannot properly explain the relationship between wish distances and genomic distances or the locations of topologically associating domains (TADs).ResultsWe have addressed one of the key issues of using Hi-C data, that is, the unclear relationship between spatial distances and the number of Hi-C contacts, which is crucial to understand significant biological functions, such as the enhancer-promoter interactions. Specifically, we developed a new method to infer this converting parameter and pairwise Euclidean distances based on the topology of the Hi-C complex network (HiCNet). The inferred distances were modeled by clustering coefficient and multiple other types of constraints. We found that our inferred distances between bead-pairs within the same TAD were apparently smaller than those distances between bead-pairs from different TADs. Our inferred distances had a higher correlation with fluorescence in situ hybridization (FISH) data, fitted the localization patterns of Xist transcripts on DNA, and better matched 156 pairs of protein-enabled long-range chromatin interactions detected by ChIA-PET. Using the inferred distances and another round of optimization, we further reconstructed 40 kb high-resolution 3D chromosomal structures of mouse male ES cells. The high-resolution structures successfully illustrate TADs and DNA loops (peaks in Hi-C contact heatmaps) that usually indicate enhancer-promoter interactions.ConclusionsWe developed a novel method to infer the wish distances between DNA bead-pairs from Hi-C contacts. High-resolution 3D structures of chromosomes were built based on the newly-inferred wish distances. This whole process has been implemented as a tool named HiCNet, which is publicly available at http://dna.cs.miami.edu/HiCNet/.

The hierarchical folding dynamics of topologically associating domains are closely related to transcriptional abnormalities in cancers

The evolution of gene regulation is considered one of the main drivers causing the astonishing morphological diversity in the animal kingdom. Gene regulation in animals heavily depends upon cis-regulatory elements, discrete pieces of DNA that interact with target promoters to regulate gene expression. In the last years, Chromosome Conformation Capture experiments (4C-seq, 5C, and HiC) in several organisms have shown that the genomes of many bilaterian animals are organized in the 3D chromatin space in compartments called topologically associated domains (TADs). The appearance of the architectural protein CTCF in the bilaterian ancestor likely facilitated the origin of this chromatin 3D organization. TADs play a critical role favoring the contact of cis-regulatory elements with their proper target promoters (that often lay within the same TAD) and preventing undesired regulatory interactions with promoters located in neighboring TADs. We propose that TAD may have had a major influence in the history of the evolution of gene regulation. They have contributed to the increment of regulatory complexity in bilaterians by allowing newly evolved cis-regulatory elements to find target promoters in a range of hundreds of kilobases. In addition, they have conditioned the mechanisms of evolution of gene regulation. These mechanisms include the appearance, removal, or relocation of TAD borders. Such architectural changes have been able to wire or unwire promoters with different sets of cis-regulatory elements in a single mutational event. We discuss the contribution of these architectural changes to the generation of critical genomic 3D structures required for new regulatory mechanisms associated to morphological novelties. WIREs Dev Biol 2017, 6:e265. doi: 10.1002/wdev.265 For further resources related to this article, please visit the WIREs website.

Studying the three-dimensional (3D) structure of a genome, including chromatin loops and Topologically Associating Domains (TADs), is essential for understanding how the genome is organized, such as gene activation, cell development, protein-protein interaction, etc. Hi-C protocol enables us to study 3D genome structure and organization. Chromatin 3D structure changes dynamically over time, and modeling these continuous changes is crucial for downstream analysis in various domains such as disease diagnosis, vaccine development, etc. The high expense and impracticality of continuous genome sequencing, particularly what evolves between two timestamps, limit the most effective genomic analysis. It is crucial to develop a straightforward and cost-efficient method for constantly generating genomic data between two timestamps in order to address these constraints. In this study, we developed HiCInterpolate, a 4D spatiotemporal interpolation architecture that accepts two timestamp Hi-C contact matrices to interpolate intermediate Hi-C contact matrices at high resolution. HiCInterpolate predicts the intermediate Hi-C contact map using a deep learning-based flow predictor, and a feature encoder and decoder architecture similar to U-Net. In addition, HiCInterpolate supports downstream analysis of multiple 3D genomic features, including A/B compartments, chromatin loops, TADs, and 3D genome structure, through an integrated analysis pipeline. Across multiple evaluation metrics, including PSNR, SSIM, GenomeDISCO, HiCRep, and LPIPS, HiCInterpolate achieved consistently strong performance. Biological validation further demonstrated preservation of key chromatin organization features, such as chromatin loops, A/B compartments, and TADs. Together, these results indicate that HiCInterpolate provides a robust computer vision-based framework for high-resolution interpolation of intermediate Hi-C contact matrices and facilitates biologically meaningful downstream analyses. HiCInterpolate is publicly available at https://github.com/OluwadareLab/HiCInterpolate .

We describe an approach for reconstructing three-dimensional (3D) structures from single-cell Hi-C data. This approach has been inspired by a method of recurrence plots and visualization tools for nonlinear time series data. Some examples are also presented.

The interphase chromatin is folded in the nucleus in a hierarchical manner, including the nucleosome, the "beads on a string" structure composed of nucleosomes, the solenoid fiber structure, the chromatin/DNA loop structure (chromatin/DNA loop), and the topologically associated domain (TAD). Among them, TAD is considered to be the basic unit of the 3D structure of chromatin because it is relatively stable and conserved in different cell types. Alu elements occupy a large proportion in the mammalian genomes. There are a wide variety of Alu elements, but their functional characterizations are limited to date. This study investigates the role of Alu elements in the assembly of 3D chromatin conformation. The evolutionary process of the Alu subfamily was explored by the distance relationship of the 3D structure of chromatin. We found that the proportion of Alu elements in high-density chromatin interaction increased with higher similarity, indicating that Alu plays an important role in the construction of chromatin 3D structure. There is a certain positive correlation between the strength of the upper interaction and the evolutionary relationship. In sum, the Alu elements with relatively close distances in the 1D sequence will also be close to each other in the 3D structure of chromatin.

Hi-C is a powerful method to investigate genome-wide, higher-order chromatin and chromosome conformations averaged from a population of cells. To expand the potential of Hi-C for single-cell analysis, we developed single-cell Hi-C. Similar to the existing "ensemble" Hi-C method, single-cell Hi-C detects proximity-dependent ligation events between cross-linked and restriction-digested chromatin fragments in cells. A major difference between the single-cell Hi-C and ensemble Hi-C protocol is that the proximity-dependent ligation is carried out in the nucleus. This allows the isolation of individual cells in which nearly the entire Hi-C procedure has been carried out, enabling the production of a Hi-C library and data from individual cells. With this new method, we studied genome conformations and found evidence for conserved topological domain organization from cell to cell, but highly variable interdomain contacts and chromosome folding genome wide. In addition, we found that the single-cell Hi-C protocol provided cleaner results with less technical noise suggesting it could be used to improve the ensemble Hi-C technique.

Eukaryotic genomes are organized within nuclei in three-dimensional space, forming structures such as loops, topologically associating domains (TADs), and chromosome territories. This 3D architecture impacts gene regulation and development, stress responses, and disease. However, current methods to infer these 3D structures from genomic data have multiple drawbacks, including varying outcomes depending on the resolution of the analysis and sequencing depth, qualitative outputs that limit statistical comparisons, and insufficient insight into structure frequency within samples. These challenges hinder rigorous comparisons of 3D properties across genomes, conditions, or species. To overcome these issues, we developed WaveTAD, a wavelet transform-based method that provides a resolution-free, probabilistic, and hierarchical description of 3D organization. WaveTAD generates TAD strengths, capturing the variable frequency of intrachromosomal interactions within samples, and shows increased accuracy and sensitivity over existing methods. We applied WaveTAD to multiple datasets from Drosophila, mouse, and humans to illustrate new biological insights that our more sensitive and quantitative approach provides, such as the widespread presence of embryonic 3D organization before zygotic genome activation, the effect of multiple CTCF units on the stability of loops and TADs, and the association between gene expression and TAD structures in COVID-19 patients or sex-specific transcription in Drosophila.

Single-cell Hi-C technology is utilized to obtain chromosome interaction information at the single-cell level and further study the differences in genome structures between different cell types. However, there are few accurate and efficient clustering methods for single-cell Hi-C data, with the following manifestations: The clustering efficacy on the dataset with a large number of cells is not very satisfactory, and it is difficult to cluster these cells with small number in the whole dataset. In this study, we propose a high-performance single-cell Hi-C clustering framework, called scHiCSC. A new smoothing method based on contact number weight is first proposed to generate cell embedding with more accurate cell features. In addition, a novel feature fusion method is proposed to further supplement the feature information of cells by fusing the chromosome structure information within cells and the distance information between cells. The experimental results show that scHiCSC has a strong generalization ability on different sizes of datasets and outperforms the existing single-cell Hi-C clustering frameworks. Moreover, scHiCSC achieves an optimal and stable clustering efficacy in the datasets with large-scale cell numbers and can cluster the cells with small number in the whole dataset. The source code of scHiCSC is freely available at https://github.com/HaoWuLab-Bioinformatics/scHiCSC.

The comprehensive annotation of regulatory elements in linear genomes is needed to elucidate the molecular mechanisms underlying chromatin loop formation in plants. Here, we characterized a novel family of conserved noncoding sequences (CNSs) in the rice (Oryza sativa) genome. These sequences, known as AT-rich pincer-like elements (APEs), are composed of 13-bp repeat unit arrays in a reverse-forward configuration. Our findings revealed that there are 611 APE copies across the japonica genome. Deletion of single APEs disrupted the long-range chromatin loops anchoring target-APE regions and moderately remodeled the profile of A/B compartments, topologically associating domains (TADs), and chromatin loops, thereby rewiring the expression of looped gene(s) including those controlling important agronomic traits. Thus, APEs function as hub motifs directly mediating chromatin looping and maintaining 3D genome integrity and stability at the levels of compartments, TADs, and loops. Moreover, neighboring genomic regions harboring numerous paired non-APE (NA) CNSs were more likely to interact with each other. This finding suggests that NA CNS pairs might play a helper role in determining loop frequency in a dose-dependent manner, likely by ensuring the pairing selectivity of anchor sites. Our study highlights the importance of APEs and NA CNSs in maintaining 3D genome structure, thereby providing the framework required to link many noncoding repetitive elements to their molecular functions in plants.

Chromatin conformation capture methods such as Hi-C have improved understanding of nuclear architecture. However, reconstruction from single-cell Hi-C (scHi-C) data is challenging due to limited DNA contacts per cell. We have previously developed the recurrence plot-based reconstruction (RPR) method for reconstructing three-dimensional (3D) genomic structure from Hi-C data even from low-coverage DNA contact information. Here we used the RPR method to analyze scHi-C data derived from early-stage F1 hybrid embryos as a proof-of-concept for understanding of global chromosomal architecture. We found that paternal and maternal chromosomes become gradually intermingled from the 1-cell to the 64-cell stage, and that discrete chromosome territories are largely established between 8-cell and 64-cell stages. We also observed Rabl-like polarization of chromosomes from the 2- to 8-cell stage, which was mostly dissolved by the 64-cell stage. We also noted transient rod-like extension and parallel chromosome alignment at the 4-cell stage. These findings indicate dynamic chromosomal reorganization before territory formation. RPR and scHi-C together capture 3D chromosomal architecture of individual cells during early embryogenesis.

In this study, by exploring chromatin conformation capture data, we show that the nuclear segregation of Topologically Associated Domains (TADs) is contributed by DNA sequence composition. GC-peaks and valleys of TADs strongly influence interchromosomal interactions and chromatin 3D structure. To gain insight on the compositional and functional constraints associated with chromatin interactions and TADs formation, we analysed intra-TAD and intra-loop GC variations. This led to the identification of clear GC-gradients, along which, the density of genes, super-enhancers, transcriptional activity, and CTCF binding sites occupancy co-vary non-randomly. Further, the analysis of DNA base composition of nucleolar aggregates and nuclear speckles showed strong sequence-dependant effects. We conjecture that dynamic DNA binding affinity and flexibility underlay the emergence of chromatin condensates, their growth is likely promoted in mechanically soft regions (GC-rich) of the lowest chromatin and nucleosome densities. As a practical perspective, the strong linear association between sequence composition and interchromosomal contacts can help define consensus chromatin interactions, which in turn may be used to study alternative states of chromatin architecture.