Frequencies and spectra of aflatoxin B1-induced mutations in liver genomes of NEIL1-deficient mice as revealed by duplex sequencing.

Increased risk for the development of hepatocellular carcinoma (HCC) is driven by a number of etiological factors including hepatitis viral infection and dietary exposures to foods contaminated with aflatoxin-producing molds. Intracellular metabolic activation of aflatoxin B1 (AFB1) to a reactive epoxide generates highly mutagenic AFB1-Fapy-dG adducts. Previously, we demonstrated that repair of AFB1-Fapy-dG adducts can be initiated by the DNA glycosylase NEIL1 and that male Neil1-/- mice were significantly more susceptible to AFB1-induced HCC relative to wild-type mice. To investigate the mechanisms underlying this enhanced carcinogenesis, WT and Neil1-/- mice were challenged with a single, 4 mg/kg dose of AFB1 and frequencies and spectra of mutations were analyzed in liver DNAs 2.5 months post-injection using duplex sequencing. The analyses of DNAs from AFB1-challenged mice revealed highly elevated mutation frequencies in the nuclear genomes of both males and females, but not the mitochondrial genomes. In both WT and Neil1-/- mice, mutation spectra were highly similar to the AFB1-specific COSMIC signature SBS24. Relative to wild-type, the NEIL1 deficiency increased AFB1-induced mutagenesis with concomitant elevated HCCs in male Neil1-/- mice. Our data establish a critical role of NEIL1 in limiting AFB1-induced mutagenesis and ultimately carcinogenesis.

Targeted accurate RNA consensus sequencing (tARC-seq) reveals mechanisms of replication error affecting SARS-CoV-2 divergence.

RNA viruses, like SARS-CoV-2, depend on their RNA-dependent RNA polymerases (RdRp) for replication, which is error prone. Monitoring replication errors is crucial for understanding the virus's evolution. Current methods lack the precision to detect rare de novo RNA mutations, particularly in low-input samples such as those from patients. Here we introduce a targeted accurate RNA consensus sequencing method (tARC-seq) to accurately determine the mutation frequency and types in SARS-CoV-2, both in cell culture and clinical samples. Our findings show an average of 2.68 × 10-5 de novo errors per cycle with a C > T bias that cannot be solely attributed to APOBEC editing. We identified hotspots and cold spots throughout the genome, correlating with high or low GC content, and pinpointed transcription regulatory sites as regions more susceptible to errors. tARC-seq captured template switching events including insertions, deletions and complex mutations. These insights shed light on the genetic diversity generation and evolutionary dynamics of SARS-CoV-2.

The multi-tissue landscape of somatic mtDNA mutations indicates tissue-specific accumulation and removal in aging.

Accumulation of somatic mutations in the mitochondrial genome (mtDNA) has long been proposed as a possible mechanism of mitochondrial and tissue dysfunction that occurs during aging. A thorough characterization of age-associated mtDNA somatic mutations has been hampered by the limited ability to detect low-frequency mutations. Here, we used Duplex Sequencing on eight tissues of an aged mouse cohort to detect >89,000 independent somatic mtDNA mutations and show significant tissue-specific increases during aging across all tissues examined which did not correlate with mitochondrial content and tissue function. G→A/C→T substitutions, indicative of replication errors and/or cytidine deamination, were the predominant mutation type across all tissues and increased with age, whereas G→T/C→A substitutions, indicative of oxidative damage, were the second most common mutation type, but did not increase with age regardless of tissue. We also show that clonal expansions of mtDNA mutations with age is tissue- and mutation type-dependent. Unexpectedly, mutations associated with oxidative damage rarely formed clones in any tissue and were significantly reduced in the hearts and kidneys of aged mice treated at late age with elamipretide or nicotinamide mononucleotide. Thus, the lack of accumulation of oxidative damage-linked mutations with age suggests a life-long dynamic clearance of either the oxidative lesions or mtDNA genomes harboring oxidative damage.

Resistance of mitochondrial DNA to cadmium and Aflatoxin B1 damage-induced germline mutation accumulation in C. elegans.

