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

URA6 mutations provide an alternative mechanism for 5-FOA resistance in Saccharomyces cerevisiae

URA3 is frequently used in the yeast community as the mutation target for 5-fluoroorotic acid (5-FOA) resistance. We identified a novel class of ura6 mutants that can grow in the presence of 5-FOA. Unlike ura3 mutants, ura6 mutants remain prototrophic and grow in the absence of uracil. In addition to 5-FOA resistance, we found that mutations to URA6 also confer resistance to 5-fluorocytosine (5-FC) and 5-fluorouracil (5-FU). In total, we identified 50 unique missense mutations across 32 residues of URA6. We found that 28 out of the 32 affected residues are located in regions conserved between Saccharomyces cerevisiae and three clinically relevant pathogenic fungi. These findings suggest that mutations to URA6 present a second target for mutation screens utilizing 5-FOA as a selection marker as well as a potential mode of resistance to the antifungal therapeutic 5-FC.

Mutational robustness and the role of buffer genes in evolvability

Organisms rely on mutations to fuel adaptive evolution. However, many mutations impose a negative effect on fitness. Cells may have therefore evolved mechanisms that affect the phenotypic effects of mutations, thus conferring mutational robustness. Specifically, so-called buffer genes are hypothesized to interact directly or indirectly with genetic variation and reduce its effect on fitness. Environmental or genetic perturbations can change the interaction between buffer genes and genetic variation, thereby unmasking the genetic variation's phenotypic effects and thus providing a source of variation for natural selection to act on. This review provides an overview of our understanding of mutational robustness and buffer genes, with the chaperone gene HSP90 as a key example. It discusses whether buffer genes merely affect standing variation or also interact with de novo mutations, how mutational robustness could influence evolution, and whether mutational robustness might be an evolved trait or rather a mere side-effect of complex genetic interactions.

Emergence of enhancers at late DNA replicating regions

Enhancers are fast-evolving genomic sequences that control spatiotemporal gene expression patterns. By examining enhancer turnover across mammalian species and in multiple tissue types, we uncover a relationship between the emergence of enhancers and genome organization as a function of germline DNA replication time. While enhancers are most abundant in euchromatic regions, enhancers emerge almost twice as often in late compared to early germline replicating regions, independent of transposable elements. Using a deep learning sequence model, we demonstrate that new enhancers are enriched for mutations that alter transcription factor (TF) binding. Recently evolved enhancers appear to be mostly neutrally evolving and enriched in eQTLs. They also show more tissue specificity than conserved enhancers, and the TFs that bind to these elements, as inferred by binding sequences, also show increased tissue-specific gene expression. We find a similar relationship with DNA replication time in cancer, suggesting that these observations may be time-invariant principles of genome evolution. Our work underscores that genome organization has a profound impact in shaping mammalian gene regulation.

DNA Damage, Genome Stability, and Adaptation: A Question of Chance or Necessity?

DNA damage causes the mutations that are the principal source of genetic variation. DNA damage detection and repair mechanisms therefore play a determining role in generating the genetic diversity on which natural selection acts. Speciation, it is commonly assumed, occurs at a rate set by the level of standing allelic diversity in a population. The process of speciation is driven by a combination of two evolutionary forces: genetic drift and ecological selection. Genetic drift takes place under the conditions of relaxed selection, and results in a balance between the rates of mutation and the rates of genetic substitution. These two processes, drift and selection, are necessarily mediated by a variety of mechanisms guaranteeing genome stability in any given species. One of the outstanding questions in evolutionary biology concerns the origin of the widely varying phylogenetic distribution of biodiversity across the Tree of Life and how the forces of drift and selection contribute to shaping that distribution. The following examines some of the molecular mechanisms underlying genome stability and the adaptive radiations that are associated with biodiversity and the widely varying species richness and evenness in the different eukaryotic lineages.

Transcription as source of genetic heterogeneity in budding yeast

Transcription presents challenges to genome stability both directly, by altering genome topology and exposing single-stranded DNA to chemical insults and nucleases, and indirectly by introducing obstacles to the DNA replication machinery. Such obstacles include the RNA polymerase holoenzyme itself, DNA-bound regulatory factors, G-quadruplexes and RNA-DNA hybrid structures known as R-loops. Here, we review the detrimental impacts of transcription on genome stability in budding yeast, as well as the mitigating effects of transcription-coupled nucleotide excision repair and of systems that maintain DNA replication fork processivity and integrity. Interactions between DNA replication and transcription have particular potential to induce mutation and structural variation, but we conclude that such interactions must have only minor effects on DNA replication by the replisome with little if any direct mutagenic outcome. However, transcription can significantly impair the fidelity of replication fork rescue mechanisms, particularly Break Induced Replication, which is used to restart collapsed replication forks when other means fail. This leads to de novo mutations, structural variation and extrachromosomal circular DNA formation that contribute to genetic heterogeneity, but only under particular conditions and in particular genetic contexts, ensuring that the bulk of the genome remains extremely stable despite the seemingly frequent interactions between transcription and DNA replication.

Enhancing cold resistance in Banana (Musa spp.) through EMS-induced mutagenesis, L-Hyp pressure selection: phenotypic alterations, biomass composition, and transcriptomic insights

BACKGROUND

The cultivation of bananas encounters substantial obstacles, particularly due to the detrimental effects of cold stress on their growth and productivity. A potential remedy that has gained attention is the utilization of ethyl mesylate (EMS)-induced mutagenesis technology, which enables the creation of a genetically varied group of banana mutants. This complex procedure entails subjecting the mutants to further stress screening utilizing L-Hyp in order to identify those exhibiting improved resistance to cold. This study conducted a comprehensive optimization of the screening conditions for EMS mutagenesis and L-Hyp, resulting in the identification of the mutant cm784, which exhibited remarkable cold resistance. Subsequent investigations further elucidated the physiological and transcriptomic responses of cm784 to low-temperature stress.

RESULTS

EMS mutagenesis had a substantial effect on banana seedlings, resulting in modifications in shoot and root traits, wherein a majority of seedlings exhibited delayed differentiation and limited elongation. Notably, mutant leaves displayed altered biomass composition, with starch content exhibiting the most pronounced variation. The application of L-Hyp pressure selection aided in the identification of cold-resistant mutants among seedling-lethal phenotypes. The mutant cm784 demonstrated enhanced cold resistance, as evidenced by improved survival rates and reduced symptoms of chilling injury. Physiological analyses demonstrated heightened activities of antioxidant enzymes and increased proline production in cm784 when subjected to cold stress. Transcriptome analysis unveiled 946 genes that were differentially expressed in cm784, with a notable enrichment in categories related to 'Carbohydrate transport and metabolism' and 'Secondary metabolites biosynthesis, transport, and catabolism'.

CONCLUSION

