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.

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.

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.

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.

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.

Clonal expansion of early to mid-life mitochondrial DNA point mutations drives mitochondrial dysfunction during human ageing.

Age-related decline in the integrity of mitochondria is an important contributor to the human ageing process. In a number of ageing stem cell populations, this decline in mitochondrial function is due to clonal expansion of individual mitochondrial DNA (mtDNA) point mutations within single cells. However the dynamics of this process and when these mtDNA mutations occur initially are poorly understood. Using human colorectal epithelium as an exemplar tissue with a well-defined stem cell population, we analysed samples from 207 healthy participants aged 17-78 years using a combination of techniques (Random Mutation Capture, Next Generation Sequencing and mitochondrial enzyme histochemistry), and show that: 1) non-pathogenic mtDNA mutations are present from early embryogenesis or may be transmitted through the germline, whereas pathogenic mtDNA mutations are detected in the somatic cells, providing evidence for purifying selection in humans, 2) pathogenic mtDNA mutations are present from early adulthood (<20 years of age), at both low levels and as clonal expansions, 3) low level mtDNA mutation frequency does not change significantly with age, suggesting that mtDNA mutation rate does not increase significantly with age, and 4) clonally expanded mtDNA mutations increase dramatically with age. These data confirm that clonal expansion of mtDNA mutations, some of which are generated very early in life, is the major driving force behind the mitochondrial dysfunction associated with ageing of the human colorectal epithelium.

In vivo generation of transplantable human hematopoietic cells from induced pluripotent stem cells.

Lineage-restricted cells can be reprogrammed to a pluripotent state known as induced pluripotent stem (iPS) cells through overexpression of 4 transcription factors. iPS cells are similar to human embryonic stem (hES) cells and have the same ability to generate all the cells of the human body, including blood cells. However, this process is extremely inefficient and to date has been unsuccessful at differentiating iPS into hematopoietic stem cells (HSCs). We hypothesized that iPS cells, injected into NOD.Cg-Prkdc(scid) Il2rg(tm1Wjl)/SzJ immunocompromised (NSG) mice could give rise to hematopoietic stem/progenitor cells (HSPCs) during teratoma formation. Here, we report a novel in vivo system in which human iPS cells differentiate within teratomas to derive functional myeloid and lymphoid cells. Similarly, HSPCs can be isolated from teratoma parenchyma and reconstitute a human immune system when transplanted into immunodeficient mice. Our data provide evidence that in vivo generation of patient customized cells is feasible, providing materials that could be useful for transplantation, human antibody generation, and drug screening applications.

Mitochondrial DNA deletions are abundant and cause functional impairment in aged human substantia nigra neurons.

Using a novel single-molecule PCR approach to quantify the total burden of mitochondrial DNA (mtDNA) molecules with deletions, we show that a high proportion of individual pigmented neurons in the aged human substantia nigra contain very high levels of mtDNA deletions. Molecules with deletions are largely clonal within each neuron; that is, they originate from a single deleted mtDNA molecule that has expanded clonally. The fraction of mtDNA deletions is significantly higher in cytochrome c oxidase (COX)-deficient neurons than in COX-positive neurons, suggesting that mtDNA deletions may be directly responsible for impaired cellular respiration.

Recombination of mitochondrial DNA in skeletal muscle of individuals with multiple mitochondrial DNA heteroplasmy.

Experimental evidence for human mitochondrial DNA (mtDNA) recombination was recently obtained in an individual with paternal inheritance of mtDNA and in an in vitro cell culture system. Whether mtDNA recombination is a common event in humans remained to be determined. To detect mtDNA recombination in human skeletal muscle, we analyzed the distribution of alleles in individuals with multiple mtDNA heteroplasmy using single-cell PCR and allele-specific PCR. In all ten individuals who carried a heteroplasmic D-loop mutation and a distantly located tRNA point mutation or a large deletion, we observed a mixture of four allelic combinations (tetraplasmy), a hallmark of recombination. Twelve of 14 individuals with closely located heteroplasmic D-loop mutation pairs contained a mixture of only three types of mitochondrial genomes (triplasmy), consistent with the absence of recombination between adjacent markers. These findings indicate that mtDNA recombination is common in human skeletal muscle.

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.

Repeats, longevity and the sources of mtDNA deletions: evidence from 'deletional spectra'.

Perfect direct repeats and, in particular, the prominent 13 bp repeat, are thought to cause mitochondrial DNA (mtDNA) deletions, which have been associated with the aging process. Accordingly, individuals lacking the 13 bp repeat are highly prevalent among centenarians and overall number of perfect repeats in mammalian mitochondrial genomes negatively correlates with species' longevity. However, detailed examination of the distribution of mtDNA deletions challenges a special role of the 13 bp repeat in generating mtDNA deletions. Instead, deletions appear to depend on long and stable, albeit imperfect, duplexes between distant mtDNA segments. Furthermore, significant dissimilarities in breakpoint distributions suggest that multiple mechanisms are involved in creating mtDNA deletions.

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.

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.

