GINS motion reveals replication fork progression is remarkably uniform throughout the yeast genome

How the Eukaryotic Replisome Achieves Rapid and Efficient DNA Replication

The eukaryotic replisome is a molecular machine that coordinates the Cdc45-MCM-GINS (CMG) replicative DNA helicase with DNA polymerases α, δ, and ε and other proteins to copy the leading- and lagging-strand templates at rates between 1 and 2 kb min-1. We have now reconstituted this sophisticated machine with purified proteins, beginning with regulated CMG assembly and activation. We show that replisome-associated factors Mrc1 and Csm3/Tof1 are crucial for in vivo rates of replisome progression. Additionally, maximal rates only occur when DNA polymerase ε catalyzes leading-strand synthesis together with its processivity factor PCNA. DNA polymerase δ can support leading-strand synthesis, but at slower rates. DNA polymerase δ is required for lagging-strand synthesis, but surprisingly also plays a role in establishing leading-strand synthesis, before DNA polymerase ε engagement. We propose that switching between these DNA polymerases also contributes to leading-strand synthesis under conditions of replicative stress.

Model-based analysis of DNA replication profiles: predicting replication fork velocity and initiation rate by profiling free-cycling cells

Eukaryotic cells initiate DNA synthesis by sequential firing of hundreds of origins. This ordered replication is described by replication profiles, which measure the DNA content within a cell population. Here, we show that replication dynamics can be deduced from replication profiles of free-cycling cells. While such profiles lack explicit temporal information, they are sensitive to fork velocity and initiation capacity through the passive replication pattern, namely the replication of origins by forks emanating elsewhere. We apply our model-based approach to a compendium of profiles that include most viable budding yeast mutants implicated in replication. Predicted changes in fork velocity or initiation capacity are verified by profiling synchronously replicating cells. Notably, most mutants implicated in late (or early) origin effects are explained by global modulation of fork velocity or initiation capacity. Our approach provides a rigorous framework for analyzing DNA replication profiles of free-cycling cells.

Pif1-family helicases cooperatively suppress widespread replication-fork arrest at tRNA genes

Saccharomyces cerevisiae encodes two distinct Pif1-family helicases – Pif1 and Rrm3 – which have been reported to play distinct roles in numerous nuclear processes. Here, we systematically characterize the roles of Pif1 helicases in replisome progression and lagging-strand synthesis in S. cerevisiae. We demonstrate that either Pif1 or Rrm3 redundantly stimulate strand-displacement by DNA polymerase δ during lagging-strand synthesis. By analyzing replisome mobility in pif1 and rrm3 mutants, we show that Rrm3, with a partially redundant contribution from Pif1, suppresses widespread terminal arrest of the replisome at tRNA genes. Although both head-on and codirectional collisions induce replication fork arrest at tRNA genes, head-on collisions arrest a higher proportion of replisomes. Consistent with this observation, we find that head-on collisions between tRNA transcription and replication are under-represented in the S. cerevisiae genome. We demonstrate that tRNA-mediated arrest is R-loop independent, and propose that replisome arrest and DNA damage are mechanistically separable.

The Eukaryotic Replication Machine

The cellular replicating machine, or "replisome," is composed of numerous different proteins. The core replication proteins in all cell types include a helicase, primase, DNA polymerases, sliding clamp, clamp loader, and single-strand binding (SSB) protein. The core eukaryotic replisome proteins evolved independently from those of bacteria and thus have distinct architectures and mechanisms of action. The core replisome proteins of the eukaryote include: an 11-subunit CMG helicase, DNA polymerase alpha-primase, leading strand DNA polymerase epsilon, lagging strand DNA polymerase delta, PCNA clamp, RFC clamp loader, and the RPA SSB protein. There are numerous other proteins that travel with eukaryotic replication forks, some of which are known to be involved in checkpoint regulation or nucleosome handling, but most have unknown functions and no bacterial analogue. Recent studies have revealed many structural and functional insights into replisome action. Also, the first structure of a replisome from any cell type has been elucidated for a eukaryote, consisting of 20 distinct proteins, with quite unexpected results. This review summarizes the current state of knowledge of the eukaryotic core replisome proteins, their structure, individual functions, and how they are organized at the replication fork as a machine.

Deciphering DNA replication dynamics in eukaryotic cell populations in relation with their averaged chromatin conformations

We propose a non-local model of DNA replication that takes into account the observed uncertainty on the position and time of replication initiation in eukaryote cell populations. By picturing replication initiation as a two-state system and considering all possible transition configurations, and by taking into account the chromatin’s fractal dimension, we derive an analytical expression for the rate of replication initiation. This model predicts with no free parameter the temporal profiles of initiation rate, replication fork density and fraction of replicated DNA, in quantitative agreement with corresponding experimental data from both S. cerevisiae and human cells and provides a quantitative estimate of initiation site redundancy. This study shows that, to a large extent, the program that regulates the dynamics of eukaryotic DNA replication is a collective phenomenon that emerges from the stochastic nature of replication origins initiation.

A Robust Workflow for Native Mass Spectrometric Analysis of Affinity-Isolated Endogenous Protein Assemblies

