Topologically associating domains are stable units of replication-timing regulation

Sex-specific DNA-replication in the early mammalian embryo

The timing of DNA replication in mammals is crucial for minimizing errors and influenced by genome usage and chromatin states. Replication timing in the newly formed mammalian embryo remains poorly understood. Here, we have investigated replication timing in mouse zygotes and 2-cell embryos, revealing that zygotes lack a conventional replication timing program, which then emerges in 2-cell embryos. This program differs from embryonic stem cells and generally correlates with transcription and genome compartmentalization of both parental genomes. However, consistent and systematic differences existed between the replication timing of the two parental genomes, including considerably later replication of maternal pericentromeric regions compared to paternal counterparts. Moreover, maternal chromatin modified by Polycomb Repressive Complexes in the oocyte, undergoes early replication, despite belonging to the typically late-replicating B-compartment of the genome. This atypical and asynchronous replication of the two parental genomes may advance our understanding of replication stress in early human embryos and trigger strategies to reduce errors and aneuploidies.

Phase-separated super-enhancers confer an innate radioresistance on genomic DNA

Recently, biomolecular condensates formed through liquid-liquid phase separation have been widely reported to regulate key intracellular processes involved in cell biology and pathogenesis. BRD4 is a nuclear protein instrumental to the establishment of phase-separated super-enhancers (SEs) to direct the transcription of important genes. We previously observed that protein droplets of BRD4 became hydrophobic as their size increase, implying an ability of SEs to limit the ionization of water molecules by irradiation. Here, we aim to establish if SEs confer radiation resistance in cancer cells. We established an in vitro DNA damage assay that measures the effect of radicals provoked by the Fenton reaction on DNA integrity. This revealed that DNA damage was markedly reduced when BRD4 underwent phase separation with DNA. Accordingly, co-focal imaging analyses revealed that SE foci and DNA damage foci are mutually exclusive in irradiated cells. Lastly, we observed that the radioresistance of cancer cells was significantly reduced when irradiation was combined with ARV-771, a BRD4 de-stabilizer. Our data revealed the existence of innately radioresistant genomic regions driven by phase separation in cancer cells. The disruption of these phase-separated components enfolding genomic DNA may represent a novel strategy to augment the effects of radiotherapy.

Replication timing alterations are associated with mutation acquisition during breast and lung cancer evolution

During each cell cycle, the process of DNA replication timing is tightly regulated to ensure the accurate duplication of the genome. The extent and significance of alterations in this process during malignant transformation have not been extensively explored. Here, we assess the impact of altered replication timing (ART) on cancer evolution by analysing replication-timing sequencing of cancer and normal cell lines and 952 whole-genome sequenced lung and breast tumours. We find that 6%-18% of the cancer genome exhibits ART, with regions with a change from early to late replication displaying an increased mutation rate and distinct mutational signatures. Whereas regions changing from late to early replication contain genes with increased expression and present a preponderance of APOBEC3-mediated mutation clusters and associated driver mutations. We demonstrate that ART occurs relatively early during cancer evolution and that ART may have a stronger correlation with mutation acquisition than alterations in chromatin structure.

Adenovirus small E1A directs activation of Alu transcription at YAP/TEAD- and AP-1-bound enhancers through interactions with the EP400 chromatin remodeler

Alu retrotransposons, which form the largest family of mobile DNA elements in the human genome, have recently come to attention as a potential source of regulatory novelties, most notably by participating in enhancer function. Even though Alu transcription by RNA polymerase III is subjected to tight epigenetic silencing, their expression has long been known to increase in response to various types of stress, including viral infection. Here we show that, in primary human fibroblasts, adenovirus small e1a triggered derepression of hundreds of individual Alus by promoting TFIIIB recruitment by Alu-bound TFIIIC. Epigenome profiling revealed an e1a-induced decrease of H3K27 acetylation and increase of H3K4 monomethylation at derepressed Alus, making them resemble poised enhancers. The enhancer nature of e1a-targeted Alus was confirmed by the enrichment, in their upstream regions, of the EP300/CBP acetyltransferase, EP400 chromatin remodeler and YAP1 and FOS transcription factors. The physical interaction of e1a with EP400 was critical for Alu derepression, which was abrogated upon EP400 ablation. Our data suggest that e1a targets a subset of enhancer Alus whose transcriptional activation, which requires EP400 and is mediated by the e1a-EP400 interaction, may participate in the manipulation of enhancer activity by adenoviruses.

PHF2 regulates genome topology and DNA replication in neural stem cells via cohesin

Cohesin plays a crucial role in the organization of topologically-associated domains (TADs), which influence gene expression and DNA replication timing. Whether epigenetic regulators may affect TADs via cohesin to mediate DNA replication remains elusive. Here, we discover that the histone demethylase PHF2 associates with RAD21, a core subunit of cohesin, to regulate DNA replication in mouse neural stem cells (NSC). PHF2 loss impairs DNA replication due to the activation of dormant replication origins in NSC. Notably, the PHF2/RAD21 co-bound genomic regions are characterized by CTCF enrichment and epigenomic features that resemble efficient, active replication origins, and can act as boundaries to separate adjacent domains. Accordingly, PHF2 loss weakens TADs and chromatin loops at the co-bound loci due to reduced RAD21 occupancy. The observed topological and DNA replication defects in PHF2 KO NSC support a cohesin-dependent mechanism. Furthermore, we demonstrate that the PHF2/RAD21 complex exerts little effect on gene regulation, and that PHF2's histone-demethylase activity is dispensable for normal DNA replication and proliferation of NSC. We propose that PHF2 may serve as a topological accessory to cohesin for cohesin localization to TADs and chromatin loops, where cohesin represses dormant replication origins directly or indirectly, to sustain DNA replication in NSC.

Identifying topologically associating domains using differential kernels

Chromatin is a polymer complex of DNA and proteins that regulates gene expression. The three-dimensional (3D) structure and organization of chromatin controls DNA transcription and replication. High-throughput chromatin conformation capture techniques generate Hi-C maps that can provide insight into the 3D structure of chromatin. Hi-C maps can be represented as a symmetric matrix [Formula: see text], where each element represents the average contact probability or number of contacts between chromatin loci i and j. Previous studies have detected topologically associating domains (TADs), or self-interacting regions in [Formula: see text] within which the contact probability is greater than that outside the region. Many algorithms have been developed to identify TADs within Hi-C maps. However, most TAD identification algorithms are unable to identify nested or overlapping TADs and for a given Hi-C map there is significant variation in the location and number of TADs identified by different methods. We develop a novel method to identify TADs, KerTAD, using a kernel-based technique from computer vision and image processing that is able to accurately identify nested and overlapping TADs. We benchmark this method against state-of-the-art TAD identification methods on both synthetic and experimental data sets. We find that the new method consistently has higher true positive rates (TPR) and lower false discovery rates (FDR) than all tested methods for both synthetic and manually annotated experimental Hi-C maps. The TPR for KerTAD is also largely insensitive to increasing noise and sparsity, in contrast to the other methods. We also find that KerTAD is consistent in the number and size of TADs identified across replicate experimental Hi-C maps for several organisms. Thus, KerTAD will improve automated TAD identification and enable researchers to better correlate changes in TADs to biological phenomena, such as enhancer-promoter interactions and disease states.

