Dynamics of DNA replication in mammalian somatic cells: nucleotide pool modulates origin choice and interorigin spacing

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.

dNTP pools determine fork progression and origin usage under replication stress

Intracellular deoxyribonucleoside triphosphate (dNTP) pools must be tightly regulated to preserve genome integrity. Indeed, alterations in dNTP pools are associated with increased mutagenesis, genomic instability and tumourigenesis. However, the mechanisms by which altered or imbalanced dNTP pools affect DNA synthesis remain poorly understood. Here, we show that changes in intracellular dNTP levels affect replication dynamics in budding yeast in different ways. Upregulation of the activity of ribonucleotide reductase (RNR) increases elongation, indicating that dNTP pools are limiting for normal DNA replication. In contrast, inhibition of RNR activity with hydroxyurea (HU) induces a sharp transition to a slow-replication mode within minutes after S-phase entry. Upregulation of RNR activity delays this transition and modulates both fork speed and origin usage under replication stress. Interestingly, we also observed that chromosomal instability (CIN) mutants have increased dNTP pools and show enhanced DNA synthesis in the presence of HU. Since upregulation of RNR promotes fork progression in the presence of DNA lesions, we propose that CIN mutants adapt to chronic replication stress by upregulating dNTP pools.

Evidence for sequential and increasing activation of replication origins along replication timing gradients in the human genome

Genome-wide replication timing studies have suggested that mammalian chromosomes consist of megabase-scale domains of coordinated origin firing separated by large originless transition regions. Here, we report a quantitative genome-wide analysis of DNA replication kinetics in several human cell types that contradicts this view. DNA combing in HeLa cells sorted into four temporal compartments of S phase shows that replication origins are spaced at 40 kb intervals and fire as small clusters whose synchrony increases during S phase and that replication fork velocity (mean 0.7 kb/min, maximum 2.0 kb/min) remains constant and narrowly distributed through S phase. However, multi-scale analysis of a genome-wide replication timing profile shows a broad distribution of replication timing gradients with practically no regions larger than 100 kb replicating at less than 2 kb/min. Therefore, HeLa cells lack large regions of unidirectional fork progression. Temporal transition regions are replicated by sequential activation of origins at a rate that increases during S phase and replication timing gradients are set by the delay and the spacing between successive origin firings rather than by the velocity of single forks. Activation of internal origins in a specific temporal transition region is directly demonstrated by DNA combing of the IGH locus in HeLa cells. Analysis of published origin maps in HeLa cells and published replication timing and DNA combing data in several other cell types corroborate these findings, with the interesting exception of embryonic stem cells where regions of unidirectional fork progression seem more abundant. These results can be explained if origins fire independently of each other but under the control of long-range chromatin structure, or if replication forks progressing from early origins stimulate initiation in nearby unreplicated DNA. These findings shed a new light on the replication timing program of mammalian genomes and provide a general model for their replication kinetics.

Eukaryotic chromosomes replicate from multiple replication origins that fire at different times in S phase. The mechanisms that specify origin position and firing time and coordinate origins to ensure complete genome duplication are unclear. Previous studies proposed either that origins are arranged in temporally coordinated groups or fire independently of each other in a stochastic manner. Here, we have performed a quantitative analysis of human genome replication kinetics using a combination of DNA combing, which reveals local patterns of origin firing and replication fork progression on single DNA molecules, and massive sequencing of newly replicated DNA, which reveals the population-averaged replication timing profile of the entire genome. We show that origins are activated synchronously in large regions of uniform replication timing but more gradually in temporal transition regions and that the rate of origin firing increases as replication progresses. Large regions of unidirectional fork progression are abundant in embryonic stem cells but rare in differentiated cells. We propose a model in which replication forks progressing from early origins stimulate initiation in nearby unreplicated DNA in a manner that explains the shape of the replication timing profile. These results provide a fundamental insight into the temporal regulation of mammalian genome replication.

Dysregulation of DNA polymerase κ recruitment to replication forks results in genomic instability

Translesion synthesis polymerases (TLS Pols) are required to tolerate DNA lesions that would otherwise cause replication arrest and cell death. Aberrant expression of these specialized Pols may be responsible for increased mutagenesis and loss of genome integrity in human cancers. The molecular events that control the usage of TLS Pols in non-pathological conditions remain largely unknown. Here, we show that aberrant recruitment of TLS Polκ to replication forks results in genomic instability and can be mediated through the loss of the deubiquitinase USP1. Moreover, artificial tethering of Polκ to proliferating cell nuclear antigen (PCNA) circumvents the need for its ubiquitin-binding domain in the promotion of genomic instability. Finally, we show that the loss of USP1 leads to a dramatic reduction of replication fork speed in a Polκ-dependent manner. We propose a mechanism whereby reversible ubiquitination of PCNA can prevent spurious TLS Pol recruitment and regulate replication fork speed to ensure the maintenance of genome integrity.

Aberrant firing of replication origins potentially explains intragenic nonrecurrent rearrangements within genes, including the human DMD gene

