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

Unique molecular mechanisms for maintenance and alteration of genetic information in the budding yeast Saccharomyces cerevisiae

The high-fidelity transmission of genetic information is crucial for the survival of organisms, the cells of which have the ability to protect DNA against endogenous and environmental agents, including reactive oxygen species (ROS), ionizing radiation, and various chemical compounds. The basis of protection mechanisms has been evolutionarily conserved from yeast to humans; however, each organism often has a specialized mode of regulation that uses different sets of machineries, particularly in lower eukaryotes. The divergence of molecular mechanisms among related organisms has provided insights into the evolution of cellular machineries to a higher architecture. Uncommon characteristics of machineries may also contribute to the development of new applications such as drugs with novel mechanisms of action. In contrast to the cellular properties for maintaining genetic information, living organisms, particularly microbes, inevitably undergo genetic alterations in order to adapt to environmental conditions. The maintenance and alteration of genetic information may be inextricably linked to each other. In this review, we describe recent findings on the unconventional molecular mechanisms of DNA damage response and DNA double-strand break (DSB) repair in the budding yeast Saccharomyces cerevisiae. We also introduce our previous research on genetic and phenotypic instabilities observed in a clonal population of clinically-derived S. cerevisiae. The molecular mechanisms of this case were associated with mutations to generate tyrosine-inserting tRNA-Tyr ochre suppressors and the position effects of mutation frequencies among eight tRNA-Tyr loci dispersed in the genome. Phenotypic variations among different strain backgrounds have also been observed by another type of nonsense suppressor, the aberrant form of the translation termination factor. Nonsense suppressors are considered to be responsible for the genome-wide translational readthrough of termination codons, including natural nonsense codons. The nonsense suppressor-mediated acquisition of phenotypic variations may be advantageous for adaptation to environmental conditions and survival during evolution.

Genome Mutational and Transcriptional Hotspots Are Traps for Duplicated Genes and Sources of Adaptations

Gene duplication generates new genetic material, which has been shown to lead to major innovations in unicellular and multicellular organisms. A whole-genome duplication occurred in the ancestor of Saccharomyces yeast species but 92% of duplicates returned to single-copy genes shortly after duplication. The persisting duplicated genes in Saccharomyces led to the origin of major metabolic innovations, which have been the source of the unique biotechnological capabilities in the Baker’s yeast Saccharomyces cerevisiae. What factors have determined the fate of duplicated genes remains unknown. Here, we report the first demonstration that the local genome mutation and transcription rates determine the fate of duplicates. We show, for the first time, a preferential location of duplicated genes in the mutational and transcriptional hotspots of S. cerevisiae genome. The mechanism of duplication matters, with whole-genome duplicates exhibiting different preservation trends compared to small-scale duplicates. Genome mutational and transcriptional hotspots are rich in duplicates with large repetitive promoter elements. Saccharomyces cerevisiae shows more tolerance to deleterious mutations in duplicates with repetitive promoter elements, which in turn exhibit higher transcriptional plasticity against environmental perturbations. Our data demonstrate that the genome traps duplicates through the accelerated regulatory and functional divergence of their gene copies providing a source of novel adaptations in yeast.

Invisible cities: segregated domains in the yeast genome with distinct structural and functional attributes

Recent advances in our understanding of the three-dimensional organization of the eukaryotic nucleus have rendered the spatial distribution of genes increasingly relevant. In a recent work (Tsochatzidou et al., Nucleic Acids Res 45:5818-5828, 2017), we proposed the existence of a functional compartmentalization of the yeast genome according to which, genes occupying the chromosomal regions at the nuclear periphery have distinct structural, functional and evolutionary characteristics compared to their centromeric-proximal counterparts. Around the same time, it was also shown that the genome of Saccharomyces cerevisiae is organized in topologically associated domains (TADs), which are largely associated with the replication timing. In this work, we proceed to investigate whether such units of three-dimensional genomic organization can be linked to transcriptional activity as a driving force for the shaping of genomic architecture. Through the application of a simple boundary-calling criterion in genome-wide 3C data, we define ~100 TAD-like domains which can be clustered in six different classes with radically different nucleosomal organizations, significant variations in transcription factor binding and uneven chromosomal distribution. Approximately ~20% of the genome is found to be confined in regions with "closed" chromatin structure around gene promoters. Most interestingly, we find both "open" and "closed" regions to be segregated, in the sense that they tend to avoid inter-chromosomal interactions. Our data further enforce the notion of a marked compartmentalization of the yeast genome in isolated territories, with implications in its function and evolution.

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.

Genetic drift, selection and the evolution of the mutation rate

As one of the few cellular traits that can be quantified across the tree of life, DNA-replication fidelity provides an excellent platform for understanding fundamental evolutionary processes. Furthermore, because mutation is the ultimate source of all genetic variation, clarifying why mutation rates vary is crucial for understanding all areas of biology. A potentially revealing hypothesis for mutation-rate evolution is that natural selection primarily operates to improve replication fidelity, with the ultimate limits to what can be achieved set by the power of random genetic drift. This drift-barrier hypothesis is consistent with comparative measures of mutation rates, provides a simple explanation for the existence of error-prone polymerases and yields a formal counter-argument to the view that selection fine-tunes gene-specific mutation rates.

Building up and breaking down: mechanisms controlling recombination during replication

The complete and faithful duplication of the genome is an essential prerequisite for proliferating cells to maintain genome integrity. This objective is greatly challenged by DNA damage encountered during replication, which causes fork stalling and in certain cases, fork breakage. DNA damage tolerance (DDT) pathways mitigate the effects on fork stability induced by replication fork stalling by mediating damage-bypass and replication fork restart. These DDT mechanisms, largely relying on homologous recombination (HR) and specialized polymerases, can however contribute to genome rearrangements and mutagenesis. There is a profound connection between replication and recombination: recombination proteins protect replication forks from nuclease-mediated degradation of the nascent DNA strands and facilitate replication completion in cells challenged by DNA damage. Moreover, in case of fork collapse and formation of double strand breaks (DSBs), the recombination factors present or recruited to the fork facilitate HR-mediated DSB repair, which is primarily error-free. Disruption of HR is inexorably linked to genome instability, but the premature activation of HR during replication often leads to genome rearrangements. Faithful replication necessitates the downregulation of HR and disruption of active RAD51 filaments at replication forks, but upon persistent fork stalling, building up of HR is critical for the reorganization of the replication fork and for filling-in of the gaps associated with discontinuous replication induced by DNA lesions. Here we summarize and reflect on our understanding of the mechanisms that either suppress recombination or locally enhance it during replication, and the principles that underlie this regulation.

