Reanalysis of mtDNA mutations of human primordial germ cells (PGCs) reveals NUMT contamination and suggests that selection in PGCs may be positive.

The resilience of the mitochondrial genome (mtDNA) to a high mutational pressure depends, in part, on negative purifying selection in the germline. A paradigm in the field has been that such selection, at least in part, takes place in primordial germ cells (PGCs). Specifically, Floros et al. (Nature Cell Biology 20: 144-51) reported an increase in the synonymity of mtDNA mutations (a sign of purifying selection) between early-stage and late-stage PGCs. We re-analyzed Floros' et al. data and determined that their mutational dataset was significantly contaminated with single nucleotide variants (SNVs) derived from a nuclear sequence of mtDNA origin (NUMT) located on chromosome 5. Contamination was caused by co-amplification of the NUMT sequence by cross-specific PCR primers. Importantly, when we removed NUMT-derived SNVs, the evidence of purifying selection was abolished. In addition to bulk PGCs, Floros et al. reported the analysis of single-cell late-stage PGCs, which were amplified with different sets of PCR primers that cannot amplify the NUMT sequence. Accordingly, there were no NUMT-derived SNVs among single PGC mutations. Interestingly, single PGC mutations show adecreaseof synonymity with increased intracellular mutant fraction. More specifically, nonsynonymous mutations show faster intracellular genetic drift towards higher mutant fraction than synonymous ones. This pattern is incompatible with predominantly negative selection. This suggests that germline selection of mtDNA mutations is a complex phenomenon and that the part of this process that takes place in PGCs may be predominantly positive. However counterintuitive, positive germline selection of detrimental mtDNA mutations has been reported previously andpotentially may be evolutionarily advantageous.

Interpreting Sequence-Levenshtein distance for determining error type and frequency between two embedded sequences of equal length.

Levenshtein distance is a commonly used edit distance metric, typically applied in language processing, and to a lesser extent, in molecular biology analysis. Biological nucleic acid sequences are often embedded in longer sequences and are subject to insertion and deletion errors that introduce frameshift during sequencing. These frameshift errors are due to string context and should not be counted as true biological errors. Sequence-Levenshtein distance is a modification to Levenshtein distance that is permissive of frameshift error without additional penalty. However, in a biological context Levenshtein distance needs to accommodate both frameshift and weighted errors, which Sequence-Levenshtein distance cannot do. Errors are weighted when they are associated with a numerical cost that corresponds to their frequency of appearance. Here, we describe a modification that allows the use of Levenshtein distance and Sequence-Levenshtein distance to appropriately accommodate penalty-free frameshift between embedded sequences and correctly weight specific error types.

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.

Are some mutations more equal than others?

A large-scale study of mutations in mitochondrial DNA has revealed a subset that do not accumulate with age.

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.

Dynamics of the most common pathogenic mtDNA variant m.3243A > G demonstrate frequency-dependency in blood and positive selection in the germline.

The A-to-G point mutation at position 3243 in the human mitochondrial genome (m.3243A > G) is the most common pathogenic mtDNA variant responsible for disease in humans. It is widely accepted that m.3243A > G levels decrease in blood with age, and an age correction representing ~ 2% annual decline is often applied to account for this change in mutation level. Here we report that recent data indicate that the dynamics of m.3243A > G are more complex and depend on the mutation level in blood in a bi-phasic way. Consequently, the traditional 2% correction, which is adequate 'on average', creates opposite predictive biases at high and low mutation levels. Unbiased age correction is needed to circumvent these drawbacks of the standard model. We propose to eliminate both biases by using an approach where age correction depends on mutation level in a biphasic way to account for the dynamics of m.3243A > G in blood. The utility of this approach was further tested in estimating germline selection of m.3243A > G. The biphasic approach permitted us to uncover patterns consistent with the possibility of positive selection for m.3243A > G. Germline selection of m.3243A > G shows an 'arching' profile by which selection is positive at intermediate mutant fractions and declines at high and low mutant fractions. We conclude that use of this biphasic approach will greatly improve the accuracy of modelling changes in mtDNA mutation frequencies in the germline and in somatic cells during aging.

Mitochondrial Pseudogenes Suggest Repeated Inter-Species Hybridization among Direct Human Ancestors.

