Exploring the common molecular basis for the universal DNA mutation bias: revival of Löwdin mutation model

Recently, numerous genome analyses revealed the existence of a universal G:C→A:T mutation bias in bacteria, fungi, plants and animals. To explore the molecular basis for this mutation bias, we examined the three well-known DNA mutation models, i.e., oxidative damage model, UV-radiation damage model and CpG hypermutation model. It was revealed that these models cannot provide a sufficient explanation to the universal mutation bias. Therefore, we resorted to a DNA mutation model proposed by L?wdin 40 years ago, which was based on inter-base double proton transfers (DPT). Since DPT is a fundamental and spontaneous chemical process and occurs much more frequently within GC pairs than AT pairs, L?wdin model offers a common explanation for the observed universal mutation bias and thus has broad biological implications.     

Serial analysis of mutation spectra (SAMS): a new approach for the determination of mutation spectra of site-specific DNA damage and their sequence dependence

Many mutations occur as a result of DNA synthesis past the site of DNA damage by DNA damage bypass polymerases. The frequency and types of mutations not only depend on the nature of the damage, but also on the sequence context, as revealed from analysis of mutation spectra of DNA exposed to mutagens. Herein we report a new method for the rapid determination of the effect of sequence context on mutagenesis called SAMS for serial analysis of mutation spectra. This technique makes use of the methodology that underlies serial analysis of gene expression (SAGE) to analyze mutations that result from DNA synthesis past a DNA lesion site-specifically embedded in a library of DNA sequences. To illustrate our technique we determined the effect of sequence context on mutations generated by DNA synthesis past a tetrahydrofuran abasic site model by the DNA damage bypass polymerase yeast polymerase eta.     

Homozygous mutation of MTPAP causes cellular radiosensitivity and persistent DNA double-strand breaks

The study of rare human syndromes characterized by radiosensitivity has been instrumental in identifying novel proteins and pathways involved in DNA damage responses to ionizing radiation. In the present study, a mutation in mitochondrial poly-A-polymerase (MTPAP), not previously recognized for its role in the DNA damage response, was identified by exome sequencing and subsequently associated with cellular radiosensitivity. Cell lines derived from two patients with the homozygous MTPAP missense mutation were radiosensitive, and this radiosensitivity could be abrogated by transfection of wild-type mtPAP cDNA into mtPAP-deficient cell lines. Further analysis of the cellular phenotype revealed delayed DNA repair, increased levels of DNA double-strand breaks, increased reactive oxygen species (ROS), and increased cell death after irradiation (IR). Pre-IR treatment of cells with the potent anti-oxidants, α-lipoic acid and n-acetylcysteine, was sufficient to abrogate the DNA repair and clonogenic survival defects. Our results firmly establish that mutation of the MTPAP gene results in a cellular phenotype of increased DNA damage, reduced repair kinetics, increased cell death by apoptosis, and reduced clonogenic survival after exposure to ionizing radiation, suggesting a pathogenesis that involves the disruption of ROS homeostasis.     

Likelihood-based estimation of microsatellite mutation rates

Microsatellites are widely used in genetic analyses, many of which require reliable estimates of microsatellite mutation rates, yet the factors determining mutation rates are uncertain. The most straightforward and conclusive method by which to study mutation is direct observation of allele transmissions in parent-child pairs, and studies of this type suggest a positive, possibly exponential, relationship between mutation rate and allele size, together with a bias toward length increase. Except for microsatellites on the Y chromosome, however, previous analyses have not made full use of available data and may have introduced bias: mutations have been identified only where child genotypes could not be generated by transmission from parents' genotypes, so that the probability that a mutation is detected depends on the distribution of allele lengths and varies with allele length. We introduce a likelihood-based approach that has two key advantages over existing methods. First, we can make formal comparisons between competing models of microsatellite evolution; second, we obtain asymptotically unbiased and efficient parameter estimates. Application to data composed of 118,866 parent-offspring transmissions of AC microsatellites supports the hypothesis that mutation rate increases exponentially with microsatellite length, with a suggestion that contractions become more likely than expansions as length increases. This would lead to a stationary distribution for allele length maintained by mutational balance. There is no evidence that contractions and expansions differ in their step size distributions.     

