The molecular basis of adaptation is a major focus of evolutionary biology, yet the dynamic process of adaptation has been explored only piecemeal. Experimental evolution of two bacteriophage lines under strong selection led to over a dozen nucleotide changes genomewide in each replicate. At least 96 percent of the amino acid substitutions appeared to be adaptive, and half the changes in one line also occurred in the other. However, the order of these changes differed between replicates, and parallel substitutions did not reflect the changes with the largest beneficial effects or indicate a common trajectory of adaptation.

The primary impediment to formulating a general theory for adaptive evolution has been the unknown distribution of fitness effects for new beneficial mutations. By applying extreme value theory, Gillespie circumvented this issue in his mutational landscape model for the adaptation of DNA sequences, and Orr recently extended Gillespie's model, generating testable predictions regarding the course of adaptive evolution. Here we provide the first empirical examination of this model, using a single-stranded DNA bacteriophage related to phiX174, and find that our data are consistent with Orr's predictions, provided that the model is adjusted to incorporate mutation bias. Orr's work suggests that there may be generalities in adaptive molecular evolution that transcend the biological details of a system, but we show that for the model to be useful as a predictive or inferential tool, some adjustments for the biology of the system will be necessary.

High temperature inhibits the growth of the wild-type bacteriophage phiX174. Three different point mutations were identified that each recovered growth at high temperature. Two affected the major capsid protein (residues F188 and F242), and one affected the internal scaffolding protein (B114). One of the major capsid mutations (F242) is located in a beta strand that contacts B114 in the procapsid during viral maturation, whereas the other capsid mutation (F188) is involved in subunit interactions at the threefold axis of symmetry. Selective coefficients of these mutations ranged from 13.9 to 3.8 in the inhibitory, hot environment, but all mutations reduced fitness at normal temperature. The selective effect of one of the mutations (F242) was evaluated at high temperature in four different genetic backgrounds and exhibited epistasis of diminishing returns: as log fitness of the background genotype increased from -0.1 to 4.1, the fitness boost provided by the F242 mutation decreased from 3.9 to 0. 8. These results support a model in which viral fitness is bounded by an upper limit and the benefit of a mutation is scaled according to the remaining opportunity for fitness improvement in the genome.

Bacteriophage genomic evolution has been largely characterized by rampant, promiscuous horizontal gene transfer involving both homologous and nonhomologous source DNA. This pattern has emerged through study of the tailed double-stranded DNA (dsDNA) phages and is based upon a sparse sampling of the enormous diversity of these phages. The single-stranded DNA phages of the family Microviridae, including X174, appear to evolve through qualitatively different mechanisms, possibly as result of their strictly lytic lifestyle and small genome size. However, this apparent difference could reflect merely a dearth of relevant data. We sought to characterize the forces that contributed to the molecular evolution of the Microviridae and to examine the genetic structure of this single family of bacteriophage by sequencing the genomes of microvirid phage isolated on a single bacterial host. Microvirids comprised 3.5% of the detectable phage in our environmental samples, and sequencing yielded 42 new microvirid genomes. Phylogenetic analysis of the genes contained in these and five previously described microvirid phages identified three distinct clades and revealed at least two horizontal transfer events between clades. All members of one clade have a block of five putative genes that are not present in any member of the other two clades. Our data indicate that horizontal transfer does contribute to the evolution of the microvirids but is both quantitatively and qualitatively different from what has been observed for the dsDNA phages.

The distribution of fitness effects for beneficial mutations is of paramount importance in determining the outcome of adaptation. It is generally assumed that fitness effects of beneficial mutations follow an exponential distribution, for example, in theoretical treatments of quantitative genetics, clonal interference, experimental evolution, and the adaptation of DNA sequences. This assumption has been justified by the statistical theory of extreme values, because the fitnesses conferred by beneficial mutations should represent samples from the extreme right tail of the fitness distribution. Yet in extreme value theory, there are three different limiting forms for right tails of distributions, and the exponential describes only those of distributions in the Gumbel domain of attraction. Using beneficial mutations from two viruses, we show for the first time that the Gumbel domain can be rejected in favor of a distribution with a right-truncated tail, thus providing evidence for an upper bound on fitness effects. Our data also violate the common assumption that smalleffect beneficial mutations greatly outnumber those of large effect, as they are consistent with a uniform distribution of beneficial effects.