The central players in most cellular events are assemblies of macromolecules. Structural and functional characterization of these assemblies requires knowledge of their subunit stoichiometry and intersubunit connectivity. One of the most direct means for acquiring such information is so-called "native mass spectrometry (MS)", wherein the masses of the intact assemblies and parts thereof are accurately determined. It is of particular interest to apply native MS to the study of endogenous protein assemblies-i.e., those wherein the component proteins are expressed at endogenous levels in their natural functional states, rather than the overexpressed (sometimes partial) constructs commonly employed in classical structural studies, whose assembly can introduce stoichiometry artifacts and other unwanted effects. To date, the application of native MS to the elucidation of endogenous protein complexes has been limited by the difficulty in obtaining pristine cell-derived assemblies at sufficiently high concentrations for effective analysis. Here, to address this challenge, we present a robust workflow that couples rapid and efficient affinity isolation of endogenous protein complexes with a sensitive native MS readout. The resulting workflow has the potential to provide a wealth of data on the stoichiometry and intersubunit connectivity of endogenous protein assemblies-information that is key to successful integrative structural elucidation of biological systems.

The Eukaryotic Mismatch Recognition Complexes Track with the Replisome during DNA Synthesis

During replication, mismatch repair proteins recognize and repair mispaired bases that escape the proofreading activity of DNA polymerase. In this work, we tested the model that the eukaryotic mismatch recognition complex tracks with the advancing replisome. Using yeast, we examined the dynamics during replication of the leading strand polymerase Polε using Pol2 and the eukaryotic mismatch recognition complex using Msh2, the invariant protein involved in mismatch recognition. Specifically, we synchronized cells and processed samples using chromatin immunoprecipitation combined with custom DNA tiling arrays (ChIP-chip). The Polε signal was not detectable in G1, but was observed at active origins and replicating DNA throughout S-phase. The Polε signal provided the resolution to track origin firing timing and efficiencies as well as replisome progression rates. By detecting Polε and Msh2 dynamics within the same strain, we established that the mismatch recognition complex binds origins and spreads to adjacent regions with the replisome. In mismatch repair defective PCNA mutants, we observed that Msh2 binds to regions of replicating DNA, but the distribution and dynamics are altered, suggesting that PCNA is not the sole determinant for the mismatch recognition complex association with replicating regions, but may influence the dynamics of movement. Using biochemical and genomic methods, we provide evidence that both MutS complexes are in the vicinity of the replisome to efficiently repair the entire spectrum of mutations during replication. Our data supports the model that the proximity of MutSα/β to the replisome for the efficient repair of the newly synthesized strand before chromatin reassembles.

During replication, errors that escape the replication machinery are identified and repaired by DNA mismatch repair proteins. A mismatch in the helix is recognized by MutS homologs and subsequent events include excision of the error-containing strand followed by re-synthesis. A critical step in this process is directing repair to the newly synthesized strand. Current data suggest that transient discontinuities in the DNA backbone, known as nicks, generated during replication serve as the strand discrimination signals. Additionally, proteins that package DNA have the capacity to block mismatch recognition and are known to rapidly assemble behind the replication fork. Thus, there must be a short window of opportunity for the mismatch recognition complexes to scan for mismatches and access the strand discrimination signals. To address these issues, we tested the model that the mismatch recognition complexes track with the replisome. We employed high resolution genomic methods to determine that during replication, the mismatch recognition complexes bind origins of replication and advances with the replisome. The findings support the hypothesis that the mismatch recognition proteins track with the DNA replication machinery to accurately survey and repair the newly synthesized strands while the DNA is unpackaged and strand specificity signals are accessible.

Single-molecule motions and interactions in live cells reveal target search dynamics in mismatch repair

MutS is responsible for initiating the correction of DNA replication errors. To understand how MutS searches for and identifies rare base-pair mismatches, we characterized the dynamic movement of MutS and the replisome in real time using superresolution microscopy and single-molecule tracking in living cells. We report that MutS dynamics are heterogeneous in cells, with one MutS population exploring the nucleoid rapidly, while another MutS population moves to and transiently dwells at the replisome region, even in the absence of appreciable mismatch formation. Analysis of MutS motion shows that the speed of MutS is correlated with its separation distance from the replisome and that MutS motion slows when it enters the replisome region. We also show that mismatch detection increases MutS speed, supporting the model for MutS sliding clamp formation after mismatch recognition. Using variants of MutS and the replication processivity clamp to impair mismatch repair, we find that MutS dynamically moves to and from the replisome before mismatch binding to scan for errors. Furthermore, a block to DNA synthesis shows that MutS is only capable of binding mismatches near the replisome. It is well-established that MutS engages in an ATPase cycle, which is necessary for signaling downstream events. We show that a variant of MutS with a nucleotide binding defect is no longer capable of dynamic movement to and from the replisome, showing that proper nucleotide binding is critical for MutS to localize to the replisome in vivo. Our results provide mechanistic insight into the trafficking and movement of MutS in live cells as it searches for mismatches.

The dynamics of eukaryotic replication initiation: origin specificity, licensing, and firing at the single-molecule level

Eukaryotic replication initiation is highly regulated and dynamic. It begins with the origin recognition complex (ORC) binding DNA sites called origins of replication. ORC, together with Cdc6 and Cdt1, mediate pre-replicative complex (pre-RC) assembly by loading a double hexamer of Mcm2-7: the core of the replicative helicase. Here, we use single-molecule imaging to directly visualize Saccharomyces cerevisiae pre-RC assembly and replisome firing in real time. We show that ORC can locate and stably bind origins within large tracts of non-origin DNA and that Cdc6 drives ordered pre-RC assembly. We further show that the dynamics of the ORC-Cdc6 interaction dictate Mcm2-7 loading specificity and that Mcm2-7 double hexamers form preferentially at a native origin sequence. Finally, we demonstrate that single Mcm2-7 hexamers propagate bidirectionally, monotonically, and processively as constituents of active replisomes.

