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.     

Estimation of DNA Sequence Context-dependent Mutation Rates Using Primate Genomic Sequences

It is understood that DNA and amino acid substitution rates are highly sequence context-dependent, e.g., C --> T substitutions in vertebrates may occur much more frequently at CpG sites and that cysteine substitution rates may depend on support of the context for participation in a disulfide bond. Furthermore, many applications rely on quantitative models of nucleotide or amino acid substitution, including phylogenetic inference and identification of amino acid sequence positions involved in functional specificity. We describe quantification of the context dependence of nucleotide substitution rates using baboon, chimpanzee, and human genomic sequence data generated by the NISC Comparative Sequencing Program. Relative mutation rates are reported for the 96 classes of mutations of the form 5' alphabetagamma 3' --> 5' alphadeltagamma 3', where alpha, beta, gamma, and delta are nucleotides and beta not equal delta, based on maximum likelihood calculations. Our results confirm that C --> T substitutions are enhanced at CpG sites compared with other transitions, relatively independent of the identity of the preceding nucleotide. While, as expected, transitions generally occur more frequently than transversions, we find that the most frequent transversions involve the C at CpG sites (CpG transversions) and that their rate is comparable to the rate of transitions at non-CpG sites. A four-class model of the rates of context-dependent evolution of primate DNA sequences, CpG transitions > non-CpG transitions approximately CpG transversions > non-CpG transversions, captures qualitative features of the mutation spectrum. We find that despite qualitative similarity of mutation rates among different genomic regions, there are statistically significant differences.     

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     

The effect of the nucleotide context on induction of mutations in hypervariable segment 1 of the human mitochondrial DNA

The distribution of unstable nucleotide positions with a higher frequency of homoplastic mutations was analyzed in hypervariable segment 1 (HVS1) of the major noncoding region of human mtDNA. Three motifs (GTAC, ACCC, CCTC) proved to be associated with a higher rate of point substitutions at unstable positions. The motifs were often arranged in direct, including tandem, repeats. Motifs CCTC and ACCC were found in extended poly(C) tracts, which form direct repeats associated with deletions and tandem duplications. The results suggested that the inconstancy of the human mitochondrial genome is to a great extent determined by context-dependent mutations.     

The Effect of the Nucleotide Context on Induction of Mutations in Hypervariable Segment 1 of Human Mitochondrial DNA

The distribution of unstable nucleotide positions with a higher frequency of homoplastic mutations was analyzed in hypervariable segment 1 (HVS1) of the major noncoding region of human mtDNA. Three motifs (GTAC, ACCC, CCTC) proved to be associated with a higher rate of point substitutions at unstable positions. The motifs were often arranged in direct, including tandem, repeats. Motifs CCTC and ACCC were found in extended poly(C) tracts, which form direct repeats associated with deletions and tandem duplications. The results suggested that the inconstancy of the human mitochondrial genome is to a great extent determined by context-dependent mutations.     

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.     

Comparison of the Mutation Spectrum of Hypervariable Segment 1 for Phylogeographical Groups of Human Mitochondrial DNA

The nucleotide sequence variation of hypervariable segment 1 (HVS1) was analyzed for mtDNAs of 88 phylogeographical clusters characteristic of African, West Eurasian, or East Eurasian populations. A significant difference in the distribution of mutations was revealed for the mitochondrial gene pools of the regional human populations. The HVS1 positions were identified whose instability is explained by strand displacement during mtDNA replication. Strand displacement was assumed to be a major mechanism of context-dependent mutagenesis associated with the regional differentiation of human populations.     

Nucleotide correspondence between protein-coding sequences of Helicobacter pylori 26695 and J99 strains

Comparison of open-reading frames (ORFs) H. pylori 26695 and J99 strains has been revealed prevalence of nucleotide replacements as transitions (more than 3%) above transversions (less than 1%). Prevalence of nucleotide transitions is caused by high speed of C : G to T : A transitions in a coding strand of DNA (3.5-5.3%) and not coding strand (2.9-3.9%). The correspondence rate of transversion (A --> C, A --> T, C --> A, C --> G, G --> C, G --> T, T --> A and T --> G) did not exceed 0.84%. The highest correspondence frequency between C and T was detected in ACGT-ATGT (28.3%) - the site of methylation by active methyltransferase M.Hpy99XI in H. pylori 26695 and J99. Thus one can speculate that predominant transition taking place in H. pylori is mutation of C into T, which is realized through cytosine methylation-deamination mechanism.     

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.     

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.     

