Multiplexed imaging of nucleome architectures in single cells of mammalian tissue

ImputeHiFI: An Imputation Method for Multiplexed DNA FISH Data by Utilizing Single-Cell Hi-C and RNA FISH Data

Although multiplexed DNA fluorescence in situ hybridization (FISH) enables tracking the spatial localization of thousands of genomic loci using probes within individual cells, the high rates of undetected probes impede the depiction of 3D chromosome structures. Current data imputation methods neither utilize single-cell Hi-C data, which elucidate 3D genome architectures using sequencing nor leverage multimodal RNA FISH data that reflect cell-type information, limiting the effectiveness of these methods in complex tissues such as the mouse brain. To this end, a novel multiplexed DNA FISH imputation method named ImputeHiFI is proposed, which fully utilizes the complementary structural information from single-cell Hi-C data and the cell type signature from RNA FISH data to obtain a high-fidelity and complete spatial location of chromatin loci. ImputeHiFI enhances cell clustering, compartment identification, and cell subtype detection at the single-cell level in the mouse brain. ImputeHiFI improves the recognition of cell-type-specific loops in three high-resolution datasets. In short, ImputeHiFI is a powerful tool capable of imputing multiplexed DNA FISH data from various resolutions and imaging protocols, facilitating studies of 3D genome structures and functions.

Tracing the evolution of single-cell cancer 3D genomes: an atlas for cancer gene discovery

Although three-dimensional (3D) genome structures are altered in cancer cells, little is known about how these changes evolve and diversify during cancer progression. Leveraging genome-wide chromatin tracing to visualize 3D genome folding directly in tissues, we generated 3D genome cancer atlases of murine lung and pancreatic adenocarcinoma. Our data reveal stereotypical, non-monotonic, and stage-specific alterations in 3D genome folding heterogeneity, compaction, and compartmentalization as cancers progress from normal to preinvasive and ultimately to invasive tumors, discovering a potential structural bottleneck in early tumor progression. Remarkably, 3D genome architectures distinguish histologic cancer states in single cells, despite considerable cell-to-cell heterogeneity. Gene-level analyses of evolutionary changes in 3D genome compartmentalization not only showed compartment-associated genes are more homogeneously regulated, but also elucidated prognostic and dependency genes in lung adenocarcinoma and a previously unappreciated role for polycomb-group protein Rnf2 in 3D genome regulation. Our results demonstrate the utility of mapping the single-cell cancer 3D genome in tissues and illuminate its potential to identify new diagnostic, prognostic, and therapeutic biomarkers in cancer.

Chromosome compartmentalization: causes, changes, consequences, and conundrums

The spatial segregation of the genome into compartments is a major feature of 3D genome organization. New data on mammalian chromosome organization across different conditions reveal important information about how and why these compartments form and change. A combination of epigenetic state, nuclear body tethering, physical forces, gene expression, and replication timing (RT) can all influence the establishment and alteration of chromosome compartments. We review the causes and implications of genomic regions undergoing a 'compartment switch' that changes their physical associations and spatial location in the nucleus. About 20-30% of genomic regions change compartment during cell differentiation or cancer progression, whereas alterations in response to a stimulus within a cell type are usually much more limited. However, even a change in 1-2% of genomic bins may have biologically relevant implications. Finally, we review the effects of compartment changes on gene regulation, DNA damage repair, replication, and the physical state of the cell.

Effect of active loop extrusion on the two-contact correlations in the interphase chromosome

The population-averaged contact maps generated by the chromosome conformation capture technique provide important information about the average frequency of contact between pairs of chromatin loci as a function of the genetic distance between them. However, these datasets do not tell us anything about the joint statistics of simultaneous contacts between genomic loci in individual cells. This kind of statistical information can be extracted using the single-cell Hi-C method, which is capable of detecting a large fraction of simultaneous contacts within a single cell, as well as through modern methods of fluorescent labeling and super-resolution imaging. Motivated by the prospect of the imminent availability of relevant experimental data, in this work, we theoretically model the joint statistics of pairs of contacts located along a line perpendicular to the main diagonal of the single-cell contact map. The analysis is performed within the framework of an ideal polymer model with quenched disorder of random loops, which, as previous studies have shown, allows us to take into account the influence of the loop extrusion process on the conformational properties of interphase chromatin.

Enhancer selectivity in space and time: from enhancer-promoter interactions to promoter activation

The primary regulators of metazoan gene expression are enhancers, originally functionally defined as DNA sequences that can activate transcription at promoters in an orientation-independent and distance-independent manner. Despite being crucial for gene regulation in animals, what mechanisms underlie enhancer selectivity for promoters, and more fundamentally, how enhancers interact with promoters and activate transcription, remain poorly understood. In this Review, we first discuss current models of enhancer-promoter interactions in space and time and how enhancers affect transcription activation. Next, we discuss different mechanisms that mediate enhancer selectivity, including repression, biochemical compatibility and regulation of 3D genome structure. Through 3D polymer simulations, we illustrate how the ability of 3D genome folding mechanisms to mediate enhancer selectivity strongly varies for different enhancer-promoter interaction mechanisms. Finally, we discuss how recent technical advances may provide new insights into mechanisms of enhancer-promoter interactions and how technical biases in methods such as Hi-C and Micro-C and imaging techniques may affect their interpretation.

Single-cell mapping of cell-type specific chromatin architecture in the central nervous system

Determining how chromatin is structured in the nucleus is critical to studying its role in gene regulation. Recent advances in the analysis of single-cell chromatin architecture have considerably improved our understanding of cell-type-specific chromosome conformation and nuclear architecture. In this review, we discuss the methods used for analysis of 3D chromatin conformation, including sequencing-based methods, imaging-based techniques, and computational approaches. We further review the application of these methods in the study of the role of chromatin topology in neural development and disorders.

Emerging Roles of Spatial Transcriptomics in Liver Research