The Ty1 LTR-retrotransposon of budding yeast, Saccharomyces cerevisiae

Long-terminal repeat (LTR)-retrotransposons generate a copy of their DNA (cDNA) by reverse transcription of their RNA genome in cytoplasmic nucleocapsids. They are widespread in the eukaryotic kingdom and are the evolutionary progenitors of retroviruses [1]. The Ty1 element of the budding yeast Saccharomyces cerevisiae was the first LTR-retrotransposon demonstrated to mobilize through an RNA intermediate, and not surprisingly, is the best studied. The depth of our knowledge of Ty1 biology stems not only from the predominance of active Ty1 elements in the S. cerevisiae genome but also the ease and breadth of genomic, biochemical and cell biology approaches available to study cellular processes in yeast. This review describes the basic structure of Ty1 and its gene products, the replication cycle, the rapidly expanding compendium of host co-factors known to influence retrotransposition and the nature of Ty1's elaborate symbiosis with its host. Our goal is to illuminate the value of Ty1 as a paradigm to explore the biology of LTR-retrotransposons in multicellular organisms, where the low frequency of retrotransposition events presents a formidable barrier to investigations of retrotransposon biology.

Dynamic loading and redistribution of the Mcm2-7 helicase complex through the cell cycle

Eukaryotic replication origins are defined by the ORC-dependent loading of the Mcm2-7 helicase complex onto chromatin in G1. Paradoxically, there is a vast excess of Mcm2-7 relative to ORC assembled onto chromatin in G1. These excess Mcm2-7 complexes exhibit little co-localization with ORC or replication foci and can function as dormant origins. We dissected the mechanisms regulating the assembly and distribution of the Mcm2-7 complex in the Drosophila genome. We found that in the absence of cyclin E/Cdk2 activity, there was a 10-fold decrease in chromatin-associated Mcm2-7 relative to the levels found at the G1/S transition. The minimal amounts of Mcm2-7 loaded in the absence of cyclin E/Cdk2 activity were strictly localized to ORC binding sites. In contrast, cyclin E/Cdk2 activity was required for maximal loading of Mcm2-7 and a dramatic genome-wide reorganization of the distribution of Mcm2-7 that is shaped by active transcription. Thus, increasing cyclin E/Cdk2 activity over the course of G1 is not only critical for Mcm2-7 loading, but also for the distribution of the Mcm2-7 helicase prior to S-phase entry.

Systems cell biology

Systems cell biology melds high-throughput experimentation with quantitative analysis and modeling to understand many critical processes that contribute to cellular organization and dynamics. Recently, there have been several advances in technology and in the application of modeling approaches that enable the exploration of the dynamic properties of cells. Merging technology and computation offers an opportunity to objectively address unsolved cellular mechanisms, and has revealed emergent properties and helped to gain a more comprehensive and fundamental understanding of cell biology.

Compromised RNA polymerase III complex assembly leads to local alterations of intergenic RNA polymerase II transcription in Saccharomyces cerevisiae

Background

Assembled RNA polymerase III (Pol III) complexes exert local effects on chromatin processes, including influencing transcription of neighboring RNA polymerase II (Pol II) transcribed genes. These properties have been designated as ‘extra-transcriptional’ effects of the Pol III complex. Previous coding sequence microarray studies using Pol III factor mutants to determine global effects of Pol III complex assembly on Pol II promoter activity revealed only modest effects that did not correlate with the proximity of Pol III complex binding sites.

Results

Given our recent results demonstrating that tDNAs block progression of intergenic Pol II transcription, we hypothesized that extra-transcriptional effects within intergenic regions were not identified in the microarray study. To reconsider global impacts of Pol III complex binding, we used RNA sequencing to compare transcriptomes of wild type versus Pol III transcription factor TFIIIC depleted mutants. The results reveal altered intergenic Pol II transcription near TFIIIC binding sites in the mutant strains, where we observe readthrough of upstream transcripts that normally terminate near these sites, 5′- and 3′-extended transcripts, and de-repression of adjacent genes and intergenic regions.

Conclusions

The results suggest that effects of assembled Pol III complexes on transcription of neighboring Pol II promoters are of greater magnitude than previously appreciated, that such effects influence expression of adjacent genes at transcriptional start site and translational levels, and may explain a function of the conserved ETC sites in yeast. The results may also be relevant to synthetic biology efforts to design a minimal yeast genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-014-0089-x) contains supplementary material, which is available to authorized users.

CMG helicase and DNA polymerase ε form a functional 15-subunit holoenzyme for eukaryotic leading-strand DNA replication

