3D Chromatin Research Articles

3D chromatin maps of a brown alga reveal U/V sex chromosome spatial organization.

Nuclear three dimensional (3D) folding of chromatin structure has been linked to gene expression regulation and correct developmental programs, but little is known about the 3D architecture of sex chromosomes within the nucleus, and how that impacts their role in sex determination. Here, we determine the sex-specific 3D organization of the model brown alga Ectocarpus chromosomes at 2 kb resolution, by mapping long-range chromosomal interactions using Hi-C coupled with Oxford Nanopore long reads. We report that Ectocarpus interphase chromatin exhibits a non-Rabl conformation, with strong contacts among telomeres and among centromeres, which feature centromere-specific LTR retrotransposons. The Ectocarpus chromosomes do not contain large local interactive domains that resemble TADs described in animals, but their 3D genome organization is largely shaped by post-translational modifications of histone proteins. We show that the sex determining region (SDR)within the U and V chromosomes are insulated and span the centromeres andwe link sex-specific chromatin dynamics and gene expression levels to the 3D chromatin structure of the U and V chromosomes. Finally, we uncover the unique conformation of a large genomic region on chromosome 6 harboring an endogenous viral element, providing insights regarding the impact of a latent giant dsDNA virus on the host genome's 3D chromosomal folding.

Chromatin structure and 3D architecture define the differential functions of PU.1 regulatory elements in blood cell lineages

The precise spatiotemporal expression of the hematopoietic ETS transcription factor PU.1, a key determinant of hematopoietic cell fates, is tightly regulated at the chromatin level. However, how chromatin signatures are linked to this dynamic expression pattern across different blood cell lineages remains uncharacterized. Here, we performed an in-depth analysis of the relationships between gene expression, chromatin structure, 3D architecture, and trans-acting factors at PU.1 cis-regulatory elements (PCREs). By identifying phylogenetically conserved DNA elements within chromatin-accessible regions in primary human blood lineages, we discovered multiple novel candidate PCREs within the upstream region of the human PU.1 locus. A subset of these elements localizes within an 8-kb-wide cluster exhibiting enhancer features, including open chromatin, demethylated DNA, enriched enhancer histone marks, present enhancer RNAs, and PU.1 occupation, presumably mediating PU.1 autoregulation. Importantly, we revealed the presence of a common 35-kb-wide CTCF-flanked insulated neighborhood that contains the PCRE cluster (PCREC), forming a chromatin territory for lineage-specific and PCRE-mediated chromatin interactions. These include functional PCRE-promoter interactions in myeloid and B cells that are absent in erythroid and T cells. By correlating chromatin structure and 3D architecture with PU.1 expression in various lineages, we were able to attribute enhancer versus silencer functions to individual elements. Our findings provide mechanistic insights into the interplay between dynamic chromatin structure and 3D architecture in the chromatin regulation of PU.1 expression. This study lays crucial groundwork for additional experimental studies that validate and dissect the role of PCREs in epigenetic regulation of normal and malignant hematopoiesis.

EMF1 functions as a 3D chromatin modulator in Arabidopsis.

EMF1 functions as a 3D chromatin modulator in Arabidopsis.

ChromaFold predicts the 3D contact map from single-cell chromatin accessibility

Identifying cell-type-specific 3D chromatin interactions between regulatory elements can help decipher gene regulation and interpret disease-associated non-coding variants. However, achieving this resolution with current 3D genomics technologies is often infeasible given limited input cell numbers. We therefore present ChromaFold, a deep learning model that predicts 3D contact maps, including regulatory interactions, from single-cell ATAC sequencing (scATAC-seq) data alone. ChromaFold uses pseudobulk chromatin accessibility, co-accessibility across metacells, and a CTCF motif track as inputs and employs a lightweight architecture to train on standard GPUs. Trained on paired scATAC-seq and Hi-C data in human samples, ChromaFold accurately predicts the 3D contact map and peak-level interactions across diverse human and mouse test cell types. Compared to leading contact map prediction models that use ATAC-seq and CTCF ChIP-seq, ChromaFold achieves state-of-the-art performance using only scATAC-seq. Finally, fine-tuning ChromaFold on paired scATAC-seq and Hi-C in a complex tissue enables deconvolution of chromatin interactions across cell subpopulations.

