Distribution of nucleotide substitutions in human mitochondrial DNA genes

To analyze the distribution pattern of nucleotide substitutions in human mitochondrial DNA (mtDNA), mutational spectra of the mitochondrial genes were reconstructed. The reconstruction procedure is based on the mutation distribution data for 47 monophyletic mtDNA clusters, to which 794 examined mtDNA sequences encoding for tRNAs, rRNAs, and mitochondrial proteins are attributed. One of specific features of mitochondrial mutational spectra revealed was homoplasy of the mutations (the mean mutation number per variable nucleotide site in the coding region varied from 1.09 to 1.43). It was established that in the mtDNA genes maximum mutational constraint fell onto the guanine bases, albeit the content of these bases in the mtDNA L-chains was minimal. Maximal bias towards parallel G to A transitions was observed for rRNA genes, with the protein- and tRNA-encoding genes ranking next. Despite the fact that the differences in the average G-nucleotides content and variability between the genes of two mtDNA segments located between the OriH and OriL were statistically significant, the results did not provide the conclusion that the G-nucleotide instability observed in the mtDNA L-spectra was determined by the mechanism of asynchronous mtDNA replication, along with the deamination of cytosines in the H-chain regions, which remained single-stranded during replication.     

Analysis of mutation mechanisms in human mitochondrial DNA

The reasons of high level of human mitochondrial DNA (mtDNA) variability remain to be largely unclear. We analyze here three probable mechanisms of mutagenesis leading to generation of mtDNA nucleotide substitutions: (1) deamination of DNA bases; (2) tautomeric migrations of protons in nitrous bases; and (3) hydrolysis of glycoside link between DNA bases and carbohydrate residue on the background of free radical damage of the mitochondrial DNA polymerase gamma. By means of quanto-chemical calculations, it was shown that the most substantiated mechanism of mutation generation is hydrolysis of N-glycoside link. This mechanism is suggestive to be more prominent on the H-strand, which remains to be single-stranded for a long time during the mtDNA replication. It was revealed also that hydrolytic deamination of adenines on the single-stranded H-strand is among of the most probable mechanisms leading to high frequency of T --> C transitions seen in the L-strand mutational spectra of the mtDNA major non-coding region.     

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.     

The role of nucleotide context in the induction of mutations in human mitochondrial DNA genes

Based on the mutations distribution patterns in the mitochondrial DNA (mtDNA) genes, context analysis of the regions, including mutable positions characterized by the appearance of more than two parallel mutations, was performed. It was demonstrated that the mechanism of dislocation mutagenesis, leading to the appearance of mismatches within the frameshift regions of either primer or template mtDNA chains during replication, accounts for the induction of 21% of unstable positions in the mtDNA genes. Context analysis showed that pyrimidine bases in the positions +1 and +2 (gYRNS, gYY, and gR consensuses, where g is mutable position) had the highest influence on the induction of mutations in G positions of the mtDNA genes. The highest effect on the mutagenesis in T positions was excreted by the bases in the positions –1 and +1 (RtY and tA consensuses, where t is mutable position). In general, these data point to the prevalence of the context-dependant mechanisms of the mutations induction in human mitochondrial genome.Translated from Genetika, Vol. 41, No. 3, 2005, pp. 385–390.Original Russian Text Copyright © 2005 by Malyarchuk.     

High mutation rates in the mitochondrial genomes of Daphnia pulex

