An Overview of Genome Organization and How We Got There: from FISH to Hi-C

Genome structural dynamics: insights from Gaussian network analysis of Hi-C data

Characterization of the spatiotemporal properties of the chromatin is essential to gaining insights into the physical bases of gene co-expression, transcriptional regulation and epigenetic modifications. The Gaussian network model (GNM) has proven in recent work to serve as a useful tool for modeling chromatin structural dynamics, using as input high-throughput chromosome conformation capture data. We focus here on the exploration of the collective dynamics of chromosomal structures at hierarchical levels of resolution, from single gene loci to topologically associating domains or entire chromosomes. The GNM permits us to identify long-range interactions between gene loci, shedding light on the role of cross-correlations between distal regions of the chromosomes in regulating gene expression. Notably, GNM analysis performed across diverse cell lines highlights the conservation of the global/cooperative movements of the chromatin across different types of cells. Variations driven by localized couplings between genomic loci, on the other hand, underlie cell differentiation, underscoring the significance of the four-dimensional properties of the genome in defining cellular identity. Finally, we demonstrate the close relation between the cell type-dependent mobility profiles of gene loci and their gene expression patterns, providing a clear demonstration of the role of chromosomal 4D features in defining cell-specific differential expression of genes.

Regulation of 3D genome organization during T cell activation

Within the three-dimensional (3D) nuclear space, the genome organizes into a series of orderly structures that impose important influences on gene regulation. T lymphocytes, crucial players in adaptive immune responses, undergo intricate transcriptional remodeling upon activation, leading to differentiation into specific effector and memory T cell subsets. Recent evidence suggests that T cell activation is accompanied by dynamic changes in genome architecture at multiple levels, providing a unique biological context to explore the functional relevance and molecular mechanisms of 3D genome organization. Here, we summarize recent advances that link the reorganization of genome architecture to the remodeling of transcriptional programs and conversion of cell fates during T cell activation and differentiation. We further discuss how various chromatin architecture regulators, including CCCTC-binding factor and several transcription factors, collectively modulate the genome architecture during this process.

Structural basis for the preservation of a subset of topologically associating domains in interphase chromosomes upon cohesin depletion

Compartment formation in interphase chromosomes is a result of spatial segregation between euchromatin and heterochromatin on a few megabase pairs (Mbp) scale. On the sub-Mbp scales, topologically associating domains (TADs) appear as interacting domains along the diagonal in the ensemble averaged Hi-C contact map. Hi-C experiments showed that most of the TADs vanish upon deleting cohesin, while the compartment structure is maintained, and perhaps even enhanced. However, closer inspection of the data reveals that a non-negligible fraction of TADs is preserved (P-TADs) after cohesin loss. Imaging experiments show that, at the single-cell level, TAD-like structures are present even without cohesin. To provide a structural basis for these findings, we first used polymer simulations to show that certain TADs with epigenetic switches across their boundaries survive after depletion of loops. More importantly, the three-dimensional structures show that many of the P-TADs have sharp physical boundaries. Informed by the simulations, we analyzed the Hi-C maps (with and without cohesin) in mouse liver and human colorectal carcinoma cell lines, which affirmed that epigenetic switches and physical boundaries (calculated using the predicted 3D structures using the data-driven HIPPS method that uses Hi-C as the input) explain the origin of the P-TADs. Single-cell structures display TAD-like features in the absence of cohesin that are remarkably similar to the findings in imaging experiments. Some P-TADs, with physical boundaries, are relevant to the retention of enhancer-promoter/promoter-promoter interactions. Overall, our study shows that preservation of a subset of TADs upon removing cohesin is a robust phenomenon that is valid across multiple cell lines.

Predicting the unpredictable: the regulatory nature and promiscuity of herbicide cross resistance

The emergence of herbicide-resistant weeds is a significant threat to modern agriculture. Cross resistance, a phenomenon where resistance to one herbicide confers resistance to another, is a particular concern owing to its unpredictability. Nontarget-site (NTS) cross resistance is especially challenging to predict, as it arises from genes that encode enzymes that do not directly involve the herbicide target site and can affect multiple herbicides. Recent advancements in genomic and structural biology techniques could provide new venues for predicting NTS resistance in weed species. In this review, we present an overview of the latest approaches that could be used. We discuss the use of genomic and epigenomics techniques such as ATAC-seq and DAP-seq to identify transcription factors and cis-regulatory elements associated with resistance traits. Enzyme/protein structure prediction and docking analysis are discussed as an initial step for predicting herbicide binding affinities with key enzymes to identify candidates for subsequent in vitro validation. We also provide example analyses that can be deployed toward elucidating cross resistance and its regulatory patterns. Ultimately, our review provides important insights into the latest scientific advancements and potential directions for predicting and managing herbicide cross resistance in weeds. © 2023 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

scENCORE: leveraging single-cell epigenetic data to predict chromatin conformation using graph embedding

Dynamic compartmentalization of eukaryotic DNA into active and repressed states enables diverse transcriptional programs to arise from a single genetic blueprint, whereas its dysregulation can be strongly linked to a broad spectrum of diseases. While single-cell Hi-C experiments allow for chromosome conformation profiling across many cells, they are still expensive and not widely available for most labs. Here, we propose an alternate approach, scENCORE, to computationally reconstruct chromatin compartments from the more affordable and widely accessible single-cell epigenetic data. First, scENCORE constructs a long-range epigenetic correlation graph to mimic chromatin interaction frequencies, where nodes and edges represent genome bins and their correlations. Then, it learns the node embeddings to cluster genome regions into A/B compartments and aligns different graphs to quantify chromatin conformation changes across conditions. Benchmarking using cell-type-matched Hi-C experiments demonstrates that scENCORE can robustly reconstruct A/B compartments in a cell-type-specific manner. Furthermore, our chromatin confirmation switching studies highlight substantial compartment-switching events that may introduce substantial regulatory and transcriptional changes in psychiatric disease. In summary, scENCORE allows accurate and cost-effective A/B compartment reconstruction to delineate higher-order chromatin structure heterogeneity in complex tissues.

