Natural DNA precursor pool asymmetry and base sequence context as determinants of replication fidelity

Periodic Variation of Mutation Rates in Bacterial Genomes Associated with Replication Timing

The causes and consequences of spatiotemporal variation in mutation rates remain to be explored in nearly all organisms. Here we examine relationships between local mutation rates and replication timing in three bacterial species whose genomes have multiple chromosomes: Vibrio fischeri, Vibrio cholerae, and Burkholderia cenocepacia. Following five mutation accumulation experiments with these bacteria conducted in the near absence of natural selection, the genomes of clones from each lineage were sequenced and analyzed to identify variation in mutation rates and spectra. In lineages lacking mismatch repair, base substitution mutation rates vary in a mirrored wave-like pattern on opposing replichores of the large chromosomes of V. fischeri and V. cholerae, where concurrently replicated regions experience similar base substitution mutation rates. The base substitution mutation rates on the small chromosome are less variable in both species but occur at similar rates to those in the concurrently replicated regions of the large chromosome. Neither nucleotide composition nor frequency of nucleotide motifs differed among regions experiencing high and low base substitution rates, which along with the inferred ~800-kb wave period suggests that the source of the periodicity is not sequence specific but rather a systematic process related to the cell cycle. These results support the notion that base substitution mutation rates are likely to vary systematically across many bacterial genomes, which exposes certain genes to elevated deleterious mutational load.

That mutation rates vary within bacterial genomes is well known, but the detailed study of these biases has been made possible only recently with contemporary sequencing methods. We applied these methods to understand how bacterial genomes with multiple chromosomes, like those of Vibrio and Burkholderia, might experience heterogeneous mutation rates because of their unusual replication and the greater genetic diversity found on smaller chromosomes. This study captured thousands of mutations and revealed wave-like rate variation that is synchronized with replication timing and not explained by sequence context. The scale of this rate variation over hundreds of kilobases of DNA strongly suggests that a temporally regulated cellular process may generate wave-like variation in mutation risk. These findings add to our understanding of how mutation risk is distributed across bacterial and likely also eukaryotic genomes, owing to their highly conserved replication and repair machinery.

Statistical Methods for Identifying Sequence Motifs Affecting Point Mutations

Mutation processes differ between types of point mutation, genomic locations, cells, and biological species. For some point mutations, specific neighboring bases are known to be mechanistically influential. Beyond these cases, numerous questions remain unresolved, including: what are the sequence motifs that affect point mutations? How large are the motifs? Are they strand symmetric? And, do they vary between samples? We present new log-linear models that allow explicit examination of these questions, along with sequence logo style visualization to enable identifying specific motifs. We demonstrate the performance of these methods by analyzing mutation processes in human germline and malignant melanoma. We recapitulate the known CpG effect, and identify novel motifs, including a highly significant motif associated with A[Formula: see text]G mutations. We show that major effects of neighbors on germline mutation lie within [Formula: see text] of the mutating base. Models are also presented for contrasting the entire mutation spectra (the distribution of the different point mutations). We show the spectra vary significantly between autosomes and X-chromosome, with a difference in T[Formula: see text]C transition dominating. Analyses of malignant melanoma confirmed reported characteristic features of this cancer, including statistically significant strand asymmetry, and markedly different neighboring influences. The methods we present are made freely available as a Python library https://bitbucket.org/pycogent3/mutationmotif.

Genome-Wide Biases in the Rate and Molecular Spectrum of Spontaneous Mutations in Vibrio cholerae and Vibrio fischeri