Non-allelic homologous recombination (NAHR), non-homologous end joining (NHEJ), and microhomology-mediated replication-dependent recombination (MMRDR) have all been put forward as mechanisms to explain DNA rearrangements associated with genomic disorders. However, many nonrecurrent rearrangements in humans remain unexplained. To further investigate the mutation mechanisms of these copy number variations (CNVs), we performed breakpoint mapping analysis for 62 clinical cases with intragenic deletions in the human DMD gene (50 cases) and other known disease-causing genes (one PCCB, one IVD, one DBT, three PAH, one STK11, one HEXB, three DBT, one HRPT1, and one EMD cases). While repetitive elements were found in only four individual cases, three involving DMD and one HEXB gene, microhomologies (2-10 bp) were observed at breakpoint junctions in 56% and insertions ranging from 1 to 48 bp were seen in 16 of the total 62 cases. Among these insertions, we observed evidence for tandem repetitions of short segments (5-20 bp) of reference sequence proximal to the breakpoints in six individual DMD cases (six repeats in one, four repeats in three, two repeats in one, and one repeat in one case), strongly indicating attempts by the replication machinery to surpass the stalled replication fork. We provide evidence of a novel template slippage event during replication rescue. With a deeper insight into the complex process of replication and its rescue during origin failure, brought forward by recent studies, we propose a hypothesis based on aberrant firing of replication origins to explain intragenic nonrecurrent rearrangements within genes, including the DMD gene.

On the opportunistic nature of transcription and replication initiation in the metazoan genome

Cellular identity and its response to external or internal signalling variations are encoded in a cell's genome as regulatory information. The genomic regions that specify this type of information are highly variable and degenerated in their sequence determinants, as it is becoming increasingly evident through the application of genome-scale methods to study gene expression. Here, we speculate that the same scenario applies to the regulatory regions controlling where DNA replication starts in the metazoan genome. We propose that replication origins cannot be defined as unique genomic features, but rather that DNA synthesis initiates opportunistically from accessible DNA sites, making cells highly robust and adaptable to environmental or developmental changes.

Wee1 controls genomic stability during replication by regulating the Mus81-Eme1 endonuclease

Wee1 is essential for normal DNA replication and for genomic stability, at least in part by inhibiting a general DNA damage response induced by the Mus81-Eme1 endonuclease.

Correct replication of the genome and protection of its integrity are essential for cell survival. In a high-throughput screen studying H2AX phosphorylation, we identified Wee1 as a regulator of genomic stability. Wee1 down-regulation not only induced H2AX phosphorylation but also triggered a general deoxyribonucleic acid (DNA) damage response (DDR) and caused a block in DNA replication, resulting in accumulation of cells in S phase. Wee1-deficient cells showed a decrease in replication fork speed, demonstrating the involvement of Wee1 in DNA replication. Inhibiting Wee1 in cells treated with short treatment of hydroxyurea enhanced the DDR, which suggests that Wee1 specifically protects the stability of stalled replication forks. Notably, the DDR induced by depletion of Wee1 critically depends on the Mus81-Eme1 endonuclease, and we found that codepletion of Mus81 and Wee1 abrogated the S phase delay. Importantly, Wee1 and Mus81 interact in vivo, suggesting direct regulation. Altogether, these results demonstrate a novel role of Wee1 in controlling Mus81 and DNA replication in human cells.

Hydroxyurea induces de novo copy number variants in human cells

Copy number variants (CNVs) are widely distributed throughout the human genome, where they contribute to genetic variation and phenotypic diversity. Spontaneous CNVs are also a major cause of genetic and developmental disorders and arise frequently in cancer cells. As with all mutation classes, genetic and environmental factors almost certainly increase the risk for new and deleterious CNVs. However, despite the importance of CNVs, there is limited understanding of these precipitating risk factors and the mechanisms responsible for a large percentage of CNVs. Here we report that low doses of hydroxyurea, an inhibitor of ribonucleotide reductase and an important drug in the treatment of sickle cell disease and other diseases induces a high frequency of de novo CNVs in cultured human cells that resemble pathogenic and aphidicolin-induced CNVs in size and breakpoint structure. These CNVs are distributed throughout the genome, with some hotspots of de novo CNV formation. Sequencing revealed that CNV breakpoint junctions are characterized by short microhomologies, blunt ends, and short insertions. These data provide direct experimental support for models of replication-error origins of CNVs and suggest that any agent or condition that leads to replication stress has the potential to induce deleterious CNVs. In addition, they point to a need for further study of the genomic consequences of the therapeutic use of hydroxyurea.

Genome-scale analysis of metazoan replication origins reveals their organization in specific but flexible sites defined by conserved features

In metazoans, thousands of DNA replication origins (Oris) are activated at each cell cycle. Their genomic organization and their genetic nature remain elusive. Here, we characterized Oris by nascent strand (NS) purification and a genome-wide analysis in Drosophila and mouse cells. We show that in both species most CpG islands (CGI) contain Oris, although methylation is nearly absent in Drosophila, indicating that this epigenetic mark is not crucial for defining the activated origin. Initiation of DNA synthesis starts at the borders of CGI, resulting in a striking bimodal distribution of NS, suggestive of a dual initiation event. Oris contain a unique nucleotide skew around NS peaks, characterized by G/T and C/A overrepresentation at the 5' and 3' of Ori sites, respectively. Repeated GC-rich elements were detected, which are good predictors of Oris, suggesting that common sequence features are part of metazoan Oris. In the heterochromatic chromosome 4 of Drosophila, Oris correlated with HP1 binding sites. At the chromosome level, regions rich in Oris are early replicating, whereas Ori-poor regions are late replicating. The genome-wide analysis was coupled with a DNA combing analysis to unravel the organization of Oris. The results indicate that Oris are in a large excess, but their activation does not occur at random. They are organized in groups of site-specific but flexible origins that define replicons, where a single origin is activated in each replicon. This organization provides both site specificity and Ori firing flexibility in each replicon, allowing possible adaptation to environmental cues and cell fates.

Nucleotide deficiency promotes genomic instability in early stages of cancer development