The present findings provide insights into the molecular mechanisms that contribute to the heightened cold resistance observed in banana mutants. These mechanisms encompass enhanced carbohydrate metabolism and secondary metabolite biosynthesis, thereby emphasizing the adaptive strategies employed to mitigate the detrimental effects induced by cold stress.

RIF1 regulates early replication timing in murine B cells

The mammalian DNA replication timing (RT) program is crucial for the proper functioning and integrity of the genome. The best-known mechanism for controlling RT is the suppression of late origins of replication in heterochromatin by RIF1. Here, we report that in antigen-activated, hypermutating murine B lymphocytes, RIF1 binds predominantly to early-replicating active chromatin and promotes early replication, but plays a minor role in regulating replication origin activity, gene expression and genome organization in B cells. Furthermore, we find that RIF1 functions in a complementary and non-epistatic manner with minichromosome maintenance (MCM) proteins to establish early RT signatures genome-wide and, specifically, to ensure the early replication of highly transcribed genes. These findings reveal additional layers of regulation within the B cell RT program, driven by the coordinated activity of RIF1 and MCM proteins.

Upper Bound on the Mutational Burden Imposed by a CRISPR-Cas9 Gene-Drive Element

CRISPR-Cas9 gene drives (CCGDs) are powerful tools for genetic control of wild populations, useful for eradication of disease vectors, conservation of endangered species and other applications. However, Cas9 alone and in a complex with gRNA can cause double-stranded DNA breaks at off-target sites, which could increase the mutational load and lead to loss of heterozygosity (LOH). These undesired effects raise potential concerns about the long-term evolutionary safety of CCGDs, but the magnitude of these effects is unknown. To estimate how the presence of a CCGD or a Cas9 alone in the genome affects the rates of LOH events and de novo mutations, we carried out a mutation accumulation experiment in yeast Saccharomyces cerevisiae. Despite its substantial statistical power, our experiment revealed no detectable effect of CCGD or Cas9 alone on the genome-wide rates of mutations or LOH events, suggesting that these rates are affected by less than 30%. Nevertheless, we found that Cas9 caused a slight but significant shift towards more interstitial and fewer terminal LOH events, and the CCGD caused a significant difference in the distribution of LOH events on Chromosome V. Taken together, our results show that these genetic elements impose a weak and likely localized additional mutational burden in the yeast model. Although the mutagenic effects of CCGDs need to be further evaluated in other systems, our results suggest that the effect of CCGDs on off-target mutation rates and genetic diversity may be acceptable.

Two independent DNA repair pathways cause mutagenesis in template switching deficient Saccharomyces cerevisiae

Upon DNA replication stress, cells utilize the postreplication repair pathway to repair single-stranded DNA and maintain genome integrity. Postreplication repair is divided into 2 branches: error-prone translesion synthesis, signaled by proliferating cell nuclear antigen (PCNA) monoubiquitination, and error-free template switching, signaled by PCNA polyubiquitination. In Saccharomyces cerevisiae, Rad5 is involved in both branches of repair during DNA replication stress. When the PCNA polyubiquitination function of Rad5 s disrupted, Rad5 recruits translesion synthesis polymerases to stalled replication forks, resulting in mutagenic repair. Details of how mutagenic repair is carried out, as well as the relationship between Rad5-mediated mutagenic repair and the canonical PCNA-mediated mutagenic repair, remain to be understood. We find that Rad5-mediated mutagenic repair requires the translesion synthesis polymerase ζ but does not require other yeast translesion polymerase activities. Furthermore, we show that Rad5-mediated mutagenic repair is independent of PCNA binding by Rev1 and so is separable from canonical mutagenic repair. In the absence of error-free template switching, both modes of mutagenic repair contribute additively to replication stress response in a replication timing-independent manner. Cellular contexts where error-free template switching is compromised are not simply laboratory phenomena, as we find that a natural variant in RAD5 is defective in PCNA polyubiquitination and therefore defective in error-free repair, resulting in Rad5- and PCNA-mediated mutagenic repair. Our results highlight the importance of Rad5 in regulating spontaneous mutagenesis and genetic diversity in S. cerevisiae through different modes of postreplication repair.

The divergence of mutation rates and spectra across the Tree of Life

Owing to advances in genome sequencing, genome stability has become one of the most scrutinized cellular traits across the Tree of Life. Despite its centrality to all things biological, the mutation rate (per nucleotide site per generation) ranges over three orders of magnitude among species and several-fold within individual phylogenetic lineages. Within all major organismal groups, mutation rates scale negatively with the effective population size of a species and with the amount of functional DNA in the genome. This relationship is most parsimoniously explained by the drift-barrier hypothesis, which postulates that natural selection typically operates to reduce mutation rates until further improvement is thwarted by the power of random genetic drift. Despite this constraint, the molecular mechanisms underlying DNA replication fidelity and repair are free to wander, provided the performance of the entire system is maintained at the prevailing level. The evolutionary flexibility of the mutation rate bears on the resolution of several prior conundrums in phylogenetic and population-genetic analysis and raises challenges for future applications in these areas.

Environmental carcinogens disproportionally mutate genes implicated in neurodevelopmental disorders

Introduction

De novo mutations contribute to a large proportion of sporadic psychiatric and developmental disorders, yet the potential role of environmental carcinogens as drivers of causal de novo mutations in neurodevelopmental disorders is poorly studied.

Methods

To explore environmental mutation vulnerability of disease-associated gene sets, we analyzed publicly available whole genome sequencing datasets of mutations in human induced pluripotent stem cell clonal lines exposed to 12 classes of environmental carcinogens, and human lung cancers from individuals living in highly polluted regions. We compared observed rates of exposure-induced mutations in disease-related gene sets with the expected rates of mutations based on control genes randomly sampled from the genome using exact binomial tests. To explore the role of sequence characteristics in mutation vulnerability, we modeled the effects of sequence length, gene expression, and percent GC content on mutation rates of entire genes and gene coding sequences using multivariate Quasi-Poisson regressions.

Results

We demonstrate that several mutagens, including radiation and polycyclic aromatic hydrocarbons, disproportionately mutate genes related to neurodevelopmental disorders including autism spectrum disorders, schizophrenia, and attention deficit hyperactivity disorder. Other disease genes including amyotrophic lateral sclerosis, Alzheimer's disease, congenital heart disease, orofacial clefts, and coronary artery disease were generally not mutated more than expected. Longer sequence length was more strongly associated with elevated mutations in entire genes compared with mutations in coding sequences. Increased expression was associated with decreased coding sequence mutation rate, but not with the mutability of entire genes. Increased GC content was associated with increased coding sequence mutation rates but decreased mutation rates in entire genes.

Discussion

Our findings support the possibility that neurodevelopmental disorder genetic etiology is partially driven by a contribution of environment-induced germ line and somatic mutations.

Identifying Targets of Selection in Laboratory Evolution Experiments

