HiCRep: assessing the reproducibility of Hi-C data using a stratum-adjusted correlation coefficient

Genome folding principles uncovered in condensin-depleted mitotic chromosomes

During mitosis, condensin activity is thought to interfere with interphase chromatin structures. To investigate genome folding principles in the absence of chromatin loop extrusion, we codepleted condensin I and condensin II, which triggered mitotic chromosome compartmentalization in ways similar to that in interphase. However, two distinct euchromatic compartments, indistinguishable in interphase, emerged upon condensin loss with different interaction preferences and dependencies on H3K27ac. Constitutive heterochromatin gradually self-aggregated and cocompartmentalized with facultative heterochromatin, contrasting with their separation during interphase. Notably, some cis-regulatory element contacts became apparent even in the absence of CTCF/cohesin-mediated structures. Heterochromatin protein 1 (HP1) proteins, which are thought to partition constitutive heterochromatin, were absent from mitotic chromosomes, suggesting, surprisingly, that constitutive heterochromatin can self-aggregate without HP1. Indeed, in cells traversing from M to G1 phase in the combined absence of HP1α, HP1β and HP1γ, constitutive heterochromatin compartments are normally re-established. In sum, condensin-deficient mitotic chromosomes illuminate forces of genome compartmentalization not identified in interphase cells.

Suv39h-catalyzed H3K9me3 is critical for euchromatic genome organization and the maintenance of gene transcription

H3K9me3-dependent heterochromatin is critical for the silencing of repeat-rich pericentromeric regions and also has key roles in repressing lineage-inappropriate protein-coding genes in differentiation and development. Here, we investigate the molecular consequences of heterochromatin loss in cells deficient in both SUV39H1 and SUV39H2 (Suv39DKO), the major mammalian histone methyltransferase enzymes that catalyze heterochromatic H3K9me3 deposition. We reveal a paradoxical repression of protein-coding genes in Suv39DKO cells, with these differentially expressed genes principally in euchromatic (Tn5-accessible, H3K4me3- and H3K27ac-marked) rather than heterochromatic (H3K9me3-marked) or polycomb (H3K27me3-marked) regions. Examination of the three-dimensional (3D) nucleome reveals that transcriptomic dysregulation occurs in euchromatic regions close to the nuclear periphery in 3D space. Moreover, this transcriptomic dysregulation is highly correlated with altered 3D genome organization in Suv39DKO cells. Together, our results suggest that the nuclear lamina-tethering of Suv39-dependent H3K9me3 domains provides an essential scaffold to support euchromatic genome organization and the maintenance of gene transcription for healthy cellular function.

A continuum of zinc finger transcription factor retention on native chromatin underlies dynamic genome organization

Transcription factor (TF) residence on chromatin translates into quantitative transcriptional or structural outcomes on genome. Commonly used formaldehyde crosslinking fixes TF-DNA interactions cumulatively and compromises the measured occupancy level. Here we mapped the occupancy level of global or individual zinc finger TFs like CTCF and MAZ, in the form of highly resolved footprints, on native chromatin. By incorporating reinforcing perturbation conditions, we established S-score, a quantitative metric to proxy the continuum of CTCF or MAZ retention across different motifs on native chromatin. The native chromatin-retained CTCF sites harbor sequence features within CTCF motifs better explained by S-score than the metrics obtained from other crosslinking or native assays. CTCF retention on native chromatin correlates with local SUMOylation level, and anti-correlates with transcriptional activity. The S-score successfully delineates the otherwise-masked differential stability of chromatin structures mediated by CTCF, or by MAZ independent of CTCF. Overall, our study established a paradigm continuum of TF retention across binding sites on native chromatin, explaining the dynamic genome organization.

Pervasive structural heterogeneity rewires glioblastoma chromosomes to sustain patient-specific transcriptional programs

Glioblastoma multiforme (GBM) encompasses brain malignancies marked by phenotypic and transcriptional heterogeneity thought to render these tumors aggressive, resistant to therapy, and inevitably recurrent. However, little is known about how the spatial organization of GBM genomes underlies this heterogeneity and its effects. Here, we compile a cohort of 28 patient-derived glioblastoma stem cell-like lines (GSCs) known to reflect the properties of their tumor-of-origin; six of these were primary-relapse tumor pairs from the same patient. We generate and analyze 5 kbp-resolution chromosome conformation capture (Hi-C) data from all GSCs to systematically map thousands of standalone and complex structural variants (SVs) and the multitude of neoloops arising as a result. By combining Hi-C, histone modification, and gene expression data with chromatin folding simulations, we explain how the pervasive, uneven, and idiosyncratic occurrence of neoloops sustains tumor-specific transcriptional programs via the formation of new enhancer-promoter contacts. We also show how even moderately recurrent neoloops can relate to patient-specific vulnerabilities. Together, our data provide a resource for dissecting GBM biology and heterogeneity, as well as for informing therapeutic approaches.

Exploring the dynamic three-dimensional chromatin architecture and transcriptional landscape in goose liver tissues underlying metabolic adaptations induced by a high-fat diet

BACKGROUND

Goose, descendants of migratory ancestors, have undergone extensive selective breeding, resulting in their remarkable ability to accumulate fat in the liver and exhibit a high tolerance for significant energy intake. As a result, goose offers an excellent model for studying obesity, metabolic disorders, and liver diseases in mammals. Although the impact of the three-dimensional arrangement of chromatin within the cell nucleus on gene expression and transcriptional regulation is widely acknowledged, the precise functions of chromatin architecture reorganization during fat deposition in goose liver tissues still need to be fully comprehended.

RESULTS

In this study, geese exhibited more pronounced changes in the liver index and triglyceride (TG) content following the consumption of the high-fat diet (HFD) than mice without significant signs of inflammation. Additionally, we performed comprehensive analyses on 10 goose liver tissues (5 HFD, 5 normal), including generating high-resolution maps of chromatin architecture, conducting whole-genome gene expression profiling, and identifying H3K27ac peaks in the livers of geese and mice subjected to the HFD. Our results unveiled a multiscale restructuring of chromatin architecture, encompassing Compartment A/B, topologically associated domains, and interactions between promoters and enhancers. The dynamism of the three-dimensional genome architecture, prompted by the HFD, assumed a pivotal role in the transcriptional regulation of crucial genes. Furthermore, we identified genes that regulate chromatin conformation changes, contributing to the metabolic adaptation process of lipid deposition and hepatic fat changes in geese in response to excessive energy intake. Moreover, we conducted a cross-species analysis comparing geese and mice exposed to the HFD, revealing unique characteristics specific to the goose liver compared to a mouse. These chromatin conformation changes help elucidate the observed characteristics of fat deposition and hepatic fat regulation in geese under conditions of excessive energy intake.

CONCLUSIONS

We examined the dynamic modifications in three-dimensional chromatin architecture and gene expression induced by an HFD in goose liver tissues. We conducted a cross-species analysis comparing that of mice. Our results contribute significant insights into the chromatin architecture of goose liver tissues, offering a novel perspective for investigating mammal liver diseases.

Unveiling Alterations of Epigenetic Modifications and Chromatin Architecture Leading to Lipid Metabolic Reprogramming during the Evolutionary Trastuzumab Adaptation of HER2-Positive Breast Cancer