Despite the great utility of mitochondrial DNA (mtDNA) sequence data in population genetics and phylogenetics, key parameters describing the process of mitochondrial mutation (e.g., the rate and spectrum of mutational change) are based on few direct estimates. Furthermore, the variation in the mtDNA mutation process within species or between lineages with contrasting reproductive strategies remains poorly understood. In this study, we directly estimate the mtDNA mutation rate and spectrum using Daphnia pulex mutation-accumulation (MA) lines derived from sexual (cyclically parthenogenetic) and asexual (obligately parthenogenetic) lineages. The nearly complete mitochondrial genome sequences of 82 sexual and 47 asexual MA lines reveal high mtDNA mutation rate of 1.37 × 10(-7) and 1.73 × 10(-7) per nucleotide per generation, respectively. The Daphnia mtDNA mutation rate is among the highest in eukaryotes, and its spectrum is dominated by insertions and deletions (70%), largely due to the presence of mutational hotspots at homopolymeric nucleotide stretches. Maximum likelihood estimates of the Daphnia mitochondrial effective population size reveal that between five and ten copies of mitochondrial genomes are transmitted per female per generation. Comparison between sexual and asexual lineages reveals no statistically different mutation rates and highly similar mutation spectra.     

Analysis of Mutation Mechanisms in Human Mitochondrial DNA

The cause of the high variability of human mitochondrial DNA (mtDNA) remains largely unknown. Three mechanisms of mutagenesis that might account for the generation of nucleotide substitutions in mtDNA have been analyzed: deamination of DNA nitrous bases caused by deamination agents, tautomeric proton migration in nitrous bases, and the hydrolysis of the glycoside bond between the nitrous base and carbohydrate residue in nucleotides against the background of the free-radical damage of DNA polymerase γ. Quantum chemical calculations demonstrated that the hydrolysis of the N-glycoside bond is the most probable mechanism; it is especially prominent in the H strand, which remains free during mtDNA replication for a relatively long time. It has also been found that hydrolytic deamination of adenine in single-stranded regions of the H strand is a possible cause of the high frequency of T → C transitions in the mutation spectra of the L-chain of the major mtDNA noncoding region.     

Direct estimation of the mitochondrial DNA mutation rate in Drosophila melanogaster

Mitochondrial DNA (mtDNA) variants are widely used in evolutionary genetics as markers for population history and to estimate divergence times among taxa. Inferences of species history are generally based on phylogenetic comparisons, which assume that molecular evolution is clock-like. Between-species comparisons have also been used to estimate the mutation rate, using sites that are thought to evolve neutrally. We directly estimated the mtDNA mutation rate by scanning the mitochondrial genome of Drosophila melanogaster lines that had undergone approximately 200 generations of spontaneous mutation accumulation (MA). We detected a total of 28 point mutations and eight insertion-deletion (indel) mutations, yielding an estimate for the single-nucleotide mutation rate of 6.2 × 10⁻⁸ per site per fly generation. Most mutations were heteroplasmic within a line, and their frequency distribution suggests that the effective number of mitochondrial genomes transmitted per female per generation is about 30. We observed repeated occurrences of some indel mutations, suggesting that indel mutational hotspots are common. Among the point mutations, there is a large excess of G→A mutations on the major strand (the sense strand for the majority of mitochondrial genes). These mutations tend to occur at nonsynonymous sites of protein-coding genes, and they are expected to be deleterious, so do not become fixed between species. The overall mtDNA mutation rate per base pair per fly generation in Drosophila is estimated to be about 10× higher than the nuclear mutation rate, but the mitochondrial major strand G→A mutation rate is about 70× higher than the nuclear rate. Silent sites are substantially more strongly biased towards A and T than nonsynonymous sites, consistent with the extreme mutation bias towards A+T. Strand-asymmetric mutation bias, coupled with selection to maintain specific nonsynonymous bases, therefore provides an explanation for the extreme base composition of the mitochondrial genome of Drosophila.     

Mitochondrial Gene Rearrangements and Partial Genome Duplications Detected by Multigene Asymmetric Compositional Bias Analysis