Chromosomal instability in early cancer stages is caused by stress on DNA replication. The molecular basis for replication perturbation in this context is currently unknown. We studied the replication dynamics in cells in which a regulator of S phase entry and cell proliferation, the Rb-E2F pathway, is aberrantly activated. Aberrant activation of this pathway by HPV-16 E6/E7 or cyclin E oncogenes significantly decreased the cellular nucleotide levels in the newly transformed cells. Exogenously supplied nucleosides rescued the replication stress and DNA damage and dramatically decreased oncogene-induced transformation. Increased transcription of nucleotide biosynthesis genes, mediated by expressing the transcription factor c-myc, increased the nucleotide pool and also rescued the replication-induced DNA damage. Our results suggest a model for early oncogenesis in which uncoordinated activation of factors regulating cell proliferation leads to insufficient nucleotides that fail to support normal replication and genome stability.

How dormant origins promote complete genome replication

Many replication origins that are licensed by loading MCM2-7 complexes in G1 are not normally used. Activation of these dormant origins during S phase provides a first line of defence for the genome if replication is inhibited. When replication forks fail, dormant origins are activated within regions of the genome currently engaged in replication. At the same time, DNA damage-response kinases activated by the stalled forks preferentially suppress the assembly of new replication factories, thereby ensuring that chromosomal regions experiencing replicative stress complete synthesis before new regions of the genome are replicated. Mice expressing reduced levels of MCM2-7 have fewer dormant origins, are cancer-prone and are genetically unstable, demonstrating the importance of dormant origins for preserving genome integrity. We review the function of dormant origins, the molecular mechanism of their regulation and their physiological implications.

Genome-scale analysis of replication timing: from bench to bioinformatics

Replication timing profiles are cell type-specific and reflect genome organization changes during differentiation. In this protocol, we describe how to analyze genome-wide replication timing (RT) in mammalian cells. Asynchronously cycling cells are pulse labeled with the nucleotide analog 5-bromo-2-deoxyuridine (BrdU) and sorted into S-phase fractions on the basis of DNA content using flow cytometry. BrdU-labeled DNA from each fraction is immunoprecipitated, amplified, differentially labeled and co-hybridized to a whole-genome comparative genomic hybridization microarray, which is currently more cost effective than high-throughput sequencing and equally capable of resolving features at the biologically relevant level of tens to hundreds of kilobases. We also present a guide to analyzing the resulting data sets based on methods we use routinely. Subjects include normalization, scaling and data quality measures, LOESS (local polynomial) smoothing of RT values, segmentation of data into domains and assignment of timing values to gene promoters. Finally, we cover clustering methods and means to relate changes in the replication program to gene expression and other genetic and epigenetic data sets. Some experience with R or similar programming languages is assumed. All together, the protocol takes ∼3 weeks per batch of samples.

Open sesame: activating dormant replication origins in the mouse immunoglobulin heavy chain (Igh) locus

Chromosomal DNA replication in mammals initiates from replication origins whose activity differs in accordance with cell type and differentiation state. In addition to origins that are active in unperturbed conditions, chromosomes also contain dormant origins that can become functional in response to certain genotoxic stress conditions. Improper regulation of origin usage can cause genomic instability leading to tumorigenesis. We review findings from recent single-molecule DNA fiber studies examining replication of the mouse immunoglobulin heavy chain (Igh) locus, in which origin activity over a 400kb region is subject to dramatic developmental regulation. Possible models are discussed to explain such differential origin usage, particularly during replication stress conditions that can activate dormant origins.

Preferential localization of human origins of DNA replication at the 5'-ends of expressed genes and at evolutionarily conserved DNA sequences

BACKGROUND

Replication of mammalian genomes requires the activation of thousands of origins which are both spatially and temporally regulated by as yet unknown mechanisms. At the most fundamental level, our knowledge about the distribution pattern of origins in each of the chromosomes, among different cell types, and whether the physiological state of the cells alters this distribution is at present very limited.

METHODOLOGY/PRINCIPAL FINDINGS

We have used standard λ-exonuclease resistant nascent DNA preparations in the size range of 0.7-1.5 kb obtained from the breast cancer cell line MCF-7 hybridized to a custom tiling array containing 50-60 nt probes evenly distributed among genic and non-genic regions covering about 1% of the human genome. A similar DNA preparation was used for high-throughput DNA sequencing. Array experiments were also performed with DNA obtained from BT-474 and H520 cell lines. By determining the sites showing nascent DNA enrichment, we have localized several thousand origins of DNA replication. Our major findings are: (a) both array and DNA sequencing assay methods produced essentially the same origin distribution profile; (b) origin distribution is largely conserved (>70%) in all cell lines tested; (c) origins are enriched at the 5'ends of expressed genes and at evolutionarily conserved intergenic sequences; and (d) ChIP on chip experiments in MCF-7 showed an enrichment of H3K4Me3 and RNA Polymerase II chromatin binding sites at origins of DNA replication.

CONCLUSIONS/SIGNIFICANCE

Our results suggest that the program for origin activation is largely conserved among different cell types. Also, our work supports recent studies connecting transcription initiation with replication, and in addition suggests that evolutionarily conserved intergenic sequences have the potential to participate in origin selection. Overall, our observations suggest that replication origin selection is a stochastic process significantly dependent upon local accessibility to replication factors.

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 origins oriGNAI3 and oriB of the mammalian AMPD2 locus nested in a region of straight DNA flanked by intrinsically bent DNA sites