The vast diversity in nucleotide composition and architecture among bacterial genomes may be partly explained by inherent biases in the rates and spectra of spontaneous mutations. Bacterial genomes with multiple chromosomes are relatively unusual but some are relevant to human health, none more so than the causative agent of cholera, Vibrio cholerae Here, we present the genome-wide mutation spectra in wild-type and mismatch repair (MMR) defective backgrounds of two Vibrio species, the low-%GC squid symbiont V. fischeri and the pathogen V. cholerae, collected under conditions that greatly minimize the efficiency of natural selection. In apparent contrast to their high diversity in nature, both wild-type V. fischeri and V. cholerae have among the lowest rates for base-substitution mutations (bpsms) and insertion-deletion mutations (indels) that have been measured, below 10^-3/genome/generation. Vibrio fischeri and V. cholerae have distinct mutation spectra, but both are AT-biased and produce a surprising number of multi-nucleotide indels. Furthermore, the loss of a functional MMR system caused the mutation spectra of these species to converge, implying that the MMR system itself contributes to species-specific mutation patterns. Bpsm and indel rates varied among genome regions, but do not explain the more rapid evolutionary rates of genes on chromosome 2, which likely result from weaker purifying selection. More generally, the very low mutation rates of Vibrio species correlate inversely with their immense population sizes and suggest that selection may not only have maximized replication fidelity but also optimized other polygenic traits relative to the constraints of genetic drift.

Pyrimidine Pool Disequilibrium Induced by a Cytidine Deaminase Deficiency Inhibits PARP-1 Activity, Leading to the Under Replication of DNA

Genome stability is jeopardized by imbalances of the dNTP pool; such imbalances affect the rate of fork progression. For example, cytidine deaminase (CDA) deficiency leads to an excess of dCTP, slowing the replication fork. We describe here a novel mechanism by which pyrimidine pool disequilibrium compromises the completion of replication and chromosome segregation: the intracellular accumulation of dCTP inhibits PARP-1 activity. CDA deficiency results in incomplete DNA replication when cells enter mitosis, leading to the formation of ultrafine anaphase bridges between sister-chromatids at "difficult-to-replicate" sites such as centromeres and fragile sites. Using molecular combing, electron microscopy and a sensitive assay involving cell imaging to quantify steady-state PAR levels, we found that DNA replication was unsuccessful due to the partial inhibition of basal PARP-1 activity, rather than slower fork speed. The stimulation of PARP-1 activity in CDA-deficient cells restores replication and, thus, chromosome segregation. Moreover, increasing intracellular dCTP levels generates under-replication-induced sister-chromatid bridges as efficiently as PARP-1 knockdown. These results have direct implications for Bloom syndrome (BS), a rare genetic disease combining susceptibility to cancer and genomic instability. BS results from mutation of the BLM gene, encoding BLM, a RecQ 3'-5' DNA helicase, a deficiency of which leads to CDA downregulation. BS cells thus have a CDA defect, resulting in a high frequency of ultrafine anaphase bridges due entirely to dCTP-dependent PARP-1 inhibition and independent of BLM status. Our study describes previously unknown pathological consequences of the distortion of dNTP pools and reveals an unexpected role for PARP-1 in preventing DNA under-replication and chromosome segregation defects.

The Rate and Molecular Spectrum of Spontaneous Mutations in the GC-Rich Multichromosome Genome of Burkholderia cenocepacia

Spontaneous mutations are ultimately essential for evolutionary change and are also the root cause of many diseases. However, until recently, both biological and technical barriers have prevented detailed analyses of mutation profiles, constraining our understanding of the mutation process to a few model organisms and leaving major gaps in our understanding of the role of genome content and structure on mutation. Here, we present a genome-wide view of the molecular mutation spectrum in Burkholderia cenocepacia, a clinically relevant pathogen with high %GC content and multiple chromosomes. We find that B. cenocepacia has low genome-wide mutation rates with insertion-deletion mutations biased toward deletions, consistent with the idea that deletion pressure reduces prokaryotic genome sizes. Unlike prior studies of other organisms, mutations in B. cenocepacia are not AT biased, which suggests that at least some genomes with high %GC content experience unusual base-substitution mutation pressure. Importantly, we also observe variation in both the rates and spectra of mutations among chromosomes and elevated G:C > T:A transversions in late-replicating regions. Thus, although some patterns of mutation appear to be highly conserved across cellular life, others vary between species and even between chromosomes of the same species, potentially influencing the evolution of nucleotide composition and genome architecture.