Mitochondrial DNA (mtDNA) is prone to mutation in aging and over evolutionary time, yet the processes that regulate the accumulation of de novo mtDNA mutations and modulate mtDNA heteroplasmy are not fully elucidated. Mitochondria lack certain DNA repair processes, which could contribute to polymerase error-induced mutations and increase susceptibility to chemical-induced mtDNA mutagenesis. We conducted error-corrected, ultra-sensitive Duplex Sequencing to investigate the effects of two known nuclear genome mutagens, cadmium and Aflatoxin B1, on germline mtDNA mutagenesis in Caenorhabditis elegans. Detection of thousands of mtDNA mutations revealed pervasive heteroplasmy in C. elegans and that mtDNA mutagenesis is dominated by C:G → A:T mutations generally attributed to oxidative damage. However, there was no effect of either exposure on mtDNA mutation frequency, spectrum, or trinucleotide context signature despite a significant increase in nuclear mutation rate after aflatoxin B1 exposure. Mitophagy-deficient mutants pink-1 and dct-1 accumulated significantly higher levels of mtDNA damage compared to wild-type C. elegans after exposures. However, there were only small differences in mtDNA mutation frequency, spectrum, or trinucleotide context signature compared to wild-type after 3050 generations, across all treatments. These findings suggest mitochondria harbor additional previously uncharacterized mechanisms that regulate mtDNA mutational processes across generations.

Myometrial oxidative stress drives MED12 mutations in leiomyoma.

BACKGROUND: More than 70% of leiomyomas (LM) harbor MED12 mutations, primarily in exon 2 at c.130-131(GG). The cause of MED12 mutations in myometrial cells remains largely unknown. We hypothesized that increased ROS promotes MED12 mutations in myometrial cells through the oxidation of guanine nucleotides followed by misrepair. METHODS: Genomic oxidative burden (8-OHdG) was evaluated in vitro and in vivo by immunohistochemistry. MED12 mutations were examined by Sanger sequencing and deep sequencing. Transcriptome examined by RNA-seq was performed in myometrium with and without LM, in primary myometrial cells treated with ROS. 8-OHdG mediated misrepair was analyzed by CRISPR/Cas9. RESULTS: Uteri with high LM burden had a significantly higher rate of MED12 mutations than uteri with low LM burden. Compelling data suggest that the uterus normally produces reactive oxidative species (ROS) in response to stress, and ROS levels in LM are elevated due to metabolic defects. We demonstrated that genomic oxidized guanine (8-OHdG) was found at a significantly higher level in the myometrium of uteri that had multiple LM compared to myometrium without LM. Transcriptome and pathway analyses detected ROS stress in myometrium with LM. Targeted replacement of guanine with 8-OHdG at MED12 c.130 by CRISPR/Cas9 significantly increased the misrepair of G>T. Exposure of primary myometrial cells to oxidative stress in vitro increased misrepair/mutations as detected by duplex sequencing. CONCLUSIONS: Together, our data identified a clear connection between increased myometrial oxidative stress and a high rate of MED12 mutations that may underlie the risk of LM development and severity in women of reproductive age.

RNA polymerase inaccuracy underlies SARS-CoV-2 variants and vaccine heterogeneity.

Both the SARS-CoV-2 virus and its mRNA vaccines depend on RNA polymerases (RNAP)1,2; however, these enzymes are inherently error-prone and can introduce variants into the RNA3. To understand SARS-CoV-2 evolution and vaccine efficacy, it is critical to identify the extent and distribution of errors introduced by the RNAPs involved in each process. Current methods lack the sensitivity and specificity to measure de novo RNA variants in low input samples like viral isolates3. Here, we determine the frequency and nature of RNA errors in both SARS-CoV-2 and its vaccine using a targeted Accurate RNA Consensus sequencing method (tARC-seq). We found that the viral RNA-dependent RNAP (RdRp) makes ~1 error every 10,000 nucleotides - higher than previous estimates4. We also observed that RNA variants are not randomly distributed across the genome but are associated with certain genomic features and genes, such as S (Spike). tARC-seq captured a number of large insertions, deletions and complex mutations that can be modeled through non-programmed RdRp template switching. This template switching feature of RdRp explains many key genetic changes observed during the evolution of different lineages worldwide, including Omicron. Further sequencing of the Pfizer-BioNTech COVID-19 vaccine revealed an RNA variant frequency of ~1 in 5,000, meaning most of the vaccine transcripts produced in vitro by T7 phage RNAP harbor a variant. These results demonstrate the extraordinary genetic diversity of viral populations and the heterogeneous nature of an mRNA vaccine fueled by RNAP inaccuracy. Along with functional studies and pandemic data, tARC-seq variant spectra can inform models to predict how SARS-CoV-2 may evolve. Finally, our results may help improve future vaccine development and study design as mRNA therapies continue to gain traction.

