Replication timing maintains the global epigenetic state in human cells

RIF1 controls replication timing in early mouse embryos independently of lamina-associated nuclear organization

Cells must duplicate their genome before they divide to ensure equal transmission of genetic information. The genome is replicated with a defined temporal order, replication timing (RT), which is cell-type specific and linked to 3D-genome organization. During mammalian development, RT is initially not well defined and becomes progressively consolidated from the 4-cell stage. However, the molecular regulators are unknown. Here, by combining loss-of-function analysis with genome-wide investigation of RT in mouse embryos, we identify Rap1 interacting factor 1 (RIF1) as a regulator of the progressive consolidation of RT. Embryos without RIF1 show DNA replication features of an early, more totipotent state. RIF1 regulates the progressive stratification of RT values and its depletion leads to global RT changes and a more heterogeneous RT program. Developmental RT changes are disentangled from changes in transcription and nuclear organization, specifically nuclear lamina association. Our work provides molecular understanding of replication and genome organization at the beginning of mammalian development.

Epigenome dynamics in early mammalian embryogenesis

During early embryonic development in mammals, the totipotency of the zygote - which is reprogrammed from the differentiated gametes - transitions to pluripotency by the blastocyst stage, coincident with the first cell fate decision. These changes in cellular potency are accompanied by large-scale alterations in the nucleus, including major transcriptional, epigenetic and architectural remodelling, and the establishment of the DNA replication programme. Advances in low-input genomics and loss-of-function methodologies tailored to the pre-implantation embryo now enable these processes to be studied at an unprecedented level of molecular detail in vivo. Such studies have provided new insights into the genome-wide landscape of epigenetic reprogramming and chromatin dynamics during this fundamental period of pre-implantation development.

Integrated promoter-capture Hi-C and Hi-C analysis reveals fine-tuned regulation of the 3D chromatin architecture in colorectal cancer

Introduction

Hi-C is a widely used technique for mapping chromosomal interactions within a 3D genomic framework, however, its resolution is often constrained by sequencing depth, making it challenging to detect fine-scale interactions. To overcome this limitation, Promoter-Capture Hi-C (PCHi-C), as it selectively enriches for promoter-associated interactions, was employed. This study integrates PCHi-C and Hi-C datasets from colorectal cancer (CRC) models investigate chromosomal interaction dynamics across various regulatory levels, from cis-regulatory elements to topologically associated domains (TADs). The primary goal is to examine how genomic structural alterations shape the epigenomic landscape in CRC and to assess their potential role in colorectal cancer susceptibility.

Methods

PCHi-C and Hi-C datasets from multiple colorectal cancer (CRC) studies were integrated to enhance the resolution of chromatin interaction mapping. The analysis focused on identifying fine-scale interactions within topologically associated domains (TADs) while incorporating histone modification landscapes (H3K27ac, H3K4me3) and transcriptomic signatures from CRC cell lines and the TCGA database. For experimental validation, ChIP-quantitative PCR was performed at the promoters of target genes using the highly malignant colorectal cell line HT29 and compared it to an embryonic cell line NT2D1.

Results

Our integrated analysis revealed significant genomic structural instability in CRC cells, closely associated with tumor-suppressive transcriptional programs. We identified nine dysregulated genes, including long non-coding RNAs (MALAT1, NEAT1, FTX, and PVT1), small nucleolar RNAs (SNORA26 and SNORA71A), and protein-coding genes (TMPRSS11D, TSPEAR, and DSG4), all of which exhibited a substantial increase in expression in CRC cell lines compared to human embryonic stem cells (hESCs). Additionally, we observed enriched activation-associated histone modifications (H3K27ac and H3K4me3) at the potential enhancer regions of these genes, indicating possible transcriptional activation. ChIP-quantitative PCRs conducted using in the highly malignant CRC cell line HT29, compared to the embryonic cell line NT2D1, further validated these findings, reinforcing the link between altered chromosomal interactions and gene dysregulation in CRC.

Discussion

This study sheds light on the dynamic 3D genome organization in CRC, highlighting critical structural changes associated with disease-associated loci. The identification of nine dysregulated genes points to potential biomarkers for colorectal cancer, with implications for diagnostic and therapeutic strategies. The combination of Hi-C and PCHi-C offers a refined approach for detecting chromosomal interactions at a higher resolution, laying the foundation for future studies on cancer-associated chromatin architecture.

Replication-dependent histone labeling dissects the physical properties of euchromatin/heterochromatin in living human cells

A string of nucleosomes, where genomic DNA is wrapped around histones, is organized in the cell as chromatin, ranging from euchromatin to heterochromatin, with distinct genome functions. Understanding physical differences between euchromatin and heterochromatin is crucial, yet specific labeling methods in living cells remain limited. Here, we have developed replication-dependent histone (Repli-Histo) labeling to mark nucleosomes in euchromatin and heterochromatin based on DNA replication timing. Using this approach, we investigated local nucleosome motion in the four known chromatin classes, from euchromatin to heterochromatin, of living human and mouse cells. The more euchromatic (earlier-replicated) and more heterochromatic (later-replicated) regions exhibit greater and lesser nucleosome motions, respectively. Notably, the motion profile in each chromatin class persists throughout interphase. Genome chromatin is essentially replicated from regions with greater nucleosome motions, although the replication timing is perturbed. Our findings, combined with computational modeling, suggest that earlier-replicated regions have more accessibility, and local chromatin motion can be a major determinant of genome-wide replication timing.