Historical Perspective of Eukaryotic DNA Replication

The replication of the genome of a eukaryotic cell is a complex process requiring the ordered assembly of multiprotein replisomes at many chromosomal sites. The process is strictly controlled during the cell cycle to ensure the complete and faithful transmission of genetic information to progeny cells. Our current understanding of the mechanisms of eukaryotic DNA replication has evolved over a period of more than 30 years through the efforts of many investigators. The aim of this perspective is to provide a brief history of the major advances during this period.

Coordination of DNA damage tolerance mechanisms with cell cycle progression in fission yeast

DNA damage tolerance (DDT) mechanisms allow cells to synthesize a new DNA strand when the template is damaged. Many mutations resulting from DNA damage in eukaryotes are generated during DDT when cells use the mutagenic translesion polymerases, Rev1 and Polζ, rather than mechanisms with higher fidelity. The coordination among DDT mechanisms is not well understood. We used live-cell imaging to study the function of DDT mechanisms throughout the cell cycle of the fission yeast Schizosaccharomyces pombe. We report that checkpoint-dependent mitotic delay provides a cellular mechanism to ensure the completion of high fidelity DDT, largely by homology-directed repair (HDR). DDT by mutagenic polymerases is suppressed during the checkpoint delay by a mechanism dependent on Rad51 recombinase. When cells pass the G2/M checkpoint and can no longer delay mitosis, they completely lose the capacity for HDR and simultaneously exhibit a requirement for Rev1 and Polζ. Thus, DDT is coordinated with the checkpoint response so that the activity of mutagenic polymerases is confined to a vulnerable period of the cell cycle when checkpoint delay and HDR are not possible.

Dosage Mutator Genes in Saccharomyces cerevisiae: A Novel Mutator Mode-of-Action of the Mph1 DNA Helicase

Mutations that cause genome instability are considered important predisposing events that contribute to initiation and progression of cancer. Genome instability arises either due to defects in genes that cause an increased mutation rate (mutator phenotype), or defects in genes that cause chromosome instability (CIN). To extend the catalog of genome instability genes, we systematically explored the effects of gene overexpression on mutation rate, using a forward-mutation screen in budding yeast. We screened ∼5100 plasmids, each overexpressing a unique single gene, and characterized the five strongest mutators, MPH1 (mutator phenotype 1), RRM3, UBP12, PIF1, and DNA2 We show that, for MPH1, the yeast homolog of Fanconi Anemia complementation group M (FANCM), the overexpression mutator phenotype is distinct from that of mph1Δ. Moreover, while four of our top hits encode DNA helicases, the overexpression of 48 other DNA helicases did not cause a mutator phenotype, suggesting this is not a general property of helicases. For Mph1 overexpression, helicase activity was not required for the mutator phenotype; in contrast Mph1 DEAH-box function was required for hypermutation. Mutagenesis by MPH1 overexpression was independent of translesion synthesis (TLS), but was suppressed by overexpression of RAD27, a conserved flap endonuclease. We propose that binding of DNA flap structures by excess Mph1 may block Rad27 action, creating a mutator phenotype that phenocopies rad27Δ. We believe this represents a novel mutator mode-of-action and opens up new prospects to understand how upregulation of DNA repair proteins may contribute to mutagenesis.

Preferential Protection of Genetic Fidelity within Open Chromatin by the Mismatch Repair Machinery

Epigenetic systems are well known for the roles they play in regulating the differential expression of the same genome in different cell types. However, epigenetic systems can also directly impact genomic integrity by protecting genetic sequences. Using an experimental evolutionary approach, we studied rates of mutation in the fission yeast Schizosaccharomyces pombe strains that lacked genes encoding several epigenetic regulators or mismatch repair components. We report that loss of a functional mismatch repair pathway in S. pombe resulted in the preferential enrichment of mutations in euchromatin, indicating that the mismatch repair machinery preferentially protected genetic fidelity in euchromatin. This preference is probably determined by differences in the accessibility of chromatin at distinct chromatin regions, which is supported by our observations that chromatin accessibility positively correlated with mutation rates in S. pombe or human cancer samples with deficiencies in mismatch repair. Importantly, such positive correlation was not observed in S. pombe strains or human cancer samples with functional mismatch repair machinery.

DNA damage tolerance by recombination: Molecular pathways and DNA structures

Replication perturbations activate DNA damage tolerance (DDT) pathways, which are crucial to promote replication completion and to prevent fork breakage, a leading cause of genome instability. One mode of DDT uses translesion synthesis polymerases, which however can also introduce mutations. The other DDT mode involves recombination-mediated mechanisms, which are generally accurate. DDT occurs prevalently postreplicatively, but in certain situations homologous recombination is needed to restart forks. Fork reversal can function to stabilize stalled forks, but may also promote error-prone outcome when used for fork restart. Recent years have witnessed important advances in our understanding of the mechanisms and DNA structures that mediate recombination-mediated damage-bypass and highlighted principles that regulate DDT pathway choice locally and temporally. In this review we summarize the current knowledge and paradoxes on recombination-mediated DDT pathways and their workings, discuss how the intermediate DNA structures may influence genome integrity, and outline key open questions for future research.

Classical Genetics with Saccharomyces cerevisiae

The budding yeast Saccharomyces cerevisiae is an outstanding experimental model organism that has been exploited since the early part of the twentieth century for studies in biochemistry and genetics. It has been the premiere experimental system for modern functional genomics and continues to make important contributions to many areas of biology. Here we discuss its many virtues as an organism for classical genetic research.

