Correction: Repair of multiple simultaneous double-strand breaks causes bursts of genome-wide clustered hypermutation.

[This corrects the article DOI: 10.1371/journal.pbio.3000464.].

Characterization of systemic genomic instability in budding yeast.

Conventional models of genome evolution are centered around the principle that mutations form independently of each other and build up slowly over time. We characterized the occurrence of bursts of genome-wide loss-of-heterozygosity (LOH) in Saccharomyces cerevisiae, providing support for an additional nonindependent and faster mode of mutation accumulation. We initially characterized a yeast clone isolated for carrying an LOH event at a specific chromosome site, and surprisingly found that it also carried multiple unselected rearrangements elsewhere in its genome. Whole-genome analysis of over 100 additional clones selected for carrying primary LOH tracts revealed that they too contained unselected structural alterations more often than control clones obtained without any selection. We also measured the rates of coincident LOH at two different chromosomes and found that double LOH formed at rates 14- to 150-fold higher than expected if the two underlying single LOH events occurred independently of each other. These results were consistent across different strain backgrounds and in mutants incapable of entering meiosis. Our results indicate that a subset of mitotic cells within a population can experience discrete episodes of systemic genomic instability, when the entire genome becomes vulnerable and multiple chromosomal alterations can form over a narrow time window. They are reminiscent of early reports from the classic yeast genetics literature, as well as recent studies in humans, both in cancer and genomic disorder contexts. The experimental model we describe provides a system to further dissect the fundamental biological processes responsible for punctuated bursts of structural genomic variation.

Long transposon-rich centromeres in an oomycete reveal divergence of centromere features in Stramenopila-Alveolata-Rhizaria lineages.

Centromeres are chromosomal regions that serve as platforms for kinetochore assembly and spindle attachments, ensuring accurate chromosome segregation during cell division. Despite functional conservation, centromere DNA sequences are diverse and often repetitive, making them challenging to assemble and identify. Here, we describe centromeres in an oomycete Phytophthora sojae by combining long-read sequencing-based genome assembly and chromatin immunoprecipitation for the centromeric histone CENP-A followed by high-throughput sequencing (ChIP-seq). P. sojae centromeres cluster at a single focus at different life stages and during nuclear division. We report an improved genome assembly of the P. sojae reference strain, which enabled identification of 15 enriched CENP-A binding regions as putative centromeres. By focusing on a subset of these regions, we demonstrate that centromeres in P. sojae are regional, spanning 211 to 356 kb. Most of these regions are transposon-rich, poorly transcribed, and lack the histone modification H3K4me2 but are embedded within regions with the heterochromatin marks H3K9me3 and H3K27me3. Strikingly, we discovered a Copia-like transposon (CoLT) that is highly enriched in the CENP-A chromatin. Similar clustered elements are also found in oomycete relatives of P. sojae, and may be applied as a criterion for prediction of oomycete centromeres. This work reveals a divergence of centromere features in oomycetes as compared to other organisms in the Stramenopila-Alveolata-Rhizaria (SAR) supergroup including diatoms and Plasmodium falciparum that have relatively short and simple regional centromeres. Identification of P. sojae centromeres in turn also advances the genome assembly.

Repair of multiple simultaneous double-strand breaks causes bursts of genome-wide clustered hypermutation.

A single cancer genome can harbor thousands of clustered mutations. Mutation signature analyses have revealed that the origin of clusters are lesions in long tracts of single-stranded (ss) DNA damaged by apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like (APOBEC) cytidine deaminases, raising questions about molecular mechanisms that generate long ssDNA vulnerable to hypermutation. Here, we show that ssDNA intermediates formed during the repair of gamma-induced bursts of double-strand breaks (DSBs) in the presence of APOBEC3A in yeast lead to multiple APOBEC-induced clusters similar to cancer. We identified three independent pathways enabling cluster formation associated with repairing bursts of DSBs: 5' to 3' bidirectional resection, unidirectional resection, and break-induced replication (BIR). Analysis of millions of mutations in APOBEC-hypermutated cancer genomes revealed that cancer tolerance to formation of hypermutable ssDNA is similar to yeast and that the predominant pattern of clustered mutagenesis is the same as in resection-defective yeast, suggesting that cluster formation in cancers is driven by a BIR-like mechanism. The phenomenon of genome-wide burst of clustered mutagenesis revealed by our study can play an important role in generating somatic hypermutation in cancers as well as in noncancerous cells.