Aicardi-Goutières syndrome: a model disease for systemic autoimmunity

Systemic autoimmunity is a complex disease process that results from a loss of immunological tolerance characterized by the inability of the immune system to discriminate self from non-self. In patients with the prototypic autoimmune disease systemic lupus erythematosus (SLE), formation of autoantibodies targeting ubiquitous nuclear antigens and subsequent deposition of immune complexes in the vascular bed induces inflammatory tissue injury that can affect virtually any organ system. Given the extraordinary genetic and phenotypic heterogeneity of SLE, one approach to the genetic dissection of complex SLE is to study monogenic diseases, for which a single gene defect is responsible. Considerable success has been achieved from the analysis of the rare monogenic disorder Aicardi-Goutières syndrome (AGS), an inflammatory encephalopathy that clinically resembles in-utero-acquired viral infection and that also shares features with SLE. Progress in understanding the cellular and molecular functions of the AGS causing genes has revealed novel pathways of the metabolism of intracellular nucleic acids, the major targets of the autoimmune attack in patients with SLE. Induction of autoimmunity initiated by immune recognition of endogenous nucleic acids originating from processes such as DNA replication/repair or endogenous retro-elements represents novel paradigms of SLE pathogenesis. These findings illustrate how investigating rare monogenic diseases can also fuel discoveries that advance our understanding of complex disease. This will not only aid the development of improved tools for SLE diagnosis and disease classification, but also the development of novel targeted therapeutic approaches.

Single-stranded nucleic acids promote SAMHD1 complex formation

SAM domain and HD domain-containing protein 1 (SAMHD1) is a dGTP-dependent triphosphohydrolase that degrades deoxyribonucleoside triphosphates (dNTPs) thereby limiting the intracellular dNTP pool. Mutations in SAMHD1 cause Aicardi-Goutières syndrome (AGS), an inflammatory encephalopathy that mimics congenital viral infection and that phenotypically overlaps with the autoimmune disease systemic lupus erythematosus. Both disorders are characterized by activation of the antiviral cytokine interferon-α initiated by immune recognition of self nucleic acids. Here we provide first direct evidence that SAMHD1 associates with endogenous nucleic acids in situ. Using fluorescence cross-correlation spectroscopy, we demonstrate that SAMHD1 specifically interacts with ssRNA and ssDNA and establish that nucleic acid-binding and formation of SAMHD1 complexes are mutually dependent. Interaction with nucleic acids and complex formation do not require the SAM domain, but are dependent on the HD domain and the C-terminal region of SAMHD1. We finally demonstrate that mutations associated with AGS exhibit both impaired nucleic acid-binding and complex formation implicating that interaction with nucleic acids is an integral aspect of SAMHD1 function.

Why genes evolve faster on secondary chromosomes in bacteria

In bacterial genomes composed of more than one chromosome, one replicon is typically larger, harbors more essential genes than the others, and is considered primary. The greater variability of secondary chromosomes among related taxa has led to the theory that they serve as an accessory genome for specific niches or conditions. By this rationale, purifying selection should be weaker on genes on secondary chromosomes because of their reduced necessity or usage. To test this hypothesis we selected bacterial genomes composed of multiple chromosomes from two genera, Burkholderia and Vibrio, and quantified the evolutionary rates (dN and dS) of all orthologs within each genus. Both evolutionary rate parameters were faster among orthologs found on secondary chromosomes than those on the primary chromosome. Further, in every bacterial genome with multiple chromosomes that we studied, genes on secondary chromosomes exhibited significantly weaker codon usage bias than those on primary chromosomes. Faster evolution and reduced codon bias could in turn result from global effects of chromosome position, as genes on secondary chromosomes experience reduced dosage and expression due to their delayed replication, or selection on specific gene attributes. These alternatives were evaluated using orthologs common to genomes with multiple chromosomes and genomes with single chromosomes. Analysis of these ortholog sets suggested that inherently fast-evolving genes tend to be sorted to secondary chromosomes when they arise; however, prolonged evolution on a secondary chromosome further accelerated substitution rates. In summary, secondary chromosomes in bacteria are evolutionary test beds where genes are weakly preserved and evolve more rapidly, likely because they are used less frequently.