The Complicated Nature of Somatic mtDNA Mutations in Aging.

Mitochondria are the main source of energy used to maintain cellular homeostasis. This aspect of mitochondrial biology underlies their putative role in age-associated tissue dysfunction. Proper functioning of the electron transport chain (ETC), which is partially encoded by the extra-nuclear mitochondrial genome (mtDNA), is key to maintaining this energy production. The acquisition of de novo somatic mutations that interrupt the function of the ETC have long been associated with aging and common diseases of the elderly. Yet, despite over 30 years of study, the exact role(s) mtDNA mutations play in driving aging and its associated pathologies remains under considerable debate. Furthermore, even fundamental aspects of age-related mtDNA mutagenesis, such as when mutations arise during aging, where and how often they occur across tissues, and the specific mechanisms that give rise to them, remain poorly understood. In this review, we address the current understanding of the somatic mtDNA mutations, with an emphasis of when, where, and how these mutations arise during aging. Additionally, we highlight current limitations in our knowledge and critically evaluate the controversies stemming from these limitations. Lastly, we highlight new and emerging technologies that offer potential ways forward in increasing our understanding of somatic mtDNA mutagenesis in the aging process.

Probing the mechanisms of two exonuclease domain mutators of DNA polymerase ϵ.

We report the properties of two mutations in the exonuclease domain of the Saccharomyces cerevisiae DNA polymerase ϵ. One, pol2-Y473F, increases the mutation rate by about 20-fold, similar to the catalytically dead pol2-D290A/E290A mutant. The other, pol2-N378K, is a stronger mutator. Both retain the ability to excise a nucleotide from double-stranded DNA, but with impaired activity. pol2-Y473F degrades DNA poorly, while pol2-N378K degrades single-stranded DNA at an elevated rate relative to double-stranded DNA. These data suggest that pol2-Y473F reduces the capacity of the enzyme to perform catalysis in the exonuclease active site, while pol2-N378K impairs partitioning to the exonuclease active site. Relative to wild-type Pol ϵ, both variants decrease the dNTP concentration required to elicit a switch between proofreading and polymerization by more than an order of magnitude. While neither mutation appears to alter the sequence specificity of polymerization, the N378K mutation stimulates polymerase activity, increasing the probability of incorporation and extension of a mismatch. Considered together, these data indicate that impairing the primer strand transfer pathway required for proofreading increases the probability of common mutations by Pol ϵ, elucidating the association of homologous mutations in human DNA polymerase ϵ with cancer.

A replication-linked mutational gradient drives somatic mutation accumulation and influences germline polymorphisms and genome composition in mitochondrial DNA.

Mutations in mitochondrial DNA (mtDNA) cause maternally inherited diseases, while somatic mutations are linked to common diseases of aging. Although mtDNA mutations impact health, the processes that give rise to them are under considerable debate. To investigate the mechanism by which de novo mutations arise, we analyzed the distribution of naturally occurring somatic mutations across the mouse and human mtDNA obtained by Duplex Sequencing. We observe distinct mutational gradients in G→A and T→C transitions delimited by the light-strand origin and the mitochondrial Control Region (mCR). The gradient increases unequally across the mtDNA with age and is lost in the absence of DNA polymerase γ proofreading activity. In addition, high-resolution analysis of the mCR shows that important regulatory elements exhibit considerable variability in mutation frequency, consistent with them being mutational 'hot-spots' or 'cold-spots'. Collectively, these patterns support genome replication via a deamination prone asymmetric strand-displacement mechanism as the fundamental driver of mutagenesis in mammalian DNA. Moreover, the distribution of mtDNA single nucleotide polymorphisms in humans and the distribution of bases in the mtDNA across vertebrate species mirror this gradient, indicating that replication-linked mutations are likely the primary source of inherited polymorphisms that, over evolutionary timescales, influences genome composition during speciation.

PolyG-DS: An ultrasensitive polyguanine tract-profiling method to detect clonal expansions and trace cell lineage.

