Calculation of 3D genome structures for comparison of chromosome conformation capture experiments with microscopy: An evaluation of single-cell Hi-C protocols

Dynamic changes in 3D chromatin structure during male gametogenesis in Arabidopsis thaliana

BACKGROUND

Chromatin higher-order structure plays an important role in genome stability maintenance and gene transcriptional regulation; however, the dynamics of the three-dimensional (3D) chromatin in male gametophytes during the two rounds of mitosis remains elusive.

RESULTS

Here, we use the optimized single-nucleus and low-input Hi-C methods to investigate changes in 3D chromatin structure in four types of male gametophyte nucleus at different stages. The reconstructed genome structures show that microspore nuclei develop towards two different directions. Although the 3D chromatin organization in generative nuclei is similar to that in microspore nuclei, vegetative nuclei lose chromosome territories, display dispersed centromeres, and switched A/B compartments, which are associated with vegetative specific gene expression. Additionally, we find that there is an active transcriptional center in sperm nuclei, emphasizing the transcription in Arabidopsis sperm is not completely inhibited despite the chromosomes being condensed.

CONCLUSIONS

Our data suggest that the special 3D structures of vegetative and sperm nuclei contribute to cell type-specific expression patterns.

CTPredictor: A comprehensive and robust framework for predicting cell types by integrating multi-scale features from single-cell Hi-C data

Single-cell Hi-C (scHi-C) has emerged as a powerful technology for deciphering cell-to-cell variability in three-dimensional (3D) chromatin organization, providing insights into genome-wide chromatin interactions and their correlation with cellular functions. Nevertheless, the accurate identification of cell types across different datasets remains a formidable challenge, hindering comprehensive investigations into genome structure. In response, we introduce CTPredictor, an innovative computational method that integrates multi-scale features to accurately predict cell types in various datasets. CTPredictor strategically incorporates three distinct feature sets, namely, small intra-domain contact probability (SICP), smoothed small intra-domain contact probability (SSICP), and smoothed bin contact probability (SBCP). The resulting fusion classification model significantly enhances the accuracy of cell type prediction based on single-cell Hi-C data (scHi-C). Rigorous benchmarking against established methods and three conventional machine learning approaches demonstrates the robust performance of CTPredictor, positioning it as an advanced tool for cell type prediction within scHi-C data. Beyond its prediction capabilities, CTPredictor holds promise in illuminating 3D genome structures and their functional significance across a wide array of biological processes.

Robust 3D modeling reveals spatiosyntenic properties of animal genomes

Animal genomes are organized into chromosomes that are remarkably conserved in their gene content, forming distinct evolutionary units (synteny). Using versatile chromosomal modeling, we infer three-dimensional topology of genomes from representative clades spanning the earliest animal diversification. We apply a partitioning approach using interaction spheres to compensate for varying quality of topological data. Using comparative genomics approaches, we test whether syntenic signal at gene pair, local, and whole chromosomal scale is reflected in the reconstructed spatial organization. We identify evolutionarily conserved three-dimensional networks at all syntenic scales revealing novel evolutionarily maintained interactors associated with known conserved local gene linkages (such as hox). We thus present evidence for evolutionary constraints that are associated with three-, rather than just two-, dimensional animal genome organization, which we term spatiosynteny. As more accurate topological data become available, together with validation approaches, spatiosynteny may become relevant in understanding the functionality behind the observed conservation of animal chromosomes.

Every gene everywhere all at once: High-precision measurement of 3D chromosome architecture with single-cell Hi-C

The three-dimensional (3D) structure of chromosomes influences essential biological processes such as gene expression, genome replication, and DNA damage repair and has been implicated in many developmental and degenerative diseases. In the past two centuries, two complementary genres of technology-microscopy, such as fluorescence in situ hybridization (FISH), and biochemistry, such as chromosome conformation capture (3C or Hi-C)-have revealed general principles of chromosome folding in the cell nucleus. However, the extraordinary complexity and cell-to-cell variability of the chromosome structure necessitate new tools with genome-wide coverage and single-cell precision. In the past decade, single-cell Hi-C emerges as a new approach that builds upon yet conceptually differs from bulk Hi-C assays. Instead of measuring population-averaged statistical properties of chromosome folding, single-cell Hi-C works as a proximity-based "biochemical microscope" that measures actual 3D structures of individual genomes, revealing features hidden in bulk Hi-C such as radial organization, multi-way interactions, and chromosome intermingling. Single-cell Hi-C has been used to study highly dynamic processes such as the cell cycle, cell-type-specific chromosome architecture ("structure types"), and structure-expression interplay, deepening our understanding of DNA organization and function.

3D genome organization in the epithelial-mesenchymal transition spectrum

Background

The plasticity along the epithelial-mesenchymal transition (EMT) spectrum has been shown to be regulated by various epigenetic repertoires. Emerging evidence of local chromatin conformation changes suggests that regulation of EMT may occur at a higher order of three-dimensional genome level.

Results

