Predictable dynamic program of timing of DNA replication in human cells

Optimized Repli-seq: improved DNA replication timing analysis by next-generation sequencing

The human genome is divided into functional units that replicate at specific times during S-phase. This temporal program is known as replication timing (RT) and is coordinated with the spatial organization of the genome and transcriptional activity. RT is also cell type-specific, dynamically regulated during development, and alterations in RT are observed in multiple diseases. Thus, the precise measure of RT is critical to understand the role of RT in gene function regulation. Distinct methods for assaying the RT program exist; however, conventional methods require thousands of cells as input, prohibiting its applicability to samples with limited cell numbers such as those from disease patients or from early developing embryos. Although single-cell RT analyses have been developed, these methods are low throughput, require generation of numerous libraries, increased sequencing costs, and produce low resolution data. Here, we developed an improved method to measure RT genome-wide that enables high-resolution analysis of low input samples. This method incorporates direct cell sorting into lysis buffer, as well as DNA fragmentation and library preparation in a single tube, resulting in higher yields, increased quality, and reproducibility with decreased costs. We also performed a systematic data processing analysis to provide standardized parameters for RT measurement. This optimized method facilitates RT analysis and will enable its application to a broad range of studies investigating the role of RT in gene expression, nuclear architecture, and disease.

Impact of Chromosomal Context on Origin Selection and the Replication Program

Eukaryotic DNA replication is regulated by conserved mechanisms that bring about a spatial and temporal organization in which distinct genomic domains are copied at characteristic times during S phase. Although this replication program has been closely linked with genome architecture, we still do not understand key aspects of how chromosomal context modulates the activity of replication origins. To address this question, we have exploited models that combine engineered genomic rearrangements with the unique replication programs of post-quiescence and pre-meiotic S phases. Our results demonstrate that large-scale inversions surprisingly do not affect cell proliferation and meiotic progression, despite inducing a restructuring of replication domains on each rearranged chromosome. Remarkably, these alterations in the organization of DNA replication are entirely due to changes in the positions of existing origins along the chromosome, as their efficiencies remain virtually unaffected genome wide. However, we identified striking alterations in origin firing proximal to the fusion points of each inversion, suggesting that the immediate chromosomal neighborhood of an origin is a crucial determinant of its activity. Interestingly, the impact of genome reorganization on replication initiation is highly comparable in the post-quiescent and pre-meiotic S phases, despite the differences in DNA metabolism in these two physiological states. Our findings therefore shed new light on how origin selection and the replication program are governed by chromosomal architecture.

Chemically Induced Chromosomal Interaction (CICI) method to study chromosome dynamics and its biological roles

Numerous intra- and inter-chromosomal contacts have been mapped in eukaryotic genomes, but it remains challenging to link these 3D structures to their regulatory functions. To establish the causal relationships between chromosome conformation and genome functions, we  develop a method, Chemically Induced Chromosomal Interaction (CICI), to selectively perturb the chromosome conformation at targeted loci. Using this method, long-distance chromosomal interactions can be induced dynamically between intra- or inter-chromosomal loci pairs, including the ones with very low Hi-C contact frequencies. Measurement of CICI formation time allows us to probe chromosome encounter dynamics between different loci pairs across the cell cycle. We also conduct two functional tests of CICI. We perturb the chromosome conformation near a DNA double-strand break and observe altered donor preference in homologous recombination; we force interactions between early and late-firing DNA replication origins and find no significant changes in replication timing. These results suggest that chromosome conformation plays a deterministic role in homology-directed DNA repair, but not in the establishment of replication timing. Overall, our study demonstrates that CICI is a powerful tool to study chromosome dynamics and 3D genome function.

Methods to selectively manipulate specific long-distance chromosomal interactions are limited. Here the authors develop a method called Chemically Induced Chromosomal Interaction (CICI) to engineer interactions and demonstrate that 3D conformation plays a causal role in establishing donor DNA preference during DNA repair.

OUP accepted manuscript

Conflicts between transcription and replication machinery are a potent source of replication stress and genome instability; however, no technique currently exists to identify endogenous genomic locations prone to transcription–replication interactions. Here, we report a novel method to identify genomic loci prone to transcription–replication interactions termed transcription–replication immunoprecipitation on nascent DNA sequencing, TRIPn-Seq. TRIPn-Seq employs the sequential immunoprecipitation of RNA polymerase 2 phosphorylated at serine 5 (RNAP2s5) followed by enrichment of nascent DNA previously labeled with bromodeoxyuridine. Using TRIPn-Seq, we mapped 1009 unique transcription–replication interactions (TRIs) in mouse primary B cells characterized by a bimodal pattern of RNAP2s5, bidirectional transcription, an enrichment of RNA:DNA hybrids, and a high probability of forming G-quadruplexes. TRIs are highly enriched at transcription start sites and map to early replicating regions. TRIs exhibit enhanced Replication Protein A association and TRI-associated genes exhibit higher replication fork termination than control transcription start sites, two marks of replication stress. TRIs colocalize with double-strand DNA breaks, are enriched for deletions, and accumulate mutations in tumors. We propose that replication stress at TRIs induces mutations potentially contributing to age-related disease, as well as tumor formation and development.

DNA polymerase zeta contributes to heterochromatin replication to prevent genome instability

The DNA polymerase zeta (Polζ) plays a critical role in bypassing DNA damage. REV3L, the catalytic subunit of Polζ, is also essential in mouse embryonic development and cell proliferation for reasons that remain incompletely understood. In this study, we reveal that REV3L protein interacts with heterochromatin components including repressive histone marks and localizes in pericentromeric regions through direct interaction with HP1 dimer. We demonstrate that Polζ/REV3L ensures progression of replication forks through difficult‐to‐replicate pericentromeric heterochromatin, thereby preventing spontaneous chromosome break formation. We also find that Rev3l‐deficient cells are compromised in the repair of heterochromatin‐associated double‐stranded breaks, eliciting deletions in late‐replicating regions. Lack of REV3L leads to further consequences that may be ascribed to heterochromatin replication and repair‐associated functions of Polζ, with a disruption of the temporal replication program at specific loci. This is correlated with changes in epigenetic landscape and transcriptional control of developmentally regulated genes. These results reveal a new function of Polζ in preventing chromosome instability during replication of heterochromatic regions.

Low Replicative Stress Triggers Cell-Type Specific Inheritable Advanced Replication Timing

DNA replication timing (RT), reflecting the temporal order of origin activation, is known as a robust and conserved cell-type specific process. Upon low replication stress, the slowing of replication forks induces well-documented RT delays associated to genetic instability, but it can also generate RT advances that are still uncharacterized. In order to characterize these advanced initiation events, we monitored the whole genome RT from six independent human cell lines treated with low doses of aphidicolin. We report that RT advances are cell-type-specific and involve large heterochromatin domains. Importantly, we found that some major late to early RT advances can be inherited by the unstressed next-cellular generation, which is a unique process that correlates with enhanced chromatin accessibility, as well as modified replication origin landscape and gene expression in daughter cells. Collectively, this work highlights how low replication stress may impact cellular identity by RT advances events at a subset of chromosomal domains.