Dynamics of replication timing during mammalian development

Recent developments in low-input genomics techniques have greatly advanced the analysis of the order in which DNA is replicated in the genome - that is, replication timing (RT) - and its interrelationships with other processes. RT correlates or anticorrelates with genomic-specific parameters such as gene expression, chromatin accessibility, histone modifications, and the 3D structure of the genome, but the significance of how they influence each other and how they relate to biological processes remains unclear. In this review I discuss the results of recent analyses of RT, the time at which it is remodeled and consolidated during embryogenesis, how it influences development and differentiation, and the regulatory mechanisms and factors involved.

Increased levels of lagging strand polymerase α in an adult stem cell lineage affect replication-coupled histone incorporation

Stem cells display asymmetric histone inheritance, while nonstem progenitor cells exhibit symmetric patterns in the Drosophila male germ line. Here, we report that components involved in lagging strand synthesis, DNA polymerases α and δ, have substantially reduced levels in stem cells compared to progenitor cells, and this promotes local asymmetry of parental histone incorporation at the replication fork. Compromising Polα genetically induces the local replication-coupled histone incorporation pattern in progenitor cells to resemble that in stem cells, seen by both nuclear localization patterns and chromatin fibers. This is recapitulated using a Polα inhibitor in a concentration-dependent manner. The local old versus new histone asymmetry is comparable between stem cells and progenitor cells at both S phase and M phase. Together, these results indicate that developmentally programmed expression of key DNA replication components is important to shape stem cell chromatin. Furthermore, manipulating one crucial DNA replication component can induce replication-coupled histone dynamics in nonstem cells to resemble those in stem cells.

Roles for the 3D genome in the cell cycle, DNA replication, and double strand break repair

Large eukaryotic genomes are packaged into the restricted area of the nucleus to protect the genetic code and provide a dedicated environment to read, copy and repair DNA. The physical organisation of the genome into chromatin loops and self-interacting domains provides the basic structural units of genome architecture. These structural arrangements are complex, multi-layered, and highly dynamic and influence how different regions of the genome interact. The role of chromatin structures during transcription via enhancer-promoter interactions is well established. Less understood is how nuclear architecture influences the plethora of chromatin transactions during DNA replication and repair. In this review, we discuss how genome architecture is regulated during the cell cycle to influence the positioning of replication origins and the coordination of DNA double strand break repair. The role of genome architecture in these cellular processes highlights its critical involvement in preserving genome integrity and cancer prevention.

KMT2C/KMT2D-dependent H3K4me1 mediates changes in DNA replication timing and origin activity during a cell fate transition

Mammalian genomes replicate in a cell-type-specific order during the S phase, correlated to transcriptional activity, histone modifications, and chromatin structure. The causal relationships between these features and DNA replication timing (RT), especially during cell fate changes, are largely unknown. Using machine learning, we quantify 21 chromatin features predicting local RT and RT changes during differentiation in embryonic stem cells (ESCs). About one-third of the genome shows RT changes during differentiation. Chromatin features accurately predict both steady-state RT and RT changes. Histone H3 lysine 4 monomethylation (H3K4me1), catalyzed by KMT2C and KMT2D (KMT2C/D), emerges as a top predictor. Loss of KMT2C/D or their enzymatic activities impairs RT changes during differentiation. This correlates with local H3K4me1 loss and reduced replication origin firing, while transcription remains largely unaffected. Our findings reveal KMT2C/D-dependent H3K4me1 as a key regulator of RT and replication initiation, a role that likely impacts diseases associated with KMT2C/D mutations.

Master transcription factor binding sites constitute the core of early replication control elements

Eukaryotic genomes replicate in a defined temporal order called the replication timing (RT) program. RT is developmentally regulated with potential to drive cell fate transitions, but mechanisms controlling RT remain elusive. We previously identified "Early Replication Control Elements" (ERCEs) necessary for early RT, domain-wide transcription, 3D chromatin architecture and compartmentalization in mouse embryonic stem cells (mESCs) but, deletions identifying ERCEs were large and encompassed many putative regulatory elements. Here, we show that ERCEs are compound elements whose RT activity can largely be accounted for by multiple sites of diverse master transcription factor binding (subERCEs), distinguished from other such sites by their long-range interactions. While deletion of subERCEs had large effects on both transcription and RT, deleting transcription start sites eliminated nearly all transcription with moderate effects on RT. Our results suggest a model in which subERCEs respond to diverse master transcription factors by functioning both as transcription enhancers and as elements that organize chromatin domains structurally and support early RT, potentially providing a feed-forward loop to drive robust epigenomic change during cell fate transitions.

Replication-transcription symbiosis in the mammalian nucleus: The art of living together

