Control of DNA replication timing in the 3D genome

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.

TULIPs decorate the three-dimensional genome of PFA ependymoma

Posterior fossa group A (PFA) ependymoma is a lethal brain cancer diagnosed in infants and young children. The lack of driver events in the PFA linear genome led us to search its 3D genome for characteristic features. Here, we reconstructed 3D genomes from diverse childhood tumor types and uncovered a global topology in PFA that is highly reminiscent of stem and progenitor cells in a variety of human tissues. A remarkable feature exclusively present in PFA are type B ultra long-range interactions in PFAs (TULIPs), regions separated by great distances along the linear genome that interact with each other in the 3D nuclear space with surprising strength. TULIPs occur in all PFA samples and recur at predictable genomic coordinates, and their formation is induced by expression of EZHIP. The universality of TULIPs across PFA samples suggests a conservation of molecular principles that could be exploited therapeutically.

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.

DNA replication in early mammalian embryos is patterned, predisposing lamina-associated regions to fragility

DNA replication in differentiated cells follows a defined program, but when and how it is established during mammalian development is not known. Here we show using single-cell sequencing, that late replicating regions are established in association with the B compartment and the nuclear lamina from the first cell cycle after fertilization on both maternal and paternal genomes. Late replicating regions contain a relative paucity of active origins and few but long genes and low G/C content. In both bovine and mouse embryos, replication timing patterns are established prior to embryonic genome activation. Chromosome breaks, which form spontaneously in bovine embryos at sites concordant with human embryos, preferentially locate to late replicating regions. In mice, late replicating regions show enhanced fragility due to a sparsity of dormant origins that can be activated under conditions of replication stress. This pattern predisposes regions with long neuronal genes to fragility and genetic change prior to separation of soma and germ cell lineages. Our studies show that the segregation of early and late replicating regions is among the first layers of genome organization established after fertilization.

Prediction of cell-type-specific cohesin-mediated chromatin loops based on chromatin state

Chromatin loop is of crucial importance for the regulation of gene transcription. Cohesin is a type of chromatin-associated protein that mediates the interaction of chromatin through the loop extrusion. Cohesin-mediated chromatin interactions have strong cell-type specificity, posing a challenge for predicting chromatin loops. Existing computational methods perform poorly in predicting cell-type-specific chromatin loops. To address this issue, we propose a random forest model to predict cell-type-specific cohesin-mediated chromatin loops based on chromatin states identified by ChromHMM and the occupancy of related factors. Our results show that chromatin state is responsible for cell-type-specificity of loops. Using only chromatin states as features, the model achieved high accuracy in predicting cell-type-specific loops between two cell types and can be applied to different cell types. Furthermore, when chromatin states are combined with the occurrence frequency of CTCF, RAD21, YY1, and H3K27ac ChIP-seq peaks, more accurate prediction can be achieved. Our feature extraction method provides novel insights into predicting cell-type-specific chromatin loops and reveals the relationship between chromatin state and chromatin loop formation.

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

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

scGHOST: identifying single-cell 3D genome subcompartments

Single-cell Hi-C (scHi-C) technologies allow for probing of genome-wide cell-to-cell variability in three-dimensional (3D) genome organization from individual cells. Computational methods have been developed to reveal single-cell 3D genome features based on scHi-C, including A/B compartments, topologically associating domains and chromatin loops. However, no method exists for annotating single-cell subcompartments, which is important for understanding chromosome spatial localization in single cells. Here we present scGHOST, a single-cell subcompartment annotation method using graph embedding with constrained random walk sampling. Applications of scGHOST to scHi-C data and contact maps derived from single-cell 3D genome imaging demonstrate reliable identification of single-cell subcompartments, offering insights into cell-to-cell variability of nuclear subcompartments. Using scHi-C data from complex tissues, scGHOST identifies cell-type-specific or allele-specific subcompartments linked to gene transcription across various cell types and developmental stages, suggesting functional implications of single-cell subcompartments. scGHOST is an effective method for annotating single-cell 3D genome subcompartments in a broad range of biological contexts.

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.

A Decade of Discovery-Eukaryotic Replisome Disassembly at Replication Termination

The eukaryotic replicative helicase (CMG complex) is assembled during DNA replication initiation in a highly regulated manner, which is described in depth by other manuscripts in this Issue. During DNA replication, the replicative helicase moves through the chromatin, unwinding DNA and facilitating nascent DNA synthesis by polymerases. Once the duplication of a replicon is complete, the CMG helicase and the remaining components of the replisome need to be removed from the chromatin. Research carried out over the last ten years has produced a breakthrough in our understanding, revealing that replication termination, and more specifically replisome disassembly, is indeed a highly regulated process. This review brings together our current understanding of these processes and highlights elements of the mechanism that are conserved or have undergone divergence throughout evolution. Finally, we discuss events beyond the classic termination of DNA replication in S-phase and go over the known mechanisms of replicative helicase removal from chromatin in these particular situations.

