Study of the DNA primary structure effect 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.     

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.     

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.     

The mutational spectra of gene p53 in different types of tumors]

The comparative analysis of the frequencies of nucleotide exchanges in mutational spectra of gene p53 (5-8 exons) between germline cancer-prone families (Li-Fraumeni syndrome), between somatic mutations in the tumors of different histogenesis and cell lines, obtained from them, was carried out. The nucleotide positions with high level of mutation events (mutational "hot spots"), typical for germline mutational spectra, tissue-specific patterns of their disappearing and appearing of new ones in solid tumors in vivo, nearly complete absence of "hot spots" in lymphomas and cell lines in vitro were revealed. The obtained results allowed to suggest, that one of the leading factor controlling hot-spots distribution in 5-8 exons of gene p53 is the specificity of cell division conditions in vivo and in vitro.     

Mutational specificity of thymine deprivation-induced mutation in the lacI gene of Escherichia coli

LacI⁻ mutants obtained following 2 and 6 h of thymine deprivation were cloned and sequenced. The mutational spectra recovered were dissimilar. After 2 h of starvation the majority of mutations were base substitutions, largely G: C→C: G transversions. Frameshift mutations but not deletions were observed. In contrast, following 6 h of starvation, with the exception of the G: C→C: G transversion, all possible base substitutions were recovered. Moreover, several deletions but no frameshift events were observed. The differences in the mutational spectra recovered after two periods of thymine deprivation highlight the role of altered nucleotide pools and the potential influence of DNA replication mechanisms.     

Mutational specificity of thymine deprivation-induced mutation in the lacI gene of Escherichia coli

LacI⁻ mutants obtained following 2 and 6 h of thymine deprivation were cloned and sequenced. The mutational spectra recovered were dissimilar. After 2 h of starvation the majority of mutations were base substitutions, largely G: C→C: G transversions. Frameshift mutations but not deletions were observed. In contrast, following 6 h of starvation, with the exception of the G: C→C: G transversion, all possible base substitutions were recovered. Moreover, several deletions but no frameshift events were observed. The differences in the mutational spectra recovered after two periods of thymine deprivation highlight the role of altered nucleotide pools and the potential influence of DNA replication mechanisms.     

The spectra of large second-step mutations are similar for two different mouse autosomes

Loss of tumor suppressor gene expression via mutations plays a critical role in cancer development, particularly when occurring in heterozygous cells. These so-called "second-step" mutational events are often large in size and arise most often from chromosome loss, mitotic recombination, or interstitial deletion. An open question in cancer research is whether different chromosomes are equally susceptible to formation of large mutations, or alternatively if the unique sequence of each chromosome will lead to chromosome-specific mutational spectra. To address this question, the spectra of second-step mutations were determined for chromosomes 8 and 11 in Aprt and Tk mutants, respectively, isolated from primary kidney clones heterozygous for both loci. The results showed that the spectra of large mutational events were essentially the same. This observation suggests that internal and external cellular environments provide the driving force for large autosomal mutational events, and that chromosome structure per se is the substrate upon which these forces act.     

Databases and software for the analysis of mutations in the human p53 gene, the human hprt gene and both the lacI and lacZ gene in transgenic rodents.

We have created databases and software applications for the analysis of DNA mutations at the humanp53gene, the humanhprtgene and both the rodent transgeniclacIandlacZlocus. The databases themselves are stand-alone dBASE files and the software for analysis of the databases runs on IBM-compatible computers. Each database has a separate software analysis program. The software created for these databases permit the filtering, ordering, report generation and display of information in the database. In addition, a significant number of routines have been developed for the analysis of single base substitutions. One method of obtaining the databases and software is via the World Wide Web (WWW). Open the following home page with a Web Browser: http://sunsite.unc.edu/dnam/mainpage.ht ml . Alternatively, the databases and programs are available via public FTP from: anonymous@sunsite.unc.edu . There is no password required to enter the system. The databases and software are found beneath the subdirectory: pub/academic/biology/dna-mutations. Two other programs are available at the site-a program for comparison of mutational spectra and a program for entry of mutational data into a relational database.     

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.     

Databases and software for the analysis of mutations in the human p53 gene, human hprt gene and both the lacI and lacZ gene in transgenic rodents.

We have created databases and software applications for the analysis of DNA mutations at the human p53 gene, the human hprt gene and both the rodent transgenic lacI and lacZ loci. The databases themselves are stand-alone dBASE files and the software for analysis of the databases runs on IBM-compatible computers with Microsoft Windows. Each database has a separate software analysis program. The software created for these databases permit the filtering, ordering, report generation and display of information in the database. In addition, a significant number of routines have been developed for the analysis of single base substitutions. One method of obtaining the databases and software is via the World Wide Web. Open the following home page with a Web Browser: http://sunsite.unc.edu/dnam/mainpage. html . Alternatively, the databases and programs are available via public FTP from: anonymous@sunsite.unc.edu. There is no password required to enter the system. The databases and software are found beneath the subdirectory: pub/academic/biology/dna-mutations. Two other programs are available at the site, a program for comparison of mutational spectra and a program for entry of mutational data into a relational database.     