Spectrum of N4-aminocytidine mutagenesis

N4-Aminocytidine, a nucleoside analog, is a potent mutagen towards phages, bacteria, Drosophila and mammalian cells in culture. In vitro, biochemical studies indicate that this reagent acts by being incorporated into DNA. To elucidate the mechanism of N4-aminocytidine mutagenesis, it is essential to identify the nature of DNA sequence alterations taking place during the mutagenesis. We have analyzed the nucleotide sequence changes in the lac promoter-lacZ alpha region of M13mp2 phage induced by treatment of phage-infected Escherichia coli with N4-aminocytidine. The sequence alterations of DNA samples from 89 mutants of the phage were determined. These mutants had single point mutations, except one mutant, in which a double point mutation was detected. Several hot spots were found: however, there are no apparent relations to particular DNA sequences regarding the locations of these spots. All the mutations are transitions; neither transversions nor deletions/insertions were found. A feature in these transitions is that the A/T to G/C and G/C to A/T changes occur at approximately equal rates. The overall picture of the mutagenesis is consistent with a scheme in which misincorporation and misreplication caused by the modified cytosine structure are the key steps in the DNA replication leading to transitions. Similar nucleotide alterations were found for the mutagenesis induced by an alkylated derivative, N'-methyl-N4-aminocytidine. N4-Aminocytidine also induced reversions of these mutants; both A/T to G/C and G/C to A/T transitions again took place.     

Context dependence, ancestral misidentification, and spurious signatures of natural selection

Population genetic analyses often use polymorphism data from one species, and orthologous genomic sequences from closely related outgroup species. These outgroup sequences are frequently used to identify ancestral alleles at segregating sites and to compare the patterns of polymorphism and divergence. Inherent in such studies is the assumption of parsimony, which posits that the ancestral state of each single nucleotide polymorphism (SNP) is the allele that matches the orthologous site in the outgroup sequence, and that all nucleotide substitutions between species have been observed. This study tests the effect of violating the parsimony assumption when mutation rates vary across sites and over time. Using a context-dependent mutation model that accounts for elevated mutation rates at CpG dinucleotides, increased propensity for transitional versus transversional mutations, as well as other directional and contextual mutation biases estimated along the human lineage, we show (using both simulations and a theoretical model) that enough unobserved substitutions could have occurred since the divergence of human and chimpanzee to cause many statistical tests to spuriously reject neutrality. Moreover, using both the chimpanzee and rhesus macaque genomes to parsimoniously identify ancestral states causes a large fraction of the data to be removed while not completely alleviating problem. By constructing a novel model of the context-dependent mutation process, we can correct polymorphism data for the effect of ancestral misidentification using a single outgroup.     

Analysis of spontaneous base substitutions generated in mutator strains of Bacillus subtilis

In the current studies, we investigated base substitutions in the Bacillus subtilis mutT, mutM, and mutY DNA error-prevention system. In the wild type strain, spontaneous mutations were mainly transitions, either G:C --> A:T or A:T --> G:C. Although both transitions and transversions were observed in mutY and mutM mutants, mutM/mutY double mutants contain strictly G:C --> T:A transversions. In the mutT strain, A:T --> C:G transversion was not observed, and over-expression of the B. subtilis mutT gene had no effect on the mutation rate in the Escherichia coli mutT strain. Using 8-oxo-dGTP-induced mutagenesis, transitions especially A:T --> G:C were predominant in the wild type and mutY strains. In contrary, transversion was high on mutY and double mutant (mutM mutY). Finally, the opuBC and yitG genes were identified from the B. subtilis chromosome as mutator genes that prevented the transition base substitutions.     

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.     

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.     

C --> T mutagenesis and gamma-radiation sensitivity due to deficiency in the Smug1 and Ung DNA glycosylases