Large-scale data-driven and physics-based models offer insights into the relationships among the structures, dynamics, and functions of chromosomes

The organized three-dimensional chromosome architecture in the cell nucleus provides scaffolding for precise regulation of gene expression. When the cell changes its identity in the cell-fate decision-making process, extensive rearrangements of chromosome structures occur accompanied by large-scale adaptations of gene expression, underscoring the importance of chromosome dynamics in shaping genome function. Over the last two decades, rapid development of experimental methods has provided unprecedented data to characterize the hierarchical structures and dynamic properties of chromosomes. In parallel, these enormous data offer valuable opportunities for developing quantitative computational models. Here, we review a variety of large-scale polymer models developed to investigate the structures and dynamics of chromosomes. Different from the underlying modeling strategies, these approaches can be classified into data-driven ('top-down') and physics-based ('bottom-up') categories. We discuss their contributions to offering valuable insights into the relationships among the structures, dynamics, and functions of chromosomes and propose the perspective of developing data integration approaches from different experimental technologies and multidisciplinary theoretical/simulation methods combined with different modeling strategies.

Leveraging polymer modeling to reconstruct chromatin connectivity from live images

Chromosomal dynamics plays a central role in a number of critical biological processes, such as transcriptional regulation, genetic recombination, and DNA replication. However, visualization of chromatin is generally limited to live imaging of a few fluorescently labeled chromosomal loci or high-resolution reconstruction of multiple loci from a single time frame. To aid in mapping the underlying chromosomal structure based on parsimonious experimental measurements, we present an exact analytical expression for the evolution of the polymer configuration based on a flexible-polymer model, and we propose an algorithm that tracks the polymer configuration from live images of chromatin marked with several fluorescent marks. Our theory identifies the resolution of microscopy needed to achieve high-accuracy tracking for a given spacing of markers, establishing the statistical confidence in the assignment of genome identity to the visualized marks. We then leverage experimental data of locus-tracking measurements to demonstrate the validity of our modeling approach and to establish a basis for the design of experiments with a desired resolution. Altogether, this work provides a computational approach founded on polymer physics that vastly improves the interpretation of in vivo measurements of biopolymer dynamics.

Fundamental insights into the correlation between chromosome configuration and transcription

Eukaryotic chromosomes exhibit a hierarchical organization that spans a spectrum of length scales, ranging from sub-regions known as loops, which typically comprise hundreds of base pairs, to much larger chromosome territories that can encompass a few mega base pairs. Chromosome conformation capture experiments that involve high-throughput sequencing methods combined with microscopy techniques have enabled a new understanding of inter- and intra-chromosomal interactions with unprecedented details. This information also provides mechanistic insights on the relationship between genome architecture and gene expression. In this article, we review the recent findings on three-dimensional interactions among chromosomes at the compartment, topologically associating domain, and loop levels and the impact of these interactions on the transcription process. We also discuss current understanding of various biophysical processes involved in multi-layer structural organization of chromosomes. Then, we discuss the relationships between gene expression and genome structure from perturbative genome-wide association studies. Furthermore, for a better understanding of how chromosome architecture and function are linked, we emphasize the role of epigenetic modifications in the regulation of gene expression. Such an understanding of the relationship between genome architecture and gene expression can provide a new perspective on the range of potential future discoveries and therapeutic research.

Model chromatin flows: numerical analysis of linear and nonlinear hydrodynamics inside a sphere

We solve a hydrodynamic model of active chromatin dynamics, within a confined geometry simulating the cell nucleus. Using both analytical and numerical methods, we describe the behavior of the chromatin polymer driven by the activity of motors having polar symmetry, both in the linear response regime as well as in the long-term, fully nonlinear regime of the flows. The introduction of a boundary induces a particular geometry in the flows of chromatin, which we describe using vector spherical harmonics, a tool which greatly simplifies both our analytical and numerical approaches. We find that the long-term behavior of this model in confinement is dominated by steady, transverse flows of chromatin which circulate around the spherical domain. These circulating flows are found to be robust to perturbations, and their characteristic size is set by the size of the domain. This gives us further insight into active chromatin dynamics in the cell nucleus, and provides a foundation for development of further, more complex models of active chromatin dynamics.

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.

The Talented LncRNAs: Meshing into Transcriptional Regulatory Networks in Cancer

As a group of diseases characterized by uncontrollable cell growth, cancer is highly multifaceted in how it overrides checkpoints controlling proliferation. Amongst the regulators of these checkpoints, long non-coding RNAs (lncRNAs) can have key roles in why natural biological processes go haywire. LncRNAs represent a large class of regulatory transcripts that can localize anywhere in cells. They were found to affect gene expression on many levels from transcription to mRNA translation and even protein stability. LncRNA participation in such control mechanisms can depend on cell context, with given transcripts sometimes acting as oncogenes or tumor suppressors. Importantly, the tissue-specificity and low expression levels of lncRNAs make them attractive therapeutic targets or biomarkers. Here, we review the various cellular processes affected by lncRNAs and outline molecular strategies they use to control gene expression, particularly in cancer and in relation to transcription factors.

Assignment of the somatic A/B compartments to chromatin domains in giant transcriptionally active lampbrush chromosomes

BACKGROUND

The three-dimensional configuration of the eukaryotic genome is an emerging area of research. Chromosome conformation capture outlined genome segregation into large scale A and B compartments corresponding mainly to transcriptionally active and repressive chromatin. It remains unknown how the compartmentalization of the genome changes in growing oocytes of animals with hypertranscriptional type of oogenesis. Such oocytes are characterized by highly elongated chromosomes, called lampbrush chromosomes, which acquire a typical chromomere-loop appearance, representing one of the classical model systems for exploring the structural and functional organization of chromatin domains.

RESULTS

