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(-6) for base substitutions and ≈2.3×10(-7) 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.

Confounders of mutation-rate estimators: selection and phenotypic lag in Thermus thermophilus.

In a recent description of the rate and character of spontaneous mutation in the hyperthermophilic bacterium Thermus thermophilus, the mutation rate was observed to be substantially lower than seen in several mesophiles. Subsequently, a report appeared indicating that this bacterium maintains an average of about 4.5 genomes per cell. This number of genomes might result in a segregation lag for the expression of a recessive mutation and might therefore lead to an underestimate of the rate of mutation. Here we describe some kinds of problems that may arise when estimating mutation rates and outline ways to adjust the rates accordingly. The emphasis is mainly on differential rates of growth of mutants versus their parents and on various kinds of phenotypic lag. We then apply these methods to the T. thermophilus data and conclude that there is as yet no reliable impact on a previously described rate.

Avoiding dangerous missense: thermophiles display especially low mutation rates.

Rates of spontaneous mutation have been estimated under optimal growth conditions for a variety of DNA-based microbes, including viruses, bacteria, and eukaryotes. When expressed as genomic mutation rates, most of the values were in the vicinity of 0.003-0.004 with a range of less than two-fold. Because the genome sizes varied by roughly 10(4)-fold, the mutation rates per average base pair varied inversely by a similar factor. Even though the commonality of the observed genomic rates remains unexplained, it implies that mutation rates in unstressed microbes reach values that can be finely tuned by evolution. An insight originating in the 1920s and maturing in the 1960s proposed that the genomic mutation rate would reflect a balance between the deleterious effect of the average mutation and the cost of further reducing the mutation rate. If this view is correct, then increasing the deleterious impact of the average mutation should be countered by reducing the genomic mutation rate. It is a common observation that many neutral or nearly neutral mutations become strongly deleterious at higher temperatures, in which case they are called temperature-sensitive mutations. Recently, the kinds and rates of spontaneous mutations were described for two microbial thermophiles, a bacterium and an archaeon. Using an updated method to extrapolate from mutation-reporter genes to whole genomes reveals that the rate of base substitutions is substantially lower in these two thermophiles than in mesophiles. This result provides the first experimental support for the concept of an evolved balance between the total genomic impact of mutations and the cost of further reducing the basal mutation rate.

The bacteriophage T4 rapid-lysis genes and their mutational proclivities.

Like most phages with double-stranded DNA, phage T4 exits the infected host cell by a lytic process requiring, at a minimum, an endolysin and a holin. Unlike most phages, T4 can sense superinfection (which signals the depletion of uninfected host cells) and responds by delaying lysis and achieving an order-of-magnitude increase in burst size using a mechanism called lysis inhibition (LIN). T4 r mutants, which are unable to conduct LIN, produce distinctly large, sharp-edged plaques. The discovery of r mutants was key to the foundations of molecular biology, in particular to discovering and characterizing genetic recombination in T4, to redefining the nature of the gene, and to exploring the mutation process at the nucleotide level of resolution. A number of r genes have been described in the past 7 decades with various degrees of clarity. Here we describe an extensive and perhaps saturating search for T4 r genes and relate the corresponding mutational spectra to the often imperfectly known physiologies of the proteins encoded by these genes. Focusing on r genes whose mutant phenotypes are largely independent of the host cell, the genes are rI (which seems to sense superinfection and signal the holin to delay lysis), rIII (of poorly defined function), rIV (same as sp and also of poorly defined function), and rV (same as t, the holin gene). We did not identify any mutations that might correspond to a putative rVI gene, and we did not focus on the famous rII genes because they appear to affect lysis only indirectly.

Templated mutagenesis in bacteriophage T4 involving imperfect direct or indirect sequence repeats.

Some mutations arise in association with a potential sequence donor that consists of an imperfect direct or reverse repeat. Many such mutations are complex; that is, they consist of multiple close sequence changes. Current models posit that the primer terminus of a replicating DNA molecule dissociates, reanneals with an ectopic template, extends briefly, and then returns to the cognate template, bringing with it a locally different sequence; alternatively, a hairpin structure may form the mutational intermediate when processed by mismatch repair. This process resembles replication repair, in which primer extension is blocked by a lesion in the template; in this case, the ectopic template is the other daughter strand, and the result is error-free bypass of the lesion. We previously showed that mutations that impair replication repair can enhance templated mutagenesis. We show here that the intensity of templated mutation can be exquisitely sensitive to its local sequence, that the donor and recipient arms of an imperfect inverse repeat can exchange roles, and that double mutants carrying two alleles, each affecting both templated mutagenesis and replication repair, can have unexpected phenotypes. We also record an instance in which the mutation rates at two particular sites change concordantly with a distant sequence change, but in a manner that appears unrelated to templated mutagenesis.

