A genome-wide comprehensive analysis of nucleosome positioning in yeast.

In eukaryotic cells, the one-dimensional DNA molecules need to be tightly packaged into the spatially constraining nucleus. Folding is achieved on its lowest level by wrapping the DNA around nucleosomes. Their arrangement regulates other nuclear processes, such as transcription and DNA repair. Despite strong efforts to study nucleosome positioning using Next Generation Sequencing (NGS) data, the mechanism of their collective arrangement along the gene body remains poorly understood. Here, we classify nucleosome distributions of protein-coding genes in Saccharomyces cerevisiae according to their profile similarity and analyse their differences using functional Principal Component Analysis. By decomposing the NGS signals into their main descriptive functions, we compared wild type and chromatin remodeler-deficient strains, keeping position-specific details preserved whilst considering the nucleosome arrangement as a whole. A correlation analysis with other genomic properties, such as gene size and length of the upstream Nucleosome Depleted Region (NDR), identified key factors that influence the nucleosome distribution. We reveal that the RSC chromatin remodeler-which is responsible for NDR maintenance-is indispensable for decoupling nucleosome arrangement within the gene from positioning outside, which interfere in rsc8-depleted conditions. Moreover, nucleosome profiles in chd1Δ strains displayed a clear correlation with RNA polymerase II presence, whereas wild type cells did not indicate a noticeable interdependence. We propose that RSC is pivotal for global nucleosome organisation, whilst Chd1 plays a key role for maintaining local arrangement.

Rif1 restrains the rate of replication origin firing in Xenopus laevis.

Metazoan genomes are duplicated by the coordinated activation of clusters of replication origins at different times during S phase, but the underlying mechanisms of this temporal program remain unclear during early development. Rif1, a key replication timing factor, inhibits origin firing by recruiting protein phosphatase 1 (PP1) to chromatin counteracting S phase kinases. We have previously described that Rif1 depletion accelerates early Xenopus laevis embryonic cell cycles. Here, we find that in the absence of Rif1, patterns of replication foci change along with the acceleration of replication cluster activation. However, initiations increase only moderately inside active clusters. Our numerical simulations suggest that the absence of Rif1 compresses the temporal program towards more homogeneity and increases the availability of limiting initiation factors. We experimentally demonstrate that Rif1 depletion increases the chromatin-binding of the S phase kinase Cdc7/Drf1, the firing factors Treslin, MTBP, Cdc45, RecQL4, and the phosphorylation of both Treslin and MTBP. We show that Rif1 globally, but not locally, restrains the replication program in early embryos, possibly by inhibiting or excluding replication factors from chromatin.

Neural network and kinetic modelling of human genome replication reveal replication origin locations and strengths.

In human and other metazoans, the determinants of replication origin location and strength are still elusive. Origins are licensed in G1 phase and fired in S phase of the cell cycle, respectively. It is debated which of these two temporally separate steps determines origin efficiency. Experiments can independently profile mean replication timing (MRT) and replication fork directionality (RFD) genome-wide. Such profiles contain information on multiple origins' properties and on fork speed. Due to possible origin inactivation by passive replication, however, observed and intrinsic origin efficiencies can markedly differ. Thus, there is a need for methods to infer intrinsic from observed origin efficiency, which is context-dependent. Here, we show that MRT and RFD data are highly consistent with each other but contain information at different spatial scales. Using neural networks, we infer an origin licensing landscape that, when inserted in an appropriate simulation framework, jointly predicts MRT and RFD data with unprecedented precision and underlies the importance of dispersive origin firing. We furthermore uncover an analytical formula that predicts intrinsic from observed origin efficiency combined with MRT data. Comparison of inferred intrinsic origin efficiencies with experimental profiles of licensed origins (ORC, MCM) and actual initiation events (Bubble-seq, SNS-seq, OK-seq, ORM) show that intrinsic origin efficiency is not solely determined by licensing efficiency. Thus, human replication origin efficiency is set at both the origin licensing and firing steps.

