Temporal profile of replication of human chromosomes

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.

Novel candidate biomarkers of origin recognition complex 1, 5 and 6 for survival surveillance in patients with hepatocellular carcinoma

Background: Hepatocellular carcinoma (HCC) has high morbidity and mortality and lacks effective biomarkers for early diagnosis and survival surveillance. Origin recognition complex (ORC), consisting of ORC1-6 isoforms, was examined to assess the potential significance of ORC isoforms for HCC prognosis. Methods: Oncomine and Gene Expression Profiling Interactive Analysis (GEPIA) databases were used to examine differential isoform expression, stage-specific expression, calculate Pearson correlations and perform survival analysis. A human protein atlas database was utilized to evaluate the protein expression of ORCs in liver tissue. The cBioPortal database was used to assess isoform mutations and the survival significance of ORCs in HCC. Cytoscape software was employed to construct gene ontologies, metabolic pathways and gene-gene interaction networks. Results: Differential expression analysis indicated that ORC1 and ORC3-6 were highly expressed in tumor tissues in the Oncomine and GEPIA databases, while ORC2 was not. All the ORCs were showed positive and statistically significant correlations with each other (all P<0.001). ORC1-2 and ORC4-6 expressions were associated with disease stages I-IV (all P<0.05), but ORC3 was not. Survival analysis found that ORC1 and ORC4-6 expressions were associated with overall survival (OS), and ORC1-3 and ORC5-6 expression were associated with recurrence-free survival (RFS; all P<0.05). In addition, low expression of these ORC genes consistently indicated better prognosis compared with high expression. Protein expression analysis revealed that ORC1 and ORC3-6 were expressed in normal liver tissues, whereas ORC2 was not. Enrichment analysis indicated that ORCs were associated with DNA metabolic process, sequence-specific DNA binding and were involved in DNA replication, cell cycle, E2F-enabled inhibition of pre-replication complex formation and G1/S transition. Conclusions: Differentially expressed ORC1, 5 and 6 are candidate biomarkers for survival prediction and recurrence surveillance in HCC.

Next-Generation Sequencing Enables Spatiotemporal Resolution of Human Centromere Replication Timing

Centromeres serve a critical function in preserving genome integrity across sequential cell divisions, by mediating symmetric chromosome segregation. The repetitive, heterochromatic nature of centromeres is thought to be inhibitory to DNA replication, but has also led to their underrepresentation in human reference genome assemblies. Consequently, centromeres have been excluded from genomic replication timing analyses, leaving their time of replication unresolved. However, the most recent human reference genome, hg38, included models of centromere sequences. To establish the experimental requirements for achieving replication timing profiles for centromeres, we sequenced G₁- and S-phase cells from five human cell lines, and aligned the sequence reads to hg38. We were able to infer DNA replication timing profiles for the centromeres in each of the five cell lines, which showed that centromere replication occurs in mid-to-late S phase. Furthermore, we found that replication timing was more variable between cell lines in the centromere regions than expected, given the distribution of variation in replication timing genome-wide. These results suggest the potential of these, and future, sequence models to enable high-resolution studies of replication in centromeres and other heterochromatic regions.

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.

Eukaryotic DNA replication: Orchestrated action of multi-subunit protein complexes

Genome duplication is an essential process to preserve genetic information between generations. The eukaryotic cell cycle is composed of functionally distinct phases: G1, S, G2, and M. One of the key replicative proteins that participate at every stage of DNA replication is the Mcm2-7 complex, a replicative helicase. In the G1 phase, inactive Mcm2-7 complexes are loaded on the replication origins by replication-initiator proteins, ORC and Cdc6. Two kinases, S-CDK and DDK, convert the inactive origin-loaded Mcm2-7 complex to an active helicase, the CMG complex in the S phase. The activated CMG complex begins DNA unwinding and recruits enzymes essential for DNA synthesis to assemble a replisome at the replication fork. After completion of DNA synthesis, the inactive CMG complex on the replicated DNA is removed from chromatin to terminate DNA replication. In this review, we will discuss the structure, function, and regulation of the molecular machines involved in each step of DNA replication.

Identification of the ORC Complex Subunits That Can Interact with the ENY2 Protein of Drosophila melanogaster

The interaction of the Drosophila ENY2 protein with the ORC complex subunits was investigated. It is found that ORC4 and ORC6 subunits directly interact with ENY2.

DNA replication origins-where do we begin?

In this review, Prioleau and MacAlpine summarize recent advances in our understanding of how primary sequence, chromatin environment, and nuclear architecture contribute to the dynamic selection and activation of replication origins across diverse cell types and developmental stages.

For more than three decades, investigators have sought to identify the precise locations where DNA replication initiates in mammalian genomes. The development of molecular and biochemical approaches to identify start sites of DNA replication (origins) based on the presence of defining and characteristic replication intermediates at specific loci led to the identification of only a handful of mammalian replication origins. The limited number of identified origins prevented a comprehensive and exhaustive search for conserved genomic features that were capable of specifying origins of DNA replication. More recently, the adaptation of origin-mapping assays to genome-wide approaches has led to the identification of tens of thousands of replication origins throughout mammalian genomes, providing an unprecedented opportunity to identify both genetic and epigenetic features that define and regulate their distribution and utilization. Here we summarize recent advances in our understanding of how primary sequence, chromatin environment, and nuclear architecture contribute to the dynamic selection and activation of replication origins across diverse cell types and developmental stages.

A chimeric protein composed of NuMA fused to the DNA binding domain of LANA is sufficient for the ori-P-dependent DNA replication

The Kaposi's sarcoma-associated herpesvirus (KSHV) genome is stably maintained in KSHV-infected PEL cell lines during cell division. We previously showed that accumulation of LANA in the nuclear matrix fraction could be important for the latent DNA replication, and that the functional significance of LANA should be its recruitment of ori-P to the nuclear matrix. Here, we investigated whether the forced localization of the LANA-DNA binding domain (DBD) to the nuclear matrix facilitated ori-P-containing plasmid replication. We demonstrated that chimeric proteins constructed by fusion of LANA DBD with the nuclear mitotic apparatus protein (NuMA), which is one of the components of the nuclear matrix, could bind with ori-P and enhance replication of an ori-P-containing plasmid, compared with that in the presence of DBD alone. These results further suggested that the ori-P recruitment to the nuclear matrix through the binding with DBD is important for latent viral DNA replication.