Similarly to life in our planet, where thousands of species inhabit the same ecosystem, the cell nucleus hosts different essential processes that share the same territory, making the interaction between them unavoidable. DNA replication and transcription are essential processes that copy and decode the information contained in our genomes, sharing -and competing for- the same chromatin template. Both activities are executed by large macromolecular machines with similar requirements to access the DNA, remodel the nucleosomes ahead of them and reassemble the chromatin make-up behind. Mechanistically, both processes cannot simultaneously act on the same DNA sequence, but emerging evidence shows that they frequently interact. Here we revise recent data on how transcription and replication occur in chromatin highlighting the symbiotic relationship between both processes, which might help explain how their activities contribute to shape the structure and function of the mammalian genome.

Transcription can be sufficient, but is not necessary, to advance replication timing

DNA replication timing (RT) is correlated with transcription during cell fate changes but there are many exceptions and our understanding of this relationship suffers from a paucity of reductionist approaches. Here, we manipulated length and strength of transcription in hybrid-genome mouse embryonic stem cells (mESCs) at a single locus upstream of the silent, late replicating, Pleiotrophin (Ptn) gene, directly comparing RT to nascent transcription rates at engineered vs. wild-type alleles. First, we inserted four reporter genes that differ only in their promoter. Two promoters transcribed the reporter gene at high rates and advanced RT. The other two transcribed at lower rates and did not advance RT. Since these promoters may prove useful in applications where effects on RT are undesirable, we confirmed the inability of one of them to advance RT at numerous ectopic sites. We next juxtaposed these same four promoters upstream of the Ptn transcription start site where they all transcribed the 96kb Ptn gene and advanced RT to different extents correlated with transcription rates. Indeed, a doxycycline-responsive promoter, which could not advance RT when induced as a small reporter gene, elicited a rapid and reversible RT advance proportional to the rate of transcription, providing direct evidence that transcription itself can advance RT. However, deletion of the Ptn promoter and enhancer, followed by directed differentiation to neural precursors, eliminated induction of transcription throughout the entire Ptn replication domain, without preventing the switch to early replication. Our results provide a solid empirical base with which to re-evaluate many decades of literature, demonstrating that length and strength of transcription is sufficient but not necessary to advance RT. Our results also provide a robust system in which to rapidly effect an RT change, permitting mechanistic studies of the role of transcription in RT and the consequences of RT changes to epigenomic remodeling.

Dynamic histone modification patterns coordinating DNA processes

Significant effort has been spent attempting to unravel the causal relationship between histone post-translational modifications and fundamental DNA processes, including transcription, replication, and repair. However, less attention has been paid to understanding the reciprocal influence-that is, how DNA processes, in turn, shape the distribution and patterns of histone modifications and how these changes convey information, both temporally and spatially, from one process to another. Here, we review how histone modifications underpin the widespread bidirectional crosstalk between different DNA processes, which allow seemingly distinct phenomena to operate as a unified whole.

Genome-wide profiling of DNA repair proteins in single cells

Accurate repair of DNA damage is critical for maintenance of genomic integrity and cellular viability. Because damage occurs non-uniformly across the genome, single-cell resolution is required for proper interrogation, but sensitive detection has remained challenging. Here, we present a comprehensive analysis of repair protein localization in single human cells using DamID and ChIC sequencing techniques. This study reports genome-wide binding profiles in response to DNA double-strand breaks induced by AsiSI, and explores variability in genomic damage locations and associated repair features in the context of spatial genome organization. By unbiasedly detecting repair factor localization, we find that repair proteins often occupy entire topologically associating domains, mimicking variability in chromatin loop anchoring. Moreover, we demonstrate the formation of multi-way chromatin hubs in response to DNA damage. Notably, larger hubs show increased coordination of repair protein binding, suggesting a preference for cooperative repair mechanisms. Together, our work offers insights into the heterogeneous processes underlying genome stability in single cells.

The chromosome folding problem and how cells solve it

Every cell must solve the problem of how to fold its genome. We describe how the folded state of chromosomes is the result of the combined activity of multiple conserved mechanisms. Homotypic affinity-driven interactions lead to spatial partitioning of active and inactive loci. Molecular motors fold chromosomes through loop extrusion. Topological features such as supercoiling and entanglements contribute to chromosome folding and its dynamics, and tethering loci to sub-nuclear structures adds additional constraints. Dramatically diverse chromosome conformations observed throughout the cell cycle and across the tree of life can be explained through differential regulation and implementation of these basic mechanisms. We propose that the first functions of chromosome folding are to mediate genome replication, compaction, and segregation and that mechanisms of folding have subsequently been co-opted for other roles, including long-range gene regulation, in different conditions, cell types, and species.

The histone chaperones ASF1 and HIRA are required for telomere length and 45S rDNA copy number homeostasis

Genome stability is significantly influenced by the precise coordination of chromatin complexes that facilitate the loading and eviction of histones from chromatin during replication, transcription, and DNA repair processes. In this study, we investigate the role of the Arabidopsis H3 histone chaperones ANTI-SILENCING FUNCTION 1 (ASF1) and HISTONE REGULATOR A (HIRA) in the maintenance of telomeres and 45S rDNA loci, genomic sites that are particularly susceptible to changes in the chromatin structure. We find that both ASF1 and HIRA are essential for telomere length regulation, as telomeres are significantly shorter in asf1a1b and hira mutants. However, these shorter telomeres remain localized around the nucleolus and exhibit a comparable relative H3 occupancy to the wild type. In addition to regulating telomere length, ASF1 and HIRA contribute to silencing 45S rRNA genes and affect their copy number. Besides, ASF1 supports global heterochromatin maintenance. Our findings also indicate that ASF1 transiently binds to the TELOMERE REPEAT BINDING 1 protein and the N terminus of telomerase in vivo, suggesting a physical link between the ASF1 histone chaperone and the telomere maintenance machinery.