Functional interplay between Mediator and RSC chromatin remodeling complex controls nucleosome-depleted region maintenance at promoters.

Chromatin organization is crucial for transcriptional regulation in eukaryotes. Mediator is an essential and conserved co-activator thought to act in concert with chromatin regulators. However, it remains largely unknown how their functions are coordinated. Here, we provide evidence in the yeast Saccharomyces cerevisiae that Mediator establishes physical contact with RSC (Remodels the Structure of Chromatin), a conserved and essential chromatin remodeling complex that is crucial for nucleosome-depleted region (NDR) formation. We determine the role of Mediator-RSC interaction in their chromatin binding, nucleosome occupancy, and transcription on a genomic scale. Mediator and RSC co-localize on wide NDRs of promoter regions, and specific Mediator mutations affect nucleosome eviction and TSS-associated +1 nucleosome stability. This work shows that Mediator contributes to RSC remodeling function to shape NDRs and maintain chromatin organization on promoter regions. It will help in our understanding of transcriptional regulation in the chromatin context relevant for severe diseases.

A quantitative modelling approach for DNA repair on a population scale.

The great advances of sequencing technologies allow the in vivo measurement of nuclear processes-such as DNA repair after UV exposure-over entire cell populations. However, data sets usually contain only a few samples over several hours, missing possibly important information in between time points. We developed a data-driven approach to analyse CPD repair kinetics over time in Saccharomyces cerevisiae. In contrast to other studies that consider sequencing signals as an average behaviour, we understand them as the superposition of signals from independent cells. By motivating repair as a stochastic process, we derive a minimal model for which the parameters can be conveniently estimated. We correlate repair parameters to a variety of genomic features that are assumed to influence repair, including transcription rate and nucleosome density. The clearest link was found for the transcription unit length, which has been unreported for budding yeast to our knowledge. The framework hence allows a comprehensive analysis of nuclear processes on a population scale.

Genomic analysis of Rad26 and Rad1-Rad10 reveals differences in their dependence on Mediator and RNA polymerase II.

Mediator is a conserved coregulator playing a key role in RNA polymerase (Pol) II transcription. Mediator also links transcription and nucleotide excision repair (NER) via a direct contact with Rad2/ERCC5(XPG) endonuclease. In this work, we analyzed the genome-wide distribution of Rad26/ERCC6(CSB) and Rad1-Rad10/ERCC4(XPF)-ERCC1, addressing the question of a potential link of these proteins with Mediator and Pol II in yeast Saccharomyces cerevisiae Our genomic analyses reveal that Rad1-Rad10 and Rad26 are present on the yeast genome in the absence of genotoxic stress, especially at highly transcribed regions, with Rad26 binding strongly correlating with that of Pol II. Moreover, we show that Rad1-Rad10 and Rad26 colocalize with Mediator at intergenic regions and physically interact with this complex. Using kin28 TFIIH mutant, we found that Mediator stabilization on core promoters leads to an increase in Rad1-Rad10 chromatin binding, whereas Rad26 occupancy follows mainly a decrease in Pol II transcription. Combined with multivariate analyses, our results show the relationships between Rad1-Rad10, Rad26, Mediator, and Pol II, modulated by the changes in binding dynamics of Mediator and Pol II transcription. In conclusion, we extend the Mediator link to Rad1-Rad10 and Rad26 NER proteins and reveal important differences in their dependence on Mediator and Pol II. Rad2 is the most dependent on Mediator, followed by Rad1-Rad10, whereas Rad26 is the most closely related to Pol II. Our work thus contributes to new concepts of the functional interplay between transcription and DNA repair machineries, which are relevant for human diseases including cancer and XP/CS syndromes.

Polo-like kinase 1 (Plk1) regulates DNA replication origin firing and interacts with Rif1 in Xenopus.

