Sci-Hi-C: A single-cell Hi-C method for mapping 3D genome organization in large number of single cells

Insights from the first chromosome-level genome assembly of the alpine gentian Gentiana straminea Maxim

Gentiana straminea Maxim. is a perennial herb and mainly distributed in the Qinghai-Tibetan Plateau. To adapt to the extreme environment, it has developed particular morphological, physiological, and genetic structures. Also, rich in iridoids, it is one of the original plants of traditional Chinese herb 'Qinjiao'. Herein, we present its first chromosome-level genome sequence assembly and compare it with the genomes of other Gentiana species to facilitate the analysis of genomic characteristics. The assembled genome size of G. straminea was 1.25 Gb, with a contig N50 of 7.5 Mb. A total of 96.08% of the genome sequences was anchored on 13 pseudochromosomes, with a scaffold N50 of 92.70 Mb. A total of 54,310 protein-coding genes were predicted, 80.25% of which were functionally annotated. Comparative genomic analyses indicated that G. straminea experienced two whole-genome duplication events after the γ whole-genome triplication with other eudicots, and it diverged from other Gentiana species at ~3.2 Mya. A total of 142 enzyme-coding genes related to iridoid biosynthesis were identified in its genome. Additionally, we identified differences in the number and expression patterns of iridoid biosynthetic pathway genes in G. straminea compared with two other Gentiana species by integrating whole-genome sequence and transcriptomic analyses.

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.

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.

GAGE-seq concurrently profiles multiscale 3D genome organization and gene expression in single cells

The organization of mammalian genomes features a complex, multiscale three-dimensional (3D) architecture, whose functional significance remains elusive because of limited single-cell technologies that can concurrently profile genome organization and transcriptional activities. Here, we introduce genome architecture and gene expression by sequencing (GAGE-seq), a scalable, robust single-cell co-assay measuring 3D genome structure and transcriptome simultaneously within the same cell. Applied to mouse brain cortex and human bone marrow CD34+ cells, GAGE-seq characterized the intricate relationships between 3D genome and gene expression, showing that multiscale 3D genome features inform cell-type-specific gene expression and link regulatory elements to target genes. Integration with spatial transcriptomic data revealed in situ 3D genome variations in mouse cortex. Observations in human hematopoiesis unveiled discordant changes between 3D genome organization and gene expression, underscoring a complex, temporal interplay at the single-cell level. GAGE-seq provides a powerful, cost-effective approach for exploring genome structure and gene expression relationships at the single-cell level across diverse biological contexts.

SnapHiC-G: identifying long-range enhancer-promoter interactions from single-cell Hi-C data via a global background model

Harnessing the power of single-cell genomics technologies, single-cell Hi-C (scHi-C) and its derived technologies provide powerful tools to measure spatial proximity between regulatory elements and their target genes in individual cells. Using a global background model, we propose SnapHiC-G, a computational method, to identify long-range enhancer-promoter interactions from scHi-C data. We applied SnapHiC-G to scHi-C datasets generated from mouse embryonic stem cells and human brain cortical cells. SnapHiC-G achieved high sensitivity in identifying long-range enhancer-promoter interactions. Moreover, SnapHiC-G can identify putative target genes for noncoding genome-wide association study (GWAS) variants, and the genetic heritability of neuropsychiatric diseases is enriched for single-nucleotide polymorphisms (SNPs) within SnapHiC-G-identified interactions in a cell-type-specific manner. In sum, SnapHiC-G is a powerful tool for characterizing cell-type-specific enhancer-promoter interactions from complex tissues and can facilitate the discovery of chromatin interactions important for gene regulation in biologically relevant cell types.

Regulating epithelial-mesenchymal plasticity from 3D genome organization

Epithelial-mesenchymal transition (EMT) is a dynamic process enabling polarized epithelial cells to acquire mesenchymal features implicated in development and carcinoma progression. As our understanding evolves, it is clear the reversible execution of EMT arises from complex epigenomic regulation involving histone modifications and 3-dimensional (3D) genome structural changes, leading to a cascade of transcriptional events. This review summarizes current knowledge on chromatin organization in EMT, with a focus on hierarchical structures of the 3D genome and chromatin accessibility changes.

Decoding polygenic diseases: advances in noncoding variant prioritization and validation

Genome-wide association studies (GWASs) provide a key foundation for elucidating the genetic underpinnings of common polygenic diseases. However, these studies have limitations in their ability to assign causality to particular genetic variants, especially those residing in the noncoding genome. Over the past decade, technological and methodological advances in both analytical and empirical prioritization of noncoding variants have enabled the identification of causative variants by leveraging orthogonal functional evidence at increasing scale. In this review, we present an overview of these approaches and describe how this workflow provides the groundwork necessary to move beyond associations toward genetically informed studies on the molecular and cellular mechanisms of polygenic disease.

Cepharanthine analogs mining and genomes of Stephania accelerate anti-coronavirus drug discovery

Cepharanthine is a secondary metabolite isolated from Stephania. It has been reported that it has anti-conronaviruses activities including severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). Here, we assemble three Stephania genomes (S. japonica, S. yunnanensis, and S. cepharantha), propose the cepharanthine biosynthetic pathway, and assess the antiviral potential of compounds involved in the pathway. Among the three genomes, S. japonica has a near telomere-to-telomere assembly with one remaining gap, and S. cepharantha and S. yunnanensis have chromosome-level assemblies. Following by biosynthetic gene mining and metabolomics analysis, we identify seven cepharanthine analogs that have broad-spectrum anti-coronavirus activities, including SARS-CoV-2, Guangxi pangolin-CoV (GX_P2V), swine acute diarrhoea syndrome coronavirus (SADS-CoV), and porcine epidemic diarrhea virus (PEDV). We also show that two other genera, Nelumbo and Thalictrum, can produce cepharanthine analogs, and thus have the potential for antiviral compound discovery. Results generated from this study could accelerate broad-spectrum anti-coronavirus drug discovery.

