A winding road to origin discovery

The CMG helicase and cancer: a tumor "engine" and weakness with missing mutations

The replicative Cdc45-MCM-GINS (CMG) helicase is a large protein complex that functions in the DNA melting and unwinding steps as a component of replisomes during DNA replication in mammalian cells. Although the CMG performs this important role in cell growth, the CMG is not a simple bystander in cell cycle events. Components of the CMG, specifically the MCM precursors, are also involved in maintaining genomic stability by regulating DNA replication fork speeds, facilitating recovery from replicative stresses, and preventing consequential DNA damage. Given these important functions, MCM/CMG complexes are highly regulated by growth factors such as TGF-ß1 and by signaling factors such as Myc, Cyclin E, and the retinoblastoma protein. Mismanagement of MCM/CMG complexes when these signaling mediators are deregulated, and in the absence of the tumor suppressor protein p53, leads to increased genomic instability and is a contributor to tumorigenic transformation and tumor heterogeneity. The goal of this review is to provide insight into the mechanisms and dynamics by which the CMG is regulated during its assembly and activation in mammalian genomes, and how errors in CMG regulation due to oncogenic changes promote tumorigenesis. Finally, and most importantly, we highlight the emerging understanding of the CMG helicase as an exploitable vulnerability and novel target for therapeutic intervention in cancer.

Genome-wide measurement of DNA replication fork directionality and quantification of DNA replication initiation and termination with Okazaki fragment sequencing

Studying the dynamics of genome replication in mammalian cells has been historically challenging. To reveal the location of replication initiation and termination in the human genome, we developed Okazaki fragment sequencing (OK-seq), a quantitative approach based on the isolation and strand-specific sequencing of Okazaki fragments, the lagging strand replication intermediates. OK-seq quantitates the proportion of leftward- and rightward-oriented forks at every genomic locus and reveals the location and efficiency of replication initiation and termination events. Here we provide the detailed experimental procedures for performing OK-seq in unperturbed cultured human cells and budding yeast and the bioinformatics pipelines for data processing and computation of replication fork directionality. Furthermore, we present the analytical approach based on a hidden Markov model, which allows automated detection of ascending, descending and flat replication fork directionality segments revealing the zones of replication initiation, termination and unidirectional fork movement across the entire genome. These tools are essential for the accurate interpretation of human and yeast replication programs. The experiments and the data processing can be accomplished within six days. Besides revealing the genome replication program in fine detail, OK-seq has been instrumental in numerous studies unravelling mechanisms of genome stability, epigenome maintenance and genome evolution.

Determination of human DNA replication origin position and efficiency reveals principles of initiation zone organisation

Replication of the human genome initiates within broad zones of ∼150 kb. The extent to which firing of individual DNA replication origins within initiation zones is spatially stochastic or localised at defined sites remains a matter of debate. A thorough characterisation of the dynamic activation of origins within initiation zones is hampered by the lack of a high-resolution map of both their position and efficiency. To address this shortcoming, we describe a modification of initiation site sequencing (ini-seq), based on density substitution. Newly replicated DNA is rendered 'heavy-light' (HL) by incorporation of BrdUTP while unreplicated DNA remains 'light-light' (LL). Replicated HL-DNA is separated from unreplicated LL-DNA by equilibrium density gradient centrifugation, then both fractions are subjected to massive parallel sequencing. This allows precise mapping of 23,905 replication origins simultaneously with an assignment of a replication initiation efficiency score to each. We show that origin firing within early initiation zones is not randomly distributed. Rather, origins are arranged hierarchically with a set of very highly efficient origins marking zone boundaries. We propose that these origins explain much of the early firing activity arising within initiation zones, helping to unify the concept of replication initiation zones with the identification of discrete replication origin sites.

The genetic architecture of DNA replication timing in human pluripotent stem cells

DNA replication follows a strict spatiotemporal program that intersects with chromatin structure but has a poorly understood genetic basis. To systematically identify genetic regulators of replication timing, we exploited inter-individual variation in human pluripotent stem cells from 349 individuals. We show that the human genome's replication program is broadly encoded in DNA and identify 1,617 cis-acting replication timing quantitative trait loci (rtQTLs) - sequence determinants of replication initiation. rtQTLs function individually, or in combinations of proximal and distal regulators, and are enriched at sites of histone H3 trimethylation of lysines 4, 9, and 36 together with histone hyperacetylation. H3 trimethylation marks are individually repressive yet synergistically associate with early replication. We identify pluripotency-related transcription factors and boundary elements as positive and negative regulators of replication timing, respectively. Taken together, human replication timing is controlled by a multi-layered mechanism with dozens of effectors working combinatorially and following principles analogous to transcription regulation.