Why many bacteria have multiple chromosomes is largely unknown, but a leading hypothesis is that secondary chromosomes evolved from plasmids and now serve as accessory genomes. We tested a key prediction of this theory that genes on secondary chromosomes should evolve faster because they are under less selective constraint. Indeed, orthologous genes shared within two groups of bacteria (Burkholderia or Vibrio) with multiple chromosomes were less conserved and evolved more rapidly when found on secondary chromosomes. Much of these patterns could stem from the tendency of secondary chromosomes to be replicated later in the cell cycle, which reduces their gene dosage, their potential for expression, and selection for their optimal translation. However, the content of secondary chromosomes appears to be predisposed to evolve faster, because these same genes still evolve more rapidly in single-chromosome genomes. In summary, the evolution of divided genomes therefore appears to allow for the long-term segregation of genome content by their rates of expression and dispensability, placing some genes at increased risk of mutational decay and greater turnover.

DNA precursor metabolism and genomic stability

Intracellular concentrations of the four deoxyribonucleoside triphosphates (dNTPs) are closely regulated, and imbalances in the four dNTP pools have genotoxic consequences. Replication errors leading to mutations can occur, for example, if one dNTP in excess drives formation of a non-Watson-Crick base pair or if it forces replicative DNA chain elongation past a mismatch before DNA polymerase can correct the error by 3' exonuclease proofreading. This review focuses on developments since 1994, when the field was last reviewed comprehensively. Emphasis is placed on the following topics: 1) novel aspects of dNTP pool regulation, 2) dNTP pool asymmetries as mutagenic determinants, 3) dNTP metabolism and hypermutagenesis of retroviral genomes, 4) dNTP metabolism and mutagenesis in the mitochondrial genome, 5) chemical modification of nucleotides as a premutagenic event, 6) relationships between dNTP metabolism, genome stability, aging, and cancer.

Natural selection and the emergence of a mutation phenotype: an update of the evolutionary synthesis considering mechanisms that affect genome variation

Most descriptions of evolution assume that all mutations are completely random with respect to their potential effects on survival. However, much like other phenotypic variations that affect the survival of the descendants, intrinsic variations in the probability, type, and location of genetic change can feel the pressure of natural selection. From site-specific recombination to changes in polymerase fidelity and repair of DNA damage, an organism's gene products affect what genetic changes occur in its genome. Through the action of natural selection on these gene products, potentially favorable mutations can become more probable than random. With examples from variation in bacterial surface proteins to the vertebrate immune response, it is clear that a great deal of genetic change is better than "random" with respect to its potential effect on survival. Indeed, some potentially useful mutations are so probable that they can be viewed as being encoded implicitly in the genome. An updated evolutionary theory includes emergence, under selective pressure, of genomic information that affects the probability of different classes of mutation, with consequences for genome survival.

Deoxyribonucleotide pool imbalance stimulates deletions in HeLa cell mitochondrial DNA