Adaptive evolution navigates a balance between chance and determinism. The stochastic processes of mutation and drift generate phenotypic variation; however, once mutations reach an appreciable frequency in the population, their fate is governed by the deterministic action of selection, enriching for favorable genotypes and purging the less-favorable ones. The net result is that replicate populations will traverse similar-but not identical-pathways to higher fitness. This parallelism in evolutionary outcomes can be leveraged to identify the genes and pathways under selection. However, distinguishing between beneficial and neutral mutations is challenging because many beneficial mutations will be lost due to drift and clonal interference, and many neutral (and even deleterious) mutations will fix by hitchhiking. Here, we review the best practices that our laboratory uses to identify genetic targets of selection from next-generation sequencing data of evolved yeast populations. The general principles for identifying the mutations driving adaptation will apply more broadly.

Using canavanine resistance to measure mutation rates in Schizosaccharomyces pombe

We constructed a panel of S. pombe strains expressing DNA polymerase ε variants associated with cancer, specifically POLES297F, POLEV411L, POLEL424V, POLES459F, and used these to compare mutation rates determined by canavanine resistance with other selective methods. Canavanine-resistance mutation rates are broadly similar to those seen with reversion of the ade-485 mutation to adenine prototrophy, but lower than 5-fluoroorotic acid (FOA)-resistance rates (inactivation of ura4+ or ura5+ genes). Inactivation of several genes has been associated with canavanine resistance in S. pombe but surprisingly whole genome sequencing showed that 8/8 spontaneous canavanine-resistant mutants have an R175C mutation in the any1/arn1 gene. This gene encodes an α-arrestin-like protein involved in mediating Pub1 ubiquitylation of target proteins, and the phenotypic resistance to canavanine by this single mutation is similar to that shown by the original "can1-1" strain, which also has the any1R175C mutation. Some of the spontaneous mutants have additional mutations in arginine transporters, suggesting that this may marginally increase resistance to canavanine. The any1R175C strain showed internalisation of the Cat1 arginine transporter as previously reported, explaining the canavanine-resistance phenotype.

Recent insights into the evolution of mutation rates in yeast

Mutation is the origin of all genetic variation, good and bad. The mutation process can evolve in response to mutations, positive or negative selection, and genetic drift, but how these forces contribute to mutation-rate variation is an unsolved problem at the heart of genetics research. Mutations can be challenging to measure, but genome sequencing and other tools have allowed for the collection of larger and more detailed datasets, particularly in the yeast-model system. We review key hypotheses for the evolution of mutation rates and describe recent advances in understanding variation in mutational properties within and among yeast species. The multidimensional spectrum of mutations is increasingly recognized as holding valuable clues about how this important process evolves.

Targeting DNA repair pathway in cancer: Mechanisms and clinical application

Over the last decades, the growing understanding on DNA damage response (DDR) pathways has broadened the therapeutic landscape in oncology. It is becoming increasingly clear that the genomic instability of cells resulted from deficient DNA damage response contributes to the occurrence of cancer. One the other hand, these defects could also be exploited as a therapeutic opportunity, which is preferentially more deleterious in tumor cells than in normal cells. An expanding repertoire of DDR-targeting agents has rapidly expanded to inhibitors of multiple members involved in DDR pathways, including PARP, ATM, ATR, CHK1, WEE1, and DNA-PK. In this review, we sought to summarize the complex network of DNA repair machinery in cancer cells and discuss the underlying mechanism for the application of DDR inhibitors in cancer. With the past preclinical evidence and ongoing clinical trials, we also provide an overview of the history and current landscape of DDR inhibitors in cancer treatment, with special focus on the combination of DDR-targeted therapies with other cancer treatment strategies.

DNA Repair in Haploid Context

DNA repair is a well-covered topic as alteration of genetic integrity underlies many pathological conditions and important transgenerational consequences. Surprisingly, the ploidy status is rarely considered although the presence of homologous chromosomes dramatically impacts the repair capacities of cells. This is especially important for the haploid gametes as they must transfer genetic information to the offspring. An understanding of the different mechanisms monitoring genetic integrity in this context is, therefore, essential as differences in repair pathways exist that differentiate the gamete’s role in transgenerational inheritance. Hence, the oocyte must have the most reliable repair capacity while sperm, produced in large numbers and from many differentiation steps, are expected to carry de novo variations. This review describes the main DNA repair pathways with a special emphasis on ploidy. Differences between Saccharomyces cerevisiae and Schizosaccharomyces pombe are especially useful to this aim as they can maintain a diploid and haploid life cycle respectively.

DNA polymerase zeta contributes to heterochromatin replication to prevent genome instability

The DNA polymerase zeta (Polζ) plays a critical role in bypassing DNA damage. REV3L, the catalytic subunit of Polζ, is also essential in mouse embryonic development and cell proliferation for reasons that remain incompletely understood. In this study, we reveal that REV3L protein interacts with heterochromatin components including repressive histone marks and localizes in pericentromeric regions through direct interaction with HP1 dimer. We demonstrate that Polζ/REV3L ensures progression of replication forks through difficult‐to‐replicate pericentromeric heterochromatin, thereby preventing spontaneous chromosome break formation. We also find that Rev3l‐deficient cells are compromised in the repair of heterochromatin‐associated double‐stranded breaks, eliciting deletions in late‐replicating regions. Lack of REV3L leads to further consequences that may be ascribed to heterochromatin replication and repair‐associated functions of Polζ, with a disruption of the temporal replication program at specific loci. This is correlated with changes in epigenetic landscape and transcriptional control of developmentally regulated genes. These results reveal a new function of Polζ in preventing chromosome instability during replication of heterochromatic regions.

A modified fluctuation assay reveals a natural mutator phenotype that drives mutation spectrum variation within Saccharomyces cerevisiae

Although studies of Saccharomyces cerevisiae have provided many insights into mutagenesis and DNA repair, most of this work has focused on a few laboratory strains. Much less is known about the phenotypic effects of natural variation within S. cerevisiae’s DNA repair pathways. Here, we use natural polymorphisms to detect historical mutation spectrum differences among several wild and domesticated S. cerevisiae strains. To determine whether these differences are likely caused by genetic mutation rate modifiers, we use a modified fluctuation assay with a CAN1 reporter to measure de novo mutation rates and spectra in 16 of the analyzed strains. We measure a 10-fold range of mutation rates and identify two strains with distinctive mutation spectra. These strains, known as AEQ and AAR, come from the panel’s ‘Mosaic beer’ clade and share an enrichment for C > A mutations that is also observed in rare variation segregating throughout the genomes of several Mosaic beer and Mixed origin strains. Both AEQ and AAR are haploid derivatives of the diploid natural isolate CBS 1782, whose rare polymorphisms are enriched for C > A as well, suggesting that the underlying mutator allele is likely active in nature. We use a plasmid complementation test to show that AAR and AEQ share a mutator allele in the DNA repair gene OGG1, which excises 8-oxoguanine lesions that can cause C > A mutations if left unrepaired.

Chromosomal Mcm2-7 distribution and the genome replication program in species from yeast to humans

