BACKGROUND
in asexual populations, adaptive evolution is fueled by novel beneficial mutations that arise and sequentially fix due to selection.
however, not all of these beneficial mutations may reach fixation; when they co-occur, competition for fixation ensues  <cit> .
those with higher selection coefficients are more likely to win the competition resulting in loss of the smaller effect mutations, a phenomenon called 'clonal interference'  <cit> .
a less severe version of this type of interference, that can decrease the efficacy of selection in sexual organisms, occurs when a locus under weak selection is tightly linked to a more strongly selected locus .
both levels of interference can significantly impact adaptive outcome by decreasing the rate of adaptation, altering levels of standing genetic variation, and favoring the evolution of sex  <cit> .
a better understanding of the effects of interference on the fate of beneficial mutations is not only important in advancing our knowledge of evolutionary processes; it also has great practical relevance in monitoring and predicting disease progression in hosts with chronic infections  <cit> .
theoretical studies of adaptation in large populations of asexual organisms show that the strength of interference is positively correlated with population size and mutation rate because mutations have a higher probability of co-occurring under these conditions  <cit> .
under the assumption that the number of available beneficial mutations is infinite, ci is predicted to: 1) slow the rate of adaptation by delaying fixation of beneficial mutations, and 2) result in a larger mean selection coefficient of fixed beneficial mutations .
kim and orr  <cit>  tested whether fixation of beneficial mutations is delayed when only a finite number of beneficial mutations are available in a given genome.
their study showed that ci can have little effect on the speed of adaptation in large populations when few beneficial mutations are available, or when the distribution of beneficial mutations has a high mean and variance.
when beneficial mutations are extremely rare, each mutation may reach fixation before a competitor arises.
likewise, when selection coefficients are very different, the transition time  may be so short that the large effect mutation reaches fixation before competitors arise.
ultimately, it is the interplay of these factors, mutation number and transition time, that determines the level of ci.
thus, our understanding of adaptation in large asexual populations could be improved by examining how parameters of the distribution of selection coefficients for beneficial mutations  influence levels of ci.
manipulating environmental conditions is a feasible experimental approach for testing effects of the distribution of beneficial mutations on levels of ci.
a more recent study of adaptation in large asexual populations distinguishes between a one-by-one fixation regime under clonal interference and a multiple mutation regime where competing genotypes may have more than one mutation  <cit> .
in contrast to the one-by-one model of clonal interference, which predicts decelerated adaptation rates and transient levels of genetic variation at large population size  and high beneficial mutation rate , the multiple mutation model predicts that the both the speed of adaptation and variation in fitness increase logarithmically with n and with μb.
consistent with predictions from the one-by-one clonal interference models  <cit> , experimental evolution studies with e. coli and vesicular stomatitis virus have shown that adaptation rates plateau, and that higher population fitness is achieved as population sizes and mutation rates increase  <cit> .
however, since these data are phenotypic, ci is not the only explanation for the observed patterns.
in larger populations, adaptive mutations will have longer transition times even in the absence of competitors, and there may be few adaptive mutations available under some conditions, especially if the first few mutations that fix leave the population stranded at a low local fitness peak.
more recent studies have provided a compelling demonstration that competition between beneficial mutations affects the probability that genotypes are lost, emphasizing the importance of multiple simultaneous mutations in determining adaptive trajectories  <cit> .
yet another recent study of adaptation in large yeast populations did not observe that the speed of adaptation decreased during adaption under conditions where multiple mutations were presumed to be pervasive, and found that variation in fitness was higher at higher nμb, supporting the multiple mutation model of adaptation  <cit> .
these more recent studies have also tested predictions of the underlying models of adaptation by observing patterns of phenotypic adaptation.
none of the empirical studies have investigated the presence of competing mutations with genomic sequence data, or examined the number of genotypes that may be competing at any given time, although holder and bull  <cit>  did observe high frequency polymorphisms during experimental evolution of bacteriophage in the family microviridae.
examining the underlying genetics of ci in large populations evolved in different environments would increase our mechanistic understanding of adaptation in large microbial populations and identify the importance of considering multiple concurrent mutations in models of adaptation.
a direct way to identify ci in experimental evolution is through genome sequence analysis.
beneficial mutations must be present simultaneously for ci to affect their fixation.
populations with mutually exclusive adaptive genotypes that coexist at relatively high frequencies indicate ci.
the bacteriophage φx <dig> is useful for identifying ci through sequence analysis and experimental evolution because it is small enough that genomic sequence data can be obtained for multiple individuals at multiple time points during the evolutionary trajectory  <cit> .
because each phage produces about  <dig> infectious progeny per  <dig> minutes, growing populations must be diluted frequently during long-term propagation.
the bottleneck effects from flask-to-flask transfers are so strong during experimental evolution of φx <dig> that there is an extremely low probability that small-effect or neutral mutations will attain appreciable frequency in the population.
when only a small number of individuals are sampled from a large phage population, there is likewise an extremely low probability of sampling an evolved genotype two or more times unless it has ascended to high frequency.
the presence of multiple high frequency genotypes in a population indicates either a selective sweep in progress or ci, depending on the nature of the differences between genotypes.
to identify the underlying genetics of ci we evolved large populations of φx <dig> and sequenced genomes of isolated genotypes from mid- and endpoint populations.
we evolved similarly large populations of phage under different levels of cacl <dig> to create benign versus harsh conditions in order to test effects of environment on the frequency of ci.
we expected that differences in cacl <dig> abundance would alter parameters of the distribution of available beneficial mutations, thus influencing levels of ci.
we assumed that by limiting the same resource to different degrees, selection would target the same trait in each treatment, but that mutations that could potentially improve the trait would be proportionately more beneficial in harsher resource conditions.
in other words, we expected that the number of available beneficial mutations  would be similar among treatments, but that the mean selection coefficient of beneficial mutations would be higher in the harsh environment since these populations would be further from their local fitness optimum.
under these assumptions, we predicted that ci would be more prevalent under higher resource levels since the smaller-effect mutations would have longer transition times, allowing for their co-occurrence.
in contrast to our assumptions, we did not find evidence of stronger selection in the harsh environment and did find evidence that there may have been fewer beneficial mutations available in the harsh environment.
thus, although our resource treatments did not have the effects we expected, they did alter the distribution of available beneficial mutations in a manner that yielded much needed insight into the effects of environment on the frequency of ci.
RESULTS
pre-adaptation to control conditions
growth rate of the lab strain increased dramatically  during the first  <dig> passages in experimental conditions, and very little  during the next  <dig>  at this point , we sequenced an isolate and found three novel mutations .
we then continued the adaptation for  <dig> more transfers  and sequenced another plaque isolate.
no additional mutations were found, and thus we concluded that the genotype was fully adapted to the control conditions such that it would not be expected to undergo further adaptive evolution under control conditions.
we measured fitness of this triple mutant at different concentrations of cacl <dig> to determine conditions where fitness was significantly decreased relative to those used in the pre-adaptation.
we used this genotype as the ancestor  to initiate the first line in each treatment.
when the control line reached the midpoint of the evolutionary trajectory, we sequenced  <dig> isolated genotypes and found that all  <dig> had an additional mutation, c788t.
in an effort to maintain a control line, we used this genotype  to initiate the second and third lines in each treatment before we knew that other adaptive mutations would be detected in the endpoint population of the first control line.
thus, the second and third lines are true replicates of each other, while the first line in each treatment was initiated from a different ancestor.
additionally, since the control lines behaved like the benign treatment, they were treated as such in the analyses.
adaptation to treatments
to examine potential differences in adaptation among the treatments with different ancestors we first conducted a two-factor anova with interaction on the relative fitness data, where ancestor and treatment were the main effects.
there was a significant difference in adaptation between the two ancestors , but no significant effect of treatment  or the interaction of ancestor and treatment , indicating that fitness increased by similar amounts among the resource treatments.
since the control and benign environments were very similar treatments in terms of both the level of cacl <dig> and the starting fitnesses of the ancestor , we pooled these two treatments to increase power for testing differences in adaptation between the benign and harsh environments.
adaptation was not significantly different between the benign and harsh treatments for anc .
on the other hand, adaptation of anc <dig> was marginally different between the two environments .
this difference was statistically significant if we analyzed the data by a nested anova that accounted for potential variation due to the different lines within treatments .
overall, the data showed that there was no significant difference in the amount of adaptation between the benign and harsh environments for anc, but that the amount of adaptation may have been slightly higher in the benign environment for anc <dig> .
due to these differences and the finding that there was no significant variation from lines within treatments, we analyzed fitness data from each ancestor separately, but pooled data for the benign treatments, and data for anc <dig> lines within each benign and harsh treatment for subsequent analyses of adaptation.
treatments
treatment: control =  <dig>  mm cacl <dig>  benign =  <dig>  mm cacl <dig>  harsh =  <dig>  mm cacl <dig>
no and n30: mean phage population sizes at the start and end of each flask growth period.
no/n30: bottleneck size; lower numbers indicate lines where beneficial mutations had higher extinction probabilities from bottlenecks.
wo: mean fitness of the ancestor in treatment conditions .
we: mean fitness of the endpoint populations in treatment conditions .
relw: mean fitness of the evolved population relative to ancestor .
s: mean selection coefficient; a measure of selection strength.
relw/we: fitness change as a proportion of total fitness; a measure of selection strength.
"harsh" treatments are shown in bold.
selection strength
selection strength is determined by parameters of the distribution of beneficial mutations.
if the mean selection coefficient  is large, then selection is strong.
thus, to evaluate differences in selection strength between treatments, we first compared the average selection coefficient of mutations by dividing the relative fitness by the average number of mutations per genome, which was weighted by the frequency of genotypes with each particular number of mutations .
this measure takes into account both the frequency and effect of observed mutations to estimate the average selection coefficient, but makes the simplifying assumption that mutations have equal selection coefficients.
we found that selection strength was significantly higher in the harsh environment for anc , but that oppositely, selection strength was stronger in the benign environment for anc <dig> .
another dimension of selection strength is the amount of fitness gain as a proportion of total fitness.
for example, if two evolved populations each increase in fitness by  <dig>  but one population has a total fitness of  <dig> and the other population has a total fitness of  <dig>  then mutations occurring in the former case cause proportionately more fitness gain which could indicate stronger selection.
thus, as a second method of evaluating selection strength, we compared the proportional fitness gain between the two treatments for each ancestor.
by this measure, selection strength was again stronger in the harsh environment for anc , but was not different between the two treatments for anc <dig> .
overall, there was evidence that selection was strongest for anc in the harsh environment.
however, since only a single mutation was found for anc in the harsh environment, it is possible that we observed this pattern by chance, and thus we did not consider this a strong result.
nature of the genetic changes
a. determining which mutations are adaptive
ci is a process concerning adaptive mutations, yet we did not measure the fitness effects for most of the individual mutations observed in this experiment.
in order to identify populations undergoing ci, we used sequence data  and evolutionary theory to infer that most, if not all, of the mutations were adaptive.
we note that the bottleneck sizes in our experiment almost completely excluded mutations with selection coefficients less than  <dig>  , suggesting that: 1) our design was biased towards sampling large-effect mutations, and 2) neutral mutations had negligible probabilities of rising to appreciable frequencies.
however, in an effort to be conservative, we considered whether each mutation showed other evidence of being adaptive .
a common criterion for identifying adaptive mutations from sequence data is parallel evolution  <cit> .
we found that ten of the  <dig> different mutations were observed in multiple lines in this study, and seven of these ten parallel mutations were observed in multiple isolates in at least one line, suggesting that they had attained high frequency.
of the two parallel mutations not observed at high frequency, one was observed six times in previous studies.
of the  <dig> mutations observed in only one line in this study, nine were observed in more than one isolate .
although it is possible that a mutation could reach high frequency by hitchhiking with a large-effect beneficial mutation, the substitution dynamics for eight of these nine high frequency mutations indicate that hitchhiking does not explain their ascent to high frequency; hitchhiking could not be ruled out for g2973a.
thus, there is strong evidence that  <dig> of the  <dig> different mutations are adaptive.
there are two mutations where evidence for their adaptive nature is not as strong – g2973a and c581t .
in the case of g2973a, this mutation has been observed previously under similar conditions , and it is in a codon adjacent to mutation c2971t, which showed strong evidence of being adaptive in the current study.
likewise, the mutation c581t was observed previously under similar conditions , and it results in a non-synonymous substitution that is two amino acids away from one that occurred independently in two lines of this experiment , as well as being two away from a site known to have a large effect on lysis regulation .
there is no evidence that the remaining four of the  <dig> mutations were adaptive .
interestingly, neither the adaptive nature of the two mutations where evidence is weaker, nor the adaptive nature of these four mutations, affects our conclusions regarding the presence of ci in any replicate.
a numbers indicate nucleotide positions that are different from ancestor; rows list the mutational position for a single genotype.
nd = not done; bold = parallel mutations.
b total number of genotypes over total number of isolates sequenced.
c ci status of mid- and endpoint populations; y means that multiple mutually exclusive genotypes are present .
y = strong ci, where multiple isolates of competing mutations were observed; y = weak ci, where only one competing isolate was observed.
d number of adaptive mutations in endpoint populations as determined from criteria in table  <dig>  data were used for mann-whitney u-test of ho: the number of adaptive mutations sampled in control/benign endpoint populations is similar to the number of adaptive mutations sampled in the harsh populations .
a genome position of nucleotide mutation.
b c = control; b = benign; h = harsh.
c gene product  and amino acid position of mutation.
d effect of the mutation at the amino acid level .
gene product c is involved in dna replication, d is the external scaffolding protein, e is the lysis protein, f is the major capsid protein, g is the minor capsid protein, h is the pilot protein, j is the dna binding protein, and k is of unknown function.
because φx <dig> has overlapping reading frames, some changes affect more than one protein.
e evidence that a particular change is adaptive includes: parallel = occurrence in multiple lines of this study; high frequency = occurrence in multiple isolates in a population in this study, indicating the change was at high frequency in the population and that there was no evidence of hitchhiking; previous  = occurrence in previous evolution experiments  ; measured = measurement of fitness effects.
none = no independent evidence that the mutation is adaptive.
of the mutations considered to be adaptive by the criteria above, hitchhiking could not be ruled out for  <dig> and  <dig>  however, both of these mutations have occurred previously under very similar experimental conditions and there are additional indicators that these mutations are adaptive:  <dig> affects an amino acid two away from another that is adaptive in this experiment , and is also two sites away from a site known to have a large effect on lysis regulation .
mutation  <dig> at h <dig> is adjacent to site h <dig> which showed strong evidence of being adaptive.
parentheses indicate synonymous substitutions.
b. number of adaptive mutations per treatment
another parameter of the distribution of beneficial mutations that can affect levels of ci is the overall frequency of adaptive mutations that are available.
if fewer adaptive mutations are available, then the waiting time is increased which could decrease levels of ci.
to investigate whether there were differences in the number of mutations available under the two environmental treatments, we compared the number of adaptive mutations identified in populations from the harsh and benign environments by a mann-whitney u-test.
for this analysis, we pooled data from both ancestors since only one population was examined for anc in the harsh condition, and we excluded mutations that showed no evidence of being adaptive.
there were significantly more adaptive mutations in endpoint populations of the control/benign treatments relative to those in the harsh treatments .
c. site of mutations
in the benign environment, adaptive mutations identified in the endpoint populations were between nucleotide positions  <dig> to  <dig>  spanning several genes, and were classified into four functional groups .
on the other hand, mutations observed in the harsh treatment were all in gene h between positions  <dig> and  <dig>  this pattern was not absolute.
the c788t change at e <dig> in anc <dig> was beneficial in the harsh environment even though it did not occur there, and a3340g at h <dig> was at high frequency in the midpoint of control line  <dig>  but lost at the end.
nevertheless, it appeared that changes affecting different viral functions dominated adaptation under benign and harsh treatments , and in the benign environment, a larger number of sites in several genes was able to contribute adaptive solutions, while in the harsh environment, solutions appeared to be more limited.
a numbers are the frequency that mutations from a particular functional category were observed in independent evolution lines under the indicated treatments.
only adaptive mutations  from endpoint populations are shown.
frequency of ci
we evaluated the frequency of ci by determining the number of endpoint populations where ci was present or absent and comparing these frequencies between treatments by a fisher exact test .
again, we pooled data from both ancestors in order to increase sample size and considered the control as a benign treatment.
we found that significantly more populations showed evidence of ci in the benign environments relative to the harsh environment .
a pops with ci: total number of populations in each treatment  showing evidence of ci ; pops with no ci: total number of populations not showing evidence of ci.
b number of endpoint populations with or without evidence of ci pooled across both ancestors.
c ho: the frequency of populations showing ci is not higher in benign conditions.
another pattern that would support the occurrence of ci under our experimental conditions is the observation of mutations in a midpoint population that were not seen in the endpoint population of the same line.
we therefore examined variation at the midpoint of some populations where the endpoint strongly indicated ci and in line  <dig> of the harsh treatment where variation  was present in the endpoint population.
ci was evident in three of the four midpoint populations from control and benign treatments, and all three of these had mutations that were not seen among the endpoint isolates.
in two of these cases, the lost genotypes were at very high frequency .
in the harsh environment, there was no evidence of ci in the midpoint population, and although there was one mutation that was not observed in the endpoint population, that mutation was a synonymous mutation with no independent evidence of being adaptive.
to examine whether mutations that were lost between the midpoint and endpoint populations could have been out-competed, we measured the selection coefficients of a mutation that showed evidence of being lost due to ci, and of the mutation that may have out-competed it .
in the midpoint population of control line  <dig>  only c788t was detected, while in the midpoint population of benign line  <dig> there were four genotypes, including c788t.
by the endpoints, c788t was present in all genotypes of control line  <dig> but was not detected in benign line  <dig>  and a1611g was detected in all genotypes in the endpoint populations of both treatments.
apparently, a1611g did not interfere with the fixation of c788t in the control line, but did interfere in the benign line.
if ci were responsible for this outcome, then in the benign line, both mutations should be adaptive and a1611g should have a significant advantage over c788t.
indeed, both mutations had a significant selective advantage over their ancestor , and the a1611g mutation had a significant selective advantage over the c788t mutation , supporting that c788t was lost due to ci by a1611g in the benign line.
frequency of multiple mutations
evidence supporting the multiple-mutations model versus the one-by-one model of clonal interference would be the occurrence of more than one change from the fixed genotype.
we found four populations where a single isolate differed from the  fixed genotype by more than one change: benign  <dig> midpoint, isolate  <dig>  1702; control  <dig> endpoint isolate  <dig>   <dig>  1216; control  <dig> midpoint, isolate  <dig>  581; and harsh  <dig> midpoint, isolate  <dig>   <dig>   in two of these cases , the changes were synonymous in gene f and there is no independent evidence to suggest that the additional change was adaptive.
in the two cases where there is evidence that the additional changes were adaptive , they were seen only in the midpoint populations.
additionally, since the control lines behaved like the benign treatment, they were treated as such in the analyses.
discussion
theoretical models and empirical data show that increased population sizes and mutation rates can slow the rate of adaptation and decrease the probability that beneficial mutations will reach fixation  <cit> .
phenotypic data from the empirical studies show that rates of adaptation do not increase linearly with population size, but it is not clear whether this pattern is actually due to competing mutations or simply because each mutation will have a longer transition time in larger populations.
it is also not clear how many genotypes might be competing at any given time, which is an important parameter for studying properties of adaptation and predicting adaptive walks.
we evolved nine lines of phage and sequenced genomes of multiple individuals from populations in these lines to examine the underlying genetics of ci.
we found that mutations with relatively large selection coefficients can be out-competed and that there were as many as three competing genotypes in many of the populations where resource was abundant, even though the regular strong bottlenecks inherent to our flask-to-flask transfer protocol strongly decreased our ability to observe most mutations.
this suggested both that beneficial mutations and ci were remarkably prevalent.
our sequence data did not reveal any ci in the harsh treatment, which emphasized that the frequency of ci can vary among similarly large populations due to environmental conditions.
adaptation under different levels of resource
the resource treatments were intended to alter parameters of the distribution of beneficial mutations in order to evaluate environmental effects on levels of ci in similarly large evolving populations.
while we found evidence that there were differences between lines in the available beneficial mutations, these differences were not as we predicted.
we expected that the ancestor would be further from its fitness optimum under harsh resource conditions, which would result in a larger mean selection coefficient for beneficial mutations.
however, there was insufficient evidence that the mean selection coefficient of beneficial mutations was higher in the harsh environment for anc, and it even appeared to be lower in the harsh treatment for anc <dig>  apparently, the c788t mutation in anc increased fitness in the harsh environment , such that selection was weaker on anc <dig> than on anc in the harsh environment.
also, although we expected that the variable resource levels would select for improvement on the same trait, the sequence data suggested that the stronger selection was on different targets between the two resource treatments, and possibly more genetic solutions were available in the benign environment.
thus, it appeared that the distribution of beneficial mutations differed in overall frequency of adaptive mutations between the resource treatments, such that more mutations were available in the benign environment.
factors affecting the frequency of ci
when considering the overall frequency of ci for both ancestors, we observed ci in most populations in the benign environments but not in any populations in the harsh environments .
while this result was consistent with our original prediction, it was not attributed to a larger mean selection coefficient in the harsh environment as we had predicted.
for example, selection was possibly weaker for anc <dig> in the harsh environment, yet no ci was observed in any of those populations.
additionally, although we found that the selection coefficient of the mutation that occurred in the harsh environment of line  <dig> was very high, it was the only mutation observed, and thus chance cannot be out-ruled as an explanation for the difference in selection strength between the treatments in line  <dig>  however, there was evidence of an overall trend of fewer adaptive mutations in the harsh condition , suggesting that the lower levels of ci in the harsh environment may have been due to the lower number of available adaptive solutions, which would increase the waiting time for beneficial mutations and decrease their co-occurrence  <cit> .
interestingly, even though the frequency of ci appeared to differ among the two resource treatments, their overall increase in fitness was similar, as were the population sizes among experimental lines.
this result emphasizes two important points: 1) that ci does not necessarily decrease rates of population adaptation, and 2) that levels of ci are not necessarily predictable from population size.
regarding adaptation rate, it is indeed possible that ci acts to increase population adaptation by maintaining higher levels of genetic variation.
for example, if a very high fitness peak can only be achieved through a specific ordering of substitutions, delayed fixation time from ci could increase the probability that primary mutations in the order are maintained long enough for the subsequent necessary mutations to arise.
in the absence of ci, fixation of alternate mutations may be too rapid, occurring before the critical primary mutations, which could result in the population being stranded on a lower fitness peak than what could have been realized through a different adaptive walk.
in this manner, the decreased rates of adaptation that are predicted to occur from competing mutations, could be outweighed by larger final fitness gains due to the more thorough exploration of sequence space that occurs in populations undergoing ci.
this could partly explain our observation of similar levels of adaptation, but different levels of ci, among the two resource treatments.
there appeared to be fewer adaptive solutions available in the harsh treatments since fewer different mutations were observed and most occurred in a single gene.
all substitutions in the endpoint populations of the harsh treatment were in gph, a capsid protein involved in guiding phage dna into host cells  <cit> , while those detected in the control and benign environments occurred in several different genes, including the gene encoding the major capsid protein , as well as in genomic regions that are putatively involved in regulation of lysis timing  <cit> .
gpf is involved in numerous processes throughout the phage life cycle including attachment to hosts, catalysis of eclipse, and assembly  <cit> .
the apparently larger number of genetic solutions in the benign environment relative to in the harsh environment suggested that: 1) selection acted on a broader range of phenotypes in the benign treatments, and/or 2) there were more genetic solutions to the types of phenotypes under selection in the benign treatments.
in the harsh environment where ca+ <dig> was very limiting, there should be strong specific selection on capsid proteins to either be able to inject dna with fewer ca+ <dig> molecules or to become more efficient at acquiring ca+ <dig> since genome injection is crucial to initiating infection.
in other words, in the harsh environment, improvement of this trait must come before other potential renovations to growth rate.
in contrast, in the benign environments where growth rates were already high at the outset, modification of any phenotype affecting growth rate  could be effective at increasing fitness.
for example, one strong source of selection inherent to all treatments is bottleneck survival.
faster lysis is a conceivable adaptive response under these conditions: phage genotypes with the highest number of free phage immediately before transfer  will have the highest probability of being transferred to the next flask, and these genotypes will go on to initiate the most individual infections in the subsequent flask.
thus, in the benign conditions adaptation could occur through any modification that produced a high number of phage by the end of each flask growth period whereas in the harsh environment selection specifically targeted injection rates before any other process in the phage life cycle, suggesting that the higher levels of ci in the benign environment could have been due to a larger number of adaptive targets under these conditions.
experimental design issues
a. bottleneck effects
since phage had higher growth rates in the benign environment, bottlenecks to a smaller population size were required to maintain population sizes similar to those in the harsh treatments .
thus, beneficial mutations that arose in the benign condition had a higher probability of extinction due to bottlenecks  <cit> , suggesting that our design was slightly biased towards finding mutations in the harsh condition.
interestingly, we found fewer adaptive mutations and lower levels of ci in the harsh treatment where the extinction probabilities of beneficial mutations due to bottlenecks were lower.
this demonstrates that the difference in ci between the two treatments may have been even more dramatic  if bottleneck sizes were equal.
thus, our approach to testing the hypothesis that ci would be stronger in benign environments was conservative.
b. growth rate differences
the combined effect of differences in growth rates for phage in the two treatments and frequent, regular, strong bottlenecking is an experimental design factor that could have increased the extinction probability of mutations in the harsh environment.
previous work has shown that phage with lower growth rates have increased burst times  <cit> , suggesting that replication is slower.
if the higher growth rates of phage in the benign environment translate to faster replication relative to that in the harsh environment, then mutations would arise earlier in flask growth.
from the work of wahl and gerrish  <cit> , we found that mutations arising in the last five minutes of flask growth have a greater than 99% chance of being lost due to sampling, regardless of their selection coefficients, while those arising in the first five minutes are likely to survive the bottleneck if their selection coefficients are high.
thus, the lower growth rates of phage in the harsh environment could have increased the extinction probability of beneficial mutations if most were produced too late during flask growth to survive the bottleneck.
if present, this effect could explain the lower levels of ci observed in the harsh condition.
however, it is also possible that the less severe bottleneck in the harsh condition, which decreased extinction probabilities , acted to negate potential effects from low growth rates.
implications of ci for predicting adaptation
in four of the five lines where we sequenced isolates from the midpoint populations, we found genotypes that were no longer present in the endpoint population, which we refer to as 'lost.'
furthermore, from measuring the selection coefficients of one of these lost individuals we found that mutations with relatively large selection coefficients can be lost.
these data demonstrated two important points that are relevant to understanding and predicting adaptation.
first, the probability that a beneficial mutation is lost due to a deterministic evolutionary phenomenon such as ci depends partly on stochasticity from the time that other mutations of similar effect arise.
in populations that are subject to strong sampling effects, it is possible for a small-effect mutation to reach fixation in the presence of a larger one, as long as the smaller-effect mutation is near fixation before the larger-effect mutation arises, such that the larger-effect mutation is more strongly affected by the stochastic sampling.
second, the fact that relatively large-effect beneficial mutations were lost in competition supports the finding that beneficial mutations may be much more frequent than believed  <cit> .
many beneficial mutations are often not seen because they can only be sampled in environments where selection coefficients are similar among mutations, or where there are no mutations with very large selection coefficients.
this emphasizes the importance in measuring adaptation in numerous environments in order to develop predictive models and understand evolutionary processes.
the co-occurrence of high frequency beneficial mutations has been previously observed during experimental evolution of microvirid phages in batch culture .
we found evidence for ci in five out of nine populations even though all mutations had high extinction probabilities, suggesting that ci can be remarkably common in large bacteriophage populations.
while in most cases competing genotypes differed from the  fixed genotype by only by a single adaptive mutation, which supports predictions of the one-by-one model of clonal interference, we found limited support for the multiple mutation model.
there were four cases where a single isolate differed from the fixed genotype by more than one change, although in two of these cases there is no evidence that the additional change is adaptive.
in no case did these genotypes carrying multiple mutations dominate in the endpoint populations.
taken together, our data and the recent study in yeast emphasize that concurrent mutations are an important determinant of adaptive trajectories in large microbial populations.
nevertheless, sequence-based models of adaptation that assume selection is strong relative to mutation , such that one beneficial mutation is fixed before the next occurs, have been considered acceptable for describing properties of adaptive walks in asexual organisms  <cit> .
our experimental data suggest that sswm conditions may often be unrealistic in large microbial populations, and thus that models of adaptation relying on this assumption may make misleading predictions regarding the properties of adaptive walks in large microbial populations.
models of adaptation that include concurrent mutations and multiple concurrent mutations on individual genotypes  <cit>  may provide more realistic frameworks for investigating the properties of adaptation in large microbial populations.
CONCLUSIONS
we verified with sequence data that ci occurs in large virus populations and that several genotypes may co-occur at very high frequencies.
we intended to test the prediction that weak selection in benign environments would result in more ci on average.
while our data showed that ci was more frequent in benign treatments, this difference was not attributed to weaker selection.
it appeared that the higher frequency of ci in the control and benign environments was partly due to a higher number of potential adaptive solutions in those environments, relative to that in the harsh treatment, which would have decreased the waiting time for beneficial mutations.
at the very least, our results showed that environmental conditions can alter the frequency of ci among similarly large populations, which demonstrates that effects of ci are not necessarily predictable from population size.
another important insight from our genetic data was that selection appeared to target different phenotypes in treatments where the same continuous resource variable was decreased to different degrees.
finally, we observed that the fate of mutations with relatively large selection coefficients depends on the time that other large-effect mutations arise, suggesting that these stochastic population-level processes are significant properties of adaptive landscapes, as predicted by theoretical models.

 in asexual populations, adaptive evolution is fueled by novel beneficial mutations that arise and sequentially fix due to selection. however, not all of these beneficial mutations may reach fixation; when they co-occur, competition for fixation ensues. those with higher selection coefficients are more likely to win the competition resulting in loss of the smaller effect mutations, a phenomenon called'clonal interference'. a less severe version of this type of interference, that can decrease the efficacy of selection in sexual organisms, occurs when a locus under weak selection is tightly linked to a more strongly selected locus. both levels of interference can significantly impact adaptive outcome by decreasing the rate of adaptation, altering levels of standing genetic variation, and favoring the evolution of sex. a better understanding of the effects of interference on the fate of beneficial mutations is not only important in advancing our knowledge of evolutionary processes; it also has great practical relevance in monitoring and predicting disease progression in hosts with chronic infections. theoretical studies of adaptation in large populations of asexual organisms show that the strength of interference is positively correlated with population size and mutation rate because mutations have a higher probability of co-occurring under these conditions. under the assumption that the number of available beneficial mutations is infinite, ci is predicted to: 1) slow the rate of adaptation by delaying fixation of beneficial mutations, and 2) result in a larger mean selection coefficient of fixed beneficial mutations. kim and orr tested whether fixation of beneficial mutations is delayed when only a finite number of beneficial mutations are available in a given genome. their study showed that ci can have little effect on the speed of adaptation in large populations when few beneficial mutations are available, or when the distribution of beneficial mutations has a high mean and variance. when beneficial mutations are extremely rare, each mutation may reach fixation before a competitor arises. likewise, when selection coefficients are very different, the transition time may be so short that the large effect mutation reaches fixation before competitors arise. ultimately, it is the interplay of these factors, mutation number and transition time, that determines the level of ci. thus, our understanding of adaptation in large asexual populations could be improved by examining how parameters of the distribution of selection coefficients for beneficial mutations influence levels of ci. manipulating environmental conditions is a feasible experimental approach for testing effects of the distribution of beneficial mutations on levels of ci. a more recent study of adaptation in large asexual populations distinguishes between a one-by-one fixation regime under clonal interference and a multiple mutation regime where competing genotypes may have more than one mutation. in contrast to the one-by-one model of clonal interference, which predicts decelerated adaptation rates and transient levels of genetic variation at large population size and high beneficial mutation rate, the multiple mutation model predicts that the both the speed of adaptation and variation in fitness increase logarithmically with n and with μb. consistent with predictions from the one-by-one clonal interference models, experimental evolution studies with e. coli and vesicular stomatitis virus have shown that adaptation rates plateau, and that higher population fitness is achieved as population sizes and mutation rates increase. however, since these data are phenotypic, ci is not the only explanation for the observed patterns. in larger populations, adaptive mutations will have longer transition times even in the absence of competitors, and there may be few adaptive mutations available under some conditions, especially if the first few mutations that fix leave the population stranded at a low local fitness peak. more recent studies have provided a compelling demonstration that competition between beneficial mutations affects the probability that genotypes are lost, emphasizing the importance of multiple simultaneous mutations in determining adaptive trajectories. yet another recent study of adaptation in large yeast populations did not observe that the speed of adaptation decreased during adaption under conditions where multiple mutations were presumed to be pervasive, and found that variation in fitness was higher at higher nμb, supporting the multiple mutation model of adaptation. these more recent studies have also tested predictions of the underlying models of adaptation by observing patterns of phenotypic adaptation. none of the empirical studies have investigated the presence of competing mutations with genomic sequence data, or examined the number of genotypes that may be competing at any given time, although holder and bull did observe high frequency polymorphisms during experimental evolution of bacteriophage in the family microviridae. examining the underlying genetics of ci in large populations evolved in different environments would increase our mechanistic understanding of adaptation in large microbial populations and identify the importance of considering multiple concurrent mutations in models of adaptation. a direct way to identify ci in experimental evolution is through genome sequence analysis. beneficial mutations must be present simultaneously for ci to affect their fixation. populations with mutually exclusive adaptive genotypes that coexist at relatively high frequencies indicate ci. the bacteriophage φx is useful for identifying ci through sequence analysis and experimental evolution because it is small enough that genomic sequence data can be obtained for multiple individuals at multiple time points during the evolutionary trajectory. because each phage produces about infectious progeny per minutes, growing populations must be diluted frequently during long-term propagation. the bottleneck effects from flask-to-flask transfers are so strong during experimental evolution of φx that there is an extremely low probability that small-effect or neutral mutations will attain appreciable frequency in the population. when only a small number of individuals are sampled from a large phage population, there is likewise an extremely low probability of sampling an evolved genotype two or more times unless it has ascended to high frequency. the presence of multiple high frequency genotypes in a population indicates either a selective sweep in progress or ci, depending on the nature of the differences between genotypes. to identify the underlying genetics of ci we evolved large populations of φx and sequenced genomes of isolated genotypes from mid- and endpoint populations. we evolved similarly large populations of phage under different levels of cacl to create benign versus harsh conditions in order to test effects of environment on the frequency of ci. we expected that differences in cacl abundance would alter parameters of the distribution of available beneficial mutations, thus influencing levels of ci. we assumed that by limiting the same resource to different degrees, selection would target the same trait in each treatment, but that mutations that could potentially improve the trait would be proportionately more beneficial in harsher resource conditions. in other words, we expected that the number of available beneficial mutations would be similar among treatments, but that the mean selection coefficient of beneficial mutations would be higher in the harsh environment since these populations would be further from their local fitness optimum. under these assumptions, we predicted that ci would be more prevalent under higher resource levels since the smaller-effect mutations would have longer transition times, allowing for their co-occurrence. in contrast to our assumptions, we did not find evidence of stronger selection in the harsh environment and did find evidence that there may have been fewer beneficial mutations available in the harsh environment. thus, although our resource treatments did not have the effects we expected, they did alter the distribution of available beneficial mutations in a manner that yielded much needed insight into the effects of environment on the frequency of ci.  pre-adaptation to control conditions growth rate of the lab strain increased dramatically during the first passages in experimental conditions, and very little during the next at this point, we sequenced an isolate and found three novel mutations. we then continued the adaptation for more transfers and sequenced another plaque isolate. no additional mutations were found, and thus we concluded that the genotype was fully adapted to the control conditions such that it would not be expected to undergo further adaptive evolution under control conditions. we measured fitness of this triple mutant at different concentrations of cacl to determine conditions where fitness was significantly decreased relative to those used in the pre-adaptation. we used this genotype as the ancestor to initiate the first line in each treatment. when the control line reached the midpoint of the evolutionary trajectory, we sequenced isolated genotypes and found that all had an additional mutation, c788t. in an effort to maintain a control line, we used this genotype to initiate the second and third lines in each treatment before we knew that other adaptive mutations would be detected in the endpoint population of the first control line. thus, the second and third lines are true replicates of each other, while the first line in each treatment was initiated from a different ancestor. additionally, since the control lines behaved like the benign treatment, they were treated as such in the analyses. adaptation to treatments to examine potential differences in adaptation among the treatments with different ancestors we first conducted a two-factor anova with interaction on the relative fitness data, where ancestor and treatment were the main effects. there was a significant difference in adaptation between the two ancestors, but no significant effect of treatment or the interaction of ancestor and treatment, indicating that fitness increased by similar amounts among the resource treatments. since the control and benign environments were very similar treatments in terms of both the level of cacl and the starting fitnesses of the ancestor, we pooled these two treatments to increase power for testing differences in adaptation between the benign and harsh environments. adaptation was not significantly different between the benign and harsh treatments for anc. on the other hand, adaptation of anc was marginally different between the two environments. this difference was statistically significant if we analyzed the data by a nested anova that accounted for potential variation due to the different lines within treatments. overall, the data showed that there was no significant difference in the amount of adaptation between the benign and harsh environments for anc, but that the amount of adaptation may have been slightly higher in the benign environment for anc. due to these differences and the finding that there was no significant variation from lines within treatments, we analyzed fitness data from each ancestor separately, but pooled data for the benign treatments, and data for anc lines within each benign and harsh treatment for subsequent analyses of adaptation. treatments treatment: control = mm cacl benign = mm cacl harsh = mm cacl no and n30: mean phage population sizes at the start and end of each flask growth period. no/n30: bottleneck size; lower numbers indicate lines where beneficial mutations had higher extinction probabilities from bottlenecks. wo: mean fitness of the ancestor in treatment conditions. we: mean fitness of the endpoint populations in treatment conditions. relw: mean fitness of the evolved population relative to ancestor. s: mean selection coefficient; a measure of selection strength. relw/we: fitness change as a proportion of total fitness; a measure of selection strength. "harsh" treatments are shown in bold. selection strength selection strength is determined by parameters of the distribution of beneficial mutations. if the mean selection coefficient is large, then selection is strong. thus, to evaluate differences in selection strength between treatments, we first compared the average selection coefficient of mutations by dividing the relative fitness by the average number of mutations per genome, which was weighted by the frequency of genotypes with each particular number of mutations. this measure takes into account both the frequency and effect of observed mutations to estimate the average selection coefficient, but makes the simplifying assumption that mutations have equal selection coefficients. we found that selection strength was significantly higher in the harsh environment for anc, but that oppositely, selection strength was stronger in the benign environment for anc. another dimension of selection strength is the amount of fitness gain as a proportion of total fitness. for example, if two evolved populations each increase in fitness by but one population has a total fitness of and the other population has a total fitness of then mutations occurring in the former case cause proportionately more fitness gain which could indicate stronger selection. thus, as a second method of evaluating selection strength, we compared the proportional fitness gain between the two treatments for each ancestor. by this measure, selection strength was again stronger in the harsh environment for anc, but was not different between the two treatments for anc. overall, there was evidence that selection was strongest for anc in the harsh environment. however, since only a single mutation was found for anc in the harsh environment, it is possible that we observed this pattern by chance, and thus we did not consider this a strong result. nature of the genetic changes a. determining which mutations are adaptive ci is a process concerning adaptive mutations, yet we did not measure the fitness effects for most of the individual mutations observed in this experiment. in order to identify populations undergoing ci, we used sequence data and evolutionary theory to infer that most, if not all, of the mutations were adaptive. we note that the bottleneck sizes in our experiment almost completely excluded mutations with selection coefficients less than, suggesting that: 1) our design was biased towards sampling large-effect mutations, and 2) neutral mutations had negligible probabilities of rising to appreciable frequencies. however, in an effort to be conservative, we considered whether each mutation showed other evidence of being adaptive. a common criterion for identifying adaptive mutations from sequence data is parallel evolution. we found that ten of the different mutations were observed in multiple lines in this study, and seven of these ten parallel mutations were observed in multiple isolates in at least one line, suggesting that they had attained high frequency. of the two parallel mutations not observed at high frequency, one was observed six times in previous studies. of the mutations observed in only one line in this study, nine were observed in more than one isolate. although it is possible that a mutation could reach high frequency by hitchhiking with a large-effect beneficial mutation, the substitution dynamics for eight of these nine high frequency mutations indicate that hitchhiking does not explain their ascent to high frequency; hitchhiking could not be ruled out for g2973a. thus, there is strong evidence that of the different mutations are adaptive. there are two mutations where evidence for their adaptive nature is not as strong – g2973a and c581t. in the case of g2973a, this mutation has been observed previously under similar conditions, and it is in a codon adjacent to mutation c2971t, which showed strong evidence of being adaptive in the current study. likewise, the mutation c581t was observed previously under similar conditions, and it results in a non-synonymous substitution that is two amino acids away from one that occurred independently in two lines of this experiment, as well as being two away from a site known to have a large effect on lysis regulation. there is no evidence that the remaining four of the mutations were adaptive. interestingly, neither the adaptive nature of the two mutations where evidence is weaker, nor the adaptive nature of these four mutations, affects our conclusions regarding the presence of ci in any replicate. a numbers indicate nucleotide positions that are different from ancestor; rows list the mutational position for a single genotype. nd = not done; bold = parallel mutations. b total number of genotypes over total number of isolates sequenced. c ci status of mid- and endpoint populations; y means that multiple mutually exclusive genotypes are present. y = strong ci, where multiple isolates of competing mutations were observed; y = weak ci, where only one competing isolate was observed. d number of adaptive mutations in endpoint populations as determined from criteria in table data were used for mann-whitney u-test of ho: the number of adaptive mutations sampled in control/benign endpoint populations is similar to the number of adaptive mutations sampled in the harsh populations. a genome position of nucleotide mutation. b c = control; b = benign; h = harsh. c gene product and amino acid position of mutation. d effect of the mutation at the amino acid level. gene product c is involved in dna replication, d is the external scaffolding protein, e is the lysis protein, f is the major capsid protein, g is the minor capsid protein, h is the pilot protein, j is the dna binding protein, and k is of unknown function. because φx has overlapping reading frames, some changes affect more than one protein. e evidence that a particular change is adaptive includes: parallel = occurrence in multiple lines of this study; high frequency = occurrence in multiple isolates in a population in this study, indicating the change was at high frequency in the population and that there was no evidence of hitchhiking; previous = occurrence in previous evolution experiments; measured = measurement of fitness effects. none = no independent evidence that the mutation is adaptive. of the mutations considered to be adaptive by the criteria above, hitchhiking could not be ruled out for and however, both of these mutations have occurred previously under very similar experimental conditions and there are additional indicators that these mutations are adaptive: affects an amino acid two away from another that is adaptive in this experiment, and is also two sites away from a site known to have a large effect on lysis regulation. mutation at h is adjacent to site h which showed strong evidence of being adaptive. parentheses indicate synonymous substitutions. b. number of adaptive mutations per treatment another parameter of the distribution of beneficial mutations that can affect levels of ci is the overall frequency of adaptive mutations that are available. if fewer adaptive mutations are available, then the waiting time is increased which could decrease levels of ci. to investigate whether there were differences in the number of mutations available under the two environmental treatments, we compared the number of adaptive mutations identified in populations from the harsh and benign environments by a mann-whitney u-test. for this analysis, we pooled data from both ancestors since only one population was examined for anc in the harsh condition, and we excluded mutations that showed no evidence of being adaptive. there were significantly more adaptive mutations in endpoint populations of the control/benign treatments relative to those in the harsh treatments. c. site of mutations in the benign environment, adaptive mutations identified in the endpoint populations were between nucleotide positions to spanning several genes, and were classified into four functional groups. on the other hand, mutations observed in the harsh treatment were all in gene h between positions and this pattern was not absolute. the c788t change at e in anc was beneficial in the harsh environment even though it did not occur there, and a3340g at h was at high frequency in the midpoint of control line but lost at the end. nevertheless, it appeared that changes affecting different viral functions dominated adaptation under benign and harsh treatments, and in the benign environment, a larger number of sites in several genes was able to contribute adaptive solutions, while in the harsh environment, solutions appeared to be more limited. a numbers are the frequency that mutations from a particular functional category were observed in independent evolution lines under the indicated treatments. only adaptive mutations from endpoint populations are shown. frequency of ci we evaluated the frequency of ci by determining the number of endpoint populations where ci was present or absent and comparing these frequencies between treatments by a fisher exact test. again, we pooled data from both ancestors in order to increase sample size and considered the control as a benign treatment. we found that significantly more populations showed evidence of ci in the benign environments relative to the harsh environment. a pops with ci: total number of populations in each treatment showing evidence of ci; pops with no ci: total number of populations not showing evidence of ci. b number of endpoint populations with or without evidence of ci pooled across both ancestors. c ho: the frequency of populations showing ci is not higher in benign conditions. another pattern that would support the occurrence of ci under our experimental conditions is the observation of mutations in a midpoint population that were not seen in the endpoint population of the same line. we therefore examined variation at the midpoint of some populations where the endpoint strongly indicated ci and in line of the harsh treatment where variation was present in the endpoint population. ci was evident in three of the four midpoint populations from control and benign treatments, and all three of these had mutations that were not seen among the endpoint isolates. in two of these cases, the lost genotypes were at very high frequency. in the harsh environment, there was no evidence of ci in the midpoint population, and although there was one mutation that was not observed in the endpoint population, that mutation was a synonymous mutation with no independent evidence of being adaptive. to examine whether mutations that were lost between the midpoint and endpoint populations could have been out-competed, we measured the selection coefficients of a mutation that showed evidence of being lost due to ci, and of the mutation that may have out-competed it. in the midpoint population of control line only c788t was detected, while in the midpoint population of benign line there were four genotypes, including c788t. by the endpoints, c788t was present in all genotypes of control line but was not detected in benign line and a1611g was detected in all genotypes in the endpoint populations of both treatments. apparently, a1611g did not interfere with the fixation of c788t in the control line, but did interfere in the benign line. if ci were responsible for this outcome, then in the benign line, both mutations should be adaptive and a1611g should have a significant advantage over c788t. indeed, both mutations had a significant selective advantage over their ancestor, and the a1611g mutation had a significant selective advantage over the c788t mutation, supporting that c788t was lost due to ci by a1611g in the benign line. frequency of multiple mutations evidence supporting the multiple-mutations model versus the one-by-one model of clonal interference would be the occurrence of more than one change from the fixed genotype. we found four populations where a single isolate differed from the fixed genotype by more than one change: benign midpoint, isolate 1702; control endpoint isolate 1216; control midpoint, isolate 581; and harsh midpoint, isolate in two of these cases, the changes were synonymous in gene f and there is no independent evidence to suggest that the additional change was adaptive. in the two cases where there is evidence that the additional changes were adaptive, they were seen only in the midpoint populations. additionally, since the control lines behaved like the benign treatment, they were treated as such in the analyses. discussion theoretical models and empirical data show that increased population sizes and mutation rates can slow the rate of adaptation and decrease the probability that beneficial mutations will reach fixation. phenotypic data from the empirical studies show that rates of adaptation do not increase linearly with population size, but it is not clear whether this pattern is actually due to competing mutations or simply because each mutation will have a longer transition time in larger populations. it is also not clear how many genotypes might be competing at any given time, which is an important parameter for studying properties of adaptation and predicting adaptive walks. we evolved nine lines of phage and sequenced genomes of multiple individuals from populations in these lines to examine the underlying genetics of ci. we found that mutations with relatively large selection coefficients can be out-competed and that there were as many as three competing genotypes in many of the populations where resource was abundant, even though the regular strong bottlenecks inherent to our flask-to-flask transfer protocol strongly decreased our ability to observe most mutations. this suggested both that beneficial mutations and ci were remarkably prevalent. our sequence data did not reveal any ci in the harsh treatment, which emphasized that the frequency of ci can vary among similarly large populations due to environmental conditions. adaptation under different levels of resource the resource treatments were intended to alter parameters of the distribution of beneficial mutations in order to evaluate environmental effects on levels of ci in similarly large evolving populations. while we found evidence that there were differences between lines in the available beneficial mutations, these differences were not as we predicted. we expected that the ancestor would be further from its fitness optimum under harsh resource conditions, which would result in a larger mean selection coefficient for beneficial mutations. however, there was insufficient evidence that the mean selection coefficient of beneficial mutations was higher in the harsh environment for anc, and it even appeared to be lower in the harsh treatment for anc apparently, the c788t mutation in anc increased fitness in the harsh environment, such that selection was weaker on anc than on anc in the harsh environment. also, although we expected that the variable resource levels would select for improvement on the same trait, the sequence data suggested that the stronger selection was on different targets between the two resource treatments, and possibly more genetic solutions were available in the benign environment. thus, it appeared that the distribution of beneficial mutations differed in overall frequency of adaptive mutations between the resource treatments, such that more mutations were available in the benign environment. factors affecting the frequency of ci when considering the overall frequency of ci for both ancestors, we observed ci in most populations in the benign environments but not in any populations in the harsh environments. while this result was consistent with our original prediction, it was not attributed to a larger mean selection coefficient in the harsh environment as we had predicted. for example, selection was possibly weaker for anc in the harsh environment, yet no ci was observed in any of those populations. additionally, although we found that the selection coefficient of the mutation that occurred in the harsh environment of line was very high, it was the only mutation observed, and thus chance cannot be out-ruled as an explanation for the difference in selection strength between the treatments in line however, there was evidence of an overall trend of fewer adaptive mutations in the harsh condition, suggesting that the lower levels of ci in the harsh environment may have been due to the lower number of available adaptive solutions, which would increase the waiting time for beneficial mutations and decrease their co-occurrence. interestingly, even though the frequency of ci appeared to differ among the two resource treatments, their overall increase in fitness was similar, as were the population sizes among experimental lines. this result emphasizes two important points: 1) that ci does not necessarily decrease rates of population adaptation, and 2) that levels of ci are not necessarily predictable from population size. regarding adaptation rate, it is indeed possible that ci acts to increase population adaptation by maintaining higher levels of genetic variation. for example, if a very high fitness peak can only be achieved through a specific ordering of substitutions, delayed fixation time from ci could increase the probability that primary mutations in the order are maintained long enough for the subsequent necessary mutations to arise. in the absence of ci, fixation of alternate mutations may be too rapid, occurring before the critical primary mutations, which could result in the population being stranded on a lower fitness peak than what could have been realized through a different adaptive walk. in this manner, the decreased rates of adaptation that are predicted to occur from competing mutations, could be outweighed by larger final fitness gains due to the more thorough exploration of sequence space that occurs in populations undergoing ci. this could partly explain our observation of similar levels of adaptation, but different levels of ci, among the two resource treatments. there appeared to be fewer adaptive solutions available in the harsh treatments since fewer different mutations were observed and most occurred in a single gene. all substitutions in the endpoint populations of the harsh treatment were in gph, a capsid protein involved in guiding phage dna into host cells, while those detected in the control and benign environments occurred in several different genes, including the gene encoding the major capsid protein, as well as in genomic regions that are putatively involved in regulation of lysis timing. gpf is involved in numerous processes throughout the phage life cycle including attachment to hosts, catalysis of eclipse, and assembly. the apparently larger number of genetic solutions in the benign environment relative to in the harsh environment suggested that: 1) selection acted on a broader range of phenotypes in the benign treatments, and/or 2) there were more genetic solutions to the types of phenotypes under selection in the benign treatments. in the harsh environment where ca+ was very limiting, there should be strong specific selection on capsid proteins to either be able to inject dna with fewer ca+ molecules or to become more efficient at acquiring ca+ since genome injection is crucial to initiating infection. in other words, in the harsh environment, improvement of this trait must come before other potential renovations to growth rate. in contrast, in the benign environments where growth rates were already high at the outset, modification of any phenotype affecting growth rate could be effective at increasing fitness. for example, one strong source of selection inherent to all treatments is bottleneck survival. faster lysis is a conceivable adaptive response under these conditions: phage genotypes with the highest number of free phage immediately before transfer will have the highest probability of being transferred to the next flask, and these genotypes will go on to initiate the most individual infections in the subsequent flask. thus, in the benign conditions adaptation could occur through any modification that produced a high number of phage by the end of each flask growth period whereas in the harsh environment selection specifically targeted injection rates before any other process in the phage life cycle, suggesting that the higher levels of ci in the benign environment could have been due to a larger number of adaptive targets under these conditions. experimental design issues a. bottleneck effects since phage had higher growth rates in the benign environment, bottlenecks to a smaller population size were required to maintain population sizes similar to those in the harsh treatments. thus, beneficial mutations that arose in the benign condition had a higher probability of extinction due to bottlenecks, suggesting that our design was slightly biased towards finding mutations in the harsh condition. interestingly, we found fewer adaptive mutations and lower levels of ci in the harsh treatment where the extinction probabilities of beneficial mutations due to bottlenecks were lower. this demonstrates that the difference in ci between the two treatments may have been even more dramatic if bottleneck sizes were equal. thus, our approach to testing the hypothesis that ci would be stronger in benign environments was conservative. b. growth rate differences the combined effect of differences in growth rates for phage in the two treatments and frequent, regular, strong bottlenecking is an experimental design factor that could have increased the extinction probability of mutations in the harsh environment. previous work has shown that phage with lower growth rates have increased burst times, suggesting that replication is slower. if the higher growth rates of phage in the benign environment translate to faster replication relative to that in the harsh environment, then mutations would arise earlier in flask growth. from the work of wahl and gerrish, we found that mutations arising in the last five minutes of flask growth have a greater than 99% chance of being lost due to sampling, regardless of their selection coefficients, while those arising in the first five minutes are likely to survive the bottleneck if their selection coefficients are high. thus, the lower growth rates of phage in the harsh environment could have increased the extinction probability of beneficial mutations if most were produced too late during flask growth to survive the bottleneck. if present, this effect could explain the lower levels of ci observed in the harsh condition. however, it is also possible that the less severe bottleneck in the harsh condition, which decreased extinction probabilities, acted to negate potential effects from low growth rates. implications of ci for predicting adaptation in four of the five lines where we sequenced isolates from the midpoint populations, we found genotypes that were no longer present in the endpoint population, which we refer to as'lost.' furthermore, from measuring the selection coefficients of one of these lost individuals we found that mutations with relatively large selection coefficients can be lost. these data demonstrated two important points that are relevant to understanding and predicting adaptation. first, the probability that a beneficial mutation is lost due to a deterministic evolutionary phenomenon such as ci depends partly on stochasticity from the time that other mutations of similar effect arise. in populations that are subject to strong sampling effects, it is possible for a small-effect mutation to reach fixation in the presence of a larger one, as long as the smaller-effect mutation is near fixation before the larger-effect mutation arises, such that the larger-effect mutation is more strongly affected by the stochastic sampling. second, the fact that relatively large-effect beneficial mutations were lost in competition supports the finding that beneficial mutations may be much more frequent than believed. many beneficial mutations are often not seen because they can only be sampled in environments where selection coefficients are similar among mutations, or where there are no mutations with very large selection coefficients. this emphasizes the importance in measuring adaptation in numerous environments in order to develop predictive models and understand evolutionary processes. the co-occurrence of high frequency beneficial mutations has been previously observed during experimental evolution of microvirid phages in batch culture. we found evidence for ci in five out of nine populations even though all mutations had high extinction probabilities, suggesting that ci can be remarkably common in large bacteriophage populations. while in most cases competing genotypes differed from the fixed genotype by only by a single adaptive mutation, which supports predictions of the one-by-one model of clonal interference, we found limited support for the multiple mutation model. there were four cases where a single isolate differed from the fixed genotype by more than one change, although in two of these cases there is no evidence that the additional change is adaptive. in no case did these genotypes carrying multiple mutations dominate in the endpoint populations. taken together, our data and the recent study in yeast emphasize that concurrent mutations are an important determinant of adaptive trajectories in large microbial populations. nevertheless, sequence-based models of adaptation that assume selection is strong relative to mutation, such that one beneficial mutation is fixed before the next occurs, have been considered acceptable for describing properties of adaptive walks in asexual organisms. our experimental data suggest that sswm conditions may often be unrealistic in large microbial populations, and thus that models of adaptation relying on this assumption may make misleading predictions regarding the properties of adaptive walks in large microbial populations. models of adaptation that include concurrent mutations and multiple concurrent mutations on individual genotypes may provide more realistic frameworks for investigating the properties of adaptation in large microbial populations.  we verified with sequence data that ci occurs in large virus populations and that several genotypes may co-occur at very high frequencies. we intended to test the prediction that weak selection in benign environments would result in more ci on average. while our data showed that ci was more frequent in benign treatments, this difference was not attributed to weaker selection. it appeared that the higher frequency of ci in the control and benign environments was partly due to a higher number of potential adaptive solutions in those environments, relative to that in the harsh treatment, which would have decreased the waiting time for beneficial mutations. at the very least, our results showed that environmental conditions can alter the frequency of ci among similarly large populations, which demonstrates that effects of ci are not necessarily predictable from population size. another important insight from our genetic data was that selection appeared to target different phenotypes in treatments where the same continuous resource variable was decreased to different degrees. finally, we observed that the fate of mutations with relatively large selection coefficients depends on the time that other large-effect mutations arise, suggesting that these stochastic population-level processes are significant properties of adaptive landscapes, as predicted by theoretical models. 
