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.     

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.     

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.     

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.     

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.     

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.     

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.     

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     

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.     

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.     

Mitochondrial dysfunction due to oxidative mitochondrial DNA damage is reduced through cooperative actions of diverse proteins

The mitochondrial genome is a significant target of exogenous and endogenous genotoxic agents; however, the determinants that govern this susceptibility and the pathways available to resist mitochondrial DNA (mtDNA) damage are not well characterized. Here we report that oxidative mtDNA damage is elevated in strains lacking Ntg1p, providing the first direct functional evidence that this mitochondrion-localized, base excision repair enzyme functions to protect mtDNA. However, ntg1 null strains did not exhibit a mitochondrial respiration-deficient (petite) phenotype, suggesting that mtDNA damage is negotiated by the cooperative actions of multiple damage resistance pathways. Null mutations in ABF2 or PIF1, two genes implicated in mtDNA maintenance and recombination, exhibit a synthetic-petite phenotype in combination with ntg1 null mutations that is accompanied by enhanced mtDNA point mutagenesis in the corresponding double-mutant strains. This phenotype was partially rescued by malonic acid, indicating that reactive oxygen species generated by the electron transport chain contribute to mitochondrial dysfunction in abf2 Delta strains. In contrast, when two other genes involved in mtDNA recombination, CCE1 and NUC1, were inactivated a strong synthetic-petite phenotype was not observed, suggesting that the effects mediated by Abf2p and Pif1p are due to novel activities of these proteins other than recombination. These results document the existence of recombination-independent mechanisms in addition to base excision repair to cope with oxidative mtDNA damage in Saccharomyces cerevisiae. Such systems are likely relevant to those operating in human cells where mtDNA recombination is less prevalent, validating yeast as a model system in which to study these important issues.     

Lineage-specific selection in human mtDNA: lack of polymorphisms in a segment of MTND5 gene in haplogroup J

Human mitochondrial DNA (mtDNA) is a nonrecombining genome that codes for 13 subunits of the mitochondrial oxidative phosphorylation system, 2 rRNAs, and 22 tRNAs. Mutations have accumulated sequentially in mtDNA lineages that diverged tens of thousands of years ago. The genes in mtDNA are subject to different functional constraints and are therefore expected to evolve at different rates, but the rank order of these rates should be the same in all lineages of a phylogeny. Previous studies have indicated, however, that specific regions of mtDNA may have experienced different histories of selection in different lineages, possibly because of lineage-specific interactions or environmental factors such as climate. We report here on a survey for lineage-specific patterns of nucleotide polymorphism in human mtDNA. We calculated molecular polymorphism indices and neutrality tests for classes of functional sites and genes in 837 human mtDNA sequences, compared the results between continent-specific mtDNA lineages, and used two sliding window methods to identify differences in the patterns of polymorphism between haplogroups. A general correlation between nucleotide position and the level of nucleotide polymorphism was identified in the coding region of the mitochondrial genome. Nucleotide diversity in the protein-coding sequence of mtDNA was generally not much higher than that found for many genes in nuclear DNA. A comparison of nonsynonymous/synonymous rate ratios in the 13 protein-coding genes suggested differences in the relative levels of selection between haplogroups, including the European haplogroup clusters. Interestingly, a segment of the MTND5 gene was found to be almost void of segregating sites and nonsynonymous mutations in haplogroup J, which has been associated with susceptibility to certain complex diseases. Our results suggest that there are haplogroup-specific differences in the intensity of selection against particular regions of the mitochondrial genome, indicating that some mutations may be non-neutral within specific phylogenetic lineages but neutral within others.     

MMS Exposure Promotes Increased MtDNA Mutagenesis in the Presence of Replication-Defective Disease-Associated DNA Polymerase γ Variants