Spatial transcriptomics, leveraging sequencing- and imaging-based techniques, has emerged as a groundbreaking technology for mapping gene expression within the complex architectures of tissues. This approach provides an in-depth understanding of cellular and molecular dynamics across various states of healthy and diseased livers. Through the integration of sophisticated bioinformatics strategies, it enables detailed exploration of cellular heterogeneity, transitions in cell states, and intricate cell-cell interactions with remarkable precision. In liver research, spatial transcriptomics has been particularly revelatory, identifying distinct zonated functions of hepatocytes that are crucial for understanding the metabolic and detoxification processes of the liver. Moreover, this technology has unveiled new insights into the pathogenesis of liver diseases, such as the role of lipid-associated macrophages in steatosis and endothelial cell signals in liver regeneration and repair. In the domain of liver cancer, spatial transcriptomics has proven instrumental in delineating intratumor heterogeneity, identifying supportive microenvironmental niches and revealing the complex interplay between tumor cells and the immune system as well as susceptibility to immune checkpoint inhibitors. In conclusion, spatial transcriptomics represents a significant advance in hepatology, promising to enhance our understanding and treatment of liver diseases.

Advances in the multimodal analysis of the 3D chromatin structure and gene regulation

Recent studies have demonstrated that the three-dimensional conformation of the chromatin plays a crucial role in gene regulation, with aberrations potentially leading to various diseases. Advanced methodologies have revealed a link between the chromatin conformation and biological function. This review divides these methodologies into sequencing-based and imaging-based methodologies, tracing their development over time. We particularly highlight innovative techniques that facilitate the simultaneous mapping of RNAs, histone modifications, and proteins within the context of the 3D architecture of chromatin. This multimodal integration substantially improves our ability to establish a robust connection between the spatial arrangement of molecular components in the nucleus and their functional roles. Achieving a comprehensive understanding of gene regulation requires capturing diverse data modalities within individual cells, enabling the direct inference of functional relationships between these components. In this context, imaging-based technologies have emerged as an especially promising approach for gathering spatial information across multiple components in the same cell.

Footprints of loop extrusion in statistics of intra-chromosomal distances: An analytically solvable model

Active loop extrusion-the process of formation of dynamically growing chromatin loops due to the motor activity of DNA-binding protein complexes-is a firmly established mechanism responsible for chromatin spatial organization at different stages of a cell cycle in eukaryotes and bacteria. The theoretical insight into the effect of loop extrusion on the experimentally measured statistics of chromatin conformation can be gained with an appropriately chosen polymer model. Here, we consider the simplest analytically solvable model of an interphase chromosome, which is treated as an ideal chain with disorder of sufficiently sparse random loops whose conformations are sampled from the equilibrium ensemble. This framework allows us to arrive at the closed-form analytical expression for the mean-squared distance between pairs of genomic loci, which is valid beyond the one-loop approximation in diagrammatic representation. In addition, we analyze the loop-induced deviation of chain conformations from the Gaussian statistics by calculating kurtosis of probability density of the pairwise separation vector. The presented results suggest the possible ways of estimating the characteristics of the loop extrusion process based on the experimental data on the scale-dependent statistics of intra-chromosomal pair-wise distances.

Multiplex DNA fluorescence in situ hybridization to analyze maternal vs. paternal C. elegans chromosomes

Recent advances in microscopy have enabled studying chromosome organization at the single-molecule level, yet little is known about inherited chromosome organization. Here we adapt single-molecule chromosome tracing to distinguish two C. elegans strains (N2 and HI) and find that while their organization is similar, the N2 chromosome influences the folding parameters of the HI chromosome, in particular the step size, across generations. Furthermore, homologous chromosomes overlap frequently, but alignment between homologous regions is rare, suggesting that transvection is unlikely. We present a powerful tool to investigate chromosome architecture and to track the parent of origin.

Single-cell multiplex chromatin and RNA interactions in aging human brain

The dynamically organized chromatin complexes often involve multiplex chromatin interactions and sometimes chromatin-associated RNA (caRNA) 1-3. Chromatin complex compositions change during cellular differentiation and aging, and are expected to be highly heterogeneous among terminally differentiated single cells 4-7. Here we introduce the Multi-Nucleic Acid Interaction Mapping in Single Cell (MUSIC) technique for concurrent profiling of multiplex chromatin interactions, gene expression, and RNA-chromatin associations within individual nuclei. Applied to 14 human frontal cortex samples from elderly donors, MUSIC delineates diverse cortical cell types and states. We observed the nuclei exhibiting fewer short-range chromatin interactions are correlated with an "older" transcriptomic signature and with Alzheimer's pathology. Furthermore, the cell type exhibiting chromatin contacts between cis expression quantitative trait loci (cis eQTLs) and a promoter tends to be the cell type where these cis eQTLs specifically affect their target gene's expression. Additionally, the female cortical cells exhibit highly heterogeneous interactions between the XIST non-coding RNA and Chromosome X, along with diverse spatial organizations of the X chromosomes. MUSIC presents a potent tool for exploring chromatin architecture and transcription at cellular resolution in complex tissues.

pyHiM: a new open-source, multi-platform software package for spatial genomics based on multiplexed DNA-FISH imaging

Genome-wide ensemble sequencing methods improved our understanding of chromatin organization in eukaryotes but lack the ability to capture single-cell heterogeneity and spatial organization. To overcome these limitations, new imaging-based methods have emerged, giving rise to the field of spatial genomics. Here, we present pyHiM, a user-friendly python toolbox specifically designed for the analysis of multiplexed DNA-FISH data and the reconstruction of chromatin traces in individual cells. pyHiM employs a modular architecture, allowing independent execution of analysis steps and customization according to sample specificity and computing resources. pyHiM aims to facilitate the democratization and standardization of spatial genomics analysis.

Longitudinal single-cell transcriptional dynamics throughout neurodegeneration in SCA1