The hypothesis that the evolution of humans involves hybridization between diverged species has been actively debated in recent years. We present the following novel evidence in support of this hypothesis: the analysis of nuclear pseudogenes of mtDNA ("NUMTs"). NUMTs are considered "mtDNA fossils" as they preserve sequences of ancient mtDNA and thus carry unique information about ancestral populations. Our comparison of a NUMT sequence shared by humans, chimpanzees, and gorillas with their mtDNAs implies that, around the time of divergence between humans and chimpanzees, our evolutionary history involved the interbreeding of individuals whose mtDNA had diverged as much as ~4.5 Myr prior. This large divergence suggests a distant interspecies hybridization. Additionally, analysis of two other NUMTs suggests that such events occur repeatedly. Our findings suggest a complex pattern of speciation in primate/human ancestors and provide one potential explanation for the mosaic nature of fossil morphology found at the emergence of the hominin lineage. A preliminary version of this manuscript was uploaded to the preprint server BioRxiv in 2017 (10.1101/134502).

3GOLD: optimized Levenshtein distance for clustering third-generation sequencing data.

BACKGROUND: Third-generation sequencing offers some advantages over next-generation sequencing predecessors, but with the caveat of harboring a much higher error rate. Clustering-related sequences is an essential task in modern biology. To accurately cluster sequences rich in errors, error type and frequency need to be accounted for. Levenshtein distance is a well-established mathematical algorithm for measuring the edit distance between words and can specifically weight insertions, deletions and substitutions. However, there are drawbacks to using Levenshtein distance in a biological context and hence has rarely been used for this purpose. We present novel modifications to the Levenshtein distance algorithm to optimize it for clustering error-rich biological sequencing data. RESULTS: We successfully introduced a bidirectional frameshift allowance with end-user determined accommodation caps combined with weighted error discrimination. Furthermore, our modifications dramatically improved the computational speed of Levenstein distance. For simulated ONT MinION and PacBio Sequel datasets, the average clustering sensitivity for 3GOLD was 41.45% (S.D. 10.39) higher than Sequence-Levenstein distance, 52.14% (S.D. 9.43) higher than Levenshtein distance, 55.93% (S.D. 8.67) higher than Starcode, 42.68% (S.D. 8.09) higher than CD-HIT-EST and 61.49% (S.D. 7.81) higher than DNACLUST. For biological ONT MinION data, 3GOLD clustering sensitivity was 27.99% higher than Sequence-Levenstein distance, 52.76% higher than Levenshtein distance, 56.39% higher than Starcode, 48% higher than CD-HIT-EST and 70.4% higher than DNACLUST. CONCLUSION: Our modifications to Levenshtein distance have improved its speed and accuracy compared to the classic Levenshtein distance, Sequence-Levenshtein distance and other commonly used clustering approaches on simulated and biological third-generation sequenced datasets. Our clustering approach is appropriate for datasets of unknown cluster centroids, such as those generated with unique molecular identifiers as well as known centroids such as barcoded datasets. A strength of our approach is high accuracy in resolving small clusters and mitigating the number of singletons.

LUCS: a high-resolution nucleic acid sequencing tool for accurate long-read analysis of individual DNA molecules.

Nucleic acid sequence analyses are fundamental to all aspects of biological research, spanning aging, mitochondrial DNA (mtDNA) and cancer, as well as microbial and viral evolution. Over the past several years, significant improvements in DNA sequencing, including consensus sequence analysis, have proven invaluable for high-throughput studies. However, all current DNA sequencing platforms have limited utility for studies of complex mixtures or of individual long molecules, the latter of which is crucial to understanding evolution and consequences of single nucleotide variants and their combinations. Here we report a new technology termed LUCS (Long-molecule UMI-driven Consensus Sequencing), in which reads from third-generation sequencing are aggregated by unique molecular identifiers (UMIs) specific for each individual DNA molecule. This enables in-silico reconstruction of highly accurate consensus reads of each DNA molecule independent of other molecules in the sample. Additionally, use of two UMIs enables detection of artificial recombinants (chimeras). As proof of concept, we show that application of LUCS to assessment of mitochondrial genomes in complex mixtures from single cells was associated with an error rate of 1X10-4 errors/nucleotide. Thus, LUCS represents a major step forward in DNA sequencing that offers high-throughput capacity and high-accuracy reads in studies of long DNA templates and nucleotide variants in heterogenous samples.

A nanoscale, multi-parametric flow cytometry-based platform to study mitochondrial heterogeneity and mitochondrial DNA dynamics.