Genome organization and stability in mammalian pre-implantation development

A largely stable genome is required for normal development, even as genetic change is an integral aspect of reproduction, genetic adaptation and evolution. Recent studies highlight a critical window of mammalian development with intrinsic DNA replication stress and genome instability in the first cell divisions after fertilization. Patterns of DNA replication and genome stability are established very early in mammals, alongside patterns of nuclear organization, and before the emergence of gene expression patterns, and prior to cell specification and germline formation. The study of DNA replication and genome stability in the mammalian embryo provides a unique cellular system due to the resetting of the epigenome to a totipotent state, and the de novo establishment of the patterns of nuclear organization, gene expression, DNA methylation, histone modifications and DNA replication. Studies on DNA replication and genome stability in the early mammalian embryo is relevant for understanding both normal and disease-causing genetic variation, and to uncover basic principles of genome regulation.

Structure and dynamics of nuclear A/B compartments and subcompartments

Mammalian chromosomes form a hierarchical structure within the cell nucleus, from chromatin loops, megabase (Mb)-sized topologically associating domains (TADs) to larger-scale A/B compartments. The molecular basis of the structures of loops and TADs has been actively studied. However, the A and B compartments, which correspond to early-replicating euchromatin and late-replicating heterochromatin, respectively, are still relatively unexplored. In this review, we focus on the A/B compartments, discuss their close relationship to DNA replication timing (RT), and introduce recent findings on the features of subcompartments revealed by detailed classification of the A/B compartments. In doing so, we speculate on the structure, potential function, and developmental dynamics of A/B compartments and subcompartments in mammalian cells.

Inferring replication timing and proliferation dynamics from single-cell DNA sequencing data

Dysregulated DNA replication is a cause and a consequence of aneuploidy in cancer, yet the interplay between copy number alterations (CNAs), replication timing (RT) and cell cycle dynamics remain understudied in aneuploid tumors. We developed a probabilistic method, PERT, for simultaneous inference of cell-specific replication and copy number states from single-cell whole genome sequencing (scWGS) data. We used PERT to investigate clone-specific RT and proliferation dynamics in  >50,000 cells obtained from aneuploid and clonally heterogeneous cell lines, xenografts and primary cancers. We observed bidirectional relationships between RT and CNAs, with CNAs affecting X-inactivation producing the largest RT shifts. Additionally, we found that clone-specific S-phase enrichment positively correlated with ground-truth proliferation rates in genomically stable but not unstable cells. Together, these results demonstrate robust computational identification of S-phase cells from scWGS data, and highlight the importance of RT and cell cycle properties in studying the genomic evolution of aneuploid tumors.

Crossing epigenetic frontiers: the intersection of novel histone modifications and diseases

Histone post-translational modifications (HPTMs), as one of the core mechanisms of epigenetic regulation, are garnering increasing attention due to their close association with the onset and progression of diseases and their potential as targeted therapeutic agents. Advances in high-throughput molecular tools and the abundance of bioinformatics data have led to the discovery of novel HPTMs which similarly affect gene expression, metabolism, and chromatin structure. Furthermore, a growing body of research has demonstrated that novel histone modifications also play crucial roles in the development and progression of various diseases, including various cancers, cardiovascular diseases, infectious diseases, psychiatric disorders, and reproductive system diseases. This review defines nine novel histone modifications: lactylation, citrullination, crotonylation, succinylation, SUMOylation, propionylation, butyrylation, 2-hydroxyisobutyrylation, and 2-hydroxybutyrylation. It comprehensively introduces the modification processes of these nine novel HPTMs, their roles in transcription, replication, DNA repair and recombination, metabolism, and chromatin structure, as well as their involvement in promoting the occurrence and development of various diseases and their clinical applications as therapeutic targets and potential biomarkers. Moreover, this review provides a detailed overview of novel HPTM inhibitors targeting various targets and their emerging strategies in the treatment of multiple diseases while offering insights into their future development prospects and challenges. Additionally, we briefly introduce novel epigenetic research techniques and their applications in the field of novel HPTM research.

Pioneer factors for DNA replication initiation in metazoans

Precise DNA replication is fundamental for genetic inheritance. In eukaryotes, replication initiates at multiple origins that are first "licensed" and subsequently "fired" to activate DNA synthesis. Despite the success in identifying origins with specific DNA motifs in Saccharomyces cerevisiae, no consensus sequence or sequences with a predictive value of replication origins have been recognized in metazoan genomes. Rather, epigenetic rules and chromatin structures are believed to play important roles in governing the selection and activation of replication origins. We propose that replication initiation is facilitated by a group of sequence-specific "replication pioneer factors," which function to increase chromatin accessibility and foster a chromatin environment that is conducive to the loading of the prereplication complex. Dysregulation of the function of these factors may lead to gene duplication, genomic instability, and ultimately the occurrence of pathological conditions such as cancer.