The rate and character of spontaneous mutation in Thermus thermophilus.

Selection of spontaneous, loss-of-function mutations at two chromosomal loci (pyrF and pyrE) enabled the first molecular-level analysis of replication fidelity in the extremely thermophilic bacterium Thermus thermophilus. Two different methods yielded similar mutation rates, and mutational spectra determined by sequencing of independent mutants revealed a variety of replication errors distributed throughout the target genes. The genomic mutation rate estimated from these targets, 0.00097 +/- 0.00052 per replication, was lower than corresponding estimates from mesophilic microorganisms, primarily because of a low rate of base substitution. However, both the rate and spectrum of spontaneous mutations in T. thermophilus resembled those of the thermoacidophilic archaeon Sulfolobus acidocaldarius, despite important molecular differences between these two thermophiles and their genomes.

The roles of Tyr391 and Tyr619 in RB69 DNA polymerase replication fidelity.

In the family-B DNA polymerase of bacteriophage RB69, the conserved aromatic palm-subdomain residues Tyr391 and Tyr619 interact with the last primer-template base-pair. Tyr619 interacts via a water-mediated hydrogen bond with the phosphate of the terminal primer nucleotide. The main-chain amide of Tyr391 interacts with the corresponding template nucleotide. A hydrogen bond has been postulated between Tyr391 and the hydroxyl group of Tyr567, a residue that plays a key role in base discrimination. This hydrogen bond may be crucial for forcing an infrequent Tyr567 rotamer conformation and, when the bond is removed, may influence fidelity. We investigated the roles of these residues in replication fidelity in vivo employing phage T4 rII reversion assays and an rI forward assay. Tyr391 was replaced by Phe, Met and Ala, and Tyr619 by Phe. The Y391A mutant, reported previously to decrease polymerase affinity for incoming nucleotides, was unable to support DNA replication in vivo, so we used an in vitro fidelity assay. Tyr391F/M replacements affect fidelity only slightly, implying that the bond with Tyr567 is not essential for fidelity. The Y391A enzyme has no mutator phenotype in vitro. The Y619F mutant displays a complex profile of impacts on fidelity but has almost the same mutational spectrum as the parental enzyme. The Y619F mutant displays reduced DNA binding, processivity, and exonuclease activity on single-stranded DNA and double-stranded DNA substrates. The Y619F substitution would disrupt the hydrogen bond network at the primer terminus and may affect the alignment of the 3' primer terminus at the polymerase active site, slowing chemistry and overall DNA synthesis.

An unusual pattern of spontaneous mutations recovered in the halophilic archaeon Haloferax volcanii.

Spontaneous mutations in the orotate:phosphoribosyl transferase (pyrE2) gene of the halophilic archaeon Haloferax volcanii were selected by 5-fluoroorotic acid plus uracil at a rate of approximately 2 x 10(-8)/cell division in fluctuation and null-fraction tests but approximately 6 x 10(-8)/cell division in mutation-accumulation tests. The corresponding genomic mutation rates were substantially lower than those observed for other mesophilic microbial DNA genomes on the basis of similar target genes. The mutational spectrum was dominated by indels adding or deleting multiples of 3 bp. Properties of the organism contributing to this unusual mutational pattern may include phenotypic lag caused by a high chromosomal copy number and efficient promotion of strand misalignments by short direct repeats.

A role for replication repair in the genesis of templated mutations.

Replication repair mediates error-free bypass of DNA damage in a series of steps that include regression of the replication fork, primer-terminus switching to use the other daughter strand as an undamaged template, primer extension, primer switching back to its cognate template with the primer terminus now having bypassed the damage, and fork rearrangement to a normal configuration. By both genetic and biochemical criteria, bacteriophage T4 catalyzes replication repair with two alternative sets of proteins, one including the gp32 SSB and the gp41 DNA helicase and the other including the UvsX recombinase. In each pathway, synthesis is conducted by the gp43 DNA polymerase. Here we show that defects in gp32, gp41 or UvsX that impair replication repair also increase mutation rates generally, but especially for templated mutations. Such templated mutations are associated with palindromic or direct repeats that are either perfect or imperfect. Models of templated mutagenesis require that the primer terminus switches to an ectopic template, but one that yields mutations instead of error-free bypass. We suggest that the proteins that conduct replication repair normally direct a blocked primer strand specifically to the other daughter strand with considerable accuracy, but that strand switching becomes promiscuous when these proteins are mutationally impaired, thus promoting templated mutations.