DNA replication in eukaryotes is asymmetric, with separate DNA polymerases (Pol) dedicated to bulk synthesis of the leading and lagging strands. Pol α/primase initiates primers on both strands that are extended by Pol ε on the leading strand and by Pol δ on the lagging strand. The CMG (Cdc45-MCM-GINS) helicase surrounds the leading strand and is proposed to recruit Pol ε for leading-strand synthesis, but to date a direct interaction between CMG and Pol ε has not been demonstrated. While purifying CMG helicase overexpressed in yeast, we detected a functional complex between CMG and native Pol ε. Using pure CMG and Pol ε, we reconstituted a stable 15-subunit CMG-Pol ε complex and showed that it is a functional polymerase-helicase on a model replication fork in vitro. On its own, the Pol2 catalytic subunit of Pol ε is inefficient in CMG-dependent replication, but addition of the Dpb2 protein subunit of Pol ε, known to bind the Psf1 protein subunit of CMG, allows stable synthesis with CMG. Dpb2 does not affect Pol δ function with CMG, and thus we propose that the connection between Dpb2 and CMG helps to stabilize Pol ε on the leading strand as part of a 15-subunit leading-strand holoenzyme we refer to as CMGE. Direct binding between Pol ε and CMG provides an explanation for specific targeting of Pol ε to the leading strand and provides clear mechanistic evidence for how strand asymmetry is maintained in eukaryotes.

Checkpoint-independent scaling of the Saccharomyces cerevisiae DNA replication program

Background

In budding yeast, perturbations that prolong S phase lead to a proportionate delay in the activation times of most origins. The DNA replication checkpoint was implicated in this scaling phenotype, as an intact checkpoint was shown to be required for the delayed activation of late origins in response to hydroxyurea treatment. In support of that, scaling is lost in cells deleted of mrc1, a mediator of the replication checkpoint signal. Mrc1p, however, also plays a role in normal replication.

Results

To examine whether the replication checkpoint is required for scaling the replication profile with S phase duration we measured the genome-wide replication profile of different MRC1 alleles that separate its checkpoint function from its role in normal replication, and further analyzed the replication profiles of S phase mutants that are checkpoint deficient. We found that the checkpoint is not required for scaling; rather the unique replication phenotype of mrc1 deleted cells is attributed to the role of Mrc1 in normal replication. This is further supported by the replication profiles of tof1Δ which functions together with Mrc1p in normal replication, and by the distinct replication profiles of specific POL2 alleles which differ in their interaction with Mrc1p.

Conclusions

We suggest that the slow fork progression in mrc1 deleted cells reduces the likelihood of passive replication leading to the activation of origins that remain mostly dormant in wild-type cells.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-014-0079-z) contains supplementary material, which is available to authorized users.

Mechanism of asymmetric polymerase assembly at the eukaryotic replication fork

Eukaryotes use distinct polymerases for leading- and lagging-strand replication, but how they target their respective strands is uncertain. We reconstituted Saccharomyces cerevisiae replication forks and found that CMG helicase selects polymerase (Pol) ɛ to the exclusion of Pol δ on the leading strand. Even if Pol δ assembles on the leading strand, Pol ɛ rapidly replaces it. Pol δ-PCNA is distributive with CMG, in contrast to its high stability on primed ssDNA. Hence CMG will not stabilize Pol δ, instead leaving the leading strand accessible for Pol ɛ and stabilizing Pol ɛ. Comparison of Pol ɛ and Pol δ on a lagging-strand model DNA reveals the opposite. Pol δ dominates over excess Pol ɛ on PCNA-primed ssDNA. Thus, PCNA strongly favors Pol δ over Pol ɛ on the lagging strand, but CMG over-rides and flips this balance in favor of Pol ɛ on the leading strand.

The precarious prokaryotic chromosome

Evolutionary selection for optimal genome preservation, replication, and expression should yield similar chromosome organizations in any type of cells. And yet, the chromosome organization is surprisingly different between eukaryotes and prokaryotes. The nuclear versus cytoplasmic accommodation of genetic material accounts for the distinct eukaryotic and prokaryotic modes of genome evolution, but it falls short of explaining the differences in the chromosome organization. I propose that the two distinct ways to organize chromosomes are driven by the differences between the global-consecutive chromosome cycle of eukaryotes and the local-concurrent chromosome cycle of prokaryotes. Specifically, progressive chromosome segregation in prokaryotes demands a single duplicon per chromosome, while other "precarious" features of the prokaryotic chromosomes can be viewed as compensations for this severe restriction.

Inferring the spatiotemporal DNA replication program from noisy data

We generalize a stochastic model of DNA replication to the case where replication-origin-initiation rates vary locally along the genome and with time. Using this generalized model, we address the inverse problem of inferring initiation rates from experimental data concerning replication in cell populations. Previous work based on curve fitting depended on arbitrarily chosen functional forms for the initiation rate, with free parameters that were constrained by the data. We introduce a nonparametric method of inference that is based on Gaussian process regression. The method replaces specific assumptions about the functional form of the initiation rate with more general prior expectations about the smoothness of variation of this rate, along the genome and in time. Using this inference method, we recover, with high precision, simulated replication schemes from noisy data that are typical of current experiments.

The dynamics of genome replication using deep sequencing

Eukaryotic genomes are replicated from multiple DNA replication origins. We present complementary deep sequencing approaches to measure origin location and activity in Saccharomyces cerevisiae. Measuring the increase in DNA copy number during a synchronous S-phase allowed the precise determination of genome replication. To map origin locations, replication forks were stalled close to their initiation sites; therefore, copy number enrichment was limited to origins. Replication timing profiles were generated from asynchronous cultures using fluorescence-activated cell sorting. Applying this technique we show that the replication profiles of haploid and diploid cells are indistinguishable, indicating that both cell types use the same cohort of origins with the same activities. Finally, increasing sequencing depth allowed the direct measure of replication dynamics from an exponentially growing culture. This is the first time this approach, called marker frequency analysis, has been successfully applied to a eukaryote. These data provide a high-resolution resource and methodological framework for studying genome biology.