The spatio-temporal program of genome replication across eukaryotes is thought to be driven both by the uneven loading of pre-replication complexes (pre-RCs) across the genome at the onset of S-phase, and by differences in the timing of activation of these complexes during S phase. To determine the degree to which distribution of pre-RC loading alone could account for chromosomal replication patterns, we mapped the binding sites of the Mcm2-7 helicase complex (MCM) in budding yeast, fission yeast, mouse and humans. We observed similar individual MCM double-hexamer (DH) footprints across the species, but notable differences in their distribution: Footprints in budding yeast were more sharply focused compared to the other three organisms, consistent with the relative sequence specificity of replication origins in S. cerevisiae. Nonetheless, with some clear exceptions, most notably the inactive X-chromosome, much of the fluctuation in replication timing along the chromosomes in all four organisms reflected uneven chromosomal distribution of pre-replication complexes.

Gene-rich regions of the genome tend to replicate earlier in S phase than do repetitive and other non-genic regions. This may be an evolutionary consequence of the fact that replication later in S phase is associated with higher frequencies of mutation and genome rearrangement. Replication timing along the chromosome is determined by 1) events prior to S-phase that specify the locations where DNA replication can be initiated, referred to as origin licensing; and 2) the timing of activation of these licensed origins during S-phase, referred to as origin firing. To determine the relative importance of these two mechanisms, here we identify both the binding sites and the abundance of a key component of the origin licensing machinery in budding yeast, fission yeast, mice, and humans, namely the replicative helicase complex. We discovered that, with a few notable exceptions, which include the inactive X chromosome in mammals, the program of replication timing can be largely explained simply on the basis of origin licensing. Our results support a model for replication timing that emphasizes stochastic firing of origins that have been licensed before S phase begins.

Daughter-strand gaps in DNA replication - substrates of lesion processing and initiators of distress signalling

Dealing with DNA lesions during genome replication is particularly challenging because damaged replication templates interfere with the progression of the replicative DNA polymerases and thereby endanger the stability of the replisome. A variety of mechanisms for the recovery of replication forks exist, but both bacteria and eukaryotic cells also have the option of continuing replication downstream of the lesion, leaving behind a daughter-strand gap in the newly synthesized DNA. In this review, we address the significance of these single-stranded DNA structures as sites of DNA damage sensing and processing at a distance from ongoing genome replication. We describe the factors controlling the emergence of daughter-strand gaps from stalled replication intermediates, the benefits and risks of their expansion and repair via translesion synthesis or recombination-mediated template switching, and the mechanisms by which they activate local as well as global replication stress signals. Our growing understanding of daughter-strand gaps not only identifies them as targets of fundamental genome maintenance mechanisms, but also suggests that proper control over their activities has important practical implications for treatment strategies and resistance mechanisms in cancer therapy.

Divergent Evolution of Mutation Rates and Biases in the Long-Term Evolution Experiment with Escherichia coli

All organisms encode enzymes that replicate, maintain, pack, recombine, and repair their genetic material. For this reason, mutation rates and biases also evolve by mutation, variation, and natural selection. By examining metagenomic time series of the Lenski long-term evolution experiment (LTEE) with Escherichia coli (Good BH, McDonald MJ, Barrick JE, Lenski RE, Desai MM. 2017. The dynamics of molecular evolution over 60,000 generations. Nature 551(7678):45–50.), we find that local mutation rate variation has evolved during the LTEE. Each LTEE population has evolved idiosyncratic differences in their rates of point mutations, indels, and mobile element insertions, due to the fixation of various hypermutator and antimutator alleles. One LTEE population, called Ara+3, shows a strong, symmetric wave pattern in its density of point mutations, radiating from the origin of replication. This pattern is largely missing from the other LTEE populations, most of which evolved missense, indel, or structural mutations in topA, fis, and dusB—loci that all affect DNA topology. The distribution of mutations in those genes over time suggests epistasis and historical contingency in the evolution of DNA topology, which may have in turn affected local mutation rates. Overall, the replicate populations of the LTEE have largely diverged in their mutation rates and biases, even though they have adapted to identical abiotic conditions.

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.

The capacity of origins to load MCM establishes replication timing patterns

Loading of the MCM replicative helicase at origins of replication is a highly regulated process that precedes DNA replication in all eukaryotes. The stoichiometry of MCM loaded at origins has been proposed to be a key determinant of when those origins initiate replication during S phase. Nevertheless, the genome-wide regulation of MCM loading stoichiometry and its direct effect on replication timing remain unclear. In order to investigate why some origins load more MCM than others, we perturbed MCM levels in budding yeast cells and, for the first time, directly measured MCM levels and replication timing in the same experiment. Reduction of MCM levels through degradation of Mcm4, one of the six obligate components of the MCM complex, slowed progression through S phase and increased sensitivity to replication stress. Reduction of MCM levels also led to differential loading at origins during G1, revealing origins that are sensitive to reductions in MCM and others that are not. Sensitive origins loaded less MCM under normal conditions and correlated with a weak ability to recruit the origin recognition complex (ORC). Moreover, reduction of MCM loading at specific origins of replication led to a delay in their replication during S phase. In contrast, overexpression of MCM had no effects on cell cycle progression, relative MCM levels at origins, or replication timing, suggesting that, under optimal growth conditions, cellular MCM levels are not limiting for MCM loading. Our results support a model in which the loading capacity of origins is the primary determinant of MCM stoichiometry in wild-type cells, but that stoichiometry is controlled by origins' ability to recruit ORC and compete for MCM when MCM becomes limiting.

Chromatin Structure and Drug Resistance in Candida spp

Anti-microbial resistance (AMR) is currently one of the most serious threats to global human health and, appropriately, research to tackle AMR garnishes significant investment and extensive attention from the scientific community. However, most of this effort focuses on antibiotics, and research into anti-fungal resistance (AFR) is vastly under-represented in comparison. Given the growing number of vulnerable, immunocompromised individuals, as well as the positive impact global warming has on fungal growth, there is an immediate urgency to tackle fungal disease, and the disturbing rise in AFR. Chromatin structure and gene expression regulation play pivotal roles in the adaptation of fungal species to anti-fungal stress, suggesting a potential therapeutic avenue to tackle AFR. In this review we discuss both the genetic and epigenetic mechanisms by which chromatin structure can dictate AFR mechanisms and will present evidence of how pathogenic yeast, specifically from the Candida genus, modify chromatin structure to promote survival in the presence of anti-fungal drugs. We also discuss the mechanisms by which anti-chromatin therapy, specifically lysine deacetylase inhibitors, influence the acquisition and phenotypic expression of AFR in Candida spp. and their potential as effective adjuvants to mitigate against AFR.

Chromatin and Nuclear Architecture: Shaping DNA Replication in 3D

