MRC1-dependent scaling of the budding yeast DNA replication timing program

Real-time assessment of mitochondrial DNA heteroplasmy dynamics at the single-cell level

Mitochondrial DNA (mtDNA) is present in multiple copies within cells and is required for mitochondrial ATP generation. Even within individual cells, mtDNA copies can differ in their sequence, a state known as heteroplasmy. The principles underlying dynamic changes in the degree of heteroplasmy remain incompletely understood, due to the inability to monitor this phenomenon in real time. Here, we employ mtDNA-based fluorescent markers, microfluidics, and automated cell tracking, to follow mtDNA variants in live heteroplasmic yeast populations at the single-cell level. This approach, in combination with direct mtDNA tracking and data-driven mathematical modeling reveals asymmetric partitioning of mtDNA copies during cell division, as well as limited mitochondrial fusion and fission frequencies, as critical driving forces for mtDNA variant segregation. Given that our approach also facilitates assessment of segregation between intact and mutant mtDNA, we anticipate that it will be instrumental in elucidating the mechanisms underlying the purifying selection of mtDNA.

Identification of different classes of genome instability suppressor genes through analysis of DNA damage response markers

Cellular pathways that detect DNA damage are useful for identifying genes that suppress DNA damage, which can cause genome instability and cancer predisposition syndromes when mutated. We identified 199 high-confidence and 530 low-confidence DNA damage-suppressing (DDS) genes in Saccharomyces cerevisiae through a whole-genome screen for mutations inducing Hug1 expression, a focused screen for mutations inducing Ddc2 foci, and data from previous screens for mutations causing Rad52 foci accumulation and Rnr3 induction. We also identified 286 high-confidence and 394 low-confidence diverse genome instability-suppressing (DGIS) genes through a whole-genome screen for mutations resulting in increased gross chromosomal rearrangements and data from previous screens for mutations causing increased genome instability as assessed in a diversity of genome instability assays. Genes that suppress both pathways (DDS+ DGIS+) prevent or repair DNA replication damage and likely include genes preventing collisions between the replication and transcription machineries. DDS+ DGIS- genes, including many transcription-related genes, likely suppress damage that is normally repaired properly or prevent inappropriate signaling, whereas DDS- DGIS+ genes, like PIF1, do not suppress damage but likely promote its proper, nonmutagenic repair. Thus, induction of DNA damage markers is not a reliable indicator of increased genome instability, and the DDS and DGIS categories define mechanistically distinct groups of genes.

Loss of mitochondrial DNA is associated with reduced DNA content variability in Saccharomyces cerevisiae

DNA content measurement by fluorescence-assisted cell sorting (FACS) provides information on cell cycle progression and DNA content variability. Saccharomyces cerevisiae mutants with DNA content variability that was reduced relative to wild-type strains had defects in mitochondrial DNA (mtDNA) maintenance and mitochondrial gene expression and were correlated with strains found to lack mtDNA ([ rho 0 ] cells) by genome sequencing and fluorescence microscopy. In contrast, mutants with increased variability had defects in cell cycle progression, which may indicate a loss of coordination between mtDNA and nuclear DNA replication. Thus, FACS measurement of DNA content variability can provide insight into cell-to-cell heterogeneity in mtDNA copy number.

Context-dependent neocentromere activity in synthetic yeast chromosome VIII

Pioneering advances in genome engineering, and specifically in genome writing, have revolutionized the field of synthetic biology, propelling us toward the creation of synthetic genomes. The Sc2.0 project aims to build the first fully synthetic eukaryotic organism by assembling the genome of Saccharomyces cerevisiae. With the completion of synthetic chromosome VIII (synVIII) described here, this goal is within reach. In addition to writing the yeast genome, we sought to manipulate an essential functional element: the point centromere. By relocating the native centromere sequence to various positions along chromosome VIII, we discovered that the minimal 118-bp CEN8 sequence is insufficient for conferring chromosomal stability at ectopic locations. Expanding the transplanted sequence to include a small segment (∼500 bp) of the CDEIII-proximal pericentromere improved chromosome stability, demonstrating that minimal centromeres display context-dependent functionality.

Replisome-cohesin interactions provided by the Tof1-Csm3 and Mrc1 cohesion establishment factors

The chromosomal cohesin complex establishes sister chromatid cohesion during S phase, which forms the basis for faithful segregation of DNA replication products during cell divisions. Cohesion establishment is defective in the absence of either of three non-essential Saccharomyces cerevisiae replication fork components Tof1-Csm3 and Mrc1. Here, we investigate how these conserved factors contribute to cohesion establishment. Tof1-Csm3 and Mrc1 serve known roles during DNA replication, including replication checkpoint signaling, securing replication fork speed, as well as recruiting topoisomerase I and the histone chaperone FACT. By modulating each of these functions independently, we rule out that one of these known replication roles explains the contribution of Tof1-Csm3 and Mrc1 to cohesion establishment. Instead, using purified components, we reveal direct and multipronged protein interactions of Tof1-Csm3 and Mrc1 with the cohesin complex. Our findings open the possibility that a series of physical interactions between replication fork components and cohesin facilitate successful establishment of sister chromatid cohesion during DNA replication.

Comprehensive analysis of DNA replication timing across 184 cell lines suggests a role for MCM10 in replication timing regulation