Neurodegeneration is a protracted process involving progressive changes in myriad cell types that ultimately results in the death of vulnerable neuronal populations. To dissect how individual cell types within a heterogeneous tissue contribute to the pathogenesis and progression of a neurodegenerative disorder, we performed longitudinal single-nucleus RNA sequencing of mouse and human spinocerebellar ataxia type 1 (SCA1) cerebellar tissue, establishing continuous dynamic trajectories of each cell population. Importantly, we defined the precise transcriptional changes that precede loss of Purkinje cells and, for the first time, identified robust early transcriptional dysregulation in unipolar brush cells and oligodendroglia. Finally, we applied a deep learning method to predict disease state accurately and identified specific features that enable accurate distinction of wild-type and SCA1 cells. Together, this work reveals new roles for diverse cerebellar cell types in SCA1 and provides a generalizable analysis framework for studying neurodegeneration.

Regulation of nuclear transcription by mitochondrial RNA in endothelial cells

Chromatin-associated RNAs (caRNAs) form a relatively poorly recognized layer of the epigenome. The caRNAs reported to date are transcribed from the nuclear genome. Here, leveraging a recently developed assay for detection of caRNAs and their genomic association, we report that mitochondrial RNAs (mtRNAs) are attached to the nuclear genome and constitute a subset of caRNA, thus termed mt-caRNA. In four human cell types analyzed, mt-caRNAs preferentially attach to promoter regions. In human endothelial cells (ECs), the level of mt-caRNA-promoter attachment changes in response to environmental stress that mimics diabetes. Suppression of a non-coding mt-caRNA in ECs attenuates stress-induced nascent RNA transcription from the nuclear genome, including that of critical genes regulating cell adhesion, and abolishes stress-induced monocyte adhesion, a hallmark of dysfunctional ECs. Finally, we report increased nuclear localization of multiple mtRNAs in the ECs of human diabetic donors, suggesting many mtRNA translocate to the nucleus in a cell stress and disease-dependent manner. These data nominate mt-caRNAs as messenger molecules responsible for mitochondrial-nuclear communication and connect the immediate product of mitochondrial transcription with the transcriptional regulation of the nuclear genome.

Hi-M: A Multiplex Oligopaint FISH Method to Capture Chromatin Conformations In Situ and Accompanying Open-Source Acquisition Software

The simultaneous observation of three-dimensional (3D) chromatin structure and transcription in single cells is critical to understand how DNA is organized inside cells and how this organization influences or is affected by other processes, such as transcription. We have recently introduced an innovative technology known as Hi-M, which enables the sequential tagging, 3D visualization, and precise localization of multiple genomic DNA regions alongside RNA expression within individual cells. In this chapter, we present a comprehensive guide outlining the creation of probes, as well as sample preparation and labeling. Finally, we provide a step-by-step guide to conduct a complete Hi-M acquisition using our open-source software package, Qudi-HiM, which controls the robotic microscope handling the entire acquisition procedure.

Interplay between epigenome and 3D chromatin structure

Epigenetic mechanisms, primarily mediated through histone and DNA modifications, play a pivotal role in orchestrating the functional identity of a cell and its response to environmental cues. Similarly, the spatial arrangement of chromatin within the threedimensional (3D) nucleus has been recognized as a significant factor influencing genomic function. Investigating the relationship between epigenetic regulation and 3D chromatin structure has revealed correlation and causality between these processes, from the global alignment of average chromatin structure with chromatin marks to the nuanced correlations at smaller scales. This review aims to dissect the biological significance and the interplay between the epigenome and 3D chromatin structure, while also exploring the underlying molecular mechanisms. By synthesizing insights from both experimental and modeling perspectives, we seek to provide a comprehensive understanding of cellular functions. [BMB Reports 2023; 56(12): 633-644].

Spatial promoter-enhancer hubs in cancer: organization, regulation, and function

Transcriptional dysregulation is a hallmark of cancer and can be driven by altered enhancer landscapes. Recent studies in genome organization have revealed that multiple enhancers and promoters can spatially coalesce to form dynamic topological assemblies, known as promoter-enhancer hubs, which strongly correlate with elevated gene expression. In this review, we discuss the structure and complexity of promoter-enhancer hubs recently identified in multiple cancer types. We further discuss underlying mechanisms driving dysregulation of promoter-enhancer hubs and speculate on their functional role in pathogenesis. Understanding the role of promoter-enhancer hubs in transcriptional dysregulation can provide insight into new therapeutic approaches to target these complex features of genome organization.

ChromaFactor: deconvolution of single-molecule chromatin organization with non-negative matrix factorization

The investigation of chromatin organization in single cells holds great promise for identifying causal relationships between genome structure and function. However, analysis of single-molecule data is hampered by extreme yet inherent heterogeneity, making it challenging to determine the contributions of individual chromatin fibers to bulk trends. To address this challenge, we propose ChromaFactor, a novel computational approach based on non-negative matrix factorization that deconvolves single-molecule chromatin organization datasets into their most salient primary components. ChromaFactor provides the ability to identify trends accounting for the maximum variance in the dataset while simultaneously describing the contribution of individual molecules to each component. Applying our approach to two single-molecule imaging datasets across different genomic scales, we find that these primary components demonstrate significant correlation with key functional phenotypes, including active transcription, enhancer-promoter distance, and genomic compartment. ChromaFactor offers a robust tool for understanding the complex interplay between chromatin structure and function on individual DNA molecules, pinpointing which subpopulations drive functional changes and fostering new insights into cellular heterogeneity and its implications for bulk genomic phenomena.

SnapFISH: a computational pipeline to identify chromatin loops from multiplexed DNA FISH data

Multiplexed DNA fluorescence in situ hybridization (FISH) imaging technologies have been developed to map the folding of chromatin fibers at tens of nanometers and up to several kilobases in resolution in single cells. However, computational methods to reliably identify chromatin loops from such imaging datasets are still lacking. Here we present a Single-Nucleus Analysis Pipeline for multiplexed DNA FISH (SnapFISH), to process the multiplexed DNA FISH data and identify chromatin loops. SnapFISH can identify known chromatin loops from mouse embryonic stem cells with high sensitivity and accuracy. In addition, SnapFISH obtains comparable results of chromatin loops across datasets generated from diverse imaging technologies. SnapFISH is freely available at https://github.com/HuMingLab/SnapFISH .