Evidence for a predominant role of oxidative damage in germline mutation in mammals

Spontaneous copying errors in replication often are assumed to be the main source of germline mutations in humans and other mammals. However, when laboratory data on context-dependent patterns of oxidative DNA damage are compared with patterns of mutation inferred from mammalian sequence evolution, the strength of the correlation suggests that damage is the main source of mutations. Analysis of damage susceptibility holds promise for improving models of mutational specificity.     

Strong and weak male mutation bias at different sites in the primate genomes: insights from the human-chimpanzee comparison

Male mutation bias is a higher mutation rate in males than in females thought to result from the greater number of germ line cell divisions in males. If errors in DNA replication cause most mutations, then the magnitude of male mutation bias, measured as the male-to-female mutation rate ratio (alpha), should reflect the relative excess of male versus female germ line cell divisions. Evolutionary rates averaged among all sites in a sequence and compared between mammalian sex chromosomes were shown to be indeed higher in males than in females. However, it is presently unknown whether individual classes of substitutions exhibit such bias. To address this issue, we investigated male mutation bias separately at non-CpG and CpG sites using human-chimpanzee whole-genome alignments. We observed strong male mutation bias at non-CpG sites: alpha in the X-autosome comparison was approximately 6-7, which was similar to the male-to-female ratio in the number of germ line cell divisions. In contrast, mutations at CpG sites exhibited weak male mutation bias: alpha in the X-autosome comparison was only approximately 2-3. This is consistent with the methylation-induced and replication-independent mechanism of CpG transitions, which constitute the majority of mutations at CpG sites. Interestingly, our study also indicated weak male mutation bias for transversions at CpG sites, implying a spontaneous mechanism largely not associated with replication. Male mutation bias was equally strong at CpG and non-CpG sites located within unmethylated "CpG islands," suggesting the replication-dependent origin of these mutations. Thus, we found that the strength of male mutation bias is nonuniform in the primate genomes. Importantly, we discovered that male mutation bias depends on the proportion of CpG sites in the loci compared. This might explain the differences in the magnitude of primate male mutation bias observed among studies.     

Why the tautomerization of the G·C Watson–Crick base pair via the DPT does not cause point mutations during DNA replication? QM and QTAIM comprehensive analysis

The ground-state tautomerization of the G·C Watson–Crick base pair by the double proton transfer (DPT) was comprehensively studied in vacuo and in the continuum with a low dielectric constant (??=?4), corresponding to a hydrophobic interface of protein–nucleic acid interactions, using DFT and MP2 levels of quantum-mechanical (QM) theory and quantum theory “Atoms in molecules” (QTAIM). Based on the sweeps of the electron-topological, geometric, polar, and energetic parameters, which describe the course of the G·C???G^?·C^? tautomerization (mutagenic tautomers of the G and C bases are marked with an asterisk) through the DPT along the intrinsic reaction coordinate (IRC), it was proved that it is, strictly speaking, a concerted asynchronous process both at the DFT and MP2 levels of theory, in which protons move with a small time gap in vacuum, while this time delay noticeably increases in the continuum with ??=?4. It was demonstrated using the conductor-like polarizable continuum model (CPCM) that the continuum with ??=?4 does not qualitatively affect the course of the tautomerization reaction. The DPT in the G·C Watson–Crick base pair occurs without any intermediates both in vacuum and in the continuum with ??=?4 at the DFT/MP2 levels of theory. The nine key points along the IRC of the G·C base pair tautomerization, which could be considered as electron-topological “fingerprints” of a concerted asynchronous process of the tautomerization via the DPT, have been identified and fully characterized. These key points have been used to define the reactant, transition state, and product regions of the DPT reaction in the G·C base pair. Analysis of the energetic characteristics of the H-bonds allows us to arrive at a definite conclusion that the middle N1H?N3/N3H?N1 and the lower N2H?O2/N2H?O2 parallel H-bonds in the G·C/G^?·C^? base pairs, respectively, are anticooperative, that is, the strengthening of the middle H-bond is accompanied by the weakening of the lower H-bond. At that point, the upper N4H?O6 and O6H?N4 H-bonds in the G·C and G^?·C^? base pairs, respectively, remain constant at the changes of the middle and the lower H-bonds at the beginning and at the ending of the G·C???G^?·C^? tautomerization. Aiming to answer the question posed in the title of the article, we established that the G^?·C^? Löwdin’s base pair satisfies all the requirements necessary to cause point mutations in DNA except its lifetime, which is much less than the period of time required for the replication machinery to forcibly dissociate a base pair into the monomers (several ns) during DNA replication. So, from the physicochemical point of view, the G^?·C^? Löwdin’s base pair cannot be considered as a source of point mutations arising during DNA replication.     

