Clusters of mutations from transient hypermutability

Targeted Hypermutation as a Survival Strategy: A Theoretical Approach

Targeted hypermutation has proven to be a useful survival strategy for bacteria under severe stress and is also used by multicellular organisms in specific instances such as the mammalian immune system. This might appear surprising, given the generally observed deleterious effects of poor replication fidelity/high mutation rate. A previous theoretical model designed to explore the role of replication fidelity in the origin of life was applied to a simulated hypermutation scenario. The results confirmed that the same model is useful for analyzing hypermutation and can predict the effects of the same parameters (survival probability, replication fidelity, mutation effect, and others) on the survival of cellular populations undergoing hypermutation as a result of severe stress.

No evidence for widespread positive selection on double substitutions within codons in primates and yeasts

Nucleotide substitutions in protein-coding genes can be divided into synonymous (S) and non-synonymous (N) ones that alter amino acids (including nonsense mutations causing stop codons). The S substitutions are expected to have little effect on function. The N substitutions almost always are affected by strong purifying selection that eliminates them from evolving populations. However, additional mutations of nearby bases can modulate the deleterious effect of single N substitutions and, thus, could be subjected to the positive selection. This effect has been demonstrated for mutations in the serine codons, stop codons and double N substitutions in prokaryotes. In all abovementioned cases, a novel technique was applied that allows elucidating the effects of selection on double substitutions considering mutational biases. Here, we applied the same technique to study double N substitutions in eukaryotic lineages of primates and yeast. We identified markedly fewer cases of purifying selection relative to prokaryotes and no evidence of codon double substitutions under positive selection. This is consistent with previous studies of serine codons in primates and yeast. In general, the obtained results strongly suggest that there are major differences between studied pro- and eukaryotes; double substitutions in primates and yeasts largely reflect mutational biases and are not hallmarks of selection. This is especially important in the context of detection of positive selection in codons because it has been suggested that multiple mutations in codons cause false inferences of lineage-specific site positive selection. It is likely that this concern is applicable to previously studied prokaryotes but not to primates and yeasts where markedly fewer double substitutions are affected by positive selection.

A Comprehensive View of Translesion Synthesis in Escherichia coli

The lesion bypass pathway, translesion synthesis (TLS), exists in essentially all organisms and is considered a pathway for postreplicative gap repair and, at the same time, for lesion tolerance. As with the saying "a trip is not over until you get back home," studying TLS only at the site of the lesion is not enough to understand the whole process of TLS. Recently, a genetic study uncovered that polymerase V (Pol V), a poorly expressed Escherichia coli TLS polymerase, is not only involved in the TLS step per se but also participates in the gap-filling reaction over several hundred nucleotides. The same study revealed that in contrast, Pol IV, another highly expressed TLS polymerase, essentially stays away from the gap-filling reaction. These observations imply fundamentally different ways these polymerases are recruited to DNA in cells. While access of Pol IV appears to be governed by mass action, efficient recruitment of Pol V involves a chaperone-like action of the RecA filament. We present a model of Pol V activation: the 3' tip of the RecA filament initially stabilizes Pol V to allow stable complex formation with a sliding β-clamp, followed by the capture of the terminal RecA monomer by Pol V, thus forming a functional Pol V complex. This activation process likely determines higher accessibility of Pol V than of Pol IV to normal DNA. Finally, we discuss the biological significance of TLS polymerases during gap-filling reactions: error-prone gap-filling synthesis may contribute as a driving force for genetic diversity, adaptive mutation, and evolution.

Nitrogen starvation reveals the mitotic potential of mutants in the S/MAPK pathways

The genetics of quiescence is an emerging field compared to that of growth, yet both states generate spontaneous mutations and genetic diversity fueling evolution. Reconciling mutation rates in dividing conditions and mutation accumulation as a function of time in non-dividing situations remains a challenge. Nitrogen-starved fission yeast cells reversibly arrest proliferation, are metabolically active and highly resistant to a variety of stresses. Here, we show that mutations in stress- and mitogen-activated protein kinase (S/MAPK) signaling pathways are enriched in aging cultures. Targeted resequencing and competition experiments indicate that these mutants arise in the first month of quiescence and expand clonally during the second month at the expense of the parental population. Reconstitution experiments show that S/MAPK modules mediate the sacrifice of many cells for the benefit of some mutants. These findings suggest that non-dividing conditions promote genetic diversity to generate a social cellular environment prone to kin selection.

Nitrogen-starved fission yeast cells survive for weeks without dividing. Here, the authors show that some of these surviving cells accumulate mutations in the stress- and mitogen-activated protein kinase pathways and outcompete their parental cells, which provide nutrients for the mutant cells.

Postzygotic Somatic Mutations in the Human Brain Expand the Threshold-Liability Model of Schizophrenia

The search for what causes schizophrenia has been onerous. This research has included extensive assessment of a variety of genetic and environmental factors using ever emerging high-resolution technologies and traditional understanding of the biology of the brain. These efforts have identified a large number of schizophrenia-associated genes, some of which are altered by mutational and epi-mutational mechanisms in a threshold liability model of schizophrenia development. The results, however, have limited predictability and the actual cause of the disease remains unknown. This current state asks for conceptualizing the problem differently in light of novel insights into the nature of mutations, the biology of the brain and the fine precision and resolution of emerging technologies. There is mounting evidence that mutations acquired during postzygotic development are more common than germline mutations. Also, the postzygotic somatic mutations including epimutations (PZMs), which often lead to somatic mosaicism, are relatively common in the mammalian brain in comparison to most other tissues and PZMs are more common in patients with neurodevelopmental mental disorders, including schizophrenia. Further, previously inaccessible, detection of PZMs is becoming feasible with the advent of novel technologies that include single-cell genomics and epigenomics and the use of exquisite experimental designs including use of monozygotic twins discordant for the disease. These developments allow us to propose a working hypothesis and expand the threshold liability model of schizophrenia that already encompasses familial genetic, epigenetic and environmental factors to include somatic de novo PZMs. Further, we offer a test for this expanded model using currently available genome sequences and methylome data on monozygotic twins discordant for schizophrenia (MZD) and their parents. The results of this analysis argue that PZMs play a significant role in the development of schizophrenia and explain extensive heterogeneity seen across patients. It also offers the potential to convincingly link PZMs to both nervous system health and disease, an area that has remained challenging to study and relatively under explored.

Crossing fitness valleys via double substitutions within codons

BACKGROUND

Single nucleotide substitutions in protein-coding genes can be divided into synonymous (S), with little fitness effect, and non-synonymous (N) ones that alter amino acids and thus generally have a greater effect. Most of the N substitutions are affected by purifying selection that eliminates them from evolving populations. However, additional mutations of nearby bases potentially could alleviate the deleterious effect of single substitutions, making them subject to positive selection. To elucidate the effects of selection on double substitutions in all codons, it is critical to differentiate selection from mutational biases.

RESULTS