A systematic analyses of different bioinformatics pipelines for genomic data and its impact on deep learning models for chromatin loop prediction

Genomic data analysis has witnessed a surge in complexity and volume, primarily driven by the advent of high-throughput technologies. In particular, studying chromatin loops and structures has become pivotal in understanding gene regulation and genome organization. This systematic investigation explores the realm of specialized bioinformatics pipelines designed specifically for the analysis of chromatin loops and structures. Our investigation incorporates two protein (CTCF and Cohesin) factor-specific loop interaction datasets from six distinct pipelines, amassing a comprehensive collection of 36 diverse datasets. Through a meticulous review of existing literature, we offer a holistic perspective on the methodologies, tools and algorithms underpinning the analysis of this multifaceted genomic feature. We illuminate the vast array of approaches deployed, encompassing pivotal aspects such as data preparation pipeline, preprocessing, statistical features and modelling techniques. Beyond this, we rigorously assess the strengths and limitations inherent in these bioinformatics pipelines, shedding light on the interplay between data quality and the performance of deep learning models, ultimately advancing our comprehension of genomic intricacies.

High resolution maps of chromatin reorganization through mouse meiosis reveal novel features of the 3D meiotic structure

When germ cells transition from the mitotic cycle into meiotic prophase I (MPI), chromosomes condense into an array of chromatin loops that are required to promote homolog pairing and genetic recombination. To identify the changes in chromosomal conformation, we isolated nuclei on a trajectory from spermatogonia to the end of MPI. At each stage along this trajectory, we built genomic interaction maps with the highest temporal and spatial resolution to date. The changes in chromatin folding coincided with a concurrent decline in mitotic cohesion and a rise in meiotic cohesin complexes. We found that the stereotypical large-scale A and B compartmentalization was lost during meiotic prophase I alongside the loss of topological associating domains (TADs). Still, local subcompartments were detected and maintained throughout meiosis. The enhanced Micro-C resolution revealed that, despite the loss of TADs, higher frequency contact sites between two loci were detectable during meiotic prophase I coinciding with CTCF bound sites. The pattern of interactions around these CTCF sites with their neighboring loci showed that CTCF sites were often anchoring the meiotic loops. Additionally, the localization of CTCF to the meiotic axes indicated that these anchors were at the base of loops. Strikingly, even in the face of the dramatic reconfiguration of interphase chromatin into a condensed loop-array, the interactions between regulatory elements remained well preserved. This establishes a potential mechanism for how the meiotic chromatin maintains active transcription within a highly structured genome. In summary, the high temporal and spatial resolution of these data revealed previously unappreciated aspects of mammalian meiotic chromatin organization.

Identifying TAD-like domains on single-cell Hi-C data by graph embedding and changepoint detection

MOTIVATION

Topologically associating domains (TADs) are fundamental building blocks of 3D genome. TAD-like domains in single cells are regarded as the underlying genesis of TADs discovered in bulk cells. Understanding the organization of TAD-like domains helps to get deeper insights into their regulatory functions. Unfortunately, it remains a challenge to identify TAD-like domains on single-cell Hi-C data due to its ultra-sparsity.

RESULTS

We propose scKTLD, an in silico tool for the identification of TAD-like domains on single-cell Hi-C data. It takes Hi-C contact matrix as the adjacency matrix for a graph, embeds the graph structures into a low-dimensional space with the help of sparse matrix factorization followed by spectral propagation, and the TAD-like domains can be identified using a kernel-based changepoint detection in the embedding space. The results tell that our scKTLD is superior to the other methods on the sparse contact matrices, including downsampled bulk Hi-C data as well as simulated and experimental single-cell Hi-C data. Besides, we demonstrated the conservation of TAD-like domain boundaries at single-cell level apart from heterogeneity within and across cell types, and found that the boundaries with higher frequency across single cells are more enriched for architectural proteins and chromatin marks, and they preferentially occur at TAD boundaries in bulk cells, especially at those with higher hierarchical levels.

AVAILABILITY AND IMPLEMENTATION

scKTLD is freely available at https://github.com/lhqxinghun/scKTLD.

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 unraveled 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 revealed 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 revealed that head-to-head interactions between replication and transcription machinery resulted in 60% DSB contribution to the head-on compared to 40% for co-directional. Our findings revealed TTR as a novel fragile class and highlighted how the linear interaction between transcription and replication impacts genome stability.

Quantity and quality of minichromosome maintenance protein complexes couple replication licensing to genome integrity

Accurate and complete replication of genetic information is a fundamental process of every cell division. The replication licensing is the first essential step that lays the foundation for error-free genome duplication. During licensing, minichromosome maintenance protein complexes, the molecular motors of DNA replication, are loaded to genomic sites called replication origins. The correct quantity and functioning of licensed origins are necessary to prevent genome instability associated with severe diseases, including cancer. Here, we delve into recent discoveries that shed light on the novel functions of licensed origins, the pathways necessary for their proper maintenance, and their implications for cancer therapies.