Polyguanine tracts (PolyGs) are short guanine homopolymer repeats that are prone to accumulating mutations when cells divide. This feature makes them especially suitable for cell lineage tracing, which has been exploited to detect and characterize precancerous and cancerous somatic evolution. PolyG genotyping, however, is challenging because of the inherent biochemical difficulties in amplifying and sequencing repetitive regions. To overcome this limitation, we developed PolyG-DS, a next-generation sequencing (NGS) method that combines the error-correction capabilities of duplex sequencing (DS) with enrichment of PolyG loci using CRISPR-Cas9-targeted genomic fragmentation. PolyG-DS markedly reduces technical artifacts by comparing the sequences derived from the complementary strands of each original DNA molecule. We demonstrate that PolyG-DS genotyping is accurate, reproducible, and highly sensitive, enabling the detection of low-frequency alleles (<0.01) in spike-in samples using a panel of only 19 PolyG markers. PolyG-DS replicated prior results based on PolyG fragment length analysis by capillary electrophoresis, and exhibited higher sensitivity for identifying clonal expansions in the nondysplastic colon of patients with ulcerative colitis. We illustrate the utility of this method for resolving the phylogenetic relationship among precancerous lesions in ulcerative colitis and for tracing the metastatic dissemination of ovarian cancer. PolyG-DS enables the study of tumor evolution without prior knowledge of tumor driver mutations and provides a tool to perform cost-effective and easily scalable ultra-accurate NGS-based PolyG genotyping for multiple applications in biology, genetics, and cancer research.

Rate volatility and asymmetric segregation diversify mutation burden in cells with mutator alleles.

Mutations that compromise mismatch repair (MMR) or DNA polymerase ε or δ exonuclease domains produce mutator phenotypes capable of fueling cancer evolution. Here, we investigate how combined defects in these pathways expands genetic heterogeneity in cells of the budding yeast, Saccharomyces cerevisiae, using a single-cell resolution approach that tallies all mutations arising from individual divisions. The distribution of replication errors present in mother cells after the initial S-phase was broader than expected for a single uniform mutation rate across all cell divisions, consistent with volatility of the mutator phenotype. The number of mismatches that then segregated to the mother and daughter cells co-varied, suggesting that each division is governed by a different underlying genome-wide mutation rate. The distribution of mutations that individual cells inherit after the second S-phase is further broadened by the sequential actions of semiconservative replication and mitotic segregation of chromosomes. Modeling suggests that this asymmetric segregation may diversify mutation burden in mutator-driven tumors.

Analyzing population structure for forensic STR markers in next generation sequencing data.

Match probabilities calculated during the evaluation of DNA evidence profiles rely on appropriate values of the population structure quantity θ. NGS-based methods will enhance forensic identification and with the transformation to such methods comes the need to facilitate NGS-based population genetics analysis. If NGS data are to be used for match probabilities there needs to be a way to accommodate population structure, which requires values for θ for those data. Such estimates have not been available. This study assesses population structure for sequence-based data using a relatively new approach applied to STR data over 27 loci in five different geographic groups. Matching proportions between individuals or groups are used to obtain locus-specific θ estimates as well as estimates per geographic group and a global measure. The results demonstrate similar effects of sequencing data on θ estimates compared to what has been seen for CE-based results.

Loss of vacuolar acidity results in iron-sulfur cluster defects and divergent homeostatic responses during aging in Saccharomyces cerevisiae.

The loss of vacuolar/lysosomal acidity is an early event during aging that has been linked to mitochondrial dysfunction. However, it is unclear how loss of vacuolar acidity results in age-related dysfunction. Through unbiased genetic screens, we determined that increased iron uptake can suppress the mitochondrial respiratory deficiency phenotype of yeast vma mutants, which have lost vacuolar acidity due to genetic disruption of the vacuolar ATPase proton pump. Yeast vma mutants exhibited nuclear localization of Aft1, which turns on the iron regulon in response to iron-sulfur cluster (ISC) deficiency. This led us to find that loss of vacuolar acidity with age in wild-type yeast causes ISC defects and a DNA damage response. Using microfluidics to investigate aging at the single-cell level, we observe grossly divergent trajectories of iron homeostasis within an isogenic and environmentally homogeneous population. One subpopulation of cells fails to mount the expected compensatory iron regulon gene expression program, and suffers progressively severe ISC deficiency with little to no activation of the iron regulon. In contrast, other cells show robust iron regulon activity with limited ISC deficiency, which allows extended passage and survival through a period of genomic instability during aging. These divergent trajectories suggest that iron regulation and ISC homeostasis represent a possible target for aging interventions.

Next-Generation Genotoxicology: Using Modern Sequencing Technologies to Assess Somatic Mutagenesis and Cancer Risk.