Next-generation forward genetic screens: uniting high-throughput perturbations with single-cell analysis

Programmable genome-engineering technologies, such as CRISPR (clustered regularly interspaced short palindromic repeats) nucleases and massively parallel CRISPR screens that capitalize on this programmability, have transformed biomedical science. These screens connect genes and noncoding genome elements to disease-relevant phenotypes, but until recently have been limited to individual phenotypes such as growth or fluorescent reporters of gene expression. By pairing massively parallel screens with high-dimensional profiling of single-cell types/states, we can now measure how individual genetic perturbations or combinations of perturbations impact the cellular transcriptome, proteome, and epigenome. We review technologies that pair CRISPR screens with single-cell multiomics and the unique opportunities afforded by extending pooled screens using deep multimodal phenotyping.

The Metasequoia genome and evolutionary relationships among redwoods

Redwood trees (Sequoioideae), including Metasequoia glyptostroboides (dawn redwood), Sequoiadendron giganteum (giant sequoia), and Sequoia sempervirens (coast redwood), are threatened and widely recognized iconic tree species. Genomic resources for redwood trees could provide clues to their evolutionary relationships. Here, we report the 8-Gb reference genome of M. glyptostroboides and a comparative analysis with two related species. More than 62% of the M. glyptostroboides genome is composed of repetitive sequences. Clade-specific bursts of long terminal repeat retrotransposons may have contributed to genomic differentiation in the three species. The chromosomal synteny between M. glyptostroboides and S. giganteum is extremely high, whereas there has been significant chromosome reorganization in S. sempervirens. Phylogenetic analysis of marker genes indicates that S. sempervirens is an autopolyploid, and more than 48% of the gene trees are incongruent with the species tree. Results of multiple analyses suggest that incomplete lineage sorting (ILS) rather than hybridization explains the inconsistent phylogeny, indicating that genetic variation among redwoods may be due to random retention of polymorphisms in ancestral populations. Functional analysis of ortholog groups indicates that gene families of ion channels, tannin biosynthesis enzymes, and transcription factors for meristem maintenance have expanded in S. giganteum and S. sempervirens, which is consistent with their extreme height. As a wetland-tolerant species, M. glyptostroboides shows a transcriptional response to flooding stress that is conserved with that of analyzed angiosperm species. Our study offers insights into redwood evolution and adaptation and provides genomic resources to aid in their conservation and management.

Concurrent profiling of multiscale 3D genome organization and gene expression in single mammalian cells

The organization of mammalian genomes within the nucleus features a complex, multiscale three-dimensional (3D) architecture. The functional significance of these 3D genome features, however, remains largely elusive due to limited single-cell technologies that can concurrently profile genome organization and transcriptional activities. Here, we report GAGE-seq, a highly scalable, robust single-cell co-assay that simultaneously measures 3D genome structure and transcriptome within the same cell. Employing GAGE-seq on mouse brain cortex and human bone marrow CD34+ cells, we comprehensively characterized the intricate relationships between 3D genome and gene expression. We found that these multiscale 3D genome features collectively inform cell type-specific gene expressions, hence contributing to defining cell identity at the single-cell level. Integration of GAGE-seq data with spatial transcriptomic data revealed in situ variations of the 3D genome in mouse cortex. Moreover, our observations of lineage commitment in normal human hematopoiesis unveiled notable discordant changes between 3D genome organization and gene expression, underscoring a complex, temporal interplay at the single-cell level that is more nuanced than previously appreciated. Together, GAGE-seq provides a powerful, cost-effective approach for interrogating genome structure and gene expression relationships at the single-cell level across diverse biological contexts.

New discoveries in the field of metabolism by applying single-cell and spatial omics

Single-cell multi-Omics (SCM-Omics) and spatial multi-Omics (SM-Omics) technologies provide state-of-the-art methods for exploring the composition and function of cell types in tissues/organs. Since its emergence in 2009, single-cell RNA sequencing (scRNA-seq) has yielded many groundbreaking new discoveries. The combination of this method with the emergence and development of SM-Omics techniques has been a pioneering strategy in neuroscience, developmental biology, and cancer research, especially for assessing tumor heterogeneity and T-cell infiltration. In recent years, the application of these methods in the study of metabolic diseases has also increased. The emerging SCM-Omics and SM-Omics approaches allow the molecular and spatial analysis of cells to explore regulatory states and determine cell fate, and thus provide promising tools for unraveling heterogeneous metabolic processes and making them amenable to intervention. Here, we review the evolution of SCM-Omics and SM-Omics technologies, and describe the progress in the application of SCM-Omics and SM-Omics in metabolism-related diseases, including obesity, diabetes, nonalcoholic fatty liver disease (NAFLD) and cardiovascular disease (CVD). We also conclude that the application of SCM-Omics and SM-Omics approaches can help resolve the molecular mechanisms underlying the pathogenesis of metabolic diseases in the body and facilitate therapeutic measures for metabolism-related diseases. This review concludes with an overview of the current status of this emerging field and the outlook for its future.

scHi-CSim: a flexible simulator that generates high-fidelity single-cell Hi-C data for benchmarking

