The tRNA-Tyr gene family of Saccharomyces cerevisiae: agents of phenotypic variation and position effects on mutation frequency

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.

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.

Muver, a computational framework for accurately calling accumulated mutations

BACKGROUND

Identification of mutations from next-generation sequencing data typically requires a balance between sensitivity and accuracy. This is particularly true of DNA insertions and deletions (indels), that can impart significant phenotypic consequences on cells but are harder to call than substitution mutations from whole genome mutation accumulation experiments. To overcome these difficulties, we present muver, a computational framework that integrates established bioinformatics tools with novel analytical methods to generate mutation calls with the extremely low false positive rates and high sensitivity required for accurate mutation rate determination and comparison.

RESULTS

Muver uses statistical comparison of ancestral and descendant allelic frequencies to identify variant loci and assigns genotypes with models that include per-sample assessments of sequencing errors by mutation type and repeat context. Muver identifies maximally parsimonious mutation pathways that connect these genotypes, differentiating potential allelic conversion events and delineating ambiguities in mutation location, type, and size. Benchmarking with a human gold standard father-son pair demonstrates muver's sensitivity and low false positive rates. In DNA mismatch repair (MMR) deficient Saccharomyces cerevisiae, muver detects multi-base deletions in homopolymers longer than the replicative polymerase footprint at rates greater than predicted for sequential single-base deletions, implying a novel multi-repeat-unit slippage mechanism.

CONCLUSIONS

Benchmarking results demonstrate the high accuracy and sensitivity achieved with muver, particularly for indels, relative to available tools. Applied to an MMR-deficient Saccharomyces cerevisiae system, muver mutation calls facilitate mechanistic insights into DNA replication fidelity.

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.

Cancer-associated isocitrate dehydrogenase mutations induce mitochondrial DNA instability

A major advance in understanding the progression and prognostic outcome of certain cancers, such as low-grade gliomas, acute myeloid leukaemia, and chondrosarcomas, has been the identification of early-occurring mutations in the NADP⁺-dependent isocitrate dehydrogenase genes IDH1 and IDH2 These mutations result in the production of the onco-metabolite D-2-hydroxyglutarate (2HG), thought to contribute to disease progression. To better understand the mechanisms of 2HG pathophysiology, we introduced the analogous glioma-associated mutations into the NADP⁺isocitrate dehydrogenase genes (IDP1, IDP2, IDP3) in Saccharomyces cerevisiae Intriguingly, expression of the mitochondrial IDP1^R148H mutant allele results in high levels of 2HG production as well as extensive mtDNA loss and respiration defects. We find no evidence for a reactive oxygen-mediated mechanism mediating this mtDNA loss. Instead, we show that 2HG production perturbs the iron sensing mechanisms as indicated by upregulation of the Aft1-controlled iron regulon and a concomitant increase in iron levels. Accordingly, iron chelation, or overexpression of a truncated AFT1 allele that dampens transcription of the iron regulon, suppresses the loss of respirative capacity. Additional suppressing factors include overexpression of the mitochondrial aldehyde dehydrogenase gene ALD5 or disruption of the retrograde response transcription factor RTG1 Furthermore, elevated α-ketoglutarate levels also suppress 2HG-mediated respiration loss; consistent with a mechanism by which 2HG contributes to mtDNA loss by acting as a toxic α-ketoglutarate analog. Our findings provide insight into the mechanisms that may contribute to 2HG oncogenicity in glioma and acute myeloid leukaemia progression, with the promise for innovative diagnostic and prognostic strategies and novel therapeutic modalities.

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.

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

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

Cytocidal amino acid starvation of Saccharomyces cerevisiae and Candida albicans acetolactate synthase (ilv2{Delta}) mutants is influenced by the carbon source and rapamycin

The isoleucine and valine biosynthetic enzyme acetolactate synthase (Ilv2p) is an attractive antifungal drug target, since the isoleucine and valine biosynthetic pathway is not present in mammals, Saccharomyces cerevisiae ilv2Delta mutants do not survive in vivo, Cryptococcus neoformans ilv2 mutants are avirulent, and both S. cerevisiae and Cr. neoformans ilv2 mutants die upon isoleucine and valine starvation. To further explore the potential of Ilv2p as an antifungal drug target, we disrupted Candida albicans ILV2, and demonstrated that Ca. albicans ilv2Delta mutants were significantly attenuated in virulence, and were also profoundly starvation-cidal, with a greater than 100-fold reduction in viability after only 4 h of isoleucine and valine starvation. As fungicidal starvation would be advantageous for drug design, we explored the basis of the starvation-cidal phenotype in both S. cerevisiae and Ca. albicans ilv2Delta mutants. Since the mutation of ILV1, required for the first step of isoleucine biosynthesis, did not suppress the ilv2Delta starvation-cidal defects in either species, the cidal phenotype was not due to alpha-ketobutyrate accumulation. We found that starvation for isoleucine alone was more deleterious in Ca. albicans than in S. cerevisiae, and starvation for valine was more deleterious than for isoleucine in both species. Interestingly, while the target of rapamycin (TOR) pathway inhibitor rapamycin further reduced S. cerevisiae ilv2Delta starvation viability, it increased Ca. albicans ilv1Delta and ilv2Delta viability. Furthermore, the recovery from starvation was dependent on the carbon source present during recovery for S. cerevisiae ilv2Delta mutants, reminiscent of isoleucine and valine starvation inducing a viable but non-culturable-like state in this species, while Ca. albicans ilv1Delta and ilv2 Delta viability was influenced by the carbon source present during starvation, supporting a role for glucose wasting in the Ca. albicans cidal phenotype.