Mutation at a distance caused by homopolymeric guanine repeats in Saccharomyces cerevisiae

Mutation provides the raw material from which natural selection shapes adaptations. The rate at which new mutations arise is therefore a key factor that determines the tempo and mode of evolution. However, an accurate assessment of the mutation rate of a given organism is difficult because mutation rate varies on a fine scale within a genome. A central challenge of evolutionary genetics is to determine the underlying causes of this variation. In earlier work, we had shown that repeat sequences not only are prone to a high rate of expansion and contraction but also can cause an increase in mutation rate (on the order of kilobases) of the sequence surrounding the repeat. We perform experiments that show that simple guanine repeats 13 bp (base pairs) in length or longer (G 13+ ) increase the substitution rate 4- to 18-fold in the downstream DNA sequence, and this correlates with DNA replication timing (R = 0.89). We show that G 13+ mutagenicity results from the interplay of both error-prone translesion synthesis and homologous recombination repair pathways. The mutagenic repeats that we study have the potential to be exploited for the artificial elevation of mutation rate in systems biology and synthetic biology applications.

Rif1: A Conserved Regulator of DNA Replication and Repair Hijacked by Telomeres in Yeasts

Rap1-interacting factor 1 (Rif1) was originally identified in the budding yeast Saccharomyces cerevisiae as a telomere-binding protein that negatively regulates telomerase-mediated telomere elongation. Although this function is conserved in the distantly related fission yeast Schizosaccharomyces pombe, recent studies, both in yeasts and in metazoans, reveal that Rif1 also functions more globally, both in the temporal control of DNA replication and in DNA repair. Rif1 proteins are large and characterized by N-terminal HEAT repeats, predicted to form an elongated alpha-helical structure. In addition, all Rif1 homologs contain two short motifs, abbreviated RVxF/SILK, that are implicated in recruitment of the PP1 (yeast Glc7) phosphatase. In yeasts the RVxF/SILK domains have been shown to play a role in control of DNA replication initiation, at least in part through targeted de-phosphorylation of proteins in the pre-Replication Complex. In human cells Rif1 is recruited to DNA double-strand breaks through an interaction with 53BP1 where it counteracts DNA resection, thus promoting repair by non-homologous end-joining. This function requires the N-terminal HEAT repeat-containing domain. Interestingly, this domain is also implicated in DNA end protection at un-capped telomeres in yeast. We conclude by discussing the deployment of Rif1 at telomeres in yeasts from both an evolutionary perspective and in light of its recently discovered global functions.

Spatially coordinated replication and minimization of expression noise constrain three-dimensional organization of yeast genome

Despite recent advances, the underlying functional constraints that shape the three-dimensional organization of eukaryotic genome are not entirely clear. Through comprehensive multivariate analyses of genome-wide datasets, we show that cis and trans interactions in yeast genome have significantly distinct functional associations. In particular, (i) the trans interactions are constrained by coordinated replication and co-varying mutation rates of early replicating domains through interactions among early origins, while cis interactions are constrained by coordination of late replication through interactions among late origins; (ii)cis and trans interactions exhibit differential preference for nucleosome occupancy; (iii)cis interactions are also constrained by the essentiality and co-fitness of interacting genes. Essential gene clusters associate with high average interaction frequency, relatively short-range interactions of low variance, and exhibit less fluctuations in chromatin conformation, marking a physically restrained state of engaged loci that, we suggest, is important to mitigate the epigenetic errors by restricting the spatial mobility of loci. Indeed, the genes with lower expression noise associate with relatively short-range interactions of lower variance and exhibit relatively higher average interaction frequency, a property that is conserved across Escherichia coli,yeast, and mESCs. Altogether, our observations highlight the coordination of replication and the minimization of expression noise, not necessarily co-expression of genes, as potent evolutionary constraints shaping the spatial organization of yeast genome.

Priming for tolerance and cohesion at replication forks

Genome duplication is coupled with DNA damage tolerance (DDT) and chromatin structural changes. Recently we reported that mutations in Primase subunits or factors that bridge Polα/Primase with the replicative helicase, Ctf4, caused abnormal usage of DDT pathways, negatively influenced sister chromatid cohesion (SCC), and associated with increased fork reversal.¹ We also found that cohesin, which is paradigmatic for SCC, facilitates recombination-mediated DDT. However, only the recombination defects of cohesin, but not of cohesion-defective Polα/Primase/Ctf4 mutants, were rescued by artificial tethering of sister chromatids. Genetic tests and electron microscopy analysis of replication intermediates made us propose that management of single-stranded DNA forming proximal to the fork is a critical determinant of chromosome and replication fork structure, and influences DDT pathway choice. Here we discuss the implications of our findings for understanding DDT regulation and cohesion establishment during replication, and outline directions to rationalize the relationship between these chromosome metabolism processes.

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.

Polymerase ζ Activity Is Linked to Replication Timing in Humans: Evidence from Mutational Signatures

Replication timing is an important determinant of germline mutation patterns, with a higher rate of point mutations in late replicating regions. Mechanisms underlying this association remain elusive. One of the suggested explanations is the activity of error-prone DNA polymerases in late-replicating regions. Polymerase zeta (pol ζ), an essential error-prone polymerase biased toward transversions, also has a tendency to produce dinucleotide mutations (DNMs), complex mutational events that simultaneously affect two adjacent nucleotides. Experimental studies have shown that pol ζ is strongly biased toward GC→AA/TT DNMs. Using primate divergence data, we show that the GC→AA/TT pol ζ mutational signature is the most frequent among DNMs, and its rate exceeds the mean rate of other DNM types by a factor of approximately 10. Unlike the overall rate of DNMs, the pol ζ signature drastically increases with the replication time in the human genome. Finally, the pol ζ signature is enriched in transcribed regions, and there is a strong prevalence of GC→TT over GC→AA DNMs on the nontemplate strand, indicating association with transcription. A recurrently occurring GC→TT DNM in HRAS and SOD1 genes causes the Costello syndrome and amyotrophic lateral sclerosis correspondently; we observe an approximately 1 kb long mutation hotspot enriched by transversions near these DNMs in both cases, suggesting a link between these diseases and pol ζ activity. This study uncovers the genomic preferences of pol ζ, shedding light on a novel cause of mutational heterogeneity along the genome.