Mitochondrial neurogastrointestinal encephalomyopathy (MNGIE) is an autosomal recessive disorder associated with multiple mutations in mitochondrial DNA, both deletions and point mutations, and mutations in the nuclear gene for thymidine phosphorylase. Spinazzola et al. (Spinazzola, A., Marti, R., Nishino, I., Andreu, A., Naini, A., Tadesse, S., Pela, I., Zammarchi, E., Donati, M., Oliver, J., and Hirano, M. (2001) J. Biol. Chem. 277, 4128-4133) showed that MNGIE patients have elevated circulating thymidine levels and they hypothesized that this generates imbalanced mitochondrial deoxyribonucleoside triphosphate (dNTP) pools, which in turn are responsible for mitochondrial (mt) DNA mutagenesis. We tested this hypothesis by culturing HeLa cells in medium supplemented with 50 microM thymidine. After 8-month growth, mtDNA in the thymidine-treated culture, but not the control, showed multiple deletions, as detected both by Southern blotting and by long extension polymerase chain reaction. After 4-h growth in thymidine-supplemented medium, we found the mitochondrial dTTP and dGTP pools to expand significantly, the dCTP pool to drop significantly, and the dATP pool to drop slightly. In whole-cell extracts, dTTP and dGTP pools also expanded, but somewhat less than in mitochondria. The dCTP pool shrank by about 50%, and the dATP pool was essentially unchanged. These results are discussed in terms of the recent report by Nishigaki et al. (Nishigaki, Y., Marti, R., Copeland, W. C., and Hirano, M. (2003) J. Clin. Invest. 111, 1913-1921) that most mitochondrial point mutations in MNGIE patients involve T --> C transitions in sequences containing two As to the 5' side of a T residue. Our finding of dTTP and dGTP elevations and dATP depletion in mitochondrial dNTP pools are consistent with a mutagenic mechanism involving T-G mispairing followed by a next-nucleotide effect involving T insertion opposite A.

Assessing the metabolic function of the MutT 8-oxodeoxyguanosine triphosphatase in Escherichia coli by nucleotide pool analysis

In Escherichia coli the mutT gene is one of several that acts to minimize mutagenesis by reactive oxygen species. The bacterial MutT protein and its mammalian homolog have been shown to catalyze in vitro the hydrolysis of the oxidized deoxyguanosine nucleotide, 8-oxo-dGTP, to its corresponding monophosphate. Thus, the protein is thought to "sanitize" the nucleotide pool by ridding the cell of a nucleotide whose incorporation into DNA would be intensely mutagenic. However, because others have shown mutT mutations to be mutagenic under some conditions of anaerobic growth, and have shown 8-oxo-dGTP to be a poor DNA polymerase substrate, there is reason to question this model. We have devised an assay for 8-oxo-dGTP in bacterial extracts. Using this assay, which involves reversed-phase high-performance liquid chromatography and electrochemical detection, we have been unable to detect 8-oxo-dGTP in extracts of three different mutT mutants of E. coli, even after growth of the bacteria in the presence of hydrogen peroxide. Our estimated upper limit for 8-oxo-dGTP content of these bacteria is about 200 molecules/cell, corresponding to a concentration of about 0.34 microm. When 8-oxo-dGTP was added at 0.34 microm to an in vitro DNA replication system primed with a DNA template that permits scoring of replication errors and with the four normal dNTPs at their estimated intracellular concentrations, there was no detectable effect upon the frequency of replication errors. These findings lead us to question the conclusion that 8-oxo-dGTP is the most significant physiological substrate for the MutT protein.

Effects of biological DNA precursor pool asymmetry upon accuracy of DNA replication in vitro

Deoxyguanosine triphosphate is underrepresented among the four common deoxyribonucleoside triphosphates (dNTPs), typically accounting for just 5-10% of the total dNTP pool. We have asked whether this pool asymmetry affects the fidelity of DNA replication, by use of an in vitro assay in which an M13 phagemid containing the Escherichia coli lacZalpha gene and an SV40 replication origin is replicated by extracts of human cells. By monitoring reversion of either a TGA or TAA codon within the lacZalpha gene, we found that replication in "biologically biased" dNTPs, representing our estimate of the concentrations in HeLa cell nuclei, is not significantly more accurate than when measured in reaction mixtures containing the four dNTPs at equimolar concentrations. However, sequence analysis of revertants revealed significantly different patterns of mispairing events leading to mutation. During replication at biased dNTP levels, mutations at the site 5' to C in the template strand for the TGA triplet were less frequent than seen in equimolar reaction mixtures, suggesting that extension from mismatches at this site is relatively slow, and proofreading efficiency high, when dGTP is the next nucleotide to be incorporated. Mismatches opposite template C, which might have been favored by the low physiological concentrations of dGTP, were not favored in our in vitro system, although one particular substitution at this site, TGA-->TTA, was strongly favored at low [dGTP]. An excess of one dNTP was found in our system to be more mutagenic than a corresponding deficiency. We also estimated dNTP concentrations in non-transformed human fibroblasts and found that in vitro replication at these levels caused significantly fewer mutations than we observed under equimolar conditions (100 microM each dNTP). This increased replication fidelity may result from increased proofreading efficiency at the lower dNTP levels; however, replication rates were decreased only slightly at these non-transformed fibroblast concentrations.