In eukaryotes, DNA replication progresses through a finely orchestrated temporal and spatial program. The 3D genome structure and nuclear architecture have recently emerged as fundamental determinants of the replication program. Factors with established roles in replication have been recognized as genome organization regulators. Exploiting paradigms from yeasts and mammals, we discuss how DNA replication is regulated in time and space through DNA-associated trans-acting factors, diffusible limiting replication initiation factors, higher-order chromatin folding, dynamic origin localization, and specific nuclear microenvironments. We present an integrated model for the regulation of DNA replication in 3D and highlight the importance of accurate spatio-temporal regulation of DNA replication in physiology and disease.

Expandable and reversible copy number amplification drives rapid adaptation to antifungal drugs

Previously, we identified long repeat sequences that are frequently associated with genome rearrangements, including copy number variation (CNV), in many diverse isolates of the human fungal pathogen Candida albicans (Todd et al., 2019). Here, we describe the rapid acquisition of novel, high copy number CNVs during adaptation to azole antifungal drugs. Single-cell karyotype analysis indicates that these CNVs appear to arise via a dicentric chromosome intermediate and breakage-fusion-bridge cycles that are repaired using multiple distinct long inverted repeat sequences. Subsequent removal of the antifungal drug can lead to a dramatic loss of the CNV and reversion to the progenitor genotype and drug susceptibility phenotype. These findings support a novel mechanism for the rapid acquisition of antifungal drug resistance and provide genomic evidence for the heterogeneity frequently observed in clinical settings.

Budding yeast complete DNA synthesis after chromosome segregation begins

To faithfully transmit genetic information, cells must replicate their entire genome before division. This is thought to be ensured by the temporal separation of replication and chromosome segregation. Here we show that in 20–40% of unperturbed yeast cells, DNA synthesis continues during anaphase, late in mitosis. High cyclin-Cdk activity inhibits DNA synthesis in metaphase, and the decrease in cyclin-Cdk activity during mitotic exit allows DNA synthesis to finish at subtelomeric and some difficult-to-replicate regions. DNA synthesis during late mitosis correlates with elevated mutation rates at subtelomeric regions, including copy number variation. Thus, yeast cells temporally overlap DNA synthesis and chromosome segregation during normal growth, possibly allowing cells to maximize population-level growth rate while simultaneously exploring greater genetic space.

In the S phase of the cell cycle, the full genome needs to be replicated before cell division occurs. Here, authors show that in budding yeast DNA synthesis is completed after chromosome segregation begins.

Checkpoint Regulation of Nuclear Tos4 Defines S Phase Arrest in Fission Yeast

From yeast to humans, the cell cycle is tightly controlled by regulatory networks that regulate cell proliferation and can be monitored by dynamic visual markers in living cells. We have observed S phase progression by monitoring nuclear accumulation of the FHA-containing DNA binding protein Tos4, which is expressed in the G1/S phase transition. We use Tos4 localization to distinguish three classes of DNA replication mutants: those that arrest with an apparent 1C DNA content and accumulate Tos4 at the restrictive temperature; those that arrest with an apparent 2C DNA content, that do not accumulate Tos4; and those that proceed into mitosis despite a 1C DNA content, again without Tos4 accumulation. Our data indicate that Tos4 localization in these conditions is responsive to checkpoint kinases, with activation of the Cds1 checkpoint kinase promoting Tos4 retention in the nucleus, and activation of the Chk1 damage checkpoint promoting its turnover. Tos4 localization therefore allows us to monitor checkpoint-dependent activation that responds to replication failure in early vs. late S phase.

Spatiotemporal regulation of PCNA ubiquitination in damage tolerance pathways

DNA is constantly exposed to a wide variety of exogenous and endogenous agents, and most DNA lesions inhibit DNA synthesis. To cope with such problems during replication, cells have molecular mechanisms to resume DNA synthesis in the presence of DNA lesions, which are known as DNA damage tolerance (DDT) pathways. The concept of ubiquitination-mediated regulation of DDT pathways in eukaryotes was established via genetic studies in the yeast Saccharomyces cerevisiae, in which two branches of the DDT pathway are regulated via ubiquitination of proliferating cell nuclear antigen (PCNA): translesion DNA synthesis (TLS) and homology-dependent repair (HDR), which are stimulated by mono- and polyubiquitination of PCNA, respectively. Over the subsequent nearly two decades, significant progress has been made in understanding the mechanisms that regulate DDT pathways in other eukaryotes. Importantly, TLS is intrinsically error-prone because of the miscoding nature of most damaged nucleotides and inaccurate replication of undamaged templates by TLS polymerases (pols), whereas HDR is theoretically error-free because the DNA synthesis is thought to be predominantly performed by pol δ, an accurate replicative DNA pol, using the undamaged sister chromatid as its template. Thus, the regulation of the choice between the TLS and HDR pathways is critical to determine the appropriate biological outcomes caused by DNA damage. In this review, we summarize our current understanding of the species-specific regulatory mechanisms of PCNA ubiquitination and how cells choose between TLS and HDR. We then provide a hypothetical model for the spatiotemporal regulation of DDT pathways in human cells.

Recurrent somatic mutations reveal new insights into consequences of mutagenic processes in cancer

The sheer size of the human genome makes it improbable that identical somatic mutations at the exact same position are observed in multiple tumours solely by chance. The scarcity of cancer driver mutations also precludes positive selection as the sole explanation. Therefore, recurrent mutations may be highly informative of characteristics of mutational processes. To explore the potential, we use recurrence as a starting point to cluster >2,500 whole genomes of a pan-cancer cohort. We describe each genome with 13 recurrence-based and 29 general mutational features. Using principal component analysis we reduce the dimensionality and create independent features. We apply hierarchical clustering to the first 18 principal components followed by k-means clustering. We show that the resulting 16 clusters capture clinically relevant cancer phenotypes. High levels of recurrent substitutions separate the clusters that we link to UV-light exposure and deregulated activity of POLE from the one representing defective mismatch repair, which shows high levels of recurrent insertions/deletions. Recurrence of both mutation types characterizes cancer genomes with somatic hypermutation of immunoglobulin genes and the cluster of genomes exposed to gastric acid. Low levels of recurrence are observed for the cluster where tobacco-smoke exposure induces mutagenesis and the one linked to increased activity of cytidine deaminases. Notably, the majority of substitutions are recurrent in a single tumour type, while recurrent insertions/deletions point to shared processes between tumour types. Recurrence also reveals susceptible sequence motifs, including TT[C>A]TTT and AAC[T>G]T for the POLE and 'gastric-acid exposure' clusters, respectively. Moreover, we refine knowledge of mutagenesis, including increased C/G deletion levels in general for lung tumours and specifically in midsize homopolymer sequence contexts for microsatellite instable tumours. Our findings are an important step towards the development of a generic cancer diagnostic test for clinical practice based on whole-genome sequencing that could replace multiple diagnostics currently in use.

Rif1 Binding and Control of Chromosome-Internal DNA Replication Origins Is Limited by Telomere Sequestration