theoretical models predict that clonal interference can slow adaptation, alter the distribution of fixed beneficial mutations, and affect disease progression by impacting within-host evolution of pathogens. the sequence data revealed that observed mutations were limited to a single portion of one gene in the harsh treatment, but spanned a different and larger region of the genome under the benign treatments, suggesting that there were more adaptive solutions to the benign treatment. additionally, the genomic sequence data showed that selection can target different types and numbers of phenotypes in environments that differ by only a single continuous variable. in large asexual populations where beneficial mutations may co-occur and recombination is absent, the fate of beneficial mutations can be significantly affected by competition. genomic sequences of multiple individuals from evolved populations were used to detect competing beneficial mutations. beneficial mutations with relatively large selection coefficients can be excluded by ci. ci may commonly determine the fate of beneficial mutations in large microbial populations, but its occurrence depends on selective conditions. while phenotypic data support that ci is a significant determinant of adaptive outcome, genetic data are needed to verify the patterns and to inform evolutionary models. ci was more frequent in benign selective conditions possibly due to a greater number of adaptive targets under this treatment. we adapted replicate populations of the bacteriophage φx under two levels of cacl to create benign and harsh environments. there were several competing genotypes in most of the populations where cacl was abundant, but no evidence of ci where cacl was scarce, even though rates of adaptation and population sizes among the treatments were similar.
 in large asexual populations where beneficial mutations may co-occur and recombination is absent, the fate of beneficial mutations can be significantly affected by competition. theoretical models predict that clonal interference can slow adaptation, alter the distribution of fixed beneficial mutations, and affect disease progression by impacting within-host evolution of pathogens. while phenotypic data support that ci is a significant determinant of adaptive outcome, genetic data are needed to verify the patterns and to inform evolutionary models. we adapted replicate populations of the bacteriophage φx under two levels of cacl to create benign and harsh environments. genomic sequences of multiple individuals from evolved populations were used to detect competing beneficial mutations.  there were several competing genotypes in most of the populations where cacl was abundant, but no evidence of ci where cacl was scarce, even though rates of adaptation and population sizes among the treatments were similar. the sequence data revealed that observed mutations were limited to a single portion of one gene in the harsh treatment, but spanned a different and larger region of the genome under the benign treatments, suggesting that there were more adaptive solutions to the benign treatment.  beneficial mutations with relatively large selection coefficients can be excluded by ci. ci may commonly determine the fate of beneficial mutations in large microbial populations, but its occurrence depends on selective conditions. ci was more frequent in benign selective conditions possibly due to a greater number of adaptive targets under this treatment. additionally, the genomic sequence data showed that selection can target different types and numbers of phenotypes in environments that differ by only a single continuous variable. 