Mitochondria are well-characterized regarding their function in both energy production and regulation of cell death; however, the heterogeneity that exists within mitochondrial populations is poorly understood. Typically analyzed as pooled samples comprised of millions of individual mitochondria, there is little information regarding potentially different functionality across subpopulations of mitochondria. Herein we present a new methodology to analyze mitochondria as individual components of a complex and heterogeneous network, using a nanoscale and multi-parametric flow cytometry-based platform. We validate the platform using multiple downstream assays, including electron microscopy, ATP generation, quantitative mass-spectrometry proteomic profiling, and mtDNA analysis at the level of single organelles. These strategies allow robust analysis and isolation of mitochondrial subpopulations to more broadly elucidate the underlying complexities of mitochondria as these organelles function collectively within a cell.

Impact of exercise on oocyte quality in the POLG mitochondrial DNA mutator mouse.

The mtDNA 'mutator' mouse, also called the 'POLG' mouse, is a well-characterized model frequently used for studies of progeroid aging. Harboring a mutation in the proofreading domain of the mitochondrial polymerase, polymerase-γ (Polg), POLG mice acquire mtDNA mutations at an accelerated rate. This results in premature mitochondrial dysfunction and a systemic aging phenotype. Previous work has demonstrated that the progeroid phenotype in POLG is attenuated following endurance exercise, the only reported intervention to extend health span and lifespan of these mice. Herein, oocyte quality was evaluated in sedentary and exercised POLG mice. In mice homozygous for the Polg mutation, litter size is dramatically reduced as compared to heterozygous Polg mice. Following ovarian hyper-stimulation, oocytes were retrieved until 9 months of age in exercised and sedentary groups, with no oocytes ovulated thereafter. Although ovulated oocyte numbers were not impacted by exercise, we did find a modest improvement in both the ovarian follicle reserve and in oocyte quality based on meiotic spindle assembly, chromosomal segregation and mitochondrial distribution at 7 months of age in exercised POLG mice as compared to sedentary counterparts. Of note, analysis of mtDNA mutational load revealed no differences between exercised and sedentary groups. Collectively, these data indicate that exercise differentially influences somatic tissues of the POLG mouse as compared to oocytes, highlighting important mechanistic differences between mitochondrial regulatory mechanisms in the soma and the germline.

Influence of Maternal Aging on Mitochondrial Heterogeneity, Inheritance, and Function in Oocytes and Preimplantation Embryos.

Contrasting the equal contribution of nuclear genetic material from maternal and paternal sources to offspring, passage of mitochondria, and thus mitochondrial DNA (mtDNA), is uniparental through the egg. Since mitochondria in eggs are ancestral to all somatic mitochondria of the next generation and to all cells of future generations, oocytes must prepare for the high energetic demands of maturation, fertilization and embryogenesis while simultaneously ensuring that their mitochondrial genomes are inherited in an undamaged state. Although significant effort has been made to understand how the mtDNA bottleneck and purifying selection act coordinately to prevent silent and unchecked spreading of invisible mtDNA mutations through the female germ line across successive generations, it is unknown if and how somatic cells of the immediate next generation are spared from inheritance of detrimental mtDNA molecules. Here, we review unique aspects of mitochondrial activity and segregation in eggs and early embryos, and how these events play into embryonic developmental competency in the face of advancing maternal age.

Data on the time of integration of the human mitochondrial pseudogenes (NUMTs) into the nuclear genome.

The data and methods presented in this article are supplementing the research article "Integration of mtDNA pseudogenes into the nuclear genome coincides with speciation of the human genus. A hypothesis", DOI: 10.1016/j.mito.2016.12.001 (Gunbin et al., 2017) [1]. Mitochondrial DNA is known to get inserted into nuclear DNA to form NUMTs, i.e. nuclear DNA pseudogenes of the mtDNA. We present here the sequences of selected NUMTs, in which time of integration can be determined with sufficient precision. We report their chromosomal positions , their position within the great ape mtDNA phylogeny, and their times of integration into the nuclear genome. The methods used to generate the data and to control their quality are also presented. The dataset is made publicly available to enable critical or extended analyzes.

Integration of mtDNA pseudogenes into the nuclear genome coincides with speciation of the human genus. A hypothesis.