Female naïve human pluripotent stem cells carry X chromosomes with Xa-like and Xi-like folding conformations

Three-dimensional (3D) genomics shows immense promise for studying X chromosome inactivation (XCI) by interrogating changes to the X chromosomes' 3D states. Here, we sought to characterize the 3D state of the X chromosome in naïve and primed human pluripotent stem cells (hPSCs). Using chromatin tracing, we analyzed X chromosome folding conformations in these cells with megabase genomic resolution. X chromosomes in female naïve hPSCs exhibit folding conformations similar to the active X chromosome (Xa) and the inactive X chromosome (Xi) in somatic cells. However, naïve X chromosomes do not exhibit the chromatin compaction typically associated with these somatic X chromosome states. In H7 naïve human embryonic stem cells, XIST accumulation observed on damaged X chromosomes demonstrates the potential for naïve hPSCs to activate XCI-related mechanisms. Overall, our findings provide insight into the X chromosome status of naïve hPSCs with a single-chromosome resolution and are critical in understanding the unique epigenetic regulation in early embryonic cells.

Spatial and temporal organization of the genome: Current state and future aims of the 4D nucleome project

The four-dimensional nucleome (4DN) consortium studies the architecture of the genome and the nucleus in space and time. We summarize progress by the consortium and highlight the development of technologies for (1) mapping genome folding and identifying roles of nuclear components and bodies, proteins, and RNA, (2) characterizing nuclear organization with time or single-cell resolution, and (3) imaging of nuclear organization. With these tools, the consortium has provided over 2,000 public datasets. Integrative computational models based on these data are starting to reveal connections between genome structure and function. We then present a forward-looking perspective and outline current aims to (1) delineate dynamics of nuclear architecture at different timescales, from minutes to weeks as cells differentiate, in populations and in single cells, (2) characterize cis-determinants and trans-modulators of genome organization, (3) test functional consequences of changes in cis- and trans-regulators, and (4) develop predictive models of genome structure and function.

Methods and applications for single-cell and spatial multi-omics

The joint analysis of the genome, epigenome, transcriptome, proteome and/or metabolome from single cells is transforming our understanding of cell biology in health and disease. In less than a decade, the field has seen tremendous technological revolutions that enable crucial new insights into the interplay between intracellular and intercellular molecular mechanisms that govern development, physiology and pathogenesis. In this Review, we highlight advances in the fast-developing field of single-cell and spatial multi-omics technologies (also known as multimodal omics approaches), and the computational strategies needed to integrate information across these molecular layers. We demonstrate their impact on fundamental cell biology and translational research, discuss current challenges and provide an outlook to the future.

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

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

Genome-wide tracing to decipher nuclear organization

Nuclear organization impacts gene expression activity and cell phenotype. Our current understanding is mainly derived from ensemble-level sequencing studies that reflect the 3D genome structure of millions of cells. These approaches have provided invaluable details on the 3D organizations of the genome and their relation to other nuclear landmarks. However, they mostly lack the ability to provide multimodal information simultaneously at the single-cell level. In recent years, cutting-edge imaging technologies have risen to the challenge of simultaneously describing multiple components of the nuclear space at the single-cell level, paving the way for a deeper understanding of the genome structure-function relationship. This review will focus on the development and utilization of such technologies to gain a multi-component view of the nucleus at single-cell resolution, dissecting the complexity and heterogeneity of nuclear organization.

High-resolution spatial multi-omics reveals cell-type specific nuclear compartments

The mammalian nucleus is compartmentalized by diverse subnuclear structures. These subnuclear structures, marked by nuclear bodies and histone modifications, are often cell-type specific and affect gene regulation and 3D genome organization1-3. Understanding nuclear organization requires identifying the molecular constituents of subnuclear structures and mapping their associations with specific genomic loci in individual cells, within complex tissues. Here, we introduce two-layer DNA seqFISH+, which allows simultaneous mapping of 100,049 genomic loci, together with nascent transcriptome for 17,856 genes and a diverse set of immunofluorescently labeled subnuclear structures all in single cells in cell lines and adult mouse cerebellum. Using these multi-omics datasets, we showed that repressive chromatin compartments are more variable by cell type than active compartments. We also discovered a single exception to this rule: an RNA polymerase II (RNAPII)-enriched compartment was associated with long, cell-type specific genes (> 200kb), in a manner distinct from nuclear speckles. Further, our analysis revealed that cell-type specific facultative and constitutive heterochromatin compartments marked by H3K27me3 and H4K20me3 are enriched at specific genes and gene clusters, respectively, and shape radial chromosomal positioning and inter-chromosomal interactions in neurons and glial cells. Together, our results provide a single-cell high-resolution multi-omics view of subnuclear compartments, associated genomic loci, and their impacts on gene regulation, directly within complex tissues.

The magic of unraveling genome architecture and function

Over the last decades, technological breakthroughs in super-resolution microscopy have allowed us to reach molecular resolution and design experiments of unprecedented complexity. Investigating how chromatin is folded in 3D, from the nucleosome level up to the entire genome, is becoming possible by "magic" (imaging genomic), i.e., the combination of imaging and genomic approaches. This offers endless opportunities to delve into the relationship between genome structure and function. Here, we review recently achieved objectives and the conceptual and technical challenges the field of genome architecture is currently undertaking. We discuss what we have learned so far and where we are heading. We elucidate how the different super-resolution microscopy approaches and, more specifically, live-cell imaging have contributed to the understanding of genome folding. Moreover, we discuss how future technical developments could address remaining open questions.

Contributions of 3D chromatin structure to cell-type-specific gene regulation

Eukaryotic genomes are organized in 3D in a multiscale manner, and different mechanisms acting at each of these scales can contribute to transcriptional regulation. However, the large single-cell variability in 3D chromatin structures represents a challenge to understand how transcription may be differentially regulated between cell types in a robust and efficient manner. Here, we describe the different mechanisms by which 3D chromatin structure was shown to contribute to cell-type-specific transcriptional regulation. Excitingly, several novel methodologies able to measure 3D chromatin conformation and transcription in single cells in their native tissue context, or to detect the dynamics of cis-regulatory interactions, are starting to allow quantitative dissection of chromatin structure noise and relate it to how transcription may be regulated between different cell types and cell states.