Here, we compared the distribution of A/B compartments in chicken somatic cells with chromatin domains in lampbrush chromosomes. We found that in lampbrush chromosomes, the extended chromatin domains, restricted by compartment boundaries in somatic cells, disintegrate into individual chromomeres. Next, we performed FISH-mapping of the genomic loci, which belong to A or B chromatin compartments as well as to A/B compartment transition regions in embryonic fibroblasts on isolated lampbrush chromosomes. We found, that in chicken lampbrush chromosomes, clusters of dense compact chromomeres bearing short lateral loops and enriched with repressive epigenetic modifications generally correspond to constitutive B compartments in somatic cells. A compartments align with lampbrush chromosome segments with smaller, less compact chromomeres, longer lateral loops, and a higher transcriptional status. Clusters of small loose chromomeres with relatively long lateral loops show no obvious correspondence with either A or B compartment identity. Some genes belonging to facultative B (sub-) compartments can be tissue-specifically transcribed during oogenesis, forming distinct lateral loops.

CONCLUSIONS

Here, we established a correspondence between the A/B compartments in somatic interphase nucleus and chromatin segments in giant lampbrush chromosomes from diplotene stage oocytes. The chromomere-loop structure of the genomic regions corresponding to interphase A and B compartments reveals the difference in how they are organized at the level of chromatin domains. The results obtained also suggest that gene-poor regions tend to be packed into chromomeres.

Theory of chromatin organization maintained by active loop extrusion

The active loop extrusion hypothesis proposes that chromatin threads through the cohesin protein complex into progressively larger loops until reaching specific boundary elements. We build upon this hypothesis and develop an analytical theory for active loop extrusion which predicts that loop formation probability is a nonmonotonic function of loop length and describes chromatin contact probabilities. We validate our model with Monte Carlo and hybrid Molecular Dynamics-Monte Carlo simulations and demonstrate that our theory recapitulates experimental chromatin conformation capture data. Our results support active loop extrusion as a mechanism for chromatin organization and provide an analytical description of chromatin organization that may be used to specifically modify chromatin contact probabilities.

Dynamic chromatin architectures provide insights into the genetics of cattle myogenesis

BACKGROUND

Sharply increased beef consumption is propelling the genetic improvement projects of beef cattle in China. Three-dimensional genome structure is confirmed to be an important layer of transcription regulation. Although genome-wide interaction data of several livestock species have already been produced, the genome structure states and its regulatory rules in cattle muscle are still limited.

RESULTS

Here we present the first 3D genome data in Longissimus dorsi muscle of fetal and adult cattle (Bos taurus). We showed that compartments, topologically associating domains (TADs), and loop undergo re-organization and the structure dynamics were consistent with transcriptomic divergence during muscle development. Furthermore, we annotated cis-regulatory elements in cattle genome during myogenesis and demonstrated the enrichments of promoter and enhancer in selection sweeps. We further validated the regulatory function of one HMGA2 intronic enhancer near a strong sweep region on primary bovine myoblast proliferation.

CONCLUSIONS

Our data provide key insights of the regulatory function of high order chromatin structure and cattle myogenic biology, which will benefit the progress of genetic improvement of beef cattle.

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

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

Exploring high-resolution chromatin interaction changes and functional enhancers of myogenic marker genes during myogenic differentiation

Skeletal muscle differentiation (myogenesis) is a complex and highly coordinated biological process regulated by a series of myogenic marker genes. Chromatin interactions between gene's promoters and their enhancers have an important role in transcriptional control. However, the high-resolution chromatin interactions of myogenic genes and their functional enhancers during myogenesis remain largely unclear. Here, we used circularized chromosome conformation capture coupled with next generation sequencing (4C-seq) to investigate eight myogenic marker genes in C2C12 myoblasts (C2C12-MBs) and C2C12 myotubes (C2C12-MTs). We revealed dynamic chromatin interactions of these marker genes during differentiation and identified 163 and 314 significant interaction sites (SISs) in C2C12-MBs and C2C12-MTs, respectively. The interacting genes of SISs in C2C12-MTs were mainly involved in muscle development, and histone modifications of the SISs changed during differentiation. Through functional genomic screening, we also identified 25 and 41 putative active enhancers in C2C12-MBs and C2C12-MTs, respectively. Using luciferase reporter assays for putative enhancers of Myog and Myh3, we identified eight activating enhancers. Furthermore, dCas9-KRAB epigenome editing and RNA-Seq revealed a role for Myog enhancers in the regulation of Myog expression and myogenic differentiation in the native genomic context. Taken together, this study lays the groundwork for understanding 3D chromatin interaction changes of myogenic genes during myogenesis and provides insights that contribute to our understanding of the role of enhancers in regulating myogenesis.

Identification of enhancers responsible for the coordinated expression of myosin heavy chain isoforms in skeletal muscle

Background

Skeletal muscles consist of fibers of differing contractility and metabolic properties, which are primarily determined by the content of myosin heavy chain (MYH) isoforms (MYH7, MYH2, MYH1, and MYH4). The regulation of Myh genes transcription depends on three-dimensional chromatin conformation interaction, but the mechanistic details remain to be determined.

Results

In this study, we characterized the interaction profiles of Myh genes using 4C-seq (circular chromosome conformation capture coupled to high-throughput sequencing). The interaction profile of Myh genes changed between fast quadriceps and slow soleus muscles. Combining chromatin immunoprecipitation-sequencing (ChIP-seq) and transposase accessible chromatin with high-throughput sequencing (ATAC-seq), we found that a 38 kb intergenic region interacting simultaneously with fast Myh genes promoters controlled the coordinated expression of fast Myh genes. We also identified four active enhancers of Myh7, and revealed that binding of MYOG and MYOD increased the activity of Myh7 enhancers.

Conclusions

This study provides new insight into the chromatin interactions that regulate Myh genes expression.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08737-9.

Dynamic regulatory module networks for inference of cell type-specific transcriptional networks