Mechanism of DNA damage tolerance

DNA damage may compromise genome integrity and lead to cell death. Cells have evolved a variety of processes to respond to DNA damage including damage repair and tolerance mechanisms, as well as damage checkpoints. The DNA damage tolerance (DDT) pathway promotes the bypass of single-stranded DNA lesions encountered by DNA polymerases during DNA replication. This prevents the stalling of DNA replication. Two mechanistically distinct DDT branches have been characterized. One is translesion synthesis (TLS) in which a replicative DNA polymerase is temporarily replaced by a specialized TLS polymerase that has the ability to replicate across DNA lesions. TLS is mechanistically simple and straightforward, but it is intrinsically error-prone. The other is the error-free template switching (TS) mechanism in which the stalled nascent strand switches from the damaged template to the undamaged newly synthesized sister strand for extension past the lesion. Error-free TS is a complex but preferable process for bypassing DNA lesions. However, our current understanding of this pathway is sketchy. An increasing number of factors are being found to participate or regulate this important mechanism, which is the focus of this editorial.

Extensive de novo mutation rate variation between individuals and across the genome of Chlamydomonas reinhardtii

Describing the process of spontaneous mutation is fundamental for understanding the genetic basis of disease, the threat posed by declining population size in conservation biology, and much of evolutionary biology. Directly studying spontaneous mutation has been difficult, however, because new mutations are rare. Mutation accumulation (MA) experiments overcome this by allowing mutations to build up over many generations in the near absence of natural selection. Here, we sequenced the genomes of 85 MA lines derived from six genetically diverse strains of the green alga Chlamydomonas reinhardtii. We identified 6843 new mutations, more than any other study of spontaneous mutation. We observed sevenfold variation in the mutation rate among strains and that mutator genotypes arose, increasing the mutation rate approximately eightfold in some replicates. We also found evidence for fine-scale heterogeneity in the mutation rate, with certain sequence motifs mutating at much higher rates, and clusters of multiple mutations occurring at closely linked sites. There was little evidence, however, for mutation rate heterogeneity between chromosomes or over large genomic regions of 200 kbp. We generated a predictive model of the mutability of sites based on their genomic properties, including local GC content, gene expression level, and local sequence context. Our model accurately predicted the average mutation rate and natural levels of genetic diversity of sites across the genome. Notably, trinucleotides vary 17-fold in rate between the most and least mutable sites. Our results uncover a rich heterogeneity in the process of spontaneous mutation both among individuals and across the genome.

Synonymous Genetic Variation in Natural Isolates of Escherichia coli Does Not Predict Where Synonymous Substitutions Occur in a Long-Term Experiment

Synonymous genetic differences vary by more than 20-fold among genes in natural isolates of Escherichia coli. One hypothesis to explain this heterogeneity is that genes with high levels of synonymous variation mutate at higher rates than genes with low synonymous variation. If so, then one would expect to observe similar mutational patterns in evolution experiments. In fact, however, the pattern of synonymous substitutions in a long-term evolution experiment with E. coli does not support this hypothesis. In particular, the extent of synonymous variation across genes in that experiment does not reflect the variation observed in natural isolates of E. coli. Instead, gene length alone predicts with high accuracy the prevalence of synonymous changes in the experimental populations. We hypothesize that patterns of synonymous variation in natural E. coli populations are instead caused by differences across genomic regions in their effective population size that, in turn, reflect different histories of recombination, horizontal gene transfer, selection, and population structure.

Control of Caenorhabditis elegans germ-line stem-cell cycling speed meets requirements of design to minimize mutation accumulation

Background

Stem cells are thought to play a critical role in minimizing the accumulation of mutations, but it is not clear which strategies they follow to fulfill that performance objective. Slow cycling of stem cells provides a simple strategy that can minimize cell pedigree depth and thereby minimize the accumulation of replication-dependent mutations. Although the power of this strategy was recognized early on, a quantitative assessment of whether and how it is employed by biological systems is missing.

Results

Here we address this problem using a simple self-renewing organ – the C. elegans gonad – whose overall organization is shared with many self-renewing organs. Computational simulations of mutation accumulation characterize a tradeoff between fast development and low mutation accumulation, and show that slow-cycling stem cells allow for an advantageous compromise to be reached. This compromise is such that worm germ-line stem cells should cycle more slowly than their differentiating counterparts, but only by a modest amount. Experimental measurements of cell cycle lengths derived using a new, quantitative technique are consistent with these predictions.

Conclusions

Our findings shed light both on design principles that underlie the role of stem cells in delaying aging and on evolutionary forces that shape stem-cell gene regulatory networks.

Electronic supplementary material

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

The Rate and Molecular Spectrum of Spontaneous Mutations in the GC-Rich Multichromosome Genome of Burkholderia cenocepacia

Spontaneous mutations are ultimately essential for evolutionary change and are also the root cause of many diseases. However, until recently, both biological and technical barriers have prevented detailed analyses of mutation profiles, constraining our understanding of the mutation process to a few model organisms and leaving major gaps in our understanding of the role of genome content and structure on mutation. Here, we present a genome-wide view of the molecular mutation spectrum in Burkholderia cenocepacia, a clinically relevant pathogen with high %GC content and multiple chromosomes. We find that B. cenocepacia has low genome-wide mutation rates with insertion-deletion mutations biased toward deletions, consistent with the idea that deletion pressure reduces prokaryotic genome sizes. Unlike prior studies of other organisms, mutations in B. cenocepacia are not AT biased, which suggests that at least some genomes with high %GC content experience unusual base-substitution mutation pressure. Importantly, we also observe variation in both the rates and spectra of mutations among chromosomes and elevated G:C > T:A transversions in late-replicating regions. Thus, although some patterns of mutation appear to be highly conserved across cellular life, others vary between species and even between chromosomes of the same species, potentially influencing the evolution of nucleotide composition and genome architecture.

