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

Mitochondrial DNA mosaicism in normal human somatic cells

Somatic cells accumulate genomic alterations with age; however, our understanding of mitochondrial DNA (mtDNA) mosaicism remains limited. Here we investigated the genomes of 2,096 clones derived from three cell types across 31 donors, identifying 6,451 mtDNA variants with heteroplasmy levels of ≳0.3%. While the majority of these variants were unique to individual clones, suggesting stochastic acquisition with age, 409 variants (6%) were shared across multiple embryonic lineages, indicating their origin from heteroplasmy in fertilized eggs. The mutational spectrum exhibited replication-strand bias, implicating mtDNA replication as a major mutational process. We evaluated the mtDNA mutation rate (5.0 × 10-8 per base pair) and a turnover frequency of 10-20 per year, which are fundamental components shaping the landscape of mtDNA mosaicism over a lifetime. The expansion of mtDNA-truncating mutations toward homoplasmy was substantially suppressed. Our findings provide comprehensive insights into the origins, dynamics and functional consequences of mtDNA mosaicism in human somatic cells.

The rate and nature of mitochondrial DNA mutations in human pedigrees

We examined the rate and nature of mitochondrial DNA (mtDNA) mutations in humans using sequence data from 64,806 contemporary Icelanders from 2,548 matrilines. Based on 116,663 mother-child transmissions, 8,199 mutations were detected, providing robust rate estimates by nucleotide type, functional impact, position, and different alleles at the same position. We thoroughly document the true extent of hypermutability in mtDNA, mainly affecting the control region but also some coding-region variants. The results reveal the impact of negative selection on viable deleterious mutations, including rapidly mutating disease-associated 3243A>G and 1555A>G and pre-natal selection that most likely occurs during the development of oocytes. Finally, we show that the fate of new mutations is determined by a drastic germline bottleneck, amounting to an average of 3 mtDNA units effectively transmitted from mother to child.

NeMu: a comprehensive pipeline for accurate reconstruction of neutral mutation spectra from evolutionary data

The recognized importance of mutational spectra in molecular evolution is yet to be fully exploited beyond human cancer studies and model organisms. The wealth of intraspecific polymorphism data in the GenBank repository, covering a broad spectrum of genes and species, presents an untapped opportunity for detailed mutational spectrum analysis. Existing methods fall short by ignoring intermediate substitutions on the inner branches of phylogenetic trees and lacking the capability for cross-species mutational comparisons. To address these challenges, we present the NeMu pipeline, available at https://nemu-pipeline.com, a tool grounded in phylogenetic principles designed to provide comprehensive and scalable analysis of mutational spectra. Utilizing extensive sequence data from numerous available genome projects, NeMu rapidly and accurately reconstructs the neutral mutational spectrum. This tool, facilitating the reconstruction of gene- and species-specific mutational spectra, contributes to a deeper understanding of evolutionary mechanisms across the broad spectrum of known species.

Disruption of recombination machinery alters the mutational landscape in plant organellar genomes

Land plant organellar genomes have extremely low rates of point mutation yet also experience high rates of recombination and genome instability. Characterizing the molecular machinery responsible for these patterns is critical for understanding the evolution of these genomes. While much progress has been made towards understanding recombination activity in land plant organellar genomes, the relationship between recombination pathways and point mutation rates remains uncertain. The organellar targeted mutS homolog MSH1 has previously been shown to suppress point mutations as well as non-allelic recombination between short repeats in Arabidopsis thaliana. We therefore implemented high-fidelity Duplex Sequencing to test if other genes that function in recombination and maintenance of genome stability also affect point mutation rates. We found small to moderate increases in the frequency of single nucleotide variants (SNVs) and indels in mitochondrial and/or plastid genomes of A. thaliana mutant lines lacking radA, recA1, or recA3. In contrast, osb2 and why2 mutants did not exhibit an increase in point mutations compared to wild type (WT) controls. In addition, we analyzed the distribution of SNVs in previously generated Duplex Sequencing data from A. thaliana organellar genomes and found unexpected strand asymmetries and large effects of flanking nucleotides on mutation rates in WT plants and msh1 mutants. Finally, using long-read Oxford Nanopore sequencing, we characterized structural variants in organellar genomes of the mutant lines and show that different short repeat sequences become recombinationally active in different mutant backgrounds. Together, these complementary sequencing approaches shed light on how recombination may impact the extraordinarily low point mutation rates in plant organellar genomes.