Mapping replication timing domains genome wide in single mammalian cells with single-cell DNA replication sequencing

Replication timing (RT) domains are stable units of chromosome structure that are regulated in the context of development and disease. Conventional genome-wide RT mapping methods require many S-phase cells for either the effective enrichment of replicating DNA through bromodeoxyuridine (BrdU) immunoprecipitation or the determination of copy-number differences during S-phase, which precludes their application to non-abundant cell types and single cells. Here, we provide a simple, cost-effective, and robust protocol for single-cell DNA replication sequencing (scRepli-seq). The scRepli-seq methodology relies on whole-genome amplification (WGA) of genomic DNA (gDNA) from single S-phase cells and next-generation sequencing (NGS)-based determination of copy-number differences that arise between replicated and unreplicated DNA. Haplotype-resolved scRepli-seq, which distinguishes pairs of homologous chromosomes within a single cell, is feasible by using single-nucleotide polymorphism (SNP)/indel information. We also provide computational pipelines for quality control, normalization, and binarization of the scRepli-seq data. The experimental portion of this protocol (before sequencing) takes 3 d.

Loss of histone H3.3 results in DNA replication defects and altered origin dynamics in C. elegans

Histone H3.3 is a replication-independent variant of histone H3 with important roles in development, differentiation, and fertility. Here, we show that loss of H3.3 results in replication defects in Caenorhabditis elegans embryos at elevated temperatures. To characterize these defects, we adapt methods to determine replication timing, map replication origins, and examine replication fork progression. Our analysis of the spatiotemporal regulation of DNA replication shows that despite the very rapid embryonic cell cycle, the genome is replicated from early and late firing origins and is partitioned into domains of early and late replication. We find that under temperature stress conditions, additional replication origins become activated. Moreover, loss of H3.3 results in altered replication fork progression around origins, which is particularly evident at stress-activated origins. These replication defects are accompanied by replication checkpoint activation, a delayed cell cycle, and increased lethality in checkpoint-compromised embryos. Our comprehensive analysis of DNA replication in C. elegans reveals the genomic location of replication origins and the dynamics of their firing, and uncovers a role of H3.3 in the regulation of replication origins under stress conditions.

Efficient, quick and easy-to-use DNA replication timing analysis with START-R suite

DNA replication must be faithful and follow a well-defined spatiotemporal program closely linked to transcriptional activity, epigenomic marks, intranuclear structures, mutation rate and cell fate determination. Among the readouts of the spatiotemporal program of DNA replication, replication timing analyses require not only complex and time-consuming experimental procedures, but also skills in bioinformatics. We developed a dedicated Shiny interactive web application, the START-R (Simple Tool for the Analysis of the Replication Timing based on R) suite, which analyzes DNA replication timing in a given organism with high-throughput data. It reduces the time required for generating and analyzing simultaneously data from several samples. It automatically detects different types of timing regions and identifies significant differences between two experimental conditions in ∼15 min. In conclusion, START-R suite allows quick, efficient and easier analyses of DNA replication timing for all organisms. This novel approach can be used by every biologist. It is now simpler to use this method in order to understand, for example, whether 'a favorite gene or protein' has an impact on replication process or, indirectly, on genomic organization (as Hi-C experiments), by comparing the replication timing profiles between wild-type and mutant cell lines.

High-resolution Repli-Seq defines the temporal choreography of initiation, elongation and termination of replication in mammalian cells

BACKGROUND

DNA replication in mammalian cells occurs in a defined temporal order during S phase, known as the replication timing (RT) programme. Replication timing is developmentally regulated and correlated with chromatin conformation and local transcriptional potential. Here, we present RT profiles of unprecedented temporal resolution in two human embryonic stem cell lines, human colon carcinoma line HCT116, and mouse embryonic stem cells and their neural progenitor derivatives.

RESULTS

Fine temporal windows revealed a remarkable degree of cell-to-cell conservation in RT, particularly at the very beginning and ends of S phase, and identified 5 temporal patterns of replication in all cell types, consistent with varying degrees of initiation efficiency. Zones of replication initiation (IZs) were detected throughout S phase and interacted in 3D space preferentially with other IZs of similar firing time. Temporal transition regions were resolved into segments of uni-directional replication punctuated at specific sites by small, inefficient IZs. Sites of convergent replication were divided into sites of termination or large constant timing regions consisting of many synchronous IZs in tandem. Developmental transitions in RT occured mainly by activating or inactivating individual IZs or occasionally by altering IZ firing time, demonstrating that IZs, rather than individual origins, are the units of developmental regulation. Finally, haplotype phasing revealed numerous regions of allele-specific and allele-independent asynchronous replication. Allele-independent asynchronous replication was correlated with the presence of previously mapped common fragile sites.

CONCLUSIONS

Altogether, these data provide a detailed temporal choreography of DNA replication in mammalian cells.

Genomic methods for measuring DNA replication dynamics

Genomic DNA replicates according to a defined temporal program in which early-replicating loci are associated with open chromatin, higher gene density, and increased gene expression levels, while late-replicating loci tend to be heterochromatic and show higher rates of genomic instability. The ability to measure DNA replication dynamics at genome scale has proven crucial for understanding the mechanisms and cellular consequences of DNA replication timing. Several methods, such as quantification of nucleotide analog incorporation and DNA copy number analyses, can accurately reconstruct the genomic replication timing profiles of various species and cell types. More recent developments have expanded the DNA replication genomic toolkit to assays that directly measure the activity of replication origins, while single-cell replication timing assays are beginning to reveal a new level of replication timing regulation. The combination of these methods, applied on a genomic scale and in multiple biological systems, promises to resolve many open questions and lead to a holistic understanding of how eukaryotic cells replicate their genomes accurately and efficiently.

Germline DNA replication timing shapes mammalian genome composition

Mammalian DNA replication is a highly organized and regulated process. Large, Mb-sized regions are replicated at defined times along S-phase. Replication Timing (RT) is thought to play a role in shaping the mammalian genome by affecting mutation rates. Previous analyses relied on somatic RT profiles. However, only germline mutations are passed on to offspring and affect genomic composition. Therefore, germ cell RT information is necessary to evaluate the influences of RT on the mammalian genome. We adapted the RT mapping technique for limited amounts of cells, and measured RT from two stages in the mouse germline - primordial germ cells (PGCs) and spermatogonial stem cells (SSCs). RT in germline cells exhibited stronger correlations to both mutation rate and recombination hotspots density than those of RT in somatic tissues, emphasizing the importance of using correct tissues-of-origin for RT profiling. Germline RT maps exhibited stronger correlations to additional genetic features including GC-content, transposable elements (SINEs and LINEs), and gene density. GC content stratification and multiple regression analysis revealed independent contributions of RT to SINE, gene, mutation, and recombination hotspot densities. Together, our results establish a central role for RT in shaping multiple levels of mammalian genome composition.