On the mutation rate of herpes simplex virus type 1.

All seven DNA-based microbes for which carefully established mutation rates and mutational spectra were previously available displayed a genomic mutation rate in the neighborhood of 0.003 per chromosome replication. The pathogenic mammalian DNA virus herpes simplex type 1 has an estimated genomic mutation rate compatible with that value.

Processivity clamp gp45 and ssDNA-binding-protein gp32 modulate the fidelity of bacteriophage RB69 DNA polymerase in a sequence-specific manner, sometimes enhancing and sometimes compromising accuracy.

Numerous studies of the impact of accessory proteins upon the fidelity of DNA synthesis have provided a complex and sometimes discordant picture. We previously described such an analysis conducted in vitro using various bacteriophage RB69 gp43 mutator DNA polymerases with or without the accessory proteins gp32 (which binds single-stranded DNA) plus gp45/44/62 (processivity clamp and its loaders). Mutations were scored at many sites in the lacZalpha mutation reporter sequence. Unexpectedly, the accessory proteins sometimes decreased and sometimes increased fidelity at a handful of specific sites. Here, we enlarge our analysis with one particular mutator polymerase compromised in both insertion accuracy and proofreading and also extend the analysis to reactions supplemented only with gp32 or only with gp45/44/62. An overall 1.56-fold increase in mutation frequencies was produced by adding single or multiple accessory proteins and was driven mainly by increased T(template)*G(primer) mispairs. Evidence was found for many additional sites where the accessory proteins influence fidelity, indicating the generality of the effect. Thus, accessory proteins contribute to the site-specific variability in mutation rates characteristically seen in mutational spectra.

Characterization of DNA synthesis catalyzed by bacteriophage T4 replication complexes reconstituted on synthetic circular substrates.

Replication complexes were reconstituted using the eight purified bacteriophage T4 replication proteins and synthetic circular 70-, 120- or 240-nt DNA substrates annealed to a leading-strand primer. To differentiate leading strands from lagging strands, the circular parts of the substrates lacked dCMP; thus, no dCTP was required for leading-strand synthesis and no dGTP for lagging-strand synthesis. The size of the substrates was crucial, the longer substrates supporting much more DNA synthesis. Leading and lagging strands were synthesized in a coupled manner. Specifically targeting leading-strand synthesis by decreasing the concentration of dGTP decreased the rate of extension of leading strands. However, blocking lagging-strand synthesis by lowering the dCTP concentration, by omitting dCTP altogether, by adding ddCTP, or with a single abasic site had no immediate effect on the rate of extension of leading strands.

Genome-wide patterns of nucleotide substitution reveal stringent functional constraints on the protein sequences of thermophiles.

To test the hypothesis that the proteins of thermophilic prokaryotes are subject to unusually stringent functional constraints, we estimated the numbers of synonymous and nonsynonymous nucleotide substitutions per site between 17,957 pairs of orthologous genes from 22 pairs of closely related species of Archaea and Bacteria. The average ratio of nonsynonymous to synonymous substitutions was significantly lower in thermophiles than in nonthermophiles, and this effect was observed in both Archaea and Bacteria. There was no evidence that this difference could be explained by factors such as nucleotide content bias. Rather, the results support the hypothesis that proteins of thermophiles are subject to unusually strong purifying selection, leading to a reduced overall level of amino acid evolution per mutational event. The results show that genome-wide patterns of sequence evolution can be influenced by natural selection exerted by a species' environment and shed light on a previous observation that relatively few of the mutations arising in a thermophilic archaeon were nucleotide substitutions in contrast to indels.

The rate and character of spontaneous mutation in an RNA virus.