Changes in transcriptional regulatory networks can significantly alter cell fate. To gain insight into transcriptional dynamics, several studies have profiled bulk multi-omic data sets with parallel transcriptomic and epigenomic measurements at different stages of a developmental process. However, integrating these data to infer cell type-specific regulatory networks is a major challenge. We present dynamic regulatory module networks (DRMNs), a novel approach to infer cell type-specific cis-regulatory networks and their dynamics. DRMN integrates expression, chromatin state, and accessibility to predict cis-regulators of context-specific expression, where context can be cell type, developmental stage, or time point, and uses multitask learning to capture network dynamics across linearly and hierarchically related contexts. We applied DRMNs to study regulatory network dynamics in three developmental processes, each showing different temporal relationships and measuring a different combination of regulatory genomic data sets: cellular reprogramming, liver dedifferentiation, and forward differentiation. DRMN identified known and novel regulators driving cell type-specific expression patterns, showing its broad applicability to examine dynamics of gene regulatory networks from linearly and hierarchically related multi-omic data sets.

LRF Promotes Indirectly Advantageous Chromatin Conformation via BGLT3-lncRNA Expression and Switch from Fetal to Adult Hemoglobin

The hemoglobin switch from fetal (HbF) to adult (HbA) has been studied intensively as an essential model for gene expression regulation, but also as a beneficial therapeutic approach for β-hemoglobinopathies, towards the objective of reactivating HbF. The transcription factor LRF (Leukemia/lymphoma-related), encoded from the ZBTB7A gene has been implicated in fetal hemoglobin silencing, though has a wide range of functions that have not been fully clarified. We thus established the LRF/ZBTB7A-overexpressing and ZBTB7A-knockdown K562 (human erythroleukemia cell line) clones to assess fetal vs. adult hemoglobin production pre- and post-induction. Transgenic K562 clones were further developed and studied under the influence of epigenetic chromatin regulators, such as DNA methyl transferase 3 (DNMT3) and Histone Deacetylase 1 (HDAC1), to evaluate LRF’s potential disturbance upon the aberrant epigenetic background and provide valuable information of the preferable epigenetic frame, in which LRF unfolds its action on the β-type globin’s expression. The ChIP-seq analysis demonstrated that LRF binds to γ-globin genes (HBG2/1) and apparently associates BCL11A for their silencing, but also during erythropoiesis induction, LRF binds the BGLT3 gene, promoting BGLT3-lncRNA production through the γ-δ intergenic region of β-type globin’s locus, triggering the transcriptional events from γ- to β-globin switch. Our findings are supported by an up-to-date looping model, which highlights chromatin alterations during erythropoiesis at late stages of gestation, to establish an “open” chromatin conformation across the γ-δ intergenic region and accomplish β-globin expression and hemoglobin switch.

Monitoring antibiotic resistance genes in wastewater environments: The challenges of filling a gap in the One-Health cycle

Antibiotic resistance (AR) is a global problem requiring international cooperation and coordinated action. Global monitoring must rely on methods available and comparable across nations to quantify AR occurrence and identify sources and reservoirs, as well as paths of AR dissemination. Numerous analytical tools that are gaining relevance in microbiology, have the potential to be applied to AR research. This review summarizes the state of the art of AR monitoring methods, considering distinct needs, objectives and available resources. Based on the overview of distinct approaches that are used or can be adapted to monitor AR, it is discussed the potential to establish reliable and useful monitoring schemes that can be implemented in distinct contexts. This discussion places the environmental monitoring within the One-Health approach, where two types of risk, dissemination across distinct environmental compartments, and transmission to humans, must be considered. The plethora of methodological approaches to monitor AR and the variable features of the monitored sites challenge the capacity of the scientific community and policy makers to reach a common understanding. However, the dialogue between different methods and the production of action-oriented data is a priority. The review aims to warm up this discussion.

Plant 3D Chromatin Organization: Important Insights from Chromosome Conformation Capture Analyses of the Last 10 Years

Over the past few decades, eukaryotic linear genomes and epigenomes have been widely and extensively studied for understanding gene expression regulation. More recently, the three-dimensional (3D) chromatin organization was found to be important for determining genome functionality, finely tuning physiological processes for appropriate cellular responses. With the development of visualization techniques and chromatin conformation capture (3C)-based techniques, increasing evidence indicates that chromosomal architecture characteristics and chromatin domains with different epigenetic modifications in the nucleus are correlated with transcriptional activities. Subsequent studies have further explored the intricate interplay between 3D genome organization and the function of interacting regions. In this review, we summarize spatial distribution patterns of chromatin, including chromatin positioning, configurations and domains, with a particular focus on the effect of a unique form of interaction between varieties of factors that shape the 3D genome conformation in plants. We further discuss the methods, advantages and limitations of various 3C-based techniques, highlighting the applications of these technologies in plants to identify chromatin domains, and address their dynamic changes and functional implications in evolution, and adaptation to development and changing environmental conditions. Moreover, the future implications and emerging research directions of 3D genome organization are discussed.

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.

Four-dimensional chromosome reconstruction elucidates the spatiotemporal reorganization of the mammalian X chromosome

Chromosomes are segmented into domains and compartments, but how these structures are spatially related in three dimensions (3D) is unclear. Here, we developed tools that directly extract 3D information from Hi-C experiments and integrate the data across time. With our "4DHiC" method, we use X chromosome inactivation (XCI) as a model to examine the time evolution of 3D chromosome architecture during large-scale changes in gene expression. Our modeling resulted in several insights. Both A/B and S1/S2 compartments divide the X chromosome into hemisphere-like structures suggestive of a spatial phase-separation. During the XCI, the X chromosome transits through A/B, S1/S2, and megadomain structures by undergoing only partial mixing to assume new structures. Interestingly, when an active X chromosome (Xa) is reorganized into an inactive X chromosome (Xi), original underlying compartment structures are not fully eliminated within the Xi superstructure. Our study affirms slow mixing dynamics in the inner chromosome core and faster dynamics near the surface where escapees reside. Once established, the Xa and Xi resemble glassy polymers where mixing no longer occurs. Finally, Xist RNA molecules initially reside within the A compartment but transition to the interface between the A and B hemispheres and then spread between hemispheres via both surface and core to establish the Xi.

Visualizing Live Chromatin Dynamics through CRISPR-Based Imaging Techniques