Quantifying DNA replication speeds in single cells by scEdU-seq

In a human cell, thousands of replication forks simultaneously coordinate duplication of the entire genome. The rate at which this process occurs might depend on the epigenetic state of the genome and vary between, or even within, cell types. To accurately measure DNA replication speeds, we developed single-cell 5-ethynyl-2'-deoxyuridine sequencing to detect nascent replicated DNA. We observed that the DNA replication speed is not constant but increases during S phase of the cell cycle. Using genetic and pharmacological perturbations we were able to alter this acceleration of replication and conclude that DNA damage inflicted by the process of transcription limits the speed of replication during early S phase. In late S phase, during which less-transcribed regions replicate, replication accelerates and approaches its maximum speed.

Differential 3D genome architecture and imprinted gene expression: cause or consequence?

Imprinted genes provide an attractive paradigm to unravel links between transcription and genome architecture. The parental allele-specific expression of these essential genes - which are clustered in chromosomal domains - is mediated by parental methylation imprints at key regulatory DNA sequences. Recent chromatin conformation capture (3C)-based studies show differential organization of topologically associating domains between the parental chromosomes at imprinted domains, in embryonic stem and differentiated cells. At several imprinted domains, differentially methylated regions show allelic binding of the insulator protein CTCF, and linked focal retention of cohesin, at the non-methylated allele only. This generates differential patterns of chromatin looping between the parental chromosomes, already in the early embryo, and thereby facilitates the allelic gene expression. Recent research evokes also the opposite scenario, in which allelic transcription contributes to the differential genome organization, similarly as reported for imprinted X chromosome inactivation. This may occur through epigenetic effects on CTCF binding, through structural effects of RNA Polymerase II, or through imprinted long non-coding RNAs that have chromatin repressive functions. The emerging picture is that epigenetically-controlled differential genome architecture precedes and facilitates imprinted gene expression during development, and that at some domains, conversely, the mono-allelic gene expression also influences genome architecture.

Replication initiation sites and zones in the mammalian genome: Where are they located and how are they defined?

Eukaryotic DNA replication is a tightly controlled process that occurs in two main steps, i.e., licensing and firing, which take place in the G1 and S phases of the cell cycle, respectively. In Saccharomyces cerevisiae, the budding yeast, replication origins contain consensus sequences that are recognized and bound by the licensing factor Orc1-6, which then recruits the replicative Mcm2-7 helicase. By contrast, mammalian initiation sites lack such consensus sequences, and the mammalian ORC does not exhibit sequence specificity. Studies performed over the past decades have identified replication initiation sites in the mammalian genome using sequencing-based assays, raising the question of whether replication initiation occurs at confined sites or in broad zones across the genome. Although recent reports have shown that the licensed MCMs in mammalian cells are broadly distributed, suggesting that ORC-dependent licensing may not determine the initiation sites/zones, they are predominantly located upstream of actively transcribed genes. This review compares the mechanism of replication initiation in yeast and mammalian cells, summarizes the sequencing-based technologies used for the identification of initiation sites/zones, and proposes a possible mechanism of initiation-site/zone selection in mammalian cells. Future directions and challenges in this field are also discussed.

Human DNA topoisomerase I poisoning causes R loop-mediated genome instability attenuated by transcription factor IIS

DNA topoisomerase I can contribute to cancer genome instability. During catalytic activity, topoisomerase I forms a transient intermediate, topoisomerase I-DNA cleavage complex (Top1cc) to allow strand rotation and duplex relaxation, which can lead to elevated levels of DNA-RNA hybrids and micronuclei. To comprehend the underlying mechanisms, we have integrated genomic data of Top1cc-triggered hybrids and DNA double-strand breaks (DSBs) shortly after Top1cc induction, revealing that Top1ccs increase hybrid levels with different mechanisms. DSBs are at highly transcribed genes in early replicating initiation zones and overlap with hybrids downstream of accumulated RNA polymerase II (RNAPII) at gene 5'-ends. A transcription factor IIS mutant impairing transcription elongation further increased RNAPII accumulation likely due to backtracking. Moreover, Top1ccs can trigger micronuclei when occurring during late G1 or early/mid S, but not during late S. As micronuclei and transcription-replication conflicts are attenuated by transcription factor IIS, our results support a role of RNAPII arrest in Top1cc-induced transcription-replication conflicts leading to DSBs and micronuclei.

RACGAP1 knockdown synergizes and enhances the effects of chemotherapeutics on ovarian cancer

Among the three most prevalent cancers affecting the female reproductive system, ovarian cancer (OV) ranks as the second most frequently diagnosed. It is important to investigate the genomic complexity of OV to develop diagnostic and therapeutic strategies. Through the utilization of bioinformatics analysis, it was determined that RacGTPase Activating Protein 1 (RACGAP1) holds significant significance in the field of OV chemotherapeutics, an aspect that has not been thoroughly explored in prior investigations. In our study, a notable increase in RACGAP1 expression was detected in ovarian cancer, demonstrating a robust association with clinicopathological features and patient prognosis. In vivo and in vitro testing revealed that RACGAP1 acts synergistically with chemotherapeutics to enhance their effects on ovarian cancer. Furthermore, an interaction between RACGAP1 and the subunit G2 of the condensin II complex, known as non-SMC condensin II complex subunit G2 (NCAPG2), has been identified. Our findings may provide new insight for improving therapeutic strategies for OV.

Mutation rate heterogeneity at the sub-gene scale due to local DNA hypomethylation