Temporal separation of replication and transcription during S-phase progression

Transcriptional events during S-phase are critical for cell cycle progression. Here, by using a nascent RNA capture assay coupled with high-throughput sequencing, we determined the temporal patterns of transcriptional events that occur during S-phase. We show that genes involved in critical S-phase-specific biological processes such as nucleosome assembly and DNA repair have temporal transcription patterns across S-phase that are not evident from total RNA levels. By comparing transcription timing with replication timing in S-phase, we show that early replicating genes show increased transcription late in S-phase whereas late replicating genes are predominantly transcribed early in S-phase. Global anti-correlation between replication and transcription timing was observed only based on nascent RNA but not total RNA. Our data provides a detailed view of ongoing transcriptional events during the S-phase of cell cycle, and supports that transcription and replication are temporally separated.

DNA replication and transcription programs respond to the same chromatin cues

DNA replication is a dynamic process that occurs in a temporal order along each of the chromosomes. A consequence of the temporally coordinated activation of replication origins is the establishment of broad domains (>100 kb) that replicate either early or late in S phase. This partitioning of the genome into early and late replication domains is important for maintaining genome stability, gene dosage, and epigenetic inheritance; however, the molecular mechanisms that define and establish these domains are poorly understood. The modENCODE Project provided an opportunity to investigate the chromatin features that define the Drosophila replication timing program in multiple cell lines. The majority of early and late replicating domains in the Drosophila genome were static across all cell lines; however, a small subset of domains was dynamic and exhibited differences in replication timing between the cell lines. Both origin selection and activation contribute to defining the DNA replication program. Our results suggest that static early and late replicating domains were defined at the level of origin selection (ORC binding) and likely mediated by chromatin accessibility. In contrast, dynamic domains exhibited low ORC densities in both cell types, suggesting that origin activation and not origin selection governs the plasticity of the DNA replication program. Finally, we show that the male-specific early replication of the X chromosome is dependent on the dosage compensation complex (DCC), suggesting that the transcription and replication programs respond to the same chromatin cues. Specifically, MOF-mediated hyperacetylation of H4K16 on the X chromosome promotes both the up-regulation of male-specific transcription and origin activation.

Next-generation sequencing and DNA replication in human cells: the future has arrived

Accurate regulation of DNA replication ensures faithful transmission of eukaryotic genomes and maintenance of genomic stability and chromatin organization. However, by itself the replication process is a threat for both DNA and chromatin integrity. This becomes particularly relevant in cancer cells, where activated oncogenes induce replication-stress, including unscheduled initiation, fork stalling and collapse and, ultimately, genomic instability. Studies addressing the relationship between (epi)genome integrity and disease have been hampered by our poor knowledge of the mechanisms regulating where and when eukaryotic replication initiates. Recently developed genome-scale methods for the analysis of DNA replication in mammals will contribute to the identification of missing links between replication, chromatin regulation and genome stability in normal and cancer cells.

A new light on DNA replication from the inactive X chromosome

While large portions of the mammalian genome are known to replicate sequentially in a distinct, tissue-specific order, recent studies suggest that the inactive X chromosome is duplicated rapidly via random, synchronous DNA synthesis at numerous adjacent regions. The rapid duplication of the inactive X chromosome was observed in high-resolution studies visualizing DNA replication patterns in the nucleus, and by allele-specific DNA sequencing studies measuring the extent of DNA synthesis. These studies conclude that inactive X chromosomes complete replication earlier than previously thought and suggest that the strict order of DNA replication detected in the majority of genomic regions is not preserved in non-transcribed, "silent" chromatin. These observations alter current concepts about the regulation of DNA replication in non-transcribed portions of the genome in general and in the inactive X-chromosome in particular.

Bubble-seq analysis of the human genome reveals distinct chromatin-mediated mechanisms for regulating early- and late-firing origins

We have devised a method for isolating virtually pure and comprehensive libraries of restriction fragments that contained replication initiation sites (bubbles) in vivo. We have now sequenced and mapped the bubble-containing fragments from GM06990, a near-normal EBV-transformed lymphoblastoid cell line, and have compared origin distributions with a comprehensive replication timing study recently published for this cell line. We find that early-firing origins, which represent ∼32% of all origins, overwhelmingly represent zones, associate only marginally with active transcription units, are localized within large domains of open chromatin, and are significantly associated with DNase I hypersensitivity. Origin "density" falls from early- to mid-S-phase, but rises again in late S-phase to levels only 17% lower than in early S-phase. Unexpectedly, late origin density calculated on the 1-Mb scale increases as a function of increasing chromatin compaction. Furthermore, the median efficiency of origins in late-replicating, heterochromatic domains is only 25% lower than in early-replicating euchromatic loci. Thus, the activation of early- and late-firing origins must be regulated by quintessentially different mechanisms. The aggregate data can be unified into a model in which initiation site selection is driven almost entirely by epigenetic factors that fashion both the long-range and local chromatin environments, with underlying DNA sequence and local transcriptional activity playing only minor roles. Importantly, the comprehensive origin map we have prepared for GM06990 overlaps moderately well with origin maps recently reported for the genomes of four different human cell lines based on the distributions of small nascent strands.

Alterations in replication timing of cancer-related genes in malignant human breast cancer cells

