An Organismal CNV Mutator Phenotype Restricted to Early Human Development.

De novo copy number variants (dnCNVs) arising at multiple loci in a personal genome have usually been considered to reflect cancer somatic genomic instabilities. We describe a multiple dnCNV (MdnCNV) phenomenon in which individuals with genomic disorders carry five to ten constitutional dnCNVs. These CNVs originate from independent formation incidences, are predominantly tandem duplications or complex gains, exhibit breakpoint junction features reminiscent of replicative repair, and show increased de novo point mutations flanking the rearrangement junctions. The active CNV mutation shower appears to be restricted to a transient perizygotic period. We propose that a defect in the CNV formation process is responsible for the "CNV-mutator state," and this state is dampened after early embryogenesis. The constitutional MdnCNV phenomenon resembles chromosomal instability in various cancers. Investigations of this phenomenon may provide unique access to understanding genomic disorders, structural variant mutagenesis, human evolution, and cancer biology.

Comprehensive Analysis of Hypermutation in Human Cancer.

We present an extensive assessment of mutation burden through sequencing analysis of >81,000 tumors from pediatric and adult patients, including tumors with hypermutation caused by chemotherapy, carcinogens, or germline alterations. Hypermutation was detected in tumor types not previously associated with high mutation burden. Replication repair deficiency was a major contributing factor. We uncovered new driver mutations in the replication-repair-associated DNA polymerases and a distinct impact of microsatellite instability and replication repair deficiency on the scale of mutation load. Unbiased clustering, based on mutational context, revealed clinically relevant subgroups regardless of the tumors' tissue of origin, highlighting similarities in evolutionary dynamics leading to hypermutation. Mutagens, such as UV light, were implicated in unexpected cancers, including sarcomas and lung tumors. The order of mutational signatures identified previous treatment and germline replication repair deficiency, which improved management of patients and families. These data will inform tumor classification, genetic testing, and clinical trial design.

Germline MBD4 deficiency causes a multi-tumor predisposition syndrome.

We report an autosomal recessive, multi-organ tumor predisposition syndrome, caused by bi-allelic loss-of-function germline variants in the base excision repair (BER) gene MBD4. We identified five individuals with bi-allelic MBD4 variants within four families and these individuals had a personal and/or family history of adenomatous colorectal polyposis, acute myeloid leukemia, and uveal melanoma. MBD4 encodes a glycosylase involved in repair of G:T mismatches resulting from deamination of 5'-methylcytosine. The colorectal adenomas from MBD4-deficient individuals showed a mutator phenotype attributable to mutational signature SBS1, consistent with the function of MBD4. MBD4-deficient polyps harbored somatic mutations in similar driver genes to sporadic colorectal tumors, although AMER1 mutations were more common and KRAS mutations less frequent. Our findings expand the role of BER deficiencies in tumor predisposition. Inclusion of MBD4 in genetic testing for polyposis and multi-tumor phenotypes is warranted to improve disease management.

Mutator-like transposable element 9A interacts with metacaspase 1 and modulates the incidence of Al-induced programmed cell death in peanut.

The toxicity of aluminum (Al) in acidic soil inhibits plant root development and reduces crop yields. In the plant response to Al toxicity, the initiation of programmed cell death (PCD) appears to be an important mechanism for the elimination of Al-damaged cells to ensure plant survival. In a previous study, the type I metacaspase AhMC1 was found to regulate the Al stress response and to be essential for Al-induced PCD. However, the mechanism by which AhMC1 is altered in the peanut response to Al stress remained unclear. Here, we show that a nuclear protein, mutator-like transposable element 9A (AhMULE9A), directly interacts with AhMC1 in vitro and in vivo. This interaction occurs in the nucleus in peanut and is weakened during Al stress. Furthermore, a conserved C2HC zinc finger domain of AhMULE9A (residues 735-751) was shown to be required for its interaction with AhMC1. Overexpression of AhMULE9A in Arabidopsis and peanut strongly inhibited root growth with a loss of root cell viability under Al treatment. Conversely, knock down of AhMULE9A in peanut significantly reduced Al uptake and Al inhibition of root growth, and alleviated the occurrence of typical hallmarks of Al-induced PCD. These findings provide novel insight into the regulation of Al-induced PCD.