The activation of eukaryotic DNA replication origins needs to be strictly controlled at multiple steps in order to faithfully duplicate the genome and to maintain its stability. How the checkpoint recovery and adaptation protein Polo-like kinase 1 (Plk1) regulates the firing of replication origins during non-challenged S phase remained an open question. Using DNA fiber analysis, we show that immunodepletion of Plk1 in the Xenopus in vitro system decreases replication fork density and initiation frequency. Numerical analyses suggest that Plk1 reduces the overall probability and synchrony of origin firing. We used quantitative chromatin proteomics and co-immunoprecipitations to demonstrate that Plk1 interacts with firing factors MTBP/Treslin/TopBP1 as well as with Rif1, a known regulator of replication timing. Phosphopeptide analysis by LC/MS/MS shows that the C-terminal domain of Rif1, which is necessary for its repressive action on origins through protein phosphatase 1 (PP1), can be phosphorylated in vitro by Plk1 on S2058 in its PP1 binding site. The phosphomimetic S2058D mutant interrupts the Rif1-PP1 interaction and modulates DNA replication. Collectively, our study provides molecular insights into how Plk1 regulates the spatio-temporal replication program and suggests that Plk1 controls origin activation at the level of large chromatin domains in vertebrates.

Organization of DNA Replication Origin Firing in Xenopus Egg Extracts: The Role of Intra-S Checkpoint.

During cell division, the duplication of the genome starts at multiple positions called replication origins. Origin firing requires the interaction of rate-limiting factors with potential origins during the S(ynthesis)-phase of the cell cycle. Origins fire as synchronous clusters which is proposed to be regulated by the intra-S checkpoint. By modelling the unchallenged, the checkpoint-inhibited and the checkpoint protein Chk1 over-expressed replication pattern of single DNA molecules from Xenopus sperm chromatin replicated in egg extracts, we demonstrate that the quantitative modelling of data requires: (1) a segmentation of the genome into regions of low and high probability of origin firing; (2) that regions with high probability of origin firing escape intra-S checkpoint regulation and (3) the variability of the rate of DNA synthesis close to replication forks is a necessary ingredient that should be taken in to account in order to describe the dynamic of replication origin firing. This model implies that the observed origin clustering emerges from the apparent synchrony of origin firing in regions with high probability of origin firing and challenge the assumption that the intra-S checkpoint is the main regulator of origin clustering.

Polo-like kinase 1 (Plk1) is a positive regulator of DNA replication in the Xenopus in vitro system.

Polo-like kinase 1 (Plk1) is a cell cycle kinase essential for mitosis progression, but also important for checkpoint recovery and adaptation in response to DNA damage and replication stress. However, although Plk1 is expressed in S phase, little is known about its function during unperturbed DNA replication. Using Xenopus laevis egg extracts, mimicking early embryonic replication, we demonstrate that Plk1 is simultaneously recruited to chromatin with pre-replication proteins where it accumulates throughout S phase. Further, we found that chromatin-bound Plk1 is phosphorylated on its activating site T201, which appears to be sensitive to dephosphorylation by protein phosphatase 2A. Extracts immunodepleted of Plk1 showed a decrease in DNA replication, rescued by wild type recombinant Plk1. Inversely, modest Plk1 overexpression accelerated DNA replication. Plk1 depletion led to an increase in Chk1 phosphorylation and to a decrease in Cdk2 activity, which strongly suggests that Plk1 could inhibit the ATR/Chk1-dependent intra-S phase checkpoint during normal S phase. In addition, we observed that phosphorylated Plk1 levels are high during the rapid, early cell cycles of Xenopus development but decrease after the mid-blastula transition when the cell cycle and the replication program slow down along with more active checkpoints. These data shed new light on the role of Plk1 as a positive regulating factor for DNA replication in early, rapidly dividing embryos.

On the Interplay of the DNA Replication Program and the Intra-S Phase Checkpoint Pathway.

