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.Translated from Genetika, Vol. 41, No. 1, 2005, pp. 93–99.Original Russian Text Copyright © 2005 by Malyarchuk.     

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.     

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.     

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     

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]     

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.     

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

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.     

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.     

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.     

Multiple Origins of a Mitochondrial Mutation Conferring Deafness

A point mutation (1555G) in the smaller ribosomal subunit of the mitochondrial DNA (mtDNA) has been associated with maternally inherited traits of hypersensitivity to streptomycin and sensorineural deafness in a number of families from China, Japan, Israel, and Africa. To determine whether this distribution was the result of a single or multiple mutational events, we carried out genetic distance analysis and phylogenetic analysis of 10 independent mtDNA D-loop sequences from Africa and Asia. The mtDNA sequence diversity was high (2.21%). Phylogenetic analysis assigned 1555G-bearing haplotypes at very divergent points in the human mtDNA evolutionary tree, and the 1555G mutations occur in many cases on race-specific mtDNA haplotypes, both facts are inconsistent with a recent introgression of the mutation into these races. The simplest interpretation of the available data is that there have been multiple origins of the 1555G mutation. The genetic distance among mtDNAs bearing the pathogenic 1555G mutation is much larger than among mtDNAs bearing either evolutionarily neutral or weakly deleterious nucleotide substitutions (such as the 4336G mutation). These results are consistent with the view that pathogenic mtDNA haplotypes such as 1555G arise on disparate mtDNA lineages which because of negative natural selection leave relatively few related descendants. The co-existence of the same mutation with deafness in individuals with very different nuclear and mitochondrial genetic backgrounds confirms the pathogenicity of the 1555G mutation.     

Mutation patterns of mitochondrial H- and L-strand DNA in closely related Cyprinid fishes

Mitochondrial genome replication is asymmetric. Replication starts from the origin of heavy (H)-strand replication, displacing the parental H-strand as it proceeds along the molecule. The H-strand remains single stranded until light (L)-strand replication is initiated from a second origin of replication. It has been suggested that single-stranded H-strand DNA is more sensitive to mutational damage, giving rise to substitutional rate differences between the two strands and among genes in mammalian mitochondrial DNA. In this study, we analyzed sequences of the cytochrome b, ND4, ND4L, and COI genes of cyprinid fishes to investigate rates and patterns of nucleotide substitution in the mitochondrial genome. To test for strand-asymmetric mutation pressure, a likelihood-ratio test was developed and applied to the cyprinid sequences. Patterns of substitution and levels of strand-asymmetric mutation pressure were largely consistent with a mutation gradient between the H- and L-strand origins of replication. Significant strand bias was observed among rates of transitional substitution. However, biological interpretation of the direction and strength of strand asymmetry for specific classes of substitutions is problematic. The problem occurs because the rate of any single class of substitution inferred from one strand is actually a sum of rates on two strands. The validity of the likelihood-ratio test is not affected by this problem.     

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.     

Rates of nuclear and cytoplasmic mitochondrial DNA sequence divergence in mammals

Differential rates of nucleotide substitution among different gene segments and between distinct evolutionary lineages is well documented among mitochondrial genes and is likely a consequence of locus-specific selective constraints that delimit mutational divergence over evolutionary time. We compared sequence variation of 18 homologous loci (15 coding genes and 3 parts of the control region) among 10 mammalian mitochondrial DNA genomes which allowed us to describe different mitochondrial evolutionary patterns and to produce an estimation of the relative order of gene divergence. The relative rates of divergence of mitochondrial DNA genes in the family Felidae were estimated by comparing their divergence from homologous counterpart genes included in nuclear mitochondrial DNA (Numt, pronounced "new might"), a genomic fossil that represents an ancient transfer of 7.9 kb of mitochondrial DNA to the nuclear genome of an ancestral species of the domestic cat (Felis catus). Phylogenetic analyses of mitochondrial (mtDNA) sequences with multiple outgroup species were conducted to date the ancestral node common to the Numt and the cytoplasmic (Cymt) mtDNA genes and to calibrate the rate of sequence divergence of mitochondrial genes relative to nuclear homologous counterparts. By setting the fastest substitution rate as strictly mutational, an empirical "selective retardation index" is computed to quantify the sum of all constraints, selective and otherwise, that limit sequence divergence of mitochondrial gene sequences over time.      

Low mutational burden of individual acquired mitochondrial DNA mutations in brain

Neurons may be particularly susceptible to oxidative damage, which has been proposed to induce somatic mutations, particularly in mitochondrial DNA (mtDNA). Therefore, acquired mtDNA mutations might preferentially accumulate in the brain and could play a role in aging and neurodegenerative disorders. Recently, a somatic T to G mtDNA mutation at noncoding nucleotide position 414 was reported in fibroblasts specifically from elderly subjects, with mutational burdens of up to 50%. We screened for this mutation in brain-derived mtDNA from 8 Alzheimer's disease patients, 27 Parkinson's disease patients, 4 multiple system atrophy patients, and 44 controls using up to three RFLP analyses. A total of 73 of these subjects were over the age of 65. The 414 mutation was absent in all cases. Next, individual mtDNA fragments from 6 elderly subjects were cloned, and a total of 70 clones were sequenced. The 414 mutation was absent in all clones, though occasional sequence variations were identified at other sites in single clones. The 414 mutation also was absent in blood (n = 6) and fibroblasts (n = 11) from elderly subjects. Our data suggest that it is rare for any one particular acquired mtDNA mutation to reach levels in the brain that are functionally significant. This does not exclude the possibility that the cumulative burden of multiple, individually rare, acquired mutations impairs mitochondrial function.     

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.