Asymmetric compositional and mutation bias between the two strands occurs in mitochondrial genomes, and an asymmetric mechanism of mtDNA replication is a potential source of this bias. Some evidence indicates that during replication the heavy strand is subject to a gradient of time spent in a single-stranded state (D _ssH) and a gradient of mutational damage. The nucleotide composition bias among genes varies with D _ssH. Consequently, partial genome duplications (PGD) will alter the skew for genes located downstream of the duplication, relatively to nascent light strand synthesis, and in the same way, gene rearrangements (GRr) will affect genes by changing their skews. We examined cases where there had been PGD or GRr and determined whether this left a trace in the form of unusual patterns of base composition. We compared the skew of genes differently located on the mtDNA genome of previously published whole mtDNA genomes from amphibians, a group that shows considerable levels of both GRr and PGD. After observing a significant correlation between AT and GC skew with D _ssH at fourfold redundant sites, we ran our analysis and detected 31.3% of the species with GRr and/or PGD. By comparing the nucleotide composition at fourfold redundant sites in normal and “abnormal” species, we found that A/C variation occurs and is associated with GRr/PGD. These results show that by analyzing the nucleotide skews of only three genes, it may be possible to predict some mitochondrial GRr and/or PGD without knowing the complete mtDNA genome sequence. [Reviewing Editor: Dr. David Pollock]     

A Genetic Bottleneck of Mitochondrial DNA During Human Lymphocyte Development

Mitochondria are essential organelles in eukaryotic cells that provide critical support for energetic and metabolic homeostasis. Although the elimination of pathogenic mitochondrial DNA (mtDNA) mutations in somatic cells has been observed, the mechanisms to maintain proper functions despite their mtDNA mutation load are poorly understood. In this study, we analyzed somatic mtDNA mutations in more than 30,000 single human peripheral and bone marrow mononuclear cells. We observed a significant overrepresentation of homoplasmic mtDNA mutations in B, T, and natural killer (NK) lymphocytes. Intriguingly, their overall mutational burden was lower than that in hematopoietic progenitors and myeloid cells. This characteristic mtDNA mutational landscape indicates a genetic bottleneck during lymphoid development, as confirmed with single-cell datasets from multiple platforms and individuals. We further demonstrated that mtDNA replication lags behind cell proliferation in both pro-B and pre-B progenitor cells, thus likely causing the genetic bottleneck by diluting mtDNA copies per cell. Through computational simulations and approximate Bayesian computation (ABC), we recapitulated this lymphocyte-specific mutational landscape and estimated the minimal mtDNA copies as <30 in T, B, and NK lineages. Our integrative analysis revealed a novel process of a lymphoid-specific mtDNA genetic bottleneck, thus illuminating a potential mechanism used by highly metabolically active immune cells to limit their mtDNA mutation load.     

Analysis of phylogenetically reconstructed mutational spectra in human mitochondrial DNA control region

Analysis of mutations in mitochondrial DNA is an important issue in population and evolutionary genetics. To study spontaneous base substitutions in human mitochondrial DNA we reconstructed the mutational spectra of the hypervariable segments I and II (HVS I and II) using published data on polymorphisms from various human populations. An excess of pyrimidine transitions was found both in HVS I and II regions. By means of classification analysis numerous mutational hotspots were revealed in these spectra. Context analysis of hotspots revealed a complex influence of neighboring bases on mutagenesis in the HVS I region. Further statistical analysis suggested that a transient misalignment dislocation mutagenesis operating in monotonous runs of nucleotides play an important role for generating base substitutions in mitochondrial DNA and define context properties of mtDNA. Our results suggest that dislocation mutagenesis in HVS I and II is a fingerprint of errors produced by DNA polymerase gamma in the course of human mitochondrial DNA replication     

Twinkle and POLG defects enhance age-dependent accumulation of mutations in the control region of mtDNA

Autosomal dominant and/or recessive progressive external ophthalmoplegia (ad/arPEO) is associated with mtDNA mutagenesis. It can be caused by mutations in three nuclear genes, encoding the adenine nucleotide translocator 1, the mitochondrial helicase Twinkle or DNA polymerase γ (POLG). How mutations in these genes result in progressive accumulation of multiple mtDNA deletions in post- mitotic tissues is still unclear. A recent hypothesis suggested that mtDNA replication infidelity could promote slipped mispairing, thereby stimulating deletion formation. This hypothesis predicts that mtDNA of ad/arPEO patients will contain frequent mutations throughout; in fact, our analysis of muscle from ad/arPEO patients revealed an age-dependent, enhanced accumulation of point mutations in addition to deletions, but specifically in the mtDNA control region. Both deleted and non-deleted mtDNA molecules showed increased point mutation levels, as did mtDNAs of patients with a single mtDNA deletion, suggesting that point mutations do not cause multiple deletions. Deletion breakpoint analysis showed frequent breakpoints around homopolymeric runs, which could be a signature of replication stalling. Therefore, we propose replication stalling as the principal cause of deletion formation.     