Secondary trastuzumab resistance represents an evolutionary adaptation of HER2-positive breast cancer during anti-HER2 treatment. Most current studies have tended to prioritize HER2 and its associated signaling pathways, often overlooking broader but seemingly less relevant cellular processes, along with their associated genetic and epigenetic mechanisms. Here, transcriptome data is not only characterized but also examined epigenomic and 3D genome architecture information in both trastuzumab-sensitive and secondary-resistant breast cancer cells. The findings reveal that the global metabolic reprogramming associated with trastuzumab resistance may stem from genome-wide alterations in both histone modifications and chromatin structure. Specifically, the transcriptional activities of key genes involved in lipid metabolism appear to be regulated by variant promoter H3K27me3 and H3K4me3 modifications, as well as promoter-enhancer interactions. These discoveries offer valuable insights into how cancer cells adapt to anti-tumor drugs and have the potential to impact future diagnostic and treatment strategies.

CTPredictor: A comprehensive and robust framework for predicting cell types by integrating multi-scale features from single-cell Hi-C data

Single-cell Hi-C (scHi-C) has emerged as a powerful technology for deciphering cell-to-cell variability in three-dimensional (3D) chromatin organization, providing insights into genome-wide chromatin interactions and their correlation with cellular functions. Nevertheless, the accurate identification of cell types across different datasets remains a formidable challenge, hindering comprehensive investigations into genome structure. In response, we introduce CTPredictor, an innovative computational method that integrates multi-scale features to accurately predict cell types in various datasets. CTPredictor strategically incorporates three distinct feature sets, namely, small intra-domain contact probability (SICP), smoothed small intra-domain contact probability (SSICP), and smoothed bin contact probability (SBCP). The resulting fusion classification model significantly enhances the accuracy of cell type prediction based on single-cell Hi-C data (scHi-C). Rigorous benchmarking against established methods and three conventional machine learning approaches demonstrates the robust performance of CTPredictor, positioning it as an advanced tool for cell type prediction within scHi-C data. Beyond its prediction capabilities, CTPredictor holds promise in illuminating 3D genome structures and their functional significance across a wide array of biological processes.

Three-dimensional chromatin landscapes in somatotroph tumour

BACKGROUND

The three-dimensional (3D) genome architecture plays a critical role inregulating gene expression. However, the specific alterations in thisarchitecture within somatotroph tumors and their implications for gene expression remain largely unexplored.

METHODS

We employed Hi-C and RNA-seq analyses to compare the 3D genomic structures of somatotroph tumors with normal pituitary tissue. This comprehensive approachenabled the characterization of A/B compartments, topologically associateddomains (TADs), and chromatin loops, integrating these with gene expression patterns.

RESULTS

We observed a decrease in both the frequency of chromosomal interactions andthe size of TADs in tumor tissue compared to normal tissue. Conversely, the number of TADs and chromatin loops was found to be increased in tumors. Integrated analysis of Hi-C and RNA-seq data demonstrated that changes inhigher-order chromat in structure were associated with alterations in gene expression. Specifically, genes in A compartments showed higher density and increased expression relative to those in B compartments. Moreover, the weakand enhanced insulation boundaries were identified, and the associated genes were enriched in the Wnt/β-Catenin signaling pathway. We identified the gainedand lost loops in tumor and integrated these differences with transcriptional changes to examine the functional relevance of the identified loops. Notably, we observed an enhanced insulation boundary and a greater number of loops in the TCF7L2 gene region within tumors, which was accompanied by an upregulation of TCF7L2 expression. Subsequently, TCF7L2 expression was confirmed through qRT-PCR, and upregulated TCF7L2 prompted cell proliferation and growth hormone (GH) secretion in vitro.

CONCLUSION

Our results provide comprehensive 3D chromatin architecture maps of somatotroph tumors and offer a valuable resource for furthering the understanding of the underlying biology and mechanisms of gene expression regulation.

TEAD2 initiates ground-state pluripotency by mediating chromatin looping

The transition of mouse embryonic stem cells (ESCs) between serum/LIF and 2i(MEK and GSK3 kinase inhibitor)/LIF culture conditions serves as a valuable model for exploring the mechanisms underlying ground and confused pluripotent states. Regulatory networks comprising core and ancillary pluripotency factors drive the gene expression programs defining stable naïve pluripotency. In our study, we systematically screened factors essential for ESC pluripotency, identifying TEAD2 as an ancillary factor maintaining ground-state pluripotency in 2i/LIF ESCs and facilitating the transition from serum/LIF to 2i/LIF ESCs. TEAD2 exhibits increased binding to chromatin in 2i/LIF ESCs, targeting active chromatin regions to regulate the expression of 2i-specific genes. In addition, TEAD2 facilitates the expression of 2i-specific genes by mediating enhancer-promoter interactions during the serum/LIF to 2i/LIF transition. Notably, deletion of Tead2 results in reduction of a specific set of enhancer-promoter interactions without significantly affecting binding of chromatin architecture proteins, CCCTC-binding factor (CTCF), and Yin Yang 1 (YY1). In summary, our findings highlight a novel prominent role of TEAD2 in orchestrating higher-order chromatin structures of 2i-specific genes to sustain ground-state pluripotency.

Three-dimensional chromatin analysis reveals Sp1 as a mediator to program and reprogram HPV-host epigenetic architecture in cervical cancer

Human papillomavirus (HPV) is predominantly associated with HPV-related cancers, however, the precise mechanisms underlying the HPV-host epigenetic architectures in HPV carcinogenesis remain elusive. Here, we employed high-throughput chromosome conformation capture (Hi-C) to comprehensively map HPV16/18-host chromatin interactions. Our study identified the transcription factor Sp1 as a pivotal mediator in programming HPV-host interactions. By targeting Sp1, the active histone modifications (H3K27ac, H3K4me1, and H3K4me3) and the HPV-host chromatin interactions are reprogrammed, which leads to the downregulation of oncogenes located near the integration sites in both HPV (E6/E7) and the host genome (KLF5/MYC). Additionally, Sp1 inhibition led to the upregulation of immune checkpoint genes by reprogramming histone modifications in host cells. Notably, humanized patient-derived xenograft (PDX-HuHSC-NSG) models demonstrated that Sp1 inhibition promoted anti-PD-1 immunotherapy via remodeling the tumor immune microenvironment in cervical cancer. Moreover, single-cell transcriptomic analysis validated the enrichment of transcription factor Sp1 in epithelial cells of cervical cancer. In summary, our findings elucidate Sp1 as a key mediator involved in the programming and reprogramming of HPV-host epigenetic architecture. Inhibiting Sp1 with plicamycin may represent a promising therapeutic option for HPV-related carcinoma.

Proximal telomeric decompaction due to telomere shortening drives FOXC1-dependent myocardial senescence