The Saccharomyces cerevisiae telomere-binding protein Rif1 plays an evolutionarily conserved role in control of DNA replication timing by promoting PP1-dependent dephosphorylation of replication initiation factors. However, ScRif1 binding outside of telomeres has never been detected, and it has thus been unclear whether Rif1 acts directly on the replication origins that it controls. Here, we show that, in unperturbed yeast cells, Rif1 primarily regulates late-replicating origins within 100 kb of a telomere. Using the chromatin endogenous cleavage ChEC-seq technique, we robustly detect Rif1 at late-replicating origins that we show are targets of its inhibitory action. Interestingly, abrogation of Rif1 telomere association by mutation of its Rap1-binding module increases Rif1 binding and origin inhibition elsewhere in the genome. Our results indicate that Rif1 inhibits replication initiation by interacting directly with origins and suggest that Rap1-dependent sequestration of Rif1 increases its effective concentration near telomeres, while limiting its action at chromosome-internal sites.

FANCD2 binding identifies conserved fragile sites at large transcribed genes in avian cells

Common Chromosomal Fragile Sites (CFSs) are specific genomic regions prone to form breaks on metaphase chromosomes in response to replication stress. Moreover, CFSs are mutational hotspots in cancer genomes, showing that the mutational mechanisms that operate at CFSs are highly active in cancer cells. Orthologs of human CFSs are found in a number of other mammals, but the extent of CFS conservation beyond the mammalian lineage is unclear. Characterization of CFSs from distantly related organisms can provide new insight into the biology underlying CFSs. Here, we have mapped CFSs in an avian cell line. We find that, overall the most significant CFSs coincide with extremely large conserved genes, from which very long transcripts are produced. However, no significant correlation between any sequence characteristics and CFSs is found. Moreover, we identified putative early replicating fragile sites (ERFSs), which is a distinct class of fragile sites and we developed a fluctuation analysis revealing high mutation rates at the CFS gene PARK2, with deletions as the most prevalent mutation. Finally, we show that avian homologs of the human CFS genes despite their fragility have resisted the general intron size reduction observed in birds suggesting that CFSs have a conserved biological function.

Post-replication repair: Rad5/HLTF regulation, activity on undamaged templates, and relationship to cancer

The eukaryotic post-replication repair (PRR) pathway allows completion of DNA replication when replication forks encounter lesions on the DNA template and are mediated by post-translational ubiquitination of the DNA sliding clamp proliferating cell nuclear antigen (PCNA). Monoubiquitinated PCNA recruits translesion synthesis (TLS) polymerases to replicate past DNA lesions in an error-prone manner while addition of K63-linked polyubiquitin chains signals for error-free template switching to the sister chromatid. Central to both branches is the E3 ubiquitin ligase and DNA helicase Rad5/helicase-like transcription factor (HLTF). Mutations in PRR pathway components lead to genomic rearrangements, cancer predisposition, and cancer progression. Recent studies have challenged the notion that the PRR pathway is involved only in DNA lesion tolerance and have shed new light on its roles in cancer progression. Molecular details of Rad5/HLTF recruitment and function at replication forks have emerged. Mounting evidence indicates that PRR is required during lesion-less replication stress, leading to TLS polymerase activity on undamaged templates. Analysis of PRR mutation status in human cancers and PRR function in cancer models indicates that down regulation of PRR activity is a viable strategy to inhibit cancer cell growth and reduce chemoresistance. Here, we review these findings, discuss how they change our views of current PRR models, and look forward to targeting the PRR pathway in the clinic.

Germline Structural Variations Are Preferential Sites of DNA Replication Timing Plasticity during Development

The DNA replication timing program is modulated throughout development and is also one of the main factors influencing the distribution of mutation rates across the genome. However, the relationship between the mutagenic influence of replication timing and its developmental plasticity remains unexplored. Here, we studied the distribution of copy number variations (CNVs) and single nucleotide polymorphisms across the zebrafish genome in relation to changes in DNA replication timing during embryonic development in this model vertebrate species. We show that CNV sites exhibit strong replication timing plasticity during development, replicating significantly early during early development but significantly late during more advanced developmental stages. Reciprocally, genomic regions that changed their replication timing during development contained a higher proportion of CNVs than developmentally constant regions. Developmentally plastic CNV sites, in particular those that become delayed in their replication timing, were enriched for the clustered protocadherins, a set of genes important for neuronal development that have undergone extensive genetic and epigenetic diversification during zebrafish evolution. In contrast, single nucleotide polymorphism sites replicated consistently early throughout embryonic development, highlighting a unique aspect of the zebrafish genome. Our results uncover a hitherto unrecognized interface between development and evolution.

Sir2 suppresses transcription-mediated displacement of Mcm2-7 replicative helicases at the ribosomal DNA repeats

Repetitive DNA sequences within eukaryotic heterochromatin are poorly transcribed and replicate late in S-phase. In Saccharomyces cerevisiae, the histone deacetylase Sir2 is required for both transcriptional silencing and late replication at the repetitive ribosomal DNA arrays (rDNA). Despite the widespread association between transcription and replication timing, it remains unclear how transcription might impinge on replication, or vice versa. Here we show that, when silencing of an RNA polymerase II (RNA Pol II)-transcribed non-coding RNA at the rDNA is disrupted by SIR2 deletion, RNA polymerase pushes and thereby relocalizes replicative Mcm2-7 helicases away from their loading sites to an adjacent region with low nucleosome occupancy, and this relocalization is associated with increased rDNA origin efficiency. Our results suggest a model in which two of the major defining features of heterochromatin, transcriptional silencing and late replication, are mechanistically linked through suppression of polymerase-mediated displacement of replication initiation complexes.

Eukaryotic genomes typically contain large regions of repetitive DNA, referred to as heterochromatin, that are both transcriptionally silent and late replicating. We provide a possible explanation for the association between transcriptional silencing and late replication. Budding yeast contains a histone deacetylase called SIR2 that was originally identified as a transcriptional repressor, but was later also found to ensure late replication of repetitive ribosomal DNA (rDNA) sequences. We show that the transcription that occurs in the absence of SIR2 directly displaces the helicase required for replication initiation at the rDNA. This work represents an important advance in understanding the interplay between transcription and replication at repetitive sequences by directly linking transcription with replication machinery in heterochromatin.

Insights into the Link between the Organization of DNA Replication and the Mutational Landscape

The generation of a complete and accurate copy of the genetic material during each cell cycle is integral to cell growth and proliferation. However, genetic diversity is essential for adaptation and evolution, and the process of DNA replication is a fundamental source of mutations. Genome alterations do not accumulate randomly, with variations in the types and frequencies of mutations that arise in different genomic regions. Intriguingly, recent studies revealed a striking link between the mutational landscape of a genome and the spatial and temporal organization of DNA replication, referred to as the replication program. In our review, we discuss how this program may contribute to shaping the profile and spectrum of genetic alterations, with implications for genome dynamics and organismal evolution in natural and pathological contexts.