Dissecting genetic and environmental mutation signatures with model organisms

Deep sequencing has impacted on cancer research by enabling routine sequencing of genomes and exomes to identify genetic changes associated with carcinogenesis. Researchers can now use the frequency, type, and context of all mutations in tumor genomes to extract mutation signatures that reflect the driving mutational processes. Identifying mutation signatures, however, may not immediately suggest a mechanism. Consequently, several recent studies have employed deep sequencing of model organisms exposed to discrete genetic or environmental perturbations. These studies exploit the simpler genomes and availability of powerful genetic tools in model organisms to analyze mutation signatures under controlled conditions, forging mechanistic links between mutational processes and signatures. We discuss the power of this approach and suggest that many such studies may be on the horizon.

The spatial and temporal organization of origin firing during the S-phase of fission yeast

Eukaryotes duplicate their genomes using multiple replication origins, but the organization of origin firing along chromosomes and during S-phase is not well understood. Using fission yeast, we report the first genome-wide analysis of the spatial and temporal organization of replication origin firing, analyzed using single DNA molecules that can approach the full length of chromosomes. At S-phase onset, origins fire randomly and sparsely throughout the chromosomes. Later in S-phase, clusters of fired origins appear embedded in the sparser regions, which form the basis of nuclear replication foci. The formation of clusters requires proper histone methylation and acetylation, and their locations are not inherited between cell cycles. The rate of origin firing increases gradually, peaking just before mid S-phase. Toward the end of S-phase, nearly all the available origins within the unreplicated regions are fired, contributing to the timely completion of genome replication. We propose that the majority of origins do not fire as a part of a deterministic program. Instead, origin firing, both individually and as clusters, should be viewed as being mostly stochastic.

A global profile of replicative polymerase usage

Three eukaryotic DNA polymerases are essential for genome replication. Polymerase (Pol) α-primase initiates each synthesis event and is rapidly replaced by processive DNA polymerases: Polɛ replicates the leading strand, whereas Polδ performs lagging-strand synthesis. However, it is not known whether this division of labor is maintained across the whole genome or how uniform it is within single replicons. Using Schizosaccharomyces pombe, we have developed a polymerase usage sequencing (Pu-seq) strategy to map polymerase usage genome wide. Pu-seq provides direct replication-origin location and efficiency data and indirect estimates of replication timing. We confirm that the division of labor is broadly maintained across an entire genome. However, our data suggest a subtle variability in the usage of the two polymerases within individual replicons. We propose that this results from occasional leading-strand initiation by Polδ followed by exchange for Polɛ.

Elevated mutation rate during meiosis in Saccharomyces cerevisiae

Mutations accumulate during all stages of growth, but only germ line mutations contribute to evolution. While meiosis contributes to evolution by reassortment of parental alleles, we show here that the process itself is inherently mutagenic. We have previously shown that the DNA synthesis associated with repair of a double-strand break is about 1000-fold less accurate than S-phase synthesis. Since the process of meiosis involves many programmed DSBs, we reasoned that this repair might also be mutagenic. Indeed, in the early 1960's Magni and Von Borstel observed elevated reversion of recessive alleles during meiosis, and found that the revertants were more likely to be associated with a crossover than non-revertants, a process that they called "the meiotic effect." Here we use a forward mutation reporter (CAN1 HIS3) placed at either a meiotic recombination coldspot or hotspot near the MAT locus on Chromosome III. We find that the increased mutation rate at CAN1 (6 to 21 -fold) correlates with the underlying recombination rate at the locus. Importantly, we show that the elevated mutation rate is fully dependent upon Spo11, the protein that introduces the meiosis specific DSBs. To examine associated recombination we selected for random spores with or without a mutation in CAN1. We find that the mutations isolated this way show an increased association with recombination (crossovers, loss of crossover interference and/or increased gene conversion tracts). Polζ appears to contribute about half of the mutations induced during meiosis, but is not the only source of mutations for the meiotic effect. We see no difference in either the spectrum or distribution of mutations between mitosis and meiosis. The correlation of hotspots with elevated mutagenesis provides a mechanism for organisms to control evolution rates in a gene specific manner.

Mutation rate, spectrum, topology, and context-dependency in the DNA mismatch repair-deficient Pseudomonas fluorescens ATCC948

High levels of genetic diversity exist among natural isolates of the bacterium Pseudomonas fluorescens, and are especially elevated around the replication terminus of the genome, where strain-specific genes are found. In an effort to understand the role of genetic variation in the evolution of Pseudomonas, we analyzed 31,106 base substitutions from 45 mutation accumulation lines of P. fluorescens ATCC948, naturally deficient for mismatch repair, yielding a base-substitution mutation rate of 2.34 × 10⁻⁸ per site per generation (SE: 0.01 × 10⁻⁸) and a small-insertion-deletion mutation rate of 1.65 × 10⁻⁹ per site per generation (SE: 0.03 × 10⁻⁹). We find that the spectrum of mutations in prophage regions, which often contain virulence factors and antibiotic resistance, is highly similar to that in the intergenic regions of the host genome. Our results show that the mutation rate varies around the chromosome, with the lowest mutation rate found near the origin of replication. Consistent with observations from other studies, we find that site-specific mutation rates are heavily influenced by the immediately flanking nucleotides, indicating that mutations are context dependent.

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.

Heterogeneous polymerase fidelity and mismatch repair bias genome variation and composition

Mutational heterogeneity must be taken into account when reconstructing evolutionary histories, calibrating molecular clocks, and predicting links between genes and disease. Selective pressures and various DNA transactions have been invoked to explain the heterogeneous distribution of genetic variation between species, within populations, and in tissue-specific tumors. To examine relationships between such heterogeneity and variations in leading- and lagging-strand replication fidelity and mismatch repair, we accumulated 40,000 spontaneous mutations in eight diploid yeast strains in the absence of selective pressure. We found that replicase error rates vary by fork direction, coding state, nucleosome proximity, and sequence context. Further, error rates and DNA mismatch repair efficiency both vary by mismatch type, responsible polymerase, replication time, and replication origin proximity. Mutation patterns implicate replication infidelity as one driver of variation in somatic and germline evolution, suggest mechanisms of mutual modulation of genome stability and composition, and predict future observations in specific cancers.

