Imputation of 3D genome structure by genetic-epigenetic interaction modeling in mice

Gene expression is known to be affected by interactions between local genetic variation and DNA accessibility, with the latter organized into three-dimensional chromatin structures. Analyses of these interactions have previously been limited, obscuring their regulatory context, and the extent to which they occur throughout the genome. Here, we undertake a genome-scale analysis of these interactions in a genetically diverse population to systematically identify global genetic-epigenetic interaction, and reveal constraints imposed by chromatin structure. We establish the extent and structure of genotype-by-epigenotype interaction using embryonic stem cells derived from Diversity Outbred mice. This mouse population segregates millions of variants from eight inbred founders, enabling precision genetic mapping with extensive genotypic and phenotypic diversity. With 176 samples profiled for genotype, gene expression, and open chromatin, we used regression modeling to infer genetic-epigenetic interactions on a genome-wide scale. Our results demonstrate that statistical interactions between genetic variants and chromatin accessibility are common throughout the genome. We found that these interactions occur within the local area of the affected gene, and that this locality corresponds to topologically associated domains (TADs). The likelihood of interaction was most strongly defined by the three-dimensional (3D) domain structure rather than linear DNA sequence. We show that stable 3D genome structure is an effective tool to guide searches for regulatory elements and, conversely, that regulatory elements in genetically diverse populations provide a means to infer 3D genome structure. We confirmed this finding with CTCF ChIP-seq that revealed strain-specific binding in the inbred founder mice. In stem cells, open chromatin participating in the most significant regression models demonstrated an enrichment for developmental genes and the TAD-forming CTCF-binding complex, providing an opportunity for statistical inference of shifting TAD boundaries operating during early development. These findings provide evidence that genetic and epigenetic factors operate within the context of 3D chromatin structure.

Spatially coordinated heterochromatinization of long synaptic genes in fragile X syndrome

Short tandem repeat (STR) instability causes transcriptional silencing in several repeat expansion disorders. In fragile X syndrome (FXS), mutation-length expansion of a CGG STR represses FMR1 via local DNA methylation. Here, we find megabase-scale H3K9me3 domains on autosomes and encompassing FMR1 on the X chromosome in FXS patient-derived iPSCs, iPSC-derived neural progenitors, EBV-transformed lymphoblasts, and brain tissue with mutation-length CGG expansion. H3K9me3 domains connect via inter-chromosomal interactions and demarcate severe misfolding of TADs and loops. They harbor long synaptic genes replicating at the end of S phase, replication-stress-induced double-strand breaks, and STRs prone to stepwise somatic instability. CRISPR engineering of the mutation-length CGG to premutation length reverses H3K9me3 on the X chromosome and multiple autosomes, refolds TADs, and restores gene expression. H3K9me3 domains can also arise in normal-length iPSCs created with perturbations linked to genome instability, suggesting their relevance beyond FXS. Our results reveal Mb-scale heterochromatinization and trans interactions among loci susceptible to instability.

Cis-regulatory effect of HPV integration is constrained by host chromatin architecture in cervical cancers

Human papillomavirus (HPV) infections are the primary drivers of cervical cancers, and often HPV DNA gets integrated into the host genome. Although the oncogenic impact of HPV encoded genes is relatively well known, the cis-regulatory effect of integrated HPV DNA on host chromatin structure and gene regulation remains less understood. We investigated genome-wide patterns of HPV integrations and associated host gene expression changes in the context of host chromatin states and topologically associating domains (TADs). HPV integrations were significantly enriched in active chromatin regions and depleted in inactive ones. Interestingly, regardless of chromatin state, genomic regions flanking HPV integrations showed transcriptional upregulation. Nevertheless, upregulation (both local and long-range) was mostly confined to TADs with integration, but not affecting adjacent TADs. Few TADs showed recurrent integrations associated with overexpression of oncogenes within them (e.g. MYC, PVT1, TP63 and ERBB2) regardless of proximity. Hi-C and 4C-seq analyses in cervical cancer cell line (HeLa) demonstrated chromatin looping interactions between integrated HPV and MYC/PVT1 regions (~ 500 kb apart), leading to allele-specific overexpression. Based on these, we propose HPV integrations can trigger multimodal oncogenic activation to promote cancer progression.

Genetic variation and domestication of horses revealed by 10 chromosome-level genomes and whole-genome resequencing

Understanding the genetic variations of the horse (Equus caballus) genome will improve breeding conservation and welfare. However, genetic variations in long segments, such as structural variants (SVs), remain understudied. We de novo assembled 10 chromosome-level three-dimensional horse genomes, each representing a distinct breed, and analysed horse SVs using a multi-assembly approach. Our findings suggest that SVs with the accumulation of mammalian-wide interspersed repeats related to long interspersed nuclear elements might be a horse-specific mechanism to modulate genome-wide gene regulatory networks. We found that olfactory receptors were commonly loss and accumulated deleterious mutations, but no purge of deleterious mutations occurred during horse domestication. We examined the potential effects of SVs on the spatial structure of chromatin via topologically associating domains (TADs). Breed-specific TADs were significantly enriched by breed-specific SVs. We identified 4199 unique breakpoint-resolved novel insertions across all chromosomes that account for 2.84 Mb sequences missing from the reference genome. Several novel insertions might have potential functional consequences, as 519 appeared to reside within 449 gene bodies. These genes are primarily involved in pathogen recognition, innate immune responses and drug metabolism. Moreover, 37 diverse horses were resequenced. Combining this with public data, we analysed 97 horses through a comparative population genomics approach to identify the genetic basis underlying breed characteristics using Thoroughbreds as a case study. We provide new scientific evidence for horse domestication, an understanding of the genetic mechanism underlying the phenotypic evolution of horses, and a comprehensive genetic variation resource for further genetic studies of horses.

Revisiting Assessment of Computational Methods for Hi-C Data Analysis