Mutations have a profound effect on human health, particularly through an increased risk of carcinogenesis and genetic disease. The strong correlation between mutagenesis and carcinogenesis has been a driving force behind genotoxicity research for more than 50 years. The stochastic and infrequent nature of mutagenesis makes it challenging to observe and to study. Indeed, decades have been spent developing increasingly sophisticated assays and methods to study these low-frequency genetic errors, in hopes of better predicting which chemicals may be carcinogens, understanding their mode of action, and informing guidelines to prevent undue human exposure. While effective, widely used genetic selection-based technologies have a number of limitations that have hampered major advancements in the field of genotoxicity. Emerging new tools, in the form of enhanced next-generation sequencing platforms and methods, are changing this paradigm. In this review, we discuss rapidly evolving sequencing tools and technologies, such as error-corrected sequencing and single cell analysis, which we anticipate will fundamentally reshape the field. In addition, we consider a variety emerging applications for these new technologies, including the detection of DNA adducts, inference of mutational processes based on genomic site and local sequence contexts, and evaluation of genome engineering fidelity, as well as other cutting-edge challenges for the next 50 years of environmental and molecular mutagenesis research. Environ. Mol. Mutagen. 61:135-151, 2020. © 2019 The Authors. Environmental and Molecular Mutagenesis published by Wiley Periodicals, Inc. on behalf of Environmental Mutagen Society.

Cancer-Associated Mutations but No Cancer: Insights into the Early Steps of Carcinogenesis and Implications for Early Cancer Detection.

Cancer is a disease of aging fueled by the accumulation of somatic mutations. While mutations in tumors are well characterized, little is known about the early mutational processes that initiate tumorigenesis. Recent advances in next-generation sequencing (NGS) have enabled the detection of mutations in normal tissue, revealing an unanticipated high level of age-related somatic mutations affecting most individuals and tissues. Surprisingly, many of these mutations are similar to mutations commonly found in tumors, suggesting an ongoing process of positive selection and clonal expansion akin to what occurs in cancer, but within normal tissue. Here we discuss some of the most important biological and clinical implications of these novel findings, with a special focus on their impact for cancer detection and prediction.

Mitochondrially-targeted APOBEC1 is a potent mtDNA mutator affecting mitochondrial function and organismal fitness in Drosophila.

Somatic mutations in the mitochondrial genome (mtDNA) have been linked to multiple disease conditions and to ageing itself. In Drosophila, knock-in of a proofreading deficient mtDNA polymerase (POLG) generates high levels of somatic point mutations and also small indels, but surprisingly limited impact on organismal longevity or fitness. Here we describe a new mtDNA mutator model based on a mitochondrially-targeted cytidine deaminase, APOBEC1. mito-APOBEC1 acts as a potent mutagen which exclusively induces C:G>T:A transitions with no indels or mtDNA depletion. In these flies, the presence of multiple non-synonymous substitutions, even at modest heteroplasmy, disrupts mitochondrial function and dramatically impacts organismal fitness. A detailed analysis of the mutation profile in the POLG and mito-APOBEC1 models reveals that mutation type (quality) rather than quantity is a critical factor in impacting organismal fitness. The specificity for transition mutations and the severe phenotypes make mito-APOBEC1 an excellent mtDNA mutator model for ageing research.

Abnormal oxidative metabolism in a quiet genomic background underlies clear cell papillary renal cell carcinoma.

While genomic sequencing routinely identifies oncogenic alterations for the majority of cancers, many tumors harbor no discernable driver lesion. Here, we describe the exceptional molecular phenotype of a genomically quiet kidney tumor, clear cell papillary renal cell carcinoma (CCPAP). In spite of a largely wild-type nuclear genome, CCPAP tumors exhibit severe depletion of mitochondrial DNA (mtDNA) and RNA and high levels of oxidative stress, reflecting a shift away from respiratory metabolism. Moreover, CCPAP tumors exhibit a distinct metabolic phenotype uniquely characterized by accumulation of the sugar alcohol sorbitol. Immunohistochemical staining of primary CCPAP tumor specimens recapitulates both the depletion of mtDNA-encoded proteins and a lipid-depleted metabolic phenotype, suggesting that the cytoplasmic clarity in CCPAP is primarily related to the presence of glycogen. These results argue for non-genetic profiling as a tool for the study of cancers of unknown driver.

Defining the impact of mutation accumulation on replicative lifespan in yeast using cancer-associated mutator phenotypes.