Emergence of a multilocus mutator genotype in mutator Escherichia coli experimental populations under repeated lethal selection.

Mutator alleles, which confer increased mutation rates, are known to spontaneously emerge and "hitchhike" to fixation in evolving asexual populations. Theory predicts that in an evolving asexual mutator population, a second mutator allele may spontaneously arise and hitchhike to fixation. Here, we describe an empirical test of the hypothesis of repeated hitchhiking. The starting population was a clonal strain of mutL-Escherichia coli whose mutation rate was 100-fold higher than wild type. We exposed the mutL- strain to a series of three antibiotics in increasing order of selective strength: fosfomycin, rifampicin, and streptomycin. Two independent replicates of the experiment were performed. As predicted, elevated mutation rates and enrichment for multilocus mutators (which bear more than one mutator allele) were observed in the end point populations of both experiments. DNA sequencing revealed an identical spontaneous 1-bp insertion in the mutator gene mutT in both end point populations. In the multilocus mutators, the causal relationship between the mutT- mutations and the increase in mutation rate was supported with mutT+ plasmid complementation tests. Surprisingly, when the experiment was repeated with the antibiotics deployed in decreasing order of selective strength, enrichment for multilocus mutators was not observed. Our data support the likelihood that the mutT- mutations rose to fixation in both populations, consistent with the hypothesis of repeated mutator hitchhiking. The escalation of mutation rates in asexual populations is relevant to multiple biological scenarios, including antibiotic resistance, host-pathogen interactions, and carcinogenesis.

Enhanced ROS Production in Mitochondria from Prematurely Aging mtDNA Mutator Mice.

An increase in mitochondrial DNA (mtDNA) mutations and an ensuing increase in mitochondrial reactive oxygen species (ROS) production have been suggested to be a cause of the aging process ("the mitochondrial hypothesis of aging"). In agreement with this, mtDNA-mutator mice accumulate a large amount of mtDNA mutations, giving rise to defective mitochondria and an accelerated aging phenotype. However, incongruously, the rates of ROS production in mtDNA mutator mitochondria have generally earlier been reported to be lower - not higher - than in wildtype, thus apparently invalidating the "mitochondrial hypothesis of aging". We have here re-examined ROS production rates in mtDNA-mutator mice mitochondria. Using traditional conditions for measuring ROS (succinate in the absence of rotenone), we indeed found lower ROS in the mtDNA-mutator mitochondria compared to wildtype. This ROS mainly results from reverse electron flow driven by the membrane potential, but the membrane potential reached in the isolated mtDNA-mutator mitochondria was 33 mV lower than that in wildtype mitochondria, due to the feedback inhibition of succinate oxidation by oxaloacetate, and to a lower oxidative capacity in the mtDNA-mutator mice, explaining the lower ROS production. In contrast, in normal forward electron flow systems (pyruvate (or glutamate) + malate or palmitoyl-CoA + carnitine), mitochondrial ROS production was higher in the mtDNA-mutator mitochondria. Particularly, even during active oxidative phosphorylation (as would be ongoing physiologically), higher ROS rates were seen in the mtDNA-mutator mitochondria than in wildtype. Thus, when examined under physiological conditions, mitochondrial ROS production rates are indeed increased in mtDNA-mutator mitochondria. While this does not prove the validity of the mitochondrial hypothesis of aging, it may no longer be said to be negated in this respect. This paper is dedicated to the memory of Professor Vladimir P. Skulachev.

Detection and characterization of constitutive replication origins defined by DNA polymerase epsilon.