Telomeres, TTAGGGn DNA repeat sequences located at the ends of eukaryotic chromosomes, play a pivotal role in aging and are targets of DNA damage response. Although we and others have demonstrated presence of short telomeres in genetic cardiomyopathic and heart failure cardiomyocytes, little is known about the role of telomere lengths in cardiomyocyte. Here, we demonstrate that in heart failure patient cardiomyocytes, telomeres are shortened compared to healthy controls. We generated isogenic human induced pluripotent stem cell derived cardiomyocytes (hiPSC-CMs) with short telomeres (sTL-CMs) and normal telomeres (nTL-CMs) as model. Compared to nTL-CMs, short telomeres result in cardiac dysfunction and expression of senescent markers. Using Hi-C and RNASeq, we observe that short telomeres induced TAD insulation decrease near telomeric ends and this correlated with a transcription upregulation in sTL-CMs. FOXC1, a key transcription factor involved in early cardiogenesis, was upregulated in sTL-CMs and its protein levels were negatively correlated with telomere lengths in heart failure patients. Overexpression of FOXC1 induced hiPSC-CM aging, mitochondrial and contractile dysfunction; knockdown of FOXC1 rescued these phenotypes. Overall, the work presented demonstrate that increased chromatin accessibility due to telomere shortening resulted in the induction of FOXC1-dependent expression network responsible for contractile dysfunction and myocardial senescence.

Building Haplotype-Resolved 3D Genome Maps of Chicken Skeletal Muscle

Haplotype-resolved 3D chromatin architecture related to allelic differences in avian skeletal muscle development has not been addressed so far, although chicken husbandry for meat consumption has been prevalent feature of cultures on every continent for more than thousands of years. Here, high-resolution Hi-C diploid maps (1.2-kb maximum resolution) are generated for skeletal muscle tissues in chicken across three developmental stages (embryonic day 15 to day 30 post-hatching). The sequence features governing spatial arrangement of chromosomes and characterize homolog pairing in the nucleus, are identified. Multi-scale characterization of chromatin reorganization between stages from myogenesis in the fetus to myofiber hypertrophy after hatching show concordant changes in transcriptional regulation by relevant signaling pathways. Further interrogation of parent-of-origin-specific chromatin conformation supported that genomic imprinting is absent in birds. This study also reveals promoter-enhancer interaction (PEI) differences between broiler and layer haplotypes in skeletal muscle development-related genes are related to genetic variation between breeds, however, only a minority of breed-specific variations likely contribute to phenotypic divergence in skeletal muscle potentially via allelic PEI rewiring. Beyond defining the haplotype-specific 3D chromatin architecture in chicken, this study provides a rich resource for investigating allelic regulatory divergence among chicken breeds.

The PTM profiling of CTCF reveals the regulation of 3D chromatin structure by O-GlcNAcylation

CCCTC-binding factor (CTCF), a ubiquitously expressed and highly conserved protein, is known to play a critical role in chromatin structure. Post-translational modifications (PTMs) diversify the functions of protein to regulate numerous cellular processes. However, the effects of PTMs on the genome-wide binding of CTCF and the organization of three-dimensional (3D) chromatin structure have not been fully understood. In this study, we uncovered the PTM profiling of CTCF and demonstrated that CTCF can be O-GlcNAcylated and arginine methylated. Functionally, we demonstrated that O-GlcNAcylation inhibits CTCF binding to chromatin. Meanwhile, deficiency of CTCF O-GlcNAcylation results in the disruption of loop domains and the alteration of chromatin loops associated with cellular development. Furthermore, the deficiency of CTCF O-GlcNAcylation increases the expression of developmental genes and negatively regulates maintenance and establishment of stem cell pluripotency. In conclusion, these results provide key insights into the role of PTMs for the 3D chromatin structure.

p53 rapidly restructures 3D chromatin organization to trigger a transcriptional response

Activation of the p53 tumor suppressor triggers a transcriptional program to control cellular response to stress. However, the molecular mechanisms by which p53 controls gene transcription are not completely understood. Here, we uncover the critical role of spatio-temporal genome architecture in this process. We demonstrate that p53 drives direct and indirect changes in genome compartments, topologically associating domains, and DNA loops prior to one hour of its activation, which escort the p53 transcriptional program. Focusing on p53-bound enhancers, we report 340 genes directly regulated by p53 over a median distance of 116 kb, with 74% of these genes not previously identified. Finally, we showcase that p53 controls transcription of distal genes through newly formed and pre-existing enhancer-promoter loops in a cohesin dependent manner. Collectively, our findings demonstrate a previously unappreciated architectural role of p53 as regulator at distinct topological layers and provide a reliable set of new p53 direct target genes that may help designs of cancer therapies.

Removing unwanted variation between samples in Hi-C experiments

Hi-C data are commonly normalized using single sample processing methods, with focus on comparisons between regions within a given contact map. Here, we aim to compare contact maps across different samples. We demonstrate that unwanted variation, of likely technical origin, is present in Hi-C data with replicates from different individuals, and that properties of this unwanted variation change across the contact map. We present band-wise normalization and batch correction, a method for normalization and batch correction of Hi-C data and show that it substantially improves comparisons across samples, including in a quantitative trait loci analysis as well as differential enrichment across cell types.

Multifaceted roles of cohesin in regulating transcriptional loops

Cohesin is required for chromatin loop formation. However, its precise role in regulating gene transcription remains largely unknown. We investigated the relationship between cohesin and RNA Polymerase II (RNAPII) using single-molecule mapping and live-cell imaging methods in human cells. Cohesin-mediated transcriptional loops were highly correlated with those of RNAPII and followed the direction of gene transcription. Depleting RAD21, a subunit of cohesin, resulted in the loss of long-range (>100 kb) loops between distal (super-)enhancers and promoters of cell-type-specific genes. By contrast, the short-range (<50 kb) loops were insensitive to RAD21 depletion and connected genes that are mostly housekeeping. This result explains why only a small fraction of genes are affected by the loss of long-range chromatin interactions due to cohesin depletion. Remarkably, RAD21 depletion appeared to up-regulate genes located in early initiation zones (EIZ) of DNA replication, and the EIZ signals were amplified drastically without RAD21. Our results revealed new mechanistic insights of cohesin's multifaceted roles in establishing transcriptional loops, preserving long-range chromatin interactions for cell-specific genes, and maintaining timely order of DNA replication.

Nuclear morphology is shaped by loop-extrusion programs

It is well established that neutrophils adopt malleable polymorphonuclear shapes to migrate through narrow interstitial tissue spaces1-3. However, how polymorphonuclear structures are assembled remains unknown4. Here we show that in neutrophil progenitors, halting loop extrusion-a motor-powered process that generates DNA loops by pulling in chromatin5-leads to the assembly of polymorphonuclear genomes. Specifically, we found that in mononuclear neutrophil progenitors, acute depletion of the loop-extrusion loading factor nipped-B-like protein (NIPBL) induced the assembly of horseshoe, banded, ringed and hypersegmented nuclear structures and led to a reduction in nuclear volume, mirroring what is observed during the differentiation of neutrophils. Depletion of NIPBL also induced cell-cycle arrest, activated a neutrophil-specific gene program and conditioned a loss of interactions across topologically associating domains to generate a chromatin architecture that resembled that of differentiated neutrophils. Removing NIPBL resulted in enrichment for mega-loops and interchromosomal hubs that contain genes associated with neutrophil-specific enhancer repertoires and an inflammatory gene program. On the basis of these observations, we propose that in neutrophil progenitors, loop-extrusion programs produce lineage-specific chromatin architectures that permit the packing of chromosomes into geometrically confined lobular structures. Our data also provide a blueprint for the assembly of polymorphonuclear structures, and point to the possibility of engineering de novo nuclear shapes to facilitate the migration of effector cells in densely populated tumorigenic environments.