Estimates of spontaneous mutation rates for RNA viruses are few and uncertain, most notably due to their dependence on tiny mutation reporter sequences that may not well represent the whole genome. We report here an estimate of the spontaneous mutation rate of tobacco mosaic virus using an 804-base cognate mutational target, the viral MP gene that encodes the movement protein (MP). Selection against newly arising mutants was countered by providing MP function from a transgene. The estimated genomic mutation rate was on the lower side of the range previously estimated for lytic animal riboviruses. We also present the first unbiased riboviral mutational spectrum. The proportion of base substitutions is the same as that in a retrovirus but is lower than that in most DNA-based organisms. Although the MP mutant frequency was 0.02-0.05, 35% of the sequenced mutants contained two or more mutations. Therefore, the mutation process in populations of TMV and perhaps of riboviruses generally differs profoundly from that in populations of DNA-based microbes and may be strongly influenced by a subpopulation of mutator polymerases.

Dissecting the fidelity of bacteriophage RB69 DNA polymerase: site-specific modulation of fidelity by polymerase accessory proteins.

Bacteriophage RB69 encodes a replicative B-family DNA polymerase (RB69 gp43) with an associated proofreading 3' exonuclease. Crystal structures have been determined for this enzyme with and without DNA substrates. We previously described the mutation rates and kinds of mutations produced in vivo by the wild-type (Pol(+) Exo(+)) enzyme, an exonuclease-deficient mutator variant (Pol(+) Exo(-)), mutator variants with substitutions at Tyr(567) in the polymerase active site (Pol(M) Exo(+)), and the double mutator Pol(M) Exo(-). Comparing the mutational spectra of the Pol(+) Exo(-) and Pol(+) Exo(+) enzymes revealed the patterns and efficiencies of proofreading, while Tyr(567) was identified as an important determinant of base-selection fidelity. Here, we sought to determine how well the fidelities of the same enzymes are reflected in vitro. Compared to their behavior in vivo, the three mutator polymerases exhibited modestly higher mutation rates in vitro and their mutational predilections were also somewhat different. Although the RB69 gp43 accessory proteins exerted little or no effect on total mutation rates in vitro, they strongly affected mutation rates at many specific sites, increasing some rates and decreasing others.

Experimental selection reveals a trade-off between fecundity and lifespan in the coliphage Qß.

Understanding virus evolution is key for improving ways to counteract virus-borne diseases. Results from comparative analyses have previously suggested a trade-off between fecundity and lifespan for viruses that infect the bacterium Escherichia coli (i.e. for coliphages), which, if confirmed, would define a particular constraint on the evolution of virus fecundity. Here, the occurrence of such a trade-off is investigated through a selection experiment using the coliphage Qß. Selection was applied for increased fecundity in three independent wild-type Qß populations, and the ability of the virions to remain viable outside the host was determined. The Qß life-history traits involved in the evolution of fecundity and the genetic changes associated with this evolution were also investigated. The results reveal that short-term evolution of increased fecundity in Qß was associated with decreased viability of phage virions. This trade-off apparently arose because fecundity increased at the expense of reducing the amount of resources (mainly time) invested per produced virion. Thus, the results also indicate that Qß fecundity may be enhanced through increases in the rates of adsorption to the host and progeny production. Finally, genomic sequencing of the evolved populations pinpointed sequences likely to be involved in the evolution of Qß fecundity.

Contrasting mutation rates from specific-locus and long-term mutation-accumulation procedures.

Until recently, the two predominant ways to estimate mutation rates were the specific-locus method and the mutation-accumulation (Bateman-Mukai) method. Both involve seeding a number of parallel lines from a small, genetically uniform population, growing as long as is feasible but not so long as to allow selection to perturb mutant frequencies, and sometimes using extreme bottlenecks to facilitate the retention of deleterious mutations. In the specific-locus method, mutations are selected according to their specific phenotypes and are confirmed by sequencing. In older versions of the mutation-accumulation method, the increase in variance of a quantitative fitness trait is measured and converted into a mutation rate. More recently, a variation on the mutation-accumulation method has become possible based on phenotype-blind genomic sequencing, which might (or might not) provide improved sampling breadth, usually at the expense of sample size. In a recent study, genomic sequencing was applied to Escherichia coli lines propagated for 40,000 generations and passaged daily via 5,000,000 cells. To mitigate the impact of selection, the only targets employed for rate calculations were putatively neutral synonymous mutations. The mutation rate estimate was about 6-fold lower than obtained previously with a robust specific-locus method. Here I argue that purifying selection acting to shape the strong codon preferences of E. coli is the probable cause of the lower estimate, rather than, for instance, a lower mutation rate in nature than in the laboratory.

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.