Mediator Subunit Med4 Enforces Metastatic Dormancy in Breast Cancer.

Long term survival of breast cancer patients is limited due to recurrence from metastatic dormant cancer cells. However, the mechanisms by which these dormant breast cancer cells survive and awaken remain poorly understood. Our unbiased genome-scale genetic screen in mice identified Med4 as a novel cancer-cell intrinsic gatekeeper in metastatic reactivation. MED4 haploinsufficiency is prevalent in metastatic breast cancer patients and correlates with poorer prognosis. Syngeneic xenograft models revealed that Med4 enforces breast cancer dormancy. Contrary to the canonical function of the Mediator complex in activating gene expression, Med4 maintains 3D chromatin compaction and enhancer landscape, by preventing enhancer priming or activation through the suppression of H3K4me1 deposition. Med4 haploinsufficiency disrupts enhancer poise and reprograms the enhancer dynamics to facilitate extracellular matrix (ECM) gene expression and integrin-mediated mechano-transduction, driving metastatic growth. Our findings establish Med4 as a key regulator of cellular dormancy and a potential biomarker for high-risk metastatic relapse.

Type II topoisomerases shape multi-scale 3D chromatin folding in regions of positive supercoils

Type II topoisomerases shape multi-scale 3D chromatin folding in regions of positive supercoils

Response eQTLs, chromatin accessibility, and 3D chromatin structure in chondrocytes provide mechanistic insight into osteoarthritis risk.

Osteoarthritis (OA) poses a significant healthcare burden with limited treatment options. While genome-wide association studies (GWAS) have identified over 100 OA-associated loci, translating these findings into therapeutic targets remains challenging. Integrating expression quantitative trait loci (eQTL), 3D chromatin structure, and other genomic approaches with OA GWAS data offers a promising approach to elucidate disease mechanisms; however, comprehensive eQTL maps in OA-relevant tissues and conditions remain scarce. We mapped gene expression, chromatin accessibility, and 3D chromatin structure in primary human articular chondrocytes in both resting and OA-mimicking conditions. We identified thousands of differentially expressed genes, including those associated with differences in sex and age. RNA-seq in chondrocytes from 101 donors across two conditions uncovered 3782 unique eGenes, including 420 that exhibited strong and significant condition-specific effects. Colocalization with OA GWAS signals revealed 13 putative OA risk genes, 10 of which have not been previously identified. Chromatin accessibility and 3D chromatin structure provided insights into the mechanisms and conditional specificity of these variants. Our findings shed light on OA pathogenesis and highlight potential targets for therapeutic development.

Mutations of PDS5 genes enhance TAD-like domain formation Arabidopsis thaliana

In eukaryotes, topologically associating domains (TADs) organize the genome into functional compartments. While TAD-like structures are common in mammals and many plants, they are challenging to detect in Arabidopsis thaliana. Here, we demonstrate that Arabidopsis PDS5 proteins play a negative role in TAD-like domain formation. Through Hi-C analysis, we show that mutations in PDS5 genes lead to the widespread emergence of enhanced TAD-like domains throughout the Arabidopsis genome, excluding pericentromeric regions. These domains exhibit increased chromatin insulation and enhanced chromatin interactions, without significant changes in gene expression or histone modifications. Our results suggest that PDS5 proteins are key regulators of genome architecture, influencing 3D chromatin organization independently of transcriptional activity. This study provides insights into the unique chromatin structure of Arabidopsis and the broader mechanisms governing plant genome folding.

Integrative 3D genomics with multi-omics analysis and functional validation of genetic regulatory mechanisms of abdominal fat deposition in chickens