Single-cell Hi-C technology provides an unprecedented opportunity to reveal chromatin structure in individual cells. However, high sequencing cost impedes the generation of biological Hi-C data with high sequencing depths and multiple replicates for downstream analysis. Here, we developed a single-cell Hi-C simulator (scHi-CSim) that generates high-fidelity data for benchmarking. scHi-CSim merges neighboring cells to overcome the sparseness of data, samples interactions in distance-stratified chromosomes to maintain the heterogeneity of single cells, and estimates the empirical distribution of restriction fragments to generate simulated data. We demonstrated that scHi-CSim can generate high-fidelity data by comparing the performance of single-cell clustering and detection of chromosomal high-order structures with raw data. Furthermore, scHi-CSim is flexible to change sequencing depth and the number of simulated replicates. We showed that increasing sequencing depth could improve the accuracy of detecting topologically associating domains. We also used scHi-CSim to generate a series of simulated datasets with different sequencing depths to benchmark scHi-C clustering methods.

Elucidating the structure and function of the nucleus-The NIH Common Fund 4D Nucleome program

Genomic architecture appears to play crucial roles in health and a variety of diseases. How nuclear structures reorganize over different timescales is elusive, partly because the tools needed to probe and perturb them are not as advanced as needed by the field. To fill this gap, the National Institutes of Health Common Fund started a program in 2015, called the 4D Nucleome (4DN), with the goal of developing and ultimately applying technologies to interrogate the structure and function of nuclear organization in space and time.

3D genome alterations in T cells associated with disease activity of systemic lupus erythematosus

OBJECTIVES

Three-dimensional (3D) genome alterations can dysregulate gene expression by rewiring physical interactions within chromosomes in a tissue-specific or cell-specific manner and lead to diseases. We aimed to elucidate the 3D genome structure and its role in gene expression networks dysregulated in systemic lupus erythematosus (SLE).

METHODS

We performed Hi-C experiments using CD4+ T cells from 7 patients with SLE and 5 age-matched and sex-matched healthy controls (HCs) combined with RNA sequencing analysis. Further integrative analyses, including transcription factor motif enrichment, SPI1 knockdown and histone modifications (H3K27ac, H3K4me1, H3K4me3), were performed for altered loop-associated gene loci in SLE.

RESULTS

We deciphered the 3D chromosome organisation in T cells of patients with SLE and found it was clearly distinct from that of HCs and closely associated with the disease activity of SLE. Importantly, we identified loops within chromosomes associated with the disease activity of SLE and differentially expressed genes and found some key histone modifications close to these loops. Moreover, we demonstrated the contribution of the transcription factor SPI1, whose motif is located in the altered loop in SLE, to the overexpression of interferon pathway gene. In addition, we identified the potential influences of genetic variations in 3D genome alterations in SLE.

CONCLUSIONS

Our results highlight the 3D genome structure alterations associated with SLE development and provide a foundation for future interrogation of the relationships between chromosome structure and gene expression control in SLE.

Characterizing cis-regulatory elements using single-cell epigenomics

Cell type-specific gene expression patterns and dynamics during development or in disease are controlled by cis-regulatory elements (CREs), such as promoters and enhancers. Distinct classes of CREs can be characterized by their epigenomic features, including DNA methylation, chromatin accessibility, combinations of histone modifications and conformation of local chromatin. Tremendous progress has been made in cataloguing CREs in the human genome using bulk transcriptomic and epigenomic methods. However, single-cell epigenomic and multi-omic technologies have the potential to provide deeper insight into cell type-specific gene regulatory programmes as well as into how they change during development, in response to environmental cues and through disease pathogenesis. Here, we highlight recent advances in single-cell epigenomic methods and analytical tools and discuss their readiness for human tissue profiling.

Disruption of the gene regulatory programme in neurodevelopmental disorders

Cortical development consists of a series of synchronised events, including fate transition of cortical progenitors, neuronal migration, specification and connectivity. It is becoming clear that gene expression programs governing these events rely on the interplay between signalling molecules, transcription factors and epigenetic mechanisms. When genetic or environmental factors disrupt expression of genes involved in important brain development processes, neurodevelopmental disorders can occur. This review aims to highlight how recent advances in technologies have helped uncover and imitate the gene regulatory mechanisms commonly disrupted in neurodevelopmental disorders.

A guide to systems-level immunomics

The immune system is highly complex and distributed throughout an organism, with hundreds to thousands of cell states existing in parallel with diverse molecular pathways interacting in a highly dynamic and coordinated fashion. Although the characterization of individual genes and molecules is of the utmost importance for understanding immune-system function, high-throughput, high-resolution omics technologies combined with sophisticated computational modeling and machine-learning approaches are creating opportunities to complement standard immunological methods with new insights into immune-system dynamics. Like systems immunology itself, immunology researchers must take advantage of these technologies and form their own diverse networks, connecting with researchers from other disciplines. This Review is an introduction and 'how-to guide' for immunologists with no particular experience in the field of omics but with the intention to learn about and apply these systems-level approaches, and for immunologists who want to make the most of interdisciplinary networks.

Single cell- and spatial 'Omics revolutionize physiology

Single cell multi- 'Omics and Spatial Transcriptomics are prominent technological highlights of recent years, and both fields still witness a ceaseless firework of novel approaches for high resolution profiling of additional omics layers. As all life processes in organs and organisms are based on the functions of their fundamental building blocks, the individual cells and their interactions, these methods are of utmost worth for the study of physiology in health and disease. Recent discoveries on embryonic development, tumor immunology, the detailed cellular composition and function of complex tissues like for example the kidney or the brain, different roles of the same cell type in different organs, the oncogenic program of individual tumor entities, or the architecture of immunopathology in infected tissue are based on single cell and spatial transcriptomics experiments. In this review, we will give a broad overview of technological concepts for single cell and spatial analysis, showing both advantages and limitations, and illustrate their impact with some particularly impressive case studies.

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.

Single-Cell Epigenomics Reveals Mechanisms of Cancer Progression