Precise estimates of mutation rate and spectrum in yeast

Mutation is the ultimate source of genetic variation. The most direct and unbiased method of studying spontaneous mutations is via mutation accumulation (MA) lines. Until recently, MA experiments were limited by the cost of sequencing and thus provided us with small numbers of mutational events and therefore imprecise estimates of rates and patterns of mutation. We used whole-genome sequencing to identify nearly 1,000 spontaneous mutation events accumulated over ∼311,000 generations in 145 diploid MA lines of the budding yeast Saccharomyces cerevisiae. MA experiments are usually assumed to have negligible levels of selection, but even mild selection will remove strongly deleterious events. We take advantage of such patterns of selection and show that mutation classes such as indels and aneuploidies (especially monosomies) are proportionately much more likely to contribute mutations of large effect. We also provide conservative estimates of indel, aneuploidy, environment-dependent dominant lethal, and recessive lethal mutation rates. To our knowledge, for the first time in yeast MA data, we identified a sufficiently large number of single-nucleotide mutations to measure context-dependent mutation rates and were able to (i) confirm strong AT bias of mutation in yeast driven by high rate of mutations from C/G to T/A and (ii) detect a higher rate of mutation at C/G nucleotides in two specific contexts consistent with cytosine methylation in S. cerevisiae.

Complex correlations: replication timing and mutational landscapes during cancer and genome evolution

A recent flurry of reports correlates replication timing (RT) with mutation rates during both evolution and cancer. Specifically, point mutations and copy number losses correlate with late replication, while copy number gains and other rearrangements correlate with early replication. In some cases, plausible mechanisms have been proposed. Point mutation rates may reflect temporal variation in repair mechanisms. Transcription-induced double-strand breaks are expected to occur in transcriptionally active early replicating chromatin. Fusion partners are generally in close proximity, and chromatin in close proximity replicates at similar times. However, temporal enrichment of copy number gains and losses remains an enigma. Moreover, many conclusions are compromised by a lack of matched RT and sequence datasets, the filtering out of developmental variation in RT, and the use of somatic cell lines to make inferences about germline evolution.

DNA bending facilitates the error-free DNA damage tolerance pathway and upholds genome integrity

DNA replication is sensitive to damage in the template. To bypass lesions and complete replication, cells activate recombination-mediated (error-free) and translesion synthesis-mediated (error-prone) DNA damage tolerance pathways. Crucial for error-free DNA damage tolerance is template switching, which depends on the formation and resolution of damage-bypass intermediates consisting of sister chromatid junctions. Here we show that a chromatin architectural pathway involving the high mobility group box protein Hmo1 channels replication-associated lesions into the error-free DNA damage tolerance pathway mediated by Rad5 and PCNA polyubiquitylation, while preventing mutagenic bypass and toxic recombination. In the process of template switching, Hmo1 also promotes sister chromatid junction formation predominantly during replication. Its C-terminal tail, implicated in chromatin bending, facilitates the formation of catenations/hemicatenations and mediates the roles of Hmo1 in DNA damage tolerance pathway choice and sister chromatid junction formation. Together, the results suggest that replication-associated topological changes involving the molecular DNA bender, Hmo1, set the stage for dedicated repair reactions that limit errors during replication and impact on genome stability.

Late-replicating CNVs as a source of new genes

Asynchronous replication of the genome has been associated with different rates of point mutation and copy number variation (CNV) in human populations. Here, our aim was to investigate whether the bias in the generation of CNV that is associated with DNA replication timing might have conditioned the birth of new protein-coding genes during evolution. We show that genes that were duplicated during primate evolution are more commonly found among the human genes located in late-replicating CNV regions. We traced the relationship between replication timing and the evolutionary age of duplicated genes. Strikingly, we found that there is a significant enrichment of evolutionary younger duplicates in late-replicating regions of the human and mouse genome. Indeed, the presence of duplicates in late-replicating regions gradually decreases as the evolutionary time since duplication extends. Our results suggest that the accumulation of recent duplications in late-replicating CNV regions is an active process influencing genome evolution.

Genome destabilizing mutator alleles drive specific mutational trajectories in Saccharomyces cerevisiae

In addition to environmental factors and intrinsic variations in base substitution rates, specific genome-destabilizing mutations can shape the mutational trajectory of genomes. How specific alleles influence the nature and position of accumulated mutations in a genomic context is largely unknown. Understanding the impact of genome-destabilizing alleles is particularly relevant to cancer genomes where biased mutational signatures are identifiable. We first created a more complete picture of cellular pathways that impact mutation rate using a primary screen to identify essential Saccharomyces cerevisiae gene mutations that cause mutator phenotypes. Drawing primarily on new alleles identified in this resource, we measure the impact of diverse mutator alleles on mutation patterns directly by whole-genome sequencing of 68 mutation-accumulation strains derived from wild-type and 11 parental mutator genotypes. The accumulated mutations differ across mutator strains, displaying base-substitution biases, allele-specific mutation hotspots, and break-associated mutation clustering. For example, in mutants of POLα and the Cdc13–Stn1–Ten1 complex, we find a distinct subtelomeric bias for mutations that we show is independent of the target sequence. Together our data suggest that specific genome-instability mutations are sufficient to drive discrete mutational signatures, some of which share properties with mutation patterns seen in tumors. Thus, in a population of cells, genome-instability mutations could influence clonal evolution by establishing discrete mutational trajectories for genomes.

DNA repair mechanisms and the bypass of DNA damage in Saccharomyces cerevisiae