Fragments of mitochondrial DNA are known to get inserted into nuclear DNA to form NUMTs, i.e. nuclear pseudogenes of the mtDNA. The insertion of a NUMT is a rare event. Hundreds of pseudogenes have been cataloged in the human genome. NUMTs are, in essence, a special type of mutation with their own internal timer, which is synchronized with an established molecular clock, the mtDNA. Thus insertion of NUMTs can be timed with respect to evolution milestones such as the emergence of new species. We asked whether NUMTs were inserted uniformly over time or preferentially during certain periods of evolution, as implied by the "punctuated evolution" model. To our surprise, the NUMT insertion times do appear nonrandom with at least one cluster positioned at around 2.8 million years ago (Ma). Interestingly, 2.8Ma closely corresponds to the time of emergence of the genus Homo, and to a well-documented period of major climate change ca. 2.9-2.5Ma. It is tempting to hypothesize that the insertion of NUMTs is related to the speciation process. NUMTs could be either "riders", i.e., their insertion could be facilitated by the overall higher genome rearrangement activity during speciation, or "drivers", i.e. they may more readily get fixed in the population due to positive selection associated with speciation. If correct, the hypothesis would support the idea that evolution of our genus may have happened in a rapid, punctuated manner.

Amelioration of premature aging in mtDNA mutator mouse by exercise: the interplay of oxidative stress, PGC-1α, p53, and DNA damage. A hypothesis.

The mtDNA mutator mouse lacks the proofreading capacity of the sole mtDNA polymerase, leading to accumulation of somatic mtDNA mutations, and a profound premature aging phenotype including elevated oxidative stress and apoptosis, and reduced mitochondrial function. We have previously reported that endurance exercise alleviates the aging phenotype in the mutator mice, reduces oxidative stress, and enhances mitochondrial biogenesis. Here we summarize our findings, with the emphasis on the central role of p53 in these adaptations. We demonstrate that mtDNA in sedentary and exercised PolG mice carry similar amounts of mutations in muscle, but in addition to that sedentary mice have more non-mutational damage, which is mitigated by exercise. It follows therefore that the profound alleviation of the mtDNA mutator phenotype in muscle by exercise may not require a reduction in mtDNA mutational load, but rather a decrease of mtDNA damage and/or oxidative stress. We further hypothesize that the observed 'alleviation without a reduction of mutational load' implies that the oxidative stress in PolG muscle is maintained, at least in part, by the 'malicious cycle', a hypothetical positive feedback potentially driven by the 'transcriptional mutagenesis', that is the conversion of chemically modified nucleotides into mutant RNA bases by the mitochondrial RNA polymerase.

Exercise-induced mitochondrial p53 repairs mtDNA mutations in mutator mice.

BACKGROUND: Human genetic disorders and transgenic mouse models have shown that mitochondrial DNA (mtDNA) mutations and telomere dysfunction instigate the aging process. Epidemiologically, exercise is associated with greater life expectancy and reduced risk of chronic diseases. While the beneficial effects of exercise are well established, the molecular mechanisms instigating these observations remain unclear. RESULTS: Endurance exercise reduces mtDNA mutation burden, alleviates multisystem pathology, and increases lifespan of the mutator mice, with proofreading deficient mitochondrial polymerase gamma (POLG1). We report evidence for a POLG1-independent mtDNA repair pathway mediated by exercise, a surprising notion as POLG1 is canonically considered to be the sole mtDNA repair enzyme. Here, we show that the tumor suppressor protein p53 translocates to mitochondria and facilitates mtDNA mutation repair and mitochondrial biogenesis in response to endurance exercise. Indeed, in mutator mice with muscle-specific deletion of p53, exercise failed to prevent mtDNA mutations, induce mitochondrial biogenesis, preserve mitochondrial morphology, reverse sarcopenia, or mitigate premature mortality. CONCLUSIONS: Our data establish a new role for p53 in exercise-mediated maintenance of the mtDNA genome and present mitochondrially targeted p53 as a novel therapeutic modality for diseases of mitochondrial etiology.

Quantitation of Mitochondrial DNA Deletions Via Restriction Digestion/Long-Range Single-Molecule PCR.

Quantification of deletions in mtDNA is a long-standing problem in mutational analysis. We describe here an approach that combines the power of single-molecule PCR of the entire mitochondrial genome with the enrichment of the deletions by restriction digestion. This approach is indispensable if information about wide range of deletion types in a sample is critical, such as in studies concerning distribution of deletion breakpoints (as opposed to approaches where fraction of a single deletion or a limited set of deletions is used as a proxy for total deletion load). Because deletions in a sample are quantified almost exhaustively, the other important application of this approach involves studies where only small amounts of tissue, such as biopsies, are available.