DNA replication in eukaryotes is achieved by the activation of multiple replication origins which needs to be precisely coordinated in space and time. This spatio-temporal replication program is regulated by many factors to maintain genome stability, which is frequently threatened through stresses of exogenous or endogenous origin. Intra-S phase checkpoints monitor the integrity of DNA synthesis and are activated when replication forks are stalled. Their activation leads to the stabilization of forks, to the delay of the replication program by the inhibition of late firing origins, and the delay of G2/M phase entry. In some cell cycles during early development these mechanisms are less efficient in order to allow rapid cell divisions. In this article, we will review our current knowledge of how the intra-S phase checkpoint regulates the replication program in budding yeast and metazoan models, including early embryos with rapid S phases. We sum up current models on how the checkpoint can inhibit origin firing in some genomic regions, but allow dormant origin activation in other regions. Finally, we discuss how numerical and theoretical models can be used to connect the multiple different actors into a global process and to extract general rules.

The eukaryotic bell-shaped temporal rate of DNA replication origin firing emanates from a balance between origin activation and passivation.

The time-dependent rate [Formula: see text] of origin firing per length of unreplicated DNA presents a universal bell shape in eukaryotes that has been interpreted as the result of a complex time-evolving interaction between origins and limiting firing factors. Here, we show that a normal diffusion of replication fork components towards localized potential replication origins (p-oris) can more simply account for the [Formula: see text] universal bell shape, as a consequence of a competition between the origin firing time and the time needed to replicate DNA separating two neighboring p-oris. We predict the [Formula: see text] maximal value to be the product of the replication fork speed with the squared p-ori density. We show that this relation is robustly observed in simulations and in experimental data for several eukaryotes. Our work underlines that fork-component recycling and potential origins localization are sufficient spatial ingredients to explain the universality of DNA replication kinetics.

Tight Chk1 Levels Control Replication Cluster Activation in Xenopus.

DNA replication in higher eukaryotes initiates at thousands of origins according to a spatio-temporal program. The ATR/Chk1 dependent replication checkpoint inhibits the activation of later firing origins. In the Xenopus in vitro system initiations are not sequence dependent and 2-5 origins are grouped in clusters that fire at different times despite a very short S phase. We have shown that the temporal program is stochastic at the level of single origins and replication clusters. It is unclear how the replication checkpoint inhibits late origins but permits origin activation in early clusters. Here, we analyze the role of Chk1 in the replication program in sperm nuclei replicating in Xenopus egg extracts by a combination of experimental and modelling approaches. After Chk1 inhibition or immunodepletion, we observed an increase of the replication extent and fork density in the presence or absence of external stress. However, overexpression of Chk1 in the absence of external replication stress inhibited DNA replication by decreasing fork densities due to lower Cdk2 kinase activity. Thus, Chk1 levels need to be tightly controlled in order to properly regulate the replication program even during normal S phase. DNA combing experiments showed that Chk1 inhibits origins outside, but not inside, already active clusters. Numerical simulations of initiation frequencies in the absence and presence of Chk1 activity are consistent with a global inhibition of origins by Chk1 at the level of clusters but need to be combined with a local repression of Chk1 action close to activated origins to fit our data.

From simple bacterial and archaeal replicons to replication N/U-domains.

The Replicon Theory proposed 50 years ago has proven to apply for replicons of the three domains of life. Here, we review our knowledge of genome organization into single and multiple replicons in bacteria, archaea and eukarya. Bacterial and archaeal replicator/initiator systems are quite specific and efficient, whereas eukaryotic replicons show degenerate specificity and efficiency, allowing for complex regulation of origin firing time. We expand on recent evidence that ~50% of the human genome is organized as ~1,500 megabase-sized replication domains with a characteristic parabolic (U-shaped) replication timing profile and linear (N-shaped) gradient of replication fork polarity. These N/U-domains correspond to self-interacting segments of the chromatin fiber bordered by open chromatin zones and replicate by cascades of origin firing initiating at their borders and propagating to their center, possibly by fork-stimulated initiation. The conserved occurrence of this replication pattern in the germline of mammals has resulted over evolutionary times in the formation of megabase-sized domains with an N-shaped nucleotide compositional skew profile due to replication-associated mutational asymmetries. Overall, these results reveal an evolutionarily conserved but developmentally plastic organization of replication that is driving mammalian genome evolution.

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.

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.

