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.     

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.     

Variation of Human Mitochondrial DNA: Distribution of Hot Spots in Hypervariable Segment I of the Major Noncoding Region

Mitochondrial DNA (mtDNA) samples belonging to fifteen phylogenetically related mtDNA types specific to the populations of Europe (H, V, J, T, U, K, I, W, and X) and Northern Asia (A, C, D, G, Y, and Z) were typed for sequence variation in hypervariable segment I (HVSI). The approach used allowed to distinguish several hypervariable sites at nucleotide positions 16093, 16129, 16189, 16311, and 16362. Identical mutations at these sites were found in 10–11 out of 15 mtDNA groups examined. Positions 16126, 16172, 16192, 16256, 16261, 16291, 16293, and 16298 appeared to be less variable, since parallel mutations at these sites were found in 6–8 European and Asian mtDNA groups. The examples of the effects of mutations in hypervariable positions at the major noncoding mtDNA region on the frequency of reverse mutations in other mtDNA regions are presented. It was shown that such effects of nucleotide context on the mutation rate could be observed in phylogenetic mtDNA networks such as cyclic structures like rhombs and cubes. Analogous structures in the networks could be seen also in the case of the appearance of recombinant mtDNA types resulted from homologous recombination between mtDNA molecules in heteroplasmic mixture. The problem of the effect of polynucleotide context on the intensity of mtDNA mutagenesis is discussed. Recombination processes along with site-directed mutagenesis caused by action of genetic factors (of nuclear genome) and/or of the environment are considered as possible mechanisms of mitochondrial genome evolution.     

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.     

Mutant copies of mitochondrial DNA in tissues and plasma of mice subjected to X-ray irradiation

Impairments of mitochondrial genome are associated with a wide spectrum of degenerative diseases, development of tumors, aging, and cell death. We studied the content of mitochondrial DNA (mtDNA) with mutations and the total content of mutations in the brain and the spleen of mice subjected to X-ray irradiation at a dose of 1–5 Gy at 8–28 days after treatment. In these mice, we studied the number of mutant copies of extracellular mtDNA (ec-mtDNA) and its total content in blood plasma. We estimated mutations in control and irradiated mice using cleavage of heteroduplexes prepared by hybridization of PCR amplicons of mtDNA (D-loop region) mediated by CEL-I endonuclease, an enzyme that specifically cleaves unpaired bases. Changes in the total number of mtDNA copies relative to nuclear DNA were assessed by real time PCR using the ND-4 and GAPDH genes, respectively. We found that the number of mutant mtDNA copies was significantly increased in the brain and the spleen of irradiated mice and reached the maximum level at the eighth day after treatment; it then decreased by the 28th day after treatment. In nuclear genes similar to mutagenesis, mutagenesis of mtDNA in the brain and spleen tissues linearly depended on irradiation dose. In contrast to mutant nuclear genes, most mutant mtDNA copies were eliminated in the brain and spleen tissues, whereas the total content of mtDNA did not change within 28 days after irradiation. Our data show that, during this period, a high level of ec-mtDNA with mutations was observed in DNA circulating in blood plasma with the maximum level found at the 14th day. We suppose that mutant mtDNA copies are eliminated from cells of animals subjected to irradiation during the posttreatment period. Higher content of ec-mtDNA in blood plasma can be considered as a potential marker of radiation damage to the body.     

A method for mutagenesis of mouse mtDNA and a resource of mouse mtDNA mutations for modeling human pathological conditions

Mutations in human mitochondrial DNA (mtDNA) can cause mitochondrial disease and have been associated with neurodegenerative disorders, cancer, diabetes and aging. Yet our progress toward delineating the precise contributions of mtDNA mutations to these conditions is impeded by the limited availability of faithful transmitochondrial animal models. Here, we report a method for the isolation of mutations in mouse mtDNA and its implementation for the generation of a collection of over 150 cell lines suitable for the production of transmitochondrial mice. This method is based on the limited mutagenesis of mtDNA by proofreading-deficient DNA-polymerase γ followed by segregation of the resulting highly heteroplasmic mtDNA population by means of intracellular cloning. Among generated cell lines, we identify nine which carry mutations affecting the same amino acid or nucleotide positions as in human disease, including a mutation in the ND4 gene responsible for 70% of Leber Hereditary Optic Neuropathies (LHON). Similar to their human counterparts, cybrids carrying the homoplasmic mouse LHON mutation demonstrated reduced respiration, reduced ATP content and elevated production of mitochondrial reactive oxygen species (ROS). The generated resource of mouse mtDNA mutants will be useful both in modeling human mitochondrial disease and in understanding the mechanisms of ROS production mediated by mutations in mtDNA.     

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.     

