Origins of human mitochondrial point mutations as DNA polymerase gamma-mediated errors

Mitochondrial mutational spectra in human cells, tissues and derived tumors for bp 10,030-10,130 are essentially identical, suggesting a predominant mutagenic role for endogenous processes. We hypothesized that errors mediated by mitochondrial DNA polymerase gamma were the primary sources of mutations. Point mutations created in this sequence by human DNA pol gamma in vitro were thus compared to the eighteen mutational hotspots, all single base substitutions, previously found in human tissues. The set of concordant hotspots accounted for 83% of these in vivo mutational events. About half of these mutations are insensitive to prolonged heating of DNA during PCR and half increase proportionally with heating time at 98 degrees C. Primary misincorporation errors and miscopying errors past thermal denaturing products such as deaminated cytosines (uracils) thus appear to be of approximately equal importance. For the sequence studied, these data support the conclusion that, endogenous error mediated by DNA pol gamma constitutes the primary source of mitochondrial point mutations in human tissues.     

Similarity of mutation spectra of the mitochondrial DNA hypervariable segment 1 in Homo and Pan species

The mutation spectrum of mtDNA hypervariable segment 1 (HVS1) was compared for east chimpanzee Pan troglodytes schweigfurthi and human. The two HVS1 had much the same nucleotide composition, and their mutation spectra were similar in major characteristics (substantial prevalence of transitions over transversions, pyrimidine transitions over purine ones, and C --> T over T --> C). DNA strand displacement (dislocation) during replication was identified as a major mechanism of context-dependent mutagenesis in human and chimpanzee mtDNAs. Nucleotide positions with mutations fitting the model of dislocation mutagenesis accounted for 21% of all variable positions in the chimpanzee HVS1. Variable motifs proved to be similar in the chimpanzee and human HVS1. Comparison of the Neanderthal and modern human HVS1 nucleotide sequences showed that most variable nucleotides are in DNA sites allowing context-dependent mutagenesis.     

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.     

Similarity of Mutation Spectra of the Mitochondrial DNA Hypervariable Segment 1 in Homo and Pan Species

The mutation spectrum of mtDNA hypervariable segment 1 (HVS1) was compared for east chimpanzee Pan troglodytes schweigfurthi and human. The two HVS1 had much the same nucleotide composition, and their mutation spectra were similar in major characteristics (substantial prevalence of transitions over transversions, pyrimidine transitions over purine ones, and C T over T C). DNA strand displacement (dislocation) during replication was identified as a major mechanism of context-dependent mutagenesis in human and chimpanzee mtDNAs. Nucleotide positions with mutations fitting the model of dislocation mutagenesis accounted for 21% of all variable positions in the chimpanzee HVS1. Variable motifs proved to be similar in the chimpanzee and human HVS1. Comparison of the Neanderthal and modern human HVS1 nucleotide sequences showed that most variable nucleotides are in DNA sites allowing context-dependent mutagenesis.     

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.     

Mutational spectrum of bleomycin in lacZ mouse kidney: a possible model for mutational spectrum of reactive oxygen species

The mutational spectrum of bleomycin was compared with the spontaneous mutational spectrum in lacZ mouse kidney. Mice were treated with four 20 mg/kg of doses of bleomycin over a two-week period, leading to a mutant fraction several times greater than that of controls. The major class of bleomycin-induced mutations consisted of small deletions, in particular -1 deletions at AT base pairs and hot spots for deletions at 5'-GTC-3' sequences. Smaller, but significant fractions of GC > AT followed by GC > TA substitutions were also observed. In untreated mice, the major class of mutations consisted of GC > AT substitutions followed by GC > TA mutations, and a much smaller fraction of deletions. Other than the specificity of bleomycin for AT base pairs and the 5'-GTC-3' hotspots, the mutational spectrum of bleomycin in mice is similar to that reported for ionizing radiation. However, bleomycin initially mediates the formation of oxidized DNA via reduction of molecular oxygen, as opposed to the radiolysis of water. In this respect mutagenesis induced by bleomycin may be more similar to that induced by endogenous reactive oxygen species (ROS) than mutagenesis induced by ionizing radiation. If bleomycin-induced mutagenesis is an appropriate model for mutagenesis induced by ROS, then, based on the difference between the mutational spectrum of bleomycin and spontaneous mutagenesis, the latter appears not to result predominantly from ROS, at least in mouse kidney.     