Epigenetic gene silencing alters the mechanisms and rate of evolutionary adaptation

Epigenetic, non-DNA sequence-based inheritance can potentially contribute to adaptation but, due to its transient nature and the difficulty involved in uncoupling it from genetic variation, it is unclear whether it has any effect on long-term evolution. However, short-term epigenetic inheritance may interact with genetic change by modifying the rate and type of adaptive mutations. Here, we test this notion in an experimental evolution set-up in yeast. We tune low, intermediate and high levels of heritable silencing of a URA3 reporter under selection by insertion at different positions within silent subtelomeric chromatin in otherwise isogenic Saccharomyces cerevisiae. Heritable silencing does not impact mutation rate but drives population size expansion and rapid epigenetic adaptation. This eventually leads to genetic assimilation of the silent phenotype by mutations that reduce or abolish URA3 expression. Moreover, at intermediate or low levels of heritable silencing we find that populations evolve more rapidly by accumulation of adaptive mutations, in part through acquisition of novel alleles that enhance gene silencing, aiding accelerated adaptation. We provide an experimental proof of concept that defines the impact and mechanisms of how short-term epigenetic inheritance can shape adaptive evolution.

Dynamics of DNA replication in a eukaryotic cell

Each genomic locus in a eukaryotic cell has a distinct average time of replication during S phase that depends on the spatial and temporal pattern of replication initiation events. Replication timing can affect genomic integrity because late replication is associated with an increased mutation rate. For most eukaryotes, the features of the genome that specify the location and timing of initiation events are unknown. To investigate these features for the fission yeast, Schizosaccharomyces pombe, we developed an integrative model to analyze large single-molecule and global genomic datasets. The model provides an accurate description of the complex dynamics of S. pombe DNA replication at high resolution. We present evidence that there are many more potential initiation sites in the S. pombe genome than previously identified and that the distribution of these sites is primarily determined by two factors: the sequence preferences of the origin recognition complex (ORC), and the interference of transcription with the assembly or stability of prereplication complexes (pre-RCs). We suggest that in addition to directly interfering with initiation, transcription has driven the evolution of the binding properties of ORC in S. pombe and other eukaryotic species to target pre-RC assembly to regions of the genome that are less likely to be transcribed.

Guidelines for DNA recombination and repair studies: Cellular assays of DNA repair pathways

Understanding the plasticity of genomes has been greatly aided by assays for recombination, repair and mutagenesis. These assays have been developed in microbial systems that provide the advantages of genetic and molecular reporters that can readily be manipulated. Cellular assays comprise genetic, molecular, and cytological reporters. The assays are powerful tools but each comes with its particular advantages and limitations. Here the most commonly used assays are reviewed, discussed, and presented as the guidelines for future studies.

Linking the organization of DNA replication with genome maintenance

The spatial and temporal organization of genome duplication, also referred to as the replication program, is defined by the distribution and the activities of the sites of replication initiation across the genome. Alterations to the replication profile are associated with cell fate changes during development and in pathologies, but the importance of undergoing S phase with distinct and specific programs remains largely unexplored. We have recently addressed this question, focusing on the interplay between the replication program and genome maintenance. In particular, we demonstrated that when cells encounter challenges to DNA synthesis, the organization of DNA replication drives the response to replication stress that is mediated by the ATR/Rad3 checkpoint pathway, thus shaping the pattern of genome instability along the chromosomes. In this review, we present the major findings of our study and discuss how they may bring new perspectives to our understanding of the biological importance of the replication program.

Superovulation Influences Methylation Reprogramming and Delays Onset of DNA Replication in Both Pronuclei of Mouse Zygotes

Although an essential component of assisted reproductive technologies, ovarian stimulation, or superovulation, may interfere with the epigenetic reprogramming machinery during early embryogenesis and gametogenesis. To investigate the possible impact of superovulation particularly on the methylation reprogramming process directly after fertilization, we performed immunofluorescence staining of pronuclear (PN) stage embryos with antibodies against 5mC and 5hmC. PN stage embryos obtained by superovulation displayed an increased incidence of abnormal methylation and hydroxymethylation patterns in both maternal and paternal pronuclear DNA. Subsequent single-cell RT-qPCR analyses of the Tet1, Tet2, and Tet3 genes revealed no significant expression differences between PN stage embryos from spontaneously and superovulated matings that could be causative for the abnormal methylation and hydroxymethylation patterns. To analyze the possible contribution of TET-independent replication-associated demethylation mechanisms, we then determined the 5mC and 5hmC levels of PN stage mouse embryos using immunofluorescence analyses after inhibition of DNA replication with aphidicolin. Inhibition of DNA replication had no effect on abnormal methylation and hydroxymethylation patterns that still persisted in the superovulated group. Interestingly, the onset of DNA replication, which was also analyzed in these experiments, was remarkably delayed in the superovulated group. Our findings imply an impact of superovulation on both replication-dependent and -independent or yet unknown demethylation mechanisms in PN stage mouse embryos. In addition, they reveal for the first time a negative effect of superovulation on the initiation of DNA replication in PN stage mouse embryos.

Identifying Drivers of Parallel Evolution: A Regression Model Approach

Parallel evolution, defined as identical changes arising in independent populations, is often attributed to similar selective pressures favoring the fixation of identical genetic changes. However, some level of parallel evolution is also expected if mutation rates are heterogeneous across regions of the genome. Theory suggests that mutation and selection can have equal impacts on patterns of parallel evolution; however, empirical studies have yet to jointly quantify the importance of these two processes. Here, we introduce several statistical models to examine the contributions of mutation and selection heterogeneity to shaping parallel evolutionary changes at the gene-level. Using this framework, we analyze published data from forty experimentally evolved Saccharomyces cerevisiae populations. We can partition the effects of a number of genomic variables into those affecting patterns of parallel evolution via effects on the rate of arising mutations, and those affecting the retention versus loss of the arising mutations (i.e., selection). Our results suggest that gene-to-gene heterogeneity in both mutation and selection, associated with gene length, recombination rate, and number of protein domains drive parallel evolution at both synonymous and nonsynonymous sites. While there are still a number of parallel changes that are not well described, we show that allowing for heterogeneous rates of mutation and selection can provide improved predictions of the prevalence and degree of parallel evolution.

Reduced intrinsic DNA curvature leads to increased mutation rate

BACKGROUND

Mutation rates vary across the genome. Many trans factors that influence mutation rates have been identified, as have specific sequence motifs at the 1-7-bp scale, but cis elements remain poorly characterized. The lack of understanding regarding why different sequences have different mutation rates hampers our ability to identify positive selection in evolution and to identify driver mutations in tumorigenesis.

RESULTS

Here, we use a combination of synthetic genes and sequences of thousands of isolated yeast colonies to show that intrinsic DNA curvature is a major cis determinant of mutation rate. Mutation rate negatively correlates with DNA curvature within genes, and a 10% decrease in curvature results in a 70% increase in mutation rate. Consistently, both yeast and humans accumulate mutations in regions with small curvature. We further show that this effect is due to differences in the intrinsic mutation rate, likely due to differences in mutagen sensitivity and not due to differences in the local activity of DNA repair.