The aim of this work was to determine whether intrinsically bent DNA sites are present at, or close to, the mammalian replication origins oriGNAI3 and oriB in the Chinese hamster AMPD2 locus. Using an electrophoretic mobility shift assay and in silico analysis, we located four intrinsically bent DNA sites (b1 to b4) in a fragment that contains the oriGNAI3 and one site (b5) proximal to oriB. The helical parameters show that each bent DNA site is curved in a left-handed superhelical writhe. A 2D projection of 3D fragment trajectories revealed that oriGNAI3 is located in a relatively straight segment flanked by bent sites b1 and b2, which map in previously identified Scaffold/Matrix Attachment Region. Sites b3 and b4 are located approximately 2 kb downstream and force the fragment into a strong closed loop structure. The b5 site is also located in an S/MAR that is found just downstream of oriB.

The dinucleotide CG as a genomic signalling module

The operon model proposed the existence of a category of proteins that control gene expression by interacting with specific DNA sequences. Since then, a large number of transcription factors recognizing a diversity of sequence motifs have been discovered. This article discusses an unusually short protein recognition sequence, 5'CG, which is read by multiple DNA binding proteins. CG exists in three distinct chemical states, two of which bind mutually exclusively to proteins that modulate chromatin structure. Non-methylated CG, which is highly concentrated at CpG island promoters, recruits enzymes that create the mark of promoter activity, trimethyl-lysine 4 of histone H3. Methylated CG, on the other hand, is a gene silencing mark and accordingly recruits enzymes that deacetylate histones. Thus, CG, despite its simplicity, has the properties of a genome-wide signalling module that adds a layer of positive or negative control over gene expression.

Chk1 inhibits replication factory activation but allows dormant origin firing in existing factories

At low levels of replication stress, Chk1 favors resolving problems at stalled replication forks over initiating origin firing in unreplicated areas of the genome.

Replication origins are licensed by loading MCM2-7 hexamers before entry into S phase. However, only ∼10% of licensed origins are normally used in S phase, with the others remaining dormant. When fork progression is inhibited, dormant origins initiate nearby to ensure that all of the DNA is eventually replicated. In apparent contrast, replicative stress activates ataxia telangiectasia and rad-3–related (ATR) and Chk1 checkpoint kinases that inhibit origin firing. In this study, we show that at low levels of replication stress, ATR/Chk1 predominantly suppresses origin initiation by inhibiting the activation of new replication factories, thereby reducing the number of active factories. At the same time, inhibition of replication fork progression allows dormant origins to initiate within existing replication factories. The inhibition of new factory activation by ATR/Chk1 therefore redirects replication toward active factories where forks are inhibited and away from regions that have yet to start replication. This minimizes the deleterious consequences of fork stalling and prevents similar problems from arising in unreplicated regions of the genome.

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.

Sister acts: coordinating DNA replication and cohesion establishment

The ring-shaped cohesin complex links sister chromatids and plays crucial roles in homologous recombination and mitotic chromosome segregation. In cycling cells, cohesin's ability to generate cohesive linkages is restricted to S phase and depends on loading and establishment factors that are intimately connected to DNA replication. Here we review how cohesin is regulated by the replication machinery, as well as recent evidence that cohesin itself influences how chromosomes are replicated.

Pre-replication complex proteins assemble at regions of low nucleosome occupancy within the Chinese hamster dihydrofolate reductase initiation zone

Genome-scale mapping of pre-replication complex proteins has not been reported in mammalian cells. Poor enrichment of these proteins at specific sites may be due to dispersed binding, poor epitope availability or cell cycle stage-specific binding. Here, we have mapped sites of biotin-tagged ORC and MCM protein binding in G1-synchronized populations of Chinese hamster cells harboring amplified copies of the dihydrofolate reductase (DHFR) locus, using avidin-affinity purification of biotinylated chromatin followed by high-density microarray analysis across the DHFR locus. We have identified several sites of significant enrichment for both complexes distributed throughout the previously identified initiation zone. Analysis of the frequency of initiations across stretched DNA fibers from the DHFR locus confirmed a broad zone of de-localized initiation activity surrounding the sites of ORC and MCM enrichment. Mapping positions of mononucleosomal DNA empirically and computing nucleosome-positioning information in silico revealed that ORC and MCM map to regions of low measured and predicted nucleosome occupancy. Our results demonstrate that specific sites of ORC and MCM enrichment can be detected within a mammalian intitiation zone, and suggest that initiation zones may be regions of generally low nucleosome occupancy where flexible nucleosome positioning permits flexible pre-RC assembly sites.

Evaluating genome-scale approaches to eukaryotic DNA replication

Mechanisms regulating where and when eukaryotic DNA replication initiates remain a mystery. Recently, genome-scale methods have been brought to bear on this problem. The identification of replication origins and their associated proteins in yeasts is a well-integrated investigative tool, but corresponding data sets from multicellular organisms are scarce. By contrast, standardized protocols for evaluating replication timing have generated informative data sets for most eukaryotic systems. Here, I summarize the genome-scale methods that are most frequently used to analyse replication in eukaryotes, the kinds of questions each method can address and the technical hurdles that must be overcome to gain a complete understanding of the nature of eukaryotic replication origins.

Chk1 promotes replication fork progression by controlling replication initiation

DNA replication starts at initiation sites termed replication origins. Metazoan cells contain many more potential origins than are activated (fired) during each S phase. Origin activation is controlled by the ATR checkpoint kinase and its downstream effector kinase Chk1, which suppresses origin firing in response to replication blocks and during normal S phase by inhibiting the cyclin-dependent kinase Cdk2. In addition to increased origin activation, cells deficient in Chk1 activity display reduced rates of replication fork progression. Here we investigate the causal relationship between increased origin firing and reduced replication fork progression. We use the Cdk inhibitor roscovitine or RNAi depletion of Cdc7 to inhibit origin firing in Chk1-inhibited or RNAi-depleted cells. We report that Cdk inhibition and depletion of Cdc7 can alleviate the slow replication fork speeds in Chk1-deficient cells. Our data suggest that increased replication initiation leads to slow replication fork progression and that Chk1 promotes replication fork progression during normal S phase by controlling replication origin activity.