The three-dimensional organization of chromatin and its time-dependent changes greatly affect virtually every cellular function, especially DNA replication, genome maintenance, transcription regulation, and cell differentiation. Sequencing-based techniques such as ChIP-seq, ATAC-seq, and Hi-C provide abundant information on how genomic elements are coupled with regulatory proteins and functionally organized into hierarchical domains through their interactions. However, visualizing the time-dependent changes of such organization in individual cells remains challenging. Recent developments of CRISPR systems for site-specific fluorescent labeling of genomic loci have provided promising strategies for visualizing chromatin dynamics in live cells. However, there are several limiting factors, including background signals, off-target binding of CRISPR, and rapid photobleaching of the fluorophores, requiring a large number of target-bound CRISPR complexes to reliably distinguish the target-specific foci from the background. Various modifications have been engineered into the CRISPR system to enhance the signal-to-background ratio and signal longevity to detect target foci more reliably and efficiently, and to reduce the required target size. In this review, we comprehensively compare the performances of recently developed CRISPR designs for improved visualization of genomic loci in terms of the reliability of target detection, the ability to detect small repeat loci, and the allowed time of live tracking. Longer observation of genomic loci allows the detailed identification of the dynamic characteristics of chromatin. The diffusion properties of chromatin found in recent studies are reviewed, which provide suggestions for the underlying biological processes.

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.

Gravitational Force-Induced 3D Chromosomal Conformational Changes Are Associated with Rapid Transcriptional Response in Human T Cells

The mechanisms underlying gravity perception in mammalian cells are unknown. We have recently discovered that the transcriptome of cells in the immune system, which is the most affected system during a spaceflight, responds rapidly and broadly to altered gravity. To pinpoint potential underlying mechanisms, we compared gene expression and three-dimensional (3D) chromosomal conformational changes in human Jurkat T cells during the short-term gravitational changes in parabolic flight and suborbital ballistic rocket flight experiments. We found that differential gene expression in gravity-responsive chromosomal regions, but not differentially regulated single genes, are highly conserved between different real altered gravity comparisons. These coupled gene expression effects in chromosomal regions could be explained by underlying chromatin structures. Based on a high-throughput chromatin conformation capture (Hi-C) analysis in altered gravity, we found that small chromosomes (chr16–22, with the exception of chr18) showed increased intra- and interchromosomal interactions in altered gravity, whereby large chromosomes showed decreased interactions. Finally, we detected a nonrandom overlap between Hi-C-identified chromosomal interacting regions and gravity-responsive chromosomal regions (GRCRs). We therefore demonstrate the first evidence that gravitational force-induced 3D chromosomal conformational changes are associated with rapid transcriptional response in human T cells. We propose a general model of cellular sensitivity to gravitational forces, where gravitational forces acting on the cellular membrane are rapidly and mechanically transduced through the cytoskeleton into the nucleus, moving chromosome territories to new conformation states and their genes into more expressive or repressive environments, finally resulting in region-specific differential gene expression.

Conformational statistics of non-equilibrium polymer loops in Rouse model with active loop extrusion

Motivated by the recent experimental observations of the DNA loop extrusion by protein motors, in this paper, we investigate the statistical properties of the growing polymer loops within the ideal chain model. The loop conformation is characterized statistically by the mean gyration radius and the pairwise contact probabilities. It turns out that a single dimensionless parameter, which is given by the ratio of the loop relaxation time over the time elapsed since the start of extrusion, controls the crossover between near-equilibrium and highly non-equilibrium asymptotics in the statistics of the extruded loop, regardless of the specific time dependence of the extrusion velocity. In addition, we show that two-sided and one-sided loop extruding motors produce the loops with almost identical properties. Our predictions are based on two rigorous semi-analytical methods accompanied by asymptotic analysis of slow and fast extrusion limits.

Multiscale modeling of genome organization with maximum entropy optimization

Three-dimensional (3D) organization of the human genome plays an essential role in all DNA-templated processes, including gene transcription, gene regulation, and DNA replication. Computational modeling can be an effective way of building high-resolution genome structures and improving our understanding of these molecular processes. However, it faces significant challenges as the human genome consists of over 6 × 109 base pairs, a system size that exceeds the capacity of traditional modeling approaches. In this perspective, we review the progress that has been made in modeling the human genome. Coarse-grained models parameterized to reproduce experimental data via the maximum entropy optimization algorithm serve as effective means to study genome organization at various length scales. They have provided insight into the principles of whole-genome organization and enabled de novo predictions of chromosome structures from epigenetic modifications. Applications of these models at a near-atomistic resolution further revealed physicochemical interactions that drive the phase separation of disordered proteins and dictate chromatin stability in situ. We conclude with an outlook on the opportunities and challenges in studying chromosome dynamics.

RUNX1-mediated alphaherpesvirus-host trans-species chromatin interaction promotes viral transcription

Like most DNA viruses, herpesviruses precisely deliver their genomes into the sophisticatedly organized nuclei of the infected host cells to initiate subsequent transcription and replication. However, it remains elusive how the viral genome specifically interacts with the host genome and hijacks host transcription machinery. Using pseudorabies virus (PRV) as model virus, we performed chromosome conformation capture assays to demonstrate a genome-wide specific trans-species chromatin interaction between the virus and host. Our data show that the PRV genome is delivered by the host DNA binding protein RUNX1 into the open chromatin and active transcription zone. This facilitates virus hijacking host RNAPII to efficiently transcribe viral genes, which is significantly inhibited by either a RUNX1 inhibitor or RNA interference. Together, these findings provide insights into the chromatin interaction between viral and host genomes and identify new areas of research to advance the understanding of herpesvirus genome transcription.

Inferring Single-Cell 3D Chromosomal Structures Based on the Lennard-Jones Potential