3D genomics and its applications in precision medicine

Three-dimensional (3D) genomics is an emerging discipline that studies the three-dimensional structure of chromatin and the three-dimensional and functions of genomes. It mainly studies the three-dimensional conformation and functional regulation of intranuclear genomes, such as DNA replication, DNA recombination, genome folding, gene expression regulation, transcription factor regulation mechanism, and the maintenance of three-dimensional conformation of genomes. Self-chromosomal conformation capture (3C) technology has been developed, and 3D genomics and related fields have developed rapidly. In addition, chromatin interaction analysis techniques developed by 3C technologies, such as paired-end tag sequencing (ChIA-PET) and whole-genome chromosome conformation capture (Hi-C), enable scientists to further study the relationship between chromatin conformation and gene regulation in different species. Thus, the spatial conformation of plant, animal, and microbial genomes, transcriptional regulation mechanisms, interaction patterns of chromosomes, and the formation mechanism of spatiotemporal specificity of genomes are revealed. With the help of new experimental technologies, the identification of key genes and signal pathways related to life activities and diseases is sustaining the rapid development of life science, agriculture, and medicine. In this paper, the concept and development of 3D genomics and its application in agricultural science, life science, and medicine are introduced, which provides a theoretical basis for the study of biological life processes.

A maximum-entropy model to predict 3D structural ensembles of chromatin from pairwise distances with applications to interphase chromosomes and structural variants

The principles that govern the organization of genomes, which are needed for an understanding of how chromosomes are packaged and function in eukaryotic cells, could be deciphered if the three-dimensional (3D) structures are known. Recently, single-cell imaging techniques have been developed to determine the 3D coordinates of genomic loci in vivo. Here, we introduce a computational method (Distance Matrix to Ensemble of Structures, DIMES), based on the maximum entropy principle, with experimental pairwise distances between loci as constraints, to generate a unique ensemble of 3D chromatin structures. Using the ensemble of structures, we quantitatively account for the distribution of pairwise distances, three-body co-localization, and higher-order interactions. The DIMES method can be applied to both small and chromosome-scale imaging data to quantify the extent of heterogeneity and fluctuations in the shapes across various length scales. We develop a perturbation method in conjunction with DIMES to predict the changes in 3D structures from structural variations. Our method also reveals quantitative differences between the 3D structures inferred from Hi-C and those measured in imaging experiments. Finally, the physical interpretation of the parameters extracted from DIMES provides insights into the origin of phase separation between euchromatin and heterochromatin domains.

The spatial organization of transcriptional control

In animals, the sequences for controlling gene expression do not concentrate just at the transcription start site of genes, but are frequently thousands to millions of base pairs distal to it. The interaction of these sequences with one another and their transcription start sites is regulated by factors that shape the three-dimensional (3D) organization of the genome within the nucleus. Over the past decade, indirect tools exploiting high-throughput DNA sequencing have helped to map this 3D organization, have identified multiple key regulators of its structure and, in the process, have substantially reshaped our view of how 3D genome architecture regulates transcription. Now, new tools for high-throughput super-resolution imaging of chromatin have directly visualized the 3D chromatin organization, settling some debates left unresolved by earlier indirect methods, challenging some earlier models of regulatory specificity and creating hypotheses about the role of chromatin structure in transcriptional regulation.

Concordance of MERFISH spatial transcriptomics with bulk and single-cell RNA sequencing

Spatial transcriptomics extends single-cell RNA sequencing (scRNA-seq) by providing spatial context for cell type identification and analysis. Imaging-based spatial technologies such as multiplexed error-robust fluorescence in situ hybridization (MERFISH) can achieve single-cell resolution, directly mapping single-cell identities to spatial positions. MERFISH produces a different data type than scRNA-seq, and a technical comparison between the two modalities is necessary to ascertain how to best integrate them. We performed MERFISH on the mouse liver and kidney and compared the resulting bulk and single-cell RNA statistics with those from the Tabula Muris Senis cell atlas and from two Visium datasets. MERFISH quantitatively reproduced the bulk RNA-seq and scRNA-seq results with improvements in overall dropout rates and sensitivity. Finally, we found that MERFISH independently resolved distinct cell types and spatial structure in both the liver and kidney. Computational integration with the Tabula Muris Senis atlas did not enhance these results. We conclude that MERFISH provides a quantitatively comparable method for single-cell gene expression and can identify cell types without the need for computational integration with scRNA-seq atlases.

Cohesin and CTCF control the dynamics of chromosome folding

In mammals, interactions between sequences within topologically associating domains enable control of gene expression across large genomic distances. Yet it is unknown how frequently such contacts occur, how long they last and how they depend on the dynamics of chromosome folding and loop extrusion activity of cohesin. By imaging chromosomal locations at high spatial and temporal resolution in living cells, we show that interactions within topologically associating domains are transient and occur frequently during the course of a cell cycle. Interactions become more frequent and longer in the presence of convergent CTCF sites, resulting in suppression of variability in chromosome folding across time. Supported by physical models of chromosome dynamics, our data suggest that CTCF-anchored loops last around 10 min. Our results show that long-range transcriptional regulation might rely on transient physical proximity, and that cohesin and CTCF stabilize highly dynamic chromosome structures, facilitating selected subsets of chromosomal interactions.

Dissecting the cosegregation probability from genome architecture mapping

Genome architecture mapping (GAM) is a recently developed methodology that offers the cosegregation probability of two genomic segments from an ensemble of thinly sliced nuclear profiles, enabling us to probe and decipher three-dimensional chromatin organization. The cosegregation probability from GAM binned at 1 Mb, which thus probes the length scale associated with the genomic separation greater than 1 Mb, is, however, not identical to the contact probability obtained from Hi-C, and its correlation with interlocus distance measured with fluorescence in situ hybridization is not so good as the contact probability. In this study, by using a polymer-based model of chromatins, we derive a theoretical expression of the cosegregation probability as well as that of the contact probability and carry out quantitative analyses of how they differ from each other. The results from our study, validated with in silico GAM analysis on three-dimensional genome structures from fluorescence in situ hybridization, suggest that to attain strong correlation with the interlocus distance, a properly normalized version of cosegregation probability needs to be calculated based on a large number of nuclear slices (n>103).