Cellular proliferation depends on the accurate and timely replication of the genome. Several genetic diseases are caused by mutations in key DNA replication genes; however, it remains unclear whether these genes influence the normal program of DNA replication timing. Similarly, the factors that regulate DNA replication dynamics are poorly understood. To systematically identify trans-acting modulators of replication timing, we profiled replication in 184 cell lines from three cell types, encompassing 60 different gene knockouts or genetic diseases. Through a rigorous approach that considers the background variability of replication timing, we concluded that most samples displayed normal replication timing. However, mutations in two genes showed consistently abnormal replication timing. The first gene was RIF1, a known modulator of replication timing. The second was MCM10, a highly conserved member of the pre-replication complex. Cells from a single patient carrying MCM10 mutations demonstrated replication timing variability comprising 46% of the genome and at different locations than RIF1 knockouts. Replication timing alterations in the mutated MCM10 cells were predominantly comprised of replication delays and initiation site gains and losses. Taken together, this study demonstrates the remarkable robustness of the human replication timing program and reveals MCM10 as a novel candidate modulator of DNA replication timing.

S Phase Duration Is Determined by Local Rate and Global Organization of Replication

The duration of the cell cycle has been extensively studied and a wide degree of variability exists between cells, tissues and organisms. However, the duration of S phase has often been neglected, due to the false assumption that S phase duration is relatively constant. In this paper, we describe the methodologies to measure S phase duration, summarize the existing knowledge about its variability and discuss the key factors that control it. The local rate of replication (LRR), which is a combination of fork rate (FR) and inter-origin distance (IOD), has a limited influence on S phase duration, partially due to the compensation between FR and IOD. On the other hand, the organization of the replication program, specifically the amount of replication domains that fire simultaneously and the degree of overlap between the firing of distinct replication timing domains, is the main determinant of S phase duration. We use these principles to explain the variation in S phase length in different tissues and conditions.

Stimulation of adaptive gene amplification by origin firing under replication fork constraint

Adaptive mutations can cause drug resistance in cancers and pathogens, and increase the tolerance of agricultural pests and diseases to chemical treatment. When and how adaptive mutations form is often hard to discern, but we have shown that adaptive copy number amplification of the copper resistance gene CUP1 occurs in response to environmental copper due to CUP1 transcriptional activation. Here we dissect the mechanism by which CUP1 transcription in budding yeast stimulates copy number variation (CNV). We show that transcriptionally stimulated CNV requires TREX-2 and Mediator, such that cells lacking TREX-2 or Mediator respond normally to copper but cannot acquire increased resistance. Mediator and TREX-2 can cause replication stress by tethering transcribed loci to nuclear pores, a process known as gene gating, and transcription at the CUP1 locus causes a TREX-2-dependent accumulation of replication forks indicative of replication fork stalling. TREX-2-dependent CUP1 gene amplification occurs by a Rad52 and Rad51-mediated homologous recombination mechanism that is enhanced by histone H3K56 acetylation and repressed by Pol32 and Pif1. CUP1 amplification is also critically dependent on late-firing replication origins present in the CUP1 repeats, and mutations that remove or inactivate these origins strongly suppress the acquisition of copper resistance. We propose that replicative stress imposed by nuclear pore association causes replication bubbles from these origins to collapse soon after activation, leaving a tract of H3K56-acetylated chromatin that promotes secondary recombination events during elongation after replication fork re-start events. The capacity for inefficient replication origins to promote copy number variation renders certain genomic regions more fragile than others, and therefore more likely to undergo adaptive evolution through de novo gene amplification.

Replication Stress, Genomic Instability, and Replication Timing: A Complex Relationship

The replication-timing program constitutes a key element of the organization and coordination of numerous nuclear processes in eukaryotes. This program is established at a crucial moment in the cell cycle and occurs simultaneously with the organization of the genome, thus indicating the vital significance of this process. With recent technological achievements of high-throughput approaches, a very strong link has been confirmed between replication timing, transcriptional activity, the epigenetic and mutational landscape, and the 3D organization of the genome. There is also a clear relationship between replication stress, replication timing, and genomic instability, but the extent to which they are mutually linked to each other is unclear. Recent evidence has shown that replication timing is affected in cancer cells, although the cause and consequence of this effect remain unknown. However, in-depth studies remain to be performed to characterize the molecular mechanisms of replication-timing regulation and clearly identify different cis- and trans-acting factors. The results of these studies will potentially facilitate the discovery of new therapeutic pathways, particularly for personalized medicine, or new biomarkers. This review focuses on the complex relationship between replication timing, replication stress, and genomic instability.

Histones on fire: the effect of Dun1 and Mrc1 on origin firing and replication of hyper-acetylated genomes

As cells replicate their DNA, there is a need to synthesize new histones with which to wrap it. Newly synthesized H3 histones that are incorporated into the assembling chromatin behind the replication fork are acetylated at lysine 56. The acetylation is removed by two deacetylases, Hst3 and Hst4. This process is tightly regulated and any perturbation leads to genomic instability and replicative stress. We recently showed that Dun1, a kinase implicated mainly in the regulation of dNTPs, is vital in cells with hyper-acetylation, to counteract Rad53's inhibition on late-firing origins of replication. Our work showed that ∆hst3 ∆hst4 cells depend on late origin firing for survival, and are unable to prevent Rad53's inhibition when Dun1 is inactive. Thus, our work describes a role for Dun1 that is independent on its known function as a regulator of dNTP levels. Here we show that Mrc1 (Claspin in mammals), a protein that moves with the replicating fork and participates in both replication and checkpoint functions, plays also an essential role in the absence of H3K56Ac deacetylation. The sum of the results shown here and in our recent publication suggests that dormant origins are also utilized in these cells, making Mrc1, which regulates firing from these origins, also essential when histone H3 is hyper-acetylated. Thus, cells suffering from hyper-acetylation of H3K56 experience replication stress caused by a combination of prone-to-collapse forks and limited replication tracts. This combination makes both Dun1 and Mrc1, each acting on different targets, essential for viability.