We addressed the evolutionary regimes of within-codon double substitutions in 37 groups of closely related prokaryotic genomes from diverse phyla by comparing the fractions of double substitutions within codons to those of the equivalent double S substitutions in adjacent codons. Under the assumption that substitutions occur one at a time, all within-codon double substitutions can be represented as "ancestral-intermediate-final" sequences (where "intermediate" refers to the first single substitution and "final" refers to the second substitution) and can be partitioned into four classes: (1) SS, S intermediate-S final; (2) SN, S intermediate-N final; (3) NS, N intermediate-S final; and (4) NN, N intermediate-N final. We found that the selective pressure on the second substitution markedly differs among these classes of double substitutions. Analogous to single S (synonymous) substitutions, SS double substitutions evolve neutrally, whereas analogous to single N (non-synonymous) substitutions, SN double substitutions are subject to purifying selection. In contrast, NS show positive selection on the second step because the original amino acid is recovered. The NN double substitutions are heterogeneous and can be subject to either purifying or positive selection, or evolve neutrally, depending on the amino acid similarity between the final or intermediate and the ancestral states.

CONCLUSIONS

The results of the present, comprehensive analysis of the evolutionary landscape of within-codon double substitutions reaffirm the largely conservative regime of protein evolution. However, the second step of a double substitution can be subject to positive selection when the first step is deleterious. Such positive selection can result in frequent crossing of valleys on the fitness landscape.

Pol V-Mediated Translesion Synthesis Elicits Localized Untargeted Mutagenesis during Post-replicative Gap Repair

In vivo, replication forks proceed beyond replication-blocking lesions by way of downstream repriming, generating daughter strand gaps that are subsequently processed by post-replicative repair pathways such as homologous recombination and translesion synthesis (TLS). The way these gaps are filled during TLS is presently unknown. The structure of gap repair synthesis was assessed by sequencing large collections of single DNA molecules that underwent specific TLS events in vivo. The higher error frequency of specialized relative to replicative polymerases allowed us to visualize gap-filling events at high resolution. Unexpectedly, the data reveal that a specialized polymerase, Pol V, synthesizes stretches of DNA both upstream and downstream of a site-specific DNA lesion. Pol V-mediated untargeted mutations are thus spread over several hundred nucleotides, strongly eliciting genetic instability on either side of a given lesion. Consequently, post-replicative gap repair may be a source of untargeted mutations critical for gene diversification in adaptation and evolution.

What is mutation? A chapter in the series: How microbes "jeopardize" the modern synthesis

Mutations drive evolution and were assumed to occur by chance: constantly, gradually, roughly uniformly in genomes, and without regard to environmental inputs, but this view is being revised by discoveries of molecular mechanisms of mutation in bacteria, now translated across the tree of life. These mechanisms reveal a picture of highly regulated mutagenesis, up-regulated temporally by stress responses and activated when cells/organisms are maladapted to their environments-when stressed-potentially accelerating adaptation. Mutation is also nonrandom in genomic space, with multiple simultaneous mutations falling in local clusters, which may allow concerted evolution-the multiple changes needed to adapt protein functions and protein machines encoded by linked genes. Molecular mechanisms of stress-inducible mutation change ideas about evolution and suggest different ways to model and address cancer development, infectious disease, and evolution generally.

Genome Plasticity of agr-Defective Staphylococcus aureus during Clinical Infection

Therapy for bacteremia caused by Staphylococcus aureus is often ineffective, even when treatment conditions are optimal according to experimental protocols. Adapted subclones, such as those bearing mutations that attenuate agr-mediated virulence activation, are associated with persistent infection and patient mortality.

Therapy for bacteremia caused by Staphylococcus aureus is often ineffective, even when treatment conditions are optimal according to experimental protocols. Adapted subclones, such as those bearing mutations that attenuate agr-mediated virulence activation, are associated with persistent infection and patient mortality. To identify additional alterations in agr-defective mutants, we sequenced and assembled the complete genomes of clone pairs from colonizing and infected sites of several patients in whom S. aureus demonstrated a within-host loss of agr function. We report that events associated with agr inactivation result in agr-defective blood and nares strain pairs that are enriched in mutations compared to pairs from wild-type controls. The random distribution of mutations between colonizing and infecting strains from the same patient, and between strains from different patients, suggests that much of the genetic complexity of agr-defective strains results from prolonged infection or therapy-induced stress. However, in one of the agr-defective infecting strains, multiple genetic changes resulted in increased virulence in a murine model of bloodstream infection, bypassing the mutation of agr and raising the possibility that some changes were selected. Expression profiling correlated the elevated virulence of this agr-defective mutant to restored expression of the agr-regulated ESAT6-like type VII secretion system, a known virulence factor. Thus, additional mutations outside the agr locus can contribute to diversification and adaptation during infection by S. aureus agr mutants associated with poor patient outcomes.

Real-time dynamics of mutagenesis reveal the chronology of DNA repair and damage tolerance responses in single cells

Evolutionary processes are driven by diverse molecular mechanisms that act in the creation and prevention of mutations. It remains unclear how these mechanisms are regulated because limitations of existing mutation assays have precluded measuring how mutation rates vary over time in single cells. Toward this goal, I detected nascent DNA mismatches as a proxy for mutagenesis and simultaneously followed gene expression dynamics in single Escherichia coli cells using microfluidics. This general microscopy-based approach revealed the real-time dynamics of mutagenesis in response to DNA alkylation damage and antibiotic treatments. It also enabled relating the creation of DNA mismatches to the chronology of the underlying molecular processes. By avoiding population averaging, I discovered cell-to-cell variation in mutagenesis that correlated with heterogeneity in the expression of alternative responses to DNA damage. Pulses of mutagenesis are shown to arise from transient DNA repair deficiency. Constitutive expression of DNA repair pathways and induction of damage tolerance by the SOS response compensate for delays in the activation of inducible DNA repair mechanisms, together providing robustness against the toxic and mutagenic effects of DNA alkylation damage.

Synthetic biology, genome mining, and combinatorial biosynthesis of NRPS-derived antibiotics: a perspective

Combinatorial biosynthesis of novel secondary metabolites derived from nonribosomal peptide synthetases (NRPSs) has been in slow development for about a quarter of a century. Progress has been hampered by the complexity of the giant multimodular multienzymes. More recently, advances have been made on understanding the chemical and structural biology of these complex megaenzymes, and on learning the design rules for engineering functional hybrid enzymes. In this perspective, I address what has been learned about successful engineering of complex lipopeptides related to daptomycin, and discuss how synthetic biology and microbial genome mining can converge to broaden the scope and enhance the speed and robustness of combinatorial biosynthesis of NRPS-derived natural products for drug discovery.

Stress-Induced Mutagenesis: Implications in Cancer and Drug Resistance

Genomic instability underlies many cancers and generates genetic variation that drives cancer initiation, progression, and therapy resistance. In contrast with classical assumptions that mutations occur purely stochastically at constant, gradual rates, microbes, plants, flies, and human cancer cells possess mechanisms of mutagenesis that are upregulated by stress responses. These generate transient, genetic-diversity bursts that can propel evolution, specifically when cells are poorly adapted to their environments-that is, when stressed. We review molecular mechanisms of stress-response-dependent (stress-induced) mutagenesis that occur from bacteria to cancer, and are activated by starvation, drugs, hypoxia, and other stressors. We discuss mutagenic DNA break repair in Escherichia coli as a model for mechanisms in cancers. The temporal regulation of mutagenesis by stress responses and spatial restriction in genomes are common themes across the tree of life. Both can accelerate evolution, including the evolution of cancers. We discuss possible anti-evolvability drugs, aimed at targeting mutagenesis and other variation generators, that could be used to delay the evolution of cancer progression and therapy resistance.

