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 A_H > G_H substitutions (H: heavy strand notation) is twice bigger in species with high versus low generation length making their mtDNA more A_H poor and G_H rich. Considering that A_H > G_H substitutions are also sensitive to the time spent single-stranded (TSSS) during asynchronous mtDNA replication we demonstrated that A_H > G_H substitution rate is a function of both species-specific generation length and position-specific TSSS. We propose that A_H > G_H is a mitochondria-specific signature of oxidative damage associated with both aging and TSSS.

Mitochondrial DNA (mtDNA) mutagenesis remains poorly understood despite its crucial role in disease, aging, and evolutionary tracing. In this study, we reconstructed a comprehensive 192-component mtDNA mutational spectrum for chordates by analyzing 118,397 synonymous mutations in the CytB gene across 1,697 species and five classes. This analysis revealed three primary forces shaping mtDNA mutagenesis: (i) symmetrical, replication-driven errors by mitochondrial polymerase (POLG), resulting in C > T and A > G mutations that are highly conserved across classes; (ii) asymmetrical, damage-driven C > T mutations on the single-stranded heavy strand with clock-like dynamics; and (iii) asymmetrical A > G mutations on the heavy strand, with dynamics suggesting sensitivity to oxidative damage. The third component, sensitive to oxidative damage, positions mtDNA mutagenesis as a promising marker for metabolic and physiological processes across various classes, species, organisms, tissues, and cells. The deconvolution of the mutational spectra into mutational signatures uncovered deficiencies in both base excision repair (BER) and mismatch repair (MMR) pathways. Further analysis of mutation hotspots, abasic sites, and mutational asymmetries underscores the critical role of single-stranded DNA damage (components ii and iii), which, uncorrected due to BER and MMR deficiencies, contributes roughly as many mutations as POLG-induced errors (component i).

To elucidate the primary factors shaping mitochondrial DNA (mtDNA) mutagenesis, we derived a comprehensive 192-component mtDNA mutational spectrum using 86,149 polymorphic synonymous mutations reconstructed from the CytB gene of 967 chordate species. The mtDNA spectrum analysis provided numerous findings on repair and mutation processes, breaking it down into three main signatures: (i) symmetrical, evenly distributed across both strands, mutations, induced by gamma DNA polymerase (about 50% of all mutations); (ii) asymmetrical, heavy-strand-specific, C>T mutations (about 30%); and (iii) asymmetrical, heavy-strand-specific A>G mutations, influenced by metabolic and age-specific factors (about 20%). We propose that both asymmetrical signatures are driven by single-strand specific damage coupled with inefficient base excision repair on the lagging (heavy) strand of mtDNA. Understanding the detailed mechanisms of this damage is crucial for developing strategies to reduce somatic mtDNA mutational load, which is vital for combating age-related diseases.

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.

Molybdenum cofactor deficiency type B (MOCODB; #252160) is an autosomal recessive metabolic disorder that has only been described in 37 affected patients. In this report, we describe the presence of an in-frame homozygous variant (c.471_477delTTTAAAAinsG) in the MOCS2 gene in an affected child, diagnosed with Ohtahara syndrome according to the clinical manifestations. The analysis of the three-dimensional structure of the protein and the amino acid substitutions suggested the pathogenicity of this mutation. To prevent transmitting this mutation to the next generation, we used preimplantation genetic testing for the monogenic disorders (PGT-M) protocol to select MOCS2 gene mutant-free embryos for transfer in an in vitro fertilization (IVF) program. As a result, a healthy child was born. Interestingly, both parents of the proband shared an identical mitochondrial (mt) DNA control region, assuming their close relationship and thus suggesting that both copies of the nuclear rare variant c.471_477delTTTAAAAinsG may have been transmitted from the same female ancestor. Our estimation of the a priori probability of meeting individuals with the same mtDNA haplotype confirms the assumption of a possible distant maternal relationship among the proband’s direct relatives.

In this work, we adapted a method for isolation of a highly purified fraction of mitochondrial DNA from muscle tissues suitable for further sample preparation of libraries without an amplification step for sequencing tasks on various NGS platforms and a method for evaluating the purity of DNA from contamination by nuclear genome regions. We optimized several techniques ^7,8 for enrichment of the mitochondrial fractions and purifying mtDNA. Here, we describe a protocol that allows getting from 80-100 mg of muscle tissues up to 1000 ng mtDNA, almost free from impurities of RNA and fragments of the nuclear genome.To assess the degree of purity of human mtDNA fraction from impurities of the nuclear genome, we adapted the PCR-screening technique ⁷ for the beta-actin gene region and AluSx-repeats in the human genome. This methodology avoids false-heteroplasmy calls (PCR biases or NUMT contamination) that occur when long-range PCR amplification is used for mtDNA enrichment.