Thermus thermophilus Argonaute Functions in the Completion of DNA Replication

In many eukaryotes, Argonaute proteins, guided by short RNA sequences, defend cells against transposons and viruses. In the eubacterium Thermus thermophilus, the DNA-guided Argonaute TtAgo defends against transformation by DNA plasmids. Here, we report that TtAgo also participates in DNA replication. In vivo, TtAgo binds 15- to 18-nt DNA guides derived from the chromosomal region where replication terminates and associates with proteins known to act in DNA replication. When gyrase, the sole T. thermophilus type II topoisomerase, is inhibited, TtAgo allows the bacterium to finish replicating its circular genome. In contrast, loss of gyrase and TtAgo activity slows growth and produces long sausage-like filaments in which the individual bacteria are linked by DNA. Finally, wild-type T. thermophilus outcompetes an otherwise isogenic strain lacking TtAgo. We propose that the primary role of TtAgo is to help T. thermophilus disentangle the catenated circular chromosomes generated by DNA replication.

Effects of cell cycle variability on lineage and population measurements of messenger RNA abundance

Many models of gene expression do not explicitly incorporate a cell cycle description. Here, we derive a theory describing how messenger RNA (mRNA) fluctuations for constitutive and bursty gene expression are influenced by stochasticity in the duration of the cell cycle and the timing of DNA replication. Analytical expressions for the moments show that omitting cell cycle duration introduces an error in the predicted mean number of mRNAs that is a monotonically decreasing function of η, which is proportional to the ratio of the mean cell cycle duration and the mRNA lifetime. By contrast, the error in the variance of the mRNA distribution is highest for intermediate values of η consistent with genome-wide measurements in many organisms. Using eukaryotic cell data, we estimate the errors in the mean and variance to be at most 3% and 25%, respectively. Furthermore, we derive an accurate negative binomial mixture approximation to the mRNA distribution. This indicates that stochasticity in the cell cycle can introduce fluctuations in mRNA numbers that are similar to the effect of bursty transcription. Finally, we show that for real experimental data, disregarding cell cycle stochasticity can introduce errors in the inference of transcription rates larger than 10%.

Genomic methods for measuring DNA replication dynamics

Genomic DNA replicates according to a defined temporal program in which early-replicating loci are associated with open chromatin, higher gene density, and increased gene expression levels, while late-replicating loci tend to be heterochromatic and show higher rates of genomic instability. The ability to measure DNA replication dynamics at genome scale has proven crucial for understanding the mechanisms and cellular consequences of DNA replication timing. Several methods, such as quantification of nucleotide analog incorporation and DNA copy number analyses, can accurately reconstruct the genomic replication timing profiles of various species and cell types. More recent developments have expanded the DNA replication genomic toolkit to assays that directly measure the activity of replication origins, while single-cell replication timing assays are beginning to reveal a new level of replication timing regulation. The combination of these methods, applied on a genomic scale and in multiple biological systems, promises to resolve many open questions and lead to a holistic understanding of how eukaryotic cells replicate their genomes accurately and efficiently.

DNA copy-number measurement of genome replication dynamics by high-throughput sequencing: the sort-seq, sync-seq and MFA-seq family

Genome replication follows a defined temporal programme that can change during cellular differentiation and disease onset. DNA replication results in an increase in DNA copy number that can be measured by high-throughput sequencing. Here we present a protocol to determine genome replication dynamics using DNA copy-number measurements. Cell populations can be obtained in three variants of the method. First, sort-seq reveals the average replication dynamics across S phase in an unperturbed cell population; FACS is used to isolate replicating and non-replicating subpopulations from asynchronous cells. Second, sync-seq measures absolute replication time at specific points during S phase using a synchronized cell population. Third, marker frequency analysis can be used to reveal the average replication dynamics using copy-number analysis in any proliferating asynchronous cell culture. These approaches have been used to reveal genome replication dynamics in prokaryotes, archaea and a wide range of eukaryotes, including yeasts and mammalian cells. We have found this approach straightforward to apply to other organisms and highlight example studies from across the three domains of life. Here we present a Saccharomyces cerevisiae version of the protocol that can be performed in 7-10 d. It requires basic molecular and cellular biology skills, as well as a basic understanding of Unix and R.

The evolutionary plasticity of chromosome metabolism allows adaptation to constitutive DNA replication stress

Many biological features are conserved and thus considered to be resistant to evolutionary change. While rapid genetic adaptation following the removal of conserved genes has been observed, we often lack a mechanistic understanding of how adaptation happens. We used the budding yeast, Saccharomyces cerevisiae, to investigate the evolutionary plasticity of chromosome metabolism, a network of evolutionary conserved modules. We experimentally evolved cells constitutively experiencing DNA replication stress caused by the absence of Ctf4, a protein that coordinates the enzymatic activities at replication forks. Parallel populations adapted to replication stress, over 1000 generations, by acquiring multiple, concerted mutations. These mutations altered conserved features of two chromosome metabolism modules, DNA replication and sister chromatid cohesion, and inactivated a third, the DNA damage checkpoint. The selected mutations define a functionally reproducible evolutionary trajectory. We suggest that the evolutionary plasticity of chromosome metabolism has implications for genome evolution in natural populations and cancer.