DNA mismatch and damage patterns revealed by single-molecule sequencing

Mutations accumulate in the genome of every cell of the body throughout life, causing cancer and other diseases1,2. Most mutations begin as nucleotide mismatches or damage in one of the two strands of the DNA before becoming double-strand mutations if unrepaired or misrepaired3,4. However, current DNA-sequencing technologies cannot accurately resolve these initial single-strand events. Here we develop a single-molecule, long-read sequencing method (Hairpin Duplex Enhanced Fidelity sequencing (HiDEF-seq)) that achieves single-molecule fidelity for base substitutions when present in either one or both DNA strands. HiDEF-seq also detects cytosine deamination-a common type of DNA damage-with single-molecule fidelity. We profiled 134 samples from diverse tissues, including from individuals with cancer predisposition syndromes, and derive from them single-strand mismatch and damage signatures. We find correspondences between these single-strand signatures and known double-strand mutational signatures, which resolves the identity of the initiating lesions. Tumours deficient in both mismatch repair and replicative polymerase proofreading show distinct single-strand mismatch patterns compared to samples that are deficient in only polymerase proofreading. We also define a single-strand damage signature for APOBEC3A. In the mitochondrial genome, our findings support a mutagenic mechanism occurring primarily during replication. As double-strand DNA mutations are only the end point of the mutation process, our approach to detect the initiating single-strand events at single-molecule resolution will enable studies of how mutations arise in a variety of contexts, especially in cancer and ageing.

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.

Mitochondrial Targeted Interventions for Aging

Changes in mitochondrial function play a critical role in the basic biology of aging and age-related disease. Mitochondria are typically thought of in the context of ATP production and oxidant production. However, it is clear that the mitochondria sit at a nexus of cell signaling where they affect metabolite, redox, and energy status, which influence many factors that contribute to the biology of aging, including stress responses, proteostasis, epigenetics, and inflammation. This has led to growing interest in identifying mitochondrial targeted interventions to delay or reverse age-related decline in function and promote healthy aging. In this review, we discuss the diverse roles of mitochondria in the cell. We then highlight some of the most promising strategies and compounds to target aging mitochondria in preclinical testing. Finally, we review the strategies and compounds that have advanced to clinical trials to test their ability to improve health in older adults.

Investigating low frequency somatic mutations in Arabidopsis with Duplex Sequencing

Mutations are the source of novel genetic diversity but can also lead to disease and maladaptation. The conventional view is that mutations occur randomly with respect to their environment-specific fitness consequences. However, intragenomic mutation rates can vary dramatically due to transcription coupled repair and based on local epigenomic modifications, which are non-uniformly distributed across genomes. One sequence feature associated with decreased mutation is higher expression level, which can vary depending on environmental cues. To understand whether the association between expression level and mutation rate creates a systematic relationship with environment-specific fitness effects, we perturbed expression through a heat treatment in Arabidopsis thaliana. We quantified gene expression to identify differentially expressed genes, which we then targeted for mutation detection using Duplex Sequencing. This approach provided a highly accurate measurement of the frequency of rare somatic mutations in vegetative plant tissues, which has been a recent source of uncertainty in plant mutation research. We included mutant lines lacking mismatch repair (MMR) and base excision repair (BER) capabilities to understand how repair mechanisms may drive biased mutation accumulation. We found wild type (WT) and BER mutant mutation frequencies to be very low (mean variant frequency 1.8×10-8 and 2.6×10-8, respectively), while MMR mutant frequencies were significantly elevated (1.13×10-6). These results show that somatic variant frequencies are extremely low in WT plants, indicating that larger datasets will be needed to address the fundamental evolutionary question as to whether environmental change leads to gene-specific changes in mutation rate.

The Fate of Oxidative Strand Breaks in Mitochondrial DNA