Ribonucleotide reductase is not limiting for mitochondrial DNA copy number in mice

Ribonucleotide reductase (RNR) is the rate-limiting enzyme in deoxyribonucleoside triphosphate (dNTP) biosynthesis, with important roles in nuclear genome maintenance. RNR is also essential for maintenance of mitochondrial DNA (mtDNA) in mammals. The mechanisms regulating mtDNA copy number in mammals are only being discovered. In budding yeast, RNR overexpression resulted in increased mtDNA levels, and rescued the disease phenotypes caused by a mutant mtDNA polymerase. This raised the question of whether mtDNA copy number increase by RNR induction could be a strategy for treating diseases with mtDNA mutations. We show here that high-level overexpression of RNR subunits (Rrm1, Rrm2 and p53R2; separately or in different combinations) in mice does not result in mtDNA copy number elevation. Instead, simultaneous expression of two RNR subunits leads to imbalanced dNTP pools and progressive mtDNA depletion in the skeletal muscle, without mtDNA mutagenesis. We also show that endogenous RNR transcripts are downregulated in response to large increases of mtDNA in mice, which is indicative of nuclear-mitochondrial crosstalk with regard to mtDNA copy number. Our results establish that RNR is not limiting for mtDNA copy number in mice, and provide new evidence for the importance of balanced dNTP pools in mtDNA maintenance in postmitotic tissues.

Nucleotide supply, not local histone acetylation, sets replication origin usage in transcribed regions

In eukaryotes, only a fraction of replication origins fire at each S phase. Local histone acetylation was proposed to control firing efficiency of origins, but conflicting results were obtained. We report that local histone acetylation does not reflect origin efficiencies along the adenosine monophosphate deaminase 2 locus in mammalian fibroblasts. Reciprocally, modulation of origin efficiency does not affect acetylation. However, treatment with a deacetylase inhibitor changes the initiation pattern. We demonstrate that this treatment alters pyrimidine biosynthesis and decreases fork speed, which recruits latent origins. Our findings reconcile results that seemed inconsistent and reveal an unsuspected effect of deacetylase inhibitors on replication dynamics.

Spi-1/PU.1 oncogene accelerates DNA replication fork elongation and promotes genetic instability in the absence of DNA breakage

The multistage process of cancer formation is driven by the progressive acquisition of somatic mutations. Replication stress creates genomic instability in mammals. Using a well-defined multistep leukemia model driven by Spi-1/PU.1 overexpression in the mouse and Spi-1/PU.1-overexpressing human leukemic cells, we investigated the relationship between DNA replication and cancer progression. Here, using DNA molecular combing and flow cytometry methods, we show that Spi-1 increases the speed of replication by acting specifically on elongation rather than enhancing origin firing. This shortens the S-phase duration. Combining data from Spi-1 knockdown in murine cells with Spi-1 overexpression in human cells, we provide evidence that inappropriate Spi-1 expression is directly responsible for the replication alteration observed. Importantly, the acceleration of replication progression coincides with an increase in the frequency of genomic mutations without inducing DNA breakage. Thus, we propose that the hitherto unsuspected role for spi-1 oncogene in promoting replication elongation and genomic mutation promotes blastic progression during leukemic development.

Single-molecule analysis reveals changes in the DNA replication program for the POU5F1 locus upon human embryonic stem cell differentiation

Human embryonic stem cells (hESCs), due to their pluripotent nature, represent a particularly relevant model system to study the relationship between the replication program and differentiation state. Here, we define the basic properties of the replication program in hESCs and compare them to the programs of hESC-derived multipotent cells (neural rosette cells) and primary differentiated cells (microvascular endothelial cells [MECs]). We characterized three genomic loci: two pluripotency regulatory genes, POU5F1 (OCT4) and NANOG, and the IGH locus, a locus that is transcriptionally active specifically in B-lineage cells. We applied a high-resolution approach to capture images of individual replicated DNA molecules. We demonstrate that for the loci studied, several basic properties of replication, including the average speed of replication forks and the average density of initiation sites, were conserved among the cells analyzed. We also demonstrate, for the first time, the presence of initiation zones in hESCs. However, significant differences were evident in other aspects of replication for the DNA segment containing the POU5F1 gene. Specifically, the locations of centers of initiation zones and the direction of replication fork progression through the POU5F1 gene were conserved in two independent hESC lines but were different in hESC-derived multipotent cells and MECs. Thus, our data identify features of the replication program characteristic of hESCs and define specific changes in replication during hESC differentiation.

Eukaryotic chromosome DNA replication: where, when, and how?

DNA replication is central to cell proliferation. Studies in the past six decades since the proposal of a semiconservative mode of DNA replication have confirmed the high degree of conservation of the basic machinery of DNA replication from prokaryotes to eukaryotes. However, the need for replication of a substantially longer segment of DNA in coordination with various internal and external signals in eukaryotic cells has led to more complex and versatile regulatory strategies. The replication program in higher eukaryotes is under a dynamic and plastic regulation within a single cell, or within the cell population, or during development. We review here various regulatory mechanisms that control the replication program in eukaryotes and discuss future directions in this dynamic field.

A winding road to origin discovery