Cancer initiation is driven by the cooperation between genetic and epigenetic aberrations that disrupt gene regulatory programs critical to maintaining specialized cellular functions. After initiation, cells acquire additional genetic and epigenetic alterations influenced by tumor-intrinsic and -extrinsic mechanisms, which increase intratumoral heterogeneity, reshape the cell's underlying gene regulatory networks and promote cancer evolution. Furthermore, environmental or therapeutic insults drive the selection of heterogeneous cell states, with implications for cancer initiation, maintenance, and drug resistance. The advancement of single-cell genomics has begun to uncover the full repertoire of chromatin and gene expression states (cell states) that exist within individual tumors. These single-cell analyses suggest that cells diversify in their regulatory states upon transformation by co-opting damage-induced and nonlineage regulatory programs that can lead to epigenomic plasticity. Here, we review these recent studies related to regulatory state changes in cancer progression and highlight the growing single-cell epigenomics toolkit poised to address unresolved questions in the field.

Cell-Type Resolved Insights into the Cis-Regulatory Genome of NAFLD

The prevalence of non-alcoholic fatty liver disease (NAFLD) is increasing rapidly, and unmet treatment can result in the development of hepatitis, fibrosis, and liver failure. There are difficulties involved in diagnosing NAFLD early and for this reason there are challenges involved in its treatment. Furthermore, no drugs are currently approved to alleviate complications, a fact which highlights the need for further insight into disease mechanisms. NAFLD pathogenesis is associated with complex cellular changes, including hepatocyte steatosis, immune cell infiltration, endothelial dysfunction, hepatic stellate cell activation, and epithelial ductular reaction. Many of these cellular changes are controlled by dramatic changes in gene expression orchestrated by the cis-regulatory genome and associated transcription factors. Thus, to understand disease mechanisms, we need extensive insights into the gene regulatory mechanisms associated with tissue remodeling. Mapping cis-regulatory regions genome-wide is a step towards this objective and several current and emerging technologies allow detection of accessible chromatin and specific histone modifications in enriched cell populations of the liver, as well as in single cells. Here, we discuss recent insights into the cis-regulatory genome in NAFLD both at the organ-level and in specific cell populations of the liver. Moreover, we highlight emerging technologies that enable single-cell resolved analysis of the cis-regulatory genome of the liver.

HiCuT: An efficient and low input method to identify protein-directed chromatin interactions

3D genome organization regulates gene expression, and disruption of these long-range (>20kB) DNA-protein interactions results in pathogenic phenotypes. Chromosome conformation methods in conjunction with chromatin immunoprecipitation were used to decipher protein-directed chromatin interactions. However, these methods required abundant starting material (>500,000 cells), sizable number of sequencing reads (>100 million reads), and elaborate data processing methods to reduce background noise, which limited their use in primary cells. Hi-C Coupled chromatin cleavage and Tagmentation (HiCuT) is a new transposase-assisted tagmentation method that generates high-resolution protein directed long-range chromatin interactions as efficiently as existing methods, HiChIP and ChIA-PET, despite using 100,000 cells (5-fold less) and 12 million sequencing reads (8-fold fewer). Moreover, HiCuT generates high resolution fragment libraries with low background signal that are easily interpreted with minimal computational processing. We used HiCuT in human primary skin cells to link previously identified single nucleotide polymorphisms (SNPs) in skin disease to candidate genes and to identify functionally relevant transcription factors in an unbiased manner. HiCuT broadens the capacity for genomic profiling in systems previously unmeasurable, including primary cells, human tissue samples, and rare cell populations, and may be a useful tool for all investigators studying human genetics and personalized epigenomics.

DNA is precisely packaged and organized within a nucleus to regulate gene expression. Altering this structure results in disease and developmental abnormalities. Current methods to probe 3D genome organization require a minimum of 500,000 cells, 100 million sequencing reads, and elaborate computational genomics skills. This limits general adoption and prevents use in small populations of cells, including primary tissues. We report a new method called Hi-C Coupled chromatin cleavage and Tagmentation (HiCuT). This method couples current 3D genome methods in conjunction with target specific tagmentation, an enzyme-assisted approach to cut and tag DNA. We benchmarked HiCuT against existing methods and found similar efficiency and specificity in assessing target specific 3D genome organization despite reducing assay requirements to 100,000 cells (5-fold less) and 12 million sequencing reads (8-fold fewer). Strikingly, HiCuT data exhibited reduced background noise and required minimal computational processing. Taken together, HiCuT broadens the capacity for 3D genome profiling in cell populations previously unmeasurable, including primary cells and human tissues, reduces sequencing costs, and lowers the need for computational expertise. HiCuT will benefit all investigators studying gene regulation and disease pathophysiology.

Whole-genome methods to define DNA and histone accessibility and long-range interactions in chromatin

Defining the genome-wide chromatin landscape has been a goal of experimentalists for decades. Here we review highlights of these efforts, from seminal experiments showing discontinuities in chromatin structure related to gene activation to extensions of these methods elucidating general features of chromatin related to gene states by exploiting deep sequencing methods. We also review chromatin conformational capture methods to identify patterns in long-range interactions between genomic loci.

Advances in profiling chromatin architecture shed light on the regulatory dynamics underlying brain disorders

Understanding the exquisitely complex nature of the three-dimensional organization of the genome and how it affects gene regulation remains a central question in biology. Recent advances in sequencing- and imaging-based approaches in decoding the three-dimensional chromatin landscape have enabled a systematic characterization of gene regulatory architecture. In this review, we outline how chromatin architecture provides a reference atlas to predict the functional consequences of non-coding variants associated with human traits and disease. High-throughput perturbation assays such as massively parallel reporter assays (MPRA) and CRISPR-based genome engineering in combination with a reference atlas opened an avenue for going beyond observational studies to experimentally validating the regulatory principles of the genome. We conclude by providing a suggested path forward by calling attention to barriers that can be addressed for a more complete understanding of the regulatory landscape of the human brain.