KDM3B inhibitors disrupt the oncogenic activity of PAX3-FOXO1 in fusion-positive rhabdomyosarcoma

Fusion-positive rhabdomyosarcoma (FP-RMS) is an aggressive pediatric sarcoma driven primarily by the PAX3-FOXO1 fusion oncogene, for which therapies targeting PAX3-FOXO1 are lacking. Here, we screen 62,643 compounds using an engineered cell line that monitors PAX3-FOXO1 transcriptional activity identifying a hitherto uncharacterized compound, P3FI-63. RNA-seq, ATAC-seq, and docking analyses implicate histone lysine demethylases (KDMs) as its targets. Enzymatic assays confirm the inhibition of multiple KDMs with the highest selectivity for KDM3B. Structural similarity search of P3FI-63 identifies P3FI-90 with improved solubility and potency. Biophysical binding of P3FI-90 to KDM3B is demonstrated using NMR and SPR. P3FI-90 suppresses the growth of FP-RMS in vitro and in vivo through downregulating PAX3-FOXO1 activity, and combined knockdown of KDM3B and KDM1A phenocopies P3FI-90 effects. Thus, we report KDM inhibitors P3FI-63 and P3FI-90 with the highest specificity for KDM3B. Their potent suppression of PAX3-FOXO1 activity indicates a possible therapeutic approach for FP-RMS and other transcriptionally addicted cancers.

Posterior inference of Hi-C contact frequency through sampling

Hi-C is one of the most widely used approaches to study three-dimensional genome conformations. Contacts captured by a Hi-C experiment are represented in a contact frequency matrix. Due to the limited sequencing depth and other factors, Hi-C contact frequency matrices are only approximations of the true interaction frequencies and are further reported without any quantification of uncertainty. Hence, downstream analyses based on Hi-C contact maps (e.g., TAD and loop annotation) are themselves point estimations. Here, we present the Hi-C interaction frequency sampler (HiCSampler) that reliably infers the posterior distribution of the interaction frequency for a given Hi-C contact map by exploiting dependencies between neighboring loci. Posterior predictive checks demonstrate that HiCSampler can infer highly predictive chromosomal interaction frequency. Summary statistics calculated by HiCSampler provide a measurement of the uncertainty for Hi-C experiments, and samples inferred by HiCSampler are ready for use by most downstream analysis tools off the shelf and permit uncertainty measurements in these analyses without modifications.

SATB2 organizes the 3D genome architecture of cognition in cortical neurons

The DNA-binding protein SATB2 is genetically linked to human intelligence. We studied its influence on the three-dimensional (3D) epigenome by mapping chromatin interactions and accessibility in control versus SATB2-deficient cortical neurons. We find that SATB2 affects the chromatin looping between enhancers and promoters of neuronal-activity-regulated genes, thus influencing their expression. It also alters A/B compartments, topologically associating domains, and frequently interacting regions. Genes linked to SATB2-dependent 3D genome changes are implicated in highly specialized neuronal functions and contribute to cognitive ability and risk for neuropsychiatric and neurodevelopmental disorders. Non-coding DNA regions with a SATB2-dependent structure are enriched for common variants associated with educational attainment, intelligence, and schizophrenia. Our data establish SATB2 as a cell-type-specific 3D genome modulator, which operates both independently and in cooperation with CCCTC-binding factor (CTCF) to set up the chromatin landscape of pyramidal neurons for cognitive processes.

Orchestrating chromosome conformation capture analysis with Bioconductor