Molecular signatures of aneuploidy-driven adaptive evolution

Alteration of normal ploidy (aneuploidy) can have a number of opposing effects, such as unbalancing protein abundances and inhibiting cell growth but also accelerating genetic diversification and rapid adaptation. The interplay of these detrimental and beneficial effects remains puzzling. Here, to understand how cells develop tolerance to aneuploidy, we subject disomic (i.e. with an extra chromosome copy) strains of yeast to long-term experimental evolution under strong selection, by forcing disomy maintenance and daily population dilution. We characterize mutations, karyotype alterations and gene expression changes, and dissect the associated molecular strategies. Cells with different extra chromosomes accumulated mutations at distinct rates and displayed diverse adaptive events. They tended to evolve towards normal ploidy through chromosomal DNA loss and gene expression changes. We identify genes with recurrent mutations and altered expression in multiple lines, revealing a variant that improves growth under genotoxic stresses. These findings support rapid evolvability of disomic strains that can be used to characterize fitness effects of mutations under different stress conditions.

Aneuploidy (abnormal chromosome number) can enable rapid adaptation to stress conditions, but it also entails fitness costs from gene imbalance. Here, the authors experimentally evolve yeast while forcing maintenance of aneuploidy to identify the mechanisms that promote tolerance of aneuploidy.

Multiple inputs ensure yeast cell size homeostasis during cell cycle progression

Coordination of cell growth with division is essential for proper cell function. In budding yeast, although some molecular mechanisms responsible for cell size control during G1 have been elucidated, the mechanism by which cell size homeostasis is established remains to be discovered. Here, we developed a new technique based on quantification of histone levels to monitor cell cycle progression in individual cells with unprecedented accuracy. Our analysis establishes the existence of a mechanism controlling bud size in G2/M that prevents premature onset of anaphase, and controls the overall size variability. While most G1 mutants do not display impaired size homeostasis, mutants in which cyclin B-Cdk regulation is altered display large size variability. Our study thus demonstrates that size homeostasis is not controlled by a G1-specific mechanism alone but is likely to be an emergent property resulting from the integration of several mechanisms that coordinate cell and bud growth with division.

G2-seq: A High Throughput Sequencing-based Technique for Identifying Late Replicating Regions of the Genome

Numerous techniques have been developed to follow the progress of DNA replication through the S phase of the cell cycle. Most of these techniques have been directed toward elucidation of the location and timing of initiation of genome duplication rather than its completion. However, it is critical that we understand regions of the genome that are last to complete replication, because these regions suffer elevated levels of chromosomal breakage and mutation, and they have been associated with both disease and aging. Here we describe how we have extended a technique that has been used to monitor replication initiation to instead identify those regions of the genome last to complete replication. This approach is based on a combination of flow cytometry and high throughput sequencing. Although this report focuses on the application of this technique to yeast, the approach can be used with any cells that can be sorted in a flow cytometer according to DNA content.

Cell-to-cell variability and robustness in S-phase duration from genome replication kinetics

Genome replication, a key process for a cell, relies on stochastic initiation by replication origins, causing a variability of replication timing from cell to cell. While stochastic models of eukaryotic replication are widely available, the link between the key parameters and overall replication timing has not been addressed systematically. We use a combined analytical and computational approach to calculate how positions and strength of many origins lead to a given cell-to-cell variability of total duration of the replication of a large region, a chromosome or the entire genome. Specifically, the total replication timing can be framed as an extreme-value problem, since it is due to the last region that replicates in each cell. Our calculations identify two regimes based on the spread between characteristic completion times of all inter-origin regions of a genome. For widely different completion times, timing is set by the single specific region that is typically the last to replicate in all cells. Conversely, when the completion time of all regions are comparable, an extreme-value estimate shows that the cell-to-cell variability of genome replication timing has universal properties. Comparison with available data shows that the replication program of three yeast species falls in this extreme-value regime.

Perturbations in the Replication Program Contribute to Genomic Instability in Cancer

Cancer and genomic instability are highly impacted by the deoxyribonucleic acid (DNA) replication program. Inaccuracies in DNA replication lead to the increased acquisition of mutations and structural variations. These inaccuracies mainly stem from loss of DNA fidelity due to replication stress or due to aberrations in the temporal organization of the replication process. Here we review the mechanisms and impact of these major sources of error to the replication program.

A role for Rtt109 in buffering gene-dosage imbalance during DNA replication

Chromatin can function as an integrator of DNA-related processes, allowing communication, for example, between DNA replication and gene transcription. Such communication is needed to overcome the gene-dosage imbalance introduced during DNA replication, when certain genes are replicated prior to others. Increased transcription of early replicating genes could alter regulatory balances. This does not occur, suggesting a mechanism that suppresses expression from newly replicated DNA. Critical to this buffering is Rtt109, which acetylates the internal K56 residue of newly synthesized histone H3 prior to incorporation onto DNA. H3K56ac distinguishes replicated from non-replicated DNA, communicating this information to the transcription machinery to ensure expression homeostasis during S phase.

Alterations in cellular metabolism triggered by URA7 or GLN3 inactivation cause imbalanced dNTP pools and increased mutagenesis