Local mutation rates in human are highly heterogeneous, with known variability at the scale of megabase-sized chromosomal domains, and, on the other extreme, at the scale of oligonucleotides. The intermediate, kilobase-scale heterogeneity in mutation risk is less well characterized. Here, by analyzing thousands of somatic genomes, we studied mutation risk gradients along gene bodies, representing a genomic scale spanning roughly 1-10 kb, hypothesizing that different mutational mechanisms are differently distributed across gene segments. The main heterogeneity concerns several kilobases at the transcription start site and further downstream into 5' ends of gene bodies; these are commonly hypomutated with several mutational signatures, most prominently the ubiquitous C > T changes at CpG dinucleotides. The width and shape of this mutational coldspot at 5' gene ends is variable across genes, and corresponds to variable interval of lowered DNA methylation depending on gene activity level and regulation. Such hypomutated loci, at 5' gene ends or elsewhere, correspond to DNA hypomethylation that can associate with various landmarks, including intragenic enhancers, Polycomb-marked regions, or chromatin loop anchor points. Tissue-specific DNA hypomethylation begets tissue-specific local hypomutation. Of note, direction of mutation risk is inverted for AID/APOBEC3 cytosine deaminase activity, whose signatures are enriched in hypomethylated regions.

Identification of key eRNAs for intervertebral disc degeneration by integrated multinomial bioinformatics analysis

BACKGROUND

Intervertebral disc degeneration (IVDD) is a common degenerative condition leading to abnormal stress distribution under load, causing intervertebral stenosis, facet joint degeneration, and foraminal stenosis. Very little is known about the molecular mechanism of eRNAs in IVDD.

METHODS

Gene expression profiles of 38 annulus disc samples composed of 27 less degenerated discs (LDs) and 11 more degenerated discs (MDs) were retrieved from the GEO database. Then, differentially expressed enhancer RNAs (DEeRNAs), differentially expressed target genes (DETGs), and differentially expressed transcription factors (DETFs), hallmark of cancer signalling pathways according to GSVA; the types and quantity of immune cells according to CIBERSORT; and immune gene sets according to ssGSEA were analysed to construct an IVDD-related eRNA network. Then, multidimensional validation was performed to explore the interactions among DEeRNAs, DETFs and DEGs in space.

RESULTS

A total of 53 components, 14 DETGs, 15 DEeRNAs, 3 DETFs, 5 immune cells, 9 hallmarks, and 7 immune gene sets, were selected to construct the regulatory network. After validation by online multidimensional databases, 21 interactive DEeRNA-DEG-DETF axes related to IVDD exacerbation were identified, among which the C1S-CTNNB1-CHD4 axis was the most significant.

CONCLUSION

Based upon the results of our study, we theorize that the C1S-CTNNB1-CHD4 axis plays a vital role in IVDD exacerbation. Specifically, C1S recruits CTNNB1 and upregulates the expression of CHD4 in IVDD, and subsequently, CHD4 suppresses glycolysis and activates oxidative phosphorylation, thus generating insoluble collagen fibre deposits and leading to the progression of IVDD. Overall, these DEeRNAs could comprise promising therapeutic targets for IVDD due to their high tissue specificity.

Linear interaction between replication and transcription shapes DNA break dynamics at recurrent DNA break Clusters

Recurrent DNA break clusters (RDCs) are replication-transcription collision hotspots; many are unique to neural progenitor cells. Through high-resolution replication sequencing and a capture-ligation assay in mouse neural progenitor cells experiencing replication stress, we unravel the replication features dictating RDC location and orientation. Most RDCs occur at the replication forks traversing timing transition regions (TTRs), where sparse replication origins connect unidirectional forks. Leftward-moving forks generate telomere-connected DNA double-strand breaks (DSBs), while rightward-moving forks lead to centromere-connected DSBs. Strand-specific mapping for DNA-bound RNA reveals co-transcriptional dual-strand DNA:RNA hybrids present at a higher density in RDC than in other actively transcribed long genes. In addition, mapping RNA polymerase activity uncovers that head-to-head interactions between replication and transcription machinery result in 60% DSB contribution to the head-on compared to 40% for co-directional. Taken together we reveal TTR as a fragile class and show how the linear interaction between transcription and replication impacts genome stability.

DNA fragility at topologically associated domain boundaries is promoted by alternative DNA secondary structure and topoisomerase II activity

CCCTC-binding factor (CTCF) binding sites are hotspots of genome instability. Although many factors have been associated with CTCF binding site fragility, no study has integrated all fragility-related factors to understand the mechanism(s) of how they work together. Using an unbiased, genome-wide approach, we found that DNA double-strand breaks (DSBs) are enriched at strong, but not weak, CTCF binding sites in five human cell types. Energetically favorable alternative DNA secondary structures underlie strong CTCF binding sites. These structures coincided with the location of topoisomerase II (TOP2) cleavage complex, suggesting that DNA secondary structure acts as a recognition sequence for TOP2 binding and cleavage at CTCF binding sites. Furthermore, CTCF knockdown significantly increased DSBs at strong CTCF binding sites and at CTCF sites that are located at topologically associated domain (TAD) boundaries. TAD boundary-associated CTCF sites that lost CTCF upon knockdown displayed increased DSBs when compared to the gained sites, and those lost sites are overrepresented with G-quadruplexes, suggesting that the structures act as boundary insulators in the absence of CTCF, and contribute to increased DSBs. These results model how alternative DNA secondary structures facilitate recruitment of TOP2 to CTCF binding sites, providing mechanistic insight into DNA fragility at CTCF binding sites.

Inference of genomic landscapes using ordered Hidden Markov Models with emission densities (oHMMed)

BACKGROUND

Genomes are inherently inhomogeneous, with features such as base composition, recombination, gene density, and gene expression varying along chromosomes. Evolutionary, biological, and biomedical analyses aim to quantify this variation, account for it during inference procedures, and ultimately determine the causal processes behind it. Since sequential observations along chromosomes are not independent, it is unsurprising that autocorrelation patterns have been observed e.g., in human base composition. In this article, we develop a class of Hidden Markov Models (HMMs) called oHMMed (ordered HMM with emission densities, the corresponding R package of the same name is available on CRAN): They identify the number of comparably homogeneous regions within autocorrelated observed sequences. These are modelled as discrete hidden states; the observed data points are realisations of continuous probability distributions with state-specific means that enable ordering of these distributions. The observed sequence is labelled according to the hidden states, permitting only neighbouring states that are also neighbours within the ordering of their associated distributions. The parameters that characterise these state-specific distributions are inferred.

RESULTS

We apply our oHMMed algorithms to the proportion of G and C bases (modelled as a mixture of normal distributions) and the number of genes (modelled as a mixture of poisson-gamma distributions) in windows along the human, mouse, and fruit fly genomes. This results in a partitioning of the genomes into regions by statistically distinguishable averages of these features, and in a characterisation of their continuous patterns of variation. In regard to the genomic G and C proportion, this latter result distinguishes oHMMed from segmentation algorithms based in isochore or compositional domain theory. We further use oHMMed to conduct a detailed analysis of variation of chromatin accessibility (ATAC-seq) and epigenetic markers H3K27ac and H3K27me3 (modelled as a mixture of poisson-gamma distributions) along the human chromosome 1 and their correlations.