Genome-wide chromatin conformation capture assays provide formidable insights into the spatial organization of genomes. However, due to the complexity of the data structure, their integration in multi-omics workflows remains challenging. We present data structures, computational methods and visualization tools available in Bioconductor to investigate Hi-C, micro-C and other 3C-related data, in R. An online book ( https://bioconductor.org/books/OHCA/ ) further provides prospective end users with a number of workflows to process, import, analyze and visualize any type of chromosome conformation capture data.

Bottom-up data integration in polymer models of chromatin organization

Cellular functions crucially depend on the precise execution of complex biochemical reactions taking place on the chromatin fiber in the tightly packed environment of the cell nucleus. Despite the availability of large datasets probing this process from multiple angles, bottom-up frameworks that allow the incorporation of the sequence-specific nature of biochemistry in a unified model of 3D chromatin structure remain scarce. Here, we propose Sequence-Enhanced Magnetic Polymer (SEMPER), a novel stochastic polymer model that naturally incorporates observational data about sequence-driven biochemical processes, such as binding of transcription factor proteins, in a 3D model of chromatin structure. We introduce a novel approximate Bayesian algorithm to quantify a posteriori the relative importance of various factors, including the polymeric nature of DNA, in determining chromatin epigenetic state, thus providing a transparent way to generate biological hypotheses. Although accurate prediction of contact frequencies (a problem already extensively studied in the literature) is not our main aim, as a by-product of the inference procedure and without additional input from the genome 3D structure, our model can predict with reasonable accuracy some notable and nontrivial conformational features of chromatin folding within the nucleus. Our work highlights the importance of introducing physically realistic statistical models for predicting chromatin states from epigenetic data and opens the way to a new class of more systematic approaches to interpreting epigenomic data.

SnapFISH-IMPUTE: an imputation method for multiplexed DNA FISH data

Chromatin spatial organization plays a crucial role in gene regulation. Recently developed and prospering multiplexed DNA FISH technologies enable direct visualization of chromatin conformation in nucleus. However, incomplete data caused by limited detection efficiency can substantially complicate and impair downstream analysis. Here, we present SnapFISH-IMPUTE that imputes missing values in multiplexed DNA FISH data. Analysis on multiple published datasets shows that the proposed method preserves the distribution of pairwise distances between imaging loci, and the imputed chromatin conformations are indistinguishable from the observed conformations. Additionally, imputation greatly improves downstream analyses such as identifying enhancer-promoter loops and clustering cells into distinct cell types. SnapFISH-IMPUTE is freely available at https://github.com/hyuyu104/SnapFISH-IMPUTE.

DiffDomain enables identification of structurally reorganized topologically associating domains

Topologically associating domains (TADs) are critical structural units in three-dimensional genome organization of mammalian genome. Dynamic reorganizations of TADs between health and disease states are associated with essential genome functions. However, computational methods for identifying reorganized TADs are still in the early stages of development. Here, we present DiffDomain, an algorithm leveraging high-dimensional random matrix theory to identify structurally reorganized TADs using high-throughput chromosome conformation capture (Hi-C) contact maps. Method comparison using multiple real Hi-C datasets reveals that DiffDomain outperforms alternative methods for false positive rates, true positive rates, and identifying a new subtype of reorganized TADs. Applying DiffDomain to Hi-C data from different cell types and disease states demonstrates its biological relevance. Identified reorganized TADs are associated with structural variations and epigenomic changes such as changes in CTCF binding sites. By applying to a single-cell Hi-C data from mouse neuronal development, DiffDomain can identify reorganized TADs between cell types with reasonable reproducibility using pseudo-bulk Hi-C data from as few as 100 cells per condition. Moreover, DiffDomain reveals differential cell-to-population variability and heterogeneous cell-to-cell variability in TADs. Therefore, DiffDomain is a statistically sound method for better comparative analysis of TADs using both Hi-C and single-cell Hi-C data.

ChIPr: accurate prediction of cohesin-mediated 3D genome organization from 2D chromatin features

The three-dimensional genome organization influences diverse nuclear processes. Here we present Chromatin Interaction Predictor (ChIPr), a suite of regression models based on deep neural networks, random forest, and gradient boosting to predict cohesin-mediated chromatin interaction strength between any two loci in the genome. The predictions of ChIPr correlate well with ChIA-PET data in four cell lines. The standard ChIPr model requires three experimental inputs: ChIP-Seq signals for RAD21, H3K27ac, and H3K27me3 but works well with just RAD21 signal. Integrative analysis reveals novel insights into the role of CTCF motif, its orientation, and CTCF binding on cohesin-mediated chromatin interactions.

Stage-specific dynamic reorganization of genome topology shapes transcriptional neighborhoods in developing human retinal organoids

We have generated a high-resolution Hi-C map of developing human retinal organoids to elucidate spatiotemporal dynamics of genomic architecture and its relationship with gene expression patterns. We demonstrate progressive stage-specific alterations in DNA topology and correlate these changes with transcription of cell-type-restricted gene markers during retinal differentiation. Temporal Hi-C reveals a shift toward A compartment for protein-coding genes and B compartment for non-coding RNAs, displaying high and low expression, respectively. Notably, retina-enriched genes are clustered near lost boundaries of topologically associated domains (TADs), and higher-order assemblages (i.e., TAD cliques) localize in active chromatin regions with binding sites for eye-field transcription factors. These genes gain chromatin contacts at their transcription start site as organoid differentiation proceeds. Our study provides a global view of chromatin architecture dynamics associated with diversification of cell types during retinal development and serves as a foundational resource for in-depth functional investigations of retinal developmental traits.

3D genome remodeling and homologous pairing during meiotic prophase of mouse oogenesis and spermatogenesis

During meiosis, the chromatin and transcriptome undergo prominent switches. Although recent studies have explored the genome reorganization during spermatogenesis, the chromatin remodeling in oogenesis and characteristics of homologous pairing remain largely elusive. We comprehensively compared chromatin structures and transcriptomes at successive substages of meiotic prophase in both female and male mice using low-input high-through chromosome conformation capture (Hi-C) and RNA sequencing (RNA-seq). Compartments and topologically associating domains (TADs) gradually disappeared and slowly recovered in both sexes. We found that homologs adopted different sex-conserved pairing strategies prior to and after the leptotene-to-zygotene transition, changing from long interspersed nuclear element (LINE)-enriched compartments B to short interspersed nuclear element (SINE)-enriched compartments A. We complemented marker genes and predicted the sex-specific meiotic sterile genes for each substage. This study provides valuable insights into the similarities and distinctions between sexes in chromosome architecture, homologous pairing, and transcriptome during meiotic prophase of both oogenesis and spermatogenesis.

Joint epigenome profiling reveals cell-type-specific gene regulatory programmes in human cortical organoids

Gene expression is regulated by multiple epigenetic mechanisms, which are coordinated in development and disease. However, current multiomics methods are frequently limited to one or two modalities at a time, making it challenging to obtain a comprehensive gene regulatory signature. Here, we describe a method-3D genome, RNA, accessibility and methylation sequencing (3DRAM-seq)-that simultaneously interrogates spatial genome organization, chromatin accessibility and DNA methylation genome-wide and at high resolution. We combine 3DRAM-seq with immunoFACS and RNA sequencing in cortical organoids to map the cell-type-specific regulatory landscape of human neural development across multiple epigenetic layers. Finally, we apply a massively parallel reporter assay to profile cell-type-specific enhancer activity in organoids and to functionally assess the role of key transcription factors for human enhancer activation and function. More broadly, 3DRAM-seq can be used to profile the multimodal epigenetic landscape in rare cell types and different tissues.

Multi-level 3D genome organization deteriorates 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) remain largely stable, 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 leads to the formation of de novo enhancer contacts and activation of the oncogene, illustrating how structural variants can interact with the 3D genome to drive oncogenic states. In summary, our findings provide evidence for the degradation of genome organization at multiple scales during breast cancer progression revealing novel relationships between genome 3D structure and oncogenic processes.

Improved quality metrics for association and reproducibility in chromatin accessibility data using mutual information

BACKGROUND

Correlation metrics are widely utilized in genomics analysis and often implemented with little regard to assumptions of normality, homoscedasticity, and independence of values. This is especially true when comparing values between replicated sequencing experiments that probe chromatin accessibility, such as assays for transposase-accessible chromatin via sequencing (ATAC-seq). Such data can possess several regions across the human genome with little to no sequencing depth and are thus non-normal with a large portion of zero values. Despite distributed use in the epigenomics field, few studies have evaluated and benchmarked how correlation and association statistics behave across ATAC-seq experiments with known differences or the effects of removing specific outliers from the data. Here, we developed a computational simulation of ATAC-seq data to elucidate the behavior of correlation statistics and to compare their accuracy under set conditions of reproducibility.

RESULTS

Using these simulations, we monitored the behavior of several correlation statistics, including the Pearson's R and Spearman's [Formula: see text] coefficients as well as Kendall's [Formula: see text] and Top-Down correlation. We also test the behavior of association measures, including the coefficient of determination R[Formula: see text], Kendall's W, and normalized mutual information. Our experiments reveal an insensitivity of most statistics, including Spearman's [Formula: see text], Kendall's [Formula: see text], and Kendall's W, to increasing differences between simulated ATAC-seq replicates. The removal of co-zeros (regions lacking mapped sequenced reads) between simulated experiments greatly improves the estimates of correlation and association. After removing co-zeros, the R[Formula: see text] coefficient and normalized mutual information display the best performance, having a closer one-to-one relationship with the known portion of shared, enhanced loci between simulated replicates. When comparing values between experimental ATAC-seq data using a random forest model, mutual information best predicts ATAC-seq replicate relationships.

CONCLUSIONS

Collectively, this study demonstrates how measures of correlation and association can behave in epigenomics experiments. We provide improved strategies for quantifying relationships in these increasingly prevalent and important chromatin accessibility assays.

Lineage-specific 3D genome organization is assembled at multiple scales by IKAROS