Computational methods for analysing multiscale 3D genome organization

Recent progress in whole-genome mapping and imaging technologies has enabled the characterization of the spatial organization and folding of the genome in the nucleus. In parallel, advanced computational methods have been developed to leverage these mapping data to reveal multiscale three-dimensional (3D) genome features and to provide a more complete view of genome structure and its connections to genome functions such as transcription. Here, we discuss how recently developed computational tools, including machine-learning-based methods and integrative structure-modelling frameworks, have led to a systematic, multiscale delineation of the connections among different scales of 3D genome organization, genomic and epigenomic features, functional nuclear components and genome function. However, approaches that more comprehensively integrate a wide variety of genomic and imaging datasets are still needed to uncover the functional role of 3D genome structure in defining cellular phenotypes in health and disease.

scENCORE: leveraging single-cell epigenetic data to predict chromatin conformation using graph embedding

Dynamic compartmentalization of eukaryotic DNA into active and repressed states enables diverse transcriptional programs to arise from a single genetic blueprint, whereas its dysregulation can be strongly linked to a broad spectrum of diseases. While single-cell Hi-C experiments allow for chromosome conformation profiling across many cells, they are still expensive and not widely available for most labs. Here, we propose an alternate approach, scENCORE, to computationally reconstruct chromatin compartments from the more affordable and widely accessible single-cell epigenetic data. First, scENCORE constructs a long-range epigenetic correlation graph to mimic chromatin interaction frequencies, where nodes and edges represent genome bins and their correlations. Then, it learns the node embeddings to cluster genome regions into A/B compartments and aligns different graphs to quantify chromatin conformation changes across conditions. Benchmarking using cell-type-matched Hi-C experiments demonstrates that scENCORE can robustly reconstruct A/B compartments in a cell-type-specific manner. Furthermore, our chromatin confirmation switching studies highlight substantial compartment-switching events that may introduce substantial regulatory and transcriptional changes in psychiatric disease. In summary, scENCORE allows accurate and cost-effective A/B compartment reconstruction to delineate higher-order chromatin structure heterogeneity in complex tissues.

Genome maintenance meets mechanobiology

Genome stability is key for healthy cells in healthy organisms, and deregulated maintenance of genome integrity is a hallmark of aging and of age-associated diseases including cancer and neurodegeneration. To maintain a stable genome, genome surveillance and repair pathways are closely intertwined with cell cycle regulation and with DNA transactions that occur during transcription and DNA replication. Coordination of these processes across different time and length scales involves dynamic changes of chromatin topology, clustering of fragile genomic regions and repair factors into nuclear repair centers, mobilization of the nuclear cytoskeleton, and activation of cell cycle checkpoints. Here, we provide a general overview of cell cycle regulation and of the processes involved in genome duplication in human cells, followed by an introduction to replication stress and to the cellular responses elicited by perturbed DNA synthesis. We discuss fragile genomic regions that experience high levels of replication stress, with a particular focus on telomere fragility caused by replication stress at the ends of linear chromosomes. Using alternative lengthening of telomeres (ALT) in cancer cells and ALT-associated PML bodies (APBs) as examples of replication stress-associated clustered DNA damage, we discuss compartmentalization of DNA repair reactions and the role of protein properties implicated in phase separation. Finally, we highlight emerging connections between DNA repair and mechanobiology and discuss how biomolecular condensates, components of the nuclear cytoskeleton, and interfaces between membrane-bound organelles and membraneless macromolecular condensates may cooperate to coordinate genome maintenance in space and time.

DNA replication in early mammalian embryos is patterned, predisposing lamina-associated regions to fragility

DNA replication in differentiated cells follows a defined program, but when and how it is established during mammalian development is not known. Here we show using single-cell sequencing, that both bovine and mouse cleavage stage embryos progress through S-phase in a defined pattern. Late replicating regions are associated with the nuclear lamina from the first cell cycle after fertilization, and contain few active origins, and few but long genes. Chromosome breaks, which form spontaneously in bovine embryos at sites concordant with human embryos, preferentially locate to late replicating regions. In mice, late replicating regions show enhanced fragility due to a sparsity of dormant origins that can be activated under conditions of replication stress. This pattern predisposes regions with long neuronal genes to fragility and genetic change prior to segregation of soma and germ line. Our studies show that the formation of early and late replicating regions is among the first layers of epigenetic regulation established on the mammalian genome after fertilization.

Insights into common fragile site instability: DNA replication challenges at DNA repeat sequences