Pervasive multinucleotide mutational events in eukaryotes

Many aspects of mutational processes are nonrandom, from the preponderance of transitions relative to transversions to the higher rate of mutation at CpG dinucleotides [1]. However, it is still often assumed that single-nucleotide mutations are independent of one another, each being caused by separate mutational events. The occurrence of multiple, closely spaced substitutions appears to violate assumptions of independence and is often interpreted as evidence for the action of adaptive natural selection [2, 3], balancing selection [4], or compensatory evolution [5, 6]. Here we provide evidence of a frequent, widespread multinucleotide mutational process active throughout eukaryotes. Genomic data from mutation-accumulation experiments, parent-offspring trios, and human polymorphisms all show that simultaneous nucleotide substitutions occur within short stretches of DNA. Regardless of species, such multinucleotide mutations (MNMs) consistently comprise ~3% of the total number of nucleotide substitutions. These results imply that previous adaptive interpretations of multiple, closely spaced substitutions may have been unwarranted and that MNMs must be considered when interpreting sequence data.     

Databases and software for the analysis of mutations in the human p53 gene, the human hprt gene and the lacZ gene in transgenic rodents.

We have created databases and software applications for the analysis of DNA mutations in the human p53 gene, the human hprt gene and the rodent transgenic lacZ locus. The databases themselves are stand-alone dBase files and the software for analysis of the databases runs on IBM- compatible computers. The software created for these databases permits filtering, ordering, report generation and display of information in the database. In addition, a significant number of routines have been developed for the analysis of single base substitutions. One method of obtaining the databases and software is via the World Wide Web (WWW). Open home page http://sunsite.unc.edu/dnam/mainpage.ht ml with a WWW browser. Alternatively, the databases and programs are available via public ftp from anonymous@sunsite.unc.edu. There is no password required to enter the system. The databases and software are found in subdirectory pub/academic/biology/dna-mutations. Two other programs are available at the WWW site, a program for comparison of mutational spectra and a program for entry of mutational data into a relational database.     

The role of nucleotide excision repair of Escherichia coli in repair of spontaneous and gamma-radiation-induced DNA damage in the lacZalpha gene

Base excision repair (BER) is a very important repair mechanism to remove oxidative DNA damage. A major oxidative DNA damage after exposure to ionizing radiation is 7,8-dihydro-8-oxoguanine (8oxoG). 8oxoG is a strong mutagenic lesion, which may cause G:C to T:A transversions if not repaired correctly. Formamidopyrimidine-DNA glycosylase (Fpg), a repair enzyme which is part of BER, is the most important enzyme to repair 8oxoG. In the past years, evidence evolved that nucleotide excision repair (NER), a repair system originally thought to repair only bulky DNA lesions, can also repair some oxidative DNA damages. Examples of DNA damages which are recognized by NER are thymine glycol and abasic sites (AP sites). The main objective of this study is to determine if NER can act as a backup system for the repair of spontaneous and gamma-radiation-induced damages when Fpg is deficient. For that purpose, the effect of a NER-deficiency on the spontaneous and gamma-radiation-induced mutation spectrum in the lacZ gene was determined, using double-stranded (ds) M13 DNA, with the lacZalpha gene inserted as mutational target sequence. Subsequently the DNA was transfected into a fpg(-)uvrA(-) Escherichia coli strain (BH420) and the mutational spectra were compared with the spectra of a fpg(-) E. coli strain (BH410) and a wild type E. coli strain (JM105), which were determined in an earlier study. Furthermore, to examine effects which are caused by UvrA-deficiency, and not by Fpg-deficiency, the spontaneous and gamma-radiation-induced mutation spectra of an E. coli strain in which only UvrA is deficient (BH430) were also determined and compared with a wild type E. coli strain (JM105). The results of this study indicate that if only UvrA is deficient, there is an increase in spontaneous G:C to T:A transversions as compared to JM105 and a decrease in A:T to G:C transitions. The gamma-radiation-induced mutation spectrum of BH420 (fpg(-)uvrA(-)) shows a significant decrease in G:C to A:T and G:C to T:A mutations, as compared to BH410 where only Fpg is deficient. Based on these results, we conclude that in our experiments NER is not acting as a backup system if Fpg is deficient. Instead, NER seems to make mistakes, leading to the formation of mutations.     

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.     