Large scale mtDNA sequencing reveals sequence and functional conservation as major determinants of homoplasmic mtDNA variant distribution

The mitochondrial DNA (mtDNA) is highly variable, containing large numbers of pathogenic mutations and neutral polymorphisms. The spectrum of homoplasmic mtDNA variation was characterized in 730 subjects and compared with known pathogenic sites. The frequency and distribution of variants in protein coding genes were inversely correlated with conservation at the amino acid level. Analysis of tRNA secondary structures indicated a preference of variants for the loops and some acceptor stem positions. This comprehensive overview of mtDNA variants distinguishes between regions and positions which are likely not critical, mainly conserved regions with pathogenic mutations and essential regions containing no mutations at all.     

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.     

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.     

Quantification of random mutations in the mitochondrial genome

Mitochondrial DNA (mtDNA) mutations contribute to the pathology of a number of age-related disorders, including Parkinson disease [A. Bender et al., Nat. Genet. 38 (2006) 515,Y. Kraytsberg et al., Nat. Genet. 38 (2006) 518], muscle-wasting [J. Wanagat, Z. Cao, P. Pathare, J.M. Aiken, FASEB J. 15 (2001) 322], and the metastatic potential of cancers [K. Ishikawa et al., Science 320 (2008) 661]. The impact of mitochondrial DNA mutations on a wide variety of human diseases has made it increasingly important to understand the mechanisms that drive mitochondrial mutagenesis. In order to provide new insight into the etiology and natural history of mtDNA mutations, we have developed an assay that can detect mitochondrial mutations in a variety of tissues and experimental settings [M. Vermulst et al., Nat. Genet. 40 (2008) 4, M. Vermulst et al., Nat. Genet. 39 (2007) 540]. This methodology, termed the Random Mutation Capture assay, relies on single-molecule amplification to detect rare mutations among millions of wild-type bases [J.H. Bielas, L.A. Loeb, Nat. Methods 2 (2005) 285], and can be used to analyze mitochondrial mutagenesis to a single base pair level in mammals.     

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     

Pathogenic mitochondrial DNA mutations are common in the general population

Mitochondrial DNA (mtDNA) mutations are a major cause of genetic disease, but their prevalence in the general population is not known. We determined the frequency of ten mitochondrial point mutations in 3168 neonatal-cord-blood samples from sequential live births, analyzing matched maternal-blood samples to estimate the de novo mutation rate. mtDNA mutations were detected in 15 offspring (0.54%, 95% CI = 0.30–0.89%). Of these live births, 0.00107% (95% CI = 0.00087–0.0127) harbored a mutation not detected in the mother's blood, providing an estimate of the de novo mutation rate. The most common mutation was m.3243A→G. m.14484T→C was only found on sub-branches of mtDNA haplogroup J. In conclusion, at least one in 200 healthy humans harbors a pathogenic mtDNA mutation that potentially causes disease in the offspring of female carriers. The exclusive detection of m.14484T→C on haplogroup J implicates the background mtDNA haplotype in mutagenesis. These findings emphasize the importance of developing new approaches to prevent transmission.     

Mutation analysis of the entire mitochondrial genome using denaturing high performance liquid chromatography

In patients with mitochondrial disease a continuously increasing number of mitochondrial DNA (mtDNA) mutations and polymorphisms have been identified. Most pathogenic mtDNA mutations are heteroplasmic, resulting in heteroduplexes after PCR amplification of mtDNA. To detect these heteroduplexes, we used the technique of denaturing high performance liquid chromatography (DHPLC). The complete mitochondrial genome was amplified in 13 fragments of 1–2 kb, digested in fragments of 90–600 bp and resolved at their optimal melting temperature. The sensitivity of the DHPLC system was high with a lowest detection of 0.5% for the A8344G mutation. The muscle mtDNA from six patients with mitochondrial disease was screened and three mutations were identified. The first patient with a limb-girdle-type myopathy carried an A3302G substitution in the tRNA^Leu(UUR) gene (70% heteroplasmy), the second patient with mitochondrial myopathy and cardiomyopathy carried a T3271C mutation in the tRNA^Leu(UUR) gene (80% heteroplasmy) and the third patient with Leigh syndrome carried a T9176C mutation in the ATPase6 gene (93% heteroplasmy). We conclude that DHPLC analysis is a sensitive and specific method to detect heteroplasmic mtDNA mutations. The entire automatic procedure can be completed within 2 days and can also be applied to exclude mtDNA involvement, providing a basis for subsequent investigation of nuclear genes.     

Variation of human mitochondrial DNA: distribution of hot spots in hypervariable segment I of the major noncoding region