Population Heterogeneity in Mutation Rate Increases the Frequency of Higher-Order Mutants and Reduces Long-Term Mutational Load

Mutation rate is a crucial evolutionary parameter that has typically been treated as a constant in population genetic analyses. However, the propensity to mutate is likely to vary among co-existing individuals within a population, due to genetic polymorphisms, heterogeneous environmental influences, and random physiological fluctuations. We review the evidence for mutation rate heterogeneity and explore its consequences by extending classic population genetic models to allow an arbitrary distribution of mutation rate among individuals, either with or without inheritance. With this general new framework, we rigorously establish the effects of heterogeneity at various evolutionary timescales. In a single generation, variation of mutation rate about the mean increases the probability of producing zero or many simultaneous mutations on a genome. Over multiple generations of mutation and selection, heterogeneity accelerates the appearance of both deleterious and beneficial multi-point mutants. At mutation-selection balance, higher-order mutant frequencies are likewise boosted, while lower-order mutants exhibit subtler effects; nonetheless, population mean fitness is always enhanced. We quantify the dependencies on moments of the mutation rate distribution and selection coefficients, and clarify the role of mutation rate inheritance. While typical methods of estimating mutation rate will recover only the population mean, analyses assuming mutation rate is fixed to this mean could underestimate the potential for multi-locus adaptation, including medically relevant evolution in pathogenic and cancerous populations. We discuss the potential to empirically parameterize mutation rate distributions, which have to date hardly been quantified.

Bioinformatics and Drug Discovery

Bioinformatic analysis can not only accelerate drug target identification and drug candidate screening and refinement, but also facilitate characterization of side effects and predict drug resistance. High-throughput data such as genomic, epigenetic, genome architecture, cistromic, transcriptomic, proteomic, and ribosome profiling data have all made significant contribution to mechanismbased drug discovery and drug repurposing. Accumulation of protein and RNA structures, as well as development of homology modeling and protein structure simulation, coupled with large structure databases of small molecules and metabolites, paved the way for more realistic protein-ligand docking experiments and more informative virtual screening. I present the conceptual framework that drives the collection of these high-throughput data, summarize the utility and potential of mining these data in drug discovery, outline a few inherent limitations in data and software mining these data, point out news ways to refine analysis of these diverse types of data, and highlight commonly used software and databases relevant to drug discovery.

Evolutionary switches between two serine codon sets are driven by selection

Serine is the only amino acid that is encoded by two disjoint codon sets so that a tandem substitution of two nucleotides is required to switch between the two sets. Previously published evidence suggests that, for the most evolutionarily conserved serines, the codon set switch occurs by simultaneous substitution of two nucleotides. Here we report a genome-wide reconstruction of the evolution of serine codons in triplets of closely related species from diverse prokaryotes and eukaryotes. The results indicate that the great majority of codon set switches proceed by two consecutive nucleotide substitutions, via a threonine or cysteine intermediate, and are driven by selection. These findings imply a strong pressure of purifying selection in protein evolution, which in the case of serine codon set switches occurs via an initial deleterious substitution quickly followed by a second, compensatory substitution. The result is frequent reversal of amino acid replacements and, at short evolutionary distances, pervasive homoplasy.

Concurrent nucleotide substitution mutations in the human genome are characterized by a significantly decreased transition/transversion ratio

There is accumulating evidence that the number of multiple-nucleotide substitutions (MNS) occurring in closely spaced sites in eukaryotic genomes is significantly higher than would be predicted from the random accumulation of independently generated single-nucleotide substitutions (SNS). Although this excess can in principle be accounted for by the concept of transient hypermutability, a general mutational signature of concurrent MNS mutations has not so far been evident. Employing a dataset (N = 449) of "concurrent" double MNS mutations causing human inherited disease, we have identified just such a mutational signature: concurrently generated double MNS mutations exhibit a >twofold lower transition/transversion ratio (termed RTs/Tv ) than independently generated de novo SNS mutations (<0.80 vs. 2.10; P = 2.69 × 10(-14) ). We replicated this novel finding through a similar analysis employing two double MNS variant datasets with differing abundances of concurrent events (150,521 variants with both substitutions on the same haplotypic lineage vs. 94,875 variants whose component substitutions were on different haplotypic lineages) plus 5,430,874 SNS variants, all being derived from the whole-genome sequencing of seven Chinese individuals. Evaluation of the newly observed mutational signature in diverse contexts provides solid support for the postulated role of translesion synthesis DNA polymerases in transient hypermutability.

Disruption of Transcriptional Coactivator Sub1 Leads to Genome-Wide Re-distribution of Clustered Mutations Induced by APOBEC in Active Yeast Genes

Mutations in genomes of species are frequently distributed non-randomly, resulting in mutation clusters, including recently discovered kataegis in tumors. DNA editing deaminases play the prominent role in the etiology of these mutations. To gain insight into the enigmatic mechanisms of localized hypermutagenesis that lead to cluster formation, we analyzed the mutational single nucleotide variations (SNV) data obtained by whole-genome sequencing of drug-resistant mutants induced in yeast diploids by AID/APOBEC deaminase and base analog 6-HAP. Deaminase from sea lamprey, PmCDA1, induced robust clusters, while 6-HAP induced a few weak ones. We found that PmCDA1, AID, and APOBEC1 deaminases preferentially mutate the beginning of the actively transcribed genes. Inactivation of transcription initiation factor Sub1 strongly reduced deaminase-induced can1 mutation frequency, but, surprisingly, did not decrease the total SNV load in genomes. However, the SNVs in the genomes of the sub1 clones were re-distributed, and the effect of mutation clustering in the regions of transcription initiation was even more pronounced. At the same time, the mutation density in the protein-coding regions was reduced, resulting in the decrease of phenotypically detected mutants. We propose that the induction of clustered mutations by deaminases involves: a) the exposure of ssDNA strands during transcription and loss of protection of ssDNA due to the depletion of ssDNA-binding proteins, such as Sub1, and b) attainment of conditions favorable for APOBEC action in subpopulation of cells, leading to enzymatic deamination within the currently expressed genes. This model is applicable to both the initial and the later stages of oncogenic transformation and explains variations in the distribution of mutations and kataegis events in different tumor cells.

Genomes of tumors are heavily enriched with mutations. Some of these mutations are distributed non-randomly, forming mutational clusters. Editing cytosine deaminases from APOBEC superfamily are responsible for the formation of many of these clusters. We have expressed APOBEC enzyme in diploid yeast cells and found that most of the mutations occur in the beginning of the active genes, where transcription starts. Clusters of mutations overlapped with promoters/transcription start sites. This is likely due to the weaker protection of ssDNA, an ultimate APOBEC deaminase enzyme target, in the beginning of the genes. This hypothesis was reinforced by the finding that inactivation of Sub1 transcription initiation factor, which is found predominantly in the regions of transcription initiation, leads to further increase in mutagenesis in the beginning of the genes. Interestingly, the total number of mutations in the genomes of Sub1-deficient clones did not change, despite the 100-fold decrease in frequency of mutants in a reporter gene. Thus, the drastic change in genome-wide distribution of mutations can be caused by inactivation of a single gene. We propose that the loss of ssDNA protection factors causes formation of mutation clusters in human cancer.