Reconstructing three-dimensional (3D) chromosomal structures based on single-cell Hi-C data is a challenging scientific problem due to the extreme sparseness of the single-cell Hi-C data. In this research, we used the Lennard-Jones potential to reconstruct both 500 kb and high-resolution 50 kb chromosomal structures based on single-cell Hi-C data. A chromosome was represented by a string of 500 kb or 50 kb DNA beads and put into a 3D cubic lattice for simulations. A 2D Gaussian function was used to impute the sparse single-cell Hi-C contact matrices. We designed a novel loss function based on the Lennard-Jones potential, in which the ε value, i.e., the well depth, was used to indicate how stable the binding of every pair of beads is. For the bead pairs that have single-cell Hi-C contacts and their neighboring bead pairs, the loss function assigns them stronger binding stability. The Metropolis-Hastings algorithm was used to try different locations for the DNA beads, and simulated annealing was used to optimize the loss function. We proved the correctness and validness of the reconstructed 3D structures by evaluating the models according to multiple criteria and comparing the models with 3D-FISH data.

Balancing cohesin eviction and retention prevents aberrant chromosomal interactions, Polycomb-mediated repression, and X-inactivation

Depletion of architectural factors globally alters chromatin structure but only modestly affects gene expression. We revisit the structure-function relationship using the inactive X chromosome (Xi) as a model. We investigate cohesin imbalances by forcing its depletion or retention using degron-tagged RAD21 (cohesin subunit) or WAPL (cohesin release factor). Cohesin loss disrupts the Xi superstructure, unveiling superloops between escapee genes with minimal effect on gene repression. By contrast, forced cohesin retention markedly affects Xi superstructure, compromises spreading of Xist RNA-Polycomb complexes, and attenuates Xi silencing. Effects are greatest at distal chromosomal ends, where looping contacts with the Xist locus are weakened. Surprisingly, cohesin loss creates an Xi superloop, and cohesin retention creates Xi megadomains on the active X chromosome. Across the genome, a proper cohesin balance protects against aberrant inter-chromosomal interactions and tempers Polycomb-mediated repression. We conclude that a balance of cohesin eviction and retention regulates X inactivation and inter-chromosomal interactions across the genome.

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.

Tracing DNA paths and RNA profiles in cultured cells and tissues with ORCA

Chromatin conformation capture (3C) methods and fluorescent in situ hybridization (FISH) microscopy have been used to investigate the spatial organization of the genome. Although powerful, both techniques have limitations. Hi-C is challenging for low cell numbers and requires very deep sequencing to achieve its high resolution. In contrast, FISH can be done on small cell numbers and capture rare cell populations, but typically targets pairs of loci at a lower resolution. Here we detail a protocol for optical reconstruction of chromatin architecture (ORCA), a microscopy approach to trace the 3D DNA path within the nuclei of fixed tissues and cultured cells with a genomic resolution as fine as 2 kb and a throughput of ~10,000 cells per experiment. ORCA can identify structural features with comparable resolution to Hi-C while providing single-cell resolution and multimodal measurements characteristic of microscopy. We describe how to use this DNA labeling in parallel with multiplexed labeling of dozens of RNAs to relate chromatin structure and gene expression in the same cells. Oligopaint probe design, primary probe making, sample collection, cryosectioning and RNA/DNA primary probe hybridization can be completed in 1.5 weeks, while automated RNA/DNA barcode hybridization and RNA/DNA imaging typically takes 2-6 d for data collection and 2-7 d for the automated steps of image analysis.

In situ genome sequencing resolves DNA sequence and structure in intact biological samples

Understanding genome organization requires integration of DNA sequence and three-dimensional spatial context; however, existing genome-wide methods lack either base pair sequence resolution or direct spatial localization. Here, we describe in situ genome sequencing (IGS), a method for simultaneously sequencing and imaging genomes within intact biological samples. We applied IGS to human fibroblasts and early mouse embryos, spatially localizing thousands of genomic loci in individual nuclei. Using these data, we characterized parent-specific changes in genome structure across embryonic stages, revealed single-cell chromatin domains in zygotes, and uncovered epigenetic memory of global chromosome positioning within individual embryos. These results demonstrate how IGS can directly connect sequence and structure across length scales from single base pairs to whole organisms.

Order and stochasticity in the folding of individual Drosophila genomes

Mammalian and Drosophila genomes are partitioned into topologically associating domains (TADs). Although this partitioning has been reported to be functionally relevant, it is unclear whether TADs represent true physical units located at the same genomic positions in each cell nucleus or emerge as an average of numerous alternative chromatin folding patterns in a cell population. Here, we use a single-nucleus Hi-C technique to construct high-resolution Hi-C maps in individual Drosophila genomes. These maps demonstrate chromatin compartmentalization at the megabase scale and partitioning of the genome into non-hierarchical TADs at the scale of 100 kb, which closely resembles the TAD profile in the bulk in situ Hi-C data. Over 40% of TAD boundaries are conserved between individual nuclei and possess a high level of active epigenetic marks. Polymer simulations demonstrate that chromatin folding is best described by the random walk model within TADs and is most suitably approximated by a crumpled globule build of Gaussian blobs at longer distances. We observe prominent cell-to-cell variability in the long-range contacts between either active genome loci or between Polycomb-bound regions, suggesting an important contribution of stochastic processes to the formation of the Drosophila 3D genome.

Genomes are partitioned into topologically associating domains (TADs). Here the authors present single-nucleus Hi-C maps in Drosophila at 10 kb resolution, demonstrating the presence of chromatin compartments in individual nuclei, and partitioning of the genome into non-hierarchical TADs at the scale of 100 kb, which resembles population TAD profiles.

Genetic approaches to revealing the principles of nuclear architecture

The spatial organization of chromosomes inside the eukaryotic nucleus is important for DNA replication, repair and gene expression. During development of multicellular organisms, different compendiums of genes are either repressed or activated to produce specific cell types. Genetic manipulation of tractable organisms is invaluable to elucidate chromosome configuration and the underlying mechanisms. Systematic inhibition of genes through RNA interference and, more recently, CRISPR/Cas9-based screens have identified new proteins with significant roles in nuclear organization. Coupling this with advances in imaging techniques, such as multiplexed DNA fluorescence in situ hybridization, and with tissue-specific genome profiling by DNA adenine methylation identification has increased our knowledge about the immense complexity and dynamics of the nucleus.

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.

Navigating the crowd: visualizing coordination between genome dynamics, structure, and transcription