Mitochondrial DNA (mtDNA) is particularly vulnerable to somatic mutagenesis. Potential mechanisms include DNA polymerase γ (POLG) errors and the effects of mutagens, such as reactive oxygen species. Here, we studied the effects of transient hydrogen peroxide (H₂O₂ pulse) on mtDNA integrity in cultured HEK 293 cells, applying Southern blotting, ultra-deep short-read and long-read sequencing. In wild-type cells, 30 min after the H₂O₂ pulse, linear mtDNA fragments appear, representing double-strand breaks (DSB) with ends characterized by short GC stretches. Intact supercoiled mtDNA species reappear within 2-6 h after treatment and are almost completely recovered after 24 h. BrdU incorporation is lower in H₂O₂-treated cells compared to non-treated cells, suggesting that fast recovery is not associated with mtDNA replication, but is driven by rapid repair of single-strand breaks (SSBs) and degradation of DSB-generated linear fragments. Genetic inactivation of mtDNA degradation in exonuclease deficient POLG p.D274A mutant cells results in the persistence of linear mtDNA fragments with no impact on the repair of SSBs. In conclusion, our data highlight the interplay between the rapid processes of SSB repair and DSB degradation and the much slower mtDNA re-synthesis after oxidative damage, which has important implications for mtDNA quality control and the potential generation of somatic mtDNA deletions.

Secondary structure of the human mitochondrial genome affects formation of deletions

BACKGROUND

Aging in postmitotic tissues is associated with clonal expansion of somatic mitochondrial deletions, the origin of which is not well understood. Such deletions are often flanked by direct nucleotide repeats, but this alone does not fully explain their distribution. Here, we hypothesized that the close proximity of direct repeats on single-stranded mitochondrial DNA (mtDNA) might play a role in the formation of deletions.

RESULTS

By analyzing human mtDNA deletions in the major arc of mtDNA, which is single-stranded during replication and is characterized by a high number of deletions, we found a non-uniform distribution with a "hot spot" where one deletion breakpoint occurred within the region of 6-9 kb and another within 13-16 kb of the mtDNA. This distribution was not explained by the presence of direct repeats, suggesting that other factors, such as the spatial proximity of these two regions, can be the cause. In silico analyses revealed that the single-stranded major arc may be organized as a large-scale hairpin-like loop with a center close to 11 kb and contacting regions between 6-9 kb and 13-16 kb, which would explain the high deletion activity in this contact zone. The direct repeats located within the contact zone, such as the well-known common repeat with a first arm at 8470-8482 bp (base pair) and a second arm at 13,447-13,459 bp, are three times more likely to cause deletions compared to direct repeats located outside of the contact zone. A comparison of age- and disease-associated deletions demonstrated that the contact zone plays a crucial role in explaining the age-associated deletions, emphasizing its importance in the rate of healthy aging.

CONCLUSIONS

Overall, we provide topological insights into the mechanism of age-associated deletion formation in human mtDNA, which could be used to predict somatic deletion burden and maximum lifespan in different human haplogroups and mammalian species.

Single-strand mismatch and damage patterns revealed by single-molecule DNA sequencing

Mutations accumulate in the genome of every cell of the body throughout life, causing cancer and other genetic diseases1-4. Almost all of these mosaic mutations begin as nucleotide mismatches or damage in only one of the two strands of the DNA prior to becoming double-strand mutations if unrepaired or misrepaired5. However, current DNA sequencing technologies cannot resolve these initial single-strand events. Here, we developed a single-molecule, long-read sequencing method that achieves single-molecule fidelity for single-base substitutions when present in either one or both strands of the DNA. It also detects single-strand cytosine deamination events, a common type of DNA damage. We profiled 110 samples from diverse tissues, including from individuals with cancer-predisposition syndromes, and define the first single-strand mismatch and damage signatures. We find correspondences between these single-strand signatures and known double-strand mutational signatures, which resolves the identity of the initiating lesions. Tumors deficient in both mismatch repair and replicative polymerase proofreading show distinct single-strand mismatch patterns compared to samples deficient in only polymerase proofreading. In the mitochondrial genome, our findings support a mutagenic mechanism occurring primarily during replication. Since the double-strand DNA mutations interrogated by prior studies are only the endpoint of the mutation process, our approach to detect the initiating single-strand events at single-molecule resolution will enable new studies of how mutations arise in a variety of contexts, especially in cancer and aging.

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.

mtUPR Modulation as a Therapeutic Target for Primary and Secondary Mitochondrial Diseases