CONCLUSIONS

Our algorithms provide a biologically assumption free approach to characterising genomic landscapes shaped by continuous, autocorrelated patterns of variation. Despite this, the resulting genome segmentation enables extraction of compositionally distinct regions for further downstream analyses.

The ORFIUS complex regulates ORC2 localization at replication origins

High-grade serous ovarian cancer (HGSC) is a lethal malignancy with elevated replication stress (RS) levels and defective RS and RS-associated DNA damage responses. Here we demonstrate that the bromodomain-containing protein BRD1 is a RS suppressing protein that forms a replication origin regulatory complex with the histone acetyltransferase HBO1, the BRCA1 tumor suppressor, and BARD1, ORigin FIring Under Stress (ORFIUS). BRD1 and HBO1 promote eventual origin firing by supporting localization of the origin licensing protein ORC2 at origins. In the absence of BRD1 and/or HBO1, both origin firing and nuclei with ORC2 foci are reduced. BRCA1 regulates BRD1, HBO1, and ORC2 localization at replication origins. In the absence of BRCA1, both origin firing and nuclei with BRD1, HBO1, and ORC2 foci are increased. In normal and non-HGSC ovarian cancer cells, the ORFIUS complex responds to ATR and CDC7 origin regulatory signaling and disengages from origins during RS. In BRCA1-mutant and sporadic HGSC cells, BRD1, HBO1, and ORC2 remain associated with replication origins, and unresponsive to RS, DNA damage, or origin regulatory kinase inhibition. ORFIUS complex dysregulation may promote HGSC cell survival by allowing for upregulated origin firing and cell cycle progression despite accumulating DNA damage, and may be a RS target.

Implications of the three-dimensional chromatin organization for genome evolution in a fungal plant pathogen

The spatial organization of eukaryotic genomes is linked to their biological functions, although it is not clear how this impacts the overall evolution of a genome. Here, we uncover the three-dimensional (3D) genome organization of the phytopathogen Verticillium dahliae, known to possess distinct genomic regions, designated adaptive genomic regions (AGRs), enriched in transposable elements and genes that mediate host infection. Short-range DNA interactions form clear topologically associating domains (TADs) with gene-rich boundaries that show reduced levels of gene expression and reduced genomic variation. Intriguingly, TADs are less clearly insulated in AGRs than in the core genome. At a global scale, the genome contains bipartite long-range interactions, particularly enriched for AGRs and more generally containing segmental duplications. Notably, the patterns observed for V. dahliae are also present in other Verticillium species. Thus, our analysis links 3D genome organization to evolutionary features conserved throughout the Verticillium genus.

Oncogenic c-Myc induces replication stress by increasing cohesins chromatin occupancy in a CTCF-dependent manner

Oncogene-induced replication stress is a crucial driver of genomic instability and one of the key events contributing to the onset and evolution of cancer. Despite its critical role in cancer, the mechanisms that generate oncogene-induced replication stress remain not fully understood. Here, we report that an oncogenic c-Myc-dependent increase in cohesins on DNA contributes to the induction of replication stress. Accumulation of cohesins on chromatin is not sufficient to cause replication stress, but also requires cohesins to accumulate at specific sites in a CTCF-dependent manner. We propose that the increased accumulation of cohesins at CTCF site interferes with the progression of replication forks, contributing to oncogene-induced replication stress. This is different from, and independent of, previously suggested mechanisms of oncogene-induced replication stress. This, together with the reported protective role of cohesins in preventing replication stress-induced DNA damage, supports a double-edge involvement of cohesins in causing and tolerating oncogene-induced replication stress.

Cell cycle gene alterations associate with a redistribution of mutation risk across chromosomal domains in human cancers

Mutations in human cells exhibit increased burden in heterochromatic, late DNA replication time (RT) chromosomal domains, with variation in mutation rates between tissues mirroring variation in heterochromatin and RT. We observed that regional mutation risk further varies between individual tumors in a manner independent of cell type, identifying three signatures of domain-scale mutagenesis in >4,000 tumor genomes. The major signature reflects remodeling of heterochromatin and of the RT program domains seen across tumors, tissues and cultured cells, and is robustly linked with higher expression of cell proliferation genes. Regional mutagenesis is associated with loss of activity of the tumor-suppressor genes RB1 and TP53, consistent with their roles in cell cycle control, with distinct mutational patterns generated by the two genes. Loss of regional heterogeneity in mutagenesis is associated with deficiencies in various DNA repair pathways. These mutation risk redistribution processes modify the mutation supply towards important genes, diverting the course of somatic evolution.

Topological reorganization and functional alteration of distinct genomic components in gallbladder cancer

Altered three-dimensional architecture of chromatin influences various genomic regulators and subsequent gene expression in human cancer. However, knowledge of the topological rearrangement of genomic hierarchical layers in cancer is largely limited. Here, by taking advantage of in situ Hi-C, RNA-sequencing, and chromatin immunoprecipitation sequencing (ChIP-seq), we investigated structural reorganization and functional changes in chromosomal compartments, topologically associated domains (TADs), and CCCTC binding factor (CTCF)-mediated loops in gallbladder cancer (GBC) tissues and cell lines. We observed that the chromosomal compartment A/B switch was correlated with CTCF binding levels and gene expression changes. Increased inter-TAD interactions with weaker TAD boundaries were identified in cancer cell lines relative to normal controls. Furthermore, the chromatin short loops and cancer unique loops associated with chromatin remodeling and epithelial-mesenchymal transition activation were enriched in cancer compared with their control counterparts. Cancer-specific enhancer-promoter loops, which contain multiple transcription factor binding motifs, acted as a central element to regulate aberrant gene expression. Depletion of individual enhancers in each loop anchor that connects with promoters led to the inhibition of their corresponding gene expressions. Collectively, our data offer the landscape of hierarchical layers of cancer genome and functional alterations that contribute to the development of GBC.

DNA damage and transcription stress

Genome damage and transcription are intimately linked. Tens to hundreds of thousands of DNA lesions arise in each cell each day, many of which can directly or indirectly impede transcription. Conversely, the process of gene expression is itself a source of endogenous DNA lesions as a result of the susceptibility of single-stranded DNA to damage, conflicts with the DNA replication machinery, and engagement by cells of topoisomerases and base excision repair enzymes to regulate the initiation and progression of gene transcription. Although such processes are tightly regulated and normally accurate, on occasion, they can become abortive and leave behind DNA breaks that can drive genome rearrangements, instability, or cell death.

DNA replication and replication stress response in the context of nuclear architecture