Studies in our laboratory over the last three decades have shown that the Chinese hamster dihydrofolate reductase (DHFR) origin of replication corresponds to a broad zone of inefficient initiation sites distributed throughout the spacer between the convergently transcribed DHFR and 2BE2121 genes. It is clear from mutational analysis that none of these sites is genetically required for controlling origin activity. However, the integrity of the promoter of the DHFR gene is needed to activate the downstream origin, while the 3' processing signals prevent invasion and inactivation of the downstream origin by transcription forks. Several other origins in metazoans have been shown to correspond to zones of inefficient sites, while a different subset appears to be similar to the fixed replicators that characterize origins in S. cerevisiae and lower organisms. These observations have led us to suggest a model in which the mammalian genome is dotted with a hierarchy of degenerate, redundant, and inefficient replicators at intervals of a kilobase or less, some of which may have evolved to be highly circumscribed and efficient. The activities of initiation sites are proposed to be largely regulated by local transcription and chromatin architecture. Recently, we and others have devised strategies for identifying active origins on a genome-wide scale in order to define their distributions between fixed and dispersive origin types and to detect relationships among origins, genes, and epigenetic markers. The global pictures emerging are suggestive but far from complete and appear to be plagued by some of the same uncertainties that have led to conflicting views of individual origins in the past (particularly DHFR). In this paper, we will trace the history of origin discovery in mammalian genomes, primarily using the well-studied DHFR origin as a model, because it has been analyzed by nearly every available origin mapping technique in several different laboratories, while many origins have been identified by only one. We will address the strengths and shortcomings of the various methods utilized to identify and characterize origins in complex genomes and will point out how we and others were sometimes led astray by false assumptions and biases, as well as insufficient information. The goal is to help guide future experiments that will provide a truly comprehensive and accurate portrait of origins and their regulation. After all, in the words of George Santayana, "Those who do not learn from history are doomed to repeat it."

Defining replication origin efficiency using DNA fiber assays

The timely duplication of eukaryotic genomes depends on the coordinated activation of thousands of replication origins distributed along the chromosomes. Origin activation follows a temporal program that is imposed by the chromosomal context and is under the control of S-phase checkpoints. Although the general mechanisms regulating DNA replication are now well-understood at the level of individual origins, little is known on the coordination of thousands of initiation events at a genome-wide level. Recent studies using DNA combing and other single-molecule assays have shown that eukaryotic genomes contain a large excess of replication origins. Most of these origins remain "dormant" in normal growth conditions but are activated when fork progression is impeded. In this review, we discuss how DNA fiber technologies have changed our view of eukaryotic replication programs and how origin redundancy contributes to the maintenance of genome integrity in eukaryotic cells.

Inhibition of histone deacetylase in cancer cells slows down replication forks, activates dormant origins, and induces DNA damage

Protein acetylation is a reversible process regulated by histone deacetylases (HDAC) that is often altered in human cancers. Suberoylanilide hydroxamic acid (SAHA) is the first HDAC inhibitor to be approved for clinical use as an anticancer agent. Given that histone acetylation is a key determinant of chromatin structure, we investigated how SAHA may affect DNA replication and integrity to gain deeper insights into the basis for its anticancer activity. Nuclear replication factories were visualized with confocal immunofluorescence microscopy and single-replicon analyses were conducted by genome-wide molecular combing after pulse labeling with two thymidine analogues. We found that pharmacologic concentrations of SAHA induce replication-mediated DNA damage with activation of histone gammaH2AX. Single DNA molecule analyses indicated slowdown in replication speed along with activation of dormant replication origins in response to SAHA. Similar results were obtained using siRNA-mediated depletion of HDAC3 expression, implicating this HDAC member as a likely target in the SAHA response. Activation of dormant origins was confirmed by molecular analyses of the beta-globin locus control region. Our findings demonstrate that SAHA produces profound alterations in DNA replication that cause DNA damage, establishing a critical link between robust chromatin acetylation and DNA replication in human cancer cells.

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.

Increased origin activity in transformed versus normal cells: identification of novel protein players involved in DNA replication and cellular transformation

Using libraries of replication origins generated previously, we identified three clones that supported the autonomous replication of their respective plasmids in transformed, but not in normal cells. Assessment of their in vivo replication activity by in situ chromosomal DNA replication assays revealed that the chromosomal loci corresponding to these clones coincided with chromosomal replication origins in all cell lines, which were more active by 2-3-fold in the transformed by comparison to the normal cells. Evaluation of pre-replication complex (pre-RC) protein abundance at these origins in transformed and normal cells by chromatin immunoprecipitation assays, using anti-ORC2, -cdc6 and -cdt1 antibodies, showed that they were bound by these pre-RC proteins in all cell lines, but a 2-3-fold higher abundance was observed in the transformed by comparison to the normal cells. Electrophoretic mobility shift assays (EMSAs) performed on the most efficiently replicating clone, using nuclear extracts from the transformed and normal cells, revealed the presence of a DNA replication complex in transformed cells, which was barely detectable in normal cells. Subsequent supershift EMSAs suggested the presence of transformation-specific complexes. Mass spectrometric analysis of these complexes revealed potential new protein players involved in DNA replication that appear to correlate with cellular transformation.

Hydroxyurea-stalled replication forks become progressively inactivated and require two different RAD51-mediated pathways for restart and repair

Faithful DNA replication is essential to all life. Hydroxyurea (HU) depletes the cells of dNTPs, which initially results in stalled replication forks that, after prolonged treatment, collapse into DSBs. Here, we report that stalled replication forks are efficiently restarted in a RAD51-dependent process that does not trigger homologous recombination (HR). The XRCC3 protein, which is required for RAD51 foci formation, is also required for replication restart of HU-stalled forks, suggesting that RAD51-mediated strand invasion supports fork restart. In contrast, replication forks collapsed by prolonged replication blocks do not restart, and global replication is rescued by new origin firing. We find that RAD51-dependent HR is triggered for repair of collapsed replication forks, without apparent restart. In conclusion, our data suggest that restart of stalled replication forks and HR repair of collapsed replication forks require two distinct RAD51-mediated pathways.