Germline Structural Variations Are Preferential Sites of DNA Replication Timing Plasticity during Development

The DNA replication timing program is modulated throughout development and is also one of the main factors influencing the distribution of mutation rates across the genome. However, the relationship between the mutagenic influence of replication timing and its developmental plasticity remains unexplored. Here, we studied the distribution of copy number variations (CNVs) and single nucleotide polymorphisms across the zebrafish genome in relation to changes in DNA replication timing during embryonic development in this model vertebrate species. We show that CNV sites exhibit strong replication timing plasticity during development, replicating significantly early during early development but significantly late during more advanced developmental stages. Reciprocally, genomic regions that changed their replication timing during development contained a higher proportion of CNVs than developmentally constant regions. Developmentally plastic CNV sites, in particular those that become delayed in their replication timing, were enriched for the clustered protocadherins, a set of genes important for neuronal development that have undergone extensive genetic and epigenetic diversification during zebrafish evolution. In contrast, single nucleotide polymorphism sites replicated consistently early throughout embryonic development, highlighting a unique aspect of the zebrafish genome. Our results uncover a hitherto unrecognized interface between development and evolution.

Mechanisms of establishment and functional significance of DNA demethylation during erythroid differentiation

Erythroid differentiation is associated with global DNA demethylation, but a complete methylome was lacking in the erythroid lineage. We have generated allele-specific base resolution methylomes of primary basophilic erythroblasts (BasoEs) and compared these with 8 other cell types. We found that DNA demethylation during differentiation from hematopoietic stem/progenitor cells (HSPCs) to BasoEs occurred predominantly in intergenic sequences and in inactive gene bodies causing the formation of partially methylated domains (PMDs) in 74% of the BasoE methylome. Moreover, differentially methylated regions (DMRs) between HSPCs and BasoEs occurred mostly in putative enhancer regions and were most often associated with GATA, EKLF, and AP1 binding motifs. Surprisingly, promoters silent in both HSPCs and BasoEs exhibited much more dramatic chromatin changes during differentiation than activated promoters. Unmethylated silent promoters were often associated with active chromatin states in highly methylated domains (HMDs) but with polycomb-repression in PMDs, indicating that silent promoters are generally regulated differently in HMDs and PMDs. We show that long PMDs replicate late, but that short PMDs replicate early and therefore that the partial methylation of DNA after replication during erythroid expansion occurs throughout S phase of the cell cycle. We propose that baseline maintenance methylation following replication decreases during erythroid differentiation resulting in PMD formation and that the presence of HMDs in the BasoE methylome results from transcription-associated DNA methylation of gene bodies. We detected ∼700 large allele-specific DMRs that were enriched in single-nucleotide polymorphisms, suggesting that primary DNA sequence might be a determinant of DNA methylation levels within PMDs.

Insights into the Link between the Organization of DNA Replication and the Mutational Landscape

The generation of a complete and accurate copy of the genetic material during each cell cycle is integral to cell growth and proliferation. However, genetic diversity is essential for adaptation and evolution, and the process of DNA replication is a fundamental source of mutations. Genome alterations do not accumulate randomly, with variations in the types and frequencies of mutations that arise in different genomic regions. Intriguingly, recent studies revealed a striking link between the mutational landscape of a genome and the spatial and temporal organization of DNA replication, referred to as the replication program. In our review, we discuss how this program may contribute to shaping the profile and spectrum of genetic alterations, with implications for genome dynamics and organismal evolution in natural and pathological contexts.

DNA Replication Timing Enters the Single-Cell Era

In mammalian cells, DNA replication timing is controlled at the level of megabase (Mb)-sized chromosomal domains and correlates well with transcription, chromatin structure, and three-dimensional (3D) genome organization. Because of these properties, DNA replication timing is an excellent entry point to explore genome regulation at various levels and a variety of studies have been carried out over the years. However, DNA replication timing studies traditionally required at least tens of thousands of cells, and it was unclear whether the replication domains detected by cell population analyses were preserved at the single-cell level. Recently, single-cell DNA replication profiling methods became available, which revealed that the Mb-sized replication domains detected by cell population analyses were actually well preserved in individual cells. In this article, we provide a brief overview of our current knowledge on DNA replication timing regulation in mammals based on cell population studies, outline the findings from single-cell DNA replication profiling, and discuss future directions and challenges.

Linking the organization of DNA replication with genome maintenance

The spatial and temporal organization of genome duplication, also referred to as the replication program, is defined by the distribution and the activities of the sites of replication initiation across the genome. Alterations to the replication profile are associated with cell fate changes during development and in pathologies, but the importance of undergoing S phase with distinct and specific programs remains largely unexplored. We have recently addressed this question, focusing on the interplay between the replication program and genome maintenance. In particular, we demonstrated that when cells encounter challenges to DNA synthesis, the organization of DNA replication drives the response to replication stress that is mediated by the ATR/Rad3 checkpoint pathway, thus shaping the pattern of genome instability along the chromosomes. In this review, we present the major findings of our study and discuss how they may bring new perspectives to our understanding of the biological importance of the replication program.

Regulation of the program of DNA replication by CDK: new findings and perspectives

Progression through the cell cycle is driven by the activities of the cyclin-dependent kinase (CDK) family of enzymes, which establish an ordered passage through the cell cycle phases. CDK activity is crucial for the cellular transitions from G1 to S and G2 to M, which are highly controlled to promote the faithful duplication of the genetic material and the transmission of the genome into daughter cells, respectively. While oscillations in CDK activity are essential for cell division, how its specific dynamics may shape cellular processes remains an open question. Recently, we have investigated the potential role of CDK in establishing the profile of replication initiation along the chromosomes, also referred to as the replication program. Our results demonstrated that the timing and level of CDK activity at G1/S provide two critical and independent inputs that modulate the pattern of origin usage. In this review, we will present the conclusions of our study and discuss the implications of our findings for cellular function and physiology.

Modeling and correct the GC bias of tumor and normal WGS data for SCNA based tumor subclonal population inferring

BACKGROUND

Somatic copy number alternations (SCNAs) can be utilized to infer tumor subclonal populations in whole genome seuqncing studies, where usually their read count ratios between tumor-normal paired samples serve as the inferring proxy. Existing SCNA based subclonal population inferring tools consider the GC bias of tumor and normal sample is of the same fature, and could be fully offset by read count ratio. However, we found that, the read count ratio on SCNA segments presents a Log linear biased pattern, which influence existing read count ratios based subclonal inferring tools performance. Currently no correction tools take into account the read ratio bias.

RESULTS

We present Pre-SCNAClonal, a tool that improving tumor subclonal population inferring by correcting GC-bias at SCNAs level. Pre-SCNAClonal first corrects GC bias using Markov chain Monte Carlo probability model, then accurately locates baseline DNA segments (not containing any SCNAs) with a hierarchy clustering model. We show Pre-SCNAClonal's superiority to exsiting GC-bias correction methods at any level of subclonal population.