Intergenic transcriptional interference is blocked by RNA polymerase III transcription factor TFIIIB in Saccharomyces cerevisiae

The major function of eukaryotic RNA polymerase III is to transcribe transfer RNA, 5S ribosomal RNA, and other small non-protein-coding RNA molecules. Assembly of the RNA polymerase III complex on chromosomal DNA requires the sequential binding of transcription factor complexes TFIIIC and TFIIIB. Recent evidence has suggested that in addition to producing RNA transcripts, chromatin-assembled RNA polymerase III complexes may mediate additional nuclear functions that include chromatin boundary, nucleosome phasing, and general genome organization activities. This study provides evidence of another such “extratranscriptional” activity of assembled RNA polymerase III complexes, which is the ability to block progression of intergenic RNA polymerase II transcription. We demonstrate that the RNA polymerase III complex bound to the tRNA gene upstream of the Saccharomyces cerevisiae ATG31 gene protects the ATG31 promoter against readthrough transcriptional interference from the upstream noncoding intergenic SUT467 transcription unit. This protection is predominately mediated by binding of the TFIIIB complex. When TFIIIB binding to this tRNA gene is weakened, an extended SUT467–ATG31 readthrough transcript is produced, resulting in compromised ATG31 translation. Since the ATG31 gene product is required for autophagy, strains expressing the readthrough transcript exhibit defective autophagy induction and reduced fitness under autophagy-inducing nitrogen starvation conditions. Given the recent discovery of widespread pervasive transcription in all forms of life, protection of neighboring genes from intergenic transcriptional interference may be a key extratranscriptional function of assembled RNA polymerase III complexes and possibly other DNA binding proteins.

High-resolution replication profiles define the stochastic nature of genome replication initiation and termination

Eukaryotic genome replication is stochastic, and each cell uses a different cohort of replication origins. We demonstrate that interpreting high-resolution Saccharomyces cerevisiae genome replication data with a mathematical model allows quantification of the stochastic nature of genome replication, including the efficiency of each origin and the distribution of termination events. Single-cell measurements support the inferred values for stochastic origin activation time. A strain, in which three origins were inactivated, confirmed that the distribution of termination events is primarily dictated by the stochastic activation time of origins. Cell-to-cell variability in origin activity ensures that termination events are widely distributed across virtually the whole genome. We propose that the heterogeneity in origin usage contributes to genome stability by limiting potentially deleterious events from accumulating at particular loci.

•

Deep sequencing reveals a high-resolution view of genome replication dynamics

•

Genome-wide modeling and single-cell imaging reveal stochastic origin activity

•

Origin activity determines the location of replication termination events

•

Termination events are widely distributed across the whole genome

Stochastic association of neighboring replicons creates replication factories in budding yeast

Inside the nucleus, DNA replication is organized at discrete sites called replication factories, consisting of DNA polymerases and other replication proteins. Replication factories play important roles in coordinating replication and in responding to replication stress. However, it remains unknown how replicons are organized for processing at each replication factory. Here we address this question using budding yeast. We analyze how individual replicons dynamically organized a replication factory using live-cell imaging and investigate how replication factories were structured using super-resolution microscopy. Surprisingly, we show that the grouping of replicons within factories is highly variable from cell to cell. Once associated, however, replicons stay together relatively stably to maintain replication factories. We derive a coherent genome-wide mathematical model showing how neighboring replicons became associated stochastically to form replication factories, which was validated by independent microscopy-based analyses. This study not only reveals the fundamental principles promoting replication factory organization in budding yeast, but also provides insight into general mechanisms by which chromosomes organize sub-nuclear structures.

DNA replication timing

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

Cost of rNTP/dNTP pool imbalance at the replication fork

The concentration of ribonucleoside triphosphates (rNTPs) in cells is far greater than the concentration of deoxyribonucleoside triphosphates (dNTPs), and this pool imbalance presents a challenge for DNA polymerases (Pols) to select their proper substrate. This report examines the effect of nucleotide pool imbalance on the rate and fidelity of the Escherichia coli replisome. We find that rNTPs decrease replication fork rate by competing with dNTPs at the active site of the C-family Pol III replicase at a step that does not require correct base-pairing. The effect of rNTPs on Pol rate generalizes to B-family eukaryotic replicases, Pols δ and ε. Imbalance of the dNTP pool also slows the replisome and thus is not specific to rNTPs. We observe a measurable frequency of rNMP incorporation that predicts one rNTP incorporated every 2.3 kb during chromosome replication. Given the frequency of rNMP incorporation, the repair of rNMPs is likely rapid. RNase HII nicks DNA at single rNMP residues to initiate replacement with dNMP. Considering that rNMPs will mark the new strand, RNase HII may direct strand-specificity for mismatch repair (MMR). How the newly synthesized strand is recognized for MMR is uncertain in eukaryotes and most bacteria, which lack a methyl-directed nicking system. Here we demonstrate that Bacillus subtilis incorporates rNMPs in vivo, that RNase HII plays a role in their removal, and the RNase HII gene deletion enhances mutagenesis, suggesting a possible role of incorporated rNMPs in MMR.