DNA polymerase ζ-dependent lesion bypass in Saccharomyces cerevisiae is accompanied by error-prone copying of long stretches of adjacent DNA

Translesion synthesis (TLS) helps cells to accomplish chromosomal replication in the presence of unrepaired DNA lesions. In eukaryotes, the bypass of most lesions involves a nucleotide insertion opposite the lesion by either a replicative or a specialized DNA polymerase, followed by extension of the resulting distorted primer terminus by DNA polymerase ζ (Polζ). The subsequent events leading to disengagement of the error-prone Polζ from the primer terminus and its replacement with an accurate replicative DNA polymerase remain largely unknown. As a first step toward understanding these events, we aimed to determine the length of DNA stretches synthesized in an error-prone manner during the Polζ-dependent lesion bypass. We developed new in vivo assays to identify the products of mutagenic TLS through a plasmid-borne tetrahydrofuran lesion and a UV-induced chromosomal lesion. We then surveyed the region downstream of the lesion site (in respect to the direction of TLS) for the presence of mutations indicative of an error-prone polymerase activity. The bypass of both lesions was associated with an approximately 300,000-fold increase in the mutation rate in the adjacent DNA segment, in comparison to the mutation rate during normal replication. The hypermutated tract extended 200 bp from the lesion in the plasmid-based assay and as far as 1 kb from the lesion in the chromosome-based assay. The mutation rate in this region was similar to the rate of errors produced by purified Polζ during copying of undamaged DNA in vitro. Further, no mutations downstream of the lesion were observed in rare TLS products recovered from Polζ-deficient cells. This led us to conclude that error-prone Polζ synthesis continues for several hundred nucleotides after the lesion bypass is completed. These results provide insight into the late steps of TLS and show that error-prone TLS tracts span a substantially larger region than previously appreciated.

Genomic instability is associated with multiple genetic diseases. Endogenous and exogenous DNA-damaging factors constitute a major source of genomic instability. Mutations occur when DNA lesions are bypassed by specialized translesion synthesis (TLS) DNA polymerases that are less accurate than the normal replicative polymerases. The discovery of the remarkable infidelity of the TLS enzymes at the turn of the century immediately suggested that their contribution to replication must be tightly restricted to sites of DNA damage to avoid excessive mutagenesis. The actual extent of error-prone synthesis that accompanies TLS in vivo has never been estimated. We describe a novel genetic approach to measure the length of DNA synthesized by TLS polymerases upon their recruitment to sites of DNA damage. We show that stretches of error-prone synthesis associated with the bypass of a single damaged nucleotide span at least 200 and sometimes up to 1,000 nucleotide-long segments, resulting in more than a 300,000-fold increase in mutagenesis in the surrounding region. We speculate that processive synthesis of long DNA stretches by error-prone polymerases could contribute to clustered mutagenesis, a phenomenon that allows for rapid genome changes without significant loss of fitness and plays an important role in tumorigenesis, the immune response and adaptation.

Error-prone polymerase activity causes multinucleotide mutations in humans

About 2% of human genetic polymorphisms have been hypothesized to arise via multinucleotide mutations (MNMs), complex events that generate SNPs at multiple sites in a single generation. MNMs have the potential to accelerate the pace at which single genes evolve and to confound studies of demography and selection that assume all SNPs arise independently. In this paper, we examine clustered mutations that are segregating in a set of 1092 human genomes, demonstrating that the signature of MNM becomes enriched as large numbers of individuals are sampled. We estimate the percentage of linked SNP pairs that were generated by simultaneous mutation as a function of the distance between affected sites and show that MNMs exhibit a high percentage of transversions relative to transitions, findings that are reproducible in data from multiple sequencing platforms and cannot be attributed to sequencing error. Among tandem mutations that occur simultaneously at adjacent sites, we find an especially skewed distribution of ancestral and derived alleles, with GC → AA, GA → TT, and their reverse complements making up 27% of the total. These mutations have been previously shown to dominate the spectrum of the error-prone polymerase Pol ζ, suggesting that low-fidelity DNA replication by Pol ζ is at least partly responsible for the MNMs that are segregating in the human population. We develop statistical estimates of MNM prevalence that can be used to correct phylogenetic and population genetic inferences for the presence of complex mutations.

Genome-wide mutation avalanches induced in diploid yeast cells by a base analog or an APOBEC deaminase

Genetic information should be accurately transmitted from cell to cell; conversely, the adaptation in evolution and disease is fueled by mutations. In the case of cancer development, multiple genetic changes happen in somatic diploid cells. Most classic studies of the molecular mechanisms of mutagenesis have been performed in haploids. We demonstrate that the parameters of the mutation process are different in diploid cell populations. The genomes of drug-resistant mutants induced in yeast diploids by base analog 6-hydroxylaminopurine (HAP) or AID/APOBEC cytosine deaminase PmCDA1 from lamprey carried a stunning load of thousands of unselected mutations. Haploid mutants contained almost an order of magnitude fewer mutations. To explain this, we propose that the distribution of induced mutation rates in the cell population is uneven. The mutants in diploids with coincidental mutations in the two copies of the reporter gene arise from a fraction of cells that are transiently hypersensitive to the mutagenic action of a given mutagen. The progeny of such cells were never recovered in haploids due to the lethality caused by the inactivation of single-copy essential genes in cells with too many induced mutations. In diploid cells, the progeny of hypersensitive cells survived, but their genomes were saturated by heterozygous mutations. The reason for the hypermutability of cells could be transient faults of the mutation prevention pathways, like sanitization of nucleotide pools for HAP or an elevated expression of the PmCDA1 gene or the temporary inability of the destruction of the deaminase. The hypothesis on spikes of mutability may explain the sudden acquisition of multiple mutational changes during evolution and carcinogenesis.

Correlation between mutation rate and genome size in riboviruses: mutation rate of bacteriophage Qβ

Genome sizes and mutation rates covary across all domains of life. In unicellular organisms and DNA viruses, they show an inverse relationship known as Drake’s rule. However, it is still unclear whether a similar relationship exists between genome sizes and mutation rates in RNA genomes. Coronaviruses, the RNA viruses with the largest genomes (∼30 kb), encode a proofreading 3′ exonuclease that allows them to increase replication fidelity. However, it is unknown whether, conversely, the RNA viruses with the smallest genomes tend to show particularly high mutation rates. To test this, we measured the mutation rate of bacteriophage Qβ, a 4.2-kb levivirus. Amber reversion-based Luria–Delbrück fluctuation tests combined with mutant sequencing gave an estimate of 1.4 × 10⁻⁴ substitutions per nucleotide per round of copying, the highest mutation rate reported for any virus using this method. This estimate was confirmed using a direct plaque sequencing approach and after reanalysis of previously published estimates for this phage. Comparison with other riboviruses (all RNA viruses except retroviruses) provided statistical support for a negative correlation between mutation rates and genome sizes. We suggest that the mutation rates of RNA viruses might be optimized for maximal adaptability and that the value of this optimum may in turn depend inversely on genome size.

Patterns and mutational signatures of tandem base substitutions causing human inherited disease