BACKGROUND: Despite the process of DNA replication being mechanistically highly conserved, the location of origins of replication (ORI) may vary from one tissue to the next, or between rounds of replication in eukaryotes, suggesting flexibility in the choice of locations to initiate replication. Lists of human ORI therefore vary widely in number and location, and there are currently no methods available to compare them. Here, we propose a method of detection of ORI based on somatic mutation patterns generated by the mutator phenotype of damaged DNA polymerase epsilon (POLE). RESULTS: We report the genome-wide localization of constitutive ORI in POLE-mutated human tumors using whole genome sequencing data. Mutations accumulated after many rounds of replication of unsynchronized dividing cell populations in tumors allow to identify constitutive origins, which we show are shared with high fidelity between individuals and tumor types. Using a Smith-Waterman-like dynamic programming approach, we compared replication origin positions obtained from multiple different methods. The comparison allowed us to define a consensus set of replication origins, identified consistently by multiple ORI detection methods. Many DNA features co-localized with the consensus set of ORI, including chromatin loop anchors, G-quadruplexes, S/MARs, and CpGs. Among all features, the H2A.Z histone exhibited the most significant association. CONCLUSIONS: Our results show that mutation-based detection of replication origins is a viable approach to determining their location and associated sequence features.

Expansion and improvement of ChinaMu by MuT-seq and chromosome-level assembly of the Mu-starter genome.

ChinaMu is the largest sequence-indexed Mutator (Mu) transposon insertional library in maize (Zea mays). In this study, we made significant improvements to the size and quality of the ChinaMu library. We developed a new Mu-tag isolation method Mu-Tn5-seq (MuT-seq). Compared to the previous method used by ChinaMu, MuT-seq recovered 1/3 more germinal insertions, while requiring only about 1/14 of the sequencing volume and 1/5 of the experimental time. Using MuT-seq, we identified 113,879 germinal insertions from 3,168 Mu-active F1 families. We also assembled a high-quality genome for the Mu-active line Mu-starter, which harbors the initial active MuDR element and was used as the pollen donor for the mutation population. Using the Mu-starter genome, we recovered 33,662 (15.6%) additional germinal insertions in 3,244 (7.4%) genes in the Mu-starter line. The Mu-starter genome also improved the assignment of 117,689 (54.5%) germinal insertions. The newly upgraded ChinaMu dataset currently contains 215,889 high-quality germinal insertions. These insertions cover 32,224 pan-genes in the Mu-starter and B73Ref5 genomes, including 23,006 (80.4%) core genes shared by the two genomes. As a test model, we investigated Mu insertions in the pentatricopeptide repeat (PPR) superfamily, discovering insertions for 92% (449/487) of PPR genes in ChinaMu, demonstrating the usefulness of ChinaMu as a functional genomics resource for maize.

Shifts in Mutation Bias Promote Mutators by Altering the Distribution of Fitness Effects.

AbstractRecent experimental evidence demonstrates that shifts in mutational biases-for example, increases in transversion frequency-can change the distribution of fitness effects of mutations (DFE). In particular, reducing or reversing a prevailing bias can increase the probability that a de novo mutation is beneficial. It has also been shown that mutator bacteria are more likely to emerge if the beneficial mutations they generate have a larger effect size than observed in the wild type. Here, we connect these two results, demonstrating that mutator strains that reduce or reverse a prevailing bias have a positively shifted DFE, which in turn can dramatically increase their emergence probability. Since changes in mutation rate and bias are often coupled through the gain and loss of DNA repair enzymes, our results predict that the invasion of mutator strains will be facilitated by shifts in mutation bias that offer improved access to previously undersampled beneficial mutations.

Characterization of POLE c.1373A > T p.(Tyr458Phe), causing high cancer risk.