Common fragile sites (CFS) are specific genomic regions prone to chromosomal instability under conditions of DNA replication stress. CFSs manifest as breaks, gaps, and constrictions on metaphase chromosomes under mild replication stress. These replication-sensitive CFS regions are preferentially unstable during cancer development, as reflected by their association with copy number variants (CNVs) frequently arise in most tumor types. Over the years, it became clear that a combination of different characteristics underlies the enhanced sensitivity of CFSs to replication stress. As of today, there is a strong evidence that the core fragility regions along CFSs overlap with actively transcribed large genes with delayed replication timing upon replication stress. Recently, the mechanistic basis for CFS instability was further extended to regions which span topologically associated domain (TAD) boundaries, generating a fragility signature composed of replication, transcription and genome organization. The presence of difficult-to-replicate AT-rich repeats was one of the early features suggested to characterize a subgroup of CFSs. These long stretches of AT-dinucleotide have the potential to fold into stable secondary structures which may impede replication fork progression, leaving the region under-replicated. Here, we focus on the molecular mechanisms underlying repeat instability at CFSs and on the proteins involved in the resolution of secondary structure impediments arising along repetitive sequence elements which are essential for the maintenance of genome stability.

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.

FAM111A regulates replication origin activation and cell fitness

FAM111A is a replisome-associated protein and dominant mutations within its trypsin-like peptidase domain are linked to severe human developmental syndrome, the Kenny-Caffey syndrome. However, FAM111A functions remain unclear. Here, we show that FAM111A facilitates efficient activation of DNA replication origins. Upon hydroxyurea treatment, FAM111A-depleted cells exhibit reduced single-stranded DNA formation and a better survival rate. Unrestrained expression of FAM111A WT and patient mutants causes accumulation of DNA damage and cell death, only when the peptidase domain remains intact. Unrestrained expression of FAM111A WT also causes increased single-stranded DNA formation that relies on S phase entry, FAM111A peptidase activity but not its binding to proliferating cell nuclear antigen. Altogether, these data unveil how FAM111A promotes DNA replication under normal conditions and becomes harmful in a disease context.

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.

Subgraph extraction and graph representation learning for single cell Hi-C imputation and clustering

Single-cell Hi-C (scHi-C) technology enables the investigation of 3D chromatin structure variability across individual cells. However, the analysis of scHi-C data is challenged by a large number of missing values. Here, we present a scHi-C data imputation model HiC-SGL, based on Subgraph extraction and graph representation learning. HiC-SGL can also learn informative low-dimensional embeddings of cells. We demonstrate that our method surpasses existing methods in terms of imputation accuracy and clustering performance by various metrics.

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.

Brd2 is dispensable for genome compartmentalization and replication timing

Replication Timing (RT) refers to the temporal order in which the genome is replicated during S phase. Early replicating regions correlate with the transcriptionally active, accessible euchromatin (A) compartment, while late replicating regions correlate with the heterochromatin (B) compartment and repressive histone marks. Previously, widespread A/B genome compartmentalization changes were reported following Brd2 depletion. Since RT and A/B compartmentalization are two of the most highly correlated chromosome properties, we evaluated the effects of Brd2 depletion on RT. We performed E/L Repli-Seq following Brd2 depletion in the previously described Brd2 conditional degron cell line and found no significant alterations in RT after Brd2 KD. This finding prompted us to re-analyze the Micro-C data from the previous publication. We report that we were unable to detect any compartmentalization changes in Brd2 depleted cells compared to DMSO control using the same data. Taken together, our findings demonstrate that Brd2 depletion alone does not affect A/B compartmentalization or RT in mouse embryonic stem cells.

Functional filter for whole-genome sequencing data identifies HHT and stress-associated non-coding SMAD4 polyadenylation site variants >5 kb from coding DNA

Despite whole-genome sequencing (WGS), many cases of single-gene disorders remain unsolved, impeding diagnosis and preventative care for people whose disease-causing variants escape detection. Since early WGS data analytic steps prioritize protein-coding sequences, to simultaneously prioritize variants in non-coding regions rich in transcribed and critical regulatory sequences, we developed GROFFFY, an analytic tool that integrates coordinates for regions with experimental evidence of functionality. Applied to WGS data from solved and unsolved hereditary hemorrhagic telangiectasia (HHT) recruits to the 100,000 Genomes Project, GROFFFY-based filtration reduced the mean number of variants/DNA from 4,867,167 to 21,486, without deleting disease-causal variants. In three unsolved cases (two related), GROFFFY identified ultra-rare deletions within the 3' untranslated region (UTR) of the tumor suppressor SMAD4, where germline loss-of-function alleles cause combined HHT and colonic polyposis (MIM: 175050). Sited >5.4 kb distal to coding DNA, the deletions did not modify or generate microRNA binding sites, but instead disrupted the sequence context of the final cleavage and polyadenylation site necessary for protein production: By iFoldRNA, an AAUAAA-adjacent 16-nucleotide deletion brought the cleavage site into inaccessible neighboring secondary structures, while a 4-nucleotide deletion unfolded the downstream RNA polymerase II roadblock. SMAD4 RNA expression differed to control-derived RNA from resting and cycloheximide-stressed peripheral blood mononuclear cells. Patterns predicted the mutational site for an unrelated HHT/polyposis-affected individual, where a complex insertion was subsequently identified. In conclusion, we describe a functional rare variant type that impacts regulatory systems based on RNA polyadenylation. Extension of coding sequence-focused gene panels is required to capture these variants.

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.