Sequential elimination of major-effect contributors identifies additional quantitative trait loci conditioning high-temperature growth in yeast

Several quantitative trait loci (QTL) mapping strategies can successfully identify major-effect loci, but often have poor success detecting loci with minor effects, potentially due to the confounding effects of major loci, epistasis, and limited sample sizes. To overcome such difficulties, we used a targeted backcross mapping strategy that genetically eliminated the effect of a previously identified major QTL underlying high-temperature growth (Htg) in yeast. This strategy facilitated the mapping of three novel QTL contributing to Htg of a clinically derived yeast strain. One QTL, which is linked to the previously identified major-effect QTL, was dissected, and NCS2 was identified as the causative gene. The interaction of the NCS2 QTL with the first major-effect QTL was background dependent, revealing a complex QTL architecture spanning these two linked loci. Such complex architecture suggests that more genes than can be predicted are likely to contribute to quantitative traits. The targeted backcrossing approach overcomes the difficulties posed by sample size, genetic linkage, and epistatic effects and facilitates identification of additional alleles with smaller contributions to complex traits.

Estimating the per-base-pair mutation rate in the yeast Saccharomyces cerevisiae

Although mutation rates are a key determinant of the rate of evolution they are difficult to measure precisely and global mutations rates (mutations per genome per generation) are often extrapolated from the per-base-pair mutation rate assuming that mutation rate is uniform across the genome. Using budding yeast, we describe an improved method for the accurate calculation of mutation rates based on the fluctuation assay. Our analysis suggests that the per-base-pair mutation rates at two genes differ significantly (3.80x10(-10) at URA3 and 6.44x10(-10) at CAN1) and we propose a definition for the effective target size of genes (the probability that a mutation inactivates the gene) that acknowledges that the mutation rate is nonuniform across the genome.

Variation in efficiency of DNA mismatch repair at different sites in the yeast genome

Evolutionary studies have suggested that mutation rates vary significantly at different positions in the eukaryotic genome. The mechanism that is responsible for this context-dependence of mutation rates is not understood. We demonstrate experimentally that frameshift mutation rates in yeast microsatellites depend on the genomic context and that this variation primarily reflects the context-dependence of the efficiency of DNA mismatch repair. We measured the stability of a 16.5-repeat polyGT tract by using a reporter gene (URA3-GT) in which the microsatellite was inserted in-frame into the yeast URA3 gene. We constructed 10 isogenic yeast strains with the reporter gene at different locations in the genome. Rates of frameshift mutations that abolished the correct reading frame of this gene were determined by fluctuation analysis. A 16-fold difference was found among these strains. We made mismatch-repair-deficient (msh2) derivatives of six of the strains. Mutation rates were elevated for all of these strains, but the differences in rates among the strains were substantially reduced. The simplest interpretation of this result is that the efficiency of DNA mismatch repair varies in different regions of the genome, perhaps reflecting some aspect of chromosome structure.

Nonsense mutations in the essential gene SUP35 of Saccharomyces cerevisiae are non-lethal

In the present work we have characterized for the first time non-lethal nonsense mutations in the essential gene SUP35, which codes for the translation termination factor eRF3 in Saccharomyces cerevisiae. The screen used was based on selection for simultaneous suppression of two auxotrophic nonsense mutations. Among 48 mutants obtained, sixteen were distinguished by the production of a reduced amount of eRF3, suggesting the appearance of nonsense mutations. Fifteen of the total mutants were sequenced, and the presence of nonsense mutations was confirmed for nine of them. Thus a substantial fraction of the sup35 mutations recovered are nonsense mutations located in different regions of SUP35, and such mutants are easily identified by the fact that they express reduced amounts of eRF3. Nonsense mutations in the SUP35 gene do not lead to a decrease in levels of SUP35 mRNA and do not influence the steady-state level of eRF1. The ability of these mutations to complement SUP35 gene disruption mutations in different genetic backgrounds and in the absence of any tRNA suppressor mutation was demonstrated. The missense mutations studied, unlike nonsense mutations, do not decrease steady-state amounts of eRF3.