Dynamics of microsatellite divergence under stepwise mutation and proportional slippage/point mutation models

Recently Kruglyak, Durrett, Schug, and Aquadro showed that microsatellite equilibrium distributions can result from a balance between polymerase slippage and point mutations. Here, we introduce an elaboration of their model that keeps track of all parts of a perfect repeat and a simplification that ignores point mutations. We develop a detailed mathematical theory for these models that exhibits properties of microsatellite distributions, such as positive skewness of allele lengths, that are consistent with data but are inconsistent with the predictions of the stepwise mutation model. We use our theoretical results to analyze the successes and failures of the genetic distances (delta(mu))(2) and D(SW) when used to date four divergences: African vs. non-African human populations, humans vs. chimpanzees, Drosophila melanogaster vs. D. simulans, and sheep vs. cattle. The influence of point mutations explains some of the problems with the last two examples, as does the fact that these genetic distances have large stochastic variance. However, we find that these two features are not enough to explain the problems of dating the human-chimpanzee split. One possible explanation of this phenomenon is that long microsatellites have a mutational bias that favors contractions over expansions.     

Comparative reactivity of mismatched and unpaired bases in relation to their type and surroundings. Chemical cleavage of DNA mismatches in mutation detection analysis

Systematic study of chemical reactivity of non-Watson–Crick base pairs depending on their type and microenvironment was performed on a model system that represents two sets of synthetic DNA duplexes with all types of mismatched and unmatched bases flanked by T·A or G·C pairs. Using comparative cleavage pattern analysis, we identified the main and additional target bases and performed quantitative study of the time course and efficacy of DNA modification caused by potassium permanganate or hydroxylamine. Potassium permanganate in combination with tetraethylammonium chloride was shown to induce DNA cleavage at all mismatched or bulged T residues, as well as at thymines of neighboring canonical pairs. Other mispaired (bulged) bases and thymine residues located on the second position from the mismatch site were not the targets for KMnO₄ attack. In contrast, hydroxylamine cleaved only heteroduplexes containing mismatched or unmatched C residues, and did not modify adjacent cytosines. However when G·C pairs flank bulged C residue, neighboring cytosines are also attacked by hydroxylamine due to defect migration. Chemical reactivity of target bases was shown to correlate strongly with the local disturbance of DNA double helix at mismatch or bulge site. With our model system, we were able to prove the absence of false-negative and false-positive results. Portion of heteroduplex reliably revealed in a mixture with corresponding homoduplex consists of 5% for bulge bases and “open” non-canonical pairs, and 10% for wobble base pairs giving minimal violations in DNA structure. This study provides a complete understanding of the principles of mutation detection methodology based on chemical cleavage of mismatches and clarifies the advantages and limitations of this approach in various biological and conformational studies of DNA.     

Joint inference of microsatellite mutation models, population history and genealogies using transdimensional Markov Chain Monte Carlo

We provide a framework for Bayesian coalescent inference from microsatellite data that enables inference of population history parameters averaged over microsatellite mutation models. To achieve this we first implemented a rich family of microsatellite mutation models and related components in the software package BEAST. BEAST is a powerful tool that performs Bayesian MCMC analysis on molecular data to make coalescent and evolutionary inferences. Our implementation permits the application of existing nonparametric methods to microsatellite data. The implemented microsatellite models are based on the replication slippage mechanism and focus on three properties of microsatellite mutation: length dependency of mutation rate, mutational bias toward expansion or contraction, and number of repeat units changed in a single mutation event. We develop a new model that facilitates microsatellite model averaging and Bayesian model selection by transdimensional MCMC. With Bayesian model averaging, the posterior distributions of population history parameters are integrated across a set of microsatellite models and thus account for model uncertainty. Simulated data are used to evaluate our method in terms of accuracy and precision of estimation and also identification of the true mutation model. Finally we apply our method to a red colobus monkey data set as an example.     