The cancer syndrome polymerase proofreading-associated polyposis results from germline mutations in the POLE and POLD1 genes. Mutations in the exonuclease domain of these genes are associated with hyper- and ultra-mutated tumors with a predominance of base substitutions resulting from faulty proofreading during DNA replication. When a new variant is identified by gene testing of POLE and POLD1, it is important to verify whether the variant is associated with PPAP or not, to guide genetic counseling of mutation carriers. In 2015, we reported the likely pathogenic (class 4) germline POLE c.1373A > T p.(Tyr458Phe) variant and we have now characterized this variant to verify that it is a class 5 pathogenic variant. For this purpose, we investigated (1) mutator phenotype in tumors from two carriers, (2) mutation frequency in cell-based mutagenesis assays, and (3) structural consequences based on protein modeling. Whole-exome sequencing of two tumors identified an ultra-mutator phenotype with a predominance of base substitutions, the majority of which are C > T. A SupF mutagenesis assay revealed increased mutation frequency in cells overexpressing the variant of interest as well as in isogenic cells encoding the variant. Moreover, exonuclease repair yeast-based assay supported defect in proofreading activity. Lastly, we present a homology model of human POLE to demonstrate structural consequences leading to pathogenic impact of the p.(Tyr458Phe) mutation. The three lines of evidence, taken together with updated co-segregation and previously published data, allow the germline variant POLE c.1373A > T p.(Tyr458Phe) to be reclassified as a class 5 variant. That means the variant is associated with PPAP.

Eukaryotic Mismatch Repair in Relation to DNA Replication.

Three processes act in series to accurately replicate the eukaryotic nuclear genome. The major replicative DNA polymerases strongly prevent mismatch formation, occasional mismatches that do form are proofread during replication, and rare mismatches that escape proofreading are corrected by mismatch repair (MMR). This review focuses on MMR in light of increasing knowledge about nuclear DNA replication enzymology and the rate and specificity with which mismatches are generated during leading- and lagging-strand replication. We consider differences in MMR efficiency in relation to mismatch recognition, signaling to direct MMR to the nascent strand, mismatch removal, and the timing of MMR. These studies are refining our understanding of relationships between generating and repairing replication errors to achieve accurate replication of both DNA strands of the nuclear genome.

Trajectory and uniqueness of mutational signatures in yeast mutators.

The acquisition of mutations plays critical roles in adaptation, evolution, senescence, and tumorigenesis. Massive genome sequencing has allowed extraction of specific features of many mutational landscapes but it remains difficult to retrospectively determine the mechanistic origin(s), selective forces, and trajectories of transient or persistent mutations and genome rearrangements. Here, we conducted a prospective reciprocal approach to inactivate 13 single or multiple evolutionary conserved genes involved in distinct genome maintenance processes and characterize de novo mutations in 274 diploid Saccharomyces cerevisiae mutation accumulation lines. This approach revealed the diversity, complexity, and ultimate uniqueness of mutational landscapes, differently composed of base substitutions, small insertions/deletions (InDels), structural variants, and/or ploidy variations. Several landscapes parallel the repertoire of mutational signatures in human cancers while others are either novel or composites of subsignatures resulting from distinct DNA damage lesions. Notably, the increase of base substitutions in the homologous recombination-deficient Rad51 mutant, specifically dependent on the Polζ translesion polymerase, yields COSMIC signature 3 observed in BRCA1/BRCA2-mutant breast cancer tumors. Furthermore, "mutome" analyses in highly polymorphic diploids and single-cell bottleneck lineages revealed a diverse spectrum of loss-of-heterozygosity (LOH) signatures characterized by interstitial and terminal chromosomal events resulting from interhomolog mitotic cross-overs. Following the appearance of heterozygous mutations, the strong stimulation of LOHs in the rad27/FEN1 and tsa1/PRDX1 backgrounds leads to fixation of homozygous mutations or their loss along the lineage. Overall, these mutomes and their trajectories provide a mechanistic framework to understand the origin and dynamics of genome variations that accumulate during clonal evolution.

Computational identification of mutator-derived lncRNA signatures of genome instability for improving the clinical outcome of cancers: a case study in breast cancer.