Combination of AID2 and BromoTag expands the utility of degron-based protein knockdowns

Acute protein knockdown is a powerful approach to dissecting protein function in dynamic cellular processes. We previously reported an improved auxin-inducible degron system, AID2, but recently noted that its ability to induce degradation of some essential replication factors, such as ORC1 and CDC6, was not enough to induce lethality. Here, we present combinational degron technologies to control two proteins or enhance target depletion. For this purpose, we initially compare PROTAC-based degrons, dTAG and BromoTag, with AID2 to reveal their key features and then demonstrate control of cohesin and condensin with AID2 and BromoTag, respectively. We develop a double-degron system with AID2 and BromoTag to enhance target depletion and accelerate depletion kinetics and demonstrate that both ORC1 and CDC6 are pivotal for MCM loading. Finally, we show that co-depletion of ORC1 and CDC6 by the double-degron system completely suppresses DNA replication, and the cells enter mitosis with single-chromatid chromosomes, indicating that DNA replication is uncoupled from cell cycle control. Our combinational degron technologies will expand the application scope for functional analyses.

Chromosome compartmentalization: causes, changes, consequences, and conundrums

The spatial segregation of the genome into compartments is a major feature of 3D genome organization. New data on mammalian chromosome organization across different conditions reveal important information about how and why these compartments form and change. A combination of epigenetic state, nuclear body tethering, physical forces, gene expression, and replication timing (RT) can all influence the establishment and alteration of chromosome compartments. We review the causes and implications of genomic regions undergoing a 'compartment switch' that changes their physical associations and spatial location in the nucleus. About 20-30% of genomic regions change compartment during cell differentiation or cancer progression, whereas alterations in response to a stimulus within a cell type are usually much more limited. However, even a change in 1-2% of genomic bins may have biologically relevant implications. Finally, we review the effects of compartment changes on gene regulation, DNA damage repair, replication, and the physical state of the cell.

DBF4, not DRF1, is the crucial regulator of CDC7 kinase at replication forks

CDC7 kinase is crucial for DNA replication initiation and is involved in fork processing and replication stress response. Human CDC7 requires the binding of either DBF4 or DRF1 for its activity. However, it is unclear whether the two regulatory subunits target CDC7 to a specific set of substrates, thus having different biological functions, or if they act redundantly. Using genome editing technology, we generated isogenic cell lines deficient in either DBF4 or DRF1: these cells are viable but present signs of genomic instability, indicating that both can independently support CDC7 for bulk DNA replication. Nonetheless, DBF4-deficient cells show altered replication efficiency, partial deficiency in MCM helicase phosphorylation, and alterations in the replication timing of discrete genomic regions. Notably, we find that CDC7 function at replication forks is entirely dependent on DBF4 and not on DRF1. Thus, DBF4 is the primary regulator of CDC7 activity, mediating most of its functions in unperturbed DNA replication and upon replication interference.

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.

Accelerating 3D genomics data analysis with Microcket

The three-dimensional (3D) organization of genome is fundamental to cell biology. To explore 3D genome, emerging high-throughput approaches have produced billions of sequencing reads, which is challenging and time-consuming to analyze. Here we present Microcket, a package for mapping and extracting interacting pairs from 3D genomics data, including Hi-C, Micro-C, and derivant protocols. Microcket utilizes a unique read-stitch strategy that takes advantage of the long read cycles in modern DNA sequencers; benchmark evaluations reveal that Microcket runs much faster than the current tools along with improved mapping efficiency, and thus shows high potential in accelerating and enhancing the biological investigations into 3D genome. Microcket is freely available at https://github.com/hellosunking/Microcket .

Reduced Levels of Lagging Strand Polymerases Shape Stem Cell Chromatin

Stem cells display asymmetric histone inheritance while non-stem progenitor cells exhibit symmetric patterns in the Drosophila male germline lineage. Here, we report that components involved in lagging strand synthesis, such as DNA polymerase α and δ (Polα and Polδ), have significantly reduced levels in stem cells compared to progenitor cells. Compromising Polα genetically induces the replication-coupled histone incorporation pattern in progenitor cells to be indistinguishable from that in stem cells, which can be recapitulated using a Polα inhibitor in a concentration-dependent manner. Furthermore, stem cell-derived chromatin fibers display a higher degree of old histone recycling by the leading strand compared to progenitor cell-derived chromatin fibers. However, upon reducing Polα levels in progenitor cells, the chromatin fibers now display asymmetric old histone recycling just like GSC-derived fibers. The old versus new histone asymmetry is comparable between stem cells and progenitor cells at both S-phase and M-phase. Together, these results indicate that developmentally programmed expression of key DNA replication components is important to shape stem cell chromatin. Furthermore, manipulating one crucial DNA replication component can induce replication-coupled histone dynamics in non-stem cells in a manner similar to that in stem cells.

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.

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.

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.

The nucleolar protein GNL3 prevents resection of stalled replication forks