Tandem base substitutions (TBSs) are multiple mutations that comprise two or more contiguous nucleotide substitutions without any net gain or loss of bases. They have recently become recognized as a distinct category of human genomic variant. However, their role in causing human inherited disease so far has not been studied methodically. Here, using data from the Human Gene Mutation Database (http://www.hgmd.org), we identified 477 events to be TBSs (doublets, 448; triplets, 16; and quadruplets to octuplets, 13). A comprehensive sequence pattern and context analysis implied the likely fundamental importance of translesion synthesis (TLS) DNA polymerases in generating these diverse TBSs but revealed that TLS polymerases may operate differently in generating TBSs of ≤ 3 bases (bypass of endogenous DNA lesions) than those of ≥ 4 bases (serial replication slippage). Moreover, GC was found to be the most frequently affected dinucleotide with GC/GC>AA/TT being the most frequent double TBS. Comparison with cancer genome mutational spectra allowed us to conclude that human germline TBSs arise predominantly through the action of endogenous mechanisms of mutagenesis rather than through exposure to exogenous mutagens. Finally, the rates of double and triple TBSs were estimated to be 0.2-1.2 × 10(-10) and 0.8-4.8 × 10(-12) per base per generation, respectively.

Variation in heterozygosity predicts variation in human substitution rates between populations, individuals and genomic regions

The "heterozygote instability" (HI) hypothesis suggests that gene conversion events focused on heterozygous sites during meiosis locally increase the mutation rate, but this hypothesis remains largely untested. As humans left Africa they lost variability, which, if HI operates, should have reduced the mutation rate in non-Africans. Relative substitution rates were quantified in diverse humans using aligned whole genome sequences from the 1,000 genomes project. Substitution rate is consistently greater in Africans than in non-Africans, but only in diploid regions of the genome, consistent with a role for heterozygosity. Analysing the same data partitioned into a series of non-overlapping 2 Mb windows reveals a strong, non-linear correlation between the amount of heterozygosity lost "out of Africa" and the difference in substitution rate between Africans and non-Africans. Putative recent mutations, derived variants that occur only once among the 80 human chromosomes sampled, occur preferentially at the centre of 2 Kb windows that have elevated heterozygosity compared both with the same region in a closely related population and with an immediately adjacent region in the same population. More than half of all substitutions appear attributable to variation in heterozygosity. This observation provides strong support for HI with implications for many branches of evolutionary biology.

Evaluating evolutionary models of stress-induced mutagenesis in bacteria

Increased mutation rates under stress allow bacterial populations to adapt rapidly to stressors, including antibiotics. Here we evaluate existing models for the evolution of stress-induced mutagenesis and present a new model arguing that it evolves as a result of a complex interplay between direct selection for increased stress tolerance, second-order selection for increased evolvability and genetic drift. Further progress in our understanding of the evolutionary biology of stress and mutagenesis will require a more detailed understanding both of the patterns of stress encountered by bacteria in nature and of the mutations that are produced under stress.

Oncogenic potential is related to activating effect of cancer single and double somatic mutations in receptor tyrosine kinases

Aberrant activation of receptor tyrosine kinases (RTKs) is a common feature of many cancer cells. It was previously suggested that the mechanisms of kinase activation in cancer might be linked to transitions between active and inactive states. Here, we estimate the effects of single and double cancer mutations on the stability of active and inactive states of the kinase domains from different RTKs. We show that singleton cancer mutations destabilize active and inactive states; however, inactive states are destabilized more than the active ones, leading to kinase activation. We show that there exists a relationship between the estimate of oncogenic potential of cancer mutation and kinase activation. Namely, more frequent mutations have a higher activating effect, which might allow us to predict the activating effect of the mutations from the mutation spectra. Independent evolutionary analysis of mutation spectra complements this observation and finds the same frequency threshold defining mutation hotspots. We analyze double mutations and report a positive epistasis and additional advantage of doublets with respect to cancer cell fitness. The activation mechanisms of double mutations differ from those of single mutations and double mutation spectrum is found to be dissimilar to the mutation spectrum of singletons.

Transient hypermutability, chromothripsis and replication-based mechanisms in the generation of concurrent clustered mutations

Clustered mutations may be broadly defined as the presence of two or more mutations within a spatially localized genomic region on a single chromosome. Known instances vary in terms of both the number and type of the component mutations, ranging from two closely spaced point mutations to tens or even hundreds of genomic rearrangements. Although clustered mutations can represent the observable net result of independent lesions sequentially acquired over multiple cell cycles, they can also be generated in a simultaneous or quasi-simultaneous manner within a single cell cycle. This review focuses on those mechanisms known to underlie the latter type. Both gene conversion and transient hypermutability are capable of generating closely spaced multiple mutations. However, a recently described phenomenon in human cancer cells, known as 'chromothripsis', has provided convincing evidence that tens to hundreds of genomic rearrangements can sometimes be generated simultaneously via a single catastrophic event. The distinctive genomic features observed in the derivative chromosomes, together with the highly characteristic junction sequences, point to non-homologous end joining (NHEJ) as being the likely underlying mutational mechanism. By contrast, replication-based mechanisms such as microhomology-mediated break-induced replication (MMBIR) which involves serial replication slippage or serial template switching probably account for those complex genomic rearrangements that comprise multiple duplications and/or triplications.

The three faces of riboviral spontaneous mutation: spectrum, mode of genome replication, and mutation rate

Riboviruses (RNA viruses without DNA replication intermediates) are the most abundant pathogens infecting animals and plants. Only a few riboviral infections can be controlled with antiviral drugs, mainly because of the rapid appearance of resistance mutations. Little reliable information is available concerning i) kinds and relative frequencies of mutations (the mutational spectrum), ii) mode of genome replication and mutation accumulation, and iii) rates of spontaneous mutation. To illuminate these issues, we developed a model in vivo system based on phage Qß infecting its natural host, Escherichia coli. The Qß RT gene encoding the Read-Through protein was used as a mutation reporter. To reduce uncertainties in mutation frequencies due to selection, the experimental Qß populations were established after a single cycle of infection and selection against RT⁻ mutants during phage growth was ameliorated by plasmid-based RT complementation in trans. The dynamics of Qß genome replication were confirmed to reflect the linear process of iterative copying (the stamping-machine mode). A total of 32 RT mutants were detected among 7,517 Qß isolates. Sequencing analysis of 45 RT mutations revealed a spectrum dominated by 39 transitions, plus 4 transversions and 2 indels. A clear template•primer mismatch bias was observed: A•C>C•A>U•G>G•U> transversion mismatches. The average mutation rate per base replication was ≈9.1×10⁻⁶ for base substitutions and ≈2.3×10⁻⁷ for indels. The estimated mutation rate per genome replication, μ_g, was ≈0.04 (or, per phage generation, ≈0.08), although secondary RT mutations arose during the growth of some RT mutants at a rate about 7-fold higher, signaling the possible impact of transitory bouts of hypermutation. These results are contrasted with those previously reported for other riboviruses to depict the current state of the art in riboviral mutagenesis.