A biokinetic model to describe consequences of inhibition/stimulation in DNA-proofreading and repair-1. Development of the model

A biokinetic model is described which deals with the mathematical consequences of the inhibition or stimulation of DNA proofreading. It demonstrates the development of the number of DNA mismatch-dependent cells (e.g. cells with a malignant phenotype), where such mismatches arise by the in situ interaction of various substances with nucleotides of the DNA. The model can test for consequences by a logic gating on an "if-then" type of analysis in relation to the separate and consecutive processes of proofreading and repair. In particular, the consequences are considered in cases where either (i) the efficacy of proofreading and repair are reduced/prevented (inhibited) or (ii) are increased by some form of stimulation. On the chosen kinetic parameters, the model is accessible to manipulation as new data arising from further investigations become available and are introduced. The model is based on recently published data which show that an increased "mutant fraction" (see note on terms) arises in DNA replication when intracellular nucleotide pools show "asymmetries" (see note on terms). Extraordinarily high mutant fractions can be predicted/have been recorded in the presence of proofreading inhibitors. The model expresses data in mathematical terms of the competition between the development of mismatch-dependent cells and those with authentic genetic information. (Feedback and metastasis-effects and those of wild-type replicates are included.) A computerized (numerical) integration of the corresponding set of differential equations is offered. (A diskette with the program CANCER.xls is available upon request.)

Deoxyribonucleoside triphosphate pool imbalances in vivo are associated with an increased retroviral mutation rate

Deoxyribonucleoside triphosphate (dNTP) pool imbalances are associated with an increase in the rate of misincorporation and hypermutation during in vitro reverse transcription reactions. However, the effects of in vivo dNTP pool imbalances on the accuracy of reverse transcription are unknown. We sought to determine the effects of in vivo dNTP pool imbalances on retroviral mutation rates and to test our hypothesis that 3'-azido-3'-deoxythymidine (AZT) increases the retroviral mutation rates through induction of dNTP pool imbalances. D17 cells were treated with thymidine, hydroxyurea (HU), or AZT, and the effects on in vivo dNTP pools were measured. Thymidine and HU treatments induced significant dNTP pool imbalances. In contrast, AZT treatment had very little effect on the dNTP pools. The effects of in vivo dNTP pool imbalances induced by thymidine and HU treatments on the retroviral mutation rates were also determined. Spleen necrosis virus (SNV)-based and murine leukemia virus (MLV)-based retroviral vectors that expressed the lacZ mutant reporter gene were used. The frequencies of inactivating mutations introduced in the lacZ gene in a single replication cycle provided a measure of the retroviral mutation rates. Treatment of D17 target cells with 500 microM thymidine increased the SNV and MLV mutant frequencies 4.7- and 4-fold, respectively. Treatment of D17 target cells with 2 mM HU increased the SNV and MLV mutant frequencies 2.1- and 2.7-fold, respectively. These results demonstrate that dNTP pool imbalances are associated with an increase in the in vivo retroviral mutation rates, but AZT treatment results in an increase in the retroviral mutation rates by a mechanism not involving alterations in dNTP pools.