The replication timing of nine genes commonly involved in cancer was investigated in the MCF10 cell lines for human breast cancer progression. Six of these nine genes are part of a constellation of tumor suppressor genes that play a major role in familial human breast cancer (TP53, ATM, PTEN, CHK2, BRCA1, and BRCA2). Three other genes are involved in a large number of human cancers including breast as either tumor suppressors (RB1 and RAD51) or as an oncogene (cMYC). Five of these nine genes (TP53, RAD51, ATM, PTEN, and cMYC) show significant differences (P < 0.05) in replication timing between MCF10A normal human breast cells and the corresponding malignant MCF10CA1a cells. These differences are specific to the malignant state of the MCF10CA1a cells since there were no significant differences in the replication timing of these genes between normal MCF10A cells and the non-malignant cancer MCF10AT1 cells. Microarray analysis further demonstrated that three of these five genes (TP53, RAD51, and cMYC) showed significant changes in gene expression (≥2-fold) between normal and malignant cells. Our findings demonstrate an alteration in the replication timing of a small subset of cancer-related genes in malignant breast cancer cells. These alterations partially correlate with the major transcriptional changes characteristic of the malignant state in these cells.

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

DNA replication timing and higher-order nuclear organization determine single-nucleotide substitution patterns in cancer genomes

Single-nucleotide substitutions are a defining characteristic of cancer genomes. Many single-nucleotide substitutions in cancer genomes arise because of errors in DNA replication, which is spatio-temporally stratified. Here we propose that DNA replication patterns help shape the mutational landscapes of normal and cancer genomes. Using data on five fully sequenced cancer types and two personal genomes, we determined that the frequency of intergenic single-nucleotide substitution is significantly higher in late DNA replication timing regions, even after controlling for a number of genomic features. Furthermore, some substitution signatures are more frequent in certain DNA replication timing zones. Finally, integrating data on higher-order nuclear organization, we found that genomic regions in close spatial proximity to late-replicating domains display similar mutation spectra as the late-replicating regions themselves. These data suggest that DNA replication timing together with higher-order genomic organization contribute to the patterns of single-nucleotide substitution in normal and cancer genomes.

Multiscale analysis of genome-wide replication timing profiles using a wavelet-based signal-processing algorithm