Oxygen accelerates the accumulation of mutations during the senescence and immortalization of murine cells in culture

Oxidative damage is a causal factor in aging and cancer, but it is still not clear how DNA damage, the cellular responses to such damage and its conversion to mutations by misrepair or misreplication contribute to these processes. Using transgenic mice carrying a lacZ mutation reporter, we have previously shown that mutations increase with age in most organs and tissues in vivo . It has also been previously shown that mouse cells respond to oxidative stress, typical of standard culture conditions, by undergoing cellular senescence. To understand better the consequences of oxidative stress, we cultured mouse embryo fibroblasts (MEFs) from lacZ mice under physiological oxygen tension (3%) or the high oxygen tension (20%) associated with standard culture, and determined the frequency and spectrum of mutations. Upon primary culture, the mutation frequency was found to increase approximately three-fold relative to the embryo. The majority of mutations were genome rearrangements. Subsequent culture in 20% oxygen resulted in senescence, followed by spontaneous immortalization. Immortalization was accompanied by an additional three-fold increase in mutations, most of which were G:C to T:A transversions, a signature mutation of oxidative DNA damage. In 3% oxygen, by contrast, MEFs did not senesce and the mutation frequency and spectrum remained similar to primary cultures. These findings demonstrate for the first time the impact of oxidative stress on the genomic integrity of murine cells during senescence and immortalization.     

Beta-thalassemia due to a T----A mutation within the ATA box

Sequence analyses of amplified DNA from a Yugoslavian patient with Hb Lepore-beta-thalassemia and from his father with a simple beta-thalassemia trait have revealed a T----A mutation within the ATA box at a position 30 base pairs upstream from the Cap site. The nucleotide substitution was confirmed through dot-blot analysis of amplified DNA with specific 32P-labeled synthetic oligonucleotide probes. The patient had a clinically severe condition; his Hb Lepore-beta-thalassemia was of the beta + type, as about 8-10% of the non-alpha chain was normal beta A. The same T----A mutation at nucleotide -30 was present on both chromosomes of a young Turkish patient who suffered from a thalassemia intermedia with a low level of Hb F (13.1%) and a relatively high beta A chain synthesis. These data are similar to those obtained for other types of beta +-thalassemia caused by comparable substitutions at positions 31, 29, and 28 base pairs upstream from the Cap site of the beta-globin gene.     

The anomalous effects of biased mutation

A model is presented in which alleles at a number of loci combine to influence the value of a quantitative trait that is subject to stabilizing selection. Mutations can occur to alleles at the loci under consideration. Some of these mutations will tend to increase the value of the trait, while others will tend to decrease it. In contrast to most previous models, we allow the mean effect of mutations to be nonzero. This means that, on average, mutations can have a bias, such that they tend to either increase or decrease the value of the trait. We find, unsurprisingly, that biased mutation moves the equilibrium mean value of the quantitative trait in the direction of the bias. What is more surprising is the behavior of the deviation of the equilibrium mean value of the trait from its optimal value. This has a nonmonotonic dependence on the degree of bias, so that increasing the degree of bias can actually bring the mean phenotype closer to the optimal phenotype. Furthermore, there is a definite maximum to the extent to which biased mutation can cause a difference between the mean phenotype and the optimum. For plausible parameter values, this maximum-possible difference is small. Typically, quantitative-genetics models assume an unconstrained model of mutation, where the expected difference in effect between a parental allele and a mutant allele is independent of the current state of the parental allele. Our results show that models of this sort can easily lead to biologically implausible consequences when mutations are biased. In particular, unconstrained mutation typically leads to a continual increase or decrease in the mean allelic effects at all trait-controlling loci. Thus at each of these loci, the mean allelic effect eventually becomes extreme. This suggests that some of the models of mutation most commonly used in quantitative genetics should be modified so as to introduce genetic constraints.     