The DNA replication process needs to be coordinated with other DNA metabolism transactions and must eventually extend to the full genome, regardless of chromatin status, gene expression, secondary structures and DNA lesions. Completeness and accuracy of DNA replication are crucial to maintain genome integrity, limiting transformation in normal cells and offering targeting opportunities for proliferating cancer cells. DNA replication is thus tightly coordinated with chromatin dynamics and 3D genome architecture, and we are only beginning to understand the underlying molecular mechanisms. While much has recently been discovered on how DNA replication initiation is organised and modulated in different genomic regions and nuclear territories-the so-called "DNA replication program"-we know much less on how the elongation of ongoing replication forks and particularly the response to replication obstacles is affected by the local nuclear organisation. Also, it is still elusive how specific components of nuclear architecture participate in the replication stress response. Here, we review known mechanisms and factors orchestrating replication initiation, and replication fork progression upon stress, focusing on recent evidence linking genome organisation and nuclear architecture with the cellular responses to replication interference, and highlighting open questions and future challenges to explore this exciting new avenue of research.

Emergence of replication timing during early mammalian development

DNA replication enables genetic inheritance across the kingdoms of life. Replication occurs with a defined temporal order known as the replication timing (RT) programme, leading to organization of the genome into early- or late-replicating regions. RT is cell-type specific, is tightly linked to the three-dimensional nuclear organization of the genome1,2 and is considered an epigenetic fingerprint3. In spite of its importance in maintaining the epigenome4, the developmental regulation of RT in mammals in vivo has not been explored. Here, using single-cell Repli-seq5, we generated genome-wide RT maps of mouse embryos from the zygote to the blastocyst stage. Our data show that RT is initially not well defined but becomes defined progressively from the 4-cell stage, coinciding with strengthening of the A and B compartments. We show that transcription contributes to the precision of the RT programme and that the difference in RT between the A and B compartments depends on RNA polymerase II at zygotic genome activation. Our data indicate that the establishment of nuclear organization precedes the acquisition of defined RT features and primes the partitioning of the genome into early- and late-replicating domains. Our work sheds light on the establishment of the epigenome at the beginning of mammalian development and reveals the organizing principles of genome organization.

A DNA methylation haplotype block landscape in human tissues and preimplantation embryos reveals regulatory elements defined by comethylation patterns

DNA methylation and associated regulatory elements play a crucial role in gene expression regulation. Previous studies have focused primarily on the distribution of mean methylation levels. Advances in whole-genome bisulfite sequencing (WGBS) have enabled the characterization of DNA methylation haplotypes (MHAPs), representing CpG sites from the same read fragment on a single chromosome, and the subsequent identification of methylation haplotype blocks (MHBs), in which adjacent CpGs on the same fragment are comethylated. Using our expert-curated WGBS data sets, we report comprehensive landscapes of MHBs in 17 representative normal somatic human tissues and during early human embryonic development. Integrative analysis reveals MHBs as a distinctive type of regulatory element characterized by comethylation patterns rather than mean methylation levels. We show the enrichment of MHBs in open chromatin regions, tissue-specific histone marks, and enhancers, including super-enhancers. Moreover, we find that MHBs tend to localize near tissue-specific genes and show an association with differential gene expression that is independent of mean methylation. Similar findings are observed in the context of human embryonic development, highlighting the dynamic nature of MHBs during early development. Collectively, our comprehensive MHB landscapes provide valuable insights into the tissue specificity and developmental dynamics of DNA methylation.

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.

RIF1 regulates early replication timing in murine B cells

The mammalian DNA replication timing (RT) program is crucial for the proper functioning and integrity of the genome. The best-known mechanism for controlling RT is the suppression of late origins of replication in heterochromatin by RIF1. Here, we report that in antigen-activated, hypermutating murine B lymphocytes, RIF1 binds predominantly to early-replicating active chromatin and promotes early replication, but plays a minor role in regulating replication origin activity, gene expression and genome organization in B cells. Furthermore, we find that RIF1 functions in a complementary and non-epistatic manner with minichromosome maintenance (MCM) proteins to establish early RT signatures genome-wide and, specifically, to ensure the early replication of highly transcribed genes. These findings reveal additional layers of regulation within the B cell RT program, driven by the coordinated activity of RIF1 and MCM proteins.

MLL3/MLL4 enzymatic activity shapes DNA replication timing

Mammalian genomes are replicated in a precise order during S phase, which is cell-type-specific1-3 and correlates with local transcriptional activity2,4-8, chromatin modifications9 and chromatin architecture1,10,11,12. However, the causal relationships between these features and the key regulators of DNA replication timing (RT) are largely unknown. Here, machine learning was applied to quantify chromatin features, including epigenetic marks, histone variants and chromatin architectural factors, best predicting local RT under steady-state and RT changes during early embryonic stem (ES) cell differentiation. About one-third of genome exhibited RT changes during the differentiation. Combined, chromatin features predicted steady-state RT and RT changes with high accuracy. Of these features, histone H3 lysine 4 monomethylation (H3K4me1) catalyzed by MLL3/4 (also known as KMT2C/D) emerged as a top predictor. Loss of Mll3/4 (but not Mll3 alone) or their enzymatic activity resulted in erasure of genome-wide RT dynamics during ES cell differentiation. Sites that normally gain H3K4me1 in a MLL3/4-dependent fashion during the transition failed to transition towards earlier RT, often with transcriptional activation unaffected. Further analysis revealed a requirement for MLL3/4 in promoting DNA replication initiation zones through MCM2 recruitment, providing a direct link for its role in regulating RT. Our results uncover MLL3/4-dependent H3K4me1 as a functional regulator of RT and highlight a causal relationship between the epigenome and RT that is largely uncoupled from transcription. These findings uncover a previously unknown role for MLL3/4-dependent chromatin functions which is likely relevant to the numerous diseases associated with MLL3/4 mutations.

Transcription induces context-dependent remodeling of chromatin architecture during differentiation

Metazoan chromosomes are organized into discrete spatial domains (TADs), believed to contribute to the regulation of transcriptional programs. Despite extensive correlation between domain organization and gene activity, a direct mechanistic link is unclear, with perturbation studies often showing little effect. To follow chromatin architecture changes during development, we used Capture Hi-C to interrogate the domains around key differentially expressed genes during mouse thymocyte maturation, uncovering specific remodeling events. Notably, one TAD boundary was broadened to accommodate RNA polymerase elongation past the border, and subdomains were formed around some activated genes without changes in CTCF binding. The ectopic induction of some genes was sufficient to recapitulate domain formation in embryonic stem cells, providing strong evidence that transcription can directly remodel chromatin structure. These results suggest that transcriptional processes drive complex chromosome folding patterns that can be important in certain genomic contexts.