Characterization of TCHQ-induced genotoxicity and mutagenesis using the pSP189 shuttle vector in mammalian cells

Tetrachlorohydroquinone (TCHQ) is a major toxic metabolite of the widely used wood preservative, pentachlorophenol (PCP), and it has also been implicated in PCP genotoxicity. However, the underlying mechanisms of genotoxicity and mutagenesis induced by TCHQ remain unclear. In this study, we examined the genotoxicity of TCHQ by using comet assays to detect DNA breakage and formation of TCHQ-DNA adducts. Then, we further verified the levels of mutagenesis by using the pSP189 shuttle vector in A549 human lung carcinoma cells. We demonstrated that TCHQ causes significant genotoxicity by inducing DNA breakage and forming DNA adducts. Additionally, DNA sequence analysis of the TCHQ-induced mutations revealed that 85.36% were single base substitutions, 9.76% were single base insertions, and 4.88% were large fragment deletions. More than 80% of the base substitutions occurred at G:C base pairs, and the mutations were G:C to C:G, G:C to T:A or G:C to A:T transversions and transitions. The most common types of mutations in A549 cells were G:C to A:T (37.14%) and A:T to C:G transitions (14.29%) and G:C to C:G (34.29%) and G:C to T:A (11.43%) transversions. We identified hotspots at nucleotides 129, 141, and 155 in the supF gene of plasmid pSP189. These mutation hotspots accounted for 63% of all single base substitutions. We conclude that TCHQ induces sequence-specific DNA mutations at high frequencies. Therefore, the safety of using this product would be carefully examined.     

Hydrogen peroxide induces G:C to TA and G:C to C:G transversions in the supF gene of Escherichia coli

A vector plasmid, pZ189, carrying an Escherichia coli supF gene as a target for mutations, was treated with a combination of hydrogen peroxide and Fe³⁺/EDTA complex and propagated in E. coli host cells that had been induced for SOS functions by ultraviolet irradiation. The mutations frequency increased by up to 30-fold over spontaneous background levels with increasing concentrations of hydrogen peroxide. The increase in mutation frequency correlated with an increase in the formation of 8-hydroxydeoxyguanosine in the pZ189 DNA. Sequence analysis of 82 independent supF mutant plasmids revealed that 70 mutants contained base substitutions, with 63 of the 70 involving a G:C base pair, and with G:C→C:G (28 cases) and G:C→T:A (26 cases) transversions predominating. Investigation of the influence of the local DNA sequence on the transversions revealed that the guanine at the center of the triplet 5′-PuGA-3′ was five times more likely to mutate after treatment with hydrogen peroxide than that at the center of 5′PyGN3′. G:C→T:A transversions presumably resulted from mispairing of an altered G (probably 8-hydroxydeoxyguanosine) with deoxyadenosine. The origin of the G:C→C:G transversions may be an as yet unidentified lesion generated by hydrogen peroxide. Mutagenic hotspots for base substitutions were found at positions 133, 160 and 168. Mutation spectra and the positions of mutagenic hotspots, when compared with a previously determined spontaneous mutagenesis spectrum, also provide information on the mechanism of spontaneous mutagenesis.     

Randomness of base substitution mutations induced in the lacI gene of Escherichia coli by ionizing radiation