In this protocol, we describe the use of the LastWave open-source signal-processing command language (http://perso.ens-lyon.fr/benjamin.audit/LastWave/) for analyzing cellular DNA replication timing profiles. LastWave makes use of a multiscale, wavelet-based signal-processing algorithm that is based on a rigorous theoretical analysis linking timing profiles to fundamental features of the cell's DNA replication program, such as the average replication fork polarity and the difference between replication origin density and termination site density. We describe the flow of signal-processing operations to obtain interactive visual analyses of DNA replication timing profiles. We focus on procedures for exploring the space-scale map of apparent replication speeds to detect peaks in the replication timing profiles that represent preferential replication initiation zones, and for delimiting U-shaped domains in the replication timing profile. In comparison with the generally adopted approach that involves genome segmentation into regions of constant timing separated by timing transition regions, the present protocol enables the recognition of more complex patterns of the spatio-temporal replication program and has a broader range of applications. Completing the full procedure should not take more than 1 h, although learning the basics of the program can take a few hours and achieving full proficiency in the use of the software may take days.

Cdc6: a multi-functional molecular switch with critical role in carcinogenesis

Research in the last decade revealed an additional role for the Replication Licensing Factor Cdc6 in transcriptional regulation. This novel function has been linked to human cancer development. Here, we summarize all the findings arguing over a role of Cdc6 as a transcriptional repressor and shed light toward new research directions for this field.

Initiation of DNA Replication in the Human Genome

Replication of the human genome relies on the presence of thousands of origins distributed along each of the chromosomes. The activation of these origins occurs in a highly regulated manner to ensure that chromosomes are faithfully duplicated only once during each cell cycle. Failure in this regulation can lead to abnormal cell proliferation, or/and genomic instability, the hallmarks of cancer cells. The mechanisms determining how, when, and where origins are activated remains still a mystery. However recent technological advances have facilitated the study of DNA replication in a genome-wide scale, and have provided a wealth of information on several features of this process. Here we present an overview of the current progress on our understanding of the initiation step of DNA replication in human cells, and its relationship to abnormal cell proliferation.

Replication origins and timing of temporal replication in budding yeast: how to solve the conundrum?

Similarly to metazoans, the budding yeast Saccharomyces cereviasiae replicates its genome with a defined timing. In this organism, well-defined, site-specific origins, are efficient and fire in almost every round of DNA replication. However, this strategy is neither conserved in the fission yeast Saccharomyces pombe, nor in Xenopus or Drosophila embryos, nor in higher eukaryotes, in which DNA replication initiates asynchronously throughout S phase at random sites. Temporal and spatial controls can contribute to the timing of replication such as Cdk activity, origin localization, epigenetic status or gene expression. However, a debate is going on to answer the question how individual origins are selected to fire in budding yeast. Two opposing theories were proposed: the "replicon paradigm" or "temporal program" vs. the "stochastic firing". Recent data support the temporal regulation of origin activation, clustering origins into temporal blocks of early and late replication. Contrarily, strong evidences suggest that stochastic processes acting on origins can generate the observed kinetics of replication without requiring a temporal order. In mammalian cells, a spatiotemporal model that accounts for a partially deterministic and partially stochastic order of DNA replication has been proposed. Is this strategy the solution to reconcile the conundrum of having both organized replication timing and stochastic origin firing also for budding yeast? In this review we discuss this possibility in the light of our recent study on the origin activation, suggesting that there might be a stochastic component in the temporal activation of the replication origins, especially under perturbed conditions.

The effect of the intra-S-phase checkpoint on origins of replication in human cells

Although many chemotherapy drugs activate the intra-S-phase checkpoint pathway to block S-phase progression, not much is known about how and where the intra-S-phase checkpoint regulates origins of replication in human chromosomes. A genomic analysis of replication in human cells in the presence of hydroxyurea (HU) revealed that only the earliest origins fire, but the forks stall within 2 kb and neighboring clusters of dormant origins are activated. The initiation events are located near expressed genes with a preference for transcription start and end sites, and when they are located in intergenic regions they are located near regulatory factor-binding regions (RFBR). The activation of clustered neo-origins by HU suggests that there are many potential replication initiation sites in permissive parts of the genome, most of which are not used in a normal S phase. Consistent with this redundancy, we see multiple sites bound to MCM3 (representative of the helicase) in the region flanking three out of three origins studied in detail. Bypass of the intra-S-phase checkpoint by caffeine activates many new origins in mid- and late-replicating parts of the genome. The intra-S-phase checkpoint suppresses origin firing after the loading of Mcm10, but before the recruitment of Cdc45 and AND-1/CTF4; i.e., after helicase loading but before helicase activation and polymerase loading. Interestingly, Cdc45 recruitment upon checkpoint bypass was accompanied by the restoration of global Cdk2 kinase activity and decrease in both global and origin-bound histone H3 Lys 4 trimethylation (H3K4me3), consistent with the suggestion that both of these factors are important for Cdc45 recruitment.

Replication timing-related and gene body-specific methylation of active human genes

Understanding how the epigenetic blueprint of the genome shapes human phenotypes requires systematic evaluation of the complex interplay between gene activity and the different layers of the epigenome. Utilizing microarray-based techniques, we explored the relationships between DNA methylation, DNA replication timing and gene expression levels across a variety of human tissues and cell lines. The analyses revealed unequal methylation levels among early- and late-replicating fractions of the genome: late-replicating DNA was hypomethylated compared with early-replicating DNA. Moreover, late-replicating regions were gradually demethylated with cell divisions, whereas the methylation of early-replicating regions was better maintained. As active genes concentrate at early-replicating regions, they are overall hypermethylated relative to inactive genes. Accordingly, we show that the previously reported positive correlation between gene-body methylation (methylation of the transcribed portion of genes) and gene expression is restricted to proliferative tissues and cell lines, whereas in tissues containing few proliferating cells, active and inactive genes have similar methylation levels. We further show that active gene bodies are hypermethylated not only compared with inactive gene bodies, but also compared with their flanking sequences. This specific hypermethylation of the active gene bodies is severely disrupted in cells of an immunodeficiency, centromeric region instability, facial anomalies (ICF) syndrome patient bearing mutated DNA methyltransferase 3B (DNMT3B). Our data show that a high methylation level is preferentially maintained in active gene bodies through independent cellular processes. Rather than serving as a distinctive mark between active and inactive genes, gene-body methylation appears to serve a vital, currently unknown function in active genes.

Developmental control of the DNA replication and transcription programs

Polyploid or polytene cells, which have more than 2C DNA content, are widespread throughout nature and present in most differentiated Drosophila tissues. These cells also can display differential replication, that is, genomic regions of increased or decreased DNA copy number relative to overall genomic ploidy. How frequently differential replication is used as a developmental strategy remains unclear. Here, we use genome-wide array-based comparative genomic hybridization (aCGH) to profile differential DNA replication in isolated and purified larval fat body and midgut tissues of Drosophila, and we compare them with recent aCGH profiles of the larval salivary gland. We identify sites of euchromatic underreplication that are common to all three tissues and others that are tissue specific. We demonstrate that both common and tissue-specific underreplicated sites are dependent on the Suppressor of Underreplication protein, SUUR. mRNA-seq profiling shows that whereas underreplicated regions are generally transcriptionally silent in the larval midgut and salivary gland, transcriptional silencing and underreplication have been uncoupled in the larval fat body. In addition to revealing the prevalence of differential replication, our results show that transcriptional silencing and underreplication can be mechanistically uncoupled.

Bubble-chip analysis of human origin distributions demonstrates on a genomic scale significant clustering into zones and significant association with transcription

We have used a novel bubble-trapping procedure to construct nearly pure and comprehensive human origin libraries from early S- and log-phase HeLa cells, and from log-phase GM06990, a karyotypically normal lymphoblastoid cell line. When hybridized to ENCODE tiling arrays, these libraries illuminated 15.3%, 16.4%, and 21.8% of the genome in the ENCODE regions, respectively. Approximately half of the origin fragments cluster into zones, and their signals are generally higher than those of isolated fragments. Interestingly, initiation events are distributed about equally between genic and intergenic template sequences. While only 13.2% and 14.0% of genes within the ENCODE regions are actually transcribed in HeLa and GM06990 cells, 54.5% and 25.6% of zonal origin fragments overlap transcribed genes, most with activating chromatin marks in their promoters. Our data suggest that cell synchronization activates a significant number of inchoate origins. In addition, HeLa and GM06990 cells activate remarkably different origin populations. Finally, there is only moderate concordance between the log-phase HeLa bubble map and published maps of small nascent strands for this cell line.

Integrating one-dimensional and three-dimensional maps of genomes

Genomes exist in vivo as complex physical structures, and their functional output (i.e. the gene expression profile of a cell) is related to their spatial organization inside the nucleus as well as to local chromatin status. Chromatin modifications and chromosome conformation are distinct in different tissues and cell types, which corresponds closely with the diversity in gene-expression patterns found in different tissues of the body. The biological processes and mechanisms driving these general correlations are currently the topic of intense study. An emerging theme is that genome compartmentalization - both along the linear length of chromosomes, and in three dimensions by the spatial colocalization of chromatin domains and genomic loci from across the genome - is a crucial parameter in regulating genome expression. In this Commentary, we propose that a full understanding of genome regulation requires integrating three different types of data: first, one-dimensional data regarding the state of local chromatin - such as patterns of protein binding along chromosomes; second, three-dimensional data that describe the population-averaged folding of chromatin inside cells and; third, single-cell observations of three-dimensional spatial colocalization of genetic loci and trans factors that reveal information about their dynamics and frequency of colocalization.

Comparative analysis of DNA replication timing reveals conserved large-scale chromosomal architecture

Recent evidence suggests that the timing of DNA replication is coordinated across megabase-scale domains in metazoan genomes, yet the importance of this aspect of genome organization is unclear. Here we show that replication timing is remarkably conserved between human and mouse, uncovering large regions that may have been governed by similar replication dynamics since these species have diverged. This conservation is both tissue-specific and independent of the genomic G+C content conservation. Moreover, we show that time of replication is globally conserved despite numerous large-scale genome rearrangements. We systematically identify rearrangement fusion points and demonstrate that replication time can be locally diverged at these loci. Conversely, rearrangements are shown to be correlated with early replication and physical chromosomal proximity. These results suggest that large chromosomal domains of coordinated replication are shuffled by evolution while conserving the large-scale nuclear architecture of the genome.

During S-phase of the cell cycle, chromosomal DNA is replicated in a complex process involving the coordinated activity of thousands of replication forks, each of which duplicates a long stretch of DNA. Recent experiments revealed that the genome is replicating as a mosaic of large-scale early and late chromosomal domains and that this high-level domain organization is correlated with genomic properties like gene density and nucleotide composition. We compared genome-wide replication time maps of compatible human and mouse cells and revealed that their organization into replication domains is highly conserved despite the numerous large-scale genome rearrangements separating the two species. Analysis of recent chromosomal interaction data shows that regions with similar time of replication are more frequently interacting with each other than expected. The data also show that evolutionary rearrangements have predominantly occurred between regions that have similar time of replication and higher-than-expected chromosomal proximity. Our data suggests that the genome, while being continuously rearranged by evolution, maintains a conserved domain organization. Whether this conservation is driven by selection, or is a consequence of the rearrangement process itself, can be resolved by enhancing the comparative approach proposed here.

Arabidopsis thaliana chromosome 4 replicates in two phases that correlate with chromatin state

DNA replication programs have been studied extensively in yeast and animal systems, where they have been shown to correlate with gene expression and certain epigenetic modifications. Despite the conservation of core DNA replication proteins, little is known about replication programs in plants. We used flow cytometry and tiling microarrays to profile DNA replication of Arabidopsis thaliana chromosome 4 (chr4) during early, mid, and late S phase. Replication profiles for early and mid S phase were similar and encompassed the majority of the euchromatin. Late S phase exhibited a distinctly different profile that includes the remaining euchromatin and essentially all of the heterochromatin. Termination zones were consistent between experiments, allowing us to define 163 putative replicons on chr4 that clustered into larger domains of predominately early or late replication. Early-replicating sequences, especially the initiation zones of early replicons, displayed a pattern of epigenetic modifications specifying an open chromatin conformation. Late replicons, and the termination zones of early replicons, showed an opposite pattern. Histone H3 acetylated on lysine 56 (H3K56ac) was enriched in early replicons, as well as the initiation zones of both early and late replicons. H3K56ac was also associated with expressed genes, but this effect was local whereas replication time correlated with H3K56ac over broad regions. The similarity of the replication profiles for early and mid S phase cells indicates that replication origin activation in euchromatin is stochastic. Replicon organization in Arabidopsis is strongly influenced by epigenetic modifications to histones and DNA. The domain organization of Arabidopsis is more similar to that in Drosophila than that in mammals, which may reflect genome size and complexity. The distinct patterns of association of H3K56ac with gene expression and early replication provide evidence that H3K56ac may be associated with initiation zones and replication origins.

During growth and development, all plants and animals must replicate their DNA. This process is regulated to ensure that all sequences are completely and accurately replicated and is limited to S phase of the cell cycle. In the cell, DNA is packaged with histone proteins into chromatin, and both DNA and histones are subject to epigenetic modifications that affect chromatin state. Euchromatin and heterochromatin are chromatin states marked by epigenetic modifications specifying open and closed conformations, respectively. Using the model plant Arabidopsis thaliana, we show that the time at which a DNA sequence replicates is influenced by the epigenetic modifications to the surrounding chromatin. DNA replication occurs in two phases, with euchromatin replicating in early and mid S phase and heterochromatin replicating late. DNA replication time has been linked to gene expression in other organisms, and this is also true in Arabidopsis because more genes are active in euchromatin when compared to heterochromatin. The earliest replicating DNA sequences are associated with acetylation of histone H3 on lysine 56 (H3K56ac). H3K56ac is also abundant in active genes, but the patterns of association of H3K56ac with gene expression and DNA replication are distinct, suggesting that H3K56ac is independently linked to both processes.

Genome-wide analysis of the replication program in mammals

Microarray technology has facilitated the research of eukaryotic DNA replication on a genome-wide scale. Recent studies have shed light on the association between time of replication and chromosome structure, on the organization principles of the replication program, and on the correlation between replication timing and transcription. In this review, we summarize various genomic measurement approaches and the biological insights achieved through applying them in the study of the mammalian replication program.

S-phase progression in mammalian cells: modelling the influence of nuclear organization

The control of DNA replication is of fundamental importance as cell proliferation demands that identical copies of the genetic material are passed to the two daughter cells that form during mitosis. These genetic copies are generated in the preceding S phase, where the entire DNA complement of the mother cell must be copied exactly once. As part of this process, it is known that different regions of mammalian genomes are replicated at specific times of a temporally defined replication programme. The key feature of this programme is that active genes in euchromatin are replicated before inactive ones in heterochromatin. This separation of S phase into periods where different classes of chromatin are duplicated is important in maintaining changes in gene expression that define individual cell types. Recent attempts to understand the structure of the S-phase timing programme have focused on the use of genome-wide strategies that inevitably use DNA isolated from large cell populations for analysis. However, this approach provides a composite view of events that occur within a population without knowledge of the cell-to-cell variability across the population. In this review, we attempt to combine information generated using genome-wide and single cell strategies in order to develop a coherent molecular understanding of S-phase progression. During this integration, we have explored how available information can be introduced into a modelling environment that best describes S-phase progression in mammalian cells.

Chk1-cyclin A/Cdk1 axis regulates origin firing programs in mammals

DNA replication is key to ensuring the complete duplication of genomic DNA prior to mitosis and is tightly regulated by both cell cycle machinery and checkpoint signals. Regulation of the S phase program occurs at several stages, affecting origin firing, replication fork elongation, fork velocity, and fork stability, all of which are dependent on S-phase-promoting kinase activity. Somatic mammalian cells use well-established origin programs by which specific regions of the genome are replicated at precise times. However, the mechanisms by which S phase kinases regulate origin firing in mammals are largely unknown. Here, we discuss recent advances in the understanding of how S phase programs are regulated in mammals at the correct regions and at the appropriate times.

Genome-wide estimation of firing efficiencies of origins of DNA replication from time-course copy number variation data

Background

DNA replication is a fundamental biological process during S phase of cell division. It is initiated from several hundreds of origins along whole chromosome with different firing efficiencies (or frequency of usage). Direct measurement of origin firing efficiency by techniques such as DNA combing are time-consuming and lack the ability to measure all origins. Recent genome-wide study of DNA replication approximated origin firing efficiency by indirectly measuring other quantities related to replication. However, these approximation methods do not reflect properties of origin firing and may lead to inappropriate estimations.

Results

In this paper, we develop a probabilistic model - Spanned Firing Time Model (SFTM) to characterize DNA replication process. The proposed model reflects current understandings about DNA replication. Origins in an individual cell may initiate replication randomly within a time window, but the population average exhibits a temporal program with some origins replicated early and the others late. By estimating DNA origin firing time and fork moving velocity from genome-wide time-course S-phase copy number variation data, we could estimate firing efficiency of all origins. The estimated firing efficiency is correlated well with the previous studies in fission and budding yeasts.

Conclusions

The new probabilistic model enables sensitive identification of origins as well as genome-wide estimation of origin firing efficiency. We have successfully estimated firing efficiencies of all origins in S.cerevisiae, S.pombe and human chromosomes 21 and 22.

Organization of DNA replication

The discovery of the DNA double helix structure half a century ago immediately suggested a mechanism for its duplication by semi-conservative copying of the nucleotide sequence into two DNA daughter strands. Shortly after, a second fundamental step toward the elucidation of the mechanism of DNA replication was taken with the isolation of the first enzyme able to polymerize DNA from a template. In the subsequent years, the basic mechanism of DNA replication and its enzymatic machinery components were elucidated, mostly through genetic approaches and in vitro biochemistry. Most recently, the spatial and temporal organization of the DNA replication process in vivo within the context of chromatin and inside the intact cell are finally beginning to be elucidated. On the one hand, recent advances in genome-wide high throughput techniques are providing a new wave of information on the progression of genome replication at high spatial resolution. On the other hand, novel super-resolution microscopy techniques are just starting to give us the first glimpses of how DNA replication is organized within the context of single intact cells with high spatial resolution. The integration of these data with time lapse microscopy analysis will give us the ability to film and dissect the replication of the genome in situ and in real time.

Impact of replication timing on non-CpG and CpG substitution rates in mammalian genomes

Neutral nucleotide substitutions occur at varying rates along genomes, and it remains a major issue to unravel the mechanisms that cause these variations and to analyze their evolutionary consequences. Here, we study the role of replication in the neutral substitution pattern. We obtained a high-resolution replication timing profile of the whole human genome by massively parallel sequencing of nascent BrdU-labeled replicating DNA. These data were compared to the neutral substitution rates along the human genome, obtained by aligning human and chimpanzee genomes using macaque and orangutan as outgroups. All substitution rates increase monotonously with replication timing even after controlling for local or regional nucleotide composition, crossover rate, distance to telomeres, and chromatin compaction. The increase in non-CpG substitution rates might result from several mechanisms including the increase in mutation-prone activities or the decrease in efficiency of DNA repair during the S phase. In contrast, the rate of C --> T transitions in CpG dinucleotides increases in later-replicating regions due to increasing DNA methylation level that reflects a negative correlation between timing and gene expression. Similar results are observed in the mouse, which indicates that replication timing is a main factor affecting nucleotide substitution dynamics at non-CpG sites and constitutes a major neutral process driving mammalian genome evolution.

Replication factory activation can be decoupled from the replication timing program by modulating Cdk levels

Cdk activity can differentially regulate the number of active replication factories, replication rates, and the rate of progression through the timing program during S phase.

In the metazoan replication timing program, clusters of replication origins located in different subchromosomal domains fire at different times during S phase. We have used Xenopus laevis egg extracts to drive an accelerated replication timing program in mammalian nuclei. Although replicative stress caused checkpoint-induced slowing of the timing program, inhibition of checkpoint kinases in an unperturbed S phase did not accelerate it. Lowering cyclin-dependent kinase (Cdk) activity slowed both replication rate and progression through the timing program, whereas raising Cdk activity increased them. Surprisingly, modest alteration of Cdk activity changed the amount of DNA synthesized during different stages of the timing program. This was associated with a change in the number of active replication factories, whereas the distribution of origins within active factories remained relatively normal. The ability of Cdks to differentially effect replication initiation, factory activation, and progression through the timing program provides new insights into the way that chromosomal DNA replication is organized during S phase.

Temporal differences in DNA replication during the S phase using single fiber analysis of normal human fibroblasts and glioblastoma T98G cells

We have recently shown that replication forks pause near origins in normal human fibroblasts (NHF1-hTERT) but not glioblastoma T98G cells. This observation led us to question whether other differences in the replication program may exist between these cell types that may relate to their genetic integrity. To identify differences, we detected immunoflourescently the sequential incorporation of the nucleotide analogs IdU and CldU into replicating DNA at the start of every hour of a synchronized S phase. We then characterized the patterns of labeled replicating DNA tracks and quantified the percentages and lengths of the tracks found at these hourly intervals. From the directionality of labeling in single extended replicating DNA fibers, tracks were categorized as single bidirectional origins, unidirectional elongations, clusters of origins firing in tandem, or merging forks (terminations). Our analysis showed that the start of S phase is enriched in single bidirectional origins in NHF1-hTERT cells, followed by an increase in clustering during mid S phase and an increase in merging forks during late S phase. Early S phase in T98G cells also largely consisted of single bidirectional origin initiations; however, an increase in clustering was delayed until an hour later, and clusters were shorter in mid/late S phase than in NHF1-hTERT cells. The spike in merging forks also did not occur until an hour later in T98G cells. Our observations suggest models to explain the temporal replication of single and clustered origins, and suggest differences in the replication program in a normal and cancer cell line.

Sequencing newly replicated DNA reveals widespread plasticity in human replication timing

Faithful transmission of genetic material to daughter cells involves a characteristic temporal order of DNA replication, which may play a significant role in the inheritance of epigenetic states. We developed a genome-scale approach--Repli Seq--to map temporally ordered replicating DNA using massively parallel sequencing and applied it to study regional variation in human DNA replication time across multiple human cell types. The method requires as few as 8,000 cytometry-fractionated cells for a single analysis, and provides high-resolution DNA replication patterns with respect to both cell-cycle time and genomic position. We find that different cell types exhibit characteristic replication signatures that reveal striking plasticity in regional replication time patterns covering at least 50% of the human genome. We also identified autosomal regions with marked biphasic replication timing that include known regions of monoallelic expression as well as many previously uncharacterized domains. Comparison with high-resolution genome-wide profiles of DNaseI sensitivity revealed that DNA replication typically initiates within foci of accessible chromatin comprising clustered DNaseI hypersensitive sites, and that replication time is better correlated with chromatin accessibility than with gene expression. The data collectively provide a unique, genome-wide picture of the epigenetic compartmentalization of the human genome and suggest that cell-lineage specification involves extensive reprogramming of replication timing patterns.

Dynamic changes in chromatin structure through post-translational modifications of histone H3 during replication origin activation

Genome duplication relies on the timely activation of multiple replication origins throughout the genome during S phase. Each origin is marked by the assembly of a multiprotein pre-replication complex (pre-RC) and the recruitment of the replicative machinery, which can gain access to replication origins on the DNA through the barrier of specific chromatin structures. Inheritance of the genetic information is further accompanied by maintenance and inheritance of the epigenetic marks, which are accomplished by the activity of histone and DNA modifying enzymes traveling with the replisome. Here, we studied the changes in the chromatin structure at the loci of three replication origins, the early activated human lamin B2 (LB2) and monkey Ors8 (mOrs8) origins and the late-activated human homologue of the latter (hOrs8), during their activation, by measuring the abundance of post-translationally modified histone H3. The data show that dynamic changes in the levels of acetylated, methylated and phosphorylated histone H3 occur during the initiation of DNA replication at these three origin loci, which differ between early- and late-firing origins as well as between human- and monkey-derived cell lines. These results suggest that specific histone modifications are associated with origin firing, temporal activation and replication fork progression and underscore the importance of species specificity.

Strategies for analyzing highly enriched IP-chip datasets

Background

Chromatin immunoprecipitation on tiling arrays (ChIP-chip) has been employed to examine features such as protein binding and histone modifications on a genome-wide scale in a variety of cell types. Array data from the latter studies typically have a high proportion of enriched probes whose signals vary considerably (due to heterogeneity in the cell population), and this makes their normalization and downstream analysis difficult.

Results

Here we present strategies for analyzing such experiments, focusing our discussion on the analysis of Bromodeoxyruridine (BrdU) immunoprecipitation on tiling array (BrdU-IP-chip) datasets. BrdU-IP-chip experiments map large, recently replicated genomic regions and have similar characteristics to histone modification/location data. To prepare such data for downstream analysis we employ a dynamic programming algorithm that identifies a set of putative unenriched probes, which we use for both within-array and between-array normalization. We also introduce a second dynamic programming algorithm that incorporates a priori knowledge to identify and quantify positive signals in these datasets.

Conclusion

Highly enriched IP-chip datasets are often difficult to analyze with traditional array normalization and analysis strategies. Here we present and test a set of analytical tools for their normalization and quantification that allows for accurate identification and analysis of enriched regions.

Predictable dynamic program of timing of DNA replication in human cells

The organization of mammalian DNA replication is poorly understood. We have produced high-resolution dynamic maps of the timing of replication in human erythroid, mesenchymal, and embryonic stem (ES) cells using TimEX, a method that relies on gaussian convolution of massive, highly redundant determinations of DNA copy-number variations during S phase to produce replication timing profiles. We first obtained timing maps of 3% of the genome using high-density oligonucleotide tiling arrays and then extended the TimEX method genome-wide using massively parallel sequencing. We show that in untransformed human cells, timing of replication is highly regulated and highly synchronous, and that many genomic segments are replicated in temporal transition regions devoid of initiation, where replication forks progress unidirectionally from origins that can be hundreds of kilobases away. Absence of initiation in one transition region is shown at the molecular level by single molecule analysis of replicated DNA (SMARD). Comparison of ES and erythroid cells replication patterns revealed that these cells replicate about 20% of their genome in different quarters of S phase. Importantly, we detected a strong inverse relationship between timing of replication and distance to the closest expressed gene. This relationship can be used to predict tissue-specific timing of replication profiles from expression data and genomic annotations. We also provide evidence that early origins of replication are preferentially located near highly expressed genes, that mid-firing origins are located near moderately expressed genes, and that late-firing origins are located far from genes.

Shifts in replication timing actively affect histone acetylation during nucleosome reassembly

The entire genome is replicated in a programmed manner, with specific regions undergoing DNA synthesis at different times in S phase. Active genes generally replicate in early S phase, while repressed genes replicate late, and for some loci this process is developmentally regulated. Using a nuclear microinjection system, we demonstrate that DNA sequences originally packaged into nucleosomes containing deacetylated histones during late S become reassembled with acetylated histones after undergoing replication in early S. Conversely, a change from early to late replication timing is accompanied by repackaging into nucleosomes containing deacetylated histones. This is carried out by differential cell-cycle-controlled acetylation and deacetylation of histones H3 and H4. These studies provide strong evidence that switches in replication timing may play a role in the regulation of nucleosome structure during development.

Metazoan origins of DNA replication: regulation through dynamic chromatin structure

DNA replication in eukaryotes is initiated at multiple replication origins distributed over the entire genome, which are normally activated once per cell cycle. Due to the complexity of the metazoan genome, the study of metazoan replication origins and their activity profiles has been less advanced than in simpler genome systems. DNA replication in eukaryotes involves many protein-protein and protein-DNA interactions, occurring in multiple stages. As in prokaryotes, control over the timing and frequency of initiation is exerted at the initiation site. A prerequisite for understanding the regulatory mechanisms of eukaryotic DNA replication is the identification and characterization of the cis-acting sequences that serve as replication origins and the trans-acting factors (proteins) that interact with them. Furthermore, in order to understand how DNA replication may become deregulated in malignant cells, the distinguishing features between normal and malignant origins of DNA replication as well as the proteins that interact with them must be determined. Based on advances that were made using simple genome model systems, several proteins involved in DNA replication have been identified. This review summarizes the current findings about metazoan origins of DNA replication and their interacting proteins as well as the role of chromatin structure in their regulation. Furthermore, progress in origin identification and isolation procedures as well as potential mechanisms to inhibit their activation in cancer development and progression are discussed.

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

Replication timing is frequently discussed superficially in terms of its relationship to transcriptional activity via chromatin structure. However, so little is known about what regulates where and when replication initiates that it has been impossible to identify mechanistic and causal relationships. Moreover, much of our knowledge base has been anecdotal, derived from analyses of a few genes in unrelated cell lines. Recent studies have revisited long-standing hypotheses using genome-wide approaches. In particular, the foundation of this field was recently shored up with incontrovertible evidence that cellular differentiation is accompanied by coordinated changes in replication timing and transcription. These changes accompany subnuclear repositioning, and take place at the level of megabase-sized domains that transcend localized changes in chromatin structure or transcription. Inferring from these results, we propose that there exists a key transition during the middle of S-phase and that changes in replication timing traversing this period are associated with subnuclear repositioning and changes in the activity of certain classes of promoters.

Human mutation rate associated with DNA replication timing

Eukaryotic DNA replication is highly stratified, with different genomic regions shown to replicate at characteristic times during S phase. Here we observe that mutation rate, as reflected in recent evolutionary divergence and human nucleotide diversity, is markedly increased in later-replicating regions of the human genome. All classes of substitutions are affected, suggesting a generalized mechanism involving replication time-dependent DNA damage. This correlation between mutation rate and regionally stratified replication timing may have substantial evolutionary implications.

Cyclin A-Cdk1 regulates the origin firing program in mammalian cells

Somatic mammalian cells possess well-established S-phase programs with specific regions of the genome replicated at precise times. The ATR-Chk1 pathway plays a central role in these programs, but the mechanism for how Chk1 regulates origin firing remains unknown. We demonstrate here the essential role of cyclin A2-Cdk1 in the regulation of late origin firing. Activity of cyclin A2-Cdk1 was hardly detected at the onset of S phase, but it was obvious at middle to late S phase under unperturbed condition. Chk1 depletion resulted in increased expression of Cdc25A, subsequent hyperactivation of cyclin A2-Cdk1, and abnormal replication at early S phase. Hence, the ectopic expression of cyclin A2-Cdk1AF (constitutively active mutant) fusion constructs resulted in abnormal origin firing, causing the premature appearance of DNA replication at late origins at early S phase. Intriguingly, inactivation of Cdk1 in temperature-sensitive Cdk1 mutant cell lines (FT210) resulted in a prolonged S phase and inefficient activation of late origin firing even at late S phase. Our results thus suggest that cyclin A2-Cdk1 is a key regulator of S-phase programs.

Analyzing origin activation patterns by copy number change experiments

Advances in microarray technology have enabled the analysis of replication dynamics on a genome-wide scale, providing deeper insight to the factors that regulate DNA replication. Studies using high-density microarrays have led to the genome-wide identification of replication origins in the budding yeast, Saccharomyces cerevisiae, and enabled the analysis of the global temporal pattern of origin activation under various conditions. We have developed a replication origin array that contains the approximately 430 potential origins in the yeast genome. By detecting the copy number change that occurs as cells progress from G1 to S phase on these arrays, we have produced origin activation patterns in wild-type cells similar to those obtained from previous studies that used whole-genome arrays. We have also applied this method to study S phase checkpoint mutants, providing insight into the genome-wide regulation of replication origin activation by S phase checkpoint kinases in the presence of replication stress. The main procedures of this technique involve arresting yeast cells in G1 and S phase, isolating and labeling genomic DNA with fluorescent dyes, and cohybridizing the DNA samples to replication origin arrays to yield copy number change data.

Introduction to molecular combing: genomics, DNA replication, and cancer

The sequencing of the human genome inaugurated a new era in both fundamental and applied genetics. At the same time, the emergence of new technologies for probing the genome has transformed the field of pharmaco-genetics and made personalized genomic profiling and high-throughput screening of new therapeutic agents all but a matter of routine. One of these technologies, molecular combing, has served to bridge the technical gap between the examination of gross chromosomal abnormalities and sequence-specific alterations. Molecular combing provides a new perspective on the structure and dynamics of the human genome at the whole genome and sub-chromosomal levels with a resolution ranging from a few kilobases up to a megabase and more. Originally developed to study genetic rearrangements and to map genes for positional cloning, recent advances have extended the spectrum of its applications to studying the real-time dynamics of the replication of the genome. Understanding how the genome is replicated is essential for elucidating the mechanisms that both maintain genome integrity and result in the instabilities leading to human genetic disease and cancer. In the following, we will examine recent discoveries and advances due to the application of molecular combing to new areas of research in the fields of molecular cytogenetics and cancer genomics.

Microarray analysis of DNA replication timing

Although all of the DNA in an eukaryotic cell replicates during the S-phase of cell cycle, there is a significant difference in the actual time in S-phase when a given chromosomal segment replicates. Methods are described here for generation of high-resolution temporal maps of DNA replication in synchronized human cells. This method does not require amplification of DNA before microarray hybridization and so avoids errors introduced during PCR. A major advantage of using this procedure is that it facilitates finer dissection of replication time in S-phase. Also, it helps delineate chromosomal regions that undergo biallelic or asynchronous replication, which otherwise are difficult to detect at a genome-wide scale by existing methods. The continuous TR50 (time of completion of 50% replication) maps of replication across chromosomal segments identify regions that undergo acute transitions in replication timing. These transition zones can play a significant role in identifying insulators that separate chromosomal domains with different chromatin modifications.