A variable fork rate affects timing of origin firing and S phase dynamics in Saccharomyces cerevisiae

Activation (in the following referred to as firing) of replication origins is a continuous and irreversible process regulated by availability of DNA replication molecules and cyclin-dependent kinase activities, which are often altered in human cancers. The temporal, progressive origin firing throughout S phase appears as a characteristic replication profile, and computational models have been developed to describe this process. Although evidence from yeast to human indicates that a range of replication fork rates is observed experimentally in order to complete a timely S phase, those models incorporate velocities that are uniform across the genome. Taking advantage of the availability of replication profiles, chromosomal position and replication timing, here we investigated how fork rate may affect origin firing in budding yeast. Our analysis suggested that patterns of origin firing can be observed from a modulation of the fork rate that strongly correlates with origin density. Replication profiles of chromosomes with a low origin density were fitted with a variable fork rate, whereas for the ones with a high origin density a constant fork rate was appropriate. This indeed supports the previously reported correlation between inter-origin distance and fork rate changes. Intriguingly, the calculated correlation between fork rate and timing of origin firing allowed the estimation of firing efficiencies for the replication origins. This approach correctly retrieved origin efficiencies previously determined for chromosome VI and provided testable prediction for other chromosomal origins. Our results gain deeper insights into the temporal coordination of genome duplication, indicating that control of the replication fork rate is required for the timely origin firing during S phase.

Quantitative, genome-wide analysis of eukaryotic replication initiation and termination

Many fundamental aspects of DNA replication, such as the exact locations where DNA synthesis is initiated and terminated, how frequently origins are used, and how fork progression is influenced by transcription, are poorly understood. Via the deep sequencing of Okazaki fragments, we comprehensively document replication fork directionality throughout the S. cerevisiae genome, which permits the systematic analysis of initiation, origin efficiency, fork progression, and termination. We show that leading-strand initiation preferentially occurs within a nucleosome-free region at replication origins. Using a strain in which late origins can be induced to fire early, we show that replication termination is a largely passive phenomenon that does not rely on cis-acting sequences or replication fork pausing. The replication profile is predominantly determined by the kinetics of origin firing, allowing us to reconstruct chromosome-wide timing profiles from an asynchronous culture.

ATR-like kinase Mec1 facilitates both chromatin accessibility at DNA replication forks and replication fork progression during replication stress

Faithful DNA replication is essential for normal cell division and differentiation. In eukaryotic cells, DNA replication takes place on chromatin. This poses the critical question as to how DNA replication can progress through chromatin, which is inhibitory to all DNA-dependent processes. Here, we developed a novel genome-wide method to measure chromatin accessibility to micrococcal nuclease (MNase) that is normalized for nucleosome density, the NCAM (normalized chromatin accessibility to MNase) assay. This method enabled us to discover that chromatin accessibility increases specifically at and ahead of DNA replication forks in normal S phase and during replication stress. We further found that Mec1, a key regulatory ATR-like kinase in the S-phase checkpoint, is required for both normal chromatin accessibility around replication forks and replication fork rate during replication stress, revealing novel functions for the kinase in replication stress response. These results suggest a possibility that Mec1 may facilitate DNA replication fork progression during replication stress by increasing chromatin accessibility around replication forks.

The spatiotemporal program of replication in the genome of Lachancea kluyveri

We generated a genome-wide replication profile in the genome of Lachancea kluyveri and assessed the relationship between replication and base composition. This species diverged from Saccharomyces cerevisiae before the ancestral whole genome duplication. The genome comprises eight chromosomes among which a chromosomal arm of 1 Mb has a G + C-content much higher than the rest of the genome. We identified 252 active replication origins in L. kluyveri and found considerable divergence in origin location with S. cerevisiae and with Lachancea waltii. Although some global features of S. cerevisiae replication are conserved: Centromeres replicate early, whereas telomeres replicate late, we found that replication origins both in L. kluyveri and L. waltii do not behave as evolutionary fragile sites. In L. kluyveri, replication timing along chromosomes alternates between regions of early and late activating origins, except for the 1 Mb GC-rich chromosomal arm. This chromosomal arm contains an origin consensus motif different from other chromosomes and is replicated early during S-phase. We showed that precocious replication results from the specific absence of late firing origins in this chromosomal arm. In addition, we found a correlation between GC-content and distance from replication origins as well as a lack of replication-associated compositional skew between leading and lagging strands specifically in this GC-rich chromosomal arm. These findings suggest that the unusual base composition in the genome of L. kluyveri could be linked to replication.

Local and global functions of Timeless and Tipin in replication fork protection

The eukaryotic cell replicates its chromosomal DNA with almost absolute fidelity in the course of every cell cycle. This accomplishment is remarkable considering that the conditions for DNA replication are rarely ideal. The replication machinery encounters a variety of obstacles on the chromosome, including damaged template DNA. In addition, a number of chromosome regions are considered to be difficult to replicate owing to DNA secondary structures and DNA binding proteins required for various transactions on the chromosome. Under these conditions, replication forks stall or break, posing grave threats to genomic integrity. How does the cell combat such stressful conditions during DNA replication? The replication fork protection complex (FPC) may help answer this question. Recent studies have demonstrated that the FPC is required for the smooth passage of replication forks at difficult-to-replicate genomic regions and plays a critical role in coordinating multiple genome maintenance processes at the replication fork.