Viral disease is a subject of major concern in public health. Diseases produced by riboviruses (RNA viruses sensu stricto) represent a special urgency, because these viruses display an exceptional capability to generate resistance mutations against antiviral drugs. Unfortunately, little is known about the rate and nature of spontaneous mutation in riboviruses. Thus, characterization of their mutation process may be helpful in the development of improved ways to counteract riboviral diseases. In this study, we investigated the mutation process in vivo of a model ribovirus, the bacteriophage Qß, focusing on three key aspects: i) the kinds and relative frequencies of mutations, ii) the mode of genome replication, and iii) the rate of spontaneous mutation. Our results, combined with other information about riboviral mutagenesis, depict a ribovirus mutation spectrum largely dominated by transitions, a predominantly linear mode of genome replication, and a mutation rate per genome replication on the order of 0.04 for bacteriophages and plant viruses but perhaps an order of magnitude higher for mammalian riboviruses.

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.

Two closely spaced nonsense mutations in the DMD gene in a Malaysian family

In Duchenne muscular dystrophy (DMD), identification of one nonsense mutation in the DMD gene has been considered an endpoint of genetic diagnosis. Here, we identified two closely spaced nonsense mutations in the DMD gene. In a Malaysian DMD patient two nonsense mutations (p.234S>X and p.249Q>X, respectively) were identified within exon 8. The proband's mother carried both mutations on one allele. Multiple mutations may explain the occasional discrepancies between genotype and phenotype in dystrophinopathy.

Gene conversion in human genetic disease

Gene conversion is a specific type of homologous recombination that involves the unidirectional transfer of genetic material from a ‘donor’ sequence to a highly homologous ‘acceptor’. We have recently reviewed the molecular mechanisms underlying gene conversion, explored the key part that this process has played in fashioning extant human genes, and performed a meta-analysis of gene-conversion events known to have caused human genetic disease. Here we shall briefly summarize some of the latest developments in the study of pathogenic gene conversion events, including (i) the emerging idea of minimal efficient sequence homology (MESH) for homologous recombination, (ii) the local DNA sequence features that appear to predispose to gene conversion, (iii) a mechanistic comparison of gene conversion and transient hypermutability, and (iv) recently reported examples of pathogenic gene conversion events.

A single-strand specific lesion drives MMS-induced hyper-mutability at a double-strand break in yeast

Localized hyper-mutability (LHM) can be important in evolution, immunity, and genetic diseases. We previously reported that single-strand DNA (ssDNA) can be an important source of damage-induced LHM in yeast. Here, we establish that the generation of LHM by methyl methanesulfonate (MMS) during repair of a chromosomal double-strand break (DSB) can result in over 0.2 mutations/kb, which is approximately 20,000-fold higher than the MMS-induced mutation density without a DSB. The MMS-induced mutations associated with DSB repair were primarily due to substitutions via translesion DNA synthesis at damaged cytosines, even though there are nearly 10 times more MMS-induced lesions at other bases. Based on this mutation bias, the promutagenic lesion dominating LHM is likely 3-methylcytosine, which is single-strand specific. Thus, the dramatic increase in mutagenesis at a DSB is concluded to result primarily from the generation of non-repairable lesions in ssDNA associated with DSB repair along with efficient induction of highly mutagenic ssDNA-specific lesions. These findings with MMS-induced LHM have broad biological implications for unrepaired damage generated in ssDNA and possibly ssRNA.

Mutational clusters generated by non-processive polymerases: A case study using DNA polymerase betain vitro

Available DNA mutational spectra reveal that the number of mutants with multiple mutations ("multiples") is usually greater than expected from a random distribution of mutations among mutants. These overloads imply the occurrence of non-random clusters of mutations, probably generated during episodes of low-fidelity DNA synthesis. Excess multiples have been reported not only for viruses, bacteria, and eukaryotic cells but also for the DNA polymerases of phages T4 and RB69 in vitro. In the simplest case of a purified polymerase, non-random clusters may be generated by a subfraction of phenotypic variants able to introduce more errors per cycle of DNA synthesis than the normal enzyme. According to this hypothesis, excess multiples are not expected with non-processive polymerases even if they harbor rare mutator variants. DNA polymerase beta (Pol beta) is a mammalian DNA-repair polymerase with very low processivity. Although several Pol beta mutational spectra have been described, there is conflicting evidence on whether or not excess multiples occur, with spectra based on the HSV-tk system tending to show excess multiples. Excess multiples generated by Pol beta or any of its mutants might imply that the excesses of multiples observed in numerous other systems, especially those with processive polymerases, could be artifactual. Here, the distributions of mutations generated by native and recombinant rat Pol beta and by the Pol beta(Y265C) mutator were analyzed in the M13mp2 lacZalpha system. Our results present no evidence for a significant excess of multiples over the expected numbers with any of the Pol beta enzymes tested in this system. The reported excess of Pol beta-generated multiples in the HSV-tk system may reflect a reduced efficiency of detection of base substitutions that cause weak phenotypes, which in turn may artifactually increase the frequency of multiples.

The rate of establishment of complex adaptations

A central problem in evolutionary theory concerns the mechanisms by which adaptations requiring multiple mutations emerge in natural populations. We develop a series of expressions that clarify the scaling of the time to establishment of complex adaptations with population size, mutation rate, magnitude of the selective disadvantage of intermediate-state alleles, and the complexity of the adaptation. In general, even in the face of deleterious intermediate steps, the time to establishment is minimized in populations with very large size. Under a broad range of conditions, the time to establishment also scales by no more than the square of the mutation rate, regardless of the number of sites contributing to the adaptive change, demonstrating that the emergence of complex adaptations is only weakly constrained by the independent acquisition of mutations at the underlying sites. Mutator alleles with deleterious side effects have only moderate effects on the rate of adaptation in large populations but can cause a quantum decrease in the time to establishment of some adaptive alleles in small populations, although probably not at a high enough rate to offset the increased deleterious mutation load. Transient hypermutability, whereby a subset of gamete-producing cells mutate at an elevated rate in a nonheritable manner, may also elevate the rate of adaptation, although the effect is modest and appears to result from a simple increase in the rate of transitions between intermediate states rather than from the saltational production of doublet mutations. Taken together, these results illustrate the plausibility of the relatively rapid emergence of specific complex adaptations by conventional population genetic mechanisms and provide insight into the relative incidences of various paths of allelic adaptation in organisms with different population genetic features.

Heterogeneous distribution of SNPs in the human genome: microsatellites as predictors of nucleotide diversity and divergence

Understanding the forces that govern the distribution of single nucleotide polymorphisms is vital for many of their applications. Here we conducted a systematic search to quantify how both SNP density and human-chimpanzee divergence vary around different repetitive sequences. We uncovered a highly complicated picture in which these quantities often differ significantly from the genome-wide average in regions extending more than 20 kb, the direction of the deviation varying with repeat number and motif. AT microsatellites in particular are potent predictors of SNP density, long (AT)(n) repeat tracts tending to be found in regions of significantly reduced SNP density and low GC content. Although the causal relationships remain difficult to determine, our results indicate a strong relationship between microsatellites and the DNA that flanks them. Our results help to explain the mixed picture that emerges from other studies and have important implications for the way in which genetic diversity is distributed in our genomes.

Closely spaced multiple mutations as potential signatures of transient hypermutability in human genes