A generic level of chromatin organization generated by the interplay between cohesin and CTCF suffices to limit promiscuous interactions between regulatory elements, but a lineage-specific chromatin assembly that supersedes these constraints is required to configure the genome to guide gene expression changes that drive faithful lineage progression. Loss-of-function approaches in B cell precursors show that IKAROS assembles interactions across megabase distances in preparation for lymphoid development. Interactions emanating from IKAROS-bound enhancers override CTCF-imposed boundaries to assemble lineage-specific regulatory units built on a backbone of smaller invariant topological domains. Gain of function in epithelial cells confirms IKAROS' ability to reconfigure chromatin architecture at multiple scales. Although the compaction of the Igκ locus required for genome editing represents a function of IKAROS unique to lymphocytes, the more general function to preconfigure the genome to support lineage-specific gene expression and suppress activation of extra-lineage genes provides a paradigm for lineage restriction.

Cellular reprogramming is driven by widespread rewiring of promoter-enhancer interactions

BACKGROUND

Long-range interactions between promoters and cis-regulatory elements, such as enhancers, play critical roles in gene regulation. However, the role of three-dimensional (3D) chromatin structure in orchestrating changes in transcriptional regulation during direct cell reprogramming is not fully understood.

RESULTS

Here, we performed integrated analyses of chromosomal architecture, epigenetics, and gene expression using Hi-C, promoter Capture Hi-C (PCHi-C), ChIP-seq, and RNA-seq during trans-differentiation of Pre-B cells into macrophages with a β-estradiol inducible C/EBPαER transgene. Within 1h of β-estradiol induction, C/EBPα translocated from the cytoplasm to the nucleus, binding to thousands of promoters and putative regulatory elements, resulting in the downregulation of Pre-B cell-specific genes and induction of macrophage-specific genes. Hi-C results were remarkably consistent throughout trans-differentiation, revealing only a small number of TAD boundary location changes, and A/B compartment switches despite significant changes in the expression of thousands of genes. PCHi-C revealed widespread changes in promoter-anchored loops with decreased interactions in parallel with decreased gene expression, and new and increased promoter-anchored interactions in parallel with increased expression of macrophage-specific genes.

CONCLUSIONS

Overall, our data demonstrate that C/EBPα-induced trans-differentiation involves few changes in genome architecture at the level of TADs and A/B compartments, in contrast with widespread reorganization of thousands of promoter-anchored loops in association with changes in gene expression and cell identity.

Genome folding principles revealed in condensin-depleted mitotic chromosomes

During mitosis, condensin activity interferes with interphase chromatin structures. Here, we generated condensin-free mitotic chromosomes to investigate genome folding principles. Co- depletion of condensin I and II, but neither alone, triggered mitotic chromosome compartmentalization in ways that differ from interphase. Two distinct euchromatic compartments, indistinguishable in interphase, rapidly emerged upon condensin loss with different interaction preferences and dependence on H3K27ac. Constitutive heterochromatin gradually self-aggregated and co-compartmentalized with the facultative heterochromatin, contrasting with their separation during interphase. While topologically associating domains (TADs) and CTCF/cohesin mediated structural loops remained undetectable, cis-regulatory element contacts became apparent, providing an explanation for their quick re-establishment during mitotic exit. HP1 proteins, which are thought to partition constitutive heterochromatin, were absent from mitotic chromosomes, suggesting, surprisingly, that constitutive heterochromatin can self-aggregate without HP1. Indeed, in cells traversing from M- to G1-phase in the combined absence of HP1α, HP1Π and HP1γ, re-established constitutive heterochromatin compartments normally. In sum, "clean-slate" condensin-deficient mitotic chromosomes illuminate mechanisms of genome compartmentalization not revealed in interphase cells.

Lineage specific 3D genome structure in the adult human brain and neurodevelopmental changes in the chromatin interactome

The human brain is a complex organ comprised of distinct cell types, and the contribution of the 3D genome to lineage specific gene expression remains poorly understood. To decipher cell type specific genome architecture, and characterize fine scale changes in the chromatin interactome across neural development, we compared the 3D genome of the human fetal cortical plate to that of neurons and glia isolated from the adult prefrontal cortex. We found that neurons have weaker genome compartmentalization compared to glia, but stronger TADs, which emerge during fetal development. Furthermore, relative to glia, the neuronal genome shifts more strongly towards repressive compartments. Neurons have differential TAD boundaries that are proximal to active promoters involved in neurodevelopmental processes. CRISPRi on CNTNAP2 in hIPSC-derived neurons reveals that transcriptional inactivation correlates with loss of insulation at the differential boundary. Finally, re-wiring of chromatin loops during neural development is associated with transcriptional and functional changes. Importantly, differential loops in the fetal cortex are associated with autism GWAS loci, suggesting a neuropsychiatric disease mechanism affecting the chromatin interactome. Furthermore, neural development involves gaining enhancer-promoter loops that upregulate genes that control synaptic activity. Altogether, our study provides multi-scale insights on the 3D genome in the human brain.

Three-dimensional genome architecture coordinates key regulators of lineage specification in mammary epithelial cells

Although lineage-specific genes have been identified in the mammary gland, little is known about the contribution of the 3D genome organization to gene regulation in the epithelium. Here, we describe the chromatin landscape of the three major epithelial subsets through integration of long- and short-range chromatin interactions, accessibility, histone modifications, and gene expression. While basal genes display exquisite lineage specificity via distal enhancers, luminal-specific genes show widespread promoter priming in basal cells. Cell specificity in luminal progenitors is largely mediated through extensive chromatin interactions with super-enhancers in gene-body regions in addition to interactions with polycomb silencer elements. Moreover, lineage-specific transcription factors appear to be controlled through cell-specific chromatin interactivity. Finally, chromatin accessibility rather than interactivity emerged as a defining feature of the activation of quiescent basal stem cells. This work provides a comprehensive resource for understanding the role of higher-order chromatin interactions in cell-fate specification and differentiation in the adult mouse mammary gland.

Exploring the Roles of RNAs in Chromatin Architecture Using Deep Learning

Recent studies have highlighted the impact of both transcription and transcripts on 3D genome organization, particularly its dynamics. Here, we propose a deep learning framework, called AkitaR, that leverages both genome sequences and genome-wide RNA-DNA interactions to investigate the roles of chromatin-associated RNAs (caRNAs) on genome folding in HFFc6 cells. In order to disentangle the cis- and trans-regulatory roles of caRNAs, we compared models with nascent transcripts, trans-located caRNAs, open chromatin data, or DNA sequence alone. Both nascent transcripts and trans-located caRNAs improved the models' predictions, especially at cell-type-specific genomic regions. Analyses of feature importance scores revealed the contribution of caRNAs at TAD boundaries, chromatin loops and nuclear sub-structures such as nuclear speckles and nucleoli to the models' predictions. Furthermore, we identified non-coding RNAs (ncRNAs) known to regulate chromatin structures, such as MALAT1 and NEAT1, as well as several novel RNAs, RNY5, RPPH1, POLG-DT and THBS1-IT, that might modulate chromatin architecture through trans-interactions in HFFc6. Our modeling also suggests that transcripts from Alus and other repetitive elements may facilitate chromatin interactions through trans R-loop formation. Our findings provide new insights and generate testable hypotheses about the roles of caRNAs in shaping chromatin organization.

Genome-wide chromatin interaction profiling reveals a vital role of super-enhancers and rearrangements in host enhancer contacts during BmNPV infection