The emerging landscape of spatial profiling technologies

Improved scale, multiplexing and resolution are establishing spatial nucleic acid and protein profiling methods as a major pillar for cellular atlas building of complex samples, from tissues to full organisms. Emerging methods yield omics measurements at resolutions covering the nano- to microscale, enabling the charting of cellular heterogeneity, complex tissue architectures and dynamic changes during development and disease. We present an overview of the developing landscape of in situ spatial genome, transcriptome and proteome technologies, exemplify their impact on cell biology and translational research, and discuss current challenges for their community-wide adoption. Among many transformative applications, we envision that spatial methods will map entire organs and enable next-generation pathology.

Nuclear speckles - a driving force in gene expression

Nuclear speckles are dynamic membraneless bodies located in the cell nucleus. They harbor RNAs and proteins, many of which are splicing factors, that together display complex biophysical properties dictating nuclear speckle formation and maintenance. Although these nuclear bodies were discovered decades ago, only recently has in-depth genomic analysis begun to unravel their essential functions in modulation of gene activity. Major advancements in genomic mapping techniques combined with microscopy approaches have enabled insights into the roles nuclear speckles may play in enhancing gene expression, and how gene positioning to specific nuclear landmarks can regulate gene expression and RNA processing. Some studies have drawn a link between nuclear speckles and disease. Certain maladies either involve nuclear speckles directly or dictate the localization and reorganization of many nuclear speckle factors. This is most striking during viral infection, as viruses alter the entire nuclear architecture and highjack host machinery. As discussed in this Review, nuclear speckles represent a fascinating target of study not only to reveal the links between gene positioning, genome subcompartments and gene activity, but also as a potential target for therapeutics.

Extracting physical characteristics of higher-order chromatin structures from 3D image data

Higher-order chromatin structures have functional impacts on gene regulation and cell identity determination. Using high-throughput sequencing (HTS)-based methods like Hi-C, active or inactive compartments and open or closed topologically associating domain (TAD) structures can be identified on a cell population level. Recently developed high-resolution three-dimensional (3D) molecular imaging techniques such as 3D electron microscopy with in situ hybridization (3D-EMSIH) and 3D structured illumination microscopy (3D-SIM) enable direct detection of physical representations of chromatin structures in a single cell. However, computational analysis of 3D image data with explainability and interpretability on functional characteristics of chromatin structures is still challenging. We developed Extracting Physical-Characteristics from Images of Chromatin Structures (EPICS), a machine-learning based computational method for processing high-resolution chromatin 3D image data. Using EPICS on images produced by 3D-EMISH or 3D-SIM techniques, we generated more direct 3D representations of higher-order chromatin structures, identified major chromatin domains, and determined the open or closed status of each domain. We identified several high-contributing features from the model as the major physical characteristics that define the open or closed chromatin domains, demonstrating the explainability and interpretability of EPICS. EPICS can be applied to the analysis of other high-resolution 3D molecular imaging data for spatial genomics studies. The R and Python codes of EPICS are available at https://github.com/zang-lab/epics.

The era of 3D and spatial genomics

Over a decade ago the advent of high-throughput chromosome conformation capture (Hi-C) sparked a new era of 3D genomics. Since then the number of methods for mapping the 3D genome has flourished, enabling an ever-increasing understanding of how DNA is packaged in the nucleus and how the spatiotemporal organization of the genome orchestrates its vital functions. More recently, the next generation of spatial genomics technologies has begun to reveal how genome sequence and 3D genome organization vary between cells in their tissue context. We summarize how the toolkit for charting genome topology has evolved over the past decade and discuss how new technological developments are advancing the field of 3D and spatial genomics.

The Impact of Space and Time on the Functional Output of the Genome

Over the past two decades, it has become clear that the multiscale spatial and temporal organization of the genome has important implications for nuclear function. This review centers on insights gained from recent advances in light microscopy on our understanding of transcription. We discuss spatial and temporal aspects that shape nuclear order and their consequences on regulatory components, focusing on genomic scales most relevant to function. The emerging picture is that spatiotemporal constraints increase the complexity in transcriptional regulation, highlighting new challenges, such as uncertainty about how information travels from molecular factors through the genome and space to generate a functional output.

Reconstruct high-resolution 3D genome structures for diverse cell-types using FLAMINGO

High-resolution reconstruction of spatial chromosome organizations from chromatin contact maps is highly demanded, but is hindered by extensive pairwise constraints, substantial missing data, and limited resolution and cell-type availabilities. Here, we present FLAMINGO, a computational method that addresses these challenges by compressing inter-dependent Hi-C interactions to delineate the underlying low-rank structures in 3D space, based on the low-rank matrix completion technique. FLAMINGO successfully generates 5 kb- and 1 kb-resolution spatial conformations for all chromosomes in the human genome across multiple cell-types, the largest resources to date. Compared to other methods using various experimental metrics, FLAMINGO consistently demonstrates superior accuracy in recapitulating observed structures with raises in scalability by orders of magnitude. The reconstructed 3D structures efficiently facilitate discoveries of higher-order multi-way interactions, imply biological interpretations of long-range QTLs, reveal geometrical properties of chromatin, and provide high-resolution references to understand structural variabilities. Importantly, FLAMINGO achieves robust predictions against high rates of missing data and significantly boosts 3D structure resolutions. Moreover, FLAMINGO shows vigorous cross cell-type structure predictions that capture cell-type specific spatial configurations via integration of 1D epigenomic signals. FLAMINGO can be widely applied to large-scale chromatin contact maps and expand high-resolution spatial genome conformations for diverse cell-types.