Human ORC/MCM density is low in active genes and correlates with replication time but does not delimit initiation zones

Eukaryotic DNA replication initiates during S phase from origins that have been licensed in the preceding G1 phase. Here, we compare ChIP-seq profiles of the licensing factors Orc2, Orc3, Mcm3, and Mcm7 with gene expression, replication timing, and fork directionality profiles obtained by RNA-seq, Repli-seq, and OK-seq. Both, the origin recognition complex (ORC) and the minichromosome maintenance complex (MCM) are significantly and homogeneously depleted from transcribed genes, enriched at gene promoters, and more abundant in early- than in late-replicating domains. Surprisingly, after controlling these variables, no difference in ORC/MCM density is detected between initiation zones, termination zones, unidirectionally replicating regions, and randomly replicating regions. Therefore, ORC/MCM density correlates with replication timing but does not solely regulate the probability of replication initiation. Interestingly, H4K20me3, a histone modification proposed to facilitate late origin licensing, was enriched in late-replicating initiation zones and gene deserts of stochastic replication fork direction. We discuss potential mechanisms specifying when and where replication initiates in human cells.

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.

How MCM loading and spreading specify eukaryotic DNA replication initiation sites

DNA replication origins strikingly differ between eukaryotic species and cell types. Origins are localized and can be highly efficient in budding yeast, are randomly located in early fly and frog embryos, which do not transcribe their genomes, and are clustered in broad (10-100 kb) non-transcribed zones, frequently abutting transcribed genes, in mammalian cells. Nonetheless, in all cases, origins are established during the G1-phase of the cell cycle by the loading of double hexamers of the Mcm 2-7 proteins (MCM DHs), the core of the replicative helicase. MCM DH activation in S-phase leads to origin unwinding, polymerase recruitment, and initiation of bidirectional DNA synthesis. Although MCM DHs are initially loaded at sites defined by the binding of the origin recognition complex (ORC), they ultimately bind chromatin in much greater numbers than ORC and only a fraction are activated in any one S-phase. Data suggest that the multiplicity and functional redundancy of MCM DHs provide robustness to the replication process and affect replication time and that MCM DHs can slide along the DNA and spread over large distances around the ORC. Recent studies further show that MCM DHs are displaced along the DNA by collision with transcription complexes but remain functional for initiation after displacement. Therefore, eukaryotic DNA replication relies on intrinsically mobile and flexible origins, a strategy fundamentally different from bacteria but conserved from yeast to human. These properties of MCM DHs likely contribute to the establishment of broad, intergenic replication initiation zones in higher eukaryotes.

The epigenetic regulation of autonomous replicons

The discovery of autonomous replicating sequences (ARSs) in Saccharomyces cerevisiae in 1979 was considered a milestone in unraveling the regulation of replication in eukaryotic cells. However, shortly afterwards it became obvious that in Saccharomyces pombe and all other higher organisms ARSs were not sufficient to initiate independent replication. Understanding the mechanisms of replication is a major challenge in modern cell biology and is also a prerequisite to developing application-oriented autonomous replicons for gene therapeutic treatments. This review will focus on the development of non-viral episomal vectors, their use in gene therapeutic applications and our current knowledge about their epigenetic regulation.

Effect of minichromosome maintenance protein 2 deficiency on the locations of DNA replication origins

Minichromosome maintenance (MCM) proteins are loaded onto chromatin during G1-phase and define potential locations of DNA replication initiation. MCM protein deficiency results in genome instability and high rates of cancer in mouse models. Here we develop a method of nascent strand capture and release and show that MCM2 deficiency reduces DNA replication initiation in gene-rich regions of the genome. DNA structural properties are shown to correlate with sequence motifs associated with replication origins and with locations that are preferentially affected by MCM2 deficiency. Reduced nascent strand density correlates with sites of recurrent focal CNVs in tumors arising in MCM2-deficient mice, consistent with a direct relationship between sites of reduced DNA replication initiation and genetic damage. Between 10% and 90% of human tumors, depending on type, carry heterozygous loss or mutation of one or more MCM2-7 genes, which is expected to compromise DNA replication origin licensing and result in elevated rates of genome damage at a subset of gene-rich locations.

The hunt for origins of DNA replication in multicellular eukaryotes

Origins of DNA replication (ORIs) occur at defined regions in the genome. Although DNA sequence defines the position of ORIs in budding yeast, the factors for ORI specification remain elusive in metazoa. Several methods have been used recently to map ORIs in metazoan genomes with the hope that features for ORI specification might emerge. These methods are reviewed here with analysis of their advantages and shortcomings. The various factors that may influence ORI selection for initiation of DNA replication are discussed.

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

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

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

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

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

Peaks cloaked in the mist: the landscape of mammalian replication origins

