3D Genome Structure Research Articles

Cohesin and CTCF emerge as building blocks of 3D genome structure.

Cohesin and CTCF emerge as building blocks of 3D genome structure.

A guide to studying 3D genome structure and dynamics in the kidney.

The human genome is tightly packed into the 3D environment of the cell nucleus. Rapidly evolving and sophisticated methods of mapping 3D genome architecture have shed light on fundamental principles of genome organization and gene regulation. The genome is physically organized on different scales, from individual genes to entire chromosomes. Nuclear landmarks such as the nuclear envelope and nucleoli have important roles in compartmentalizing the genome within the nucleus. Genome activity (for example, gene transcription) is also functionally partitioned within this 3D organization. Rather than being static, the 3D organization of the genome is tightly regulated over various time scales. These dynamic changes in genome structure over time represent the fourth dimension of the genome. Innovative methods have been used to map the dynamic regulation of genome structure during important cellular processes including organism development, responses to stimuli, cell division and senescence. Furthermore, disruptions to the 4D genome have been linked to various diseases, including of the kidney. As tools and approaches to studying the 4D genome become more readily available, future studies that apply these methods to study kidney biology will provide insights into kidney function in health and disease.

Analysis and Visualization of Multiple Hi-C and Micro-C Data with CustardPy.

Three-dimensional (3D) genome structure plays crucial roles in biological processes and disease pathogenesis. Hi-C and Micro-C, well-established methods for 3D genome analysis, can identify a variety of 3D genome structures. However, selecting appropriate pipelines and tools for the analysis and setting up the required computing environment can sometimes pose challenges. To address this, we have introduced CustardPy, a Docker-based pipeline specifically designed for 3D genome analysis. CustardPy is designed to compare and evaluate multiple samples and wraps several existing tools to cover the entire workflow from FASTQ mapping to visualization. In this chapter, we demonstrate how to analyze and visualize Hi-C data using CustardPy and introduce several 3D genome features observed in Hi-C data.

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

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

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

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

ARGV: 3D genome structure exploration using augmented reality

Over the past two decades, scientists have increasingly realized the importance of the three-dimensional (3D) genome organization in regulating cellular activity. Hi-C and related experiments yield 2D contact matrices that can be used to infer 3D models of chromosome structure. Visualizing and analyzing genomes in 3D space remains challenging. Here, we present ARGV, an augmented reality 3D Genome Viewer. ARGV contains more than 350 pre-computed and annotated genome structures inferred from Hi-C and imaging data. It offers interactive and collaborative visualization of genomes in 3D space, using standard mobile phones or tablets. A user study comparing ARGV to existing tools demonstrates its benefits.

Deciphering Pre-existing and Induced 3D Genome Architecture Changes involved in Constricted Melanoma Migration.

Metastatic cancer cells traverse constricted spaces that exert forces on their nucleus and the genomic contents within. Cancerous tumors are highly heterogeneous and not all cells within them can achieve such a feat. Here, we investigated what initial genome architecture characteristics favor the constricted migratory ability of cancer cells and which arise only after passage through multiple constrictions. We identified a cell surface protein (ITGB4) whose expression correlates with increased initial constricted migration ability in human melanoma A375 cells. Sorting out this subpopulation allowed us to identify cellular and nuclear features that pre-exist and favor migration, as well as alterations that only appear after cells have passed through constrictions. We identified specific genomic regions that experienced altered genome spatial compartment profiles only after constricted migration. Our study reveals 3D genome structure contributions to both selection and induction mechanisms of cell fate change during cancer metastasis.

Loss of multi-level 3D genome organization during breast cancer progression.

Breast cancer entails intricate alterations in genome organization and expression. However, how three-dimensional (3D) chromatin structure changes in the progression from a normal to a breast cancer malignant state remains unknown. To address this, we conducted an analysis combining Hi-C data with lamina-associated domains (LADs), epigenomic marks, and gene expression in an in vitro model of breast cancer progression. Our results reveal that while the fundamental properties of topologically associating domains (TADs) are overall maintained, significant changes occur in the organization of compartments and subcompartments. These changes are closely correlated with alterations in the expression of oncogenic genes. We also observe a restructuring of TAD-TAD interactions, coinciding with a loss of spatial compartmentalization and radial positioning of the 3D genome. Notably, we identify a previously unrecognized interchromosomal insertion event, wherein a locus on chromosome 8 housing the MYC oncogene is inserted into a highly active subcompartment on chromosome 10. This insertion is accompanied by the formation of de novo enhancer contacts and activation of MYC , illustrating how structural genomic variants can alter the 3D genome to drive oncogenic states. In summary, our findings provide evidence for the loss of genome organization at multiple scales during breast cancer progression revealing novel relationships between genome 3D structure and oncogenic processes.

Reversible acetylation of HDAC8 regulates cell cycle