The lacI system of Escherichia coli provides a method for monitoring mutational events at a large number of sites. Using this system, we have previously determined the mutational spectra for gamma-ray and beta-particle emissions resulting from the decay of tritium. Analysis of these mutational spectra reveals that base substitution mutations induced by ionizing radiation are distributed nearly randomly throughout the lacI gene and include all detectable substitution events. The distribution of ionizing radiation-induced mutagenesis is similar to the low frequency of occurrence mutational events induced by other SOS-dependent mutagens. The lack of an apparent nonrandom or high frequency of occurrence component seen with other SOS-dependent mutagens can be best explained as the result of the random interaction of ionizing radiation with the DNA bases leading to production of a variety of base substitutions.     

Spectrum of spontaneous mutations in a cDNA of the human hprt gene integrated in chromosomal DNA

Summary Altered sequences were determined of 52 independent spontaneous mutations occuring in a cDNA of the human hypoxanthine phosphoribosyltransferase (hprt) gene, which was integrated into chromosomal DNA of the mouse cell as a part of the retroviral shuttle vector. Spontaneous mutations comprised a variety of events: base substitutions, frameshifts, deletions, duplications, and complex mutational events, and were distributed randomly over the coding region of the gene. Frameshifts were the most frequent mutational event (38%), and base substitutions were the next most frequent (25%), followed by deletions (19%). Frameshift and deletion mutations commonly occurred preferentially at sites flanked by short direct repeats. Short inverted repeats were frequently found to be associated with duplication and complex mutational events. Analysis of the sequence alterations in the mutant genes suggests that misalignment mutagenesis represents an important molecular mechanism for the generation of spontaneous mutations in eukaryotic cells.     

Mutations induced by ultraviolet light

The different ultraviolet (UV) wavelength components, UVA (320-400 nm), UVB (280-320 nm), and UVC (200-280 nm), have distinct mutagenic properties. A hallmark of UVC and UVB mutagenesis is the high frequency of transition mutations at dipyrimidine sequences containing cytosine. In human skin cancers, about 35% of all mutations in the p53 gene are transitions at dipyrimidines within the sequence 5'-TCG and 5'-CCG, and these are localized at several mutational hotspots. Since 5'-CG sequences are methylated along the p53 coding sequence in human cells, these mutations may be derived from sunlight-induced pyrimidine dimers forming at sequences that contain 5-methylcytosine. Cyclobutane pyrimidine dimers (CPDs) form preferentially at dipyrimidines containing 5-methylcytosine when cells are irradiated with UVB or sunlight. In order to define the contribution of 5-methylcytosine to sunlight-induced mutations, the lacI and cII transgenes in mouse fibroblasts were used as mutational targets. After 254 nm UVC irradiation, only 6-9% of the base substitutions were at dipyrimidines containing 5-methylcytosine. However, 24-32% of the solar light-induced mutations were at dipyrimidines that contain 5-methylcytosine and most of these mutations were transitions. Thus, CPDs forming preferentially at dipyrimidines with 5-methylcytosine are responsible for a considerable fraction of the mutations induced by sunlight in mammalian cells. Using mouse cell lines harboring photoproduct-specific photolyases and mutational reporter genes, we showed that CPDs (rather than 6-4 photoproducts or other lesions) are responsible for the great majority of UVB-induced mutations. An important component of UVB mutagenesis is the deamination of cytosine and 5-methylcytosine within CPDs. The mutational specificity of long-wave UVA (340-400 nm) is distinct from that of the shorter wavelength UV and is characterized mainly by G to T transversions presumably arising through mechanisms involving oxidized DNA bases. We also discuss the role of DNA damage-tolerant DNA polymerases in UV lesion bypass and mutagenesis.     

Modulation of an ultraviolet mutational hotspot in a shuttle vector Xeroderma cells.