Mitochondrial dysfunction is a key pathological event in many diseases. Its role in energy production, calcium homeostasis, apoptosis regulation, and reactive oxygen species (ROS) balance render mitochondria essential for cell survival and fitness. However, there are no effective treatments for most primary and secondary mitochondrial diseases to this day. Therefore, new therapeutic approaches, such as the modulation of the mitochondrial unfolded protein response (mtUPR), are being explored. mtUPRs englobe several compensatory processes related to proteostasis and antioxidant system mechanisms. mtUPR activation, through an overcompensation for mild intracellular stress, promotes cell homeostasis and improves lifespan and disease alterations in biological models of mitochondrial dysfunction in age-related diseases, cardiopathies, metabolic disorders, and primary mitochondrial diseases. Although mtUPR activation is a promising therapeutic option for many pathological conditions, its activation could promote tumor progression in cancer patients, and its overactivation could lead to non-desired side effects, such as the increased heteroplasmy of mitochondrial DNA mutations. In this review, we present the most recent data about mtUPR modulation as a therapeutic approach, its role in diseases, and its potential negative consequences in specific pathological situations.

Modulating mitochondrial DNA mutations: factors shaping heteroplasmy in the germ line and somatic cells

Until recently it was thought that most humans only harbor one type of mitochondrial DNA (mtDNA), however, deep sequencing and single-cell analysis has shown the converse - that mixed populations of mtDNA (heteroplasmy) are the norm. This is important because heteroplasmy levels can change dramatically during transmission in the female germ line, leading to high levels causing severe mitochondrial diseases. There is also emerging evidence that low level mtDNA mutations contribute to common late onset diseases such as neurodegenerative disorders and cardiometabolic diseases because the inherited mutation levels can change within developing organs and non-dividing cells over time. Initial predictions suggested that the segregation of mtDNA heteroplasmy was largely stochastic, with an equal tendency for levels to increase or decrease. However, transgenic animal work and single-cell analysis have shown this not to be the case during germ-line transmission and in somatic tissues during life. Mutation levels in specific mtDNA regions can increase or decrease in different contexts and the underlying molecular mechanisms are starting to be unraveled. In this review we provide a synthesis of recent literature on the mechanisms of selection for and against mtDNA variants. We identify the most pertinent gaps in our understanding and suggest ways these could be addressed using state of the art techniques.

A mitochondria-specific mutational signature of aging: increased rate of A > G substitutions on the heavy strand

The mutational spectrum of the mitochondrial DNA (mtDNA) does not resemble any of the known mutational signatures of the nuclear genome and variation in mtDNA mutational spectra between different organisms is still incomprehensible. Since mitochondria are responsible for aerobic respiration, it is expected that mtDNA mutational spectrum is affected by oxidative damage. Assuming that oxidative damage increases with age, we analyse mtDNA mutagenesis of different species in regards to their generation length. Analysing, (i) dozens of thousands of somatic mtDNA mutations in samples of different ages (ii) 70053 polymorphic synonymous mtDNA substitutions reconstructed in 424 mammalian species with different generation lengths and (iii) synonymous nucleotide content of 650 complete mitochondrial genomes of mammalian species we observed that the frequency of AH > GH substitutions (H: heavy strand notation) is twice bigger in species with high versus low generation length making their mtDNA more AH poor and GH rich. Considering that AH > GH substitutions are also sensitive to the time spent single-stranded (TSSS) during asynchronous mtDNA replication we demonstrated that AH > GH substitution rate is a function of both species-specific generation length and position-specific TSSS. We propose that AH > GH is a mitochondria-specific signature of oxidative damage associated with both aging and TSSS.

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.

Mitochondrial DNA: Consensuses and Controversies

In the course of its short history, mitochondrial DNA (mtDNA) has made a long journey from obscurity to the forefront of research on major biological processes. mtDNA alterations have been found in all major disease groups, and their significance remains the subject of intense research. Despite remarkable progress, our understanding of the major aspects of mtDNA biology, such as its replication, damage, repair, transcription, maintenance, etc., is frustratingly limited. The path to better understanding mtDNA and its role in cells, however, remains torturous and not without errors, which sometimes leave a long trail of controversy behind them. This review aims to provide a brief summary of our current knowledge of mtDNA and highlight some of the controversies that require attention from the mitochondrial research community.

Mutational profiling of mtDNA control region reveals tumor-specific evolutionary selection involved in mitochondrial dysfunction

Background