Timing of replication is a determinant of neutral substitution rates but does not explain slow Y chromosome evolution in rodents

Mutation rates, assayed as substitution rates of putatively neutral sites, are highly variable around mammalian genomes: There is heterogeneity between genes, between autosomes, and between X, Y, and autosomes. The differences between X, Y, and autosomes are typically assumed to reflect the greater number of cell divisions in the male germ-line. Such an effect can neither account for within-autosome differences nor does it predict the differences between X, Y, and autosome observed in rodents. It has recently been proposed that in primates, the time during S-phase when a gene is replicated is an important determinant of neutral rates of evolution. Here we ask 1) whether we can replicate this result in rodents, 2) whether different autosomes replicate on average at different times, and 3) whether this might explain differences in their substitution rates. Finally we ask 4) whether X, Y, and autosome replicate at different times and 5) whether any difference might explain why the number of replication events alone cannot explain their substitution rates. We find that, as in primates, autosomal intronic rates of evolution increase significantly during S-phase. Different autosomes do have different average replication times, and together with rearrangement, this is a significant predictor of between-autosome differences in substitution rate. Although we find that autosomal, X-, and Y-linked genes replicate at different times, it is paradoxical that the Y-linked genes replicate latest, and replicate more often, but are not especially fast evolving. These results support the hypothesis that replication timing is an important source of substitution rate heterogeneity.

Analysis of replication factories in human cells by super-resolution light microscopy

Background

DNA replication in human cells is performed in discrete sub-nuclear locations known as replication foci or factories. These factories form in the nucleus during S phase and are sites of DNA synthesis and high local concentrations of enzymes required for chromatin replication. Why these structures are required, and how they are organised internally has yet to be identified. It has been difficult to analyse the structure of these factories as they are small in size and thus below the resolution limit of the standard confocal microscope. We have used stimulated emission depletion (STED) microscopy, which improves on the resolving power of the confocal microscope, to probe the structure of these factories at sub-diffraction limit resolution.

Results

Using immunofluorescent imaging of PCNA (proliferating cell nuclear antigen) and RPA (replication protein A) we show that factories are smaller in size (approximately 150 nm diameter), and greater in number (up to 1400 in an early S- phase nucleus), than is determined by confocal imaging. The replication inhibitor hydroxyurea caused an approximately 40% reduction in number and a 30% increase in diameter of replication factories, changes that were not clearly identified by standard confocal imaging.

Conclusions

These measurements for replication factory size now approach the dimensions suggested by electron microscopy. This agreement between these two methods, that use very different sample preparation and imaging conditions, suggests that we have arrived at a true measurement for the size of these structures. The number of individual factories present in a single nucleus that we measure using this system is greater than has been previously reported. This analysis therefore suggests that each replication factory contains fewer active replication forks than previously envisaged.

An interfering activity against lymphocytic choriomeningitis virus replication associated with enhanced mutagenesis

Previous studies have documented that, in the presence of the mutagenic base analogue 5-fluorouracil (FU), lymphocytic choriomeningitis virus (LCMV) that persisted in BHK-21 cells decreased its infectivity to a larger extent than intracellular viral RNA levels, prior to virus extinction. This observation, together with in silico simulations, led to the proposal of the lethal defection model of virus extinction. This model suggests the participation of defective-interfering genomes in the loss of infectivity by increased mutagenesis. Since LCMV naturally produces defective-interfering particles, it was important to show that a capacity to interfere is produced in association with FU treatment. Here, we document that BHK-21 cells persistently infected with LCMV grown in the presence of FU, but not in its absence, generated an interfering activity that suppressed LCMV infectivity. Interference was specific for LCMV and was sensitive to UV irradiation and its activity was dose- and time-dependent. The interfering preparations produced positive LCMV immunofluorescence and viral particles seen by electron microscopy when used to infect cells, despite some preparations being devoid of detectable infectivity. Interference did not involve significant increases of mutant spectrum complexity, as predicted by the lethal defection model. The results provide support for a specific interference associated with LCMV when the virus replicates in the presence of FU. The excess of interference relative to that observed in the absence of FU is necessary for LCMV extinction.

The cell cycle is a redox cycle: linking phase-specific targets to cell fate

Reactive oxygen species (ROS) regulate the strength and duration of signaling through redox-dependent signal transduction pathways via the cyclic oxidation/reduction of cysteine residues in kinases, phosphatases, and other regulatory factors. Signaling circuits may be segregated in organelles or other subcellular domains with distinct redox states, permitting them to respond independently to changes in the oxidation state of two major thiol reductants, glutathione and thioredoxin. Studies in yeast, and in complex eukaryotes, show that oscillations in oxygen consumption, energy metabolism, and redox state are intimately integrated with cell cycle progression. Because signaling pathways play specific roles in different phases of the cell cycle and the hierarchy of redox-dependent regulatory checkpoints changes during cell cycle progression, the effects of ROS on cell fate vary during the cell cycle. In G1, ROS stimulate mitogenic pathways that control the activity of cyclin-dependent kinases (CDKs) and phosphorylation of the retinoblastoma protein (pRB), thereby regulating S-phase entry. In response to oxidative stress, Nrf2 and Foxo3a promote cell survival by inducing the expression of antioxidant enzymes and factors involved in cell cycle withdrawal, such as the cyclin-dependent kinase inhibitor (CKI) p27. In S phase, ROS induce S-phase arrest via PP2A-dependent dephosphorylation of pRB. In precancerous cells, unconstrained mitogenic signaling by activated oncogenes induces replication stress in S phase, which activates the DNA-damage response and induces cell senescence. A number of studies suggest that interactions of ROS with the G1 CDK/CKI network play a fundamental role in senescence, which is considered a barrier to tumorigenesis. Adaptive responses and loss of checkpoint proteins such as p53 and p16(INK4a) allow tumor cells to tolerate constitutive mitogenic signaling and enhanced production of ROS, leading to altered redox status in many fully transformed cells. Alterations in oxidant and energy metabolism of cancer cells have emerged as fertile ground for new therapeutic targets. The present challenge is to identify redox-dependent targets relevant to each cell cycle phase, to understand how these targets control fate decisions, and to describe the mechanisms that link metabolism to cell cycle progression.