Identifying quantitatively differential chromosomal compartmentalization changes and their biological significance from Hi-C data using DARIC

BACKGROUND

Chromosomal compartmentalization plays a critical role in maintaining proper transcriptional programs in cell differentiation and oncogenesis. However, currently the prevalent method for comparative analysis of compartmentalization landscapes between different cell types is limited to the qualitative switched compartments.

RESULTS

To identify genomic regions with quantitatively differential compartmentalization changes from genome-wide chromatin conformation data like Hi-C, we developed a computational framework named DARIC. DARIC includes three modules: compartmentalization quantification, normalization, and differential analysis. Comparing DARIC with the conventional compartment switching analysis reveals substantial regions characterized by quantitatively significant compartmentalization changes without switching. These changes are accompanied by changes in gene expression, chromatin accessibility, H3K27ac intensity, as well as the interactions with nuclear lamina proteins and nuclear positioning, highlighting the functional importance of such quantitative changes in gene regulation. We applied DARIC to dissect the quantitative compartmentalization changes during human cardiomyocyte differentiation and identified two distinct mechanisms for gene activation based on the association with compartmentalization changes. Using the quantitative compartmentalization measurement module from DARIC, we further dissected the compartment variability landscape in the human genome by analyzing a compendium of 32 Hi-C datasets from 4DN. We discovered an interesting correlation between compartmentalization variability and sub-compartments.

CONCLUSIONS

DARIC is a useful tool for analyzing quantitative compartmentalization changes and mining novel biological insights from increasing Hi-C data. Our results demonstrate the functional significance of quantitative compartmentalization changes in gene regulation, and provide new insights into the relationship between compartmentalization variability and sub-compartments in the human genome.

A CRISPR-based rapid DNA repositioning strategy and the early intranuclear life of HSV-1

The relative positions of viral DNA genomes to the host intranuclear environment play critical roles in determining virus fate. Recent advances in the application of chromosome conformation capture-based sequencing analysis (3 C technologies) have revealed valuable aspects of the spatiotemporal interplay of viral genomes with host chromosomes. However, to elucidate the causal relationship between the subnuclear localization of viral genomes and the pathogenic outcome of an infection, manipulative tools are needed. Rapid repositioning of viral DNAs to specific subnuclear compartments amid infection is a powerful approach to synchronize and interrogate this dynamically changing process in space and time. Herein, we report an inducible CRISPR-based two-component platform that relocates extrachromosomal DNA pieces (5 kb to 170 kb) to the nuclear periphery in minutes (CRISPR-nuPin). Based on this strategy, investigations of herpes simplex virus 1 (HSV-1), a prototypical member of the human herpesvirus family, revealed unprecedently reported insights into the early intranuclear life of the pathogen: (I) Viral genomes tethered to the nuclear periphery upon entry, compared with those freely infecting the nucleus, were wrapped around histones with increased suppressive modifications and subjected to stronger transcriptional silencing and prominent growth inhibition. (II) Relocating HSV-1 genomes at 1 hr post infection significantly promoted the transcription of viral genes, termed an 'Escaping' effect. (III) Early accumulation of ICP0 was a sufficient but not necessary condition for 'Escaping'. (IV) Subnuclear localization was only critical during early infection. Importantly, the CRISPR-nuPin tactic, in principle, is applicable to many other DNA viruses.

G-quadruplex DNA contributes to RNA polymerase II-mediated 3D chromatin architecture

High-order chromatin organization plays an important role in biological processes and disease development. Previous studies revealed a widespread occurrence of guanine quadruplex (G4) structures in the human genome, with enrichment in gene regulatory regions, especially in promoters. However, it remains unclear whether G4 structures contribute to RNA polymerase II (RNAPII)-mediated long-range DNA interactions and transcription activity. In this study, we conducted an intuitive overlapping analysis of previously published RNAPII ChIA-PET (chromatin interaction analysis with paired-end tag) and BG4 ChIP-seq (chromatin immunoprecipitation followed by sequencing using a G4 structure-specific antibody) data. We observed a strong positive correlation between RNAPII-linked DNA loops and G4 structures in chromatin. Additionally, our RNAPII HiChIP-seq (in situ Hi-C followed by ChIP-seq) results showed that treatment of HepG2 cells with pyridostatin (PDS), a small-molecule G4-binding ligand, could diminish RNAPII-linked long-range DNA contacts, with more pronounced diminutions being observed for those contacts involving G4 structure loci. RNA sequencing data revealed that PDS treatment modulates the expression of not only genes with G4 structures in their promoters, but also those with promoters being connected with distal G4s through RNAPII-linked long-range DNA interactions. Together, our data substantiate the function of DNA G4s in RNAPII-associated DNA looping and transcription regulation.

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.