As influential regulatory elements in the genome, enhancers control gene expression under specific cellular conditions, and such connections are dynamic under different conditions. However, because of the lack of a genome-wide enhancer-gene connection map, the roles and regulatory pattern of enhancers were poorly investigated in insects, and the dynamic changes of enhancer contacts and functions under different conditions remain elusive. Here, combining Hi-C, ATAC-seq, and H3K27ac ChIP-seq data, we generate the genome-wide enhancer-gene map of silkworm and identify super-enhancers with a role in regulating the expression of vital genes related to cell state maintenance through a sophisticated interaction network. Additionally, a radical attenuation of chromatin interactions is found after infection of Bombyx mori nucleopolyhedrovirus (BmNPV), the main pathogen of silkworm, which directly rearranges the enhancer contacts. Such a rearrangement disturbs the intrinsic enhancer-gene connections in several antiviral genes, resulting in reduced expression of these genes, which accelerates viral infection. Overall, our results reveal the regulatory role of super-enhancers and shed new light on the mechanisms and dynamic changes of the genome-wide enhancer regulatory network.

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 found 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.

HiC4D: forecasting spatiotemporal Hi-C data with residual ConvLSTM

The Hi-C experiments have been extensively used for the studies of genomic structures. In the last few years, spatiotemporal Hi-C has largely contributed to the investigation of genome dynamic reorganization. However, computationally modeling and forecasting spatiotemporal Hi-C data still have not been seen in the literature. We present HiC4D for dealing with the problem of forecasting spatiotemporal Hi-C data. We designed and benchmarked a novel network and named it residual ConvLSTM (ResConvLSTM), which is a combination of residual network and convolutional long short-term memory (ConvLSTM). We evaluated our new ResConvLSTM networks and compared them with the other five methods, including a naïve network (NaiveNet) that we designed as a baseline method and four outstanding video-prediction methods from the literature: ConvLSTM, spatiotemporal LSTM (ST-LSTM), self-attention LSTM (SA-LSTM) and simple video prediction (SimVP). We used eight different spatiotemporal Hi-C datasets for the blind test, including two from mouse embryogenesis, one from somatic cell nuclear transfer (SCNT) embryos, three embryogenesis datasets from different species and two non-embryogenesis datasets. Our evaluation results indicate that our ResConvLSTM networks almost always outperform the other methods on the eight blind-test datasets in terms of accurately predicting the Hi-C contact matrices at future time-steps. Our benchmarks also indicate that all of the methods that we benchmarked can successfully recover the boundaries of topologically associating domains called on the experimental Hi-C contact matrices. Taken together, our benchmarks suggest that HiC4D is an effective tool for predicting spatiotemporal Hi-C data. HiC4D is publicly available at both http://dna.cs.miami.edu/HiC4D/ and https://github.com/zwang-bioinformatics/HiC4D/.

Context-dependent perturbations in chromatin folding and the transcriptome by cohesin and related factors

Cohesin regulates gene expression through context-specific chromatin folding mechanisms such as enhancer-promoter looping and topologically associating domain (TAD) formation by cooperating with factors such as cohesin loaders and the insulation factor CTCF. We developed a computational workflow to explore how three-dimensional (3D) structure and gene expression are regulated collectively or individually by cohesin and related factors. The main component is CustardPy, by which multi-omics datasets are compared systematically. To validate our methodology, we generated 3D genome, transcriptome, and epigenome data before and after depletion of cohesin and related factors and compared the effects of depletion. We observed diverse effects on the 3D genome and transcriptome, and gene expression changes were correlated with the splitting of TADs caused by cohesin loss. We also observed variations in long-range interactions across TADs, which correlated with their epigenomic states. These computational tools and datasets will be valuable for 3D genome and epigenome studies.

PfMORC protein regulates chromatin accessibility and transcriptional repression in the human malaria parasite, P. falciparum

The environmental challenges the human malaria parasite, Plasmodium falciparum, faces during its progression into its various lifecycle stages warrant the use of effective and highly regulated access to chromatin for transcriptional regulation. Microrchidia (MORC) proteins have been implicated in DNA compaction and gene silencing across plant and animal kingdoms. Accumulating evidence has shed light into the role MORC protein plays as a transcriptional switch in apicomplexan parasites. In this study, using CRISPR/Cas9 genome editing tool along with complementary molecular and genomics approaches, we demonstrate that PfMORC not only modulates chromatin structure and heterochromatin formation throughout the parasite erythrocytic cycle, but is also essential to the parasite survival. Chromatin immunoprecipitation followed by deep sequencing (ChIP-seq) experiments suggest that PfMORC binds to not only sub-telomeric regions and genes involved in antigenic variation but is also most likely a key modulator of stage transition. Protein knockdown experiments followed by chromatin conformation capture (Hi-C) studies indicate that downregulation of PfMORC induces the collapse of the parasite heterochromatin structure leading to its death. All together these findings confirm that PfMORC plays a crucial role in chromatin structure and gene regulation, validating this factor as a strong candidate for novel antimalarial strategies.

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.

Simultaneous profiling of chromatin architecture and transcription in single cells

The three-dimensional structure of chromatin plays a crucial role in development and disease, both of which are associated with transcriptional changes. However, given the heterogeneity in single-cell chromatin architecture and transcription, the regulatory relationship between the three-dimensional chromatin structure and gene expression is difficult to explain based on bulk cell populations. Here we develop a single-cell, multimodal, omics method allowing the simultaneous detection of chromatin architecture and messenger RNA expression by sequencing (single-cell transcriptome sequencing (scCARE-seq)). Applying scCARE-seq to examine chromatin architecture and transcription from 2i to serum single mouse embryonic stem cells, we observe improved separation of cell clusters compared with single-cell chromatin conformation capture. In addition, after defining the cell-cycle phase of each cell through chromatin architecture extracted by scCARE-seq, we find that periodic changes in chromatin architecture occur in parallel with transcription during the cell cycle. These findings highlight the potential of scCARE-seq to facilitate comprehensive analyses that may boost our understanding of chromatin architecture and transcription in the same single cell.

Evaluating the role of the nuclear microenvironment in gene function by population-based modeling

The nuclear folding of chromosomes relative to nuclear bodies is an integral part of gene function. Here, we demonstrate that population-based modeling-from ensemble Hi-C data-provides a detailed description of the nuclear microenvironment of genes and its role in gene function. We define the microenvironment by the subnuclear positions of genomic regions with respect to nuclear bodies, local chromatin compaction, and preferences in chromatin compartmentalization. These structural descriptors are determined in single-cell models, thereby revealing the structural variability between cells. We demonstrate that the microenvironment of a genomic region is linked to its functional potential in gene transcription, replication, and chromatin compartmentalization. Some chromatin regions feature a strong preference for a single microenvironment, due to association with specific nuclear bodies in most cells. Other chromatin shows high structural variability, which is a strong indicator of functional heterogeneity. Moreover, we identify specialized nuclear microenvironments, which distinguish chromatin in different functional states and reveal a key role of nuclear speckles in chromosome organization. We demonstrate that our method produces highly predictive three-dimensional genome structures, which accurately reproduce data from a variety of orthogonal experiments, thus considerably expanding the range of Hi-C data analysis.

High-throughput Oligopaint screen identifies druggable 3D genome regulators