Epigenetic-structural changes in X chromosomes promote Xic pairing during early differentiation of mouse embryonic stem cells

X chromosome inactivation center (Xic) pairing occurs during the differentiation of embryonic stem (ES) cells from female mouse embryos, and is related to X chromosome inactivation, the circadian clock, intra-nucleus architecture, and metabolism. However, the mechanisms underlying the identification and approach of X chromosome pairs in the crowded nucleus are unclear. To elucidate the driving force of Xic pairing, we developed a coarse-grained molecular dynamics model of intranuclear chromosomes in ES cells and in cells 2 days after the onset of differentiation (2-day cells) by considering intrachromosomal epigenetic-structural feature-dependent mechanics. The analysis of the experimental data showed that X-chromosomes exhibit the rearrangement of their distributions of open/closed chromatin regions on their surfaces during cell differentiation. By simulating models where the excluded volume effects of closed chromatin regions are stronger than those of open chromatin regions, such rearrangement of open/closed chromatin regions on X-chromosome surfaces promoted the mutual approach of the Xic pair. These findings suggested that local intrachromosomal epigenetic features may contribute to the regulation of cell species-dependent differences in intranuclear architecture.

A Shift in Paradigms: Spatial Genomics Approaches to Reveal Single-Cell Principles of Genome Organization

The genome tridimensional (3D) organization and its role towards the regulation of key cell processes such as transcription is currently a main question in biology. Interphase chromosomes are spatially segregated into "territories," epigenetically-defined large domains of chromatin that interact to form "compartments" with common transcriptional status, and insulator-flanked domains called "topologically associating domains" (TADs). Moreover, chromatin organizes around nuclear structures such as lamina, speckles, or the nucleolus to acquire a higher-order genome organization. Due to recent technological advances, the different hierarchies are being solved. Particularly, advances in microscopy technologies are shedding light on the genome structure at multiple levels. Intriguingly, more and more reports point to high variability and stochasticity at the single-cell level. However, the functional consequences of such variability in genome conformation are still unsolved. Here, I will discuss the implication of the cell-to-cell heterogeneity at the different scales in the context of newly developed imaging approaches, particularly multiplexed Fluorescence in situ hybridization methods that enabled "chromatin tracing." Extensions of these methods are now combining spatial information of dozens to thousands of genomic loci with the localization of nuclear features such as the nucleolus, nuclear speckles, or even histone modifications, creating the fast-moving field of "spatial genomics." As our view of genome organization shifts the focus from ensemble to single-cell, new insights to fundamental questions begin to emerge.

TAD-like single-cell domain structures exist on both active and inactive X chromosomes and persist under epigenetic perturbations

BACKGROUND

Topologically associating domains (TADs) are important building blocks of three-dimensional genome architectures. The formation of TADs has been shown to depend on cohesin in a loop-extrusion mechanism. Recently, advances in an image-based spatial genomics technique known as chromatin tracing lead to the discovery of cohesin-independent TAD-like structures, also known as single-cell domains, which are highly variant self-interacting chromatin domains with boundaries that occasionally overlap with TAD boundaries but tend to differ among single cells and among single chromosome copies. Recent computational modeling studies suggest that epigenetic interactions may underlie the formation of the single-cell domains.

RESULTS

Here we use chromatin tracing to visualize in female human cells the fine-scale chromatin folding of inactive and active X chromosomes, which are known to have distinct global epigenetic landscapes and distinct population-averaged TAD profiles, with inactive X chromosomes largely devoid of TADs and cohesin. We show that both inactive and active X chromosomes possess highly variant single-cell domains across the same genomic region despite the fact that only active X chromosomes show clear TAD structures at the population level. These X chromosome single-cell domains exist in distinct cell lines. Perturbations of major epigenetic components and transcription mostly do not affect the frequency or strength of the single-cell domains. Increased chromatin compaction of inactive X chromosomes occurs at a length scale above that of the single-cell domains.

CONCLUSIONS

In sum, this study suggests that single-cell domains are genome architecture building blocks independent of the tested major epigenetic components.

FISH Going Meso-Scale: A Microscopic Search for Chromatin Domains

The intimate relationships between genome structure and function direct efforts toward deciphering three-dimensional chromatin organization within the interphase nuclei at different genomic length scales. For decades, major insights into chromatin structure at the level of large-scale euchromatin and heterochromatin compartments, chromosome territories, and subchromosomal regions resulted from the evolution of light microscopy and fluorescence in situ hybridization. Studies of nanoscale nucleosomal chromatin organization benefited from a variety of electron microscopy techniques. Recent breakthroughs in the investigation of mesoscale chromatin structures have emerged from chromatin conformation capture methods (C-methods). Chromatin has been found to form hierarchical domains with high frequency of local interactions from loop domains to topologically associating domains and compartments. During the last decade, advances in super-resolution light microscopy made these levels of chromatin folding amenable for microscopic examination. Here we are reviewing recent developments in FISH-based approaches for detection, quantitative measurements, and validation of contact chromatin domains deduced from C-based data. We specifically focus on the design and application of Oligopaint probes, which marked the latest progress in the imaging of chromatin domains. Vivid examples of chromatin domain FISH-visualization by means of conventional, super-resolution light and electron microscopy in different model organisms are provided.

Single-cell nuclear architecture across cell types in the mouse brain

Diverse cell types in tissues have distinct gene expression programs, chromatin states, and nuclear architectures. To correlate such multimodal information across thousands of single cells in mouse brain tissue sections, we use integrated spatial genomics, imaging thousands of genomic loci along with RNAs and epigenetic markers simultaneously in individual cells. We reveal that cell type–specific association and scaffolding of DNA loci around nuclear bodies organize the nuclear architecture and correlate with differential expression levels in different cell types. At the submegabase level, active and inactive X chromosomes access similar domain structures in single cells despite distinct epigenetic and expression states. This work represents a major step forward in linking single-cell three-dimensional nuclear architecture, gene expression, and epigenetic modifications in a native tissue context.

Spatial transcriptome profiling by MERFISH reveals fetal liver hematopoietic stem cell niche architecture