Mutational process in protein-coding genes of human mitochondrial genome in context of evolution of Homo genus

The human mitochondrial genome, although small in size, shows a high level of variation that differs across nucleotide groups. In this work, mutation rates in mtDNA were compared in species of the Homo genus, including humans, Neanderthals, Denisova hominins, and other primate species. It was found that more than half (56.5%) of the polymorphisms in protein-coding genes of human mtDNA are actually reverse mutations to the pre-H. sapiens state of the mitochondrial genome. Among hypervariable nucleotide positions, only a small portion of mutations are specific to H. sapiens, while the majority of mutations (both nucleotide and amino acid substitutions) result in a loss of Homo-specific variants of polymorphisms. Most commonly, polymorphism variants specific to H. sapiens arise as a result of unique forward mutations and disappear mainly due to multiple reverse mutations, including those in mutational hot spots.     

The MYC dualism in growth and death

Investigations of intraindividual sequence diversity in mtDNA are a key step in exploring the linkage between somatic mutations in mtDNA and mitochondrial genome evolution. This paper reports a directional cloning procedure enabling the isolation of multiple copies of the D-loop region of the mitochondrial genome from the fish Ameiurus nebulosus. Sequence analysis of 708 D-loop molecules revealed eight mutants, an average intraindividual mutation frequency of 1.12%. Three different types of mutations were detected but each derived from a single mutational event. By contrasting the spectrum of nucleotide variation at multiple biological levels, one can investigate the effects of spontaneous mutations on genome evolution. Such hierarchical analysis suggested shifts in the type and distribution of mtDNA (mitochondrial DNA) mutations at different biological levels, indicating the need to recognize three different rates of mtDNA sequence change from the cellular to population level.     

Deficiency of the ADP-forming succinyl-CoA synthase activity is associated with encephalomyopathy and mitochondrial DNA depletion

The mitochondrial DNA (mtDNA) depletion syndrome is a quantitative defect of mtDNA resulting from dysfunction of one of several nuclear-encoded factors responsible for maintenance of mitochondrial deoxyribonucleoside triphosphate (dNTP) pools or replication of mtDNA. Markedly decreased succinyl-CoA synthetase activity due to a deleterious mutation in SUCLA2, the gene encoding the beta subunit of the ADP-forming succinyl-CoA synthetase ligase, was found in muscle mitochondria of patients with encephalomyopathy and mtDNA depletion. Succinyl-CoA synthetase is invariably in a complex with mitochondrial nucleotide diphosphate kinase; hence, we propose that a defect in the last step of mitochondrial dNTP salvage is a novel cause of the mtDNA depletion syndrome.     

The role of nucleotide context in the induction of mutations in human mitochondrial DNA genes

Based on the mutations distribution patterns in the mitochondrial DNA (mtDNA) genes, context analysis of the regions, including mutable positions characterized by the appearance of more than two parallel mutations, was performed. It was demonstrated that the mechanism of dislocation mutagenesis, leading to the appearance of mismatches within the frameshift regions of either primer or template mtDNA chains during replication, accounts for the induction of 21% of unstable positions in the mtDNA genes. Context analysis showed that pyrimidine bases in the positions +1 and +2 (gYRNS, gYY, and gR consensuses, where g is mutable position) had the highest influence on the induction of mutations in G positions of the mtDNA genes. The highest effect on the mutagenesis in T positions was excreted by the bases in the positions -1 and +1 (RyT and tA consensuses, where t is mutable position). In general, these data point to the prevalence of the context-dependant mechanisms of the mutations induction in human mitochondrial genome.     