Data from diverse organisms suggests that transient hypermutability is a general mutational mechanism with the potential to generate multiple synchronous mutations, a phenomenon probably best exemplified by closely spaced multiple mutations (CSMMs). Here we have attempted to extend the concept of transient hypermutability from somatic cells to the germline, using human inherited disease-causing multiple mutations as a model system. Employing stringent criteria for data inclusion, we have retrospectively identified numerous potential examples of pathogenic CSMMs that exhibit marked similarities to the CSMMs reported in other systems. These examples include (1) eight multiple mutations, each comprising three or more components within a sequence tract of <100 bp; (2) three possible instances of "mutation showers"; and (3) numerous highly informative "homocoordinate" mutations. Using the proportion of CpG substitution as a crude indicator of the relative likelihood of transient hypermutability, we present evidence to suggest that CSMMs comprising at least one pair of mutations separated by < or =100 bp may constitute signatures of transient hypermutability in human genes. Although this analysis extends the generality of the concept of transient hypermutability and provides new insights into what may be considered a novel mechanism of mutagenesis underlying human inherited disease, it has raised serious concerns regarding current practices in mutation screening.

A cancer-associated DNA polymerase delta variant modeled in yeast causes a catastrophic increase in genomic instability

Accurate DNA synthesis by the replicative DNA polymerases alpha, delta, and epsilon is critical for genome stability in eukaryotes. In humans, over 20 SNPs were reported that result in amino-acid changes in Poldelta or Polepsilon. In addition, Poldelta variants were found in colon-cancer cell lines and in sporadic colorectal carcinomas. Using the yeast-model system, we examined the functional consequences of two cancer-associated Poldelta mutations and four polymorphisms affecting well-conserved regions of Poldelta or Polepsilon. We show that the R696W substitution in Poldelta (analog of the R689W change in the human cancer-cell line DLD-1) is lethal in haploid and homozygous diploid yeast. The cell death results from a catastrophic increase in spontaneous mutagenesis attributed to low-fidelity DNA synthesis by Poldelta-R696W. Heterozygotes survive, and the mutation rate depends on the relative expression level of wild-type versus mutant alleles. Based on these observations, we propose that the mutation rate in heterozygous human cells could be regulated by transient changes in gene expression leading to a temporary excess of Poldelta-R689W. The similarities between the mutational spectra of the yeast strains producing Poldelta-R696W and DLD-1 cells suggest that the altered Poldelta could be responsible for a significant proportion of spontaneous mutations in this cancer cell line. These results suggest that the highly error-prone Poldelta-R689W could contribute to cancer initiation and/or progression in humans.

Genomic mutation rates: what high-throughput methods can tell us

High-throughput DNA analyses are increasingly being used to detect rare mutations in moderately sized genomes. These methods have yielded genome mutation rates that are markedly higher than those obtained using pre-genomic strategies. Recent work in a variety of organisms has shown that mutation rate is strongly affected by sequence context and genome position. These observations suggest that high-throughput DNA analyses will ultimately allow researchers to identify trans-acting factors and cis sequences that underlie mutation rate variation. Such work should provide insights on how mutation rate variability can impact genome organization and disease progression.

Low-fidelity DNA synthesis by the L979F mutator derivative of Saccharomyces cerevisiae DNA polymerase zeta

To probe Pol ζ functions in vivo via its error signature, here we report the properties of Saccharomyces cerevisiae Pol ζ in which phenyalanine was substituted for the conserved Leu-979 in the catalytic (Rev3) subunit. We show that purified L979F Pol ζ is 30% as active as wild-type Pol ζ when replicating undamaged DNA. L979F Pol ζ shares with wild-type Pol ζ the ability to perform moderately processive DNA synthesis. When copying undamaged DNA, L979F Pol ζ is error-prone compared to wild-type Pol ζ, providing a biochemical rationale for the observed mutator phenotype of rev3-L979F yeast strains. Errors generated by L979F Pol ζ in vitro include single-base insertions, deletions and substitutions, with the highest error rates involving stable misincorporation of dAMP and dGMP. L979F Pol ζ also generates multiple errors in close proximity to each other. The frequency of these events far exceeds that expected for independent single changes, indicating that the first error increases the probability of additional errors within 10 nucleotides. Thus L979F Pol ζ, and perhaps wild-type Pol ζ, which also generates clustered mutations at a lower but significant rate, performs short patches of processive, error-prone DNA synthesis. This may explain the origin of some multiple clustered mutations observed in vivo.

Discovery of proteomic code with mRNA assisted protein folding

The 3x redundancy of the Genetic Code is usually explained as a necessity to increase the mutation-resistance of the genetic information. However recent bioinformatical observations indicate that the redundant Genetic Code contains more biological information than previously known and which is additional to the 64/20 definition of amino acids. It might define the physico-chemical and structural properties of amino acids, the codon boundaries, the amino acid co-locations (interactions) in the coded proteins and the free folding energy of mRNAs. This additional information, which seems to be necessary to determine the 3D structure of coding nucleic acids as well as the coded proteins, is known as the Proteomic Code and mRNA Assisted Protein Folding.

Hypermutability of damaged single-strand DNA formed at double-strand breaks and uncapped telomeres in yeast Saccharomyces cerevisiae

The major DNA repair pathways operate on damage in double-strand DNA because they use the intact strand as a template after damage removal. Therefore, lesions in transient single-strand stretches of chromosomal DNA are expected to be especially threatening to genome stability. To test this hypothesis, we designed systems in budding yeast that could generate many kilobases of persistent single-strand DNA next to double-strand breaks or uncapped telomeres. The systems allowed controlled restoration to the double-strand state after applying DNA damage. We found that lesions induced by UV-light and methyl methanesulfonate can be tolerated in long single-strand regions and are hypermutagenic. The hypermutability required PCNA monoubiquitination and was largely attributable to translesion synthesis by the error-prone DNA polymerase ζ. In support of multiple lesions in single-strand DNA being a source of hypermutability, analysis of the UV-induced mutants revealed strong strand-specific bias and unexpectedly high frequency of alleles with widely separated multiple mutations scattered over several kilobases. Hypermutability and multiple mutations associated with lesions in transient stretches of long single-strand DNA may be a source of carcinogenesis and provide selective advantage in adaptive evolution.

A variety of error avoidance mechanisms assure low mutation rates across the genome. Genetic defects in DNA replication or repair can lead to genome-wide increase in mutation frequency that may result in cancer predisposition and genetic disease. Transient localized hypermutability drastically differs in its biological consequences from genome-wide mutators. Since genome-wide hypermutability can cause reduced fitness due to accumulation of dysfunctional alleles, mutators are under negative selection pressure. By contrast, there would be less selection against temporary hypermutability within limited genomic regions, suggesting a special role in adaptive evolution and carcinogenesis. Mechanisms of transient hypermutability are poorly understood. Long stretches of single-strand DNA have been implicated but not demonstrated as a source of localized transient hypermutability. Using sophisticated yeast genetic systems that we developed, we found that transient stretches of chromosomal single-strand DNA at double-strand breaks and that telomeres can tolerate multiple lesions and are highly prone to damage-induced mutations, including a very unusual class of widely spaced multiple mutations. The hypermutability relied on error prone translesion DNA synthesis. Our work demonstrates a simple in vivo mechanism for localized transient hypermutability extending over several kilobases that can result in widely spaced multiple mutations without severe mutation load in the rest of the genome.

Epidemiology of doublet/multiplet mutations in lung cancers: evidence that a subset arises by chronocoordinate events

Background