Eukaryotic DNA replication fidelity relies on the concerted action of DNA polymerase nucleotide selectivity, proofreading activity, and DNA mismatch repair (MMR). Nucleotide selectivity and proofreading are affected by the balance and concentration of deoxyribonucleotide (dNTP) pools, which are strictly regulated by ribonucleotide reductase (RNR). Mutations preventing DNA polymerase proofreading activity or MMR function cause mutator phenotypes and consequently increased cancer susceptibility. To identify genes not previously linked to high-fidelity DNA replication, we conducted a genome-wide screen in Saccharomyces cerevisiae using DNA polymerase active-site mutants as a "sensitized mutator background." Among the genes identified in our screen, three metabolism-related genes (GLN3, URA7, and SHM2) have not been previously associated to the suppression of mutations. Loss of either the transcription factor Gln3 or inactivation of the CTP synthetase Ura7 both resulted in the activation of the DNA damage response and imbalanced dNTP pools. Importantly, these dNTP imbalances are strongly mutagenic in genetic backgrounds where DNA polymerase function or MMR activity is partially compromised. Previous reports have shown that dNTP pool imbalances can be caused by mutations altering the allosteric regulation of enzymes involved in dNTP biosynthesis (e.g., RNR or dCMP deaminase). Here, we provide evidence that mutations affecting genes involved in RNR substrate production can cause dNTP imbalances, which cannot be compensated by RNR or other enzymatic activities. Moreover, Gln3 inactivation links nutrient deprivation to increased mutagenesis. Our results suggest that similar genetic interactions could drive mutator phenotypes in cancer cells.

Checkpoint-Independent Regulation of Origin Firing by Mrc1 through Interaction with Hsk1 Kinase

Mrc1 is a conserved checkpoint mediator protein that transduces the replication stress signal to the downstream effector kinase. The loss of mrc1 checkpoint activity results in the aberrant activation of late/dormant origins in the presence of hydroxyurea. Mrc1 was also suggested to regulate orders of early origin firing in a checkpoint-independent manner, but its mechanism was unknown. Here we identify HBS (Hsk1 bypass segment) on Mrc1. An ΔHBS mutant does not activate late/dormant origin firing in the presence of hydroxyurea but causes the precocious and enhanced activation of weak early-firing origins during normal S-phase progression and bypasses the requirement for Hsk1 for growth. This may be caused by the disruption of intramolecular binding between HBS and NTHBS (N-terminal target of HBS). Hsk1 binds to Mrc1 through HBS and phosphorylates a segment adjacent to NTHBS, disrupting the intramolecular interaction. We propose that Mrc1 exerts a "brake" on initiation (through intramolecular interactions) and that this brake can be released (upon the loss of intramolecular interactions) by either the Hsk1-mediated phosphorylation of Mrc1 or the deletion of HBS (or a phosphomimic mutation of putative Hsk1 target serine/threonine), which can bypass the function of Hsk1 for growth. The brake mechanism may explain the checkpoint-independent regulation of early origin firing in fission yeast.

Coupling phenotypic persistence to DNA damage increases genetic diversity in severe stress

Mutation rate balances the need to protect genome integrity with the advantage of evolutionary innovations. Microorganisms increase their mutation rate when stressed, perhaps addressing the growing need for evolutionary innovation. Such a strategy, however, is only beneficial under moderate stresses that allow cells to divide and realize their mutagenic potential. In contrast, severe stresses rapidly kill the majority of the population with the exception of a small minority of cells that are in a phenotypically distinct state termed persistence. Although persisters were discovered many decades ago, the stochastic event triggering persistence is poorly understood. We report that spontaneous DNA damage triggers persistence in Saccharomyces cerevisiae by activating the general stress response, providing protection against a range of harsh stress and drug environments. We further show that the persister subpopulation carries an increased load of genetic variants in the form of insertions, deletions or large structural variations, which are unrelated to their stress survival. This coupling of DNA damage to phenotypic persistence may increase genetic diversity specifically in severe stress conditions, where diversity is beneficial but the ability to generate de novo mutations is limited.

SIR2 suppresses replication gaps and genome instability by balancing replication between repetitive and unique sequences

Replication gaps that persist into mitosis likely represent important threats to genome stability, but experimental identification of these gaps has proved challenging. We have developed a technique that allows us to explore the dynamics by which genome replication is completed before mitosis. Using this approach, we demonstrate that excessive allocation of replication resources to origins within repetitive regions, induced by SIR2 deletion, leads to persistent replication gaps and genome instability. Conversely, the weakening of replication origins in repetitive regions suppresses these gaps. Given known age- and cancer-associated changes in chromatin accessibility at repetitive sequences, we suggest that replication gaps resulting from misallocation of replication resources underlie age- and disease-associated genome instability.

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.

The Genomic Landscape of Position Effects on Protein Expression Level and Noise in Yeast

Position effect, the influence of the chromosomal location of a gene on its activity, is a fundamental property of the genome. By placing a GFP gene cassette at 482 different locations across all chromosomes in budding yeast, we quantified the position effects on protein expression level and noise at the genomic scale. The position effects are significant, altering the mean protein expression level by up to 15 times and expression noise by up to 20 times. DNA replication timing, 3D chromosomal conformation, and several histone modifications are major covariates of position effects. Essential genes are enriched in genomic regions with inherently low expression noise, supporting the hypothesis that chromosomal clustering of essential genes results from selection against their expressional stochasticity. Position effects exhibit significant interactions with promoters. Together, our results suggest that position effects have shaped the evolution of chromosome organization and should inform future genome engineering efforts.