Mitochondrial DNA (mtDNA) samples belonging to fifteen phylogenetically related mtDNA types specific to the populations of Europe (H, V, J, T, U, K, I, W, and X) and Northern Asia (A, C, D, G, Y, and Z) were typed for sequence variation in hypervariable segment I (HVSI). The approach used allowed to distinguish several hypervariable sites at nucleotide positions 16093, 16129, 16189, 16311, and 16362. Identical mutations at these sites were found in 10-11 out of 15 mtDNA groups examined. Positions 16126, 16172, 16192, 16256, 16261, 16291, 16293, and 16298 appeared to be less variable, since parallel mutations at these sites were found in 6-8 European and Asian mtDNA groups. The examples of the effects of mutations in hypervariable positions at the major noncoding mtDNA region on the frequency of reverse mutations in other mtDNA regions are presented. It was shown that such effects of nucleotide context on the mutation rate could be observed in phylogenetic mtDNA networks such as cyclic structures like rhombs and cubes. Analogous structures in the networks could be seen also in the case of the appearance of recombinant mtDNA types resulted from homologous recombination between mtDNA molecules in heteroplasmic mixture. The problem of the effect of polynucleotide context on the intensity of mtDNA mutagenesis is discussed. Recombination processes along with site-directed mutagenesis caused by action of genetic factors (of nuclear genome) and/or of the environment are considered as possible mechanisms of mitochondrial genome evolution.     

The effect of the Primary Structure of DNA on Induction of Mutations by Alkylating Agents

Analysis and comparison of mutation spectra is one of the major tasks of molecular biology, since mutation spectra often reveal important properties of various mutagens and proteins involved in the repair/replication systems. Mutability is known to vary significantly along the nucleotide sequence. Mutations are abundant at certain positions (mutation hotspots). In this work, we applied regression analysis based on the basic logic patterns to understand the role of the nucleotide sequence context in mutation induction. The spectra of mutations induced by various alkylating agents were studied. The nucleotide bases at positions –2, –1, +1 and +2 were shown to have the most significant effect in G : C A : T replacements.     

Analysis of the mitochondrial DNA from patients with Wolfram (DIDMOAD) syndrome

Wolfram or DIDMOAD (Diabetes Insipidus, Diabetes Mellitus, Optic Atrophy and Deafness) syndrome, which has long been known as an autosomal-recessive disorder, has recently been proposed to be a mitochondrial-mediated disease with either a nuclear or a mitochondrial genetic background. The phenotypic characteristics of the syndrome resemble those found in other mitochondrial (mt)DNA mediated disorders such as Leber's hereditary optic neuropathy (LHON) or maternally inherited diabetes and deafness (MIDD). Therefore, we looked for respective mtDNA alterations in blood samples from 7 patients with DIDMOAD syndrome using SSCP-analysis of PCR-amplified fragments, encompassing all mitochondrial ND and tRNA genes, followed by direct sequencing. Subsequently, we compared mtDNA variants identified in this disease group with those detected in a group of LHON patients (n = 17) and in a group of 69 healthy controls. We found that 4/7 (57%) DIDMOAD patients harbored a specific set of point mutations in tRNA and ND genes including the so-called class II or secondary LHON mutations at nucleotide positions (nps) 4216 and 4917 (haplogroup B). In contrast, LHON-patients were frequently (10/17, 59%) found in association with another cluster of mtDNA variants including the secondary LHON mutations at nps 4216 and 13708 and further mtDNA polymorphisms in ND genes (haplogroup A), overlapping with haplogroup B only by variants at nps 4216 and 11251. The frequencies of both haplogroups were significantly lower in the control group versus the respective disease groups. We propose that haplogroup B represents a susceptibility factor for DIDMOAD which, by interaction with further exogeneous or genetic factors, might increase the risk for disease. (Mol Cell Biochem 174: 209–213, 1997)     

Evolution of Ig DNA sequence to target specific base positions within codons for somatic hypermutation

Ig variable (V) region genes are subjected to a somatic hypermutation process as B lymphocytes participate in immune reactions to protein Ags. Although little is known regarding the mechanism of mutagenesis, a consistent hierarchy of trinucleotide target preferences is evident. Analysis of trinucleotide regional distributions predicted and we now empirically confirm the surprising finding that the framework 2 region of kappa V region genes is highly mutable despite its importance to the structural integrity and function of the Ab molecule. Interestingly, much of this mutability appears to be focused on the third codon position where synonymous substitutions are most likely to occur. We also observed a trend for high predicted mutability for codon positions 1 and 2 in complementarity-determining regions. Consequently, amino acid replacements should occur at a higher rate in complementarity-determining regions than in framework regions due to the distribution and subsequent targeting of microsequences by the mutation mechanism. Our results reveal a subtle tier of V region gene evolution in which DNA sequence has been molded to direct mutations to specific base positions within codons in a manner that minimizes damage and maximizes the benefits of the somatic hypermutation process.