CONCLUSIONS

Pre-SCNAClonal could be run independently as well as serving as pre-processing/gc-correction step in conjuntion with exsiting SCNA-based subclonal inferring tools.

Correlative live and super-resolution imaging reveals the dynamic structure of replication domains

Combining fluorescence labeling with live-cell confocal and correlative super-resolution microscopy, Xiang et al. characterize biophysical parameters defining the internal organization, spacing, and mechanical coupling of replication domains.

Chromosome organization in higher eukaryotes controls gene expression, DNA replication, and DNA repair. Genome mapping has revealed the functional units of chromatin at the submegabase scale as self-interacting regions called topologically associating domains (TADs) and showed they correspond to replication domains (RDs). A quantitative structural and dynamic description of RD behavior in the nucleus is, however, missing because visualization of dynamic subdiffraction-sized RDs remains challenging. Using fluorescence labeling of RDs combined with correlative live and super-resolution microscopy in situ, we determined biophysical parameters to characterize the internal organization, spacing, and mechanical coupling of RDs. We found that RDs are typically 150 nm in size and contain four co-replicating regions spaced 60 nm apart. Spatially neighboring RDs are spaced 300 nm apart and connected by highly flexible linker regions that couple their motion only <550 nm. Our pipeline allows a robust quantitative characterization of chromosome structure in situ and provides important biophysical parameters to understand general principles of chromatin organization.

Genome-Wide Analysis of the Arabidopsis Replication Timing Program

Eukaryotes use a temporally regulated process, known as the replication timing program, to ensure that their genomes are fully and accurately duplicated during S phase. Replication timing programs are predictive of genomic features and activity and are considered to be functional readouts of chromatin organization. Although replication timing programs have been described for yeast and animal systems, much less is known about the temporal regulation of plant DNA replication or its relationship to genome sequence and chromatin structure. We used the thymidine analog, 5-ethynyl-2'-deoxyuridine, in combination with flow sorting and Repli-Seq to describe, at high-resolution, the genome-wide replication timing program for Arabidopsis (Arabidopsis thaliana) Col-0 suspension cells. We identified genomic regions that replicate predominantly during early, mid, and late S phase, and correlated these regions with genomic features and with data for chromatin state, accessibility, and long-distance interaction. Arabidopsis chromosome arms tend to replicate early while pericentromeric regions replicate late. Early and mid-replicating regions are gene-rich and predominantly euchromatic, while late regions are rich in transposable elements and primarily heterochromatic. However, the distribution of chromatin states across the different times is complex, with each replication time corresponding to a mixture of states. Early and mid-replicating sequences interact with each other and not with late sequences, but early regions are more accessible than mid regions. The replication timing program in Arabidopsis reflects a bipartite genomic organization with early/mid-replicating regions and late regions forming separate, noninteracting compartments. The temporal order of DNA replication within the early/mid compartment may be modulated largely by chromatin accessibility.

CDK activity provides temporal and quantitative cues for organizing genome duplication

In eukaryotes, the spatial and temporal organization of genome duplication gives rise to distinctive profiles of replication origin usage along the chromosomes. While it has become increasingly clear that these programs are important for cellular physiology, the mechanisms by which they are determined and modulated remain elusive. Replication initiation requires the function of cyclin-dependent kinases (CDKs), which associate with various cyclin partners to drive cell proliferation. Surprisingly, although we possess detailed knowledge of the CDK regulators and targets that are crucial for origin activation, little is known about whether CDKs play a critical role in establishing the genome-wide pattern of origin selection. We have addressed this question in the fission yeast, taking advantage of a simplified cell cycle network in which cell proliferation is driven by a single cyclin-CDK module. This system allows us to precisely control CDK activity in vivo using chemical genetics. First, in contrast to previous reports, our results clearly show that distinct cyclin-CDK pairs are not essential for regulating specific subsets of origins and for establishing a normal replication program. Importantly, we then demonstrate that the timing at which CDK activity reaches the S phase threshold is critical for the organization of replication in distinct efficiency domains, while the level of CDK activity at the onset of S phase is a dose-dependent modulator of overall origin efficiencies. Our study therefore implicates these different aspects of CDK regulation as versatile mechanisms for shaping the architecture of DNA replication across the genome.

The duplication of the genetic material is a highly conserved and tightly regulated process that is essential for cell proliferation. DNA synthesis initiates at sites called origins that are distributed throughout the genome. Replication origins are not all used equivalently, and their patterns of activation along the chromosomes give rise to specific profiles, or programs, of DNA replication. These programs change during development and in response to external stimuli, and we have previously shown that they have important consequences for cellular function. However, we still do not understand the mechanisms by which cells establish different replication patterns. Here we investigate the role of the cyclin-dependent kinase (CDK) family of proteins, whose activities are critical for cell cycle progression, in regulating the organization of genome duplication. Taking advantage of a system that allows us to precisely modulate CDK activity levels in living cells, we demonstrate that both the temporal and quantitative controls of CDK function are crucial for determining distinct programs of DNA replication. Our work therefore uncovers a fundamental link between CDK activity, a central input in a variety of cellular and developmental processes, and the architecture of genome duplication.

Genomic Analysis of the DNA Replication Timing Program during Mitotic S Phase in Maize (Zea mays) Root Tips

All plants and animals must replicate their DNA, using a regulated process to ensure that their genomes are completely and accurately replicated. DNA replication timing programs have been extensively studied in yeast and animal systems, but much less is known about the replication programs of plants. We report a novel adaptation of the "Repli-seq" assay for use in intact root tips of maize (Zea mays) that includes several different cell lineages and present whole-genome replication timing profiles from cells in early, mid, and late S phase of the mitotic cell cycle. Maize root tips have a complex replication timing program, including regions of distinct early, mid, and late S replication that each constitute between 20 and 24% of the genome, as well as other loci corresponding to ∼32% of the genome that exhibit replication activity in two different time windows. Analyses of genomic, transcriptional, and chromatin features of the euchromatic portion of the maize genome provide evidence for a gradient of early replicating, open chromatin that transitions gradually to less open and less transcriptionally active chromatin replicating in mid S phase. Our genomic level analysis also demonstrated that the centromere core replicates in mid S, before heavily compacted classical heterochromatin, including pericentromeres and knobs, which replicate during late S phase.

Distinct epigenetic features of differentiation-regulated replication origins

BACKGROUND

Eukaryotic genome duplication starts at discrete sequences (replication origins) that coordinate cell cycle progression, ensure genomic stability and modulate gene expression. Origins share some sequence features, but their activity also responds to changes in transcription and cellular differentiation status.

RESULTS