Detecting ATM-dependent chromatin modification in DNA damage response

Loss of function or mutation of the ataxia-telangiectasia mutated gene product (ATM) results in inherited genetic disorders characterized by neurodegeneration, immunodeficiency, and cancer. Ataxia-telangiectasia mutated (ATM) gene product belongs to the PI3K-like protein kinase (PIKKs) family and is functionally implicated in mitogenic signal transduction, chromosome condensation, meiotic recombination, cell-cycle control, and telomere maintenance. The ATM protein kinase is primarily activated in response to DNA double strand breaks (DSBs), the most deleterious form of DNA damage produced by ionizing radiation (IR) or radiomimetic drugs. It is detected at DNA damage sites, where ATM autophosphorylation causes dissociation of the inactive homodimeric form to the activated monomeric form. Interestingly, heat shock can activate ATM independent of the presence of DNA strand breaks. ATM is an integral part of the sensory machinery that detects DSBs during meiosis, mitosis, or DNA breaks mediated by free radicals. These DNA lesions can trigger higher order chromatin reorganization fuelled by posttranslational modifications of histones and histone binding proteins. Our group, and others, have shown that ATM activation is tightly regulated by chromatin modifications. This review summarizes the multiple approaches used to discern the role of ATM and other associated proteins in chromatin modification in response to DNA damage.

Nascent RNA folding mitigates transcription-associated mutagenesis

Transcription is mutagenic, in part because the R-loop formed by the binding of the nascent RNA with its DNA template exposes the nontemplate DNA strand to mutagens and primes unscheduled error-prone DNA synthesis. We hypothesize that strong folding of nascent RNA weakens R-loops and hence decreases mutagenesis. By a yeast forward mutation assay, we show that strengthening RNA folding and reducing R-loop formation by synonymous changes in a reporter gene can lower mutation rate by >80%. This effect is diminished after the overexpression of the gene encoding RNase H1 that degrades the RNA in a DNA–RNA hybrid, indicating that the effect is R-loop-dependent. Analysis of genomic data of yeast mutation accumulation lines and human neutral polymorphisms confirms the generality of these findings. This mechanism for local protection of genome integrity is of special importance to highly expressed genes because of their frequent transcription and strong RNA folding, the latter also improves translational fidelity. As a result, strengthening RNA folding simultaneously curtails genotypic and phenotypic mutations.

Replication timing is regulated by the number of MCMs loaded at origins

Replication timing is a crucial aspect of genome regulation that is strongly correlated with chromatin structure, gene expression, DNA repair, and genome evolution. Replication timing is determined by the timing of replication origin firing, which involves activation of MCM helicase complexes loaded at replication origins. Nonetheless, how the timing of such origin firing is regulated remains mysterious. Here, we show that the number of MCMs loaded at origins regulates replication timing. We show for the first time in vivo that multiple MCMs are loaded at origins. Because early origins have more MCMs loaded, they are, on average, more likely to fire early in S phase. Our results provide a mechanistic explanation for the observed heterogeneity in origin firing and help to explain how defined replication timing profiles emerge from stochastic origin firing. These results establish a framework in which further mechanistic studies on replication timing, such as the strong effect of heterochromatin, can be pursued.

Systematic identification of cell size regulators in budding yeast

Cell size is determined by a complex interplay between growth and division, involving multiple cellular pathways. To identify systematically processes affecting size control in G1 in budding yeast, we imaged and analyzed the cell cycle of millions of individual cells representing 591 mutants implicated in size control. Quantitative metric distinguished mutants affecting the mechanism of size control from the majority of mutants that have a perturbed size due to indirect effects modulating cell growth. Overall, we identified 17 negative and dozens positive size control regulators, with the negative regulators forming a small network centered on elements of mitotic exit network. Some elements of the translation machinery affected size control with a notable distinction between the deletions of parts of small and large ribosomal subunit: parts of small ribosomal subunit tended to regulate size control, while parts of the large subunit affected cell growth. Analysis of small cells revealed additional size control mechanism that functions in G2/M, complementing the primary size control in G1. Our study provides new insights about size control mechanisms in budding yeast.

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.

DNA replication timing: Coordinating genome stability with genome regulation on the X chromosome and beyond

Recent studies based on next-generation DNA sequencing have revealed that the female inactive X chromosome is replicated in a rapid, unorganized manner, and undergoes increased rates of mutation. These observations link the organization of DNA replication timing to gene regulation on one hand, and to the generation of mutations on the other hand. More generally, the exceptional biology of the inactive X chromosome highlights general principles of genome replication. Cells may control replication timing by a combination of intrinsic replication origin properties, local chromatin states and global levels of replication factors, leading to a functional separation between the activity of genes and their mutation.

Telomere length kinetics assay (TELKA) sorts the telomere length maintenance (tlm) mutants into functional groups