Single-cell measurement of higher-order 3D genome organization with scSPRITE

Although three-dimensional (3D) genome organization is central to many aspects of nuclear function, it has been difficult to measure at the single-cell level. To address this, we developed 'single-cell split-pool recognition of interactions by tag extension' (scSPRITE). scSPRITE uses split-and-pool barcoding to tag DNA fragments in the same nucleus and their 3D spatial arrangement. Because scSPRITE measures multiway DNA contacts, it generates higher-resolution maps within an individual cell than can be achieved by proximity ligation. We applied scSPRITE to thousands of mouse embryonic stem cells and detected known genome structures, including chromosome territories, active and inactive compartments, and topologically associating domains (TADs) as well as long-range inter-chromosomal structures organized around various nuclear bodies. We observe that these structures exhibit different levels of heterogeneity across the population, with TADs representing dynamic units of genome organization across cells. We expect that scSPRITE will be a critical tool for studying genome structure within heterogeneous populations.

Taxus yunnanensis genome offers insights into gymnosperm phylogeny and taxol production

Taxol, a natural product derived from Taxus, is one of the most effective natural anticancer drugs and the biosynthetic pathway of Taxol is the basis of heterologous bio-production. Here, we report a high-quality genome assembly and annotation of Taxus yunnanensis based on 10.7 Gb sequences assembled into 12 chromosomes with contig N50 and scaffold N50 of 2.89 Mb and 966.80 Mb, respectively. Phylogenomic analyses show that T. yunnanensis is most closely related to Sequoiadendron giganteum among the sampled taxa, with an estimated divergence time of 133.4-213.0 MYA. As with most gymnosperms, and unlike most angiosperms, there is no evidence of a recent whole-genome duplication in T. yunnanensis. Repetitive sequences, especially long terminal repeat retrotransposons, are prevalent in the T. yunnanensis genome, contributing to its large genome size. We further integrated genomic and transcriptomic data to unveil clusters of genes involved in Taxol synthesis, located on the chromosome 12, while gene families encoding hydroxylase in the Taxol pathway exhibited significant expansion. Our study contributes to the further elucidation of gymnosperm relationships and the Taxol biosynthetic pathway.

Single-cell landscape of nuclear configuration and gene expression during stem cell differentiation and X inactivation

BACKGROUND

Mammalian development is associated with extensive changes in gene expression, chromatin accessibility, and nuclear structure. Here, we follow such changes associated with mouse embryonic stem cell differentiation and X inactivation by integrating, for the first time, allele-specific data from these three modalities obtained by high-throughput single-cell RNA-seq, ATAC-seq, and Hi-C.

RESULTS

Allele-specific contact decay profiles obtained by single-cell Hi-C clearly show that the inactive X chromosome has a unique profile in differentiated cells that have undergone X inactivation. Loss of this inactive X-specific structure at mitosis is followed by its reappearance during the cell cycle, suggesting a "bookmark" mechanism. Differentiation of embryonic stem cells to follow the onset of X inactivation is associated with changes in contact decay profiles that occur in parallel on both the X chromosomes and autosomes. Single-cell RNA-seq and ATAC-seq show evidence of a delay in female versus male cells, due to the presence of two active X chromosomes at early stages of differentiation. The onset of the inactive X-specific structure in single cells occurs later than gene silencing, consistent with the idea that chromatin compaction is a late event of X inactivation. Single-cell Hi-C highlights evidence of discrete changes in nuclear structure characterized by the acquisition of very long-range contacts throughout the nucleus. Novel computational approaches allow for the effective alignment of single-cell gene expression, chromatin accessibility, and 3D chromosome structure.

CONCLUSIONS

Based on trajectory analyses, three distinct nuclear structure states are detected reflecting discrete and profound simultaneous changes not only to the structure of the X chromosomes, but also to that of autosomes during differentiation. Our study reveals that long-range structural changes to chromosomes appear as discrete events, unlike progressive changes in gene expression and chromatin accessibility.

tagHi-C Reveals 3D Chromatin Architecture Dynamics during Mouse Hematopoiesis

Spatiotemporal chromatin reorganization during hematopoietic differentiation has not been comprehensively characterized, mainly because of the large numbers of starting cells required for current chromatin conformation capture approaches. Here, we introduce a low-input tagmentation-based Hi-C (tagHi-C) method to capture the chromatin structures of hundreds of cells. Using tagHi-C, we are able to map the spatiotemporal dynamics of chromatin structure in ten primary hematopoietic stem, progenitor, and differentiated cell populations from mouse bone marrow. Our results reveal that changes in compartment dynamics and the Rabl configuration occur during hematopoietic cell differentiation. We identify gene-body-associating domains (GADs) as general structures for highly expressed genes. Moreover, we extend the body of knowledge regarding genes influenced by genome-wide association study (GWAS) loci through spatial chromatin looping. Our study provides the tagHi-C method for studying the three-dimensional (3D) genome of a small number of cells and maps the comprehensive 3D chromatin landscape of bone marrow hematopoietic cells.

The 3D nuclear conformation of the major histocompatibility complex changes upon cell activation both in porcine and human macrophages

BACKGROUND