Mitochondrial DNA (mtDNA) mutations alter mitochondrial function in oxidative metabolism and play an important role in tumorigenesis. A series of studies have demonstrated that the mtDNA control region (mtCTR), which is essential for mtDNA replication and transcription, represents a mutational hotspot in human tumors. However, a comprehensive pan-cancer evolutionary pattern analysis of mtCTR mutations is urgently needed.

Methods

We generated a comprehensive combined dataset containing 10026 mtDNA somatic mutations from 4664 patients, covering 20 tumor types based on public and private next-generation sequencing data.

Findings

Our results demonstrated a significantly higher and much more variable mutation rate in mtCTR than in the coding region across different tumor types. Moreover, our data showed a remarkable distributional bias of tumor somatic mutations between the hypervariable segment (HVS) and non-HVS, with a significantly higher mutation density and average mutation sites in HVS. Importantly, the tumor-specific mutational pattern between mtCTR HVS and non-HVS was identified, which was classified into three evolutionary selection types (relaxed, moderate, and strict constraint types). Analysis of substitution patterns revealed that the prevalence of C_H > T_H in non-HVS greatly contributed to the mutational selection pattern of mtCTR across different tumor types. Furthermore, we found that the mutational pattern of mtCTR in the four tumor types was clearly associated with mitochondrial biogenesis, mitochondrial oxidative metabolism, and the overall survival of patients.

Interpretation

Our results suggest that somatic mutations in mtCTR may be shaped by tumor-specific selective pressure and are involved in tumorigenesis.

Fundings

National Natural Science Foundation of China [grants 82020108023, 81830070, 81872302], and Autonomous Project of State Key Laboratory of Cancer Biology, China [grants CBSKL2019ZZ06, CBSKL2019ZZ27].

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.

The Mitochondrial Genome in Aging and Disease and the Future of Mitochondrial Therapeutics

Mitochondria are intracellular organelles that utilize nutrients to generate energy in the form of ATP by oxidative phosphorylation. Mitochondrial DNA (mtDNA) in humans is a 16,569 base pair double-stranded circular DNA that encodes for 13 vital proteins of the electron transport chain. Our understanding of the mitochondrial genome’s transcription, translation, and maintenance is still emerging, and human pathologies caused by mtDNA dysfunction are widely observed. Additionally, a correlation between declining mitochondrial DNA quality and copy number with organelle dysfunction in aging is well-documented in the literature. Despite tremendous advancements in nuclear gene-editing technologies and their value in translational avenues, our ability to edit mitochondrial DNA is still limited. In this review, we discuss the current therapeutic landscape in addressing the various pathologies that result from mtDNA mutations. We further evaluate existing gene therapy efforts, particularly allotopic expression and its potential to become an indispensable tool for restoring mitochondrial health in disease and aging.

Mitochondrial mutations in Caenorhabditis elegans show signatures of oxidative damage and an AT-bias

Rapid mutation rates are typical of mitochondrial genomes (mtDNAs) in animals, but it is not clear why. The difficulty of obtaining measurements of mtDNA mutation that are not biased by natural selection has stymied efforts to distinguish between competing hypotheses about the causes of high mtDNA mutation rates. Several studies which have measured mtDNA mutations in nematodes have yielded small datasets with conflicting conclusions about the relative abundance of different substitution classes (i.e., the mutation spectrum). We therefore leveraged Duplex Sequencing, a high-fidelity DNA sequencing technique, to characterize de novo mtDNA mutations in Caenorhabditis elegans. This approach detected nearly an order of magnitude more mtDNA mutations than documented in any previous nematode mutation study. Despite an existing extreme AT bias in the C. elegans mtDNA (75.6% AT), we found that a significant majority of mutations increase genomic AT content. Compared to some prior studies in nematodes and other animals, the mutation spectrum reported here contains an abundance of CG→AT transversions, supporting the hypothesis that oxidative damage may be a driver of mtDNA mutations in nematodes. Furthermore, we found an excess of G→T and C→T changes on the coding DNA strand relative to the template strand, consistent with increased exposure to oxidative damage. Analysis of the distribution of mutations across the mtDNA revealed significant variation among protein-coding genes and as well as among neighboring nucleotides. This high-resolution view of mitochondrial mutations in C. elegans highlights the value of this system for understanding relationships among oxidative damage, replication error, and mtDNA mutation.