Genome-wide systematic screens in yeast have uncovered a large gene network (the telomere length maintenance network or TLM), encompassing more than 400 genes, which acts coordinatively to maintain telomere length. Identifying the genes was an important first stage; the next challenge is to decipher their mechanism of action and to organize then into functional groups or pathways. Here we present a new telomere-length measuring program, TelQuant, and a novel assay, telomere length kinetics assay, and use them to organize tlm mutants into functional classes. Our results show that a mutant defective for the relatively unknown MET7 gene has the same telomeric kinetics as mutants defective for the ribonucleotide reductase subunit Rnr1, in charge of the limiting step in dNTP synthesis, or for the Ku heterodimer, a well-established telomere complex. We confirm the epistatic relationship between the mutants and show that physical interactions exist between Rnr1 and Met7. We also show that Met7 and the Ku heterodimer affect dNTP formation, and play a role in non-homologous end joining. Thus, our telomere kinetics assay uncovers new functional groups, as well as complex genetic interactions between tlm mutants.

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.

Time to be versatile: regulation of the replication timing program in budding yeast

Eukaryotic replication origins are activated at different times during the S phase of the cell cycle, following a temporal program that is stably transmitted to daughter cells. Although the mechanisms that control initiation at the level of individual origins are now well understood, much less is known on how cells coordinate replication at hundreds of origins distributed on the chromosomes. In this review, we discuss recent advances shedding new light on how this complex process is regulated in the budding yeast Saccharomyces cerevisiae. The picture that emerges from these studies is that replication timing is regulated in cis by mechanisms modulating the chromatin structure and the subnuclear organization of origins. These mechanisms do not affect the licensing of replication origins but determine their ability to compete for limiting initiation factors, which are recycled from early to late origins throughout the length of the S phase.

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.

Distinct SUMO ligases cooperate with Esc2 and Slx5 to suppress duplication-mediated genome rearrangements

Suppression of duplication-mediated gross chromosomal rearrangements (GCRs) is essential to maintain genome integrity in eukaryotes. Here we report that SUMO ligase Mms21 has a strong role in suppressing GCRs in Saccharomyces cerevisiae, while Siz1 and Siz2 have weaker and partially redundant roles. Understanding the functions of these enzymes has been hampered by a paucity of knowledge of their substrate specificity in vivo. Using a new quantitative SUMO-proteomics technology, we found that Siz1 and Siz2 redundantly control the abundances of most sumoylated substrates, while Mms21 more specifically regulates sumoylation of RNA polymerase-I and the SMC-family proteins. Interestingly, Esc2, a SUMO-like domain-containing protein, specifically promotes the accumulation of sumoylated Mms21-specific substrates and functions with Mms21 to suppress GCRs. On the other hand, the Slx5-Slx8 complex, a SUMO-targeted ubiquitin ligase, suppresses the accumulation of sumoylated Mms21-specific substrates. Thus, distinct SUMO ligases work in concert with Esc2 and Slx5-Slx8 to control substrate specificity and sumoylation homeostasis to prevent GCRs.

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.

Saccharomyces cerevisiae deletion strains with complex DNA content profiles

To identify Saccharomyces cerevisiae genes required for the proper timing of cell cycle transitions, we previously reported a systematic examination of the DNA content of homozygous diploid deletion strains. However, deletion strains with complex DNA content profiles were not examined in that study. Here, we report S. cerevisiae genes that when deleted give rise to DNA content profiles consistent with roles of the corresponding gene products during DNA replication. We also identified a set of genes whose deletion leads to increased DNA content, consistent with defects in mitosis, cytokinesis, or cell separation. Finally, we examined known interactions between the gene products of each group, placing these gene products in functional networks. Taken together, the data we present further validate the roles of the corresponding gene products in these processes, facilitating efforts to delineate gene function critical for genome replication, maintenance, and segregation.

DNA replication stress-induced loss of reproductive capacity in S. cerevisiae and its inhibition by caloric restriction

In many organisms, attenuation of growth signaling by caloric restriction or mutational inactivation of growth signaling pathways extends lifespan and protects against cancer and other age-related diseases. The focus of many efforts to understand these effects has been on the induction of oxidative stress defenses that inhibit cellular senescence and cell death. Here we show that in the model organism S. cerevisiae, growth signaling induces entry of cells in stationary phase into S phase in parallel with loss of reproductive capacity, which is enhanced by elevated concentrations of glucose. Overexpression of RNR1 encoding a ribonucleotide reductase subunit required for the synthesis of deoxynucleotide triphosphates and DNA replication suppresses the accelerated loss of reproductive capacity of cells cultured in high glucose. The reduced reproductive capacity of these cells is also suppressed by excess threonine, which buffers dNTP pools when ribonucleotide reductase activity is limiting. Caloric restriction or inactivation of the AKT homolog Sch9p inhibits senescence and death in stationary phase cells caused by the DNA replication inhibitor hydroxyurea or by inactivation of the DNA replication and repair proteins Sgs1p or Rad27p. Inhibition of DNA replication stress represents a novel mechanism by which caloric restriction promotes longevity in S. cerevisiae. A similar mechanism may promote longevity and inhibit cancer and other age-related diseases in humans.

Differential relationship of DNA replication timing to different forms of human mutation and variation

Human genetic variation is distributed nonrandomly across the genome, though the principles governing its distribution are only partially known. DNA replication creates opportunities for mutation, and the timing of DNA replication correlates with the density of SNPs across the human genome. To enable deeper investigation of how DNA replication timing relates to human mutation and variation, we generated a high-resolution map of the human genome's replication timing program and analyzed its relationship to point mutations, copy number variations, and the meiotic recombination hotspots utilized by males and females. DNA replication timing associated with point mutations far more strongly than predicted from earlier analyses and showed a stronger relationship to transversion than transition mutations. Structural mutations arising from recombination-based mechanisms and recombination hotspots used more extensively by females were enriched in early-replicating parts of the genome, though these relationships appeared to relate more strongly to the genomic distribution of causative sequence features. These results indicate differential and sex-specific relationship of DNA replication timing to different forms of mutation and recombination.

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.