Mutation signatures specific to DNA alkylating agents in yeast and cancers.

Alkylation is one of the most ubiquitous forms of DNA lesions. However, the motif preferences and substrates for the activity of the major types of alkylating agents defined by their nucleophilic substitution reactions (SN1 and SN2) are still unclear. Utilizing yeast strains engineered for large-scale production of single-stranded DNA (ssDNA), we probed the substrate specificity, mutation spectra and signatures associated with DNA alkylating agents. We determined that SN1-type agents preferably mutagenize double-stranded DNA (dsDNA), and the mutation signature characteristic of the activity of SN1-type agents was conserved across yeast, mice and human cancers. Conversely, SN2-type agents preferably mutagenize ssDNA in yeast. Moreover, the spectra and signatures derived from yeast were detectable in lung cancers, head and neck cancers and tumors from patients exposed to SN2-type alkylating chemicals. The estimates of mutation loads associated with the SN2-type alkylation signature were higher in lung tumors from smokers than never-smokers, pointing toward the mutagenic activity of the SN2-type alkylating carcinogens in cigarettes. In summary, our analysis of mutations in yeast strains treated with alkylating agents, as well as in whole-exome and whole-genome-sequenced tumors identified signatures highly specific to alkylation mutagenesis and indicate the pervasive nature of alkylation-induced mutagenesis in cancers.

Repair of base damage within break-induced replication intermediates promotes kataegis associated with chromosome rearrangements.

Break induced replication (BIR) is a double strand break repair pathway that can promote genetic instabilities similar to those observed in cancer. Instead of a replication fork, BIR is driven by a migration bubble where asynchronous synthesis between leading and lagging strands leads to accumulation of single-stranded DNA (ssDNA) that promotes mutation. However, the details of the mechanism of mutagenesis, including the identity of the participating proteins, remain unknown. Using yeast as a model, we demonstrate that mutagenic ssDNA is formed at multiple positions along the BIR track and that Pol ζ is responsible for the majority of both spontaneous and damage-induced base substitutions during BIR. We also report that BIR creates a potent substrate for APOBEC3A (A3A) cytidine deaminase that can promote formation of mutation clusters along the entire track of BIR. Finally, we demonstrate that uracil glycosylase initiates the bypass of DNA damage induced by A3A in the context of BIR without formation of base substitutions, but instead this pathway frequently leads to chromosomal rearrangements. Together, the expression of A3A during BIR in yeast recapitulates the main features of APOBEC-induced kataegis in human cancers, suggesting that BIR might represent an important source of these hyper-mutagenic events.

Genome-wide maps of alkylation damage, repair, and mutagenesis in yeast reveal mechanisms of mutational heterogeneity.

DNA base damage is an important contributor to genome instability, but how the formation and repair of these lesions is affected by the genomic landscape and contributes to mutagenesis is unknown. Here, we describe genome-wide maps of DNA base damage, repair, and mutagenesis at single nucleotide resolution in yeast treated with the alkylating agent methyl methanesulfonate (MMS). Analysis of these maps revealed that base excision repair (BER) of alkylation damage is significantly modulated by chromatin, with faster repair in nucleosome-depleted regions, and slower repair and higher mutation density within strongly positioned nucleosomes. Both the translational and rotational settings of lesions within nucleosomes significantly influence BER efficiency; moreover, this effect is asymmetric relative to the nucleosome dyad axis and is regulated by histone modifications. Our data also indicate that MMS-induced mutations at adenine nucleotides are significantly enriched on the nontranscribed strand (NTS) of yeast genes, particularly in BER-deficient strains, due to higher damage formation on the NTS and transcription-coupled repair of the transcribed strand (TS). These findings reveal the influence of chromatin on repair and mutagenesis of base lesions on a genome-wide scale and suggest a novel mechanism for transcription-associated mutation asymmetry, which is frequently observed in human cancers.

Clustered mutations in yeast and in human cancers can arise from damaged long single-strand DNA regions.

Mutations are typically perceived as random, independent events. We describe here nonrandom clustered mutations in yeast and in human cancers. Genome sequencing of yeast grown under chronic alkylation damage identified mutation clusters that extend up to 200 kb. A predominance of "strand-coordinated" changes of either cytosines or guanines in the same strand, mutation patterns, and genetic controls indicated that simultaneous mutations were generated by base alkylation in abnormally long single-strand DNA (ssDNA) formed at double-strand breaks (DSBs) and replication forks. Significantly, we found mutation clusters with analogous features in sequenced human cancers. Strand-coordinated clusters of mutated cytosines or guanines often resided near chromosome rearrangement breakpoints and were highly enriched with a motif targeted by APOBEC family cytosine-deaminases, which strongly prefer ssDNA. These data indicate that hypermutation via multiple simultaneous changes in randomly formed ssDNA is a general phenomenon that may be an important mechanism producing rapid genetic variation.