Faithful DNA replication requires specific proteins that protect replication forks and so prevent the formation of DNA lesions that may damage the genome. Identification of new proteins involved in this process is essential to understand how DNA lesions accumulate in cancer cells and how they tolerate them. Here, we show that human GNL3/nucleostemin, a GTP-binding protein localized mostly in the nucleolus and highly expressed in cancer cells, prevents nuclease-dependent resection of nascent DNA in response to replication stress. We demonstrate that inhibiting origin firing reduces resection. This suggests that the heightened replication origin activation observed upon GNL3 depletion largely drives the observed DNA resection probably due to the exhaustion of the available RPA pool. We show that GNL3 and DNA replication initiation factor ORC2 interact in the nucleolus and that the concentration of GNL3 in the nucleolus is required to limit DNA resection. We propose that the control of origin firing by GNL3 through the sequestration of ORC2 in the nucleolus is critical to prevent nascent DNA resection in response to replication stress.

Diffusion controls local versus dispersed inheritance of histones during replication and shapes epigenomic architecture

The dynamics of inheritance of histones and their associated modifications across cell divisions can have major consequences on maintenance of the cellular epigenomic state. Recent experiments contradict the long-held notion that histone inheritance during replication is always local, suggesting that active and repressed regions of the genome exhibit fundamentally different histone dynamics independent of transcription-coupled turnover. Here we develop a stochastic model of histone dynamics at the replication fork and demonstrate that differential diffusivity of histones in active versus repressed chromatin is sufficient to quantitatively explain these recent experiments. Further, we use the model to predict patterns in histone mark similarity between pairs of genomic loci that should be developed as a result of diffusion, but cannot originate from either PRC2 mediated mark spreading or transcriptional processes. Interestingly, using a combination of CHIP-seq, replication timing and Hi-C datasets we demonstrate that all the computationally predicted patterns are consistently observed for both active and repressive histone marks in two different cell lines. While direct evidence for histone diffusion remains controversial, our results suggest that dislodged histones in euchromatin and facultative heterochromatin may exhibit some level of diffusion within "Diffusion-Accessible-Domains" (DADs), leading to redistribution of epigenetic marks within and across chromosomes. Preservation of the epigenomic state across cell divisions therefore might be achieved not by passing on strict positional information of histone marks, but by maintaining the marks in somewhat larger DADs of the genome.

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.

Epigenetic Regulation in Response to CO2 Fluctuation in Marine Microalga Nannochloropsis oceanica

Microalgae often undergo different CO2 experiment in their habitat. To adapt to low CO2, carbon concentrating mechanism (CCM) could be launched in majority of microalgae and CCM are regulated at RNA level are well known. However, epigenetic modifications and their potential regulation of the transcription of masked genes at the genome level in response to CO2 fluctuation remain unclear. Here epigenetic regulation in response to CO2 fluctuation and epigenome-association with phenotypic plasticity of CCM are firstly uncovered in marine microalga Nannochloropsis oceanica IMET1. The result showed that lysine butyrylation (Kbu) and histone H3K9m2 modifications were present in N. oceanica IMET1. Moreover, Kbu modification positively regulated gene expression. In response to CO2 fluctuation, there were 5,438 and 1,106 genes regulated by Kbu and H3K9m2 in Nannochloropsis, respectively. Gained or lost histone methylations were closely associated with activating or repressing gene expressions. Differential modifications were mainly enriched in carbon fixation, photorespiration, photosynthesis, and lipid metabolism etc. Massive genome-wide epigenetic reprogramming was observed after N. oceanica cells shifted from high CO2 to low CO2. Particularly, we firstly noted that the transcription of the key low CO2 responsive carbonic anhydrase (CA5), a key component involved in CCM stress signaling, was potentially regulated by bivalent Kbu-H3K9m2 modifications in microalgae. This study provides novel insights into the relationship between gene transcription and epigenetic modification in Nannochloropsis, which will lay foundation on genetic improvement of CCM at epigenetic level.

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.

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.

ChArmTelo Enables Large-Scale Chromosome Arm-Level Telomere Analysis across Human Populations and in Cancer Patients

Telomeres are structures protecting chromosome ends. However, a scalable and cost-effective method to investigate chromosome arm-level (ChArm) telomeres (Telos) in large-scale projects is still lacking, hindering intensive investigation of high-resolution telomeres across cancers and other diseases. Here, ChArmTelo, the first computational toolbox to analyze telomeres at chromosome arm level in human and other animal species, using 10X linked-read and similar technologies, is presented. ChArmTelo currently consists of two algorithms, TeloEM and TeloKnow, for arm-level telomere length (TL) analysis. The algorithms are demonstrated by comprehensive analysis of chromosome arm-level telomere lengths (chArmTLs) in nearly 400 whole genome sequencing samples (WGS) from human populations and animals, including healthy and cancer samples. Notably, considerable performance improvement contributed by using the latest complete telomere-to-telomere reference genome (CHM13v2), compared to hg38, is shown. ChArmTelo reveals population-specific chArmTL differences and liver cancer signatures of chArmTLs and that DNA replication origin disruption may contribute to cancer by affecting TLs. Importantly, ChArmTelo can be readily applied to tens of thousands of cancer and healthy samples with published WGS data.

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.