HDAC8, a member of class I HDACs, plays a pivotal role in cell cycle regulation by deacetylating the cohesin subunit SMC3. While cyclins and CDKs are well-established cell cycle regulators, our knowledge of other regulators remains limited. Here we reveal the acetylation of K202 in HDAC8 as a key cell cycle regulator responsive to stress. K202 acetylation in HDAC8, primarily catalyzed by Tip60, restricts HDAC8 activity, leading to increased SMC3 acetylation and cell cycle arrest. Furthermore, cells expressing the mutant form of HDAC8 mimicking K202 acetylation display significant alterations in gene expression, potentially linked to changes in 3D genome structure, including enhanced chromatid loop interactions. K202 acetylation impairs cell cycle progression by disrupting the expression of cell cycle-related genes and sister chromatid cohesion, resulting in G2/M phase arrest. These findings indicate the reversible acetylation of HDAC8 as a cell cycle regulator, expanding our understanding of stress-responsive cell cycle dynamics.

Hi-C-guided many-polymer model to decipher 3D genome organization

We propose a high-throughput chromosome conformation capture data-based many-polymer model that allows us to generate an ensemble of multi-scale genome structures. We demonstrate the efficacy of our model by validating the generated structures against experimental measurements and employ them to address key questions regarding genome organization. Our model first confirms a significant correlation between chromosome size and nuclear positioning. Specifically, smaller chromosomes are distributed at the core region, whereas larger chromosomes are at the periphery, interacting with the nuclear envelope. The spatial distribution of A- and B-type compartments, which is nontrivial to infer from the corresponding high-throughput chromosome conformation capture maps alone, can also be elucidated using our model, accounting for an issue such as the effect of chromatin-lamina interaction on the compartmentalization of conventional and inverted nuclei. In accordance with imaging data, the overall shape of the 3D genome structures generated from our model displays significant variation. As a case study, we apply our method to the yellow fever mosquito genome, finding that the predicted morphology displays, on average, a more globular shape than the previously suggested spindle-like organization and that our prediction better aligns with the fluorescence in situ hybridization data. Our model has great potential to be extended to investigate many outstanding issues concerning 3D genome organization.

Reprogramming of 3D genome structure underlying HSPC development in zebrafish

BackgroundDevelopment of hematopoietic stem and progenitor cells (HSPC) is a multi-staged complex process that conserved between zebrafish and mammals. Understanding the mechanism underlying HSPC development is a holy grail of hematopoietic biology, which is helpful for HSPC clinical application. Chromatin conformation plays important roles in transcriptional regulation and cell fate decision; however, its dynamic and role in HSPC development is poorly investigated.MethodsWe performed chromatin structure and multi-omics dissection across different stages of HSPC developmental trajectory in zebrafish for the first time, including Hi-C, RNA-seq, ATAC-seq, H3K4me3 and H3K27ac ChIP-seq.ResultsThe chromatin organization of zebrafish HSPC resemble mammalian cells with similar hierarchical structure. We revealed the multi-scale reorganization of chromatin structure and its influence on transcriptional regulation and transition of cell fate during HSPC development. Nascent HSPC is featured by loose conformation with obscure structure at all layers. Notably, PU.1 was identified as a potential factor mediating formation of promoter-involved loops and regulating gene expression of HSPC.ConclusionsOur results provided a global view of chromatin structure dynamics associated with development of zebrafish HSPC and discovered key transcription factors involved in HSPC chromatin interactions, which will provide new insights into the epigenetic regulatory mechanisms underlying vertebrate HSPC fate decision.

A three-dimensional genomics view for speciation research

Genomic information is folded in a three-dimensional (3D) structure, a rarely explored evolutionary driver of speciation. Technological advances now enable the study of 3D genome structures (3DGSs) across the Tree of Life. At the onset of 3D speciation genomics, we discuss the putative roles of 3DGSs in speciation.

Estrogen-induced chromatin looping changes identify a subset of functional regulatory elements.

Transcriptional enhancers can regulate individual or multiple genes through long-range three-dimensional (3D) genome interactions, and these interactions are commonly altered in cancer. Yet, the functional relationship between changes in 3D interactions associated with regulatory regions and differential gene expression appears context-dependent. In this study, we used HiChiP to capture changes in 3D genome interactions between active regulatory regions of endometrial cancer cells in response to estrogen treatment and uncovered significant differential long-range interactions that are strongly enriched for estrogen receptor α (ER) bound sites (ERBS). The ERBS anchoring differential loops with either a gene's promoter or distal regions were correlated with larger transcriptional responses to estrogen compared to ERBS not involved in differential interactions. To functionally test this observation, CRISPR-based Enhancer-i was used to deactivate specific ERBS, which revealed a wide range of effects on the transcriptional response to estrogen. However, these effects are only subtly and not significantly stronger for ERBS in differential loops. In addition, we observed an enrichment of 3D interactions between the promoters of estrogen up-regulated genes and found that looped promoters can work together cooperatively. Overall, our work suggests that changes in 3D genome structure upon estrogen treatment identify some functionally important regulatory regions; however, these changes aren't required for a transcriptional response to E2 in endometrial cancer cells.