The most common genetic change in aerobic organisms is a C:G to T:A mutation. C --> T transitions can arise through spontaneous hydrolytic deamination of cytosine to give a miscoding uracil residue. This is also a frequent DNA lesion induced by oxidative damage, through exposure to agents such as ionizing radiation, or from endogenous sources that are implicated in the aetiology of degenerative diseases, ageing and cancer. The Ung and Smug1 enzymes excise uracil from DNA to effect repair in mammalian cells, and gene-targeted Ung(-/-) mice exhibit a moderate increase in genome-wide spontaneous mutagenesis. Here, we report that stable siRNA-mediated silencing of Smug1 in mouse embryo fibroblasts also generates a mutator phenotype. However, an additive 10-fold increase in spontaneous C:G to T:A transitions in cells deficient in both Smug1 and Ung demonstrates that these enzymes have distinct and nonredundant roles in suppressing C --> T mutability at non-CpG sites. Such cells are also hypersensitive to ionizing radiation, and reveal a role of Smug1 in the repair of lesions generated by oxidation of cytosine.     

Mutation spectra induced by 1-nitropyrene 4,5-oxide and 1-nitropyrene 9,10-oxide in the supF gene of human XP-A fibroblasts

1-Nitropyrene 4,5-oxide and 1-nitropyrene 9,10-oxide are oxidative metabolites that are responsible for the mutagenicity of 1-nitropyrene. In this study, the mutation spectra induced by oxidative metabolites in human cells were determined using a shuttle vector assay. The mutation frequencies induced by 1-nitropyrene 9,10-oxide were 2-3 times higher than those induced by 1-nitropyrene 4,5-oxide. The base substitutions induced by 1-nitropyrene 4,5-oxide were G --> A transitions, G --> C transversions, and G --> T transversions. In the case of 1-nitropyrene 9,10-oxide, G --> A transitions, G --> T transversions, A --> G transitions and G --> C transversions were observed. Most base substitution mutations induced by oxidative metabolites occurred at the guanine sites in the supF gene. These sequence-specific hot spots were commonly identified as 5'-GA sequences for both metabolites. On the other hand, the sequence-specific hot spots at the adenine sites were identified as 5'-CAC sequences for 1-nitropyrene 9,10-oxide. These results suggest that the oxidative metabolites of 1-nitropyrene induce sequence-specific DNA mutations at the guanine and adenine sites at high frequency.     

Mutagenic properties of the 8-amino-2'-deoxyguanosine DNA adduct in mammalian cells.

The DNA adduct 8-amino-2'-deoxyguanosine (8-amino-dG) is found in liver DNA of rats treated with the hepatocarcinogen 2-nitropropane. Site-specifically modified oligodeoxynucleotides were used to explore the mutagenic potential of 8-amino-dG in simian kidney (COS-7) cells. Oligodeoxynucleotides (5'-TCCTCCTX1G2CCTCTC and 5'-TCCTCCTG1X2CCTCTC, X = dG or 8-amino-dG) with the lesion positioned at codon 60 or 61 of the non-coding strand of the human c-Ha- ras1 gene were inserted into single-stranded phagemid vectors and transfected into COS-7 cells. The progeny plasmid obtained was used to transform Escherichia coli DH10B. Transformants were analyzed by oligodeoxynucleotide hybridization and DNA sequencing to establish the mutation frequency and spectrum produced by the modified base. The correct base, dCMP, was incorporated preferentially opposite 8-amino-dG at X1and X2. When 8-amino-dG was at X1, targeted GNH2-->T transversions were detected, along with smaller numbers of GNH2-->A transitions and GNH2-->C transversions. When the adduct was at X2, only GNH2-->T transversions were observed. The frequencies of targeted mutation at X1and X2were 2.7 and 1.7%, respectively. Mutation frequency and mutagenic spectrum were sequence context dependent. In addition, non-targeted G-->T transversions, accompanied by some G-->A transitions, were detected 5' to 8-amino-dG when the lesion was at X2. We conclude that 8-amino-dG is a mutagenic lesion, generating G-->T and G-->C transversions and G-->A transitions in mammalian cells.