Conservation of replication timing reveals global and local regulation of replication origin activity

DNA replication initiates from defined locations called replication origins; some origins are highly active, whereas others are dormant and rarely used. Origins also differ in their activation time, resulting in particular genomic regions replicating at characteristic times and in a defined temporal order. Here we report the comparison of genome replication in four budding yeast species: Saccharomyces cerevisiae, S. paradoxus, S. arboricolus, and S. bayanus. First, we find that the locations of active origins are predominantly conserved between species, whereas dormant origins are poorly conserved. Second, we generated genome-wide replication profiles for each of these species and discovered that the temporal order of genome replication is highly conserved. Therefore, active origins are not only conserved in location, but also in activation time. Only a minority of these conserved origins show differences in activation time between these species. To gain insight as to the mechanisms by which origin activation time is regulated we generated replication profiles for a S. cerevisiae/S. bayanus hybrid strain and find that there are both local and global regulators of origin function.

Separation of DNA replication from the assembly of break-competent meiotic chromosomes

The meiotic cell division reduces the chromosome number from diploid to haploid to form gametes for sexual reproduction. Although much progress has been made in understanding meiotic recombination and the two meiotic divisions, the processes leading up to recombination, including the prolonged pre-meiotic S phase (meiS) and the assembly of meiotic chromosome axes, remain poorly defined. We have used genome-wide approaches in Saccharomyces cerevisiae to measure the kinetics of pre-meiotic DNA replication and to investigate the interdependencies between replication and axis formation. We found that replication initiation was delayed for a large number of origins in meiS compared to mitosis and that meiotic cells were far more sensitive to replication inhibition, most likely due to the starvation conditions required for meiotic induction. Moreover, replication initiation was delayed even in the absence of chromosome axes, indicating replication timing is independent of the process of axis assembly. Finally, we found that cells were able to install axis components and initiate recombination on unreplicated DNA. Thus, although pre-meiotic DNA replication and meiotic chromosome axis formation occur concurrently, they are not strictly coupled. The functional separation of these processes reveals a modular method of building meiotic chromosomes and predicts that any crosstalk between these modules must occur through superimposed regulatory mechanisms.

Replication timing and its emergence from stochastic processes

The temporal organization of DNA replication has puzzled cell biologists since before the mechanism of replication was understood. The realization that replication timing correlates with important features, such as transcription, chromatin structure and genome evolution, and is misregulated in cancer and aging has only deepened the fascination. Many ideas about replication timing have been proposed, but most have been short on mechanistic detail. However, recent work has begun to elucidate basic principles of replication timing. In particular, mathematical modeling of replication kinetics in several systems has shown that the reproducible replication timing patterns seen in population studies can be explained by stochastic origin firing at the single-cell level. This work suggests that replication timing need not be controlled by a hierarchical mechanism that imposes replication timing from a central regulator, but instead results from simple rules that affect individual origins.

Multiple pathways can bypass the essential role of fission yeast Hsk1 kinase in DNA replication initiation

A number of different genetic backgrounds and growth conditions bypass DNA replication defects caused by the absence of yeast Hsk1 kinase, demonstrating the plasticity of the eukaryotic DNA replication program.

Cdc7/Hsk1 is a conserved kinase required for initiation of DNA replication that potentially regulates timing and locations of replication origin firing. Here, we show that viability of fission yeast hsk1Δ cells can be restored by loss of mrc1, which is required for maintenance of replication fork integrity, by cds1Δ, or by a checkpoint-deficient mutant of mrc1. In these mutants, normally inactive origins are activated in the presence of hydroxyurea and binding of Cdc45 to MCM is stimulated. mrc1Δ bypasses hsk1Δ more efficiently because of its checkpoint-independent inhibitory functions. Unexpectedly, hsk1Δ is viable at 37°C. More DNA is synthesized, and some dormant origins fire in the presence of hydroxyurea at 37°C. Furthermore, hsk1Δ bypass strains grow poorly at 25°C compared with higher temperatures. Our results show that Hsk1 functions for DNA replication can be bypassed by different genetic backgrounds as well as under varied physiological conditions, providing additional evidence for plasticity of the replication program in eukaryotes.

Genetic variation and DNA replication timing, or why is there late replicating DNA?

Mutation rates vary significantly within the genome and across species. Recent studies revealed a long suspected replication-timing effect on mutation rate, but the mechanisms that regulate the increase in mutation rate as the genome is replicated remain unclear. Evidence is emerging, however, that DNA repair systems, in general, are less efficient in late replicating heterochromatic regions compared to early replicating euchromatic regions of the genome. At the same time, mutation rates in both vertebrates and invertebrates have been shown to vary with generation time (GT). GT is correlated with genome size, which suggests a possible nucleotypic effect on species-specific mutation rates. These and other observations all converge on a role for DNA replication checkpoints in modulating generation times and mutation rates during the DNA synthetic phase (S phase) of the cell cycle. The following will examine the potential role of the intra-S checkpoint in regulating cell cycle times (GT) and mutation rates in eukaryotes. This article was published online on August 5, 2011. An error was subsequently identified. This notice is included in the online and print versions to indicate that both have been corrected October 4, 2011.

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

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