To identify chromatin states and histone modifications that locally mark replication origins, we profiled origin distributions in eight human cell lines representing embryonic and differentiated cell types. Consistent with a role of chromatin structure in determining origin activity, we found that cancer and non-cancer cells of similar lineages exhibited highly similar replication origin distributions. Surprisingly, our study revealed that DNase hypersensitivity, which often correlates with early replication at large-scale chromatin domains, did not emerge as a strong local determinant of origin activity. Instead, we found that two distinct sets of chromatin modifications exhibited strong local associations with two discrete groups of replication origins. The first origin group consisted of about 40,000 regions that actively initiated replication in all cell types and preferentially colocalized with unmethylated CpGs and with the euchromatin markers, H3K4me3 and H3K9Ac. The second group included origins that were consistently active in cells of a single type or lineage and preferentially colocalized with the heterochromatin marker, H3K9me3. Shared origins replicated throughout the S-phase of the cell cycle, whereas cell-type-specific origins preferentially replicated during late S-phase.

CONCLUSIONS

These observations are in line with the hypothesis that differentiation-associated changes in chromatin and gene expression affect the activation of specific replication origins.

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

DNA replication is essential for faithful transmission of genetic information and is intimately tied to chromosome structure and function. Genome duplication occurs in a defined temporal order known as the replication-timing (RT) program, which is regulated during the cell cycle and development in discrete units referred to as replication domains (RDs). RDs correspond to topologically-associating domains (TADs) and are spatio-temporally compartmentalized in the nucleus. While improvements in experimental tools have begun to reveal glimpses of causality, they have also unveiled complex context-dependent relationships that challenge long recognized correlations of RT to chromatin organization and gene regulation. In particular, RDs/TADs that switch RT during development march to the beat of a different drummer.

How to rule the nucleus: divide et impera

Genome-wide molecular studies have provided new insights into the organization of nuclear chromatin by revealing the presence of chromatin domains of differing transcriptional activity, frequency of cis-interactions, proximity to scaffolding structures and replication timing. These studies have not only brought our understanding of genome function to a new level, but also offered functional insight for many phenomena observed in microscopic studies. In this review, we discuss the major principles of nuclear organization based on the spatial segregation of euchromatin and heterochromatin, as well as the dynamic genome rearrangements occurring during cell differentiation and development. We hope to unite the existing molecular and microscopic data on genome organization to get a holistic view of the nucleus, and propose a model, in which repeat repertoire together with scaffolding structures blueprint the functional nuclear architecture.

Large-Scale Chromatin Structure-Function Relationships during the Cell Cycle and Development: Insights from Replication Timing

Chromosome architecture has received a lot of attention since the recent development of genome-scale methods to measure chromatin interactions (Hi-C), enabling the first sequence-based models of chromosome tertiary structure. A view has emerged of chromosomes as a string of structural units (topologically associating domains; TADs) whose boundaries persist through the cell cycle and development. TADs with similar chromatin states tend to aggregate, forming spatially segregated chromatin compartments. However, high-resolution Hi-C has revealed substructure within TADs (subTADs) that poses a challenge for models that attribute significance to structural units at any given scale. More than 20 years ago, the DNA replication field independently identified stable structural (and functional) units of chromosomes (replication foci) as well as spatially segregated chromatin compartments (early and late foci), but lacked the means to link these units to genomic map units. Genome-wide studies of replication timing (RT) have now merged these two disciplines by identifying individual units of replication regulation (replication domains; RDs) that correspond to TADs and are arranged in 3D to form spatiotemporally segregated subnuclear compartments. Furthermore, classifying RDs/TADs by their constitutive versus developmentally regulated RT has revealed distinct classes of chromatin organization, providing unexpected insight into the relationship between large-scale chromosome structure and function.

Chromatin at the nuclear periphery and the regulation of genome functions

Chromatin is not randomly organized in the nucleus, and its spatial organization participates in the regulation of genome functions. However, this spatial organization is also not entirely fixed and modifications of chromatin architecture are implicated in physiological processes such as differentiation or senescence. One of the most striking features of chromatin architecture is the concentration of heterochromatin at the nuclear periphery. A closer examination of the association of chromatin at the nuclear periphery reveals that heterochromatin accumulates at the nuclear lamina, whereas nuclear pores are usually devoid of heterochromatin. After summarizing the current techniques used to study the attachment of chromatin at the nuclear lamina or the nuclear pores, we review the mechanisms underlying these attachments, their plasticity and their consequences on the regulation of gene expression, DNA repair and replication.

Dynamic changes in replication timing and gene expression during lineage specification of human pluripotent stem cells

Duplication of the genome in mammalian cells occurs in a defined temporal order referred to as its replication-timing (RT) program. RT changes dynamically during development, regulated in units of 400-800 kb referred to as replication domains (RDs). Changes in RT are generally coordinated with transcriptional competence and changes in subnuclear position. We generated genome-wide RT profiles for 26 distinct human cell types, including embryonic stem cell (hESC)-derived, primary cells and established cell lines representing intermediate stages of endoderm, mesoderm, ectoderm, and neural crest (NC) development. We identified clusters of RDs that replicate at unique times in each stage (RT signatures) and confirmed global consolidation of the genome into larger synchronously replicating segments during differentiation. Surprisingly, transcriptome data revealed that the well-accepted correlation between early replication and transcriptional activity was restricted to RT-constitutive genes, whereas two-thirds of the genes that switched RT during differentiation were strongly expressed when late replicating in one or more cell types. Closer inspection revealed that transcription of this class of genes was frequently restricted to the lineage in which the RT switch occurred, but was induced prior to a late-to-early RT switch and/or down-regulated after an early-to-late RT switch. Analysis of transcriptional regulatory networks showed that this class of genes contains strong regulators of genes that were only expressed when early replicating. These results provide intriguing new insight into the complex relationship between transcription and RT regulation during human development.

Allele-specific analysis of DNA replication origins in mammalian cells

The mechanisms that control the location and timing of firing of replication origins are poorly understood. Using a novel functional genomic approach based on the analysis of SNPs and indels in phased human genomes, we observe that replication asynchrony is associated with small cumulative variations in the initiation efficiency of multiple origins between the chromosome homologues, rather than with the activation of dormant origins. Allele-specific measurements demonstrate that the presence of G-quadruplex-forming sequences does not correlate with the efficiency of initiation. Sequence analysis reveals that the origins are highly enriched in sequences with profoundly asymmetric G/C and A/T nucleotide distributions and are almost completely depleted of antiparallel triplex-forming sequences. We therefore propose that although G4-forming sequences are abundant in replication origins, an asymmetry in nucleotide distribution, which increases the propensity of origins to unwind and adopt non-B DNA structure, rather than the ability to form G4, is directly associated with origin activity.

Structural organization of human replication timing domains

Recent analysis of genome-wide epigenetic modification data, mean replication timing (MRT) profiles and chromosome conformation data in mammals have provided increasing evidence that flexibility in replication origin usage is regulated locally by the epigenetic landscape and over larger genomic distances by the 3D chromatin architecture. Here, we review the recent results establishing some link between replication domains and chromatin structural domains in pluripotent and various differentiated cell types in human. We reconcile the originally proposed dichotomic picture of early and late constant timing regions that replicate by multiple rather synchronous origins in separated nuclear compartments of open and closed chromatins, with the U-shaped MRT domains bordered by "master" replication origins specified by a localized (∼200-300 kb) zone of open and transcriptionally active chromatin from which a replication wave likely initiates and propagates toward the domain center via a cascade of origin firing. We discuss the relationships between these MRT domains, topologically associated domains and lamina-associated domains. This review sheds a new light on the epigenetically regulated global chromatin reorganization that underlies the loss of pluripotency and the determination of differentiation properties.