Where and when to start: Regulating DNA replication origin activity in eukaryotic genomes

In eukaryotic genomes, hundreds to thousands of potential start sites of DNA replication named origins are dispersed across each of the linear chromosomes. During S-phase, only a subset of origins is selected in a stochastic manner to assemble bidirectional replication forks and initiate DNA synthesis. Despite substantial progress in our understanding of this complex process, a comprehensive 'identity code' that defines origins based on specific nucleotide sequences, DNA structural features, the local chromatin environment, or 3D genome architecture is still missing. In this article, we review the genetic and epigenetic features of replication origins in yeast and metazoan chromosomes and highlight recent insights into how this flexibility in origin usage contributes to nuclear organization, cell growth, differentiation, and genome stability.

Three-dimensional nuclear organisation and the DNA replication timing program

In eukaryotic cells, genome duplication is temporally organised according to a program referred to as the replication-timing (RT) program. The RT of individual genomic domains strikingly parallels the three-dimensional architecture of their chromatin contacts and subnuclear distribution. However, it is unclear whether this correspondence is coincidental or whether it indicates a causal and regulatory relationship. In either case, the nature of the molecular mechanisms ensuring this spatio-temporal coordination is still unknown. Here, we review recent evidence that begins to uncover the existence of a shared molecular machinery at the core of the spatio-temporal co-regulation of DNA replication and genome architecture. Finally, we discuss the outstanding, key question of the biological role of their coordination.

Transcription-Replication Conflicts as a Source of Genome Instability

Transcription and replication both require large macromolecular complexes to act on a DNA template, yet these machineries cannot simultaneously act on the same DNA sequence. Conflicts between the replication and transcription machineries (transcription-replication conflicts, or TRCs) are widespread in both prokaryotes and eukaryotes and have the capacity to both cause DNA damage and compromise complete, faithful replication of the genome. This review will highlight recent studies investigating the genomic locations of TRCs and the mechanisms by which they may be prevented, mitigated, or resolved. We address work from both model organisms and mammalian systems but predominantly focus on multicellular eukaryotes owing to the additional complexities inherent in the coordination of replication and transcription in the context of cell type-specific gene expression and higher-order chromatin organization.

The Response of the Replication Apparatus to Leading Template Strand Blocks

Duplication of the genome requires the replication apparatus to overcome a variety of impediments, including covalent DNA adducts, the most challenging of which is on the leading template strand. Replisomes consist of two functional units, a helicase to unwind DNA and polymerases to synthesize it. The helicase is a multi-protein complex that encircles the leading template strand and makes the first contact with a leading strand adduct. The size of the channel in the helicase would appear to preclude transit by large adducts such as DNA: protein complexes (DPC). Here we discuss some of the extensively studied pathways that support replication restart after replisome encounters with leading template strand adducts. We also call attention to recent work that highlights the tolerance of the helicase for adducts ostensibly too large to pass through the central channel.

Improving Hi-C contact matrices using genome graphs

Three-dimensional chromosome structure plays an important role in fundamental genomic functions. Hi-C, a high-throughput, sequencing-based technique, has drastically expanded our comprehension of 3D chromosome structures. The first step of Hi-C analysis pipeline involves mapping sequencing reads from Hi-C to linear reference genomes. However, the linear reference genome does not incorporate genetic variation information, which can lead to incorrect read alignments, especially when analyzing samples with substantial genomic differences from the reference such as cancer samples. Using genome graphs as the reference facilitates more accurate mapping of reads, however, new algorithms are required for inferring linear genomes from Hi-C reads mapped on genome graphs and constructing corresponding Hi-C contact matrices, which is a prerequisite for the subsequent steps of the Hi-C analysis such as identifying topologically associated domains and calling chromatin loops. We introduce the problem of genome sequence inference from Hi-C data mediated by genome graphs. We formalize this problem, show the hardness of solving this problem, and introduce a novel heuristic algorithm specifically tailored to this problem. We provide a theoretical analysis to evaluate the efficacy of our algorithm. Finally, our empirical experiments indicate that the linear genomes inferred from our method lead to the creation of improved Hi-C contact matrices. These enhanced matrices show a reduction in erroneous patterns caused by structural variations and are more effective in accurately capturing the structures of topologically associated domains.

Perturb-tracing enables high-content screening of multiscale 3D genome regulators

Three-dimensional (3D) genome organization becomes altered during development, aging, and disease1-23, but the factors regulating chromatin topology are incompletely understood and currently no technology can efficiently screen for new regulators of multiscale chromatin organization. Here, we developed an image-based high-content screening platform (Perturb-tracing) that combines pooled CRISPR screen, a new cellular barcode readout method (BARC-FISH), and chromatin tracing. We performed a loss-of-function screen in human cells, and visualized alterations to their genome organization from 13,000 imaging target-perturbation combinations, alongside perturbation-paired barcode readout in the same single cells. Using 1.4 million 3D positions along chromosome traces, we discovered tens of new regulators of chromatin folding at different length scales, ranging from chromatin domains and compartments to chromosome territory. A subset of the regulators exhibited 3D genome effects associated with loop-extrusion and A-B compartmentalization mechanisms, while others were largely unrelated to these known 3D genome mechanisms. We found that the ATP-dependent helicase CHD7, the loss of which causes the congenital neural crest syndrome CHARGE24 and a chromatin remodeler previously shown to promote local chromatin openness25-27, counter-intuitively compacts chromatin over long range in different genomic contexts and cell backgrounds including neural crest cells, and globally represses gene expression. The DNA compaction effect of CHD7 is independent of its chromatin remodeling activity and does not require other protein partners. Finally, we identified new regulators of nuclear architectures and found a functional link between chromatin compaction and nuclear shape. Altogether, our method enables scalable, high-content identification of chromatin and nuclear topology regulators that will stimulate new insights into the 3D genome functions, such as global gene and nuclear regulation, in health and disease.

Quantifying the large-scale chromosome structural dynamics during the mitosis-to-G1 phase transition of cell cycle

Cell cycle is known to be regulated by the underlying gene network. Chromosomes, which serve as the scaffold for gene expressions, undergo significant structural reorganizations during mitosis. Understanding the mechanism of the cell cycle from the chromosome structural perspective remains a grand challenge. In this study, we applied an integrated theoretical approach to investigate large-scale chromosome structural dynamics during the mitosis-to-G1 phase transition. We observed that the chromosome structural expansion and adaptation of the structural asphericity do not occur synchronously and attributed this behaviour to the unique unloading sequence of the two types of condensins. Furthermore, we observed that the coherent motions between the chromosomal loci are primarily enhanced within the topologically associating domains (TADs) as cells progress to the G1 phase, suggesting that TADs can be considered as both structural and dynamical units for organizing the three-dimensional chromosome. Our analysis also reveals that the quantified pathways of chromosome structural reorganization during the mitosis-to-G1 phase transition exhibit high stochasticity at the single-cell level and show nonlinear behaviours in changing TADs and contacts formed at the long-range regions. Our findings offer valuable insights into large-scale chromosome structural dynamics after mitosis.