We perform Hi-C analysis and combine ChIP-seq data across cancer cell lines representing different EMT states. We demonstrate that the epithelial and mesenchymal genes are regulated distinctively. We find that EMT genes are regulated within their topologically associated domains (TADs), with only a subset of mesenchymal genes being influenced by A/B compartment switches, indicating topological remodeling is required in the transcriptional regulation of these genes. At the TAD level, epithelial and mesenchymal genes are associated with different regulatory trajectories. The epithelial gene-residing TADs are enriched with H3K27me3 marks in the mesenchymal-like states. The mesenchymal gene-residing TADs, which do not show enrichment of H3K27me3 in epithelial-like states, exhibit increased interaction frequencies with regulatory elements in the mesenchymal-like states.

Conclusions

We propose a novel workflow coupling immunofluorescence and dielectrophoresis to unravel EMT heterogeneity at single-cell resolution. The predicted three-dimensional structures of chromosome 10, harboring Vimentin, identify cell clusters of different states. Our results pioneer a novel avenue to decipher the complexities underlying the regulation of EMT and may infer the barriers of plasticity in the 3D genome context.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02687-x.

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

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

Single-cell Hi-C data analysis: safety in numbers

Over the past decade, genome-wide assays for chromatin interactions in single cells have enabled the study of individual nuclei at unprecedented resolution and throughput. Current chromosome conformation capture techniques survey contacts for up to tens of thousands of individual cells, improving our understanding of genome function in 3D. However, these methods recover a small fraction of all contacts in single cells, requiring specialised processing of sparse interactome data. In this review, we highlight recent advances in methods for the interpretation of single-cell genomic contacts. After discussing the strengths and limitations of these methods, we outline frontiers for future development in this rapidly moving field.

scHiCStackL: a stacking ensemble learning-based method for single-cell Hi-C classification using cell embedding

Single-cell Hi-C data are a common data source for studying the differences in the three-dimensional structure of cell chromosomes. The development of single-cell Hi-C technology makes it possible to obtain batches of single-cell Hi-C data. How to quickly and effectively discriminate cell types has become one hot research field. However, the existing computational methods to predict cell types based on Hi-C data are found to be low in accuracy. Therefore, we propose a high accuracy cell classification algorithm, called scHiCStackL, based on single-cell Hi-C data. In our work, we first improve the existing data preprocessing method for single-cell Hi-C data, which allows the generated cell embedding better to represent cells. Then, we construct a two-layer stacking ensemble model for classifying cells. Experimental results show that the cell embedding generated by our data preprocessing method increases by 0.23, 1.22, 1.46 and 1.61$\%$ comparing with the cell embedding generated by the previously published method scHiCluster, in terms of the Acc, MCC, F1 and Precision confidence intervals, respectively, on the task of classifying human cells in the ML1 and ML3 datasets. When using the two-layer stacking ensemble framework with the cell embedding, scHiCStackL improves by 13.33, 19, 19.27 and 14.5 over the scHiCluster, in terms of the Acc, ARI, NMI and F1 confidence intervals, respectively. In summary, scHiCStackL achieves superior performance in predicting cell types using the single-cell Hi-C data. The webserver and source code of scHiCStackL are freely available at http://hww.sdu.edu.cn:8002/scHiCStackL/ and https://github.com/HaoWuLab-Bioinformatics/scHiCStackL, respectively.

Si-C is a method for inferring super-resolution intact genome structure from single-cell Hi-C data

There is a strong demand for methods that can efficiently reconstruct valid super-resolution intact genome 3D structures from sparse and noise single-cell Hi-C data. Here, we develop Single-Cell Chromosome Conformation Calculator (Si-C) within the Bayesian theory framework and apply this approach to reconstruct intact genome 3D structures from single-cell Hi-C data of eight G1-phase haploid mouse ES cells. The inferred 100-kb and 10-kb structures consistently reproduce the known conserved features of chromatin organization revealed by independent imaging experiments. The analysis of the 10-kb resolution 3D structures reveals cell-to-cell varying domain structures in individual cells and hyperfine structures in domains, such as loops. An average of 0.2 contact reads per divided bin is sufficient for Si-C to obtain reliable structures. The valid super-resolution structures constructed by Si-C demonstrate the potential for visualizing and investigating interactions between all chromatin loci at the genome scale in individual cells.

Constructing valid super-resolution intact genome 3D structures from single-cell Hi-C data is essential in investigating chromosome folding. Here the authors develop a method that makes it possible to visualize and investigate chromosome folding in individual cells at the genome scale

Comparison of the Hi-C, GAM and SPRITE methods using polymer models of chromatin

Hi-C, split-pool recognition of interactions by tag extension (SPRITE) and genome architecture mapping (GAM) are powerful technologies utilized to probe chromatin interactions genome wide, but how faithfully they capture three-dimensional (3D) contacts and how they perform relative to each other is unclear, as no benchmark exists. Here, we compare these methods in silico in a simplified, yet controlled, framework against known 3D structures of polymer models of murine and human loci, which can recapitulate Hi-C, GAM and SPRITE experiments and multiplexed fluorescence in situ hybridization (FISH) single-molecule conformations. We find that in silico Hi-C, GAM and SPRITE bulk data are faithful to the reference 3D structures whereas single-cell data reflect strong variability among single molecules. The minimal number of cells required in replicate experiments to return statistically similar contacts is different across the technologies, being lowest in SPRITE and highest in GAM under the same conditions. Noise-to-signal levels follow an inverse power law with detection efficiency and grow with genomic distance differently among the three methods, being lowest in GAM for genomic separations >1 Mb.