Ubiquitous human 'master' origins of replication are encoded in the DNA sequence via a local enrichment in nucleosome excluding energy barriers

As the elementary building block of eukaryotic chromatin, the nucleosome is at the heart of the compromise between the necessity of compacting DNA in the cell nucleus and the required accessibility to regulatory proteins. The recent availability of genome-wide experimental maps of nucleosome positions for many different organisms and cell types has provided an unprecedented opportunity to elucidate to what extent the DNA sequence conditions the primary structure of chromatin and in turn participates in the chromatin-mediated regulation of nuclear functions, such as gene expression and DNA replication. In this study, we use in vivo and in vitro genome-wide nucleosome occupancy data together with the set of nucleosome-free regions (NFRs) predicted by a physical model of nucleosome formation based on sequence-dependent bending properties of the DNA double-helix, to investigate the role of intrinsic nucleosome occupancy in the regulation of the replication spatio-temporal programme in human. We focus our analysis on the so-called replication U/N-domains that were shown to cover about half of the human genome in the germline (skew-N domains) as well as in embryonic stem cells, somatic and HeLa cells (mean replication timing U-domains). The 'master' origins of replication (MaOris) that border these megabase-sized U/N-domains were found to be specified by a few hundred kb wide regions that are hyper-sensitive to DNase I cleavage, hypomethylated, and enriched in epigenetic marks involved in transcription regulation, the hallmarks of localized open chromatin structures. Here we show that replication U/N-domain borders that are conserved in all considered cell lines have an environment highly enriched in nucleosome-excluding-energy barriers, suggesting that these ubiquitous MaOris have been selected during evolution. In contrast, MaOris that are cell-type-specific are mainly regulated epigenetically and are no longer favoured by a local abundance of intrinsic NFRs encoded in the DNA sequence. At the smaller few hundred bp scale of gene promoters, CpG-rich promoters of housekeeping genes found nearby ubiquitous MaOris as well as CpG-poor promoters of tissue-specific genes found nearby cell-type-specific MaOris, both correspond to in vivo NFRs that are not coded as nucleosome-excluding-energy barriers. Whereas the former promoters are likely to correspond to high occupancy transcription factor binding regions, the latter are an illustration that gene regulation in human is typically cell-type-specific.

Embryonic stem cell specific "master" replication origins at the heart of the loss of pluripotency

Epigenetic regulation of the replication program during mammalian cell differentiation remains poorly understood. We performed an integrative analysis of eleven genome-wide epigenetic profiles at 100 kb resolution of Mean Replication Timing (MRT) data in six human cell lines. Compared to the organization in four chromatin states shared by the five somatic cell lines, embryonic stem cell (ESC) line H1 displays (i) a gene-poor but highly dynamic chromatin state (EC4) associated to histone variant H2AZ rather than a HP1-associated heterochromatin state (C4) and (ii) a mid-S accessible chromatin state with bivalent gene marks instead of a polycomb-repressed heterochromatin state. Plastic MRT regions (≲ 20% of the genome) are predominantly localized at the borders of U-shaped timing domains. Whereas somatic-specific U-domain borders are gene-dense GC-rich regions, 31.6% of H1-specific U-domain borders are early EC4 regions enriched in pluripotency transcription factors NANOG and OCT4 despite being GC poor and gene deserts. Silencing of these ESC-specific “master” replication initiation zones during differentiation corresponds to a loss of H2AZ and an enrichment in H3K9me3 mark characteristic of late replicating C4 heterochromatin. These results shed a new light on the epigenetically regulated global chromatin reorganization that underlies the loss of pluripotency and lineage commitment.

During development, embryonic stem cell (ESC) enter a program of cell differentiation eventually leading to all the necessary differentiated cell types. Understanding the mechanisms responsible for the underlying modifications of the gene expression program is of fundamental importance, as it will likely have strong impact on the development of regenerative medicine. We show that besides some epigenetic regulation, ubiquitous master replication origins at replication timing U-domain borders shared by 6 human cell types are transcriptionally active open chromatin regions specified by a local enrichment in nucleosome free regions encoded in the DNA sequence suggesting that they have been selected during evolution. In contrast, ESC specific master replication origins bear a unique epigenetic signature (enrichment in CTCF, H2AZ, NANOG, OCT4, …) likely contributing to maintain ESC chromatin in a highly dynamic and accessible state that is refractory to polycomb and HP1 heterochromatin spreading. These ESC specific master origins thus appear as key genomic regions where epigenetic control of chromatin organization is at play to maintain pluripotency of stem cell lineages and to guide lineage commitment to somatic cell types.

Genome architecture and its roles in human copy number variation

Besides single-nucleotide variants in the human genome, large-scale genomic variants, such as copy number variations (CNVs), are being increasingly discovered as a genetic source of human diversity and the pathogenic factors of diseases. Recent experimental findings have shed light on the links between different genome architectures and CNV mutagenesis. In this review, we summarize various genomic features and discuss their contributions to CNV formation. Genomic repeats, including both low-copy and high-copy repeats, play important roles in CNV instability, which was initially known as DNA recombination events. Furthermore, it has been found that human genomic repeats can also induce DNA replication errors and consequently result in CNV mutations. Some recent studies showed that DNA replication timing, which reflects the high-order information of genomic organization, is involved in human CNV mutations. Our review highlights that genome architecture, from DNA sequence to high-order genomic organization, is an important molecular factor in CNV mutagenesis and human genomic instability.

Best practices for mapping replication origins in eukaryotic chromosomes

Understanding the regulatory principles ensuring complete DNA replication in each cell division is critical for deciphering the mechanisms that maintain genomic stability. Recent advances in genome sequencing technology facilitated complete mapping of DNA replication sites and helped move the field from observing replication patterns at a handful of single loci to analyzing replication patterns genome-wide. These advances address issues, such as the relationship between replication initiation events, transcription, and chromatin modifications, and identify potential replication origin consensus sequences. This unit summarizes the technological and fundamental aspects of replication profiling and briefly discusses novel insights emerging from mining large datasets, published in the last 3 years, and also describes DNA replication dynamics on a whole-genome scale.

Large replication skew domains delimit GC-poor gene deserts in human