Chickens are the most abundant agricultural animals globally, with controlling abdominal fat deposition being a key objective in poultry breeding. While GWAS can identify genetic variants associated with abdominal fat deposition, the precise roles and mechanisms of these variants remain largely unclear. Here, we use male chickens from two lines divergently selected for abdominal fat deposition as experimental models. Through the integration of genomic, epigenomic, 3D genomic, and transcriptomic data, we build a comprehensive chromatin 3D regulatory network map to identify the genetic regulatory mechanisms that influence abdominal fat deposition in chickens. Notably, we find that the rs734209466 variant functions as an allele-specific enhancer, remotely enhancing the transcription of IGFBP2 and IGFBP5 by the binding transcription factor IRF4. This interaction influences the differentiation and proliferation of preadipocytes, which ultimately affects phenotype. This work presents a detailed genetic regulatory map for chicken abdominal fat deposition, offering molecular targets for selective breeding.

Targeting the 3D genome by anthracyclines for chemotherapeutic effects.

The chromatins are folded into three-dimensional (3D) structures inside cells, which coordinates the regulation of gene transcription by the non-coding regulatory elements. Aberrant chromatin 3D folding has been shown in many diseases, such as acute myeloid leukemia (AML), and may contribute to tumorigenesis. The anthracycline topoisomerase II inhibitors can induce histone eviction and DNA damage. We performed genome-wide high-resolution mapping of the chemotherapeutic effects of various clinically used anthracycline drugs. ATAC-seq was used to profile the histone eviction effects of different anthracyclines. TOP2A ChIP-seq was used to profile the potential DNA damage regions. Integrated analyses show that different anthracyclines have distinct target selectivity on epigenomic regions, based on their respective ATAC-seq and ChIP-seq profiles. We identified the underlying molecular mechanism that unique anthracycline variants selectively target chromatin looping anchors via disrupting CTCF binding, suggesting an additional potential therapeutic effect on the 3D genome. We further performed Hi-C experiments, and data from K562 cells treated with the selective anthracycline drugs indicate that the 3D chromatin organization is disrupted. Furthermore, AML patients receiving anthracycline drugs showed altered chromatin structures around potential looping anchors, which linked to distinct clinical outcomes. Our data indicate that anthracyclines are potent and selective epigenomic targeting drugs and can target the 3D genome for anticancer therapy, which could be used for personalized medicine to treat tumors with aberrant 3D chromatin structures.

Drought-responsive dynamics of H3K9ac-marked 3D chromatin interactions are integrated by OsbZIP23-associated super-enhancer-like promoter regions in rice

BackgroundIn response to drought stress (DS), plants undergo complex processes that entail significant transcriptome reprogramming. However, the intricate relationship between the dynamic alterations in the three-dimensional (3D) genome and the modulation of gene co-expression in drought responses remains a relatively unexplored area.ResultsIn this study, we reconstruct high-resolution 3D genome maps based on genomic regions marked by H3K9ac, an active histone modification that dynamically responds to soil water variations in rice. We discover a genome-wide disconnection of 3D genome contact upon DS with over 10,000 chromatin loops lost, which are partially recovered in the subsequent re-watering. Loops integrating promoter–promoter interactions (PPI) contribute to gene expression in addition to basal H3K9ac modifications. Moreover, H3K9ac-marked promoter regions with high affinities in mediating PPIs, termed as super-promoter regions (SPRs), integrate spatially clustered PPIs in a super-enhancer-like manner. Interestingly, the knockout mutation of OsbZIP23, a well-defined DS-responsive transcription factor, leads to the disassociation of over 80% DS-specific PPIs and decreased expression of the corresponding genes under DS. As a case study, we show how OsbZIP23 integrates the PPI cluster formation and the co-expression of four dehydrin genes, RAB16A–D, through targeting the RAB16C SPR in a stress signaling-dependent manner.ConclusionsOur high-resolution 3D genome maps unveil the principles and details of dynamic genome folding in response to water supply variations and illustrate OsbZIP23 as an indispensable integrator of the yet unique 3D genome organization that is essential for gene co-expression under DS in rice.