The human genome functions as a three-dimensional chromatin polymer, driven by a complex collection of chromosome interactions1-3. Although the molecular rules governing these interactions are being quickly elucidated, relatively few proteins regulating this process have been identified. Here, to address this gap, we developed high-throughput DNA or RNA labelling with optimized Oligopaints (HiDRO)-an automated imaging pipeline that enables the quantitative measurement of chromatin interactions in single cells across thousands of samples. By screening the human druggable genome, we identified more than 300 factors that influence genome folding during interphase. Among these, 43 genes were validated as either increasing or decreasing interactions between topologically associating domains. Our findings show that genetic or chemical inhibition of the ubiquitous kinase GSK3A leads to increased long-range chromatin looping interactions in a genome-wide and cohesin-dependent manner. These results demonstrate the importance of GSK3A signalling in nuclear architecture and the use of HiDRO for identifying mechanisms of spatial genome organization.

Single-cell Hi-C data enhancement with deep residual and generative adversarial networks

MOTIVATION

The spatial genome organization of a eukaryotic cell is important for its function. The development of single-cell technologies for probing the 3D genome conformation, especially single-cell chromosome conformation capture techniques, has enabled us to understand genome function better than before. However, due to extreme sparsity and high noise associated with single-cell Hi-C data, it is still difficult to study genome structure and function using the HiC-data of one single cell.

RESULTS

In this work, we developed a deep learning method ScHiCEDRN based on deep residual networks and generative adversarial networks for the imputation and enhancement of Hi-C data of a single cell. In terms of both image evaluation and Hi-C reproducibility metrics, ScHiCEDRN outperforms the four deep learning methods (DeepHiC, HiCPlus, HiCSR, and Loopenhance) on enhancing the raw single-cell Hi-C data of human and Drosophila. The experiments also show that it can generate single-cell Hi-C data more suitable for identifying topologically associating domain boundaries and reconstructing 3D chromosome structures than the existing methods. Moreover, ScHiCEDRN's performance generalizes well across different single cells and cell types, and it can be applied to improving population Hi-C data.

AVAILABILITY AND IMPLEMENTATION

The source code of ScHiCEDRN is available at the GitHub repository: https://github.com/BioinfoMachineLearning/ScHiCEDRN.

Multi-omic analysis reveals dynamic changes of three-dimensional chromatin architecture during T cell differentiation

Cell differentiation results in widespread changes in transcriptional programs as well as multi-level remodeling of three-dimensional genome architecture. Nonetheless, few synthetically investigate the chromatin higher-order landscapes in different T helper (Th) cells. Using RNA-Seq, ATAC-Seq and Hi-C assays, we characterize dynamic changes in chromatin organization at different levels during Naive CD4⁺ T cells differentiation into T helper 17 (Th17) and T helper 1 (Th1) cells. Upon differentiation, we observe decreased short-range and increased extra-long-range chromatin interactions. Although there is no apparent global switch in the A/B compartments, Th cells display the weaker compartmentalization. A portion of topologically associated domains are rearranged. Furthermore, we identify cell-type specific enhancer-promoter loops, many of which are associated with functional genes in Th cells, such as Rorc facilitating Th17 differentiation and Hif1a responding to intracellular oxygen levels in Th1. Taken together, these results uncover the general patterns of chromatin reorganization and epigenetic landscapes of gene regulation during T helper cell differentiation.

Multidimensional profiling reveals GATA1-modulated stage-specific chromatin states and functional associations during human erythropoiesis

Mammalian erythroid development can be divided into three stages: hematopoietic stem and progenitor cell (HSPC), erythroid progenitor (Ery-Pro), and erythroid precursor (Ery-Pre). However, the mechanisms by which the 3D genome changes to establish the stage-specific transcription programs that are critical for erythropoiesis remain unclear. Here, we analyze the chromatin landscape at multiple levels in defined populations from primary human erythroid culture. While compartments and topologically associating domains remain largely unchanged, ∼50% of H3K27Ac-marked enhancers are dynamic in HSPC versus Ery-Pre. The enhancer anchors of enhancer-promoter loops are enriched for occupancy of respective stage-specific transcription factors (TFs), indicating these TFs orchestrate the enhancer connectome rewiring. The master TF of erythropoiesis, GATA1, is found to occupy most erythroid gene promoters at the Ery-Pro stage, and mediate conspicuous local rewiring through acquiring binding at the distal regions in Ery-Pre, promoting productive erythroid transcription output. Knocking out GATA1 binding sites precisely abrogates local rewiring and corresponding gene expression. Interestingly, knocking down GATA1 can transiently revert the cell state to an earlier stage and prolong the window of progenitor state. This study reveals mechanistic insights underlying chromatin rearrangements during development by integrating multidimensional chromatin landscape analyses to associate with transcription output and cellular states.

Improving comparative analyses of Hi-C data via contrastive self-supervised learning

Hi-C is a widely applied chromosome conformation capture (3C)-based technique, which has produced a large number of genomic contact maps with high sequencing depths for a wide range of cell types, enabling comprehensive analyses of the relationships between biological functionalities (e.g. gene regulation and expression) and the three-dimensional genome structure. Comparative analyses play significant roles in Hi-C data studies, which are designed to make comparisons between Hi-C contact maps, thus evaluating the consistency of replicate Hi-C experiments (i.e. reproducibility measurement) and detecting statistically differential interacting regions with biological significance (i.e. differential chromatin interaction detection). However, due to the complex and hierarchical nature of Hi-C contact maps, it remains challenging to conduct systematic and reliable comparative analyses of Hi-C data. Here, we proposed sslHiC, a contrastive self-supervised representation learning framework, for precisely modeling the multi-level features of chromosome conformation and automatically producing informative feature embeddings for genomic loci and their interactions to facilitate comparative analyses of Hi-C contact maps. Comprehensive computational experiments on both simulated and real datasets demonstrated that our method consistently outperformed the state-of-the-art baseline methods in providing reliable measurements of reproducibility and detecting differential interactions with biological meanings.

Schizophrenia Risk Mapping and Functional Engineering of the 3D Genome in Three Neuronal Subtypes

Common variants associated with schizophrenia are concentrated in non-coding regulatory sequences, but their precise target genes are context-dependent and impacted by cell-type-specific three-dimensional spatial chromatin organization. Here, we map long-range chromosomal conformations in isogenic human dopaminergic, GABAergic, and glutamatergic neurons to track developmentally programmed shifts in the regulatory activity of schizophrenia risk loci. Massive repressive compartmentalization, concomitant with the emergence of hundreds of neuron-specific multi-valent chromatin architectural stripes, occurs during neuronal differentiation, with genes interconnected to genetic risk loci through these long-range chromatin structures differing in their biological roles from genes more proximal to sequences conferring heritable risk. Chemically induced CRISPR-guided chromosomal loop-engineering for the proximal risk gene SNAP91 and distal risk gene BHLHE22 profoundly alters synaptic development and functional activity. Our findings highlight the large-scale cell-type-specific reorganization of chromosomal conformations at schizophrenia risk loci during neurodevelopment and establish a causal link between risk-associated gene-regulatory loop structures and neuronal function.