Besides their large-scale organization in isochores, mammalian genomes display megabase-sized regions, spanning both genes and intergenes, where the strand nucleotide composition asymmetry decreases linearly, possibly due to replication activity. These so-called skew-N domains cover about a third of the human genome and are bordered by two skew upward jumps that were hypothesized to compose a subset of "master" replication origins active in the germline. Skew-N domains were shown to exhibit a particular gene organization. Genes with CpG-rich promoters likely expressed in the germline are over represented near the master replication origins, with large genes being co-oriented with replication fork progression, which suggests some coordination of replication and transcription. In this study, we describe another skew structure that covers ∼13% of the human genome and that is bordered by putative master replication origins similar to the ones flanking skew-N domains. These skew-split-N domains have a shape reminiscent of a N, but split in half, leaving in the center a region of null skew whose length increases with domain size. These central regions (median size ∼860 kb) have a homogeneous composition, i.e. both a null and constant skew and a constant and low GC content. They correspond to heterochromatin gene deserts found in low-GC isochores with an average gene density of 0.81 promoters/Mb as compared to 7.73 promoters/Mb genome wide. The analysis of epigenetic marks and replication timing data confirms that, in these late replicating heterochomatic regions, the initiation of replication is likely to be random. This contrasts with the transcriptionally active euchromatin state found around the bordering well positioned master replication origins. Altogether skew-N domains and skew-split-N domains cover about 50% of the human genome.

DNA replication timing: Coordinating genome stability with genome regulation on the X chromosome and beyond

Recent studies based on next-generation DNA sequencing have revealed that the female inactive X chromosome is replicated in a rapid, unorganized manner, and undergoes increased rates of mutation. These observations link the organization of DNA replication timing to gene regulation on one hand, and to the generation of mutations on the other hand. More generally, the exceptional biology of the inactive X chromosome highlights general principles of genome replication. Cells may control replication timing by a combination of intrinsic replication origin properties, local chromatin states and global levels of replication factors, leading to a functional separation between the activity of genes and their mutation.