Mathematical modeling of genome replication

Eukaryotic DNA replication is initiated from multiple sites on the chromosome, but little is known about the global and local regulation of replication. We present a mathematical model for the spatial dynamics of DNA replication, which offers insight into the kinetics of replication in different types of organisms. Most biological experiments involve average quantities over large cell populations (typically >10(7) cells) and therefore can mask the cell-to-cell variability present in the system. Although the model is formulated in terms of a population of cells, using mathematical analysis we show that one can obtain signatures of stochasticity in individual cells from averaged quantities. This work generalizes the result by Retkute et al. [Phys. Rev. Lett. 107, 068103 (2011)] to a broader set of parameter regimes.

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.

Physical tethering and volume exclusion determine higher-order genome organization in budding yeast

In this paper we show that tethering of heterochromatic regions to nuclear landmarks and random encounters of chromosomes in the confined nuclear volume are sufficient to explain the higher-order organization of the budding yeast genome. We have quantitatively characterized the contact patterns and nuclear territories that emerge when chromosomes are allowed to behave as constrained but otherwise randomly configured flexible polymer chains in the nucleus. Remarkably, this constrained random encounter model explains in a statistical manner the experimental hallmarks of the S. cerevisiae genome organization, including (1) the folding patterns of individual chromosomes; (2) the highly enriched interactions between specific chromatin regions and chromosomes; (3) the emergence, shape, and position of gene territories; (4) the mean distances between pairs of telomeres; and (5) even the co-location of functionally related gene loci, including early replication start sites and tRNA genes. Therefore, most aspects of the yeast genome organization can be explained without calling on biochemically mediated chromatin interactions. Such interactions may modulate the pre-existing propensity for co-localization but seem not to be the cause for the observed higher-order organization. The fact that geometrical constraints alone yield a highly organized genome structure, on which different functional elements are specifically distributed, has strong implications for the folding principles of the genome and the evolution of its function.

Genome-wide localization of replication factors

Chromatin Immunoprecipitation (ChIP) is a powerful tool for the identification and characterization of protein-DNA interactions in vivo. ChIP has been utilized to study diverse nuclear processes such as transcription regulation, chromatin modification, DNA recombination and DNA replication at specific loci and, more recently, across the entire genome. Advances in genomic approaches, and whole genome sequencing in particular, have made it possible and affordable to comprehensively identify specific protein binding sites throughout the genomes of higher eukaryotes. The dynamic nature of the DNA replication program and the transient occupancy of many replication factors throughout the cell cycle present additional challenges that may not pertain to the mapping of site specific transcription factors. Here we discuss the specific considerations that need to be addressed in the application of ChIP to the genome-wide location analysis of protein factors involved in DNA replication.

Optimal placement of origins for DNA replication

DNA replication is an essential process in biology and its timing must be robust so that cells can divide properly. Random fluctuations in the formation of replication starting points, called origins, and the subsequent activation of proteins lead to variations in the replication time. We analyze these stochastic properties of DNA and derive the positions of origins corresponding to the minimum replication time. We show that under some conditions the minimization of replication time leads to the grouping of origins, and relate this to experimental data in a number of species showing origin grouping.

The GINS complex: structure and function

Eukaryotic chromosomal DNA replication is controlled by a highly ordered series of steps involving multiple proteins at replication origins. The eukaryotic GINS complex is essential for the establishment of DNA replication forks and replisome progression. GINS is one of the core components of the eukaryotic replicative helicase, the CMG (Cdc45-MCM-GINS) complex, which unwinds duplex DNA ahead of the moving replication fork. Eukaryotic GINS also links with other key proteins at the fork to maintain an active replisome progression complex. Archaeal GINS homologues play a central role in chromosome replication by associating with other replisome components. This chapter focuses on the molecular events related with DNA replication initiation, and summarizes our current understanding of the function, structure and evolution of the GINS complex in eukaryotes and archaea.

OriDB, the DNA replication origin database updated and extended