The eukaryotic genome is hierarchically structured yet highly dynamic. Regulating transcription in this environment demands a high level of coordination to permit many proteins to interact with chromatin fiber at appropriate sites in a timely manner. We describe how recent advances in quantitative imaging techniques overcome caveats of sequencing-based methods (Hi-C and related) by enabling direct visualization of transcription factors and chromatin at high resolution, from single genes to the whole nucleus. We discuss the contribution of fluorescence imaging to deciphering the principles underlying this coordination within the crowded nuclear space in living cells and discuss challenges ahead.

Chromatin Remodelers in the 3D Nuclear Compartment

Chromatin remodeling complexes (CRCs) use ATP hydrolysis to maintain correct expression profiles, chromatin stability, and inherited epigenetic states. More than 20 CRCs have been described to date, which encompass four large families defined by their ATPase subunits. These complexes and their subunits are conserved from yeast to humans through evolution. Their activities depend on their catalytic subunits which through ATP hydrolysis provide the energy necessary to fulfill cellular functions such as gene transcription, DNA repair, and transposon silencing. These activities take place at the first levels of chromatin compaction, and CRCs have been recognized as essential elements of chromatin dynamics. Recent studies have demonstrated an important role for these complexes in the maintenance of higher order chromatin structure. In this review, we present an overview of the organization of the genome within the cell nucleus, the different levels of chromatin compaction, and importance of the architectural proteins, and discuss the role of CRCs and how their functions contribute to the dynamics of the 3D genome organization.

The Interchromatin Compartment Participates in the Structural and Functional Organization of the Cell Nucleus

This article focuses on the role of the interchromatin compartment (IC) in shaping nuclear landscapes. The IC is connected with nuclear pore complexes (NPCs) and harbors splicing speckles and nuclear bodies. It is postulated that the IC provides routes for imported transcription factors to target sites, for export routes of mRNA as ribonucleoproteins toward NPCs, as well as for the intranuclear passage of regulatory RNAs from sites of transcription to remote functional sites (IC hypothesis). IC channels are lined by less-compacted euchromatin, called the perichromatin region (PR). The PR and IC together form the active nuclear compartment (ANC). The ANC is co-aligned with the inactive nuclear compartment (INC), comprising more compacted heterochromatin. It is postulated that the INC is accessible for individual transcription factors, but inaccessible for larger macromolecular aggregates (limited accessibility hypothesis). This functional nuclear organization depends on still unexplored movements of genes and regulatory sequences between the two compartments.

Characterizing chromatin folding coordinate and landscape with deep learning

Genome organization is critical for setting up the spatial environment of gene transcription, and substantial progress has been made towards its high-resolution characterization. The underlying molecular mechanism for its establishment is much less understood. We applied a deep-learning approach, variational autoencoder (VAE), to analyze the fluctuation and heterogeneity of chromatin structures revealed by single-cell imaging and to identify a reaction coordinate for chromatin folding. This coordinate connects the seemingly random structures observed in individual cohesin-depleted cells as intermediate states along a folding pathway that leads to the formation of topologically associating domains (TAD). We showed that folding into wild-type-like structures remain energetically favorable in cohesin-depleted cells, potentially as a result of the phase separation between the two chromatin segments with active and repressive histone marks. The energetic stabilization, however, is not strong enough to overcome the entropic penalty, leading to the formation of only partially folded structures and the disappearance of TADs from contact maps upon averaging. Our study suggests that machine learning techniques, when combined with rigorous statistical mechanical analysis, are powerful tools for analyzing structural ensembles of chromatin.

Hi-C chromosome conformation capture sequencing of avian genomes using the BGISEQ-500 platform

BACKGROUND

Hi-C experiments couple DNA-DNA proximity with next-generation sequencing to yield an unbiased description of genome-wide interactions. Previous methods describing Hi-C experiments have focused on the industry-standard Illumina sequencing. With new next-generation sequencing platforms such as BGISEQ-500 becoming more widely available, protocol adaptations to fit platform-specific requirements are useful to give increased choice to researchers who routinely generate sequencing data.

RESULTS

We describe an in situ Hi-C protocol adapted to be compatible with the BGISEQ-500 high-throughput sequencing platform. Using zebra finch (Taeniopygia guttata) as a biological sample, we demonstrate how Hi-C libraries can be constructed to generate informative data using the BGISEQ-500 platform, following circularization and DNA nanoball generation. Our protocol is a modification of an Illumina-compatible method, based around blunt-end ligations in library construction, using un-barcoded, distally overhanging double-stranded adapters, followed by amplification using indexed primers. The resulting libraries are ready for circularization and subsequent sequencing on the BGISEQ series of platforms and yield data similar to what can be expected using Illumina-compatible approaches.

CONCLUSIONS

Our straightforward modification to an Illumina-compatible in situHi-C protocol enables data generation on the BGISEQ series of platforms, thus expanding the options available for researchers who wish to utilize the powerful Hi-C techniques in their research.

Functional annotation of human long noncoding RNAs via molecular phenotyping

Long noncoding RNAs (lncRNAs) constitute the majority of transcripts in the mammalian genomes, and yet, their functions remain largely unknown. As part of the FANTOM6 project, we systematically knocked down the expression of 285 lncRNAs in human dermal fibroblasts and quantified cellular growth, morphological changes, and transcriptomic responses using Capped Analysis of Gene Expression (CAGE). Antisense oligonucleotides targeting the same lncRNAs exhibited global concordance, and the molecular phenotype, measured by CAGE, recapitulated the observed cellular phenotypes while providing additional insights on the affected genes and pathways. Here, we disseminate the largest-to-date lncRNA knockdown data set with molecular phenotyping (over 1000 CAGE deep-sequencing libraries) for further exploration and highlight functional roles for ZNF213-AS1 and lnc-KHDC3L-2.

The contribution of non-coding regulatory elements to cardiovascular disease