The Impact of Environmental and Endogenous Damage on Somatic Mutation Load in Human Skin Fibroblasts.

Accumulation of somatic changes, due to environmental and endogenous lesions, in the human genome is associated with aging and cancer. Understanding the impacts of these processes on mutagenesis is fundamental to understanding the etiology, and improving the prognosis and prevention of cancers and other genetic diseases. Previous methods relying on either the generation of induced pluripotent stem cells, or sequencing of single-cell genomes were inherently error-prone and did not allow independent validation of the mutations. In the current study we eliminated these potential sources of error by high coverage genome sequencing of single-cell derived clonal fibroblast lineages, obtained after minimal propagation in culture, prepared from skin biopsies of two healthy adult humans. We report here accurate measurement of genome-wide magnitude and spectra of mutations accrued in skin fibroblasts of healthy adult humans. We found that every cell contains at least one chromosomal rearrangement and 600­13,000 base substitutions. The spectra and correlation of base substitutions with epigenomic features resemble many cancers. Moreover, because biopsies were taken from body parts differing by sun exposure, we can delineate the precise contributions of environmental and endogenous factors to the accrual of genetic changes within the same individual. We show here that UV-induced and endogenous DNA damage can have a comparable impact on the somatic mutation loads in skin fibroblasts. Trial Registration: ClinicalTrials.gov NCT01087307.

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.

Variant calling in low-coverage whole genome sequencing of a Native American population sample.

BACKGROUND: The reduction in the cost of sequencing a human genome has led to the use of genotype sampling strategies in order to impute and infer the presence of sequence variants that can then be tested for associations with traits of interest. Low-coverage Whole Genome Sequencing (WGS) is a sampling strategy that overcomes some of the deficiencies seen in fixed content SNP array studies. Linkage-disequilibrium (LD) aware variant callers, such as the program Thunder, may provide a calling rate and accuracy that makes a low-coverage sequencing strategy viable. RESULTS: We examined the performance of an LD-aware variant calling strategy in a population of 708 low-coverage whole genome sequences from a community sample of Native Americans. We assessed variant calling through a comparison of the sequencing results to genotypes measured in 641 of the same subjects using a fixed content first generation exome array. The comparison was made using the variant calling routines GATK Unified Genotyper program and the LD-aware variant caller Thunder. Thunder was found to improve concordance in a coverage dependent fashion, while correctly calling nearly all of the common variants as well as a high percentage of the rare variants present in the sample. CONCLUSIONS: Low-coverage WGS is a strategy that appears to collect genetic information intermediate in scope between fixed content genotyping arrays and deep-coverage WGS. Our data suggests that low-coverage WGS is a viable strategy with a greater chance of discovering novel variants and associations than fixed content arrays for large sample association analyses.

The generation of oxidative stress-induced rearrangements in Saccharomyces cerevisiae mtDNA is dependent on the Nuc1 (EndoG/ExoG) nuclease and is enhanced by inactivation of the MRX complex.

Oxidative stress is known to enhance the frequency of two major types of alterations in the mitochondrial genome of Saccharomyces cerevisiae: point mutations and large deletions resulting in the generation of respiration-deficient petite rho mutants. We investigated the effect of antimycin A, a well-known agent inducing oxidative stress, on the stability of mtDNA. We show that antimycin enhances exclusively the generation of respiration-deficient petite mutants and this is accompanied by a significant increase in the level of reactive oxygen species (ROS) and in a marked drop of cellular ATP. Whole mitochondrial genome sequencing revealed that mtDNAs of antimycin-induced petite mutants are deleted for most of the wild-type sequence and usually contain one of the active origins of mtDNA replication: ori1, ori2 ori3 or ori5. We show that the frequency of antimycin-induced rho mutants is significantly elevated in mutants deleted either for the RAD50 or XRS2 gene, both encoding the components of the MRX complex, which is known to be involved in the repair of double strand breaks (DSBs) in DNA. Furthermore, enhanced frequency of rho mutants in cultures of antimycin-treated cells lacking Rad50 was further increased by the simultaneous absence of the Ogg1 glycosylase, an important enzyme functioning in mtBER. We demonstrate also that rad50Δ and xrs2Δ deletion mutants display a considerable reduction in the frequency of allelic mitochondrial recombination, suggesting that it is the deficiency in homologous recombination which is responsible for enhanced rearrangements of mtDNA in antimycin-treated cells of these mutants. Finally, we show that the generation of large-scale mtDNA deletions induced by antimycin is markedly decreased in a nuc1Δ mutant lacking the activity of the Nuc1 nuclease, an ortholog of the mammalian mitochondrial nucleases EndoG and ExoG. This result indicates that the nuclease plays an important role in processing of oxidative stress-induced lesions in the mitochondrial genome.