Evidence strongly suggests that spontaneous doublet mutations in normal mouse tissues generally arise from chronocoordinate events. These chronocoordinate mutations sometimes reflect “mutation showers”, which are multiple chronocoordinate mutations spanning many kilobases. However, little is known about mutagenesis of doublet and multiplet mutations (domuplets) in human cancer. Lung cancer accounts for about 25% of all cancer deaths. Herein, we analyze the epidemiology of domuplets in the EGFR and TP53 genes in lung cancer. The EGFR gene is an oncogene in which doublets are generally driver plus driver mutations, while the TP53 gene is a tumor suppressor gene with a more typical situation in which doublets derive from a driver and passenger mutation.

Methodology/Principal Findings

EGFR mutations identified by sequencing were collected from 66 published papers and our updated EGFR mutation database (www.egfr.org). TP53 mutations were collected from IARC version 12 (www-p53.iarc.fr). For EGFR and TP53 doublets, no clearly significant differences in race, ethnicity, gender and smoking status were observed. Doublets in the EGFR and TP53 genes in human lung cancer are elevated about eight- and three-fold, respectively, relative to spontaneous doublets in mouse (6% and 2.3% versus 0.7%).

Conclusions/Significance

Although no one characteristic is definitive, the aggregate properties of doublet and multiplet mutations in lung cancer are consistent with a subset derived from chronocoordinate events in the EGFR gene: i) the eight frameshift doublets (present in 0.5% of all patients with EGFR mutations) are clustered and produce a net in-frame change; ii) about 32% of doublets are very closely spaced (≤30 nt); and iii) multiplets contain two or more closely spaced mutations. TP53 mutations in lung cancer are very closely spaced (≤30 nt) in 33% of doublets, and multiplets generally contain two or more very closely spaced mutations. Work in model systems is necessary to confirm the significance of chronocoordinate events in lung and other cancers.

Mutability and importance of a hypermutable cell subpopulation that produces stress-induced mutants in Escherichia coli

In bacterial, yeast, and human cells, stress-induced mutation mechanisms are induced in growth-limiting environments and produce non-adaptive and adaptive mutations. These mechanisms may accelerate evolution specifically when cells are maladapted to their environments, i.e., when they are are stressed. One mechanism of stress-induced mutagenesis in Escherichia coli occurs by error-prone DNA double-strand break (DSB) repair. This mechanism was linked previously to a differentiated subpopulation of cells with a transiently elevated mutation rate, a hypermutable cell subpopulation (HMS). The HMS could be important, producing essentially all stress-induced mutants. Alternatively, the HMS was proposed to produce only a minority of stress-induced mutants, i.e., it was proposed to be peripheral. We characterize three aspects of the HMS. First, using improved mutation-detection methods, we estimate the number of mutations per genome of HMS-derived cells and find that it is compatible with fitness after the HMS state. This implies that these mutants are not necessarily an evolutionary dead end, and could contribute to adaptive evolution. Second, we show that stress-induced Lac(+) mutants, with and without evidence of descent from the HMS, have similar Lac(+) mutation sequences. This provides evidence that HMS-descended and most stress-induced mutants form via a common mechanism. Third, mutation-stimulating DSBs introduced via I-SceI endonuclease in vivo do not promote Lac(+) mutation independently of the HMS. This and the previous finding support the hypothesis that the HMS underlies most stress-induced mutants, not just a minority of them, i.e., it is important. We consider a model in which HMS differentiation is controlled by stress responses. Differentiation of an HMS potentially limits the risks of mutagenesis in cell clones.

High-precision, whole-genome sequencing of laboratory strains facilitates genetic studies

Whole-genome sequencing is a powerful technique for obtaining the reference sequence information of multiple organisms. Its use can be dramatically expanded to rapidly identify genomic variations, which can be linked with phenotypes to obtain biological insights. We explored these potential applications using the emerging next-generation sequencing platform Solexa Genome Analyzer, and the well-characterized model bacterium Bacillus subtilis. Combining sequencing with experimental verification, we first improved the accuracy of the published sequence of the B. subtilis reference strain 168, then obtained sequences of multiple related laboratory strains and different isolates of each strain. This provides a framework for comparing the divergence between different laboratory strains and between their individual isolates. We also demonstrated the power of Solexa sequencing by using its results to predict a defect in the citrate signal transduction pathway of a common laboratory strain, which we verified experimentally. Finally, we examined the molecular nature of spontaneously generated mutations that suppress the growth defect caused by deletion of the stringent response mediator relA. Using whole-genome sequencing, we rapidly mapped these suppressor mutations to two small homologs of relA. Interestingly, stable suppressor strains had mutations in both genes, with each mutation alone partially relieving the relA growth defect. This supports an intriguing three-locus interaction module that is not easily identifiable through traditional suppressor mapping. We conclude that whole-genome sequencing can drastically accelerate the identification of suppressor mutations and complex genetic interactions, and it can be applied as a standard tool to investigate the genetic traits of model organisms.

In this manuscript, we describe novel applications of the newly developed Solexa sequencing technology. We aim to provide insights into the following questions: (1) Can whole-genome sequencing, while rapidly surveying mega-bases of genome information, also reliably identify variations at the base-pair resolution? (2) Can it be used to identify the differences between isolates of the same laboratory strain and between different laboratory strains? (3) Can it be used as a genetic tool to predict phenotypes and identify suppressors? To this end, we performed whole-genome shotgun sequencing of several related strains of the widely studied model bacterium Bacillus subtilis, we identified genomic variations that potentially underlie strain-specific phenotypes, which occur frequently in biological studies, and we found multiple suppressor mutations within a single strain that are difficult to discern through traditional methods. We conclude that whole-genome sequencing can be directly used to guide genetic studies.

Rates of evolutionary change in viruses: patterns and determinants

Understanding the factors that determine the rate at which genomes generate and fix mutations provides important insights into key evolutionary mechanisms. We review our current knowledge of the rates of mutation and substitution, as well as their determinants, in RNA viruses, DNA viruses and retroviruses. We show that the high rate of nucleotide substitution in RNA viruses is matched by some DNA viruses, suggesting that evolutionary rates in viruses are explained by diverse aspects of viral biology, such as genomic architecture and replication speed, and not simply by polymerase fidelity.

Double complex mutations involving F8 and FUNDC2 caused by distinct break-induced replication

Genomic rearrangements are a well-recognized cause of genetic disease and can be formed by a variety of mechanisms. We report a complex rearrangement causing severe hemophilia A, identified and further characterized using a range of PCR-based methods, and confirmed using array-comparative genomic hybridization (array-CGH). This rearrangement consists of a 15.5-kb deletion/16-bp insertion located 0.6 kb from a 28.1-kb deletion/263-kb insertion at Xq28 and is one of the most complex rearrangements described at a DNA sequence level. We propose that the rearrangement was generated by distinct but linked cellular responses to double strand breakage, namely break-induced replication (BIR) and a novel model of break-induced serial replication slippage (SRS). The copy number of several genes is affected by this rearrangement, with deletion of part of the Factor VIII gene (F8, causing hemophilia A) and the FUNDC2 gene, and duplication of the TMEM185A, HSFX1, MAGEA9, and MAGEA11 genes. As the patient exhibits no clinically detectable phenotype other than hemophilia A, it appears that the biological effects of the other genes involved are not dosage-dependent. This investigation has provided novel insights into processes of DNA repair including BIR and the first description of SRS during repair in a pathological context.