The performances of algorithms for Hi-C data preprocessing, the identification of topologically associating domains, and the detection of chromatin interactions and promoter-enhancer interactions have been mostly evaluated using semi-quantitative or synthetic data approaches, without utilizing the most recent methods, since 2017. In this study, we comprehensively evaluated 24 popular state-of-the-art methods for the complete end-to-end pipeline of Hi-C data analysis, using manually curated or experimentally validated benchmark datasets, including a CRISPR dataset for promoter-enhancer interaction validation. Our results indicate that, although no single method exhibited superior performance in all situations, HiC-Pro, DomainCaller, and Fit-Hi-C2 showed relatively balanced performances of most evaluation metrics for preprocessing, topologically associating domain identification, and chromatin interaction/promoter-enhancer interaction detection, respectively. The comprehensive comparison presented in this manuscript provides a reference for researchers to choose Hi-C analysis tools that best suit their needs.

Glucocorticoid stimulation induces regionalized gene responses within topologically associating domains

Transcription-factor binding to cis-regulatory regions regulates the gene expression program of a cell, but occupancy is often a poor predictor of the gene response. Here, we show that glucocorticoid stimulation led to the reorganization of transcriptional coregulators MED1 and BRD4 within topologically associating domains (TADs), resulting in active or repressive gene environments. Indeed, we observed a bias toward the activation or repression of a TAD when their activities were defined by the number of regions gaining and losing MED1 and BRD4 following dexamethasone (Dex) stimulation. Variations in Dex-responsive genes at the RNA levels were consistent with the redistribution of MED1 and BRD4 at the associated cis-regulatory regions. Interestingly, Dex-responsive genes without the differential recruitment of MED1 and BRD4 or binding by the glucocorticoid receptor were found within TADs, which gained or lost MED1 and BRD4, suggesting a role of the surrounding environment in gene regulation. However, the amplitude of the response of Dex-regulated genes was higher when the differential recruitment of the glucocorticoid receptor and transcriptional coregulators was observed, reaffirming the role of transcription factor-driven gene regulation and attributing a lesser role to the TAD environment. These results support a model where a signal-induced transcription factor induces a regionalized effect throughout the TAD, redefining the notion of direct and indirect effects of transcription factors on target genes.

CATAD: exploring topologically associating domains from an insight of core-attachment structure

Identifying topologically associating domains (TADs), which are considered as the basic units of chromosome structure and function, can facilitate the exploration of the 3D-structure of chromosomes. Methods have been proposed to identify TADs by detecting the boundaries of TADs or identifying the closely interacted regions as TADs, while the possible inner structure of TADs is seldom investigated. In this study, we assume that a TAD is composed of a core and its surrounding attachments, and propose a method, named CATAD, to identify TADs based on the core-attachment structure model. In CATAD, the cores of TADs are identified based on the local density and cosine similarity, and the surrounding attachments are determined based on boundary insulation. CATAD was applied to the Hi-C data of two human cell lines and two mouse cell lines, and the results show that the boundaries of TADs identified by CATAD are significantly enriched by structural proteins, histone modifications, transcription start sites and enzymes. Furthermore, CATAD outperforms other methods in many cases, in terms of the average peak, boundary tagged ratio and fold change. In addition, CATAD is robust and rarely affected by the different resolutions of Hi-C matrices. Conclusively, identifying TADs based on the core-attachment structure is useful, which may inspire researchers to explore TADs from the angles of possible spatial structures and formation process.

Chromatin modules and their implication in genomic organization and gene regulation

Regulation of gene expression is a complex but highly guided process. While genomic technologies and computational approaches have allowed high-throughput mapping of cis-regulatory elements (CREs) and their interactions in 3D, their precise role in regulating gene expression remains obscure. Recent complementary observations revealed that interactions between CREs frequently result in the formation of small-scale functional modules within topologically associating domains. Such chromatin modules likely emerge from a complex interplay between regulatory machineries assembled at CREs, including site-specific binding of transcription factors. Here, we review the methods that allow identifying chromatin modules, summarize possible mechanisms that steer CRE interactions within these modules, and discuss outstanding challenges to uncover how chromatin modules fit in our current understanding of the functional 3D genome.

Mode and Tempo of 3D Genome Evolution in Drosophila

Topologically associating domains (TADs) are thought to play an important role in preventing gene misexpression by spatially constraining enhancer-promoter contacts. The deleterious nature of gene misexpression implies that TADs should, therefore, be conserved among related species. Several early studies comparing chromosome conformation between species reported high levels of TAD conservation; however, more recent studies have questioned these results. Furthermore, recent work suggests that TAD reorganization is not associated with extensive changes in gene expression. Here, we investigate the evolutionary conservation of TADs among 11 species of Drosophila. We use Hi-C data to identify TADs in each species and employ a comparative phylogenetic approach to derive empirical estimates of the rate of TAD evolution. Surprisingly, we find that TADs evolve rapidly. However, we also find that the rate of evolution depends on the chromatin state of the TAD, with TADs enriched for developmentally regulated chromatin evolving significantly slower than TADs enriched for broadly expressed, active chromatin. We also find that, after controlling for differences in chromatin state, highly conserved TADs do not exhibit higher levels of gene expression constraint. These results suggest that, in general, most TADs evolve rapidly and their divergence is not associated with widespread changes in gene expression. However, higher levels of evolutionary conservation and gene expression constraints in TADs enriched for developmentally regulated chromatin suggest that these TAD subtypes may be more important for regulating gene expression, likely due to the larger number of long-distance enhancer-promoter contacts associated with developmental genes.

Evolution of the ancestral mammalian karyotype and syntenic regions

Computational reconstruction of ancestral mammalian karyotypes revealed a comprehensive picture of the chromosome rearrangements that occurred over the evolutionary history of mammals. Ancient gene order, in some cases extending to full chromosomes, was found conserved for more than 300 My, demonstrating strong evolutionary constraint against rearrangements in some regions. Conserved segments of chromosomes are enriched for genes that control developmental processes. Therefore, Darwinian selection likely maintains ancient gene combinations while allowing for genomic innovations within or near chromosomal sites that break and rearrange over evolutionary time. The revealed relationship between the three-dimensional structure of chromosomes and the evolutionary stability of chromosome segments provides additional insights into the mechanisms of chromosome evolution and diseases associated with genome rearrangements.