Altered chromatin topologies caused by balanced chromosomal translocation lead to central iris hypoplasia

Despite the advent of genomic sequencing, molecular diagnosis remains unsolved in approximately half of patients with Mendelian disorders, largely due to unclarified functions of noncoding regions and the difficulty in identifying complex structural variations. In this study, we map a unique form of central iris hypoplasia in a large family to 6q15-q23.3 and 18p11.31-q12.1 using a genome-wide linkage scan. Long-read sequencing reveals a balanced translocation t(6;18)(q22.31;p11.22) with intergenic breakpoints. By performing Hi-C on induced pluripotent stem cells from a patient, we identify two chromatin topologically associating domains spanning across the breakpoints. These alterations lead the ectopic chromatin interactions between APCDD1 on chromosome 18 and enhancers on chromosome 6, resulting in upregulation of APCDD1. Notably, APCDD1 is specifically localized in the iris of human eyes. Our findings demonstrate that noncoding structural variations can lead to Mendelian diseases by disrupting the 3D genome structure and resulting in altered gene expression.

A Constitutive Heterochromatic Region Shapes Genome Organization and Impacts Gene Expression in Neurospora crassa.

Genome organization is essential for proper function, including gene expression. In metazoan genome organization, chromatin loops and Topologically Associated Domains (TADs) facilitate local gene clustering, while chromosomes form distinct nuclear territories characterized by compartmentalization of silent heterochromatin at the nuclear periphery and active euchromatin in the nucleus center. A similar hierarchical organization occurs in the fungus Neurospora crassa where its seven chromosomes form a Rabl conformation, where heterochromatic centromeres and telomeres independently cluster at the nuclear membrane, while interspersed heterochromatic loci in Neurospora aggregate across megabases of linear genomic distance for forming TAD-like structures. However, the role of individual heterochromatic loci in normal genome organization and function is unknown. Here, we examined the genome organization of a Neurospora strain harboring a ~47.4 kilobase facultative (temporarily silent) heterochromatic region deletion, as well as the genome organization of a strain deleted of a 110.6 kilobase permanently silent constitutive heterochromatic region. While the facultative heterochromatin deletion had little effect on local chromatin structure, the constitutive heterochromatin deletion altered local TAD-like structures, gene expression, and the predicted 3D genome structure by qualitatively repositioning genes into the nucleus center. Our work elucidates the role of individual heterochromatic regions for genome organization and function.

Stratifying TAD boundaries pinpoints focal genomic regions of regulation, damage, and repair.

Advances in chromatin mapping have exposed the complex chromatin hierarchical organization in mammals, including topologically associating domains (TADs) and their substructures, yet the functional implications of this hierarchy in gene regulation and disease progression are not fully elucidated. Our study delves into the phenomenon of shared TAD boundaries, which are pivotal in maintaining the hierarchical chromatin structure and regulating gene activity. By integrating high-resolution Hi-C data, chromatin accessibility, and DNA double-strand breaks (DSBs) data from various cell lines, we systematically explore the complex regulatory landscape at high-level TAD boundaries. Our findings indicate that these boundaries are not only key architectural elements but also vibrant hubs, enriched with functionally crucial genes and complex transcription factor binding site-clustered regions. Moreover, they exhibit a pronounced enrichment of DSBs, suggesting a nuanced interplay between transcriptional regulation and genomic stability. Our research provides novel insights into the intricate relationship between the 3D genome structure, gene regulation, and DNA repair mechanisms, highlighting the role of shared TAD boundaries in maintaining genomic integrity and resilience against perturbations. The implications of our findings extend to understanding the complexities of genomic diseases and open new avenues for therapeutic interventions targeting the structural and functional integrity of TAD boundaries.

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

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

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.

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.

Simultaneous single-cell three-dimensional genome and gene expression profiling uncovers dynamic enhancer connectivity underlying olfactory receptor choice

The simultaneous measurement of three-dimensional (3D) genome structure and gene expression of individual cells is critical for understanding a genome’s structure–function relationship, yet this is challenging for existing methods. Here we present ‘Linking mRNA to Chromatin Architecture (LiMCA)’, which jointly profiles the 3D genome and transcriptome with exceptional sensitivity and from low-input materials. Combining LiMCA and our high-resolution scATAC-seq assay, METATAC, we successfully characterized chromatin accessibility, as well as paired 3D genome structures and gene expression information, of individual developing olfactory sensory neurons. We expanded the repertoire of known olfactory receptor (OR) enhancers and discovered unexpected rules of their dynamics: OR genes and their enhancers are most accessible during early differentiation. Furthermore, we revealed the dynamic spatial relationship between ORs and enhancers behind stepwise OR expression. These findings offer valuable insights into how 3D connectivity of ORs and enhancers dynamically orchestrate the ‘one neuron–one receptor’ selection process.