7690 Acute hormonal signaling-induced 3D genome organization controls transcriptional dynamics during adipocyte thermogenesis

Abstract Disclosure: Y. Zhang: None. R. Zheng: None. T. Tsuji: None. Y. Xing: None. C. Wang: None. J. Darcy: None. K. Chen: None. Y. Tseng: None. Brown adipose tissue (BAT) maintains body temperature by dissipating energy as heat, presenting a promising target against obesity and its associated metabolic disorders. When exposed to cold, the β3-adrenergic receptor (β3-AR) in brown adipocytes is stimulated, initiating a sequential signaling cascade and activating a thermogenic gene program essential for heat generation. Despite remarkable progress has been made towards understanding the engagement of transcriptional and epigenetic regulators during thermogenesis in response to β3-AR stimulation, the precise impact of three-dimensional (3D) genome organization, an increasingly important modulator of gene expression, on the activation of thermogenic gene program remains elusive. In this study, we employed high-resolution Micro-C in murine brown adipocytes to interrogate the rapidity of 3D genome reorganization in response to β3-AR hormonal signaling and its impact on transcriptional dynamics during thermogenesis. Our findings unveiled substantial rewiring of chromatin loops within 4 hours of β3-AR stimulation. These dynamic changes in chromatin loops, particularly enhancer-promoter (E-P) interactions, were coupled with gene activation during thermogenesis. Furthermore, we uncovered the mechanisms by which β3-AR hormonal signaling mediates 3D chromatin reorganization to regulate the thermogenic gene program. We found that β3-AR hormonal signaling induced the phosphorylation of p18Hamlet, a subunit of the SRCAP complex, facilitating the deposition of the histone variant H2A.Z into regulatory DNA regions. This nucleosome incorporation of H2A.Z further mediated nucleosome rearrangement, favoring the formation of nucleosome-depleted regions (NDRs) and enhancing the accessibility of regulatory DNA regions. As a result, H2A.Z chromatin incorporation facilitated interactions between cis-regulatory regions and gene promoters, such as E-P interactions, thereby activating thermogenic gene expression in response to β3-AR stimulation. Loss of H2A.Z disrupted loop dynamics, impeding gene activation and BAT thermogenic capacity. The conserved role of H2A.Z in human thermogenesis was validated, and we identified adiposity-related SNPs in human homologous loop anchors that could potentially be targeted for treating obesity. Taken together, our results unveil a previously underappreciated aspect of acute hormonal signaling-induced chromatin dynamics in orchestrating the thermogenic program in brown adipocytes. These findings not only enhance our understanding of molecular regulation of thermogenesis but also shed light on the role of rapid and dynamic chromatin reconfiguration in regulating the cellular response to external stimuli. Our study bridges a critical gap between 3D genome organization and gene activation in the thermogenic program mediated by acute β3-AR hormonal signaling. Presentation: 6/2/2024

Machine Learning Reveals the Diversity of Human 3D Chromatin Contact Patterns.

Understanding variation in chromatin contact patterns across diverse humans is critical for interpreting noncoding variants and their effects on gene expression and phenotypes. However, experimental determination of chromatin contact patterns across large samples is prohibitively expensive. To overcome this challenge, we develop and validate a machine learning method to quantify the variation in 3D chromatin contacts at 2 kilobase resolution from genome sequence alone. We apply this approach to thousands of human genomes from the 1000 Genomes Project and the inferred hominin ancestral genome. While patterns of 3D contact divergence genome wide are qualitatively similar to patterns of sequence divergence, we find substantial differences in 3D divergence and sequence divergence in local 1 megabase genomic windows. In particular, we identify 392 windows with significantly greater 3D divergence than expected from sequence. Moreover, for 31% of genomic windows, a single individual has a rare divergent 3D contact map pattern. Using in silico mutagenesis, we find that most single nucleotide sequence changes do not result in changes to 3D chromatin contacts. However, in windows with substantial 3D divergence just one or a few variants can lead to divergent 3D chromatin contacts without the individuals carrying those variants having high sequence divergence. In summary, inferring 3D chromatin contact maps across human populations reveals variable contact patterns. We anticipate that these genetically diverse maps of 3D chromatin contact will provide a reference for future work on the function and evolution of 3D chromatin contact variation across human populations.