Ultraviolet mutagenesis of the shuttle vector plasmid pZ189 in Xeroderma Pigmentosum cells yields a mutational pattern marked by hotspots at photoproduct sites on both strands of the supF marker gene. In order to test the influence of strand orientation on the appearance of hotspots the mutagenesis study was repeated on a vector with the supF gene in the inverted orientation. We recovered a pattern the same as that in the earlier work and conclude that the nature of the DNA polymerase involved in the replication of specific strands is not a primary determinant of hotspot occurrence in this system. One of the hotspots lies in an 8 base palindrome while the corresponding site on the other strand was not a hotspot. These results were obtained with calcium phosphate transfection of the UV treated vector. When DEAE dextran was used as a transfection agent both sites in the palindrome were hotspots. In a mixing experiment the calcium phosphate pattern was recovered. Our data suggest that the sequence determinants of mutational probability at these two sites lie outside the 8 bases of the palindrome and that mutagenesis at one, but not the other, site is sensitive to perturbation of cellular calcium levels.     

DNA damage and mutations produced by chloroacetaldehyde in a CpG-methylated target gene

Chloroacetaldehyde (CAA) is a metabolite of the human carcinogen vinyl chloride. CAA produces several types of DNA adducts including the exocyclic base adducts 3,N(4)-ethenocytosine, 1,N(6)-ethenoadenine, N(2),3-ethenoguanine, and 1,N(2)-ethenoguanine. Adducts of CAA with 5-methylcytosine have not yet been characterized. Here we have analyzed the mutational spectra produced by CAA in the supF gene of the pSP189 shuttle vector when present in either an unmethylated or CpG-methylated state. The vectors were replicated in human nucleotide excision repair-deficient XP-A fibroblasts. The mutational spectra obtained with the unmethylated and methylated supF target genes were generally similar with a preponderance of C/G to T/A transitions and C/G to A/T transversions. CAA-induced DNA adducts were mapped along the supF gene by using thermostable thymine DNA glycosylase (TDG) in conjunction with ligation-mediated PCR or by a Taq polymerase stop assay. Prominent CAA-induced TDG-sensitive sites were seen at several CpG positions but were independent of methylation. Methylated CpG sites were sites of CAA-induced mutations but were not the major mutational hotspots. Taq polymerase arrest sites were observed at numerous sequence positions in the supF gene and reflected the rather broad distributions of mutations along the sequence. We conclude that methylated CpG sites are not preferential targets for chloroacetaldehyde-induced mutagenesis.     

The distribution of UV damage in the lacI gene of Escherichia coli: correlation with mutation spectrum.

We have determined the UV (254 nm) damage distribution in the transcribed and non-transcribed strands of the i-d region of the Escherichia coli lacI gene. The locations of replication blocking lesions were revealed as termination sites of T7 DNA polymerase and/or T4 DNA polymerase 3'-5' exonuclease. Termination products, i.e. both cyclobutane pyrimidine dimers and 6-4 photoproducts, were resolved and analysed on an automated DNA sequencer. These two major photoproducts are not randomly distributed along the gene, but occur in clusters, in longer runs of pyrimidines. We also have compared the UV damage distribution with the previously reported UV-induced base substitutions in the same region. Mutational hotspots, in both repair-deficient and repair-proficient strains of E. coli, are all located in stretches of pyrimidines, and thus correlate well with the distribution of photolesions. One mutational hotspot in the wild-type strain may reflect the high frequency of closely opposed lesions. To explain the other mutational hotspots, we propose that the repair of UV lesions is impaired due to the local conformation of the DNA, which might deviate from the B-form. This hypothesis is supported by the excess of mutational hotspots in pyrimidine runs in the Uvr+ strain compared to Uvr-. Runs of pyrimidines thus represent both damage- and mutation-prone sequences following UV treatment.     

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.     

Theoretical analysis of mutation hotspots and their DNA sequence context specificity