The crucial role of the major histocompatibility complex (MHC) for the immune response to infectious diseases is well-known, but no information is available on the 3D nuclear organization of this gene-dense region in immune cells, whereas nuclear architecture is known to play an essential role on genome function regulation. We analyzed the spatial arrangement of the three MHC regions (class I, III and II) in macrophages using 3D-FISH. Since this complex presents major differences in humans and pigs with, notably, the presence of the centromere between class III and class II regions in pigs, the analysis was implemented in both species to determine the impact of this organization on the 3D conformation of the MHC. The expression level of the three genes selected to represent each MHC region was assessed by quantitative real-time PCR. Resting and lipopolysaccharide (LPS)-activated states were investigated to ascertain whether a response to a pathogen modifies their expression level and their 3D organization.

RESULTS

While the three MHC regions occupy an intermediate radial position in porcine macrophages, the class I region was clearly more peripheral in humans. The BAC center-to-center distances allowed us to propose a 3D nuclear organization of the MHC in each species. LPS/IFNγ activation induces a significant decompaction of the chromatin between class I and class III regions in pigs and between class I and class II regions in humans. We detected a strong overexpression of TNFα (class III region) in both species. Moreover, a single nucleus analysis revealed that the two alleles can have either the same or a different compaction pattern. In addition, macrophage activation leads to an increase in alleles that present a decompacted pattern in humans and pigs.

CONCLUSIONS

The data presented demonstrate that: (i) the MHC harbors a different 3D organization in humans and pigs; (ii) LPS/IFNγ activation induces chromatin decompaction, but it is not the same area affected in the two species. These findings were supported by the application of an original computation method based on the geometrical distribution of the three target genes. Finally, the position of the centromere inside the swine MHC could influence chromatin reorganization during the activation process.

The 3D Genome Structure of Single Cells

The spatial organization of the genome in the cell nucleus is pivotal to cell function. However, how the 3D genome organization and its dynamics influence cellular phenotypes remains poorly understood. The very recent development of single-cell technologies for probing the 3D genome, especially single-cell Hi-C (scHi-C), has ushered in a new era of unveiling cell-to-cell variability of 3D genome features at an unprecedented resolution. Here, we review recent developments in computational approaches to the analysis of scHi-C, including data processing, dimensionality reduction, imputation for enhancing data quality, and the revealing of 3D genome features at single-cell resolution. While much progress has been made in computational method development to analyze single-cell 3D genomes, substantial future work is needed to improve data interpretation and multimodal data integration, which are critical to reveal fundamental connections between genome structure and function among heterogeneous cell populations in various biological contexts.

Chromatin Profiling Techniques: Exploring the Chromatin Environment and Its Contributions to Complex Traits

The genetic architecture of complex traits is multifactorial. Genome-wide association studies (GWASs) have identified risk loci for complex traits and diseases that are disproportionately located at the non-coding regions of the genome. On the other hand, we have just begun to understand the regulatory roles of the non-coding genome, making it challenging to precisely interpret the functions of non-coding variants associated with complex diseases. Additionally, the epigenome plays an active role in mediating cellular responses to fluctuations of sensory or environmental stimuli. However, it remains unclear how exactly non-coding elements associate with epigenetic modifications to regulate gene expression changes and mediate phenotypic outcomes. Therefore, finer interrogations of the human epigenomic landscape in associating with non-coding variants are warranted. Recently, chromatin-profiling techniques have vastly improved our understanding of the numerous functions mediated by the epigenome and DNA structure. Here, we review various chromatin-profiling techniques, such as assays of chromatin accessibility, nucleosome distribution, histone modifications, and chromatin topology, and discuss their applications in unraveling the brain epigenome and etiology of complex traits at tissue homogenate and single-cell resolution. These techniques have elucidated compositional and structural organizing principles of the chromatin environment. Taken together, we believe that high-resolution epigenomic and DNA structure profiling will be one of the best ways to elucidate how non-coding genetic variations impact complex diseases, ultimately allowing us to pinpoint cell-type targets with therapeutic potential.

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.

Single Cell Technologies: Beyond Microfluidics

Single-cell RNA-sequencing (scRNA-seq) has been widely adopted in recent years due to standardized protocols and automation, reliability, and standardized bioinformatic pipelines. The most widely adopted platform is the 10× Genomics solution. Although powerful, this system is limited by its high cost, moderate throughput, and the inability to customize due to fixed kit components. This study will cover new approaches that do not rely on microfluidics and thus have low entry costs, are highly customizable, and are within the reach of any laboratory possessing molecular biology expertise.

Hi-C as a molecular rangefinder to examine genomic rearrangements

The mammalian genome is highly packed into the nucleus. Over the past decade, the development of Hi-C has contributed significantly to our understanding of the three-dimensional (3D) chromatin structure, uncovering the principles and functions of higher-order chromatin organizations. Recent studies have repositioned its property in spatial proximity measurement to address challenging problems in genome analyses including genome assembly, haplotype phasing, and the detection of genomic rearrangements. In particular, the power of Hi-C in detecting large-scale structural variations (SVs) in the cancer genome has been demonstrated, which is challenging to be addressed solely with short-read-based whole-genome sequencing analyses. In this review, we first provide a comprehensive view of Hi-C as an intuitive and effective SV detection tool. Then, we introduce recently developed bioinformatics tools utilizing Hi-C to investigate genomic rearrangements. Finally, we discuss the potential application of single-cell Hi-C to address the heterogeneity of genomic rearrangements and sub-population identification in the cancer genome.

Understanding the Adult Mammalian Heart at Single-Cell RNA-Seq Resolution