Chromatin and gene expression changes during female Drosophila germline stem cell development illuminate the biology of highly potent stem cells

Highly potent animal stem cells either self renew or launch complex differentiation programs, using mechanisms that are only partly understood. Drosophila female germline stem cells (GSCs) perpetuate without change over evolutionary time and generate cystoblast daughters that develop into nurse cells and oocytes. Cystoblasts initiate differentiation by generating a transient syncytial state, the germline cyst, and by increasing pericentromeric H3K9me3 modification, actions likely to suppress transposable element activity. Relatively open GSC chromatin is further restricted by Polycomb repression of testis or somatic cell-expressed genes briefly active in early female germ cells. Subsequently, Neijre/CBP and Myc help upregulate growth and reprogram GSC metabolism by altering mitochondrial transmembrane transport, gluconeogenesis, and other processes. In all these respects GSC differentiation resembles development of the totipotent zygote. We propose that the totipotent stem cell state was shaped by the need to resist transposon activity over evolutionary timescales.

Single-cell DNA replication dynamics in genomically unstable cancers

Dysregulated DNA replication is both a cause and a consequence of aneuploidy, yet the dynamics of DNA replication in aneuploid cell populations remains understudied. We developed a new method, PERT, for inferring cell-specific DNA replication states from single-cell whole genome sequencing, and investigated clone-specific DNA replication dynamics in >50,000 cells obtained from a collection of aneuploid and clonally heterogeneous cell lines, xenografts and primary cancer tissues. Clone replication timing (RT) profiles correlated with future copy number changes in serially passaged cell lines. Cell type was the strongest determinant of RT heterogeneity, while whole genome doubling and mutational process were associated with accumulation of late S-phase cells and weaker RT associations. Copy number changes affecting chromosome X had striking impact on RT, with loss of the inactive X allele shifting replication earlier, and loss of inactive Xq resulting in reactivation of Xp. Finally, analysis of time series xenografts illustrate how cell cycle distributions approximate clone proliferation, recapitulating expected relationships between proliferation and fitness in treatment-naive and chemotherapeutic contexts.

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.

Spatial and temporal organization of the genome: Current state and future aims of the 4D nucleome project

The four-dimensional nucleome (4DN) consortium studies the architecture of the genome and the nucleus in space and time. We summarize progress by the consortium and highlight the development of technologies for (1) mapping genome folding and identifying roles of nuclear components and bodies, proteins, and RNA, (2) characterizing nuclear organization with time or single-cell resolution, and (3) imaging of nuclear organization. With these tools, the consortium has provided over 2,000 public datasets. Integrative computational models based on these data are starting to reveal connections between genome structure and function. We then present a forward-looking perspective and outline current aims to (1) delineate dynamics of nuclear architecture at different timescales, from minutes to weeks as cells differentiate, in populations and in single cells, (2) characterize cis-determinants and trans-modulators of genome organization, (3) test functional consequences of changes in cis- and trans-regulators, and (4) develop predictive models of genome structure and function.

Regulatory principles and mechanisms governing the onset of random X-chromosome inactivation

X-chromosome inactivation (XCI) has evolved in mammals to compensate for the difference in X-chromosomal dosage between the sexes. In placental mammals, XCI is initiated during early embryonic development through upregulation of the long noncoding RNA Xist from one randomly chosen X chromosome in each female cell. The Xist locus must thus integrate both X-linked and developmental trans-regulatory factors in a dosage-dependent manner. Furthermore, the two alleles must coordinate to ensure inactivation of exactly one X chromosome per cell. In this review, we summarize the regulatory principles that govern the onset of XCI. We go on to provide an overview over the factors that have been implicated in Xist regulation and discuss recent advances in our understanding of how Xist's cis-regulatory landscape integrates information in a precise fashion.

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.

Epigenetic plasticity safeguards heterochromatin configuration in mammals

Heterochromatin is a key architectural feature of eukaryotic chromosomes critical for cell type-specific gene expression and genome stability. In the mammalian nucleus, heterochromatin segregates from transcriptionally active genomic regions and exists in large, condensed, and inactive nuclear compartments. However, the mechanisms underlying the spatial organization of heterochromatin need to be better understood. Histone H3 lysine 9 trimethylation (H3K9me3) and lysine 27 trimethylation (H3K27me3) are two major epigenetic modifications that enrich constitutive and facultative heterochromatin, respectively. Mammals have at least five H3K9 methyltransferases (SUV39H1, SUV39H2, SETDB1, G9a and GLP) and two H3K27 methyltransferases (EZH1 and EZH2). In this study, we addressed the role of H3K9 and H3K27 methylation in heterochromatin organization using a combination of mutant cells for five H3K9 methyltransferases and an EZH1/2 dual inhibitor, DS3201. We showed that H3K27me3, which is normally segregated from H3K9me3, was redistributed to regions targeted by H3K9me3 after the loss of H3K9 methylation and that the loss of both H3K9 and H3K27 methylation resulted in impaired condensation and spatial organization of heterochromatin. Our data demonstrate that the H3K27me3 pathway safeguards heterochromatin organization after the loss of H3K9 methylation in mammalian cells.

Chronic treatment with ATR and CHK1 inhibitors does not substantially increase the mutational burden of human cells