Advancements and future directions in single-cell Hi-C based 3D chromatin modeling

Advancements and future directions in single-cell Hi-C based 3D chromatin modeling

Multiscale molecular modeling of chromatin with MultiMM: From nucleosomes to the whole genome

Multiscale molecular modeling of chromatin with MultiMM: From nucleosomes to the whole genome

De novo genome assemblies of two cryptodiran turtles with ZZ/ZW and XX/XY sex chromosomes provide insights into patterns of genome reshuffling and uncover novel 3D genome folding in amniotes.

Understanding the evolution of chromatin conformation among species is fundamental to elucidate the architecture and plasticity of genomes. Nonrandom interactions of linearly distant loci regulate gene function in species-specific patterns, affecting genome function, evolution, and, ultimately, speciation. Yet, data from nonmodel organisms are scarce. To capture the macroevolutionary diversity of vertebrate chromatin conformation, here we generate de novo genome assemblies for two cryptodiran (hidden-neck) turtles via Illumina sequencing, chromosome conformation capture, and RNA-seq: Apalone spinifera (ZZ/ZW, 2n = 66) and Staurotypus triporcatus (XX/XY, 2n = 54). We detected differences in the three-dimensional (3D) chromatin structure in turtles compared to other amniotes beyond the fusion/fission events detected in the linear genomes. Namely, whole-genome comparisons revealed distinct trends of chromosome rearrangements in turtles: (1) a low rate of genome reshuffling in Apalone (Trionychidae) whose karyotype is highly conserved when compared to chicken (likely ancestral for turtles), and (2) a moderate rate of fusions/fissions in Staurotypus (Kinosternidae) and Trachemys scripta (Emydidae). Furthermore, we identified a chromosome folding pattern that enables "centromere-telomere interactions" previously undetected in turtles. The combined turtle pattern of "centromere-telomere interactions" (discovered here) plus "centromere clustering" (previously reported in sauropsids) is novel for amniotes and it counters previous hypotheses about amniote 3D chromatin structure. We hypothesize that the divergent pattern found in turtles originated from an amniote ancestral state defined by a nuclear configuration with extensive associations among microchromosomes that were preserved upon the reshuffling of the linear genome.

The chromatin tapestry as a framework for neurodevelopment.

The neuronal nucleus houses a meticulously organized genome. Within this structure, genetic material is not simply compacted but arranged into a precise and functional 3D chromatin landscape essential for cellular regulation. This mini-review highlights the importance of this chromatin landscape in healthy neurodevelopment, as well as the diseases that occur with aberrant chromatin architecture. We discuss insights into the fundamental mechanistic relationship between histone modifications, DNA methylation, and genome organization. We then discuss findings that reveal how these epigenetic features change throughout normal neurodevelopment. Finally, we highlight single-gene neurodevelopmental disorders that illustrate the interdependence of epigenetic features, showing how disruptions in DNA methylation or genome architecture can ripple across the entire epigenome. As such, we emphasize the importance of measuring multiple chromatin architectural aspects, as the disruption of one mechanism can likely impact others in the intricate epigenetic network. This mini-review underscores the vast gaps in our understanding of chromatin structure in neurodevelopmental diseases and the substantial research needed to understand the interplay between chromatin features and neurodevelopment.

Integration of scHi-C and scRNA-seq data defines distinct 3D-regulated and biological-context dependent cell subpopulations