Y RNA functions at the initiation step of mammalian chromosomal DNA replication

Non-coding Y RNAs have recently been identified as essential novel factors for chromosomal DNA replication in mammalian cell nuclei, but mechanistic details of their function have not been defined. Here, we identify the execution point for Y RNA function during chromosomal DNA replication in a mammalian cell-free system. We determined the effect of degradation of Y3 RNA on replication origin activation and on fork progression rates at single-molecule resolution by DNA combing and nascent-strand analysis. Degradation of Y3 RNA inhibits the establishment of new DNA replication forks at the G1- to S-phase transition and during S phase. This inhibition is negated by addition of exogenous Y1 RNA. By contrast, progression rates of DNA replication forks are not affected by degradation of Y3 RNA or supplementation with exogenous Y1 RNA. These data indicate that Y RNAs are required for the establishment, but not for the elongation, of chromosomal DNA replication forks in mammalian cell nuclei. We conclude that the execution point for non-coding Y RNA function is the activation of chromosomal DNA replication origins.

Topoisomerase I suppresses genomic instability by preventing interference between replication and transcription

Topoisomerase I (Top1) is a key enzyme in functioning at the interface between DNA replication, transcription and mRNA maturation. Here, we show that Top1 suppresses genomic instability in mammalian cells by preventing a conflict between transcription and DNA replication. Using DNA combing and ChIP (chromatin immunoprecipitation)-on-chip, we found that Top1-deficient cells accumulate stalled replication forks and chromosome breaks in S phase, and that breaks occur preferentially at gene-rich regions of the genome. Notably, these phenotypes were suppressed by preventing the formation of RNA-DNA hybrids (R-loops) during transcription. Moreover, these defects could be mimicked by depletion of the splicing factor ASF/SF2 (alternative splicing factor/splicing factor 2), which interacts functionally with Top1. Taken together, these data indicate that Top1 prevents replication fork collapse by suppressing the formation of R-loops in an ASF/SF2-dependent manner. We propose that interference between replication and transcription represents a major source of spontaneous replication stress, which could drive genomic instability during the early stages of tumorigenesis.

The role of the Fanconi anemia network in the response to DNA replication stress

Fanconi anemia is a genetically heterogeneous disorder associated with chromosome instability and a highly elevated risk for developing cancer. The mutated genes encode proteins involved in the cellular response to DNA replication stress. Fanconi anemia proteins are extensively connected with DNA caretaker proteins, and appear to function as a hub for the coordination of DNA repair with DNA replication and cell cycle progression. At a molecular level, however, the raison d'être of Fanconi anemia proteins still remains largely elusive. The thirteen Fanconi anemia proteins identified to date have not been embraced into a single and defined biological process. To help put the Fanconi anemia puzzle into perspective, we begin this review with a summary of the strategies employed by prokaryotes and eukaryotes to tolerate obstacles to the progression of replication forks. We then summarize what we know about Fanconi anemia with an emphasis on biochemical aspects, and discuss how the Fanconi anemia network, a late acquisition in evolution, may function to permit the faithful and complete duplication of our very large vertebrate chromosomes.

DNA replication initiation

DNA replication is fundamental to cellular life on earth, and replication initiation provides the primary point of control over this process. Replication initiation in all organisms involves the interaction of initiator proteins with one or more origins of replication in the DNA, with subsequent regulated assembly of two replisome complexes at each origin, melting of the DNA, and primed initiation of DNA synthesis on leading and lagging strands. Archaea and Eukarya share homologous systems for DNA replication initiation, but differ in the complexity of these; Bacteria appear to have analogous, rather than homologous, mechanisms for replication initiation. This chapter provides an overview of current knowledge of initiation of chromosomal DNA replication in the three domains of life.

Mammalian telomeres resemble fragile sites and require TRF1 for efficient replication

Telomeres protect chromosome ends through the interaction of telomeric repeats with shelterin, a protein complex that represses DNA damage signaling and DNA repair reactions. The telomeric repeats are maintained by telomerase, which solves the end replication problem. We report that the TTAGGG repeat arrays of mammalian telomeres pose a challenge to the DNA replication machinery, giving rise to replication-dependent defects that resemble those of aphidicolin-induced common fragile sites. Gene deletion experiments showed that efficient duplication of telomeres requires the shelterin component TRF1. Without TRF1, telomeres activate the ATR kinase in S phase and show a fragile-site phenotype in metaphase. Single-molecule analysis of replicating telomeres showed that TRF1 promotes efficient replication of TTAGGG repeats and prevents fork stalling. Two helicases implicated in the removal of G4 DNA structures, BLM and RTEL1, were required to repress the fragile-telomere phenotype. These results identify a second telomere replication problem that is solved by the shelterin component TRF1.

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.