Decrypting the rearrangements that drive mammalian chromosome evolution is critical to understanding the molecular bases of speciation, adaptation, and disease susceptibility. Using 8 scaffolded and 26 chromosome-scale genome assemblies representing 23/26 mammal orders, we computationally reconstructed ancestral karyotypes and syntenic relationships at 16 nodes along the mammalian phylogeny. Three different reference genomes (human, sloth, and cattle) representing phylogenetically distinct mammalian superorders were used to assess reference bias in the reconstructed ancestral karyotypes and to expand the number of clades with reconstructed genomes. The mammalian ancestor likely had 19 pairs of autosomes, with nine of the smallest chromosomes shared with the common ancestor of all amniotes (three still conserved in extant mammals), demonstrating a striking conservation of synteny for ∼320 My of vertebrate evolution. The numbers and types of chromosome rearrangements were classified for transitions between the ancestral mammalian karyotype, descendent ancestors, and extant species. For example, 94 inversions, 16 fissions, and 14 fusions that occurred over 53 My differentiated the therian from the descendent eutherian ancestor. The highest breakpoint rate was observed between the mammalian and therian ancestors (3.9 breakpoints/My). Reconstructed mammalian ancestor chromosomes were found to have distinct evolutionary histories reflected in their rates and types of rearrangements. The distributions of genes, repetitive elements, topologically associating domains, and actively transcribed regions in multispecies homologous synteny blocks and evolutionary breakpoint regions indicate that purifying selection acted over millions of years of vertebrate evolution to maintain syntenic relationships of developmentally important genes and regulatory landscapes of gene-dense chromosomes.

Repression and 3D-restructuring resolves regulatory conflicts in evolutionarily rearranged genomes

Regulatory landscapes drive complex developmental gene expression, but it remains unclear how their integrity is maintained when incorporating novel genes and functions during evolution. Here, we investigated how a placental mammal-specific gene, Zfp42, emerged in an ancient vertebrate topologically associated domain (TAD) without adopting or disrupting the conserved expression of its gene, Fat1. In ESCs, physical TAD partitioning separates Zfp42 and Fat1 with distinct local enhancers that drive their independent expression. This separation is driven by chromatin activity and not CTCF/cohesin. In contrast, in embryonic limbs, inactive Zfp42 shares Fat1's intact TAD without responding to active Fat1 enhancers. However, neither Fat1 enhancer-incompatibility nor nuclear envelope-attachment account for Zfp42's unresponsiveness. Rather, Zfp42's promoter is rendered inert to enhancers by context-dependent DNA methylation. Thus, diverse mechanisms enabled the integration of independent Zfp42 regulation in the Fat1 locus. Critically, such regulatory complexity appears common in evolution as, genome wide, most TADs contain multiple independently expressed genes.

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.

Opportunity to improve livestock traits using 3D genomics

The advent of high-throughput chromosome conformation capture and sequencing (Hi-C) has enabled researchers to probe the 3D architecture of the mammalian genome in a genome-wide manner. Simultaneously, advances in epigenomic assays, such as chromatin immunoprecipitation and sequencing (ChIP-seq) and DNase-seq, have enabled researchers to study cis-regulatory interactions and chromatin accessibility across the same genome-wide scale. The use of these data has revealed many unique insights into gene regulation and disease pathomechanisms in several model organisms. With the advent of these high-throughput sequencing technologies, there has been an ever-increasing number of datasets available for study; however, this is often limited to model organisms. Livestock species play critical roles in the economies of developing and developed nations alike. Despite this, they are greatly underrepresented in the 3D genomics space; Hi-C and related technologies have the potential to revolutionise livestock breeding by enabling a more comprehensive understanding of how production traits are controlled. The growth in human and model organism Hi-C data has seen a surge in the availability of computational tools for use in 3D genomics, with some tools using machine learning techniques to predict features and improve dataset quality. In this review, we provide an overview of the 3D genome and discuss the status of 3D genomics in livestock before delving into advancing the field by drawing inspiration from research in human and mouse. We end by offering future directions for livestock research in the field of 3D genomics.

Structural Variations of the 3D Genome Architecture in Cervical Cancer Development

Human papillomavirus (HPV) integration is the major contributor to cervical cancer (CC) development by inducing structural variations (SVs) in the human genome. SVs are directly associated with the three-dimensional (3D) genome structure leading to cancer development. The detection of SVs is not a trivial task, and several genome-wide techniques have greatly helped in the identification of SVs in the cancerous genome. However, in cervical cancer, precise prediction of SVs mainly translocations and their effects on 3D-genome and gene expression still need to be explored. Here, we have used high-throughput chromosome conformation capture (Hi-C) data of cervical cancer to detect the SVs, especially the translocations, and validated it through whole-genome sequencing (WGS) data. We found that the cervical cancer 3D-genome architecture rearranges itself as compared to that in the normal tissue, and 24% of the total genome switches their A/B compartments. Moreover, translocation detection from Hi-C data showed the presence of high-resolution t(4;7) (q13.1; q31.32) and t(1;16) (q21.2; q22.1) translocations, which disrupted the expression of the genes located at and nearby positions. Enrichment analysis suggested that the disrupted genes were mainly involved in controlling cervical cancer-related pathways. In summary, we detect the novel SVs through Hi-C data and unfold the association among genome-reorganization, translocations, and gene expression regulation. The results help understand the underlying pathogenicity mechanism of SVs in cervical cancer development and identify the targeted therapeutics against cervical cancer.

3D genome organisation in Drosophila

Ever since Thomas Hunt Morgan's discovery of the chromosomal basis of inheritance by using Drosophila melanogaster as a model organism, the fruit fly has remained an essential model system in studies of genome biology, including chromatin organisation. Very much as in vertebrates, in Drosophila, the genome is organised in territories, compartments and topologically associating domains (TADs). However, these domains might be formed through a slightly different mechanism than in vertebrates due to the presence of a large and potentially redundant set of insulator proteins and the minor role of dCTCF in TAD boundary formation. Here, we review the different levels of chromatin organisation in Drosophila and discuss mechanisms and factors that might be involved in TAD formation. The dynamics of TADs and enhancer-promoter interactions in the context of transcription are covered in the light of currently conflicting results. Finally, we illustrate the value of polymer modelling approaches to infer the principles governing the three-dimensional organisation of the Drosophila genome.

Topologically Associating Domains and Regulatory Landscapes in Development, Evolution and Disease