OriDB (http://www.oridb.org/) is a database containing collated genome-wide mapping studies of confirmed and predicted replication origin sites. The original database collated and curated Saccharomyces cerevisiae origin mapping studies. Here, we report that the OriDB database and web site have been revamped to improve user accessibility to curated data sets, to greatly increase the number of curated origin mapping studies, and to include the collation of replication origin sites in the fission yeast Schizosaccharomyces pombe. The revised database structure underlies these improvements and will facilitate further expansion in the future. The updated OriDB for S. cerevisiae is available at http://cerevisiae.oridb.org/ and for S. pombe at http://pombe.oridb.org/.

Do replication forks control late origin firing in Saccharomyces cerevisiae?

Recent studies of eukaryotic DNA replication timing profiles suggest that the time-dependent rate of origin firing, I(t), has a universal shape, which ensures a reproducible replication completion time. However, measurements of I(t) are based on population averages, which may bias the shape of the I(t) because of imperfect cell synchrony and cell-to-cell variability. Here, we measure the population-averaged I(t) profile from synchronized Saccharomyces cerevisiae cells using DNA combing and we extract the single-cell I(t) profile using numerical deconvolution. The single cell I(t) and the population-averaged I(t) extracted from DNA combing and replication timing profiles are similar, indicating a genome scale invariance of the replication process, and excluding cell-to-cell variability in replication time as an explanation for the shape of I(t). The single cell I(t) correlates with fork density in wild-type cells, which is specifically loosened in late S phase in the clb5Δ mutant. A previously proposed numerical model that reproduces the wild-type I(t) profile, could also describe the clb5Δ mutant I(t) once modified to incorporate the decline in CDK activity and the looser dependency of initiation on fork density in the absence of Clb5p. Overall, these results suggest that the replication forks emanating from early fired origins facilitate origin firing in later-replicating regions.

Extra-transcriptional functions of RNA Polymerase III complexes: TFIIIC as a potential global chromatin bookmark

RNA polymerase III (Pol III) is one of three eukaryotic transcription complexes, and was identified as the complex responsible for production of transfer RNA and a limited number of other small RNAs. Pol III transcription at tRNA genes (tDNAs) requires the binding of two transcription factor complexes, TFIIIC and TFIIIB. Recent evidence points to a larger role for the Pol III transcription system in various other nuclear processes, including effects on nucleosome positioning, global genome and sub-nuclear organization, and direct effects on RNA polymerase II (Pol II) transcription. These effects are perhaps mediated by recruitment of a host of other chromatin proteins, including Pol II transcription factors and chromatin enzymes. Extra-TFIIIC sites (ETC sites) are chromosomal locations bound by TFIIIC without the rest of the Pol III complex, and bound TFIIIC alone is also able to mediate additional functions. These so called "extra-transcriptional effects" of the Pol III system are reviewed here, and a model is put forth suggesting that the TFIIIC transcription factor may act as a stably bound, global "bookmark" within chromatin to establish, maintain, or demarcate chromatin states as cells divide or change gene expression patterns.

Dynamics of DNA replication in yeast

We present a mathematical model for the spatial dynamics of DNA replication. Using this model we determine the probability distribution for the time at which each chromosomal position is replicated. From this we show, contrary to previous reports, that mean replication time curves cannot be used to directly determine origin parameters. We demonstrate that the stochastic nature of replication dynamics leaves a clear signature in experimentally measured population average data, and we show that the width of the activation time probability distribution can be inferred from this data. Our results compare favorably with experimental measurements in Saccharomyces cerevisae.

Mutation rates across budding yeast chromosome VI are correlated with replication timing

Previous experimental studies suggest that the mutation rate is nonuniform across the yeast genome. To characterize this variation across the genome more precisely, we measured the mutation rate of the URA3 gene integrated at 43 different locations tiled across Chromosome VI. We show that mutation rate varies 6-fold across a single chromosome, that this variation is correlated with replication timing, and we propose a model to explain this variation that relies on the temporal separation of two processes for replicating past damaged DNA: error-free DNA damage tolerance and translesion synthesis. This model is supported by the observation that eliminating translesion synthesis decreases this variation.

The torsional state of DNA within the chromosome

Virtually all processes of the genome biology affect or are affected by the torsional state of DNA. Torsional energy associated with an altered twist facilitates or hinders the melting of the double helix, its molecular interactions, and its spatial folding in the form of supercoils. Yet, understanding how the torsional state of DNA is modulated remains a challenging task due to the multiplicity of cellular factors involved in the generation, transmission, and dissipation of DNA twisting forces. Here, an overview of the implication of DNA topoisomerases, DNA revolving motors, and other DNA interactions that determine local levels of torsional stress in bacterial and eukaryotic chromosomes is provided. Particular emphasis is made on the experimental approaches being developed to assess the torsional state of intracellular DNA and its organization into topological domains.

Genome stability control by checkpoint regulation of tRNA gene transcription

The RNA polymerase III pre-initiation complex (PIC) assembled on yeast tRNA genes naturally causes replication fork pausing that contributes to genome instability. Mechanistic coupling of the fork pausing activity of tRNA genes to replication has long been considered likely, but only recently demonstrated. In contrast to the expectation that this coupling might occur by a passive mechanism such as direct disruption of transcription factor-DNA complexes by a component of the replisome, it turns out that disassembly of the RNA polymerase III PIC is actively controlled by the replication stress checkpoint signal transduction pathway. This advance supports a new model in which checkpoint-dependent disassembly of the transcription machinery at tRNA genes is a vital component of an overall system of genome stability control that also targets replication and DNA repair proteins.

Modeling genome-wide replication kinetics reveals a mechanism for regulation of replication timing

Microarrays are powerful tools to probe genome-wide replication kinetics. The rich data sets that result contain more information than has been extracted by current methods of analysis. In this paper, we present an analytical model that incorporates probabilistic initiation of origins and passive replication. Using the model, we performed least-squares fits to a set of recently published time course microarray data on Saccharomyces cerevisiae. We extracted the distribution of firing times for each origin and found that the later an origin fires on average, the greater the variation in firing times. To explain this trend, we propose a model where earlier-firing origins have more initiator complexes loaded and a more accessible chromatin environment. The model demonstrates how initiation can be stochastic and yet occur at defined times during S phase, without an explicit timing program. Furthermore, we hypothesize that the initiators in this model correspond to loaded minichromosome maintenance complexes. This model is the first to suggest a detailed, testable, biochemically plausible mechanism for the regulation of replication timing in eukaryotes.