DNA replication stress (RS) entails the frequent slow down and arrest of replication forks by a variety of conditions that hinder accurate and processive genome duplication. Elevated RS leads to genome instability, replication catastrophe and eventually cell death. RS is particularly prevalent in cancer cells and its exacerbation to unsustainable levels by chemotherapeutic agents remains a cornerstone of cancer treatments. The adverse consequences of RS are normally prevented by the ATR and CHK1 checkpoint kinases that stabilize stressed forks, suppress origin firing and promote cell cycle arrest when replication is perturbed. Specific inhibitors of these kinases have been developed and shown to potentiate RS and cell death in multiple in vitro cancer settings. Ongoing clinical trials are now probing their efficacy against various cancer types, either as single agents or in combination with mainstay chemotherapeutics. Despite their promise as valuable additions to the anti-cancer pharmacopoeia, we still lack a genome-wide view of the potential mutagenicity of these new drugs. To investigate this question, we performed chronic long-term treatments of TP53-depleted human cancer cells with ATR and CHK1 inhibitors (ATRi, AZD6738/ceralasertib and CHK1i, MK8776/SCH-900776). ATR or CHK1 inhibition did not significantly increase the mutational burden of cells, nor generate specific mutational signatures. Indeed, no notable changes in the numbers of base substitutions, short insertions/deletions and larger scale rearrangements were observed despite induction of replication-associated DNA breaks during treatments. Interestingly, ATR inhibition did induce a slight increase in closely-spaced mutations, a feature previously attributed to translesion synthesis DNA polymerases. The results suggest that ATRi and CHK1i do not have substantial mutagenic effects in vitro when used as standalone agents.

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.

PARP-dependent and NAT10-independent acetylation of N4-cytidine in RNA appears in UV-damaged chromatin

RNA modifications have been known for many years, but their function has not been fully elucidated yet. For instance, the regulatory role of acetylation on N4-cytidine (ac4C) in RNA can be explored not only in terms of RNA stability and mRNA translation but also in DNA repair. Here, we observe a high level of ac4C RNA at DNA lesions in interphase cells and irradiated cells in telophase. Ac4C RNA appears in the damaged genome from 2 to 45 min after microirradiation. However, RNA cytidine acetyltransferase NAT10 did not accumulate to damaged sites, and NAT10 depletion did not affect the pronounced recruitment of ac4C RNA to DNA lesions. This process was not dependent on the G1, S, and G2 cell cycle phases. In addition, we observed that the PARP inhibitor, olaparib, prevents the recruitment of ac4C RNA to damaged chromatin. Our data imply that the acetylation of N4-cytidine, especially in small RNAs, has an important role in mediating DNA damage repair. Ac4C RNA likely causes de-condensation of chromatin in the vicinity of DNA lesions, making it accessible for other DNA repair factors involved in the DNA damage response. Alternatively, RNA modifications, including ac4C, could be direct markers of damaged RNAs.

Linking genome structures to functions by simultaneous single-cell Hi-C and RNA-seq

Much progress has been made recently in single-cell chromosome conformation capture technologies. However, a method that allows simultaneous profiling of chromatin architecture and gene expression has not been reported. Here, we developed an assay named "Hi-C and RNA-seq employed simultaneously" (HiRES) and performed it on thousands of single cells from developing mouse embryos. Single-cell three-dimensional genome structures, despite being heavily determined by the cell cycle and developmental stages, gradually diverged in a cell type-specific manner as development progressed. By comparing the pseudotemporal dynamics of chromatin interactions with gene expression, we found a widespread chromatin rewiring that occurred before transcription activation. Our results demonstrate that the establishment of specific chromatin interactions is tightly related to transcriptional control and cell functions during lineage specification.

Drug-induced loss of imprinting revealed using bioluminescent reporters of Cdkn1c

Genomic imprinting is an epigenetically mediated mechanism that regulates allelic expression of genes based upon parent-of-origin and provides a paradigm for studying epigenetic silencing and release. Here, bioluminescent reporters for the maternally-expressed imprinted gene Cdkn1c are used to examine the capacity of chromatin-modifying drugs to reverse paternal Cdkn1c silencing. Exposure of reporter mouse embryonic stem cells (mESCs) to 5-Azacytidine, HDAC inhibitors, BET inhibitors or GSK-J4 (KDM6A/B inhibitor) relieved repression of paternal Cdkn1c, either selectively or by inducing biallelic effects. Treatment of reporter fibroblasts with HDAC inhibitors or GSK-J4 resulted in similar paternal Cdkn1c activation, whereas BET inhibitor-induced loss of imprinting was specific to mESCs. Changes in allelic expression were generally not sustained in dividing cultures upon drug removal, indicating that the underlying epigenetic memory of silencing was maintained. In contrast, Cdkn1c de-repression by GSK-J4 was retained in both mESCs and fibroblasts following inhibitor removal, although this impact may be linked to cellular stress and DNA damage. Taken together, these data introduce bioluminescent reporter cells as tools for studying epigenetic silencing and disruption, and demonstrate that Cdkn1c imprinting requires distinct and cell-type specific chromatin features and modifying enzymes to enact and propagate a memory of silencing.