Metal-induced point defects in DNA: model and mechanisms

The aim of this work was to study the role of H3O+ and transition-metal (TM) ions in keto-enol and amino-imino tautomeric transitions in DNA base pairs and depurination. In this regard, we discuss the thermodynamic model of ion-DNA interactions and UV display of double-proton transfer (DPT) in GC. The probabilities and energies of rare tautomeric forms of GC pairs in DNA induced by H3O+ and TMwere determined being in the range from0.02 (forMg2+) to 1 ( forCu2+), and from 0 kcal/m (for Cu2+) to 2.3 kcal/m (for Mg2+), respectively. It was shown that 3'ACC5'/5'TGG3' site of DNA double helix, which corresponds to the only triplet 5'UGG3' of RNA that codes the most valuable amino acid tryptophan, is a good target for TM ions to attack. It was also shown that the only way to obtain the tryptophan-coding 5'UGG3' triplet in RNA via transition-type G --> A point mutation caused by TM ions is their interaction with the site of a DNA double helix, which corresponds to 5'CGG3' triplet of RNA that codes arginine.     

A single aspartate mutation in the conserved catalytic site of Rev3L generates a hypomorphic phenotype in vivo and in vitro

Rev3, the catalytic subunit of yeast DNA polymerase ζ, is required for UV resistance and UV-induced mutagenesis, while its mammalian ortholog, REV3L, plays further vital roles in cell proliferation and embryonic development. To assess the contribution of REV3L catalytic activity to its in vivo function, we generated mutant mouse strains in which one or two Ala residues were substituted to the Asp of the invariant catalytic YGDTDS motif. The simultaneous mutation of both Asp (ATA) phenocopies the Rev3l knockout, which proves that the catalytic activity is mandatory for the vital functions of Rev3L, as reported recently. Surprisingly, although the mutation of the first Asp severely impairs the enzymatic activity of other B-family DNA polymerases, the corresponding mutation of Rev3 (ATD) is hypomorphic in yeast and mouse, as it does not affect viability and proliferation and moderately impacts UVC-induced cell death and mutagenesis. Interestingly, Rev3l hypomorphic mutant mice display a distinct, albeit modest, alteration of the immunoglobulin gene mutation spectrum at G-C base pairs, further documenting its role in this process.     

Nitrative and oxidative DNA damage caused by K-ras mutation in mice

Ras mutation is important for carcinogenesis. Carcinogenesis consists of multi-step process with mutations in several genes. We investigated the role of DNA damage in carcinogenesis initiated by K-ras mutation, using conditional transgenic mice. Immunohistochemical analysis revealed that mutagenic 8-nitroguanine and 8-oxo-7,8-dihydro-2′-deoxyguanosine (8-oxodG) were apparently formed in adenocarcinoma caused by mutated K-ras. 8-Nitroguanine was co-localized with iNOS, eNOS, NF-κB, IKK, MAPK, MEK, and mutated K-ras, suggesting that oncogenic K-ras causes additional DNA damage via signaling pathway involving these molecules. It is noteworthy that K-ras mutation mediates not only cell over-proliferation but also the accumulation of mutagenic DNA lesions, leading to carcinogenesis.     

Sperm competition can drive a male-biased mutation rate

A pattern of male-biased mutation has been found in a wide range of species. The standard explanation for this bias is that there are greater numbers of mitotic cell divisions in the history of the average sperm, compared to the average egg, and that mutations typically result from errors made during replication. However, this fails to provide an ultimate evolutionary explanation for why the male germline would tolerate more mutations that are typically deleterious. One possibility is that if there is a tradeoff between producing large numbers of sperm and expending energetic resources in maintaining a lower mutation rate, sperm competition would select for males that produce larger numbers of sperm despite a higher resulting mutation rate. Here I describe a model that jointly considers the fitness consequences of deleterious mutation and mating success in the face of sperm competition. I show that a moderate level of sperm competition can account for the observation that the male germline tolerates a higher mutation rate than the female germline.