Animal genomes are folded in topologically associating domains (TADs) that have been linked to the regulation of the genes they contain by constraining regulatory interactions between cis-regulatory elements and promoters. Therefore, TADs are proposed as structural scaffolds for the establishment of regulatory landscapes (RLs). In this review, we discuss recent advances in the connection between TADs and gene regulation, their relationship with gene RLs and their dynamics during development and differentiation. Moreover, we describe how restructuring TADs may lead to pathological conditions, which explains their high evolutionary conservation, but at the same time it provides a substrate for the emergence of evolutionary innovations that lay at the origin of vertebrates and other phylogenetic clades.

Utilization of Host Cell Chromosome Conformation by Viral Pathogens: Knowing When to Hold and When to Fold

Though viruses have their own genomes, many depend on the nuclear environment of their hosts for replication and survival. A substantial body of work has therefore been devoted to understanding how viral and eukaryotic genomes interact. Recent advances in chromosome conformation capture technologies have provided unprecedented opportunities to visualize how mammalian genomes are organized and, by extension, how packaging of nuclear DNA impacts cellular processes. Recent studies have indicated that some viruses, upon entry into host cell nuclei, produce factors that alter host chromatin topology, and thus, impact the 3D organization of the host genome. Additionally, a variety of distinct viruses utilize host genome architectural factors to advance various aspects of their life cycles. Indeed, human gammaherpesviruses, known for establishing long-term reservoirs of latent infection in B lymphocytes, utilize 3D principles of genome folding to package their DNA and establish latency in host cells. This manipulation of host epigenetic machinery by latent viral genomes is etiologically linked to the onset of B cell oncogenesis. Small DNA viruses, by contrast, are tethered to distinct cellular sites that support virus expression and replication. Here, we briefly review the recent findings on how viruses and host genomes spatially communicate, and how this impacts virus-induced pathology.

A TAD Skeptic: Is 3D Genome Topology Conserved?

The notion that topologically associating domains (TADs) are highly conserved across species is prevalent in the field of 3D genomics. However, what exactly is meant by 'highly conserved' and what are the actual comparative data that support this notion? To address these questions, we performed a historical review of the relevant literature and retraced numerous citation chains to reveal the primary data that were used as the basis for the widely accepted conclusion that TADs are highly conserved across evolution. A thorough review of the available evidence suggests the answer may be more complex than what is commonly presented.

Epigenetics of Muscle- and Brain-Specific Expression of KLHL Family Genes

KLHL and the related KBTBD genes encode components of the Cullin-E3 ubiquitin ligase complex and typically target tissue-specific proteins for degradation, thereby affecting differentiation, homeostasis, metabolism, cell signaling, and the oxidative stress response. Despite their importance in cell function and disease (especially, KLHL40, KLHL41, KBTBD13, KEAP1, and ENC1), previous studies of epigenetic factors that affect transcription were predominantly limited to promoter DNA methylation. Using diverse tissue and cell culture whole-genome profiles, we examined 17 KLHL or KBTBD genes preferentially expressed in skeletal muscle or brain to identify tissue-specific enhancer and promoter chromatin, open chromatin (DNaseI hypersensitivity), and DNA hypomethylation. Sixteen of the 17 genes displayed muscle- or brain-specific enhancer chromatin in their gene bodies, and most exhibited specific intergenic enhancer chromatin as well. Seven genes were embedded in super-enhancers (particularly strong, tissue-specific clusters of enhancers). The enhancer chromatin regions typically displayed foci of DNA hypomethylation at peaks of open chromatin. In addition, we found evidence for an intragenic enhancer in one gene upregulating expression of its neighboring gene, specifically for KLHL40/HHATL and KLHL38/FBXO32 gene pairs. Many KLHL/KBTBD genes had tissue-specific promoter chromatin at their 5' ends, but surprisingly, two (KBTBD11 and KLHL31) had constitutively unmethylated promoter chromatin in their 3' exons that overlaps a retrotransposed KLHL gene. Our findings demonstrate the importance of expanding epigenetic analyses beyond the 5' ends of genes in studies of normal and abnormal gene regulation.

Dysregulation of chromatin organization in pediatric and adult brain tumors: oncoepigenomic contributions to tumorigenesis and cancer stem cell properties

The three-dimensional (3D) organization of the genome is a crucial enabler of cell fate, identity, and function. In this review, we will focus on the emerging role of altered 3D genome organization in the etiology of disease, with a special emphasis on brain cancers. We discuss how different genetic alterations can converge to disrupt the epigenome in childhood and adult brain tumors, by causing aberrant DNA methylation and by affecting the amounts and genomic distribution of histone post-translational modifications. We also highlight examples that illustrate how epigenomic alterations have the potential to affect 3D genome architecture in brain tumors. Finally, we will propose the concept of "epigenomic erosion" to explain the transition from stem-like cells to differentiated cells in hierarchically organized brain cancers.

Transcription of carbonyl reductase 1 is regulated by DNA topoisomerase II beta

DNA topoisomerase II beta (TOP2B) has a role in transcriptional regulation. Here, to further investigate transcriptional regulation by TOP2B, we used RNA-sequencing and real-time PCR to analyse the differential gene expression profiles of wild-type and two independent TOP2B-null pre-B Nalm-6 cell lines, one generated by targeted insertion and the other using CRISPR-Cas9 gene editing. We identified carbonyl reductase 1 (CBR1) among the most significantly downregulated genes in these TOP2B-null cells. Reduced CBR1 expression was accompanied by loss of binding of the transcription factors USF2 and MAX to the CBR1 promoter. We describe possible mechanisms by which loss of TOP2B results in CBR1 downregulation. To our knowledge, this is the first report of a link between TOP2B and CBR1.

Hi-C Identifies Complex Genomic Rearrangements and TAD-Shuffling in Developmental Diseases

Genome-wide analysis methods, such as array comparative genomic hybridization (CGH) and whole-genome sequencing (WGS), have greatly advanced the identification of structural variants (SVs) in the human genome. However, even with standard high-throughput sequencing techniques, complex rearrangements with multiple breakpoints are often difficult to resolve, and predicting their effects on gene expression and phenotype remains a challenge. Here, we address these problems by using high-throughput chromosome conformation capture (Hi-C) generated from cultured cells of nine individuals with developmental disorders (DDs). Three individuals had previously been identified as harboring duplications at the SOX9 locus and six had been identified with translocations. Hi-C resolved the positions of the duplications and was instructive in interpreting their distinct pathogenic effects, including the formation of new topologically associating domains (neo-TADs). Hi-C was very sensitive in detecting translocations, and it revealed previously unrecognized complex rearrangements at the breakpoints. In several cases, we observed the formation of fused-TADs promoting ectopic enhancer-promoter interactions that were likely to be involved in the disease pathology. In summary, we show that Hi-C is a sensible method for the detection of complex SVs in a clinical setting. The results help interpret the possible pathogenic effects of the SVs in individuals with DDs.