During the last decade, extensive efforts have been made to comprehend cardiac cell genetic and functional diversity. Such knowledge allows for the definition of the cardiac cellular interactome as a reasonable strategy to increase our understanding of the normal and pathologic heart. Previous experimental approaches including cell lineage tracing, flow cytometry, and bulk RNA-Seq have often tackled the analysis of cardiac cell diversity as based on the assumption that cell types can be identified by the expression of a single gene. More recently, however, the emergence of single-cell RNA-Seq technology has led us to explore the diversity of individual cells, enabling the cardiovascular research community to redefine cardiac cell subpopulations and identify relevant ones, and even novel cell types, through their cell-specific transcriptomic signatures in an unbiased manner. These findings are changing our understanding of cell composition and in consequence the identification of potential therapeutic targets for different cardiac diseases. In this review, we provide an overview of the continuously changing cardiac cellular landscape, traveling from the pre-single-cell RNA-Seq times to the single cell-RNA-Seq revolution, and discuss the utilities and limitations of this technology.

3D reconstruction of genomic regions from sparse interaction data

Chromosome conformation capture (3C) technologies measure the interaction frequency between pairs of chromatin regions within the nucleus in a cell or a population of cells. Some of these 3C technologies retrieve interactions involving non-contiguous sets of loci, resulting in sparse interaction matrices. One of such 3C technologies is Promoter Capture Hi-C (pcHi-C) that is tailored to probe only interactions involving gene promoters. As such, pcHi-C provides sparse interaction matrices that are suitable to characterize short- and long-range enhancer-promoter interactions. Here, we introduce a new method to reconstruct the chromatin structural (3D) organization from sparse 3C-based datasets such as pcHi-C. Our method allows for data normalization, detection of significant interactions and reconstruction of the full 3D organization of the genomic region despite of the data sparseness. Specifically, it builds, with as low as the 2-3% of the data from the matrix, reliable 3D models of similar accuracy of those based on dense interaction matrices. Furthermore, the method is sensitive enough to detect cell-type-specific 3D organizational features such as the formation of different networks of active gene communities.

X-chromosome regulation and sex differences in brain anatomy

Humans show reproducible sex-differences in cognition and psychopathology that may be contributed to by influences of gonadal sex-steroids and/or sex-chromosomes on regional brain development. Gonadal sex-steroids are well known to play a major role in sexual differentiation of the vertebrate brain, but far less is known regarding the role of sex-chromosomes. Our review focuses on this latter issue by bridging together two literatures that have to date been largely disconnected. We first consider "bottom-up" genetic and molecular studies focused on sex-chromosome gene content and regulation. This literature nominates specific sex-chromosome genes that could drive developmental sex-differences by virtue of their sex-biased expression and their functions within the brain. We then consider the complementary "top down" view, from magnetic resonance imaging studies that map sex- and sex chromosome effects on regional brain anatomy, and link these maps to regional gene-expression within the brain. By connecting these top-down and bottom-up approaches, we emphasize the potential role of X-linked genes in driving sex-biased brain development and outline key goals for future work in this field.

Hi-C implementation of genome structure for in silico models of radiation-induced DNA damage

Developments in the genome organisation field has resulted in the recent methodology to infer spatial conformations of the genome directly from experimentally measured genome contacts (Hi-C data). This provides a detailed description of both intra- and inter-chromosomal arrangements. Chromosomal intermingling is an important driver for radiation-induced DNA mis-repair. Which is a key biological endpoint of relevance to the fields of cancer therapy (radiotherapy), public health (biodosimetry) and space travel. For the first time, we leverage these methods of inferring genome organisation and couple them to nano-dosimetric radiation track structure modelling to predict quantities and distribution of DNA damage within cell-type specific geometries. These nano-dosimetric simulations are highly dependent on geometry and are benefited from the inclusion of experimentally driven chromosome conformations. We show how the changes in Hi-C contract maps impact the inferred geometries resulting in significant differences in chromosomal intermingling. We demonstrate how these differences propagate through to significant changes in the distribution of DNA damage throughout the cell nucleus, suggesting implications for DNA repair fidelity and subsequent cell fate. We suggest that differences in the geometric clustering for the chromosomes between the cell-types are a plausible factor leading to changes in cellular radiosensitivity. Furthermore, we investigate changes in cell shape, such as flattening, and show that this greatly impacts the distribution of DNA damage. This should be considered when comparing in vitro results to in vivo systems. The effect may be especially important when attempting to translate radiosensitivity measurements at the experimental in vitro level to the patient or human level.

Novel orthogonal methods to uncover the complexity and diversity of nuclear architecture

Recent years have seen a vast expansion of knowledge on three-dimensional (3D) genome organization. The majority of studies on chromosome topology consists of pairwise interaction data of bulk populations of cells and therefore conceals heterogenic and more complex folding patterns. Here, we discuss novel methodologies to study the variation in genome topologies between different cells and techniques that allow analysis of complex, multi-way interactions. These technologies will aid the interpretation of genome-wide chromosome conformation data and provide strategies to further dissect the interplay between genome architecture and transcription regulation.

Data integration for inferring context-specific gene regulatory networks

Transcriptional regulatory networks control context-specific gene expression patterns and play important roles in normal and disease processes. Advances in genomics are rapidly increasing our ability to measure different components of the regulation machinery at the single-cell and bulk population level. An important challenge is to combine different types of regulatory genomic measurements to construct a more complete picture of gene regulatory networks across different disease, environmental, and developmental contexts. In this review, we focus on recent computational methods that integrate regulatory genomic data sets to infer context specificity and dynamics in regulatory networks.

Epigenomic technologies for precision oncology

Epigenetic patterns in a cell control the expression of genes and consequently determine the phenotype of a cell. Cancer cells possess altered epigenomes which include aberrant patterns of DNA methylation, histone tail modifications, nucleosome positioning and of the three-dimensional chromatin organization within a nucleus. These altered epigenetic patterns are potential useful biomarkers to detect cancer cells and to classify tumor types. In addition, the cancer epigenome dictates the response of a cancer cell to therapeutic intervention and, therefore its knowledge, will allow to predict response to different therapeutic approaches. Here we review the current state-of-the-art technologies that have been developed to decipher epigenetic patterns on the genomic level and discuss how these methods are potentially useful for precision oncology.