Metabolic Potential, Ecology and Presence of Associated Bacteria Is Reflected in Genomic Diversity of Mucoromycotina.

Mucoromycotina are often considered mainly in pathogenic context but their biology remains understudied. We describe the genomes of six Mucoromycotina fungi representing distant saprotrophic lineages within the subphylum (i.e., Umbelopsidales and Mucorales). We selected two Umbelopsis isolates from soil (i.e., U. isabellina, U. vinacea), two soil-derived Mucor isolates (i.e., M. circinatus, M. plumbeus), and two Mucorales representatives with extended proteolytic activity (i.e., Thamnidium elegans and Mucor saturninus). We complement computational genome annotation with experimental characteristics of their digestive capabilities, cell wall carbohydrate composition, and extensive total lipid profiles. These traits inferred from genome composition, e.g., in terms of identified encoded enzymes, are in accordance with experimental results. Finally, we link the presence of associated bacteria with observed characteristics. Thamnidium elegans genome harbors an additional, complete genome of an associated bacterium classified to Paenibacillus sp. This fungus displays multiple altered traits compared to the remaining isolates, regardless of their evolutionary distance. For instance, it has expanded carbon assimilation capabilities, e.g., efficiently degrades carboxylic acids, and has a higher diacylglycerol:triacylglycerol ratio and skewed phospholipid composition which suggests a more rigid cellular membrane. The bacterium can complement the host enzymatic capabilities, alter the fungal metabolism, cell membrane composition but does not change the composition of the cell wall of the fungus. Comparison of early-diverging Umbelopsidales with evolutionary younger Mucorales points at several subtle differences particularly in their carbon source preferences and encoded carbohydrate repertoire. Nevertheless, all tested Mucoromycotina share features including the ability to produce 18:3 gamma-linoleic acid, use TAG as the storage lipid and have fucose as a cell wall component.

The Shu complex prevents mutagenesis and cytotoxicity of single-strand specific alkylation lesions.

Three-methyl cytosine (3meC) are toxic DNA lesions, blocking base pairing. Bacteria and humans express members of the AlkB enzymes family, which directly remove 3meC. However, other organisms, including budding yeast, lack this class of enzymes. It remains an unanswered evolutionary question as to how yeast repairs 3meC, particularly in single-stranded DNA. The yeast Shu complex, a conserved homologous recombination factor, aids in preventing replication-associated mutagenesis from DNA base damaging agents such as methyl methanesulfonate (MMS). We found that MMS-treated Shu complex-deficient cells exhibit a genome-wide increase in A:T and G:C substitutions mutations. The G:C substitutions displayed transcriptional and replicational asymmetries consistent with mutations resulting from 3meC. Ectopic expression of a human AlkB homolog in Shu-deficient yeast rescues MMS-induced growth defects and increased mutagenesis. Thus, our work identifies a novel homologous recombination-based mechanism mediated by the Shu complex for coping with alkylation adducts.

Msh1p counteracts oxidative lesion-induced instability of mtDNA and stimulates mitochondrial recombination in Saccharomyces cerevisiae.