Comparative analysis of the hypervariable segment 1 mutational spectra in human mitochondrial DNA phylogeographic groups

Variability of the mtDNA hypervariable segment 1 (HVS 1) nucleotide sequences belonging to 88 phylogeographic clusters characteristic for human populations of Africa, West and East Eurasia was analyzed. Statistically significant differences between distribution of mutations in mitochondrial gene pools of the human continental groups were revealed. The list of the HVS 1 nucleotide positions characterizing by instability explained by the model of mtDNA strands dislocation during the replication process is suggested. It was shown that DNA strands dislocation during mtDNA replication is one of the key mechanisms of the context-dependent mtDNA mutagenesis during the regional differentiation of human populations.     

The complete mitochondrial DNA sequence of the horseshoe crab Limulus polyphemus

We determined the complete 14,985-nt sequence of the mitochondrial DNA of the horseshoe crab Limulus polyphemus (Arthropoda: Xiphosura). This mtDNA encodes the 13 protein, 2 rRNA, and 22 tRNA genes typical for metazoans. The arrangement of these genes and about half of the sequence was reported previously; however, the sequence contained a large number of errors, which are corrected here. The two strands of Limulus mtDNA have significantly different nucleotide compositions. The strand encoding most mitochondrial proteins has 1. 25 times as many A's as T's and 2.33 times as many C's as G's. This nucleotide bias correlates with the biases in amino acid content and synonymous codon usage in proteins encoded by different strands and with the number of non-Watson-Crick base pairs in the stem regions of encoded tRNAs. The sizes of most mitochondrial protein genes in Limulus are either identical to or slightly smaller than those of their Drosophila counterparts. The usage of the initiation and termination codons in these genes seems to follow patterns that are conserved among most arthropod and some other metazoan mitochondrial genomes. The noncoding region of Limulus mtDNA contains a potential stem-loop structure, and we found a similar structure in the noncoding region of the published mtDNA of the prostriate tick Ixodes hexagonus. A simulation study was designed to evaluate the significance of these secondary structures; it revealed that they are statistically significant. No significant, comparable structure can be identified for the metastriate ticks Rhipicephalus sanguineus and Boophilus microplus. The latter two animals also share a mitochondrial gene rearrangement and an unusual structure of mt-tRNA(C) that is exactly the same association of changes as previously reported for a group of lizards. This suggests that the changes observed are not independent and that the stem-loop structure found in the noncoding regions of Limulus and Ixodes mtDNA may play the same role as that between trnN and trnC in vertebrates, i.e., the role of lagging strand origin of replication.     

Somatic mtDNA Mutation Spectra in the Aging Human Putamen

The accumulation of heteroplasmic mitochondrial DNA (mtDNA) deletions and single nucleotide variants (SNVs) is a well-accepted facet of the biology of aging, yet comprehensive mutation spectra have not been described. To address this, we have used next generation sequencing of mtDNA-enriched libraries (Mito-Seq) to investigate mtDNA mutation spectra of putamen from young and aged donors. Frequencies of the “common” deletion and other “major arc” deletions were significantly increased in the aged cohort with the fold increase in the frequency of the common deletion exceeding that of major arc deletions. SNVs also increased with age with the highest rate of accumulation in the non-coding control region which contains elements necessary for translation and replication. Examination of predicted amino acid changes revealed a skew towards pathogenic SNVs in the coding region driven by mutation bias. Levels of the pathogenic m.3243A>G tRNA mutation were also found to increase with age. Novel multimeric tandem duplications that resemble murine control region multimers and yeast ρ⁻ mtDNAs, were identified in both young and aged specimens. Clonal ∼50 bp deletions in the control region were found at high frequencies in aged specimens. Our results reveal the complex manner in which the mitochondrial genome alters with age and provides a foundation for studies of other tissues and disease states.