DNA repair mechanisms are critical for maintaining the integrity of genomic DNA, and their loss is associated with cancer predisposition syndromes. Studies in Saccharomyces cerevisiae have played a central role in elucidating the highly conserved mechanisms that promote eukaryotic genome stability. This review will focus on repair mechanisms that involve excision of a single strand from duplex DNA with the intact, complementary strand serving as a template to fill the resulting gap. These mechanisms are of two general types: those that remove damage from DNA and those that repair errors made during DNA synthesis. The major DNA-damage repair pathways are base excision repair and nucleotide excision repair, which, in the most simple terms, are distinguished by the extent of single-strand DNA removed together with the lesion. Mistakes made by DNA polymerases are corrected by the mismatch repair pathway, which also corrects mismatches generated when single strands of non-identical duplexes are exchanged during homologous recombination. In addition to the true repair pathways, the postreplication repair pathway allows lesions or structural aberrations that block replicative DNA polymerases to be tolerated. There are two bypass mechanisms: an error-free mechanism that involves a switch to an undamaged template for synthesis past the lesion and an error-prone mechanism that utilizes specialized translesion synthesis DNA polymerases to directly synthesize DNA across the lesion. A high level of functional redundancy exists among the pathways that deal with lesions, which minimizes the detrimental effects of endogenous and exogenous DNA damage.

Three wise centromere functions: see no error, hear no break, speak no delay

The main function of the centromere is to promote kinetochore assembly for spindle microtubule attachment. Two additional functions of the centromere, however, are becoming increasingly clear: facilitation of robust sister-chromatid cohesion at pericentromeres and advancement of replication of centromeric regions. The combination of these three centromere functions ensures correct chromosome segregation during mitosis. Here, we review the mechanisms of the kinetochore-microtubule interaction, focusing on sister-kinetochore bi-orientation (or chromosome bi-orientation). We also discuss the biological importance of robust pericentromeric cohesion and early centromere replication, as well as the mechanisms orchestrating these two functions at the microtubule attachment site.

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

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

•

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

•

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

•

Origin activity determines the location of replication termination events

•

Termination events are widely distributed across the whole genome

Replisome stall events have shaped the distribution of replication origins in the genomes of yeasts

During S phase, the entire genome must be precisely duplicated, with no sections of DNA left unreplicated. Here, we develop a simple mathematical model to describe the probability of replication failing due to the irreversible stalling of replication forks. We show that the probability of complete genome replication is maximized if replication origins are evenly spaced, the largest inter-origin distances are minimized, and the end-most origins are positioned close to chromosome ends. We show that origin positions in the yeast Saccharomyces cerevisiae genome conform to all three predictions thereby maximizing the probability of complete replication if replication forks stall. Origin positions in four other yeasts-Kluyveromyces lactis, Lachancea kluyveri, Lachancea waltii and Schizosaccharomyces pombe-also conform to these predictions. Equating failure rates at chromosome ends with those in chromosome interiors gives a mean per nucleotide fork stall rate of ∼5 × 10(-8), which is consistent with experimental estimates. Using this value in our theoretical predictions gives replication failure rates that are consistent with data from replication origin knockout experiments. Our theory also predicts that significantly larger genomes, such as those of mammals, will experience a much greater probability of replication failure genome-wide, and therefore will likely require additional compensatory mechanisms.

A Link between ORC-origin binding mechanisms and origin activation time revealed in budding yeast

Eukaryotic DNA replication origins are selected in G1-phase when the origin recognition complex (ORC) binds chromosomal positions and triggers molecular events culminating in the initiation of DNA replication (a.k.a. origin firing) during S-phase. Each chromosome uses multiple origins for its duplication, and each origin fires at a characteristic time during S-phase, creating a cell-type specific genome replication pattern relevant to differentiation and genome stability. It is unclear whether ORC-origin interactions are relevant to origin activation time. We applied a novel genome-wide strategy to classify origins in the model eukaryote Saccharomyces cerevisiae based on the types of molecular interactions used for ORC-origin binding. Specifically, origins were classified as DNA-dependent when the strength of ORC-origin binding in vivo could be explained by the affinity of ORC for origin DNA in vitro, and, conversely, as ‘chromatin-dependent’ when the ORC-DNA interaction in vitro was insufficient to explain the strength of ORC-origin binding in vivo. These two origin classes differed in terms of nucleosome architecture and dependence on origin-flanking sequences in plasmid replication assays, consistent with local features of chromatin promoting ORC binding at ‘chromatin-dependent’ origins. Finally, the ‘chromatin-dependent’ class was enriched for origins that fire early in S-phase, while the DNA-dependent class was enriched for later firing origins. Conversely, the latest firing origins showed a positive association with the ORC-origin DNA paradigm for normal levels of ORC binding, whereas the earliest firing origins did not. These data reveal a novel association between ORC-origin binding mechanisms and the regulation of origin activation time.

Cell division requires the duplication of chromosomes, protein-DNA complexes harboring genetic information. Specific chromosomal positions, origins, initiate this duplication. Multiple origins are required for accurate, efficient duplication—an insufficient number leads to mistakes in the genetic material and pathologies such as cancer. Origins are chosen when the origin recognition complex (ORC) binds to them. The molecular interactions controlling this binding remain unclear. Understanding these interactions will lead to new ways to control cell division, which could aid in treatments of disease. Experiments were performed in the eukaryotic microbe budding yeast to define the types of molecular interactions ORC uses to bind origins. Yeasts are useful for these studies because chromosome duplication and structure are well conserved from yeast to humans. While ORC-DNA interactions were important, interactions between ORC and chromosomal proteins played a role. In addition, different origins relied on different types of molecular interactions with ORC. Finally, ORC-protein interactions but not ORC-DNA interactions were associated with enhanced origin function during chromosome-duplication, revealing an unanticipated link between the types of molecular interactions ORC uses to select an origin and the ultimate function of that origin. These results have implications for interfering with ORC-origin interactions to control cell division.

Mutation rates, spectra, and genome-wide distribution of spontaneous mutations in mismatch repair deficient yeast