Global Regulatory DNA Potentiation by SMARCA4 Propagates to Selective Gene Expression Programs via Domain-Level Remodeling

The human genome encodes millions of regulatory elements, of which only a small fraction are active within a given cell type. Little is known about the global impact of chromatin remodelers on regulatory DNA landscapes and how this translates to gene expression. We use precision genome engineering to reawaken homozygously inactivated SMARCA4, a central ATPase of the human SWI/SNF chromatin remodeling complex, in lung adenocarcinoma cells. Here, we combine DNase I hypersensitivity, histone modification, and transcriptional profiling to show that SMARCA4 dramatically increases both the number and magnitude of accessible chromatin sites genome-wide, chiefly by unmasking sites of low regulatory factor occupancy. By contrast, transcriptional changes are concentrated within well-demarcated remodeling domains wherein expression of specific genes is gated by both distal element activation and promoter chromatin configuration. Our results provide a perspective on how global chromatin remodeling activity is translated to gene expression via regulatory DNA.

Supramolecular multilayer organization of chromosomes: possible functional roles of planar chromatin in gene expression and DNA replication and repair

Experimental evidence indicates that the chromatin filament is self-organized into a multilayer planar structure that is densely stacked in metaphase and unstacked in interphase. This chromatin organization is unexpected, but it is shown that diverse supramolecular assemblies, including dinoflagellate chromosomes, are multilayered. The mechanical strength of planar chromatin protects the genome integrity, even when double-strand breaks are produced. Here, it is hypothesized that the chromatin filament in the loops and topologically associating domains is folded within the thin layers of the multilaminar chromosomes. It is also proposed that multilayer chromatin has two states: inactive when layers are stacked and active when layers are unstacked. Importantly, the well-defined topology of planar chromatin may facilitate DNA replication without entanglements and DNA repair by homologous recombination.

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

The tridimensional (3D) organization of mammalian genomes combines structures from different length scales. Within this organization, Topologically Associating Domains (TADs) are visible in Hi-C heat maps at the sub-megabase scale. The integrity of TADs is important for correct gene expression, but in a context-dependent and variable manner. The correct structure and function of TADs require the binding of the CTCF protein at both borders, which appears to block an active and dynamic mechanism of "Cohesin-mediated loop extrusion." As a result, mammalian TADs appear as so-called "loop domains" in Hi-C data, which are the focus of this review. Here, we present a reanalysis of TADs from three "golden-standard" mammalian Hi-C data sets. Despite the prominent presence of TADs in Hi-C heat maps from all studies, we find consistently that regions within these domains are only moderately insulated from their surroundings. Moreover, single-cell Hi-C and superresolution microscopy have revealed that the structure of TADs and the position of their borders can vary from cell to cell. The function of TADs as units of gene regulation may thus require additional aspects, potentially incorporating the mechanism of loop extrusion as well. Recent developments in single-cell and multi-contact genomics and superresolution microscopy assays will be instrumental to link TAD formation and structure to their function in transcriptional regulation.

The Genomic Code: A Pervasive Encoding/Molding of Chromatin Structures and a Solution of the "Non-Coding DNA" Mystery

Recent investigations have revealed 1) that the isochores of the human genome group into two super-families characterized by two different long-range 3D structures, and 2) that these structures, essentially based on the distribution and topology of short sequences, mold primary chromatin domains (and define nucleosome binding). More specifically, GC-poor, gene-poor isochores are low-heterogeneity sequences with oligo-A spikes that mold the lamina-associated domains (LADs), whereas GC-rich, gene-rich isochores are characterized by single or multiple GC peaks that mold the topologically associating domains (TADs). The formation of these "primary TADs" may be followed by extrusion under the action of cohesin and CTCF. Finally, the genomic code, which is responsible for the pervasive encoding and molding of primary chromatin domains (LADs and primary TADs, namely the "gene spaces"/"spatial compartments") resolves the longstanding problems of "non-coding DNA," "junk DNA," and "selfish DNA" leading to a new vision of the genome as shaped by DNA sequences.

MeCP2 Represses Enhancers through Chromosome Topology-Associated DNA Methylation

The genomes of mammalian neurons contain uniquely high levels of non-CG DNA methylation that can be bound by the Rett syndrome protein, MeCP2, to regulate gene expression. How patterns of non-CG methylation are established in neurons and the mechanism by which this methylation works with MeCP2 to control gene expression is unclear. Here, we find that genes repressed by MeCP2 are often located within megabase-scale regions of high non-CG methylation that correspond with topologically associating domains of chromatin folding. MeCP2 represses enhancers found in these domains that are enriched for non-CG and CG methylation, with the strongest repression occurring for enhancers located within MeCP2-repressed genes. These alterations in enhancer activity provide a mechanism for how MeCP2 disruption in disease can lead to widespread changes in gene expression. Hence, we find that DNA topology can shape non-CG DNA methylation across the genome to dictate MeCP2-mediated enhancer regulation in the brain.

HiCeekR: A Novel Shiny App for Hi-C Data Analysis