Emerging evidence revealed the critical roles of long non-coding RNAs (lncRNAs) in maintaining genomic instability. However, identification of genome instability-associated lncRNAs and their clinical significance in cancers remain largely unexplored. Here, we developed a mutator hypothesis-derived computational frame combining lncRNA expression profiles and somatic mutation profiles in a tumor genome and identified 128 novel genomic instability-associated lncRNAs in breast cancer as a case study. We then identified a genome instability-derived two lncRNA-based gene signature (GILncSig) that stratified patients into high- and low-risk groups with significantly different outcome and was further validated in multiple independent patient cohorts. Furthermore, the GILncSig correlated with genomic mutation rate in both ovarian cancer and breast cancer, indicating its potential as a measurement of the degree of genome instability. The GILncSig was able to divide TP53 wide-type patients into two risk groups, with the low-risk group showing significantly improved outcome and the high-risk group showing no significant difference compared with those with TP53 mutation. In summary, this study provided a critical approach and resource for further studies examining the role of lncRNAs in genome instability and introduced a potential new avenue for identifying genomic instability-associated cancer biomarkers.

The Barley Chloroplast Mutator (cpm) Mutant: All Roads Lead to the Msh1 Gene.

The barley chloroplast mutator (cpm) is a nuclear gene mutant that induces a wide spectrum of cytoplasmically inherited chlorophyll deficiencies. Plastome instability of cpm seedlings was determined by identification of a particular landscape of polymorphisms that suggests failures in a plastome mismatch repair (MMR) protein. In Arabidopsis, MSH genes encode proteins that are in charge of mismatch repair and have anti-recombination activity. In this work, barley homologs of these genes were identified, and their sequences were analyzed in control and cpm mutant seedlings. A substitution, leading to a premature stop codon and a truncated MSH1 protein, was identified in the Msh1 gene of cpm plants. The relationship between this mutation and the presence of chlorophyll deficiencies was established in progenies from crosses and backcrosses. These results strongly suggest that the mutation identified in the Msh1 gene of the cpm mutant is responsible for the observed plastome instabilities. Interestingly, comparison of mutant phenotypes and molecular changes induced by the barley cpm mutant with those of Arabidopsis MSH1 mutants revealed marked differences.

Cloning of Maize TED Transposon into Escherichia coli Reveals the Polychromatic Sequence Landscape of Refractorily Propagated Plasmids.

MuDR, the founder member of the Mutator superfamily and its MURA transcripts, has been identified as toxic sequences to Escherichia coli (E. coli), which heavily hindered the elucidation of the biochemical features of MURA transposase and confined the broader application of the Mutator system in other organisms. To harness less constrained systems as alternatives, we attempted to clone TED and Jittery, two recently isolated autonomous Mutator-like elements (MULEs) from maize, respectively. Their full-length transcripts and genomic copies are successfully cloned when the incubation time for bacteria to recover from heat shock is extended appropriately prior to plating. However, during their proliferation in E. coli, TED transformed plasmids are unstable, as evidenced by derivatives from which frameshift, deletion mutations, or IS transposon insertions are readily detected. Our results suggest that neither leaky expression of the transposase nor the presence of terminal inverse repeats (TIRs) are responsible for the cloning barriers, which were once ascribed to the presence of the Shine-Dalgarno-like sequence. Instead, the internal sequence of TED (from 1250 to 2845 bp), especially the exons in this region, was the most likely causer. The findings provide novel insights into the property and function of the Mutator superfamily and shed light on the dissection of toxic effects on cloning from MULEs.

Mutations Beget More Mutations-Rapid Evolution of Mutation Rate in Response to the Risk of Runaway Accumulation.