Replication of mammalian genomes starts at sites termed replication origins, which historically have been difficult to locate as a result of large genome sizes, limited power of genetic identification schemes, and rareness and fragility of initiation intermediates. However, origins are now mapped by the thousands using microarrays and sequencing techniques. Independent studies show modest concordance, suggesting that mammalian origins can form at any DNA sequence but are suppressed by read-through transcription or that they can overlap the 5' end or even the entire gene. These results require a critical reevaluation of whether origins form at specific DNA elements and/or epigenetic signals or require no such determinants.

Isolation and sequencing of active origins of DNA replication by nascent strand capture and release (NSCR)

Nascent strand capture and release (NSCR) is a method for isolation of short nascent strands to identify origins of DNA replication. The protocol provided involves isolation of total DNA, denaturation, size fractionation on a sucrose gradient, 5'-biotinylation of the appropriate size nucleic acids, binding to a streptavidin coated magnetic beads, intensive washing, and specific release of only the RNA-containing chimeric nascent strand DNA using ribonuclease I (RNase I). The method has been applied to mammalian cells derived from proliferative tissues and cell culture but could be used for any system where DNA replication is primed by a small RNA resulting in chimeric RNA-DNA molecules.

Exonuclease mutations in DNA polymerase epsilon reveal replication strand specific mutation patterns and human origins of replication

Tumors with somatic mutations in the proofreading exonuclease domain of DNA polymerase epsilon (POLE-exo*) exhibit a novel mutator phenotype, with markedly elevated TCT→TAT and TCG→TTG mutations and overall mutation frequencies often exceeding 100 mutations/Mb. Here, we identify POLE-exo* tumors in numerous cancers and classify them into two groups, A and B, according to their mutational properties. Group A mutants are found only in POLE, whereas Group B mutants are found in POLE and POLD1 and appear to be nonfunctional. In Group A, cell-free polymerase assays confirm that mutations in the exonuclease domain result in high mutation frequencies with a preference for C→A mutation. We describe the patterns of amino acid substitutions caused by POLE-exo* and compare them to other tumor types. The nucleotide preference of POLE-exo* leads to increased frequencies of recurrent nonsense mutations in key tumor suppressors such as TP53, ATM, and PIK3R1. We further demonstrate that strand-specific mutation patterns arise from some of these POLE-exo* mutants during genome duplication. This is the first direct proof of leading strand-specific replication by human POLE, which has only been demonstrated in yeast so far. Taken together, the extremely high mutation frequency and strand specificity of mutations provide a unique identifier of eukaryotic origins of replication.

Human genome replication proceeds through four chromatin states

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

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

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

Trapping DNA replication origins from the human genome

Synthesis of chromosomal DNA is initiated from multiple origins of replication in higher eukaryotes; however, little is known about these origins’ structures. We isolated the origin-derived nascent DNAs from a human repair-deficient cell line by blocking the replication forks near the origins using two different origin-trapping methods (i.e., UV- or chemical crosslinker-treatment and cell synchronization in early S phase using DNA replication inhibitors). Single-stranded DNAs (of 0.5–3 kb) that accumulated after such treatments were labeled with bromodeoxyuridine (BrdU). BrdU-labeled DNA was immunopurified after fractionation by alkaline sucrose density gradient centrifugation and cloned by complementary-strand synthesis and PCR amplification. Competitive PCR revealed an increased abundance of DNA derived from known replication origins (c-myc and lamin B2 genes) in the nascent DNA fractions from the UV-treated or crosslinked cells. Nucleotide sequences of 85 and 208 kb were obtained from the two libraries (I and II) prepared from the UV-treated log-phase cells and early S phase arrested cells, respectively. The libraries differed from each other in their G+C composition and replication-related motif contents, suggesting that differences existed between the origin fragments isolated by the two different origin-trapping methods. The replication activities for seven out of 12 putative origin loci from the early-S phase cells were shown by competitive PCR. We mapped 117 (library I) and 172 (library II) putative origin loci to the human genome; approximately 60% and 50% of these loci were assigned to the G-band and intragenic regions, respectively. Analyses of the flanking sequences of the mapped loci suggested that the putative origin loci tended to associate with genes (including conserved sites) and DNase I hypersensitive sites; however, poor correlations were found between such loci and the CpG islands, transcription start sites, and K27-acetylated histone H3 peaks.

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.

Genome-wide mapping of human DNA-replication origins: levels of transcription at ORC1 sites regulate origin selection and replication timing