DNA mismatch repair is a highly conserved DNA repair pathway. In humans, germline mutations in hMSH2 or hMLH1, key components of mismatch repair, have been associated with Lynch syndrome, a leading cause of inherited cancer mortality. Current estimates of the mutation rate and the mutational spectra in mismatch repair defective cells are primarily limited to a small number of individual reporter loci. Here we use the yeast Saccharomyces cerevisiae to generate a genome-wide view of the rates, spectra, and distribution of mutation in the absence of mismatch repair. We performed mutation accumulation assays and next generation sequencing on 19 strains, including 16 msh2 missense variants implicated in Lynch cancer syndrome. The mutation rate for DNA mismatch repair null strains was approximately 1 mutation per genome per generation, 225-fold greater than the wild-type rate. The mutations were distributed randomly throughout the genome, independent of replication timing. The mutation spectra included insertions/deletions at homopolymeric runs (87.7%) and at larger microsatellites (5.9%), as well as transitions (4.5%) and transversions (1.9%). Additionally, repeat regions with proximal repeats are more likely to be mutated. A bias toward deletions at homopolymers and insertions at (AT)_n microsatellites suggests a different mechanism for mismatch generation at these sites. Interestingly, 5% of the single base pair substitutions might represent double-slippage events that occurred at the junction of immediately adjacent repeats, resulting in a shift in the repeat boundary. These data suggest a closer scrutiny of tumor suppressors with homopolymeric runs with proximal repeats as the potential drivers of oncogenesis in mismatch repair defective cells.

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.

The yeast environmental stress response regulates mutagenesis induced by proteotoxic stress

Conditions of chronic stress are associated with genetic instability in many organisms, but the roles of stress responses in mutagenesis have so far been elucidated only in bacteria. Here, we present data demonstrating that the environmental stress response (ESR) in yeast functions in mutagenesis induced by proteotoxic stress. We show that the drug canavanine causes proteotoxic stress, activates the ESR, and induces mutagenesis at several loci in an ESR-dependent manner. Canavanine-induced mutagenesis also involves translesion DNA polymerases Rev1 and Polζ and non-homologous end joining factor Ku. Furthermore, under conditions of chronic sub-lethal canavanine stress, deletions of Rev1, Polζ, and Ku-encoding genes exhibit genetic interactions with ESR mutants indicative of ESR regulating these mutagenic DNA repair processes. Analyses of mutagenesis induced by several different stresses showed that the ESR specifically modulates mutagenesis induced by proteotoxic stress. Together, these results document the first known example of an involvement of a eukaryotic stress response pathway in mutagenesis and have important implications for mechanisms of evolution, carcinogenesis, and emergence of drug-resistant pathogens and chemotherapy-resistant tumors.

Cellular capability to mutate its DNA plays an important role in evolution and impinges on medical issues, including acquisition of mutator phenotypes by cancer cells and emergence of drug-resistant pathogens. Whether and how the environment affects rates of mutation has been studied predominantly in the context of environmental agents that damage DNA (e.g. UV and γ-rays). However, it has been observed that conditions of chronic non-DNA-damaging stress (e.g. starvation or heat shock) also increase mutagenesis. It has been shown that in bacteria, activation of the general stress response activates a pro-mutagenic pathway and thus promotes mutagenesis during periods of stress. However, in eukaryotes, so far there has been no evidence of a stress response regulating mutagenesis. In this manuscript we demonstrate that in budding yeast, a model eukaryote, the general environmental stress response (ESR) regulates mutagenesis induced by proteotoxic stress (accumulation of unfolded proteins) at several loci. We also identify two pro-mutagenic DNA metabolic pathways that contribute to this mutagenesis and present genetic data showing that the ESR regulates these pathways. Together, these data advance our understanding of how cellular sensing and responding to environmental cues affect cellular capability for mutagenesis.

Estimating selection on synonymous codon usage from noisy experimental data

A key goal in molecular evolution is to extract mechanistic insights from signatures of selection. A case study is codon usage, where despite many recent advances and hypotheses, two longstanding problems remain: the relative contribution of selection and mutation in determining codon frequencies and the relative contribution of translational speed and accuracy to selection. The relevant targets of selection--the rate of translation and of mistranslation of a codon per unit time in the cell--can only be related to mechanistic properties of the translational apparatus if the number of transcripts per cell is known, requiring use of gene expression measurements. Perhaps surprisingly, different gene-expression data sets yield markedly different estimates of selection. We show that this is largely due to measurement noise, notably due to differences between studies rather than instrument error or biological variability. We develop an analytical framework that explicitly models noise in expression in the context of the population-genetic model. Estimates of mutation and selection strength in budding yeast produced by this method are robust to the expression data set used and are substantially higher than estimates using a noise-blind approach. We introduce per-gene selection estimates that correlate well with previous scoring systems, such as the codon adaptation index, while now carrying an evolutionary interpretation. On average, selection for codon usage in budding yeast is weak, yet our estimates show that genes range from virtually unselected to average per-codon selection coefficients above the inverse population size. Our analytical framework may be generally useful for distinguishing biological signals from measurement noise in other applications that depend upon measurements of gene expression.

DNA replication timing, genome stability and cancer: late and/or delayed DNA replication timing is associated with increased genomic instability

Normal cellular division requires that the genome be faithfully replicated to ensure that unaltered genomic information is passed from one generation to the next. DNA replication initiates from thousands of origins scattered throughout the genome every cell cycle; however, not all origins initiate replication at the same time. A vast amount of work over the years indicates that different origins along each eukaryotic chromosome are activated in early, middle or late S phase. This temporal control of DNA replication is referred to as the replication-timing program. The replication-timing program represents a very stable epigenetic feature of chromosomes. Recent evidence has indicated that the replication-timing program can influence the spatial distribution of mutagenic events such that certain regions of the genome experience increased spontaneous mutagenesis compared to surrounding regions. This influence has helped shape the genomes of humans and other multicellular organisms and can affect the distribution of mutations in somatic cells. It is also becoming clear that the replication-timing program is deregulated in many disease states, including cancer. Aberrant DNA replication timing is associated with changes in gene expression, changes in epigenetic modifications and an increased frequency of structural rearrangements. Furthermore, certain replication timing changes can directly lead to overt genomic instability and may explain unique mutational signatures that are present in cells that have undergone the recently described processes of "chromothripsis" and "kataegis". In this review, we will discuss how the normal replication timing program, as well as how alterations to this program, can contribute to the evolution of the genomic landscape in normal and cancerous cells.