Repli-seq Sample Preparation using Cell Sorting with Cell-Permeant Dyes

Replication timing is significantly correlated with gene expression and chromatin organization, changes dynamically during cell differentiation, and is altered in diseased states. Genome-wide analysis of replication timing is performed in actively replicating cells by Repli-seq. Current methods for Repli-seq require cells to be fixed in large numbers. This is a barrier for sample types that are sensitive to fixation or are in very limited numbers. In this article, we outline optimized methods to process live cells and intact nuclei for Repli-seq. Our protocol enables the processing of a smaller number of cells per sample and reduces processing time and sample loss while obtaining high-quality data. Further, for samples that tend to form clumps and are difficult to dissociate into a single-cell suspension, we also outline methods for isolation, staining, and processing of nuclei for Repli-seq. The Repli-seq data obtained from live cells and intact nuclei are comparable to those obtained from the standard protocols. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol: Live cell isolation and staining Alternate Protocol: Nuclei isolation and staining.

High-speed AFM imaging reveals DNA capture and loop extrusion dynamics by cohesin-NIPBL

3D chromatin organization plays a critical role in regulating gene expression, DNA replication, recombination, and repair. While initially discovered for its role in sister chromatid cohesion, emerging evidence suggests that the cohesin complex (SMC1, SMC3, RAD21, and SA1/SA2), facilitated by NIPBL, mediates topologically associating domains and chromatin loops through DNA loop extrusion. However, information on how conformational changes of cohesin-NIPBL drive its loading onto DNA, initiation, and growth of DNA loops is still lacking. In this study, high-speed atomic force microscopy imaging reveals that cohesin-NIPBL captures DNA through arm extension, assisted by feet (shorter protrusions), and followed by transfer of DNA to its lower compartment (SMC heads, RAD21, SA1, and NIPBL). While binding at the lower compartment, arm extension leads to the capture of a second DNA segment and the initiation of a DNA loop that is independent of ATP hydrolysis. The feet are likely contributed by the C-terminal domains of SA1 and NIPBL and can transiently bind to DNA to facilitate the loading of the cohesin complex onto DNA. Furthermore, high-speed atomic force microscopy imaging reveals distinct forward and reverse DNA loop extrusion steps by cohesin-NIPBL. These results advance our understanding of cohesin by establishing direct experimental evidence for a multistep DNA-binding mechanism mediated by dynamic protein conformational changes.

Choreography of lamina-associated domains: structure meets dynamics

Lamina-associated domains are large regions of heterochromatin positioned at the nuclear periphery. These domains have been implicated in gene repression, especially in the context of development. In mammals, LAD organization is dependent on nuclear lamins, inner nuclear membrane proteins, and chromatin state. In addition, chromatin readers and modifier proteins have been implicated in this organization, potentially serving as molecular tethers that interact with both nuclear envelope proteins and chromatin. More recent studies have focused on teasing apart the rules that govern dynamic LAD organization and how LAD organization, in turn, relates to gene regulation and overall 3D genome organization. This review highlights recent studies in mammalian cells uncovering factors that instruct the choreography of LAD organization, re-organization, and dynamics at the nuclear lamina, including LAD dynamics in interphase and through mitotic exit, when LAD organization is re-established, as well as intra-LAD subdomain variations.

Coinhibition of topoisomerase 1 and BRD4-mediated pause release selectively kills pancreatic cancer via readthrough transcription

Pancreatic carcinoma lacks effective therapeutic strategies resulting in poor prognosis. Transcriptional dysregulation due to alterations in KRAS and MYC affects initiation, development, and survival of this tumor type. Using patient-derived xenografts of KRAS- and MYC-driven pancreatic carcinoma, we show that coinhibition of topoisomerase 1 (TOP1) and bromodomain-containing protein 4 (BRD4) synergistically induces tumor regression by targeting promoter pause release. Comparing the nascent transcriptome with the recruitment of elongation and termination factors, we found that coinhibition of TOP1 and BRD4 disrupts recruitment of transcription termination factors. Thus, RNA polymerases transcribe downstream of genes for hundreds of kilobases leading to readthrough transcription. This occurs during replication, perturbing replisome progression and inducing DNA damage. The synergistic effect of TOP1 + BRD4 inhibition is specific to cancer cells leaving normal cells unaffected, highlighting the tumor's vulnerability to transcriptional defects. This preclinical study provides a mechanistic understanding of the benefit of combining TOP1 and BRD4 inhibitors to treat pancreatic carcinomas addicted to oncogenic drivers of transcription and replication.

Multi-feature clustering of CTCF binding creates robustness for loop extrusion blocking and Topologically Associating Domain boundaries

Topologically Associating Domains (TADs) separate vertebrate genomes into insulated regulatory neighborhoods that focus genome-associated processes. TADs are formed by Cohesin-mediated loop extrusion, with many TAD boundaries consisting of clustered binding sites of the CTCF insulator protein. Here we determine how this clustering of CTCF binding contributes to the blocking of loop extrusion and the insulation between TADs. We identify enrichment of three features of CTCF binding at strong TAD boundaries, consisting of strongly bound and closely spaced CTCF binding peaks, with a further enrichment of DNA-binding motifs within these peaks. Using multi-contact Nano-C analysis in cells with normal and perturbed CTCF binding, we establish that individual CTCF binding sites contribute to the blocking of loop extrusion, but in an incomplete manner. When clustered, individual CTCF binding sites thus create a stepwise insulation between neighboring TADs. Based on these results, we propose a model whereby multiple instances of temporal loop extrusion blocking create strong insulation between TADs.

The location and development of Replicon Cluster Domains in early replicating DNA