CONCLUSIONS

Our study establishes a framework for understanding the cis properties of DNA sequence in modulating the local mutation rate and identifies a novel causal source of non-uniform mutation rates across the genome.

Mutational signature distribution varies with DNA replication timing and strand asymmetry

BACKGROUND

DNA replication plays an important role in mutagenesis, yet little is known about how it interacts with other mutagenic processes. Here, we use somatic mutation signatures-each representing a mutagenic process-derived from 3056 patients spanning 19 cancer types to quantify the strand asymmetry of mutational signatures around replication origins and between early and late replicating regions.

RESULTS

We observe that most of the detected mutational signatures are significantly correlated with the timing or direction of DNA replication. The properties of these associations are distinct for different signatures and shed new light on several mutagenic processes. For example, our results suggest that oxidative damage to the nucleotide pool substantially contributes to the mutational landscape of esophageal adenocarcinoma.

CONCLUSIONS

Together, our results indicate an interaction between DNA replication, the associated damage repair, and most mutagenic processes.

Periodic Variation of Mutation Rates in Bacterial Genomes Associated with Replication Timing

The causes and consequences of spatiotemporal variation in mutation rates remain to be explored in nearly all organisms. Here we examine relationships between local mutation rates and replication timing in three bacterial species whose genomes have multiple chromosomes: Vibrio fischeri, Vibrio cholerae, and Burkholderia cenocepacia. Following five mutation accumulation experiments with these bacteria conducted in the near absence of natural selection, the genomes of clones from each lineage were sequenced and analyzed to identify variation in mutation rates and spectra. In lineages lacking mismatch repair, base substitution mutation rates vary in a mirrored wave-like pattern on opposing replichores of the large chromosomes of V. fischeri and V. cholerae, where concurrently replicated regions experience similar base substitution mutation rates. The base substitution mutation rates on the small chromosome are less variable in both species but occur at similar rates to those in the concurrently replicated regions of the large chromosome. Neither nucleotide composition nor frequency of nucleotide motifs differed among regions experiencing high and low base substitution rates, which along with the inferred ~800-kb wave period suggests that the source of the periodicity is not sequence specific but rather a systematic process related to the cell cycle. These results support the notion that base substitution mutation rates are likely to vary systematically across many bacterial genomes, which exposes certain genes to elevated deleterious mutational load.

That mutation rates vary within bacterial genomes is well known, but the detailed study of these biases has been made possible only recently with contemporary sequencing methods. We applied these methods to understand how bacterial genomes with multiple chromosomes, like those of Vibrio and Burkholderia, might experience heterogeneous mutation rates because of their unusual replication and the greater genetic diversity found on smaller chromosomes. This study captured thousands of mutations and revealed wave-like rate variation that is synchronized with replication timing and not explained by sequence context. The scale of this rate variation over hundreds of kilobases of DNA strongly suggests that a temporally regulated cellular process may generate wave-like variation in mutation risk. These findings add to our understanding of how mutation risk is distributed across bacterial and likely also eukaryotic genomes, owing to their highly conserved replication and repair machinery.

GC content elevates mutation and recombination rates in the yeast Saccharomyces cerevisiae

The chromosomes of many eukaryotes have regions of high GC content interspersed with regions of low GC content. In the yeast Saccharomyces cerevisiae, high-GC regions are often associated with high levels of meiotic recombination. In this study, we constructed URA3 genes that differ substantially in their base composition [URA3-AT (31% GC), URA3-WT (43% GC), and URA3-GC (63% GC)] but encode proteins with the same amino acid sequence. The strain with URA3-GC had an approximately sevenfold elevated rate of ura3 mutations compared with the strains with URA3-WT or URA3-AT About half of these mutations were single-base substitutions and were dependent on the error-prone DNA polymerase ζ. About 30% were deletions or duplications between short (5-10 base) direct repeats resulting from DNA polymerase slippage. The URA3-GC gene also had elevated rates of meiotic and mitotic recombination relative to the URA3-AT or URA3-WT genes. Thus, base composition has a substantial effect on the basic parameters of genome stability and evolution.

The organization of genome duplication is a critical determinant of the landscape of genome maintenance

Genome duplication is essential for cell proliferation, and the mechanisms regulating its execution are highly conserved. These processes give rise to a spatiotemporal organization of replication initiation across the genome, referred to as the replication program. Despite the identification of such programs in diverse eukaryotic organisms, their biological importance for cellular physiology remains largely unexplored. We address this fundamental question in the context of genome maintenance, taking advantage of the inappropriate origin firing that occurs when fission yeast cells lacking the Rad3/ATR checkpoint kinase are subjected to replication stress. Using this model, we demonstrate that the replication program quantitatively dictates the extent of origin de-regulation and the clustered localization of these events. Furthermore, our results uncover an accumulation of abnormal levels of single-stranded DNA (ssDNA) and the Rad52 repair protein at de-regulated origins. We show that these loci constitute a defining source of the overall ssDNA and Rad52 hotspots in the genome, generating a signature pattern of instability along the chromosomes. We then induce a genome-wide reprogramming of origin usage and evaluate its consequences in our experimental system. This leads to a complete redistribution of the sites of both inappropriate initiation and associated Rad52 recruitment. We therefore conclude that the organization of genome duplication governs the checkpoint control of origin-associated hotspots of instability and plays an integral role in shaping the landscape of genome maintenance.

The genome-wide rate and spectrum of spontaneous mutations differ between haploid and diploid yeast

By altering the dynamics of DNA replication and repair, alternative ploidy states may experience different rates and types of new mutations, leading to divergent evolutionary outcomes. We report a direct comparison of the genome-wide spectrum of spontaneous mutations arising in haploids and diploids following a mutation-accumulation experiment in the budding yeast Saccharomyces cerevisiae Characterizing the number, types, locations, and effects of thousands of mutations revealed that haploids were more prone to single-nucleotide mutations (SNMs) and mitochondrial mutations, while larger structural changes were more common in diploids. Mutations were more likely to be detrimental in diploids, even after accounting for the large impact of structural changes, contrary to the prediction that mutations would have weaker effects, due to masking, in diploids. Haploidy is expected to reduce the opportunity for conservative DNA repair involving homologous chromosomes, increasing the insertion-deletion rate, but we found little support for this idea. Instead, haploids were more susceptible to SNMs in late-replicating genomic regions, resulting in a ploidy difference in the spectrum of substitutions. In diploids, we detect mutation rate variation among chromosomes in association with centromere location, a finding that is supported by published polymorphism data. Diploids are not simply doubled haploids; instead, our results predict that the spectrum of spontaneous mutations will substantially shape the dynamics of genome evolution in haploid and diploid populations.

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.