Mitochondrial DNA (mtDNA) encodes proteins essential for ATP production. Mutant variants of the mtDNA polymerase cause mutagenesis that contributes to aging, genetic diseases, and sensitivity to environmental agents. We interrogated mtDNA replication in Saccharomyces cerevisiae strains with disease-associated mutations affecting conserved regions of the mtDNA polymerase, Mip1, in the presence of the wild type Mip1. Mutant frequency arising from mtDNA base substitutions that confer erythromycin resistance and deletions between 21-nucleotide direct repeats was determined. Previously, increased mutagenesis was observed in strains encoding mutant variants that were insufficient to maintain mtDNA and that were not expected to reduce polymerase fidelity or exonuclease proofreading. Increased mutagenesis could be explained by mutant variants stalling the replication fork, thereby predisposing the template DNA to irreparable damage that is bypassed with poor fidelity. This hypothesis suggests that the exogenous base-alkylating agent, methyl methanesulfonate (MMS), would further increase mtDNA mutagenesis. Mitochondrial mutagenesis associated with MMS exposure was increased up to 30-fold in mip1 mutants containing disease-associated alterations that affect polymerase activity. Disrupting exonuclease activity of mutant variants was not associated with increased spontaneous mutagenesis compared with exonuclease-proficient alleles, suggesting that most or all of the mtDNA was replicated by wild type Mip1. A novel subset of C to G transversions was responsible for about half of the mutants arising after MMS exposure implicating error-prone bypass of methylated cytosines as the predominant mutational mechanism. Exposure to MMS does not disrupt exonuclease activity that suppresses deletions between 21-nucleotide direct repeats, suggesting the MMS-induce mutagenesis is not explained by inactivated exonuclease activity. Further, trace amounts of CdCl₂ inhibit mtDNA replication but suppresses MMS-induced mutagenesis. These results suggest a novel mechanism wherein mutations that lead to hypermutation by DNA base-damaging agents and associate with mitochondrial disease may contribute to previously unexplained phenomena, such as the wide variation of age of disease onset and acquired mitochondrial toxicities.     

Variation of the mitochondrial genome in the evolution of Drosophila

The evidence on mitochondrial genome variation and its role in evolution of the genus Drosophila are reviewed. The mitochondrial genome is represented by a circular double-stranded DNA molecule 16 to 19 kb in length. The genome contains no introns involved in recombination. The entire mitochondrial genome can be arbitrarily divided into three parts: (1) protein-coding genes; (2) genes encoding rRNA and tRNA; and (3) the noncoding regulatory region (A + T region). The selective importance of mutations within different mtDNA regions is therefore unequal. In Drosophila, the content of the A + T pairs in mtDNA is extremely low and a pattern of nucleotide substitution is characterized by a low transition/transversion ratio (and a low threshold of mutation saturation). The deletions and duplications are of common occurrence in the mitochondrial genome. However, this genome lacks such characteristic for the nuclear genome aberrations as the inversions and transpositions. The phenomena of introgression and heteroplasmy provide an opportunity to study the adaptive role of the mitochondrial genome and its role in speciation. Analysis of evidence concerning mtDNA variation in different species of the genus Drosophila made it possible to ascertain data on phylogenetic relationships among species obtained by studying nuclear genome variation. In some species, mtDNA variation may serve as a reliable marker for population differentiation within a species, although evidence on the population dynamics of the mtDNA variation is very scarce.     

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.     

Information content of amino acid residues in putative helix VIII of the lac permease from Escherichia coli

Mutants in putative helix VIII of lactose permease that retain the ability to accumulate lactose were created by cassette mutagenesis. A mutagenic insert encoding amino acid residues 259-278 was synthesized chemically by using reagents contaminated with 1% each of the other three bases and ligated into a KpnI/BclI site in the lacY gene in plasmid pGEM-4. Mutants that retain transport activity were selected by transforming a strain of Escherichia coli containing a wild-type lacZ gene, but deleted in lacY, with the mutant library and identifying colonies that transport lactose on indicator plates. Sequencing of the mutated region in lacY in 129 positive colonies reveals 43 single amino acid mutations at 26 sites and 26 multiple mutations. The variable amino acid positions are largely on one side of the putative alpha-helix, a stripe opposite Glu269. This mutable stripe of low information content is probably in contact with the membrane phospholipids.     

Variation of the Mitochondral Genome in the Evolution of Drosophila

The evidence on mitochondrial genome variation and its role in evolution of the genus Drosophila are reviewed. The mitochondrial genome is represented by a circular double-stranded DNA molecule 16 to 19 kb in length. Mitochondrial genes lack introns and recombination. The entire mitochondrial genome can be arbitrarily divided into three parts: (1) protein-coding genes; (2) genes encoding rRNA and tRNA; and (3) the noncoding regulatory region (A + T region). The selective importance of mutations within different mtDNA regions is therefore unequal. In Drosophila, the content of the A + T pairs in mtDNA is extremely high and a pattern of nucleotide substitution is characterized by a low transition/transversion ratio (and a low threshold of mutation saturation). The deletions and duplications are of common occurrence in the mitochondrial genome. However, this genome lacks such characteristic for the nuclear genome aberrations as inversions and transpositions. The phenomena of introgression and heteroplasmy provide an opportunity to study the adaptive role of the mitochondrial genome and its role in speciation. Analysis of evidence concerning mtDNA variation in different species of the genus Drosophilamade it possible to ascertain data on phylogenetic relationships among species obtained by studying nuclear genome variation. In some species, mtDNA variation may serve as a reliable marker for population differentiation within a species, although evidence on the population dynamics of the mtDNA variation is very scarce.