Background: It has been known for many years that in metazoan cells, replication origins are organised into clusters where origins within each cluster fire near-synchronously. Despite clusters being a fundamental organising principle of metazoan DNA replication, the genomic location of origin clusters has not been documented. Methods: We synchronised human U2OS by thymidine block and release followed by L-mimosine block and release to create a population of cells progressing into S phase with a high degree of synchrony. At different times after release into S phase, cells were pulsed with EdU; the EdU-labelled DNA was then pulled down, sequenced and mapped onto the human genome. Results: The early replicating DNA showed features at a range of scales. Wavelet analysis showed that the major feature of the early replicating DNA was at a size of 500 kb, consistent with clusters of replication origins. Over the first two hours of S phase, these Replicon Cluster Domains broadened in width, consistent with their being enlarged by the progression of replication forks at their outer boundaries. The total replication signal associated with each Replicon Cluster Domain varied considerably, and this variation was reproducible and conserved over time. We provide evidence that this variability in replication signal was at least in part caused by Replicon Cluster Domains being activated at different times in different cells in the population. We also provide evidence that adjacent clusters had a statistical preference for being activated in sequence across a group, consistent with the 'domino' model of replication focus activation order observed by microscopy. Conclusions: We show that early replicating DNA is organised into Replicon Cluster Domains that behave as expected of replicon clusters observed by DNA fibre analysis. The coordinated activation of different Replicon Cluster Domains can generate the replication timing programme by which the genome is duplicated.

CTCF and R-loops are boundaries of cohesin-mediated DNA looping

Cohesin and CCCTC-binding factor (CTCF) are key regulatory proteins of three-dimensional (3D) genome organization. Cohesin extrudes DNA loops that are anchored by CTCF in a polar orientation. Here, we present direct evidence that CTCF binding polarity controls cohesin-mediated DNA looping. Using single-molecule imaging, we demonstrate that a critical N-terminal motif of CTCF blocks cohesin translocation and DNA looping. The cryo-EM structure of the cohesin-CTCF complex reveals that this CTCF motif ahead of zinc fingers can only reach its binding site on the STAG1 cohesin subunit when the N terminus of CTCF faces cohesin. Remarkably, a C-terminally oriented CTCF accelerates DNA compaction by cohesin. DNA-bound Cas9 and Cas12a ribonucleoproteins are also polar cohesin barriers, indicating that stalling may be intrinsic to cohesin itself. Finally, we show that RNA-DNA hybrids (R-loops) block cohesin-mediated DNA compaction in vitro and are enriched with cohesin subunits in vivo, likely forming TAD boundaries.

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.

Histone H1 facilitates restoration of H3K27me3 during DNA replication by chromatin compaction

During cell renewal, epigenetic information needs to be precisely restored to maintain cell identity and genome integrity following DNA replication. The histone mark H3K27me3 is essential for the formation of facultative heterochromatin and the repression of developmental genes in embryonic stem cells. However, how the restoration of H3K27me3 is precisely achieved following DNA replication is still poorly understood. Here we employ ChOR-seq (Chromatin Occupancy after Replication) to monitor the dynamic re-establishment of H3K27me3 on nascent DNA during DNA replication. We find that the restoration rate of H3K27me3 is highly correlated with dense chromatin states. In addition, we reveal that the linker histone H1 facilitates the rapid post-replication restoration of H3K27me3 on repressed genes and the restoration rate of H3K27me3 on nascent DNA is greatly compromised after partial depletion of H1. Finally, our in vitro biochemical experiments demonstrate that H1 facilitates the propagation of H3K27me3 by PRC2 through compacting chromatin. Collectively, our results indicate that H1-mediated chromatin compaction facilitates the propagation and restoration of H3K27me3 after DNA replication.

Hi-TrAC detects active sub-TADs and reveals internal organizations of super-enhancers

The spatial folding of eukaryotic genome plays a key role in genome function. We report here that our recently developed method, Hi-TrAC, which specializes in detecting chromatin loops among accessible genomic regions, can detect active sub-TADs with a median size of 100 kb, most of which harbor one or two cell specifically expressed genes and regulatory elements such as super-enhancers organized into nested interaction domains. These active sub-TADs are characterized by highly enriched histone mark H3K4me1 and chromatin-binding proteins, including Cohesin complex. Deletion of selected sub-TAD boundaries have different impacts, such as decreased chromatin interaction and gene expression within the sub-TADs or compromised insulation between the sub-TADs, depending on the specific chromatin environment. We show that knocking down core subunit of the Cohesin complex using shRNAs in human cells or decreasing the H3K4me1 modification by deleting the H3K4 methyltransferase Mll4 gene in mouse Th17 cells disrupted the sub-TADs structure. Our data also suggest that super-enhancers exist as an equilibrium globule structure, while inaccessible chromatin regions exist as a fractal globule structure. In summary, Hi-TrAC serves as a highly sensitive and inexpensive approach to study dynamic changes of active sub-TADs, providing more explicit insights into delicate genome structures and functions.

Oncogenic IDH mutations increase heterochromatin-related replication stress without impacting homologous recombination

Oncogenic mutations in isocitrate dehydrogenases 1 and 2 (IDH1/2) produce 2-hydroxyglutarate (2HG), which inhibits dioxygenases that modulate chromatin dynamics. The effects of 2HG have been reported to sensitize IDH tumors to poly-(ADP-ribose) polymerase (PARP) inhibitors. However, unlike PARP-inhibitor-sensitive BRCA1/2 tumors, which exhibit impaired homologous recombination, IDH-mutant tumors have a silent mutational profile and lack signatures associated with impaired homologous recombination. Instead, 2HG-producing IDH mutations lead to a heterochromatin-dependent slowing of DNA replication accompanied by increased replication stress and DNA double-strand breaks. This replicative stress manifests as replication fork slowing, but the breaks are repaired without a significant increase in mutation burden. Faithful resolution of replicative stress in IDH-mutant cells is dependent on poly-(ADP-ribosylation). Treatment with PARP inhibitors increases DNA replication but results in incomplete DNA repair. These findings demonstrate a role for PARP in the replication of heterochromatin and further validate PARP as a therapeutic target in IDH-mutant tumors.

A spatial genome aligner for resolving chromatin architectures from multiplexed DNA FISH

Multiplexed fluorescence in situ hybridization (FISH) is a widely used approach for analyzing three-dimensional genome organization, but it is challenging to derive chromosomal conformations from noisy fluorescence signals, and tracing chromatin is not straightforward. Here we report a spatial genome aligner that parses true chromatin signal from noise by aligning signals to a DNA polymer model. Using genomic distances separating imaged loci, our aligner estimates spatial distances expected to separate loci on a polymer in three-dimensional space. Our aligner then evaluates the physical probability observed signals belonging to these loci are connected, thereby tracing chromatin structures. We demonstrate that this spatial genome aligner can efficiently model chromosome architectures from DNA FISH data across multiple scales and be used to predict chromosome ploidies de novo in interphase cells. Reprocessing of previous whole-genome chromosome tracing data with this method indicates the spatial aggregation of sister chromatids in S/G2 phase cells in asynchronous mouse embryonic stem cells and provides evidence for extranumerary chromosomes that remain tightly paired in postmitotic neurons of the adult mouse cortex.