The proximity of the mitochondrial genome to the respiratory chain, a major source of ROS (radical oxygen species), makes mtDNA more vulnerable to oxidative damage than nuclear DNA. Mitochondrial BER (base excision repair) is generally considered to be the main pathway involved in the prevention of oxidative lesion-induced mutations in mtDNA. However, we previously demonstrated that the increased frequency of mitochondrial Oli(r) mutants in an ogg1Delta strain, lacking the activity of a crucial mtBER glycosylase, is reduced in the presence of plasmids encoding Msh1p, the mitochondrial homologue of the bacterial mismatch protein MutS. This finding suggested that Msh1p might be involved in the prevention of mitochondrial mutagenesis induced by oxidative stress. Here we show that a double mutant carrying the msh1-R813W allele, encoding a variant of the protein defective in the ATP hydrolysis activity, combined with deletion of SOD2, encoding the mitochondrial superoxide dismutase, displays a synergistic effect on the frequency of Oli(r) mutants, indicating that Msh1p prevents generation of oxidative lesion-induced mitochondrial mutations. We also show that double mutants carrying the msh1-R813W allele, combined with deletion of either OGG1 or APN1, the latter resulting in deficiency of the Apn1 endonuclease, exhibit a synergistic effect on the frequency of respiration-defective mutants having gross rearrangements of the mitochondrial genome. This suggests that Msh1p, Ogg1p and Apn1p play overlapping functions in maintaining the stability of mtDNA. In addition, we demonstrate, using a novel ARG8(m) recombination assay, that a surplus of Msh1p results in enhanced mitochondrial recombination. Interestingly, the mutant forms of the protein, msh1p-R813W and msh1p-G776D, fail to stimulate recombination. We postulate that the Msh1p-enhanced homologous recombination may play an important role in the prevention of oxidative lesion-induced rearrangements of the mitochondrial genome.

Congenital Midline Cervical Cleft: First Report and Genetic Analysis of Two Related Patients.

OBJECTIVES: Congenital midline cervical cleft (CMCC) is a rare congenital anterior neck anatomical anomaly. We present the case of two related patients (grandchild and maternal grandmother) who were both born with a congenital midline cervical cleft along with genetic analysis. METHODS: Clinical examination of both patients and surgical excision of the grandchild was performed. Genetic analysis with exome sequencing (ES) was conducted for both patients. RESULTS: Genetic analysis with exome sequencing (ES) revealed apparently novel single nucleotide variants in 66 genes present in both proband and grandmother. Five of these variants are predicted to cause frameshifting in the coding region of the respective genes and truncated proteins (OVGP1, TYW1B, ZAN, SSPO, FOLR3). Two of these genes (TYW1B and SSPO) have homozygous indel mutations in both patients. CONCLUSIONS: To our knowledge, this is the first case of two related patients with a congenital midline cervical cleft. The results of our genetic analysis reveal potential relevance to CMCC development.

Atypical UV Photoproducts Induce Non-canonical Mutation Classes Associated with Driver Mutations in Melanoma.

Somatic mutations in skin cancers and other ultraviolet (UV)-exposed cells are typified by C>T and CC>TT substitutions at dipyrimidine sequences; however, many oncogenic "driver" mutations in melanoma do not fit this UV signature. Here, we use genome sequencing to characterize mutations in yeast repeatedly irradiated with UV light. Analysis of ~50,000 UV-induced mutations reveals abundant non-canonical mutations, including T>C, T>A, and AC>TT substitutions. These mutations display transcriptional asymmetry that is modulated by nucleotide excision repair (NER), indicating that they are caused by UV photoproducts. Using a sequencing method called UV DNA endonuclease sequencing (UVDE-seq), we confirm the existence of an atypical thymine-adenine photoproduct likely responsible for UV-induced T>A substitutions. Similar non-canonical mutations are present in skin cancers, which also display transcriptional asymmetry and dependence on NER. These include multiple driver mutations, most prominently the recurrent BRAF V600E and V600K substitutions, suggesting that mutations arising from rare, atypical UV photoproducts may play a role in melanomagenesis.

Inactivation of the 20S proteasome maturase, Ump1p, leads to the instability of mtDNA in Saccharomyces cerevisiae.

The proteasome plays fundamental roles in the removal of oxidized proteins and in normal degradation of short-lived proteins. Increasing evidence suggests that the proteasome may be an important factor in both oxidative stress response and cellular aging. Moreover, it was recently reported that proteasome inhibition leads to mitochondrial dysfunction. In this study, we have investigated whether proteasome impairment, caused by deletion of UMP1, a gene necessary for the 20S proteasome biogenesis, may influence the stability of the yeast mitochondrial genome. Here we show that an ump1Delta mutant displays enhanced mitochondrial point mutagenesis, measured by the frequency of oligomycin-resistant (Oli(r)) and erythromycin-resistant (Ery(r)) mutants, compared to that of the isogenic wild-type strain. Deletion of UMP1 significantly increases also the frequency of respiration-defective mutants having gross rearrangements of the mitochondrial genome. We show that this mitochondrial mutator phenotype of the ump1Delta strain is considerably reduced in the presence of a plasmid encoding Msh1p, the mitochondrial homologue of the bacterial mismatch protein MutS, which was shown previously to counteract oxidative lesion-induced instability of mtDNA. In search of the mechanism underlying the decreased stability of mtDNA in the ump1Delta deletion mutant, we have determined the level of reactive oxygen species (ROS) in the mutant cells and have found that they are exposed to endogenous oxidative stress. Furthermore, we show also that both cellular and intramitochondrial levels of Msh1p are significantly reduced in the mutant cells compared to the wild-type cells. We conclude, therefore, that both an increased ROS production and a markedly decreased level of Msh1p, a protein crucial for the repair of mtDNA, lead in S. cerevisiae cells with impaired proteasome activity to the increased instability of their mitochondrial genome.