Different rates of DNA replication at early versus late S-phase sections: multiscale modeling of stochastic events related to DNA content/EdU (5-ethynyl-2'deoxyuridine) incorporation distributions

Mathematical modeling allows relating molecular events to single-cell characteristics assessed by multiparameter cytometry. In the present study we labeled newly synthesized DNA in A549 human lung carcinoma cells with 15-120 min pulses of EdU. All DNA was stained with DAPI and cellular fluorescence was measured by laser scanning cytometry. The frequency of cells in the ascending (left) side of the "horseshoe"-shaped EdU/DAPI bivariate distributions reports the rate of DNA replication at the time of entrance to S phase while their frequency in the descending (right) side is a marker of DNA replication rate at the time of transition from S to G2 phase. To understand the connection between molecular-scale events and scatterplot asymmetry, we developed a multiscale stochastic model, which simulates DNA replication and cell cycle progression of individual cells and produces in silico EdU/DAPI scatterplots. For each S-phase cell the time points at which replication origins are fired are modeled by a non-homogeneous Poisson Process (NHPP). Shifted gamma distributions are assumed for durations of cell cycle phases (G1, S and G2 M), Depending on the rate of DNA synthesis being an increasing or decreasing function, simulated EdU/DAPI bivariate graphs show predominance of cells in left (early-S) or right (late-S) side of the horseshoe distribution. Assuming NHPP rate estimated from independent experiments, simulated EdU/DAPI graphs are nearly indistinguishable from those experimentally observed. This finding proves consistency between the S-phase DNA-replication rate based on molecular-scale analyses, and cell population kinetics ascertained from EdU/DAPI scatterplots and demonstrates that DNA replication rate at entrance to S is relatively slow compared with its rather abrupt termination during S to G2 transition. Our approach opens a possibility of similar modeling to study the effect of anticancer drugs on DNA replication/cell cycle progression and also to quantify other kinetic events that can be measured during S-phase.

Allele-specific genome-wide profiling in human primary erythroblasts reveal replication program organization

We have developed a new approach to characterize allele-specific timing of DNA replication genome-wide in human primary basophilic erythroblasts. We show that the two chromosome homologs replicate at the same time in about 88% of the genome and that large structural variants are preferentially associated with asynchronous replication. We identified about 600 megabase-sized asynchronously replicated domains in two tested individuals. The longest asynchronously replicated domains are enriched in imprinted genes suggesting that structural variants and parental imprinting are two causes of replication asynchrony in the human genome. Biased chromosome X inactivation in one of the two individuals tested was another source of detectable replication asynchrony. Analysis of high-resolution TimEX profiles revealed small variations termed timing ripples, which were undetected in previous, lower resolution analyses. Timing ripples reflect highly reproducible, variations of the timing of replication in the 100 kb-range that exist within the well-characterized megabase-sized replication timing domains. These ripples correspond to clusters of origins of replication that we detected using novel nascent strands DNA profiling methods. Analysis of the distribution of replication origins revealed dramatic differences in initiation of replication frequencies during S phase and a strong association, in both synchronous and asynchronous regions, between origins of replication and three genomic features: G-quadruplexes, CpG Islands and transcription start sites. The frequency of initiation in asynchronous regions was similar in the two homologs. Asynchronous regions were richer in origins of replication than synchronous regions.

DNA replication in mammalian cells proceeds according to a distinct order. Genes that are expressed tend to replicate before genes that are not expressed. We report here that we have developed a method to measure the timing of replication of the maternal and paternal chromosomes separately. We found that the paternal and maternal chromosomes replicate at exactly the same time in the large majority of the genome and that the 12% of the genome that replicated asynchronously was enriched in imprinted genes and in structural variants. Previous experiments have shown that chromosomes could be divided into replication timing domains that are a few hundred thousand to a few megabases in size. We show here that these domains can be divided into sub-domains defined by ripples in the timing profile. These ripples corresponded to clusters of origins of replication. Finally, we show that the frequency of initiation in asynchronous regions was similar in the two homologs.

The programme of DNA replication: beyond genome duplication

The accurate duplication and transmission of genetic information is critical for cell growth and proliferation, and this is ensured in part by the multi-layered regulation of DNA synthesis. One of the key steps in this process is the selection and activation of the sites of replication initiation, or origins, across the genome. Interestingly, origin usage changes during development and in different pathologies, suggesting an integral interplay between the establishment of replication initiation along the chromosomes and cellular function. The present review discusses how the spatiotemporal organization of replication origin activation may play crucial roles in the control of biological events.

General and specific replication profiles are detected in normal human cells by genome-wide and single-locus molecular combing

Mammalian genomes are replicated under a flexible program, with random use of origins and variable fork rates, and many details of the process must be still unraveled. Molecular combing provides a set of direct data regarding the replication profile of eukaryotic cells: fork rates; organization of the replication clusters; proportion of unidirectional forks; and fork dynamics. In this study the replication profiles of different primary and immortalized non-cancer human cells (lymphocytes, lymphoblastoid cells, fibroblasts) were evaluated at the whole-genome level or within reference genomic regions harboring coding genes. It emerged that these different cell types are characterized by specific replication profiles. In primary fibroblasts, a remarkable fraction of the mammalian genome was found to be replicated by unidirectional forks, and interestingly, the proportion of unidirectional forks further increased in the replicating genome along the population divisions. A second difference concerned in the proportion of paused replication forks, again more frequent in primary fibroblasts than in PBL/lymphoblastoid cells. We concluded that these patterns, whose relevance could escape when genomic methods are applied, represent normal replication features. In single-locus analyses, unidirectional and paused replication forks were highly represented in all genomic regions considered with respect to the average estimates referring to the whole-genome. In addition, fork rates were significantly lower than whole-genome estimates. Instead, when considering the specificities of each genomic region investigated (early to late replication, normal or fragile site) no further differentiating features of replication profiles were detected. These data, representing the integration of genome-wide and single-locus analyses, highlight a large heterogeneity of replication profiles among cell types and within the genome, which should be considered for the correct use of replication datasets.

Human genome replication proceeds through four chromatin states

Advances in genomic studies have led to significant progress in understanding the epigenetically controlled interplay between chromatin structure and nuclear functions. Epigenetic modifications were shown to play a key role in transcription regulation and genome activity during development and differentiation or in response to the environment. Paradoxically, the molecular mechanisms that regulate the initiation and the maintenance of the spatio-temporal replication program in higher eukaryotes, and in particular their links to epigenetic modifications, still remain elusive. By integrative analysis of the genome-wide distributions of thirteen epigenetic marks in the human cell line K562, at the 100 kb resolution of corresponding mean replication timing (MRT) data, we identify four major groups of chromatin marks with shared features. These states have different MRT, namely from early to late replicating, replication proceeds though a transcriptionally active euchromatin state (C1), a repressive type of chromatin (C2) associated with polycomb complexes, a silent state (C3) not enriched in any available marks, and a gene poor HP1-associated heterochromatin state (C4). When mapping these chromatin states inside the megabase-sized U-domains (U-shaped MRT profile) covering about 50% of the human genome, we reveal that the associated replication fork polarity gradient corresponds to a directional path across the four chromatin states, from C1 at U-domains borders followed by C2, C3 and C4 at centers. Analysis of the other genome half is consistent with early and late replication loci occurring in separate compartments, the former correspond to gene-rich, high-GC domains of intermingled chromatin states C1 and C2, whereas the latter correspond to gene-poor, low-GC domains of alternating chromatin states C3 and C4 or long C4 domains. This new segmentation sheds a new light on the epigenetic regulation of the spatio-temporal replication program in human and provides a framework for further studies in different cell types, in both health and disease.

Previous studies revealed spatially coherent and biological-meaningful chromatin mark combinations in human cells. Here, we analyze thirteen epigenetic mark maps in the human cell line K562 at 100 kb resolution of MRT data. The complexity of epigenetic data is reduced to four chromatin states that display remarkable similarities with those reported in fly, worm and plants. These states have different MRT: (C1) is transcriptionally active, early replicating, enriched in CTCF; (C2) is Polycomb repressed, mid-S replicating; (C3) lacks of marks and replicates late and (C4) is a late-replicating gene-poor HP1 repressed heterochromatin state. When mapping these states inside the 876 replication U-domains of K562, the replication fork polarity gradient observed in these U-domains comes along with a remarkable epigenetic organization from C1 at U-domain borders to C2, C3 and ultimately C4 at centers. The remaining genome half displays early replicating, gene rich and high GC domains of intermingled C1 and C2 states segregating from late replicating, gene poor and low GC domains of concatenated C3 and/or C4 states. This constitutes the first evidence of epigenetic compartmentalization of the human genome into replication domains likely corresponding to autonomous units in the 3D chromatin architecture.

DNA replication timing

Patterns of replication within eukaryotic genomes correlate with gene expression, chromatin structure, and genome evolution. Recent advances in genome-scale mapping of replication kinetics have allowed these correlations to be explored in many species, cell types, and growth conditions, and these large data sets have allowed quantitative and computational analyses. One striking new correlation to emerge from these analyses is between replication timing and the three-dimensional structure of chromosomes. This correlation, which is significantly stronger than with any single histone modification or chromosome-binding protein, suggests that replication timing is controlled at the level of chromosomal domains. This conclusion dovetails with parallel work on the heterogeneity of origin firing and the competition between origins for limiting activators to suggest a model in which the stochastic probability of individual origin firing is modulated by chromosomal domain structure to produce patterns of replication. Whether these patterns have inherent biological functions or simply reflect higher-order genome structure is an open question.

[Regulation of DNA replication timing]

Although distinct chromatin types have been long known to replicate at different timepoints of S phase, fine replication control has only recently become considered as an epigenetic phenomenon. It is now clear that in course of differentiation significant changes in genome replication timing occur, and these changes are intimately linked with the changes in transcriptional activity and nuclear architecture. Temporally coordinate replication is organized spatially into discrete units having specific chromosomal organization and function. Even though the functional aspects of such tight control of replication timing remain to be explored, one can confidently consider the replication program as yet another fundamental feature characteristic of the given differentiation state. The present review touches upon the molecular mechanisms of spatial and temporal control of replication timing, involving individual replication origins as well as large chromatin domains.

The spatial organization of the human genome

In vivo, the human genome functions as a complex, folded, three-dimensional chromatin polymer. Understanding how the human genome is spatially organized and folded inside the cell nucleus is therefore central to understanding how genes are regulated in normal development and dysregulated in disease. Established light microscopy-based approaches and more recent molecular chromosome conformation capture methods are now combining to give us unprecedented insight into this fascinating aspect of human genomics.

The replication domain model: regulating replicon firing in the context of large-scale chromosome architecture

The "Replicon Theory" of Jacob, Brenner, and Cuzin has reliably served as the paradigm for regulating the sites where individual replicons initiate replication. Concurrent with the replicon model was Taylor's demonstration that plant and animal chromosomes replicate segmentally in a defined temporal sequence, via cytologically defined units too large to be accounted for by a single replicon. Instead, there seemed to be a program to choreograph when chromosome units replicate during S phase, executed by initiation at clusters of individual replicons within each segment. Here, we summarize recent molecular evidence for the existence of such units, now known as "replication domains", and discuss how the organization of large chromosomes into structural units has added additional layers of regulation to the original replicon model.

Transcription profiling during the cell cycle shows that a subset of Polycomb-targeted genes is upregulated during DNA replication

Genome-wide gene expression analyses of the human somatic cell cycle have indicated that the set of cycling genes differ between primary and cancer cells. By identifying genes that have cell cycle dependent expression in HaCaT human keratinocytes and comparing these with previously identified cell cycle genes, we have identified three distinct groups of cell cycle genes. First, housekeeping genes enriched for known cell cycle functions; second, cell type-specific genes enriched for HaCaT-specific functions; and third, Polycomb-regulated genes. These Polycomb-regulated genes are specifically upregulated during DNA replication, and consistent with being epigenetically silenced in other cell cycle phases, these genes have lower expression than other cell cycle genes. We also find similar patterns in foreskin fibroblasts, indicating that replication-dependent expression of Polycomb-silenced genes is a prevalent but unrecognized regulatory mechanism.