Mutation frequencies vary significantly along nucleotide sequences such that mutations often concentrate at certain positions called hotspots. Mutation hotspots in DNA reflect intrinsic properties of the mutation process, such as sequence specificity, that manifests itself at the level of interaction between mutagens, DNA, and the action of the repair and replication machineries. The hotspots might also reflect structural and functional features of the respective DNA sequences. When mutations in a gene are identified using a particular experimental system, resulting hotspots could reflect the properties of the gene product and the mutant selection scheme. Analysis of the nucleotide sequence context of hotspots can provide information on the molecular mechanisms of mutagenesis. However, the determinants of mutation frequency and specificity are complex, and there are many analytical methods for their study. Here we review computational approaches for analyzing mutation spectra (distribution of mutations along the target genes) that include many mutable (detectable) positions. The following methods are reviewed: derivation of a consensus sequence, application of regression approaches to correlate nucleotide sequence features with mutation frequency, mutation hotspot prediction, analysis of oligonucleotide composition of regions containing mutations, pairwise comparison of mutation spectra, analysis of multiple spectra, and analysis of "context-free" characteristics. The advantages and pitfalls of these methods are discussed and illustrated by examples from the literature. The most reliable analyses were obtained when several methods were combined and information from theoretical analysis and experimental observations was considered simultaneously. Simple, robust approaches should be used with small samples of mutations, whereas combinations of simple and complex approaches may be required for large samples. We discuss several well-documented studies where analysis of mutation spectra has substantially contributed to the current understanding of molecular mechanisms of mutagenesis. The nucleotide sequence context of mutational hotspots is a fingerprint of interactions between DNA and DNA repair, replication, and modification enzymes, and the analysis of hotspot context provides evidence of such interactions.     

Distinct functional contributions of three potential secondary structures in the phage G4 origin of complementary DNA strand synthesis

Three potential secondary structures, stem-loops I, II, and III, are contained in the phage G4 origin of complementary DNA strand synthesis, G4oric, and are believed to be involved in its recognition by dnaG-encoded primase and the synthesis of primer RNA. In a previous publication [Sakai et al., Gene 71 (1988) 323-330], we suggested that base pairing between the loops of stem-loops I, and II, and/or II and III, might play a role in G4oric function. To test this hypothesis, site-directed mutagenesis was used to construct mutants which carried base substitutions in loops I, II and III that destroyed possible interloop base pairing. These mutations, however, did not seriously affect G4oric activity. This indicates that base pairing between the loops is not essential for G4oric functional activity, and also that base substitutions which do not affect the secondary structure of stem-loops I, II and III, do not affect G4oric activity. To complete an analysis of the effects of altering the structure of the G4oric stem-loops, insertions were made into stem-loop III. In contrast to stem-loops I and II, all insertions into stem-loop III destroyed in vivo G4oric activity.     

Guanine Holes Are Prominent Targets for Mutation in Cancer and Inherited Disease

Single base substitutions constitute the most frequent type of human gene mutation and are a leading cause of cancer and inherited disease. These alterations occur non-randomly in DNA, being strongly influenced by the local nucleotide sequence context. However, the molecular mechanisms underlying such sequence context-dependent mutagenesis are not fully understood. Using bioinformatics, computational and molecular modeling analyses, we have determined the frequencies of mutation at G•C bp in the context of all 64 5′-NGNN-3′ motifs that contain the mutation at the second position. Twenty-four datasets were employed, comprising >530,000 somatic single base substitutions from 21 cancer genomes, >77,000 germline single-base substitutions causing or associated with human inherited disease and 16.7 million benign germline single-nucleotide variants. In several cancer types, the number of mutated motifs correlated both with the free energies of base stacking and the energies required for abstracting an electron from the target guanines (ionization potentials). Similar correlations were also evident for the pathological missense and nonsense germline mutations, but only when the target guanines were located on the non-transcribed DNA strand. Likewise, pathogenic splicing mutations predominantly affected positions in which a purine was located on the non-transcribed DNA strand. Novel candidate driver mutations and tissue-specific mutational patterns were also identified in the cancer datasets. We conclude that electron transfer reactions within the DNA molecule contribute to sequence context-dependent mutagenesis, involving both somatic driver and passenger mutations in cancer, as well as germline alterations causing or associated with inherited disease.