Principal curve approaches for inferring 3D chromatin architecture

Three-dimensional (3D) genome spatial organization is critical for numerous cellular processes, including transcription, while certain conformation-driven structural alterations are frequently oncogenic. Genome architecture had been notoriously difficult to elucidate, but the advent of the suite of chromatin conformation capture assays, notably Hi-C, has transformed understanding of chromatin structure and provided downstream biological insights. Although many findings have flowed from direct analysis of the pairwise proximity data produced by these assays, there is added value in generating corresponding 3D reconstructions deriving from superposing genomic features on the reconstruction. Accordingly, many methods for inferring 3D architecture from proximity data have been advanced. However, none of these approaches exploit the fact that single chromosome solutions constitute a one-dimensional (1D) curve in 3D. Rather, this aspect has either been addressed by imposition of constraints, which is both computationally burdensome and cell type specific, or ignored with contiguity imposed after the fact. Here, we target finding a 1D curve by extending principal curve methodology to the metric scaling problem. We illustrate how this approach yields a sequence of candidate solutions, indexed by an underlying smoothness or degrees-of-freedom parameter, and propose methods for selection from this sequence. We apply the methodology to Hi-C data obtained on IMR90 cells and so are positioned to evaluate reconstruction accuracy by referencing orthogonal imaging data. The results indicate the utility and reproducibility of our principal curve approach in the face of underlying structural variation.

GiniQC: a measure for quantifying noise in single-cell Hi-C data

SUMMARY

Single-cell Hi-C (scHi-C) allows the study of cell-to-cell variability in chromatin structure and dynamics. However, the high level of noise inherent in current scHi-C protocols necessitates careful assessment of data quality before biological conclusions can be drawn. Here, we present GiniQC, which quantifies unevenness in the distribution of inter-chromosomal reads in the scHi-C contact matrix to measure the level of noise. Our examples show the utility of GiniQC in assessing the quality of scHi-C data as a complement to existing quality control measures. We also demonstrate how GiniQC can help inform the impact of various data processing steps on data quality.

AVAILABILITY AND IMPLEMENTATION

Source code and documentation are freely available at https://github.com/4dn-dcic/GiniQC.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The Dynamic 3D Genome in Gametogenesis and Early Embryonic Development

During gametogenesis and early embryonic development, the chromatin architecture changes dramatically, and both the transcriptomic and epigenomic landscape are comprehensively reprogrammed. Understanding these processes is the holy grail in developmental biology and a key step towards evolution. The 3D conformation of chromatin plays a central role in the organization and function of nuclei. Recently, the dynamics of chromatin structures have been profiled in many model and non-model systems, from insects to mammals, resulting in an interesting comparison. In this review, we first introduce the research methods of 3D chromatin structure with low-input material suitable for embryonic study. Then, the dynamics of 3D chromatin architectures during gametogenesis and early embryonic development is summarized and compared between species. Finally, we discuss the possible mechanisms for triggering the formation of genome 3D conformation in early development.

Single-Cell RNA Sequencing-Based Computational Analysis to Describe Disease Heterogeneity

The trillions of cells in the human body can be viewed as elementary but essential biological units that achieve different body states, but the low resolution of previous cell isolation and measurement approaches limits our understanding of the cell-specific molecular profiles. The recent establishment and rapid growth of single-cell sequencing technology has facilitated the identification of molecular profiles of heterogeneous cells, especially on the transcription level of single cells [single-cell RNA sequencing (scRNA-seq)]. As a novel method, the robustness of scRNA-seq under changing conditions will determine its practical potential in major research programs and clinical applications. In this review, we first briefly presented the scRNA-seq-related methods from the point of view of experiments and computation. Then, we compared several state-of-the-art scRNA-seq analysis frameworks mainly by analyzing their performance robustness on independent scRNA-seq datasets for the same complex disease. Finally, we elaborated on our hypothesis on consensus scRNA-seq analysis and summarized the potential indicative and predictive roles of individual cells in understanding disease heterogeneity by single-cell technologies.

The Pediatric Cell Atlas: Defining the Growth Phase of Human Development at Single-Cell Resolution

Single-cell gene expression analyses of mammalian tissues have uncovered profound stage-specific molecular regulatory phenomena that have changed the understanding of unique cell types and signaling pathways critical for lineage determination, morphogenesis, and growth. We discuss here the case for a Pediatric Cell Atlas as part of the Human Cell Atlas consortium to provide single-cell profiles and spatial characterization of gene expression across human tissues and organs. Such data will complement adult and developmentally focused HCA projects to provide a rich cytogenomic framework for understanding not only pediatric health and disease but also environmental and genetic impacts across the human lifespan.