Mutations accumulate within somatic cells and have been proposed to contribute to aging. It is unclear what level of mutation burden may be required to consistently reduce cellular lifespan. Human cancers driven by a mutator phenotype represent an intriguing model to test this hypothesis, since they carry the highest mutation burdens of any human cell. However, it remains technically challenging to measure the replicative lifespan of individual mammalian cells. Here, we modeled the consequences of cancer-related mutator phenotypes on lifespan using yeast defective for mismatch repair (MMR) and/or leading strand (Polε) or lagging strand (Polδ) DNA polymerase proofreading. Only haploid mutator cells with significant lifetime mutation accumulation (MA) exhibited shorter lifespans. Diploid strains, derived by mating haploids of various genotypes, carried variable numbers of fixed mutations and a range of mutator phenotypes. Some diploid strains with fewer than two mutations per megabase displayed a 25% decrease in lifespan, suggesting that moderate numbers of random heterozygous mutations can increase mortality rate. As mutation rates and burdens climbed, lifespan steadily eroded. Strong diploid mutator phenotypes produced a form of genetic anticipation with regard to aging, where the longer a lineage persisted, the shorter lived cells became. Using MA lines, we established a relationship between mutation burden and lifespan, as well as population doubling time. Our observations define a threshold of random mutation burden that consistently decreases cellular longevity in diploid yeast cells. Many human cancers carry comparable mutation burdens, suggesting that while cancers appear immortal, individual cancer cells may suffer diminished lifespan due to accrued mutation burden.

Mitochondrial DNA Mutations are Associated with Ulcerative Colitis Preneoplasia but Tend to be Negatively Selected in Cancer.

The role of mitochondrial DNA (mtDNA) mutations in cancer remains controversial. Ulcerative colitis is an inflammatory bowel disease that increases the risk of colorectal cancer and involves mitochondrial dysfunction, making it an ideal model to study the role of mtDNA in tumorigenesis. Our goal was to comprehensively characterize mtDNA mutations in ulcerative colitis tumorigenesis using Duplex Sequencing, an ultra-accurate next-generation sequencing method. We analyzed 46 colon biopsies from non-ulcerative colitis control patients and ulcerative colitis patients with and without cancer, including biopsies at all stages of dysplastic progression. mtDNA was sequenced at a median depth of 1,364x. Mutations were classified by mutant allele frequency: clonal > 0.95, subclonal 0.01-0.95, and very low frequency (VLF) < 0.01. We identified 208 clonal and subclonal mutations and 56,764 VLF mutations. Mutations were randomly distributed across the mitochondrial genome. Clonal and subclonal mutations increased in number and pathogenicity in early dysplasia, but decreased in number and pathogenicity in cancer. Most clonal, subclonal, and VLF mutations were C>T transitions in the heavy strand of mtDNA, which likely arise from DNA replication errors. A subset of VLF mutations were C>A transversions, which are probably due to oxidative damage. VLF transitions and indels were less abundant in the non-D-loop region and decreased with progression. Our results indicate that mtDNA mutations are frequent in ulcerative colitis preneoplasia but negatively selected in cancers. IMPLICATIONS: While mtDNA mutations might contribute to early ulcerative colitis tumorigenesis, they appear to be selected against in cancer, suggesting that functional mitochondria might be required for malignant transformation in ulcerative colitis.

Deleterious mitochondrial DNA point mutations are overrepresented in Drosophila expressing a proofreading-defective DNA polymerase γ.

Mitochondrial DNA (mtDNA) mutations cause severe maternally inherited syndromes and the accumulation of somatic mtDNA mutations is implicated in aging and common diseases. However, the mechanisms that influence the frequency and pathogenicity of mtDNA mutations are poorly understood. To address this matter, we created a Drosophila mtDNA mutator strain expressing a proofreading-deficient form of the mitochondrial DNA polymerase. Mutator flies have a dramatically increased somatic mtDNA mutation frequency that correlates with the dosage of the proofreading-deficient polymerase. Mutator flies also exhibit mitochondrial dysfunction, shortened lifespan, a progressive locomotor deficit, and loss of dopaminergic neurons. Surprisingly, the frequency of nonsynonymous, pathogenic, and conserved-site mutations in mutator flies exceeded predictions of a neutral mutational model, indicating the existence of a positive selection mechanism that favors deleterious mtDNA variants. We propose from these findings that deleterious mtDNA mutations are overrepresented because they selectively evade quality control surveillance or because they are amplified through compensatory mitochondrial biogenesis.