Both cell autonomous and non-autonomous processes modulate the association between replication timing and mutation rate

Cancer somatic mutations are the product of multiple mutational and repair processes, some of which are tightly associated with DNA replication. Mutation rates (MR) are known to be higher in late replication timing (RT) regions, but different processes can affect this association. Systematic analysis of the mutational landscape of 2787 tumors from 32 tumor types revealed that approximately one third of the tumor samples show weak association between replication timing and mutation rate. Further analyses revealed that those samples have unique mutational signatures and are enriched with mutations in genes involved in DNA replication, DNA repair and chromatin structure. Surprisingly, analysis of differentially expressed genes between weak and strong RT-MR association groups revealed that tumors with weak association are enriched with genes associated with cell-cell communication and the immune system, suggesting a non-autonomous response to DNA damage.

Dimeric G-quadruplex motifs-induced NFRs determine strong replication origins in vertebrates

Replication of vertebrate genomes is tightly regulated to ensure accurate duplication, but our understanding of the interplay between genetic and epigenetic factors in this regulation remains incomplete. Here, we investigated the involvement of three elements enriched at gene promoters and replication origins: guanine-rich motifs potentially forming G-quadruplexes (pG4s), nucleosome-free regions (NFRs), and the histone variant H2A.Z, in the firing of origins of replication in vertebrates. We show that two pG4s on the same DNA strand (dimeric pG4s) are sufficient to induce the assembly of an efficient minimal replication origin without inducing transcription in avian DT40 cells. Dimeric pG4s in replication origins are associated with formation of an NFR next to precisely-positioned nucleosomes enriched in H2A.Z on this minimal origin and genome-wide. Thus, our data suggest that dimeric pG4s are important for the organization and duplication of vertebrate genomes. It supports the hypothesis that a nucleosome close to an NFR is a shared signal for the formation of replication origins in eukaryotes.

Global effects of identity and aging on the human sperm methylome

BACKGROUND

As the average age of fatherhood increases worldwide, so too does the need for understanding effects of aging in male germline cells. Molecular change, including epigenomic alterations, may impact offspring. Age-associated change to DNA cytosine methylation in the cytosine-guanine (CpG) context is a hallmark of aging tissues, including sperm. Prior studies have led to accurate models that predict a man's age based on specific methylation features in the DNA of sperm, but the relationship between aging and global DNA methylation in sperm remains opaque. Further clarification requires a more complete survey of the methylome with assessment of variability within and between individuals.

RESULTS

We collected sperm methylome data in a longitudinal study of ten healthy fertile men. We used whole-genome bisulfite sequencing of samples collected 10 to 18 years apart from each donor. We found that, overall, variability between donors far exceeds age-associated variation. After controlling for donor identity, we see significant age-dependent genome-wide change to the methylome. Notably, trends of change with age depend on genomic location or annotation, with contrasting signatures that correlate with gene density and proximity to centromeres and promoter regions.

CONCLUSIONS

We uncovered epigenetic signatures that reflect a stable process which begins in early adulthood, progressing steadily through most of the male lifespan, and warrants consideration in any future study of the aging sperm epigenome.

Rif1 restrains the rate of replication origin firing in Xenopus laevis

Metazoan genomes are duplicated by the coordinated activation of clusters of replication origins at different times during S phase, but the underlying mechanisms of this temporal program remain unclear during early development. Rif1, a key replication timing factor, inhibits origin firing by recruiting protein phosphatase 1 (PP1) to chromatin counteracting S phase kinases. We have previously described that Rif1 depletion accelerates early Xenopus laevis embryonic cell cycles. Here, we find that in the absence of Rif1, patterns of replication foci change along with the acceleration of replication cluster activation. However, initiations increase only moderately inside active clusters. Our numerical simulations suggest that the absence of Rif1 compresses the temporal program towards more homogeneity and increases the availability of limiting initiation factors. We experimentally demonstrate that Rif1 depletion increases the chromatin-binding of the S phase kinase Cdc7/Drf1, the firing factors Treslin, MTBP, Cdc45, RecQL4, and the phosphorylation of both Treslin and MTBP. We show that Rif1 globally, but not locally, restrains the replication program in early embryos, possibly by inhibiting or excluding replication factors from chromatin.

Concurrent profiling of multiscale 3D genome organization and gene expression in single mammalian cells