We report the genome-wide mapping of ORC1 binding sites in mammals, by chromatin immunoprecipitation and parallel sequencing (ChIP-seq). ORC1 binding sites in HeLa cells were validated as active DNA replication origins (ORIs) using Repli-seq, a method that allows identification of ORI-containing regions by parallel sequencing of temporally ordered replicating DNA. ORC1 sites were universally associated with transcription start sites (TSSs) of coding or noncoding RNAs (ncRNAs). Transcription levels at the ORC1 sites directly correlated with replication timing, suggesting the existence of two classes of ORIs: those associated with moderate/high transcription levels (≥1 RNA copy/cell), firing in early S and mapping to the TSSs of coding RNAs; and those associated with low transcription levels (<1 RNA copy/cell), firing throughout the entire S and mapping to TSSs of ncRNAs. These findings are compatible with a scenario whereby TSS expression levels influence the efficiency of ORC1 recruitment at G₁ and the probability of firing during S.

Open chromatin structures regulate the efficiencies of pre-RC formation and replication initiation in Epstein-Barr virus

Studies of EBV replication origins demonstrate an excess of pre-replication complexes that are formed at flexible MNase-sensitive sites in the genome.

Whether or not metazoan replication initiates at random or specific but flexible sites is an unsolved question. The lack of sequence specificity in origin recognition complex (ORC) DNA binding complicates genome-scale chromatin immunoprecipitation (ChIP)-based studies. Epstein-Barr virus (EBV) persists as chromatinized minichromosomes that are replicated by the host replication machinery. We used EBV to investigate the link between zones of pre-replication complex (pre-RC) assembly, replication initiation, and micrococcal nuclease (MNase) sensitivity at different cell cycle stages in a genome-wide fashion. The dyad symmetry element (DS) of EBV’s latent origin, a well-established and very efficient pre-RC assembly region, served as an internal control. We identified 64 pre-RC zones that correlate spatially with 57 short nascent strand (SNS) zones. MNase experiments revealed that pre-RC and SNS zones were linked to regions of increased MNase sensitivity, which is a marker of origin strength. Interestingly, although spatially correlated, pre-RC and SNS zones were characterized by different features. We propose that pre-RCs are formed at flexible but distinct sites, from which only a few are activated per single genome and cell cycle.

Inferring where and when replication initiates from genome-wide replication timing data

Based on an analogy between DNA replication and one dimensional nucleation-and-growth processes, various attempts to infer the local initiation rate I(x,t) of DNA replication origins from replication timing data have been developed in the framework of phase transition kinetics theories. These works have all used curve-fit strategies to estimate I(x,t) from genome-wide replication timing data. Here, we show how to invert analytically the Kolmogorov-Johnson-Mehl-Avrami model and extract I(x,t) directly. Tests on both simulated and experimental budding-yeast data confirm the location and firing-time distribution of replication origins.

Replication fork polarity gradients revealed by megabase-sized U-shaped replication timing domains in human cell lines

In higher eukaryotes, replication program specification in different cell types remains to be fully understood. We show for seven human cell lines that about half of the genome is divided in domains that display a characteristic U-shaped replication timing profile with early initiation zones at borders and late replication at centers. Significant overlap is observed between U-domains of different cell lines and also with germline replication domains exhibiting a N-shaped nucleotide compositional skew. From the demonstration that the average fork polarity is directly reflected by both the compositional skew and the derivative of the replication timing profile, we argue that the fact that this derivative displays a N-shape in U-domains sustains the existence of large-scale gradients of replication fork polarity in somatic and germline cells. Analysis of chromatin interaction (Hi-C) and chromatin marker data reveals that U-domains correspond to high-order chromatin structural units. We discuss possible models for replication origin activation within U/N-domains. The compartmentalization of the genome into replication U/N-domains provides new insights on the organization of the replication program in the human genome.

Integrative analysis of gene amplification in Drosophila follicle cells: parameters of origin activation and repression

In metazoans, how replication origins are specified and subsequently activated is not well understood. Drosophila amplicons in follicle cells (DAFCs) are genomic regions that undergo rereplication to increase DNA copy number. We identified all DAFCs by comparative genomic hybridization, uncovering two new amplicons in addition to four known previously. The complete identification of all DAFCs enabled us to investigate these in vivo replicons with respect to parameters of transcription, localization of the origin recognition complex (ORC), and histone acetylation, yielding important insights into gene amplification as a metazoan replication model. Significantly, ORC is bound across domains spanning 10 or more kilobases at the DAFC rather than at a specific site. Additionally, ORC is bound at many regions that do not undergo amplification, and, in contrast to cell culture, these regions do not correlate with high gene expression. As a developmental strategy, gene amplification is not the predominant means of achieving high expression levels, even in cells capable of amplification. Intriguingly, we found that, in some strains, a new amplicon, DAFC-22B, does not amplify, a consequence of distant repression of ORC binding and origin activation. This repression is alleviated when a fragment containing the origin is placed in different genomic contexts.

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.

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.

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.