The hematopoietic stem cell (HSC) niche has been extensively studied in bone marrow, yet a more systematic investigation into the microenvironment regulation of hematopoiesis in fetal liver is necessary. Here we investigate the spatial organization and transcriptional profile of individual cells in both wild type (WT) and Tet2^−/− fetal livers, by multiplexed error robust fluorescence in situ hybridization. We find that specific pairs of fetal liver cell types are preferentially positioned next to each other. Ligand-receptor signaling molecule pairs such as Kitl and Kit are enriched in neighboring cell types. The majority of HSCs are in direct contact with endothelial cells (ECs) in both WT and Tet2^−/− fetal livers. Loss of Tet2 increases the number of HSCs, and upregulates Wnt and Notch signaling genes in the HSC niche. Two subtypes of ECs, arterial ECs and sinusoidal ECs, and other cell types contribute distinct signaling molecules to the HSC niche. Collectively, this study provides a comprehensive picture and bioinformatic foundation for HSC spatial regulation in fetal liver.

Deep learning connects DNA traces to transcription to reveal predictive features beyond enhancer-promoter contact

Chromatin architecture plays an important role in gene regulation. Recent advances in super-resolution microscopy have made it possible to measure chromatin 3D structure and transcription in thousands of single cells. However, leveraging these complex data sets with a computationally unbiased method has been challenging. Here, we present a deep learning-based approach to better understand to what degree chromatin structure relates to transcriptional state of individual cells. Furthermore, we explore methods to "unpack the black box" to determine in an unbiased manner which structural features of chromatin regulation are most important for gene expression state. We apply this approach to an Optical Reconstruction of Chromatin Architecture dataset of the Bithorax gene cluster in Drosophila and show it outperforms previous contact-focused methods in predicting expression state from 3D structure. We find the structural information is distributed across the domain, overlapping and extending beyond domains identified by prior genetic analyses. Individual enhancer-promoter interactions are a minor contributor to predictions of activity.

Advances in higher-order chromatin architecture: the move towards 4D genome

In eukaryotes, the genome is hierarchically packed inside the nucleus, which facilitates physical contact between cis-regulatory elements (CREs), such as enhancers and promoters. Accumulating evidence highlights the critical role of higher-order chromatin structure in precise regulation of spatiotemporal gene expression under diverse biological contexts including lineage commitment and cell activation by external stimulus. Genomics and imaging-based technologies, such as Hi-C and DNA fluorescence in situ hybridization (FISH), have revealed the key principles of genome folding, while newly developed tools focus on improvement in resolution, throughput and modality at single-cell and population levels, and challenge the knowledge obtained through conventional approaches. In this review, we discuss recent advances in our understanding of principles of higher-order chromosome conformation and technologies to investigate 4D chromatin interactions.

The two waves in single-cell 3D genomics

For decades, biochemical methods for the analysis of genome structure and function provided cell-population-averaged data that allowed general principles and tendencies to be disclosed. Microscopy-based studies, which immanently involve single-cell analysis, did not provide sufficient spatial resolution to investigate the particularly small details of 3D genome folding. Nevertheless, these studies demonstrated that mutual positions of chromosome territories within cell nuclei and individual genomic loci within chromosomal territories can vary significantly in individual cells. The development of new technologies in biochemistry and the advent of super-resolution microscopy in the last decade have made possible the full-scale study of 3D genome organization in individual cells. Maps of the 3D genome build based on C-data and super-resolution microscopy are highly consistent and, therefore, biologically relevant. The internal structures of individual chromosomes, loci, and topologically associating domains (TADs) are resolved as well as cell-cycle dynamics. 3D modeling allows one to investigate the physical mechanisms underlying genome folding. Finally, joint profiling of genome topology and epigenetic features will allow 3D genomics to handle complex cell-to-cell heterogeneity. In this review, we summarize the present state of studies into 3D genome organization in individual cells, analyze the technical problems of single-cell studies, and outline perspectives of 3D genomics.

Chromatin tracing and multiplexed imaging of nucleome architectures (MINA) and RNAs in single mammalian cells and tissue

The genome is hierarchically organized into several 3D architectures, including chromatin loops, domains, compartments and regions associated with nuclear lamina and nucleoli. Changes in these architectures have been associated with normal development, aging and a wide range of diseases. Despite its critical importance, understanding how the genome is spatially organized in single cells, how organization varies in different cell types in mammalian tissue and how organization affects gene expression remains a major challenge. Previous approaches have been limited by a lack of capacity to directly trace chromatin folding in 3D and to simultaneously measure genomic organization in relation to other nuclear components and gene expression in the same single cells. We have developed an image-based 3D genomics technique termed 'chromatin tracing', which enables direct 3D tracing of chromatin folding along individual chromosomes in single cells. More recently, we also developed multiplexed imaging of nucleome architectures (MINA), which enables simultaneous measurements of multiscale chromatin folding, associations of genomic regions with nuclear lamina and nucleoli and copy numbers of numerous RNA species in the same single cells in mammalian tissue. Here, we provide detailed protocols for chromatin tracing in cell lines and MINA in mammalian tissue, which take 3-4 d for experimental work and 2-3 d for data analysis. We expect these developments to be broadly applicable and to affect many lines of research on 3D genomics by depicting multiscale genomic architectures associated with gene expression, in different types of cells and tissue undergoing different biological processes.

Interactions Between Nucleosomes: From Atomistic Simulation to Polymer Model

The organization of genomes in space and time dimension plays an important role in gene expression and regulation. Chromatin folding occurs in a dynamic, structured way that is subject to biophysical rules and biological processes. Nucleosomes are the basic unit of chromatin in living cells, and here we report on the effective interactions between two nucleosomes in physiological conditions using explicit-solvent all-atom simulations. Free energy landscapes derived from umbrella sampling simulations agree well with recent experimental and simulation results. Our simulations reveal the atomistic details of the interactions between nucleosomes in solution and can be used for constructing the coarse-grained model for chromatin in a bottom-up manner.