The rapidity with which the mutation rate evolves could greatly impact evolutionary patterns. Nevertheless, most studies simply assume a constant rate in the time scale of interest (Kimura 1983; Drake 1991; Kumar 2005; Li 2007; Lynch 2010). In contrast, recent studies of somatic mutations suggest that the mutation rate may vary by several orders of magnitude within a lifetime (Kandoth et al. 2013; Lawrence et al. 2013). To resolve the discrepancy, we now propose a runaway model, applicable to both the germline and soma, whereby mutator mutations form a positive-feedback loop. In this loop, any mutator mutation would increase the rate of acquiring the next mutator, thus triggering a runaway escalation in mutation rate. The process can be initiated more readily if there are many weak mutators than a few strong ones. Interestingly, even a small increase in the mutation rate at birth could trigger the runaway process, resulting in unfit progeny. In slowly reproducing species, the need to minimize the risk of this uncontrolled accumulation would thus favor setting the mutation rate low. In comparison, species that starts and ends reproduction sooner do not face the risk and may set the baseline mutation rate higher. The mutation rate would evolve in response to the risk of runaway mutation, in particular, when the generation time changes. A rapidly evolving mutation rate may shed new lights on many evolutionary phenomena (Elango et al. 2006; Thomas et al. 2010, 2018; Langergraber et al. 2012; Besenbacher et al. 2019).

Calculated identification of mutator-derived lncRNA signatures of genomic instability to predict the clinical outcome of muscle-invasive bladder cancer.

BACKGROUND: Muscle-invasive bladder cancer (MIBC) is one of the most important type of bladder cancer, with a high morbidity and mortality rate. Studies have found that long non-coding RNA (lncRNA) plays a key role in maintaining genomic instability. However, Identification of lncRNAs related to genomic instability (GIlncRNAs) and their clinical significance in cancers have not been extensively studied yet. METHODS: Here, we downloaded the lncRNA expression profiles, somatic mutation profiles and clinical related data in MIBC patients from The Cancer Genome Atlas (TCGA) database. A lncRNA computational framework was used to find differentially expressed GIlncRNAs. Multivariate Cox regression analysis was used to construct a genomic instability-related lncRNA signature (GIlncSig). Univariate and multivariate Cox analyses were used to assess the independent prognostic for the GIlncSig and other key clinical factors. RESULTS: We found 43 differentially expressed GIlncRNAs and constructed the GIlncSig with 6 GIlncRNAs in the training cohort. The patients were divided into two risk groups. The overall survival of patients in the high-risk group was lower than that in the low-risk group (P < 0.001), which were further verified in the testing cohort and the entire TCGA cohort. Univariate and multivariate Cox regression showed that the GIlncSig was an independent prognostic factor. In addition, the GIlncSig correlated with the genomic mutation rate of MIBC, indicating its potential as a measure of the degree of genomic instability. The GIlncSig was able to divide FGFR3 wild- and mutant-type patients into two risk groups, and effectively enhanced the prediction effect. CONCLUSION: Our study introduced an important reference for further research on the role of GIlncRNAs, and provided prognostic indicators and potential biological therapy targets for MIBC.

Mutation frequency of Escherichia coli isolated from river water: potential role in the development of antimicrobial resistance.

Bacteria acquire genetic variations that help them to adapt to stressful environmental conditions, and these changes may be associated with the development of antimicrobial resistance. In this study, we investigated the mutation frequencies of 270 isolates of Escherichia coli from river water, which represents a relatively unstressful environment. As we predicted, mutation frequencies of the E. coli isolates ranged from <1 × 10-11 to 6.3 × 10-8 (median, 1.7 × 10-9), and a strong mutator (≥ 4 × 10-7) was not detected. To better understand the role of mutation frequency in the development of antimicrobial resistance, we assessed antimicrobial sensitivity after exposure of the E. coli isolates to subinhibitory concentrations of ciprofloxacin, as a surrogate for stress. We found that antimicrobial resistance increased in bacteria with a low mutation frequency after exposure, and the relative increase in antimicrobial resistance generally increased, depending on the mutation frequency. Thus, mutation frequency may contribute to the development of antimicrobial resistance of bacteria in natural environments.

Tyr82 Amino Acid Mutation in PB1 Polymerase Induces an Influenza Virus Mutator Phenotype.