The organization of mammalian genomes within the nucleus features a complex, multiscale three-dimensional (3D) architecture. The functional significance of these 3D genome features, however, remains largely elusive due to limited single-cell technologies that can concurrently profile genome organization and transcriptional activities. Here, we report GAGE-seq, a highly scalable, robust single-cell co-assay that simultaneously measures 3D genome structure and transcriptome within the same cell. Employing GAGE-seq on mouse brain cortex and human bone marrow CD34+ cells, we comprehensively characterized the intricate relationships between 3D genome and gene expression. We found that these multiscale 3D genome features collectively inform cell type-specific gene expressions, hence contributing to defining cell identity at the single-cell level. Integration of GAGE-seq data with spatial transcriptomic data revealed in situ variations of the 3D genome in mouse cortex. Moreover, our observations of lineage commitment in normal human hematopoiesis unveiled notable discordant changes between 3D genome organization and gene expression, underscoring a complex, temporal interplay at the single-cell level that is more nuanced than previously appreciated. Together, GAGE-seq provides a powerful, cost-effective approach for interrogating genome structure and gene expression relationships at the single-cell level across diverse biological contexts.

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.

Single cell Hi-C identifies plastic chromosome conformations underlying the gastrulation enhancer landscape

Embryonic development involves massive proliferation and differentiation of cell lineages. This must be supported by chromosome replication and epigenetic reprogramming, but how proliferation and cell fate acquisition are balanced in this process is not well understood. Here we use single cell Hi-C to map chromosomal conformations in post-gastrulation mouse embryo cells and study their distributions and correlations with matching embryonic transcriptional atlases. We find that embryonic chromosomes show a remarkably strong cell cycle signature. Despite that, replication timing, chromosome compartment structure, topological associated domains (TADs) and promoter-enhancer contacts are shown to be variable between distinct epigenetic states. About 10% of the nuclei are identified as primitive erythrocytes, showing exceptionally compact and organized compartment structure. The remaining cells are broadly associated with ectoderm and mesoderm identities, showing only mild differentiation of TADs and compartment structures, but more specific localized contacts in hundreds of ectoderm and mesoderm promoter-enhancer pairs. The data suggest that while fully committed embryonic lineages can rapidly acquire specific chromosomal conformations, most embryonic cells are showing plastic signatures driven by complex and intermixed enhancer landscapes.

Flexible Attachment and Detachment of Centromeres and Telomeres to and from Chromosomes

Accurate transmission of genomic information across multiple cell divisions and generations, without any losses or errors, is fundamental to all living organisms. To achieve this goal, eukaryotes devised chromosomes. Eukaryotic genomes are represented by multiple linear chromosomes in the nucleus, each carrying a centromere in the middle, a telomere at both ends, and multiple origins of replication along the chromosome arms. Although all three of these DNA elements are indispensable for chromosome function, centromeres and telomeres possess the potential to detach from the original chromosome and attach to new chromosomal positions, as evident from the events of telomere fusion, centromere inactivation, telomere healing, and neocentromere formation. These events seem to occur spontaneously in nature but have not yet been elucidated clearly, because they are relatively infrequent and sometimes detrimental. To address this issue, experimental setups have been developed using model organisms such as yeast. In this article, we review some of the key experiments that provide clues as to the extent to which these paradoxical and elusive features of chromosomally indispensable elements may become valuable in the natural context.

Determining chromatin architecture with Micro Capture-C

Micro Capture-C (MCC) is a chromatin conformation capture (3C) method for visualizing reproducible three-dimensional contacts of specified regions of the genome at base pair resolution. These methods are an established family of techniques that use proximity ligation to assay the topology of chromatin. MCC can generate data at substantially higher resolution than previous techniques through multiple refinements of the 3C method. Using a sequence agnostic nuclease, the maintenance of cellular integrity and full sequencing of the ligation junctions, MCC achieves subnucleosomal levels of resolution, which can be used to reveal transcription factor binding sites analogous to DNAse I footprinting. Gene dense regions, close-range enhancer-promoter contacts, individual enhancers within super-enhancers and multiple other types of loci or regulatory regions that were previously challenging to assay with conventional 3C techniques, are readily observed using MCC. MCC requires training in common molecular biology techniques and bioinformatics to perform the experiment and analyze the data. The protocol can be expected to be completed in a 3 week timeframe for experienced molecular biologists.

scGHOST: Identifying single-cell 3D genome subcompartments

New single-cell Hi-C (scHi-C) technologies enable probing of the genome-wide cell-to-cell variability in 3D genome organization from individual cells. Several computational methods have been developed to reveal single-cell 3D genome features based on scHi-C data, including A/B compartments, topologically-associating domains, and chromatin loops. However, no scHi-C analysis method currently exists for annotating single-cell subcompartments, which are crucial for providing a more refined view of large-scale chromosome spatial localization in single cells. Here, we present scGhost, a single-cell subcompartment annotation method based on graph embedding with constrained random walk sampling. Applications of scGhost to scHi-C data and single-cell 3D genome imaging data demonstrate the reliable identification of single-cell subcompartments and offer new insights into cell-to-cell variability of nuclear subcompartments. Using scHi-C data from the human prefrontal cortex, scGhost identifies cell type-specific subcompartments that are strongly connected to cell type-specific gene expression, suggesting the functional implications of single-cell subcompartments. Overall, scGhost is an effective new method for single-cell 3D genome subcompartment annotation based on scHi-C data for a broad range of biological contexts.