The High-throughput Chromosome Conformation Capture (Hi-C) technique combines the power of the Next Generation Sequencing technologies with chromosome conformation capture approach to study the 3D chromatin organization at the genome-wide scale. Although such a technique is quite recent, many tools are already available for pre-processing and analyzing Hi-C data, allowing to identify chromatin loops, topological associating domains and A/B compartments. However, only a few of them provide an exhaustive analysis pipeline or allow to easily integrate and visualize other omic layers. Moreover, most of the available tools are designed for expert users, who have great confidence with command-line applications. In this paper, we present HiCeekR (https://github.com/lucidif/HiCeekR), a novel R Graphical User Interface (GUI) that allows researchers to easily perform a complete Hi-C data analysis. With the aid of the Shiny libraries, it integrates several R/Bioconductor packages for Hi-C data analysis and visualization, guiding the user during the entire process. Here, we describe its architecture and functionalities, then illustrate its capabilities using a publicly available dataset.

X Chromosome Domain Architecture Regulates Caenorhabditis elegans Lifespan but Not Dosage Compensation

Mechanisms establishing higher-order chromosome structures and their roles in gene regulation are elusive. We analyzed chromosome architecture during nematode X chromosome dosage compensation, which represses transcription via a dosage-compensation condensin complex (DCC) that binds hermaphrodite Xs and establishes megabase-sized topologically associating domains (TADs). We show that DCC binding at high-occupancy sites (rex sites) defines eight TAD boundaries. Single rex deletions disrupted boundaries, and single insertions created new boundaries, demonstrating that a rex site is necessary and sufficient to define DCC-dependent boundary locations. Deleting eight rex sites (8rexΔ) recapitulated TAD structure of DCC mutants, permitting analysis when chromosome-wide domain architecture was disrupted but most DCC binding remained. 8rexΔ animals exhibited no changes in X expression and lacked dosage-compensation mutant phenotypes. Hence, TAD boundaries are neither the cause nor the consequence of DCC-mediated gene repression. Abrogating TAD structure did, however, reduce thermotolerance, accelerate aging, and shorten lifespan, implicating chromosome architecture in stress responses and aging.

NETWORK MODELLING OF TOPOLOGICAL DOMAINS USING HI-C DATA

Chromosome conformation capture experiments such as Hi-C are used to map the three-dimensional spatial organization of genomes. One specific feature of the 3D organization is known as topologically associating domains (TADs), which are densely interacting, contiguous chromatin regions playing important roles in regulating gene expression. A few algorithms have been proposed to detect TADs. In particular, the structure of Hi-C data naturally inspires application of community detection methods. However, one of the drawbacks of community detection is that most methods take exchangeability of the nodes in the network for granted; whereas the nodes in this case, that is, the positions on the chromosomes, are not exchangeable. We propose a network model for detecting TADs using Hi-C data that takes into account this nonexchangeability. in addition, our model explicitly makes use of cell-type specific CTCF binding sites as biological covariates and can be used to identify conserved TADs across multiple cell types. The model leads to a likelihood objective that can be efficiently optimized via relaxation. We also prove that when suitably initialized, this model finds the underlying TAD structure with high probability. using simulated data, we show the advantages of our method and the caveats of popular community detection methods, such as spectral clustering, in this application. Applying our method to real Hi-C data, we demonstrate the domains identified have desirable epigenetic features and compare them across different cell types.

The Secret Life of Chromosome Loops upon DNA Double-Strand Break

DNA double-strand breaks (DSBs) are harmful lesions that severely challenge genomic integrity, and recent evidence suggests that DSBs occur more frequently on the genome than previously thought. These lesions activate a complex and multilayered response called the DNA damage response, which allows to coordinate their repair with the cell cycle progression. While the mechanistic details of repair processes have been narrowed, thanks to several decades of intense studies, our knowledge of the impact of DSB on chromatin composition and chromosome architecture is still very sparse. However, the recent development of various tools to induce DSB at annotated loci, compatible with next-generation sequencing-based approaches, is opening a new framework to tackle these questions. Here we discuss the influence of initial and DSB-induced chromatin conformation and the strong potential of 3C-based technologies to decipher the contribution of chromosome architecture during DSB repair.

Gain of CTCF-Anchored Chromatin Loops Marks the Exit from Naive Pluripotency

The genome of pluripotent stem cells adopts a unique three-dimensional architecture featuring weakly condensed heterochromatin and large nucleosome-free regions. Yet, it is unknown whether structural loops and contact domains display characteristics that distinguish embryonic stem cells (ESCs) from differentiated cell types. We used genome-wide chromosome conformation capture and super-resolution imaging to determine nuclear organization in mouse ESC and neural stem cell (NSC) derivatives. We found that loss of pluripotency is accompanied by widespread gain of structural loops. This general architectural change correlates with enhanced binding of CTCF and cohesins and more pronounced insulation of contacts across chromatin boundaries in lineage-committed cells. Reprogramming NSCs to pluripotency restores the unique features of ESC domain topology. Domains defined by the anchors of loops established upon differentiation are enriched for developmental genes. Chromatin loop formation is a pervasive structural alteration to the genome that accompanies exit from pluripotency and delineates the spatial segregation of developmentally regulated genes.

Molecular Lesions of Insulator CTCF and Its Paralogue CTCFL (BORIS) in Cancer: An Analysis from Published Genomic Studies

CTCF (CCCTC-binding factor) is a transcription regulator with hundreds of binding sites in the human genome. It has a main function as an insulator protein, defining together with cohesins the boundaries of areas of the genome called topologically associating domains (TADs). TADs contain regulatory elements such as enhancers which function as regulators of the transcription of genes inside the boundaries of the TAD while they are restricted from regulating genes outside these boundaries. This paper will examine the most common genetic lesions of CTCF as well as its related protein CTCFL (CTCF-like also called BORIS) in cancer using publicly available data from published genomic studies. Cancer types where abnormalities in the two genes are more common will be examined for possible associations with underlying repair defects or other prevalent genetic lesions. The putative functional effects in CTCF and CTCFL lesions will also be explored.

Disease-Associated Short Tandem Repeats Co-localize with Chromatin Domain Boundaries

More than 25 inherited human disorders are caused by the unstable expansion of repetitive DNA sequences termed short tandem repeats (STRs). A fundamental unresolved question is why some STRs are susceptible to pathologic expansion, whereas thousands of repeat tracts across the human genome are relatively stable. Here, we discover that nearly all disease-associated STRs (daSTRs) are located at boundaries demarcating 3D chromatin domains. We identify a subset of boundaries with markedly higher CpG island density compared to the rest of the genome. daSTRs specifically localize to ultra-high-density CpG island boundaries, suggesting they might be hotspots for epigenetic misregulation or topological disruption linked to STR expansion. Fragile X syndrome patients exhibit severe boundary disruption in a manner that correlates with local loss of CTCF occupancy and the degree of FMR1 silencing. Our data uncover higher-order chromatin architecture as a new dimension in understanding repeat expansion disorders.

Parvovirus minute virus of mice interacts with sites of cellular DNA damage to establish and amplify its lytic infection

We have developed a generally adaptable, novel high-throughput Viral Chromosome Conformation Capture assay (V3C-seq) for use in trans that allows genome-wide identification of the direct interactions of a lytic virus genome with distinct regions of the cellular chromosome. Upon infection, we found that the parvovirus Minute Virus of Mice (MVM) genome initially associated with sites of cellular DNA damage that in mock-infected cells also exhibited DNA damage as cells progressed through S-phase. As infection proceeded, new DNA damage sites were induced, and virus subsequently also associated with these. Sites of association identified biochemically were confirmed microscopically and MVM could be targeted specifically to artificially induced sites of DNA damage. Thus, MVM established replication at cellular DNA damage sites, which provide replication and expression machinery, and as cellular DNA damage accrued, virus spread additionally to newly damaged sites to amplify infection. MVM-associated sites overlap significantly with previously identified topologically-associated domains (TADs).

Viruses are small infectious particles that can only reproduce with the help of a host. Once they are inside their victim, they hijack the cells’ genetic material and reprogram it to become a virus factory that produces more virus particles. Parvoviruses, for example, are among the simplest of viruses and need all resources a cell has to offer to successfully replicate.

This process often takes place at so-called replication centers that contain these necessary factors. It was previously thought that parvoviruses set up such centers randomly, and gather the required molecules such as proteins to these sites. However, it was not well understood how they do this.

Now, Majumder et al. have developed a new method that enabled them to study in detail how parvoviruses gain access to the resources of the cell they need to initiate and amplify replication. The results show that parvoviruses set up their replication centers at sites on the host DNA that are already rich in proteins needed to repair and then replicate damaged DNA. Some of these sites already exist in the cell’s genetic material as a consequence of naturally occurring processes, but others are created during infection by the virus. These findings may have important implications for how other viruses may establish their replication.

Viruses, including parvoviruses, are important pathogens. Like many microbes, viruses can be beneficial for our health and environment. Others, however, can be harmful. A clearer understanding of how viruses establish and amplify an infection may provide new treatment opportunities.

SMCHD1 Merges Chromosome Compartments and Assists Formation of Super-Structures on the Inactive X

Mammalian chromosomes are partitioned into A/B compartments and topologically associated domains (TADs). The inactive X (Xi) chromosome, however, adopts a distinct conformation without evident compartments or TADs. Here, through exploration of an architectural protein, structural-maintenance-of-chromosomes hinge domain containing 1 (SMCHD1), we probe how the Xi is reconfigured during X chromosome inactivation. A/B compartments are first fused into "S1" and "S2" compartments, coinciding with Xist spreading into gene-rich domains. SMCHD1 then binds S1/S2 compartments and merges them to create a compartment-less architecture. Contrary to current views, TADs remain on the Xi but in an attenuated state. Ablating SMCHD1 results in a persistent S1/S2 organization and strengthening of TADs. Furthermore, loss of SMCHD1 causes regional defects in Xist spreading and erosion of heterochromatic silencing. We present a stepwise model for Xi folding, where SMCHD1 attenuates a hidden layer of Xi architecture to facilitate Xist spreading.

Principles of Chromosome Architecture Revealed by Hi-C

Chromosomes are folded and compacted in interphase nuclei, but the molecular basis of this folding is poorly understood. Chromosome conformation capture methods, such as Hi-C, combine chemical crosslinking of chromatin with fragmentation, DNA ligation, and high-throughput DNA sequencing to detect neighboring loci genome-wide. Hi-C has revealed the segregation of chromatin into active and inactive compartments and the folding of DNA into self-associating domains and loops. Depletion of CTCF, cohesin, or cohesin-associated proteins was recently shown to affect the majority of domains and loops in a manner that is consistent with a model of DNA folding through extrusion of chromatin loops. Compartmentation was not dependent on CTCF or cohesin. Hi-C contact maps represent the superimposition of CTCF/cohesin-dependent and -independent folding states.

Chromothripsis and ring chromosome 22: a paradigm of genomic complexity in the Phelan-McDermid syndrome (22q13 deletion syndrome)

Introduction

Phelan-McDermid syndrome (PMS) is caused by SHANK3 haploinsufficiency. Its wide phenotypic variation is attributed partly to the type and size of 22q13 genomic lesion (deletion, unbalanced translocation, ring chromosome), partly to additional undefined factors. We investigated a child with severe global neurodevelopmental delay (NDD) compatible with her distal 22q13 deletion, complicated by bilateral perisylvian polymicrogyria (BPP) and urticarial rashes, unreported in PMS.

Methods

Following the cytogenetic and array-comparative genomic hybridization (CGH) detection of a r(22) with SHANK3 deletion and two upstream duplications, whole-genome sequencing (WGS) in blood and whole-exome sequencing (WES) in blood and saliva were performed to highlight potential chromothripsis/chromoanagenesis events and any possible BPP-associated variants, even in low-level mosaicism.

Results

WGS confirmed the deletion and highlighted inversion and displaced order of eight fragments, three of them duplicated. The microhomology-mediated insertion of partial Alu-elements at one breakpoint junction disrupted the topological associating domain joining NFAM1 to the transcriptional coregulator TCF20. WES failed to detect BPP-associated variants.

Conclusions

Although we were unable to highlight the molecular basis of BPP, our data suggest that SHANK3 haploinsufficiency and TCF20 misregulation, both associated with intellectual disability, contributed to the patient’s NDD, while NFAM1 interruption likely caused her skin rashes, as previously reported. We provide the first example of chromoanasynthesis in a constitutional ring chromosome and reinforce the growing evidence that chromosomal rearrangements may be more complex than estimated by conventional diagnostic approaches and affect the phenotype by global alteration of the topological chromatin organisation rather than simply by deletion or duplication of dosage-sensitive genes.

Dynamic Rewiring of Promoter-Anchored Chromatin Loops during Adipocyte Differentiation

Interactions between transcriptional promoters and their distal regulatory elements play an important role in transcriptional regulation; however, the extent to which these interactions are subject to rapid modulations in response to signals is unknown. Here, we use promoter capture Hi-C to demonstrate a rapid reorganization of promoter-anchored chromatin loops within 4 hr after inducing differentiation of 3T3-L1 preadipocytes. The establishment of new promoter-enhancer loops is tightly coupled to activation of poised (histone H3 lysine 4 mono- and dimethylated) enhancers, as evidenced by the acquisition of histone H3 lysine 27 acetylation and the binding of MED1, SMC1, and P300 proteins to these regions, as well as to activation of target genes. Intriguingly, formation of loops connecting activated enhancers and promoters is also associated with extensive recruitment of corepressors such as NCoR and HDACs, indicating that this class of coregulators may play a previously unrecognized role during enhancer activation.

Chromatin Architecture Emerges during Zygotic Genome Activation Independent of Transcription

Chromatin architecture is fundamental in regulating gene expression. To investigate when spatial genome organization is first established during development, we examined chromatin conformation during Drosophila embryogenesis and observed the emergence of chromatin architecture within a tight time window that coincides with the onset of transcription activation in the zygote. Prior to zygotic genome activation, the genome is mostly unstructured. Early expressed genes serve as nucleation sites for topologically associating domain (TAD) boundaries. Activation of gene expression coincides with the establishment of TADs throughout the genome and co-localization of housekeeping gene clusters, which remain stable in subsequent stages of development. However, the appearance of TAD boundaries is independent of transcription and requires the transcription factor Zelda for locus-specific TAD boundary insulation. These results offer insight into when spatial organization of the genome emerges and identify a key factor that helps trigger this architecture.

Unraveling the mechanisms of chromatin fibril packaging

Recent data indicate that eukaryotic chromosomes are organized into Topologically Associating Domains (TADs); however, the mechanisms underlying TAD formation remain obscure. Based on the results of Hi-C analysis performed on 4 Drosophila melanogaster cell lines, we have proposed that specific properties of nucleosomes in active and repressed chromatin play a key role in the formation of TADs. Our computer simulations showed that the ability of "inactive" nucleosomes to stick to each other and the lack of such ability in "active" nucleosomes is sufficient for spatial segregation of these types of chromatin, which is revealed in the Hi-C analysis as TAD/inter-TAD partitioning. However, some Drosophila and mammalian TADs contain both active and inactive chromatin, a fact that does not fit this model. Herein, we present additional arguments for the model by postulating that transcriptionally active chromatin is extruded on the surface of a TAD, and discuss the possible impact of this organization on the enhancer-promoter communication and on the segregation of TADs.

Gene Regulatory Elements, Major Drivers of Human Disease

Gene expression changes, the driving forces for cellular diversity in multicellular organisms, are regulated by a diverse set of gene regulatory elements that direct transcription in specific cells. Mutations in these elements, ranging from chromosomal aberrations to single-nucleotide polymorphisms, are a major cause of human disease. However, we currently have a very limited understanding of how regulatory element genotypes lead to specific phenotypes. In this review, we discuss the various methods of regulatory element identification, the different types of mutations they harbor, and their impact on human disease. We highlight how these variations can affect transcription of multiple genes in gene regulatory networks. In addition, we describe how novel technologies, such as massively parallel reporter assays and CRISPR/Cas9 genome editing, are beginning to provide a better understanding of the functional roles that these elements have and how their alteration can lead to specific phenotypes.

Making Sense of the Tangle: Insights into Chromatin Folding and Gene Regulation

Proximity ligation assays such as circularized chromosome conformation capture and high-throughput chromosome capture assays have shed light on the structural organization of the interphase genome. Functional topologically associating domains (TADs) that constitute the building blocks of genomic organization are disrupted and reconstructed during the cell cycle. Epigenetic memory, as well as the sequence of chromosomes, regulate TAD reconstitution. Sub-TAD domains that are invariant across cell types have been identified, and contacts between these domains, rather than looping, are speculated to drive chromatin folding. Replication domains are established simultaneously with TADs during the cell cycle and the two correlate well in terms of characteristic features, such as lamin association and histone modifications. CCCTC-binding factor (CTCF) and cohesin cooperate across different cell types to regulate genes and genome organization. CTCF elements that demarcate TAD boundaries are commonly disrupted in cancer and promote oncogene activation. Chromatin looping facilitates interactions between distant promoters and enhancers, and the resulting enhanceosome complex promotes gene expression. Deciphering the chromatin tangle requires comprehensive integrative analyses of DNA- and protein-dependent factors that regulate genomic organization.

Nanoscale spatial organization of the HoxD gene cluster in distinct transcriptional states

Chromatin condensation plays an important role in the regulation of gene expression. Recently, it was shown that the transcriptional activation of Hoxd genes during vertebrate digit development involves modifications in 3D interactions within and around the HoxD gene cluster. This reorganization follows a global transition from one set of regulatory contacts to another, between two topologically associating domains (TADs) located on either side of the HoxD locus. Here, we use 3D DNA FISH to assess the spatial organization of chromatin at and around the HoxD gene cluster and report that although the two TADs are tightly associated, they appear as spatially distinct units. We measured the relative position of genes within the cluster and found that they segregate over long distances, suggesting that a physical elongation of the HoxD cluster can occur. We analyzed this possibility by super-resolution imaging (STORM) and found that tissues with distinct transcriptional activity exhibit differing degrees of elongation. We also observed that the morphological change of the HoxD cluster in developing digits is associated with its position at the boundary between the two TADs. Such variations in the fine-scale architecture of the gene cluster suggest causal links among its spatial configuration, transcriptional activation, and the flanking chromatin context.

Topology, structures, and energy landscapes of human chromosomes

Chromosome conformation capture experiments provide a rich set of data concerning the spatial organization of the genome. We use these data along with a maximum entropy approach to derive a least-biased effective energy landscape for the chromosome. Simulations of the ensemble of chromosome conformations based on the resulting information theoretic landscape not only accurately reproduce experimental contact probabilities, but also provide a picture of chromosome dynamics and topology. The topology of the simulated chromosomes is probed by computing the distribution of their knot invariants. The simulated chromosome structures are largely free of knots. Topologically associating domains are shown to be crucial for establishing these knotless structures. The simulated chromosome conformations exhibit a tendency to form fibril-like structures like those observed via light microscopy. The topologically associating domains of the interphase chromosome exhibit multistability with varying liquid crystalline ordering that may allow discrete unfolding events and the landscape is locally funneled toward "ideal" chromosome structures that represent hierarchical fibrils of fibrils.