Epistatic interactions between genes and individual mutations are major determinants of the evolutionary properties of genetic systems and have therefore been well documented, but few quantitative data exist on epistatic interactions between beneficial mutations, presumably because such mutations are so much rarer than deleterious ones. We explored epistasis for beneficial mutations by constructing genotypes with pairs of mutations that had been previously identified as beneficial to the ssDNA bacteriophage ID11 and by measuring the effects of these mutations alone and in combination. We constructed 18 of the 36 possible double mutants for the nine available beneficial mutations. We found that epistatic interactions between beneficial mutations were all antagonistic—the effects of the double mutations were less than the sums of the effects of their component single mutations. We found a number of cases of decompensatory interactions, an extreme form of antagonistic epistasis in which the second mutation is actually deleterious in the presence of the first. In the vast majority of cases, recombination uniting two beneficial mutations into the same genome would not be favored by selection, as the recombinant could not outcompete its constituent single mutations. In an attempt to understand these results, we developed a simple model in which the phenotypic effects of mutations are completely additive and epistatic interactions arise as a result of the form of the phenotype-to-fitness mapping. We found that a model with an intermediate phenotypic optimum and additive phenotypic effects provided a good explanation for our data and the observed patterns of epistatic interactions.

Adaptation in haploid organisms has been extensively modeled but little tested. Using a microvirid bacteriophage (ID11), we conducted serial passage adaptations at two bottleneck sizes (10 4 and 10 6 ), followed by fitness assays and whole-genome sequencing of 631 individual isolates. Extensive genetic variation was observed including 22 beneficial, several nearly neutral, and several deleterious mutations. In the three large bottleneck lines, up to eight different haplotypes were observed in samples of 23 genomes from the final time point. The small bottleneck lines were less diverse. The small bottleneck lines appeared to operate near the transition between isolated selective sweeps and conditions of complex dynamics (e.g., clonal interference). The large bottleneck lines exhibited extensive interference and less stochasticity, with multiple beneficial mutations establishing on a variety of backgrounds. Several leapfrog events occurred. The distribution of first-step adaptive mutations differed significantly from the distribution of second-steps, and a surprisingly large number of second-step beneficial mutations were observed on a highly fit first-step background. Furthermore, few first-step mutations appeared as second-steps and second-steps had substantially smaller selection coefficients. Collectively, the results indicate that the fitness landscape falls between the extremes of smooth and fully uncorrelated, violating the assumptions of many current mutational landscape models.

Theories of adaptive molecular evolution have recently experienced significant expansion, and their predictions and assumptions have begun to be subjected to rigorous empirical testing. However, these theories focus largely on predicting the first event in adaptive evolution, the fixation of a single beneficial mutation. To address longterm adaptation it is necessary to include new assumptions, but empirical data are needed for guidance. To empirically characterize the general properties of adaptive walks, eight recently isolated relatives of the single-stranded DNA (ssDNA) bacteriophage uX174 (family Microviridae) were adapted to identical selective conditions. Three of the eight genotypes were adapted in replicate, for a total of 11 adaptive walks. We measured fitness improvement and identified the genetic changes underlying the observed adaptation. Nearly all phages were evolvable; nine of the 11 lineages showed a significant increase in fitness. However, fitness plateaued quickly, and adaptation was achieved through only three substitutions on average. Parallel evolution was rampant, both across replicates of the same genotype as well as across different genotypes, yet adaptation of replicates never proceeded through the exact same set of mutations. Despite this, final fitnesses did not vary significantly among replicates. Final fitnesses did vary significantly across genotypes but not across phylogenetic groupings of genotypes. A positive correlation was found between the number of substitutions in an adaptive walk and the magnitude of fitness improvement, but no correlation was found between starting and ending fitness. These results provide an empirical framework for future adaptation theory.