DNA replication induces compositional biases in yeast.

Asymmetries intrinsic to the process of DNA replication are expected to cause differences in the substitution patterns of the leading and the lagging strands and to induce compositional biases. These biases have been detected in the majority of eubacterial genomes but rarely in eukaryotes. Only in the human genome, the activity of a minority of replication origins seems to generate compositional biases. In this work, we provide evidence for replication-associated GC and TA skews in the genomes of two yeast species, Saccharomyces cerevisiae and Kluyveromyces lactis, whereas the data for the Schizosaccharomyces pombe genome are less conclusive. In contrast with the genomes of Homo sapiens and of the majority of eubacteria, the leading strand is enriched in cytosine and adenine in both S. cerevisiae and K. lactis. We observed significant variations across the interorigin intervals of several substitution rates in the S. cerevisiae lineage since its divergence from S. paradoxus. We also found that the S. cerevisiae genome is far from compositional equilibrium and that its present compositional biases are due to substitution rates operating before its divergence from S. paradoxus. Finally, we observed that replication and transcription tend to be cooriented in the S. cerevisiae genome, especially for genes encoding subunits of protein complexes. Taken together, our results suggest that replication-related compositional biases may be a feature of many eukaryotic genomes despite the stochastic nature of the firing of replication origins in these genomes.

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.

Giant yeast cells with nonrecyclable ribonucleotide reductase.

Ribonucleotide reductase (RNR) catalyzes the reduction of ribonucleotides to deoxyribonucleotides and thereby provides the precursors required for DNA synthesis and repair. In an attempt to test cell resistance to a permanent replicational stress, we constructed a mutant Saccharomyces cerevisiae strain containing exclusively nonrecyclable catalytic subunits of RNR that become inactivated following the reduction of one ribonucleoside diphosphate. In this rnr1C883A rnr3Δ mutant, the synthesis of each deoxyribonucleotide thus requires the production of one Rnr1C883A protein, which means that 26 million Rnr1C883A proteins (half the protein complement of a wild-type cell) have to be produced during each cell cycle. rnr1C883A rnr3Δ cells grow under constant replicational stress, as evidenced by the constitutive activation of the checkpoint effector Rad53, and their S phase is considerably extended compared to the wild type. rnr1C883A rnr3Δ mutants also display additional abnormalities such as a median cell volume increased by a factor of 8, and the presence of massive inclusion bodies. However, they exhibit a good plating efficiency and can be propagated indefinitely. rnr1C883A rnr3Δ cells, which can be used as a protein overexpression system, thus illustrate the robustness of S. cerevisiae to multiple physiological parameters.

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

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

Mathematical modelling of eukaryotic DNA replication.

Eukaryotic DNA replication is a complex process. Replication starts at thousand origins that are activated at different times in S phase and terminates when converging replication forks meet. Potential origins are much more abundant than actually fire within a given S phase. The choice of replication origins and their time of activation is never exactly the same in any two cells. Individual origins show different efficiencies and different firing time probability distributions, conferring stochasticity to the DNA replication process. High-throughput microarray and sequencing techniques are providing increasingly huge datasets on the population-averaged spatiotemporal patterns of DNA replication in several organisms. On the other hand, single-molecule replication mapping techniques such as DNA combing provide unique information about cell-to-cell variability in DNA replication patterns. Mathematical modelling is required to fully comprehend the complexity of the chromosome replication process and to correctly interpret these data. Mathematical analysis and computer simulations have been recently used to model and interpret genome-wide replication data in the yeast Saccharomyces cerevisiae and Schizosaccharomyces pombe, in Xenopus egg extracts and in mammalian cells. These works reveal how stochasticity in origin usage confers robustness and reliability to the DNA replication process.