Transversion-enriched sequence saturation mutagenesis (SeSaM-Tv+): a random mutagenesis method with consecutive nucleotide exchanges that complements the bias of error-prone PCR

The sequence saturation mutagenesis (SeSaM) method has been advanced to a random mutagenesis method with adjustable mutational biases. SeSaM offers, for example, a bias that is complementary to error-prone (ep) PCR and is enriched in transversions (SeSaM-Tv(+)). dNTP alpha S and three degenerate bases (P, K and I) are used to control mutational bias flexibly. After quantifying incorporation rates of dPTP, dKTP and dITP by terminal transferase using a luciferase-based assay and investigating the read and/or write activities of eight DNA polymerases, a transversion-enriched protocol has been developed. In a mutant library generated using dGTP alpha S and dPTP, transversion frequencies of 16.22-22.58% (G-->T) and 6.38-9.69% (G-->C) were achieved. These mutational spectra are complementary and occur twice as frequently in comparison to standard epPCR methods employing Taq DNA polymerase. For generating more complex mutant libraries, the occurrence of consecutive nucleotide exchanges was increased by 10(5)-10(6)-fold compared to epPCR. Finally, 16.7% of all sequenced mutants contained consecutive nucleotide exchanges composed mainly of a transversion followed by a transition.     

Strange goings-on in the mouse germ line

It is a conventional paradigm that mutagens lead to changes in nucleotide sequence when the cell attempts to repair or replicate lesions in DNA (such as adducts or strand breaks) that have been produced by the mutagens or their metabolites. The resulting changes are located at (or very near) the sites of the initial damage. This is the underlying theory behind mutational spectra work, but how general is it in vivo? Work with ionising radiation has shown that there are interesting things going on in the mouse germ line that do not fall within the conventional paradigm. Mutations occur at certain sites remote from initial DNA damage and in greater than expected number. Bryn Bridges discusses some recent papers on mutational changes in the germ line of mice following exposure to chemical mutagens that suggest that such phenomena may not be confined to radiation.     

Combinatorial codons: a computer program to approximate amino acid probabilities with biased nucleotide usage

Using techniques from optimization theory, we have developed a computer program that approximates a desired probability distribution for amino acids by imposing a probability distribution on the four nucleotides in each of the three codon positions. These base probabilities allow for the generation of biased codons for use in mutational studies and in the design of biologically encoded libraries. The dependencies between codons in the genetic code often makes the exact generation of the desired probability distribution for amino acids impossible. Compromises are often necessary. The program, therefore, not only solves for the "optimal" approximation to the desired distribution (where the definition of "optimal" is influenced by several types of parameters entered by the user), but also solves for a number of "sub-optimal" solutions that are classified into families of similar solutions. A representative of each family is presented to the program user, who can then choose the type of approximation that is best for the intended application. The Combinatorial Codons program is available for use over the web from http://www.wi.mit.edu/kim/computing.html.     

Thermophilic prokaryotes have characteristic patterns of codon usage, amino acid composition and nucleotide content

A number of recent studies have shown that thermophilic prokaryotes have distinguishable patterns of both synonymous codon usage and amino acid composition, indicating the action of natural selection related to thermophily. On the other hand, several other studies of whole genomes have illustrated that nucleotide bias can have dramatic effects on synonymous codon usage and also on the amino acid composition of the encoded proteins. This raises the possibility that the thermophile-specific patterns observed at both the codon and protein levels are merely reflections of a single underlying effect at the level of nucleotide composition. Moreover, such an effect at the nucleotide level might be due entirely to mutational bias. In this study, we have compared the genomes of thermophiles and mesophiles at three levels: nucleotide content, codon usage and amino acid composition. Our results indicate that the genomes of thermophiles are distinguishable from mesophiles at all three levels and that the codon and amino acid frequency differences cannot be explained simply by the patterns of nucleotide composition. At the nucleotide level, we see a consistent tendency for the frequency of adenine to increase at all codon positions within the thermophiles. Thermophiles are also distinguished by their pattern of synonymous codon usage for several amino acids, particularly arginine and isoleucine. At the protein level, the most dramatic effect is a two-fold decrease in the frequency of glutamine residues among thermophiles. These results indicate that adaptation to growth at high temperature requires a coordinated set of evolutionary changes affecting (i) mRNA thermostability, (ii) stability of codon-anticodon interactions and (iii) increased thermostability of the protein products. We conclude that elevated growth temperature imposes selective constraints at all three molecular levels: nucleotide content, codon usage and amino acid composition. In addition to these multiple selective effects, however, the genomes of both thermophiles and mesophiles are often subject to superimposed large changes in composition due to mutational bias.