A Case Study of Genomic Instability in an Industrial Strain of Saccharomyces cerevisiae.

The Saccharomyces cerevisiae strain JAY270/PE2 is a highly efficient biocatalyst used in the production of bioethanol from sugarcane feedstock. This strain is heterothallic and diploid, and its genome is characterized by abundant structural and nucleotide polymorphisms between homologous chromosomes. One of the reasons it is favored by many distilleries is that its cells do not normally aggregate, a trait that facilitates cell recycling during batch-fed fermentations. However, long-term propagation makes the yeast population vulnerable to the effects of genomic instability, which may trigger the appearance of undesirable phenotypes such as cellular aggregation. In pure cultures of JAY270, we identified the recurrent appearance of mutants displaying a mother-daughter cell separation defect resulting in rough colonies in agar media and fast sedimentation in liquid culture. We investigated the genetic basis of the colony morphology phenotype and found that JAY270 is heterozygous for a frameshift mutation in the ACE2 gene (ACE2/ace2-A7), which encodes a transcriptional regulator of mother-daughter cell separation. All spontaneous rough colony JAY270-derived isolates analyzed carried copy-neutral loss-of-heterozygosity (LOH) at the region of chromosome XII where ACE2 is located (ace2-A7/ace2-A7). We specifically measured LOH rates at the ACE2 locus, and at three additional chromosomal regions in JAY270 and in a conventional homozygous diploid laboratory strain. This direct comparison showed that LOH rates at all sites were quite similar between the two strain backgrounds. In this case study of genomic instability in an industrial strain, we showed that the JAY270 genome is dynamic and that structural changes to its chromosomes can lead to new phenotypes. However, our analysis also indicated that the inherent level of genomic instability in this industrial strain is normal relative to a laboratory strain. Our work provides an important frame of reference to contextualize the interpretation of instability processes observed in the complex genomes of industrial yeast strains.

APOBEC3B cytidine deaminase targets the non-transcribed strand of tRNA genes in yeast.

Variations in mutation rates across the genome have been demonstrated both in model organisms and in cancers. This phenomenon is largely driven by the damage specificity of diverse mutagens and the differences in DNA repair efficiency in given genomic contexts. Here, we demonstrate that the single-strand DNA-specific cytidine deaminase APOBEC3B (A3B) damages tRNA genes at a 1000-fold higher efficiency than other non-tRNA genomic regions in budding yeast. We found that A3B-induced lesions in tRNA genes were predominantly located on the non-transcribed strand, while no transcriptional strand bias was observed in protein coding genes. Furthermore, tRNA gene mutations were exacerbated in cells where RNaseH expression was completely abolished (Δrnh1Δrnh35). These data suggest a transcription-dependent mechanism for A3B-induced tRNA gene hypermutation. Interestingly, in strains proficient in DNA repair, only 1% of the abasic sites formed upon excision of A3B-deaminated cytosines were not repaired leading to mutations in tRNA genes, while 18% of these lesions failed to be repaired in the remainder of the genome. A3B-induced mutagenesis in tRNA genes was found to be efficiently suppressed by the redundant activities of both base excision repair (BER) and the error-free DNA damage bypass pathway. On the other hand, deficiencies in BER did not have a profound effect on A3B-induced mutations in CAN1, the reporter for protein coding genes. We hypothesize that differences in the mechanisms underlying ssDNA formation at tRNA genes and other genomic loci are the key determinants of the choice of the repair pathways and consequently the efficiency of DNA damage repair in these regions. Overall, our results indicate that tRNA genes are highly susceptible to ssDNA-specific DNA damaging agents. However, increased DNA repair efficacy in tRNA genes can prevent their hypermutation and maintain both genome and proteome homeostasis.