Cardiovascular disease collectively accounts for a quarter of deaths worldwide. Genome-wide association studies across a range of cardiovascular traits and pathologies have highlighted the prevalence of common non-coding genetic variants within candidate loci. Here, we review genetic, epigenomic and molecular approaches to investigate the contribution of non-coding regulatory elements in cardiovascular biology. We then discuss recent insights on the emerging role of non-coding variation in predisposition to cardiovascular disease, with a focus on novel mechanistic examples from functional genomics studies. Lastly, we consider the clinical significance of these findings at present, and some of the current challenges facing the field.

Programmable Chromosome Painting with Oligopaints

Current methods for chromosome painting via fluorescence in situ hybridization (FISH) are costly, time-consuming, and limited in complexity. In contrast to conventional sources of probe, Oligopaints are computationally designed, synthesized on microarrays, and amplified by PCR. This approach allows for precise control over the sequences they target, which can range from a few kilobases to entire chromosomes with the same basic protocol. We have utilized the flexibility and scalability of Oligopaints to generate low-cost and renewable chromosome paints for Drosophila, mouse, and human chromosomes. These Oligopaint libraries can be customized to label any genomic feature(s) in a chromosome-wide manner. Additionally, this method is compatible with sequential FISH to label entire genomes with a single denaturation step. Here, we outline a protocol and considerations to scale the Oligopaint technology for fluorescent labeling of whole chromosomes.

HIFI: estimating DNA-DNA interaction frequency from Hi-C data at restriction-fragment resolution

Hi-C is a popular technique to map three-dimensional chromosome conformation. In principle, Hi-C's resolution is only limited by the size of restriction fragments. However, insufficient sequencing depth forces researchers to artificially reduce the resolution of Hi-C matrices at a loss of biological interpretability. We present the Hi-C Interaction Frequency Inference (HIFI) algorithms that accurately estimate restriction-fragment resolution Hi-C matrices by exploiting dependencies between neighboring fragments. Cross-validation experiments and comparisons to 5C data and known regulatory interactions demonstrate HIFI's superiority to existing approaches. In addition, HIFI's restriction-fragment resolution reveals a new role for active regulatory regions in structuring topologically associating domains.

Centric Fusions behind the Karyotype Evolution of Neotropical Nannostomus Pencilfishes (Characiforme, Lebiasinidae): First Insights from a Molecular Cytogenetic Perspective

Lebiasinidae is a Neotropical freshwater family widely distributed throughout South and Central America. Due to their often very small body size, Lebiasinidae species are cytogenetically challenging and hence largely underexplored. However, the available but limited karyotype data already suggested a high interspecific variability in the diploid chromosome number (2n), which is pronounced in the speciose genus Nannostomus, a popular taxon in ornamental fish trade due to its remarkable body coloration. Aiming to more deeply examine the karyotype diversification in Nannostomus, we combined conventional cytogenetics (Giemsa-staining and C-banding) with the chromosomal mapping of tandemly repeated 5S and 18S rDNA clusters and with interspecific comparative genomic hybridization (CGH) to investigate genomes of four representative Nannostomus species: N. beckfordi, N. eques, N. marginatus, and N. unifasciatus. Our data showed a remarkable variability in 2n, ranging from 2n = 22 in N. unifasciatus (karyotype composed exclusively of metacentrics/submetacentrics) to 2n = 44 in N. beckfordi (karyotype composed entirely of acrocentrics). On the other hand, patterns of 18S and 5S rDNA distribution in the analyzed karyotypes remained rather conservative, with only two 18S and two to four 5S rDNA sites. In view of the mostly unchanged number of chromosome arms (FN = 44) in all but one species (N. eques; FN = 36), and with respect to the current phylogenetic hypothesis, we propose Robertsonian translocations to be a significant contributor to the karyotype differentiation in (at least herein studied) Nannostomus species. Interspecific comparative genome hybridization (CGH) using whole genomic DNAs mapped against the chromosome background of N. beckfordi found a moderate divergence in the repetitive DNA content among the species’ genomes. Collectively, our data suggest that the karyotype differentiation in Nannostomus has been largely driven by major structural rearrangements, accompanied by only low to moderate dynamics of repetitive DNA at the sub-chromosomal level. Possible mechanisms and factors behind the elevated tolerance to such a rate of karyotype change in Nannostomus are discussed.

An RB-Condensin II Complex Mediates Long-Range Chromosome Interactions and Influences Expression at Divergently Paired Genes

Interphase chromosomes are organized into topologically associated domains in order to establish and maintain integrity of transcriptional programs that remain poorly understood. Here, we show that condensin II and TFIIIC are recruited to bidirectionally transcribed promoters by a mechanism that is dependent on the retinoblastoma (RB) protein. Long-range chromosome contacts are disrupted by loss of condensin II loading, which leads to altered expression at bidirectional gene pairs. This study demonstrates that mammalian condensin II functions to organize long-range chromosome contacts and regulate transcription at specific genes. In addition, RB dependence of condensin II suggests that widespread misregulation of chromosome contacts and transcriptional alterations are a consequence of RB mutation.

In silico prediction of high-resolution Hi-C interaction matrices

The three-dimensional (3D) organization of the genome plays an important role in gene regulation bringing distal sequence elements in 3D proximity to genes hundreds of kilobases away. Hi-C is a powerful genome-wide technique to study 3D genome organization. Owing to experimental costs, high resolution Hi-C datasets are limited to a few cell lines. Computational prediction of Hi-C counts can offer a scalable and inexpensive approach to examine 3D genome organization across multiple cellular contexts. Here we present HiC-Reg, an approach to predict contact counts from one-dimensional regulatory signals. HiC-Reg predictions identify topologically associating domains and significant interactions that are enriched for CCCTC-binding factor (CTCF) bidirectional motifs and interactions identified from complementary sources. CTCF and chromatin marks, especially repressive and elongation marks, are most important for HiC-Reg’s predictive performance. Taken together, HiC-Reg provides a powerful framework to generate high-resolution profiles of contact counts that can be used to study individual locus level interactions and higher-order organizational units of the genome.

Existing computational approaches to predict long-range regulatory interactions do not fully exploit high-resolution Hi-C datasets. Here the authors present a Random Forests regression-based approach to predict high-resolution Hi-C counts using one-dimensional regulatory genomic signals.