An integration of 3D chromatin structure and gene expression at single-cell resolution has yet been demonstrated. Here, we develop a computational method, a multiomic data integration (MUDI) algorithm, which integrates scHi-C and scRNA-seq data to precisely define the 3D-regulated and biological-context dependent cell subpopulations or topologically integrated subpopulations (TISPs). We demonstrate its algorithmic utility on the publicly available and newly generated scHi-C and scRNA-seq data. We then test and apply MUDI in a breast cancer cell model system to demonstrate its biological-context dependent utility. We find the newly defined topologically conserved associating domain (CAD) is the characteristic single-cell 3D chromatin structure and better characterizes chromatin domains in single-cell resolution. We further identify 20 TISPs uniquely characterizing 3D-regulated breast cancer cellular states. We reveal two of TISPs are remarkably resemble to high cycling breast cancer persister cells and chromatin modifying enzymes might be functional regulators to drive the alteration of the 3D chromatin structures. Our comprehensive integration of scHi-C and scRNA-seq data in cancer cells at single-cell resolution provides mechanistic insights into 3D-regulated heterogeneity of developing drug-tolerant cancer cells.

3D genome contributes to MHC-II neoantigen prediction.

Reliable and ultra-fast DNA and RNA sequencing have been achieved with the emergence of high-throughput sequencing technology. When combining the results of DNA and RNA sequencing for tumor cells of cancer patients, neoantigens that potentially stimulate the immune response of either CD4+ or CD8+ T cells can be identified. However, due to the abundance of somatic mutations and the high polymorphic nature of human leukocyte antigen (HLA) it is challenging to accurately predict the neoantigens. Moreover, comparing to HLA-I presented peptides, the HLA-II presented peptides are more variable in length, making the prediction of HLA-II loaded neoantigens even harder. A number of computational approaches have been proposed to address this issue but none of them considers the DNA origin of the neoantigens from the perspective of 3D genome. Here we investigate the DNA origins of the immune-positive and non-negative HLA-II neoantigens in the context of 3D genome and discovered that the chromatin 3D architecture plays an important role in more effective HLA-II neoantigen prediction. We believe that the 3D genome information will help to increase the precision of HLA-II neoantigen discovery and eventually benefit precision and personalized medicine in cancer immunotherapy.

Peculiar k-mer Spectra Are Correlated with 3D Contact Frequencies and Breakpoint Regions in the Human Genome.

It is widely accepted that the 3D chromatin organization in human cell nuclei is not random and recent investigations point towards an interactive relation of epigenetic functioning and chromatin (re-)organization. Although chromatin organization seems to be the result of self-organization of the entirety of all molecules available in the cell nucleus, a general question remains open as to what extent chromatin organization might additionally be predetermined by the DNA sequence and, if so, if there are characteristic differences that distinguish typical regions involved in dysfunction-related aberrations from normal ones, since typical DNA breakpoint regions involved in disease-related chromosome aberrations are not randomly distributed along the DNA sequence. Highly conserved k-mer patterns in intronic and intergenic regions have been reported in eukaryotic genomes. In this article, we search and analyze regions deviating from average spectra (ReDFAS) of k-mer word frequencies in the human genome. This includes all assembled regions, e.g., telomeric, centromeric, genic as well as intergenic regions. A positive correlation between k-mer spectra and 3D contact frequencies, obtained exemplarily from given Hi-C datasets, has been found indicating a relation of ReDFAS to chromatin organization and interactions. We also searched and found correlations of known functional annotations, e.g., genes correlating with ReDFAS. Selected regions known to contain typical breakpoints on chromosomes 9 and 5 that are involved in cancer-related chromosomal aberrations appear to be enriched in ReDFAS. Since transposable elements like ALUs are often assigned as major players in 3D genome organization, we also studied their impact on our examples but could not find a correlation between ALU regions and breakpoints comparable to ReDFAS. Our findings might show that ReDFAS are associated with instable regions of the genome and regions with many chromatin contacts which is in line with current research indicating that chromatin loop anchor points lead to genomic instability.