In various positive-sense single-stranded RNA viruses, a low-fidelity viral RNA-dependent RNA polymerase (RdRp) confers attenuated phenotypes by increasing the mutation frequency. We report a negative-sense single-stranded RNA virus RdRp mutant strain with a mutator phenotype. Based on structural data of RdRp, rational targeting of key residues, and screening of fidelity variants, we isolated a novel low-fidelity mutator strain of influenza virus that harbors a Tyr82-to-Cys (Y82C) single-amino-acid substitution in the PB1 polymerase subunit. The purified PB1-Y82C polymerase indeed showed an increased frequency of misincorporation compared with the wild-type PB1 in an in vitro biochemical assay. To further investigate the effects of position 82 on PB1 polymerase fidelity, we substituted various amino acids at this position. As a result, we isolated various novel mutators other than PB1-Y82C with higher mutation frequencies. The structural model of influenza virus polymerase complex suggested that the Tyr82 residue, which is located at the nucleoside triphosphate entrance tunnel, may influence a fidelity checkpoint. Interestingly, although the PB1-Y82C variant replicated with wild-type PB1-like kinetics in tissue culture, the 50% lethal dose of the PB1-Y82C mutant was 10 times lower than that of wild-type PB1 in embryonated chicken eggs. In conclusion, our data indicate that the Tyr82 residue of PB1 has a crucial role in regulating polymerase fidelity of influenza virus and is closely related to attenuated pathogenic phenotypes in vivoIMPORTANCE Influenza A virus rapidly acquires antigenic changes and antiviral drug resistance, which limit the effectiveness of vaccines and drug treatments, primarily owing to its high rate of evolution. Virus populations formed by quasispecies can contain resistance mutations even before a selective pressure is applied. To study the effects of the viral mutation spectrum and quasispecies, high- and low-fidelity variants have been isolated for several RNA viruses. Here, we report the discovery of a low-fidelity RdRp variant of influenza A virus that contains a substitution at Tyr82 in PB1. Viruses containing the PB1-Y82C substitution showed growth kinetics and viral RNA synthesis levels similar to those of the wild-type virus in cell culture; however, they had significantly attenuated phenotypes in a chicken egg infection experiment. These data demonstrated that decreased RdRp fidelity attenuates influenza A virus in vivo, which is a desirable feature for the development of safer live attenuated vaccine candidates.

Hypermutation signature reveals a slippage and realignment model of translesion synthesis by Rev3 polymerase in cisplatin-treated yeast.

Gene-gene or gene-drug interactions are typically quantified using fitness as a readout because the data are continuous and easily measured in high throughput. However, to what extent fitness captures the range of other phenotypes that show synergistic effects is usually unknown. Using Saccharomyces cerevisiae and focusing on a matrix of DNA repair mutants and genotoxic drugs, we quantify 76 gene-drug interactions based on both mutation rate and fitness and find that these parameters are not connected. Independent of fitness defects, we identified six cases of synthetic hypermutation, where the combined effect of the drug and mutant on mutation rate was greater than predicted. One example occurred when yeast lacking RAD1 were exposed to cisplatin, and we characterized this interaction using whole-genome sequencing. Our sequencing results indicate mutagenesis by cisplatin in rad1Δ cells appeared to depend almost entirely on interstrand cross-links at GpCpN motifs. Interestingly, our data suggest that the following base on the template strand dictates the addition of the mutated base. This result differs from cisplatin mutation signatures in XPF-deficient Caenorhabditis elegans and supports a model in which translesion synthesis polymerases perform a slippage and realignment extension across from the damaged base. Accordingly, DNA polymerase ζ activity was essential for mutagenesis in cisplatin-treated rad1Δ cells. Together these data reveal the potential to gain new mechanistic insights from nonfitness measures of gene-drug interactions and extend the use of mutation accumulation and whole-genome sequencing analysis to define DNA repair mechanisms.