Replication timing and transcriptional control: beyond cause and effect - part IV

Decades of work on the spatiotemporal organization of mammalian DNA replication timing (RT) continues to unveil novel correlations with aspects of transcription and chromatin organization but, until recently, mechanisms regulating RT and the biological significance of the RT program had been indistinct. We now know that the RT program is both influenced by and necessary to maintain chromatin structure, forming an epigenetic positive feedback loop. Moreover, the discovery of specific cis-acting elements regulating mammalian RT at both the domain and the whole-chromosome level has revealed multiple cell-type-specific and developmentally regulated mechanisms of RT control. We review recent evidence for diverse mechanisms employed by different cell types to regulate their RT programs and the biological significance of RT regulation during development.

Profiling Chromatin Accessibility on Replicated DNA with repli-ATAC-Seq

Open or accessible chromatin typifies euchromatic regions and helps define cell type-specific transcription programs. DNA replication massively disorders chromatin composition and structure, and how accessible regions are affected by and recover from this disruption has been unclear. Here, we present repli-ATAC-seq, a protocol to profile accessible chromatin genome-wide on replicated DNA starting from 100,000 cells. In this method, replicated DNA is labeled with a short 5-ethynyl-2'-deoxyuridine (EdU) pulse in cultured cells and isolated from a population of tagmented fragments for amplification and next-generation sequencing. Repli-ATAC-seq provides high-resolution information on chromatin dynamics after DNA replication and reveals new insights into the interplay between DNA replication, transcription, and the chromatin landscape.

Deep generative modeling and clustering of single cell Hi-C data

Deciphering 3D genome conformation is important for understanding gene regulation and cellular function at a spatial level. The recent advances of single cell Hi-C technologies have enabled the profiling of the 3D architecture of DNA within individual cell, which allows us to study the cell-to-cell variability of 3D chromatin organization. Computational approaches are in urgent need to comprehensively analyze the sparse and heterogeneous single cell Hi-C data. Here, we proposed scDEC-Hi-C, a new framework for single cell Hi-C analysis with deep generative neural networks. scDEC-Hi-C outperforms existing methods in terms of single cell Hi-C data clustering and imputation. Moreover, the generative power of scDEC-Hi-C could help unveil the differences of chromatin architecture across cell types. We expect that scDEC-Hi-C could shed light on deepening our understanding of the complex mechanism underlying the formation of chromatin contacts.

Double-strand break repair and mis-repair in 3D

DNA double-strand breaks (DSBs) are lesions that arise frequently from exposure to damaging agents as well as from ongoing physiological DNA transactions. Mis-repair of DSBs leads to rearrangements and structural variations in chromosomes, including insertions, deletions, and translocations implicated in disease. The DNA damage response (DDR) limits pathologic mutations and large-scale chromosome rearrangements. DSB repair initiates in 2D at DNA lesions with the stepwise recruitment of repair proteins and local chromatin remodeling which facilitates break accessibility. More complex structures are then formed via protein assembly into nanodomains and via genome folding into chromatin loops. Subsequently, 3D reorganization of DSBs is guided by clustering forces which drive the assembly of repair domains harboring multiple lesions. These domains are further stabilized and insulated into condensates via liquid-liquid phase-separation. Here, we discuss the benefits and risks associated with this 3D reorganization of the broken genome.

Spatiotemporal Epigenetic Control of the Histone Gene Chromatin Landscape during the Cell Cycle

Higher-order genomic organization supports the activation of histone genes in response to cell cycle regulatory cues that epigenetically mediates stringent control of transcription at the G1/S-phase transition. Histone locus bodies (HLBs) are dynamic, non-membranous, phase-separated nuclear domains where the regulatory machinery for histone gene expression is organized and assembled to support spatiotemporal epigenetic control of histone genes. HLBs provide molecular hubs that support synthesis and processing of DNA replication-dependent histone mRNAs. These regulatory microenvironments support long-range genomic interactions among non-contiguous histone genes within a single topologically associating domain (TAD). HLBs respond to activation of the cyclin E/CDK2/NPAT/HINFP pathway at the G1/S transition. HINFP and its coactivator NPAT form a complex within HLBs that controls histone mRNA transcription to support histone protein synthesis and packaging of newly replicated DNA. Loss of HINFP compromises H4 gene expression and chromatin formation, which may result in DNA damage and impede cell cycle progression. HLBs provide a paradigm for higher-order genomic organization of a subnuclear domain that executes an obligatory cell cycle-controlled function in response to cyclin E/CDK2 signaling. Understanding the coordinately and spatiotemporally organized regulatory programs in focally defined nuclear domains provides insight into molecular infrastructure for responsiveness to cell signaling pathways that mediate biological control of growth, differentiation phenotype, and are compromised in cancer.