Mitochondrial DNA mutations and aging: devils in the details?

Although several lines of evidence support a role for accumulating somatic mitochondrial DNA (mtDNA) mutations in the etiology of aging, it remains unclear if they are a major cause of age-related deterioration and death. Mouse models that harbor elevated mtDNA mutation frequencies age prematurely; these findings were thought to provide conclusive evidence for a causal role of such mutations in aging. Yet, the presence of several conflicting reports has sparked controversy in the field and this is further aggravated by discrepancies in the estimates of mtDNA mutant fractions, which disagree by orders of magnitude. Here, we briefly review the evidence and some of the unresolved questions surrounding a causative role for accumulating mtDNA mutations in aging.

Clonally expanded mitochondrial DNA deletions within the choroid plexus in multiple sclerosis.

Mitochondrial DNA deletions (∆-mtDNA) have been implicated in the pathogenesis of Alzheimer's disease (AD), multiple sclerosis (MS) and Parkinson's disease (PD), as well as ageing. Clonal expansion of ∆-mtDNA is the process by which a mutant mtDNA molecule increases to high levels within a single cell containing both wild-type and mutant mtDNA. Unlike in AD and PD, the diffuse inflammatory process in MS involves the choroid plexus, and mitochondria are exposed to reactive oxygen and nitrogen species over a prolonged period. We determined the extent of respiratory enzyme deficiency and ∆-mtDNA at a single cell level within choroid plexus epithelial cells in MS as well as in AD, PD and controls. The respiratory enzyme-deficient (lacking complex IV and with intact complex II activity) cells were more prevalent within the choroid plexus in AD, MS and PD compared with controls. The main catalytic subunit of complex IV (subunit-I of cytochrome c oxidase) was lacking in significantly more respiratory enzyme-deficient cells in MS compared with AD, PD and controls. The single cell analysis showed a fourfold increase in the percentage of respiratory enzyme-deficient choroid plexus epithelial cells harbouring clonally expanded ∆-mtDNA in MS. Our findings establish clonal expansion of ∆-mtDNA as a feature relatively more prominent within the choroid plexus epithelium in MS than AD, PD or controls. We propose clonal expansion of ∆-mtDNA as a molecular link between inflammation and part of a delayed cellular energy failure in MS.

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).

Mitochondrial DNA deletions in mice in men: substantia nigra is much less affected in the mouse.

Mitochondrial DNA (mtDNA) deletions have been reported to accumulate to high levels in substantia nigra of older humans, and these mutations are suspected of causing age-related degeneration in this area. We have compared levels of mtDNA deletions in humans and mice and report here that levels of deletions in the mouse are very significantly lower than in humans. While human mtDNA from substantia nigra contained more than 5% of deleted molecules, mouse substantia nigra contained less than 0.5%. These results imply that mtDNA deletions are unlikely to play any significant role in of murine substantia nigra aging and further call for caution in using mouse models in studies of the role of mtDNA deletions in aging and neurodegeneration. On a more general note, these results support the view that critical targets of the various aging processes may differ significantly between distant species.

Distinct patterns of mitochondrial genome diversity in bonobos (Pan paniscus) and humans.

BACKGROUND: We have analyzed the complete mitochondrial genomes of 22 Pan paniscus (bonobo, pygmy chimpanzee) individuals to assess the detailed mitochondrial DNA (mtDNA) phylogeny of this close relative of Homo sapiens. RESULTS: We identified three major clades among bonobos that separated approximately 540,000 years ago, as suggested by Bayesian analysis. Incidentally, we discovered that the current reference sequence for bonobo likely is a hybrid of the mitochondrial genomes of two distant individuals. When comparing spectra of polymorphic mtDNA sites in bonobos and humans, we observed two major differences: (i) Of all 31 bonobo mtDNA homoplasies, i.e. nucleotide changes that occurred independently on separate branches of the phylogenetic tree, 13 were not homoplasic in humans. This indicates that at least a part of the unstable sites of the mitochondrial genome is species-specific and difficult to be explained on the basis of a mutational hotspot concept. (ii) A comparison of the ratios of non-synonymous to synonymous changes (dN/dS) among polymorphic positions in bonobos and in 4902 Homo sapiens mitochondrial genomes revealed a remarkable difference in the strength of purifying selection in the mitochondrial genes of the F0F1-ATPase complex. While in bonobos this complex showed a similar low value as complexes I and IV, human haplogroups displayed 2.2 to 7.6 times increased dN/dS ratios when compared to bonobos. CONCLUSIONS: Some variants of mitochondrially encoded subunits of the ATPase complex in humans very likely decrease the efficiency of energy conversion leading to production of extra heat. Thus, we hypothesize that the species-specific release of evolutionary constraints for the mitochondrial genes of the proton-translocating ATPase is a consequence of altered heat homeostasis in modern humans.