Genomic Cytometry and New Modalities for Deep Single-Cell Interrogation

In the past few years, the rapid development of single-cell analysis techniques has allowed for increasingly in-depth analysis of DNA, RNA, protein, and epigenetic states, at the level of the individual cell. This unprecedented characterization ability has been enabled through the combination of cytometry, microfluidics, genomics, and informatics. Although traditionally discrete, when properly integrated, these fields create the synergistic field of Genomic Cytometry. In this review, we look at the individual methods that together gave rise to the broad field of Genomic Cytometry. We further outline the basic concepts that drive the field and provide a framework to understand this increasingly complex, technology-intensive space. Thus, we introduce Genomic Cytometry as an emerging field and propose that synergistic rationalization of disparate modalities of cytometry, microfluidics, genomics, and informatics under one banner will enable massive leaps forward in the understanding of complex biology. © 2020 International Society for Advancement of Cytometry.

Transcription factors: building hubs in the 3D space

The hierarchical three-dimensional folding of the mammalian genome constitutes an important regulatory layer of gene expression and cell fate control during processes such as development and tumorigenesis. Accumulating evidence supports the existence of complex topological assemblies in which multiple genes and regulatory elements are frequently interacting with each other in the 3D nucleus. Here, we will discuss the nature, organizational principles, and potential function of such assemblies, including the recently reported enhancer “hubs,” “cliques,” and FIREs (frequently interacting regions) as well as multi-contact hubs. We will also review recent studies that investigate the role of transcription factors (TFs) in driving the topological genome reorganization and hub formation in the context of cell fate transitions and cancer. Finally, we will highlight technological advances that enabled these studies, current limitations, and future directions necessary to advance our understating in the field.

Use of Single-Cell -Omic Technologies to Study the Gastrointestinal Tract and Diseases, From Single Cell Identities to Patient Features

Single cells are the building blocks of tissue systems that determine organ phenotypes, behaviors, and functions. Understanding the differences between cell types and their activities might provide us with insights into normal tissue physiology, development of disease, and new therapeutic strategies. Although -omic level single-cell technologies are a relatively recent development that have been used only in research settings, these approaches might eventually be used in the clinic. We review the prospects of applying single-cell genome, transcriptome, epigenome, proteome, and metabolome analyses to gastroenterology and hepatology research. Combining data from multi-omic platforms coupled to rapid technological development could lead to new diagnostic, prognostic, and therapeutic approaches.

Midbrain Dopaminergic Neuron Development at the Single Cell Level: In vivo and in Stem Cells

Parkinson's disease (PD) is a progressive neurodegenerative disorder that predominantly affects dopaminergic (DA) neurons of the substantia nigra. Current treatment options for PD are symptomatic and typically involve the replacement of DA neurotransmission by DA drugs, which relieve the patients of some of their motor symptoms. However, by the time of diagnosis, patients have already lost about 70% of their substantia nigra DA neurons and these drugs offer only temporary relief. Therefore, cell replacement therapy has garnered much interest as a potential treatment option for PD. Early studies using human fetal tissue for transplantation in PD patients provided proof of principle for cell replacement therapy, but they also highlighted the ethical and practical difficulties associated with using human fetal tissue as a cell source. In recent years, advancements in stem cell research have made human pluripotent stem cells (hPSCs) an attractive source of material for cell replacement therapy. Studies on how DA neurons are specified and differentiated in the developing mouse midbrain have allowed us to recapitulate many of the positional and temporal cues needed to generate DA neurons in vitro. However, little is known about the developmental programs that govern human DA neuron development. With the advent of single-cell RNA sequencing (scRNA-seq) and bioinformatics, it has become possible to analyze precious human samples with unprecedented detail and extract valuable high-quality information from large data sets. This technology has allowed the systematic classification of cell types present in the human developing midbrain along with their gene expression patterns. By studying human development in such an unbiased manner, we can begin to elucidate human DA neuron development and determine how much it differs from our knowledge of the rodent brain. Importantly, this molecular description of the function of human cells has become and will increasingly be a reference to define, evaluate, and engineer cell types for PD cell replacement therapy and disease modeling.

IDR2D identifies reproducible genomic interactions

Chromatin interaction data from protocols such as ChIA-PET, HiChIP and Hi-C provide valuable insights into genome organization and gene regulation, but can include spurious interactions that do not reflect underlying genome biology. We introduce an extension of the Irreproducible Discovery Rate (IDR) method called IDR2D that identifies replicable interactions shared by chromatin interaction experiments. IDR2D provides a principled set of interactions and eliminates artifacts from single experiments. The method is available as a Bioconductor package for the R community, as well as an online service at https://idr2d.mit.edu.

Radial Organization in the Mammalian Nucleus

In eukaryotic cells, most of the genetic material is contained within a highly specialized organelle-the nucleus. A large body of evidence indicates that, within the nucleus, chromatinized DNA is spatially organized at multiple length scales. The higher-order organization of chromatin is crucial for proper execution of multiple genome functions, including DNA replication and transcription. Here, we review our current knowledge on the spatial organization of chromatin in the nucleus of mammalian cells, focusing in particular on how chromatin is radially arranged with respect to the nuclear lamina. We then discuss the possible mechanisms by which the radial organization of chromatin in the cell nucleus is established. Lastly, we propose a unifying model of nuclear spatial organization, and suggest novel approaches to test it.