- Research article
- Open Access
A common, non-optimal phenotypic endpoint in experimental adaptations of bacteriophage lysis time
BMC Evolutionary Biology volume 12, Article number: 37 (2012)
Optimality models of evolution, which ignore genetic details and focus on natural selection, are widely used but sometimes criticized as oversimplifications. Their utility for quantitatively predicting phenotypic evolution can be tested experimentally. One such model predicts optimal bacteriophage lysis interval, how long a virus should produce progeny before lysing its host bacterium to release them. The genetic basis of this life history trait is well studied in many easily propagated phages, making it possible to test the model across a variety of environments and taxa.
We adapted two related small single-stranded DNA phages, ΦX174 and ST-1, to various conditions. The model predicted the evolution of the lysis interval in response to host density and other environmental factors. In all cases the initial phages lysed later than predicted. The ΦX174 lysis interval did not evolve detectably when the phage was adapted to normal hosts, indicating complete failure of optimality predictions. ΦX174 grown on slyD-defective hosts which initially entirely prevented lysis readily recovered to a lysis interval similar to that attained on normal hosts. Finally, the lysis interval still evolved to the same endpoint when the environment was altered to delay optimal lysis interval. ST-1 lysis interval evolved to be ~2 min shorter, qualitatively in accord with predictions. However, there were no changes in the single known lysis gene. Part of ST-1's total lysis time evolution consisted of an earlier start to progeny production, an unpredicted phenotypic response outside the boundaries of the optimality model.
The consistent failure of the optimality model suggests that constraint and genetic details affect quantitative and even qualitative success of optimality predictions. Several features of ST-1 adaptation show that lysis time is best understood as an output of multiple traits, rather than in isolation.
Models that describe evolutionary processes must strike a balance between simplicity and detail. At one end of the spectrum are optimality models, which predict phenotypic evolution based on natural selection in the context of ecology, often neglecting genetics except in the form of simple 'trade-offs' between phenotypic alternatives. The appeal of such models is in their apparent generality, enabling prediction across a wide variety of organisms and environments with minimal input about genetic parameters. The weakness of optimality models is argued to lie in their neglect of genetic details, the very matrix on which the evolutionary process rests . If these details lead to constraint or unexpected pleiotropic interactions between phenotypes, optimality model predictions may mislead.
There is a long history of disagreement about the importance of genetic details in understanding phenotypic evolution . There seems to be little resolution. Optimality models are still widely applied with apparent success for some biological systems [3–10], but studies that account for genetic details seem to gain insight from their incorporation [11–14]. As yet, however, there are few studies that both enlist information about genetic details and conduct quantitative tests of optimality model predictions. Such direct comparisons of the two approaches are essential in resolving the merits of genetic details in understanding evolution.
Here we use experimental evolution of bacteriophage to test an optimality model for a life history trait, lysis time. Environmental conditions are easily controlled when growing bacteriophages, making it simple to parameterize and test quantitative predictions about optimality in reasonably well-understood genetic systems. The lysis time model has been tested on three phages, λ, T7, and T4 [15–17], but the present work extends these tests to a pair of phages with a very different lysis mechanism.
Lysis time optimality model
Most bacteriophages lyse their hosts to enable progeny to escape into the environment (Figure 1). Lysis of the infected cell releases all phage progeny at once and so is the phage equivalent of the age of reproduction. Furthermore, this life history trait has major effects on phage fitness. The infected cell typically produces phages indefinitely until lysis, so a long lysis interval increases progeny but extends generation time, and vice versa.
A simple optimality model predicts which lysis interval will maximize fitness simply as a function of the intrinsic growth rate of the phage population [22, 23]. The optimality model assumes linear accumulation of phages within the cell and parallels a theory of optimal foraging in patchy environments originally developed by Charnov . It leads to an approximate lysis time optimum of
[16, 23]. The optimal lysis time () minus the eclipse time (E , when functional phage progeny begin to accumulate inside the cell) should evolve to equal the inverse of r (the intrinsic rate of increase, a measure of growth rate). is defined here as the lysis interval. All viable phage progeny are produced inside the cell during this interval, at an approximately constant rate (Figure 1). Since this equation applies at the evolutionary equilibrium, the optimal lysis interval changes when the environment changes or the phage evolves in response to a change in the environment.
The simple model presented here ignores the effects of non-genetic lysis variation. Lysis asynchrony (in which genetically identical phages lyse at different times due to phenotypic noise) is important in ΦX174 -like phages and affects the optimal lysis interval. A more complete model that accounts for this parameter was used for optimality predictions throughout (see equations 2 and 3 in Methods for details).
This model can be explored with experimental evolution [15–17, 25]. As the genetic basis of lysis is well studied and apparently simple for many phages, lysis time is an ideal system for exploring how optimality interacts with genetic details in a life history trait.
Lysis in isometric phages
Only two assumptions about lysis mechanisms are required by the optimality model; both are met in ΦX174 bacteriophage. First, average lysis time is evolvable, i.e. mutations exist that can change it [18, 26]. Second, phage progeny accumulate linearly within the infected cell until the end of the lysis interval [18, 27]. By similarity of genome organizations, we assume ST-1 also meets these assumptions.
ST-1 and ΦX174 are isometric, icosahedral phages of the family microviridae with small single-stranded DNA genomes. ΦX174 has 11 genes, while ST-1 encodes homologs to those 11 plus an extra region with putative genes, like certain other isometric phages . In both phages, only one gene (gene E) is directly implicated in lysis [29, 30], and it overlaps, out of frame, the highly conserved gene D  which encodes the external scaffolding protein .
The gene E protein, GpE, causes lysis by interfering with cell wall synthesis during cell division [32, 33], a simpler mechanism than that seen in larger phages . All mutations in gene E known to affect lysis do so by altering expression or by eliminating lytic activity altogether [27, 35]. GpE noncompetitively inhibits MraY , and without this enzyme the cell expands and eventually lyses rather than dividing .
ΦX174 requires the host factor slyD for lysis [29, 38, 39]. Mutations in gene E that increase gpE expression can restore lysis [35, 40]. In addition, when these mutant alleles were expressed from a λ phage construct they hastened mean lysis time in the presence of slyD , suggesting a potential for quantitative variation in ΦX174 lysis time.
We adapted these two phages across serial flask transfers for hundreds of generations to identify the mechanisms and generality of evolution to an optimal lysis time. The universal outcome for ΦX174 adaptations was a lysis interval indistinguishable from that of wild-type, despite the predictions of the optimality model that a shorter lysis interval should evolve. In ST-1, the lysis interval shortened somewhat but did not reach the predicted optimum.
We adapted phages under a variety of environmental conditions to determine if they evolved to lyse at the predicted optimal time. In every replicate, the lysis interval remained longer than the optimality model predicted. Changes in the single known lysis protein occurred only when phages were adapted on hosts that eliminated wild-type phage lysis entirely.
Lysis time is the sum of eclipse and lysis interval (before and after phage progeny start to accumulate in the cell, respectively) and the model assumes that eclipse time is fixed and that the optimal lysis time is achieved by changes in the lysis interval. We thus address primarily the lysis interval.
ΦX174 phenotypic adaptation to E. coli C: Failure to evolve shorter lysis time
We adapted ΦX174 to E. coli C twice under similar conditions, in which the predicted optimal lysis interval was shorter than the initial lysis interval. Thus we expected the lysis interval to shorten over adaptation. One adaptation started with wild-type phage, while the other started with a recombinant population to increase initial variation. Through this latter mechanism, ΦX174+ was provided with substitutions from two lines preadapted to conditions differing slightly from those used here. Higher genetic variance in the population should facilitate adaptation. Preliminary analysis showed that total lysis time was broadly similar between the two adaptations, so an isolate from the recombinant population (ΦX174C) was chosen for further characterization.
Following 34 hours of adaptation, ΦX174C fitness increased from 20.8 to 23.4 db/hr, an increase of 2.6 db/hr (Table 1, p < 0.03 by 1-tailed t-test). No other phenotypes measured changed significantly, although eclipse time was possibly shorter by ~1 min (p < 0.07 by 2-tailed t-test, Table 1). The lysis interval (11.4 min) did not evolve to the predicted optimum of 5.5 min (calculated using lysis time variance as in equation 3, p < 0.0002 by the parametric model).
The failure of the mean lysis interval to adapt closer to the optimum over ~120 generations means that even qualitative predictions of evolution failed for ΦX174C. Within the framework of our model, selection for earlier lysis did exist; from projections based on measurements of life history traits, the evolved phage was 4.3 db/hr away from the fitness it might have reached at optimality. This potential fitness advantage is large compared to the fitness difference of 2.6 db/hr between the initial and evolved phages, although measurements of small differences in absolute fitness may not be reliable .
ΦX174 phenotypic adaptation to hosts blocking lysis: Failure to evolve optimally despite evolution
The slyD gene of E. coli is needed for lysis by wild-type ΦX174. However, in agreement with previous work , ΦX174 mutants that formed plaques on slyD-defective hosts were easily isolated. By definition, any mutant that lyses this host at all lyses earlier than the wild-type phage. We were interested in the lysis time of these mutant phages and of their descendants that were further adapted to the inhibitory host. A mix of two plaques was used to initiate an extended adaptation. After 35.5 hours of adaptation, the resultant phage (ΦX174ΔslyD) had a lysis interval of 11.2 min, similar to that of ΦX174C (11.4 min, Table 1). Not only are the lysis intervals indistinguishable between the evolved phages on their respective hosts, but the fitnesses are similar: 23.4 db/hr for ΦX174C versus 21.2 for ΦX174ΔslyD (p < 0.006). Again, the virus failed to reach its predicted optimum (5.5 min) by almost 6 min (p < 0.0052 by parametric test).
The mechanism of adaptation to slyD-deficient hosts is increased expression of the phage lysis gene E, which in a λ phage construct results in earlier lysis even in the presence of slyD . By this criterion, a phage adapted to grow on the slyD host should lyse earlier when grown on normal hosts. Matching this expectation, ΦX174ΔslyD grown on normal hosts had a lysis interval of 8.9 min, shorter than the interval on slyD hosts (11.2 min, p < 0.025 by 2-tailed t-test).
On normal hosts, ΦX174ΔslyD lysis interval was also somewhat shorter than that of ΦX174C, although this difference was only marginally significant (Table 1, 8.9 min rather than 11.4 min, p < 0.1 by 2-tailed t-test). If this difference is real it would suggest that ΦX174C lyses later than predicted due to pleiotropic costs, rather than absolute constraint.
Selecting for later lysis
If our model was wrong because it calculated an inaccurate optimum rather than an unattainable optimum, then a change in the growth conditions to select for later lysis should result in slower lysis after adaptation, because we know the virus can attain that phenotypic state. This can be tested by growing the phage at a sufficiently lower host density to favor longer lysis interval. The same starting population of phages used for ΦX174ΔslyD was grown on hosts at 2 × 106 cells/ml, which greatly lowered fitness. Even after 25.6 hr of adaptation, the final phage (ΦX174ΔslyDlow) had a much lower fitness (absolute growth rate) under these conditions than ΦX174ΔslyD had when grown at 1 × 108 cells/ml (Table 1, 10.6 rather than 20.4 db/hr). However, lysis interval evolved to a similar endpoint (10.0 min ± 0.11 min standard error for ΦX174ΔslyDlow as compared to 11.2 ± 0.54 min for ΦX174ΔslyD (p < 0.13 by 2-tailed t-test). Although the optimal lysis time for ΦX174ΔslyDlow (8.7 min) was considerably later than that of ΦX174ΔslyD (5.5 min), ΦX174ΔslyDlow still failed to attain its optimum (p < 0.0046 by parametric test).
ST-1 phenotypic adaptation: Failure to evolve optimally despite evolution
ST-1 initially had an eclipse time of 7.7 min and a lysis interval of 10.4 min. The phage adapted to laboratory conditions (ST-1K-12) had a significantly earlier eclipse time (p < 0.01) and shorter lysis intervals (p < 0.005), each by ~2 minutes, resulting in lysis approximately 4 min earlier than ST-1+ (Table 1, Figure 2; see Additional file 1 for representative curves). Attachment rate was high throughout and did not change (Table 1), suggesting this adaptation was not driven by attachment evolution. Although the burst size estimate was somewhat smaller after adaptation, this difference was not significant (p < 0.17 by 2-tailed t-test).
The lysis interval of the evolved phage ST-1K-12 (8.0 min) did not reach the predicted optimum of 4.7 min (p < 0.004 by parametric test). This deviation has a predicted fitness cost of 4.4 db/hr. Thus, ST-1 lysis interval appears to have evolved qualitatively as expected, but did not reach the optimum despite selection to do so.
Molecular evolution in both phages
Neither ΦX174C nor ST-1K-12 carried any coding or apparent regulatory changes in gene E, which generates the only protein known to play a role in ΦX174 lysis. Thus there is no obvious genetic cause for the decrease in lysis interval observed in ST-1. ΦX174C contained seven substitutions relative to ΦX174+ (Table 2), two of which were not present in the preadapted founding population. ST-1K-12 acquired four substitutions, two during the initial adaptation to purely permissive hosts and two more in phase II when nonpermissive hosts were also present (Table 3).
While many of the genetic changes of ΦX174C were present in the founding lines, two mutations were new. This shows that there was time for new substitutions to fix, and thus the potential for evolution closer to the optimum that did not rely on changes from the preadaptations. All 7 of the ΦX174C changes together only increased the fitness of the wild-type phage by 2.6 db/hr, suggesting that any single change moving lysis time near the optimum was not available.
In the two adaptations to slyD-deficient hosts, however, every change observed was in the lysis gene E (and overlapping gene D, Table 3). Two of these changes have been previously observed to restore lytic capacity on these hosts . As the adaptations were started from the same population, no conclusion can be drawn from the similarity of genetic changes between the adaptations. However, they do demonstrate that gene E evolution changes lysis time under some conditions. If ΦX174ΔslyD does lyse faster on regular hosts than ΦX174C does, as seems possible (see previous section) the lack of these gene E changes in the other lines suggests that there is some cost to the changes that prevent them from spreading under normal conditions despite optimality predictions.
This study provides an experimental test of a quantitative optimality model in a phage system. Although optimality models are widely used to explain nature, there have been few attempts to actually evolve organisms to an optimum. Is it possible to predict phenotypic evolution using a simple model based largely on natural selection?
Optimality models may fail for a variety of reasons. An evolutionary constraint might prevent any adaptation at all from the starting, non-optimal state. For example, in a study of optimal foraging evolution, comparisons between ΦX174 and the related phage G4 demonstrated that selection would favor indiscriminate foraging under certain conditions, yet an absolute genetic constraint prevented adaptation of G4 to infect a new host . However, another optimal foraging study using T7 showed that adaptation was possible for that phage . Another form of constraint would be genes with pleiotropic effects on traits not included in the model. In Φ6, expansions in host range were often associated with decreased attachment to the original host  while in T7 a similar mechanism caused a failure of some optimal foraging predictions . There is thus precedence for constraint on adaptation of simple phenotypes in phage genomes, but the extent of constraint is variable.
Phage lysis interval, and thus lysis time, is likely to be under strong and fluctuating selection in the wild, thus providing an ecologically important model system for the evolution of a life history trait. The lysis time optimality model is based on a trade-off between progeny number and generation time. It predicts how long a phage should wait after eclipse to lyse its host based on the growth rate of the phage, assuming that phage progeny accumulate linearly within the host after eclipse. The model has been previously tested in three phages: λ, T4 and T7. In all of those phages, lysis is thought to require a 2-component lysis mechanism; a lysin that destroys the cell wall and a holin that destroys the cell membrane allowing the lysin to access the cell wall. The holin triggers lysis with great precision in at least some phages .
The optimality model was tested in phage λ and T4 by competing lysis mutants with wild-type phages. In both cases, as predicted by the model, different environments selected for different mutants. [15, 16, 25]. In T7, quantitative optimality predictions were tested by allowing lysis to evolve without previously introduced genetic changes known to affect lysis. Adaptation had mixed success; selection for late lysis did not generate the optimal phenotype. Model failure was attributed to a violation of the model assumption that phage particles increase linearly inside the cell if the normal lysis mechanism is suppressed .
This prior work motivated the work here. How would the model fare in a one-component system such as ΦX174 and ST-1? ΦX174 satisfies the model assumption of linear accumulation . In addition, mutations in the lysis gene (E) have been found that hasten lysis time under some conditions . There was thus no a priori reason that the mean lysis interval for the population could not evolve to the optimum.
Under selection for a shorter lysis interval (according to the optimality model) ST-1 evolved earlier lysis by 4 min (although 2 min of this was caused by earlier eclipse), while ΦX174 did not; no changes in the ΦX174 lysis gene evolved in any adaptation to normal hosts. In none of the adaptations on either phage did lysis interval reach the predicted optimum. Quantitatively optimal evolution is not expected in a single step, because a large change is not likely to carry the phenotype directly to the optimum [47, 48]. However, it is surprising how far from the optimum the phages in this study remained after adaptations involving a number of genetic changes.
While precise control of the phenotype appears difficult, it is clear that rapid evolution is possible in a different part of the phenotype space. When ΦX174+ is grown on hosts lacking the cofactor slyD, it does not lyse detectably, yet the adapted phage attains an lysis interval on the slyD-defective hosts indistinguishable from that of a virus adapted to normal hosts. This is similar to the pattern seen in T7, which can recover from complete genetic ablation of lysis but cannot reach the predicted optimum [17, 49]. Interestingly, despite the ability of ΦX174 to recover lysis on slyD-defective hosts, genetic details are relevant here as well; the coding mutations seen in our adaptations were the same as those previously found in another study , demonstrating a limited number of paths to maximal recovery and thus the possibility for constraint.
It could be that the phages would have reached the optimum if allowed to adapt for longer. However, repeated evolution of ΦX174 to a common phenotypic endpoint and the observation that phenotypic evolution was rapid over other areas of phenotype space suggest that this was not the primary limit on adaptation.
Instead, constraint on lysis interval, either absolute or due to pleiotropic effects on other phenotypes, seems the most likely explanation for failure to reach the predicted optimum. This may be caused by the overlap of gene E with gene D, which is highly conserved , or by the relatively crude ΦX174 lysis mechanism . If evolution were possible but the quantitative optimality prediction itself was inaccurate, introducing major changes that would alter the optimum would lead to a different results, yet this had no effect.
A large part of the ST-1 lysis time adaptation was driven by changes in eclipse time. Phenotypes that are outside the scope of a model may frequently play a role in evolution. When optimality models are tested in nature, researchers often form conclusions after considering only the life history trait thought to be under selection. However, most life history traits, like lysis time, result from the interplay of a number of factors, which may mask the nature of adaptation.
This complexity is also evident on the genetic level. For example, ST-1 evolved to lyse closer to the optimum (shortening both its post-infection lysis time and its lysis interval) without changes in gene E, the only gene implicated in the widely accepted 1-component model for lysis in ΦX174-like phages . Other genes (A, F and H) that did evolve in ST-1 may play some heretofore unsuspected role in mediating lysis. However, it is far more likely that, as previously observed in T7 phage, genes with no direct role in cell lysis affect lysis timing by altering timing of lysis gene expression or by indirect effects on other stages of the life cycle rather than lysis mechanisms per se [17, 20, 49]. Similar interactions may make optimality predictions in other organisms more difficult.
Optimality models can be used in two ways: as a means of producing accurate phenotypic predictions or as null hypotheses. As null hypotheses, optimality models help us detect and explore violations of model assumptions. Here, for example, the quantitative optimality model was strictly required in order to identify constraint. However, it may be that even in well studied systems, genetic details make it difficult to use simple models to quantifiably predict the evolution of life history traits across taxa. The constraint we observed, the difference in model success between replicates from different species, and the evolution of traits outside the scope of the model are all consistent with this outlook.
Cell and phage lines
ΦX174 was grown on E. coli C cells. The ΦX174 adaptation was initiated with a wild-type ΦX174 allowed to recombine with two preadapted ΦX174 phages (by cross-streaking all phages to allow coinfection) in order to introduce genetic variation. One isolate, ΦX174 Af , had been adapted to E. coli C under the same conditions used here except that phage generations were restarted approximately every 35 minutes through addition of chloroform. The other, ΦX174pif- , was adapted without frequent chloroform treatment but on a different host, IJ1862 (E. coli C/F'128 pifA15 lacZΔM15Tn10). This recombinant population was used deliberately to provide potentially adaptive substitutions to the wild-type phage and thus maximize the potential evolution. The source of all changes in the resultant evolved line was determined by sequencing. ΦX174ΔslyD and ΦX174ΔslyDlow were adapted to slyD-defective host CCX1 (E. coli C slyD zhd::Tn10).
ST-1 phage [51, 52] was adapted on E. coli K-12 strain BW25113 (rrnB3 ΔlacZ4787 hsdR514 Δ(araBAD)567 Δ(rhaBAD)568 rph-1) with an additional deletion of lacY, part of the Keio collection . Another Keio collection mutant (Δrep), identical except that it instead lacked the rep helicase protein, made up 95% of the cells in the second phase of ST-1 adaptation. The Δrep host was nonpermissive, as ST-1 infects it but does not produce viable phage offspring , and thus its presence lowered overall fitness.
Adaptation protocols were similar to that described in other work  with some modifications. Specifically, 10 ml of LB media (10 g NaCl, 10 g Bacto tryptone, 5 g Bacto yeast extract), supplemented with 5 mM CaCl2 to improve adsorption of both Φ174 and ST-1 to their hosts, was pre-warmed to 37°C in an orbital water bath (200 rpm). Cells from recently thawed 20% glycerol-LB stocks stored at -80°C were then added and allowed to grow for an hour such that the cell population was 1-2 × 108 cells/ml (1-2 × 106 cells/ml for ΦX174ΔslyDlow) and in exponential phase when phages were added.
Phages were incubated in this cell culture for only 20-30 min, minimizing changes in host density over time, which would confound the optimality model predictions. Following incubation, an aliquot of the infected culture - which included both free phage and infected cells - was transferred to a new flask in which cells were at the requisite density (having grown one hour). The amount transferred was such that roughly 105 phages total (no less than 104 and no more than 107) were added to the next flask. In addition, at the end of each flask, a sample of the completed flask was chloroform-treated and stored, preserving free phage and killing remaining hosts. The chloroformed sample from the final flask of the day was then used as the starting population for the first flask of the next day.
ΦX174 was adapted for 34 hours (excluding the pre-adaptation time of the starting populations) to yield ΦX174C, with a second adaptation of ΦX174+ over 35 hours to get ΦX174C#2. Both adaptations to slyD-defective hosts started from a 1:1 mixture of two mutant plaques. ΦX174ΔslyD was adapted for 35.5 hr, while ΦX174ΔslyDlow was adapted for 26.6 hr. ST-1 was adapted for 26.3 hours on purely permissive hosts in phase I of adaptation, then on a mix of 95% nonpermissive and 5% permissive hosts for 41.5 hours in phase II of adaptation (resulting in ST-1K-12). The latter conditions did not change the direction of selection on lysis time, but merely facilitated longer, more convenient passages. The nonpermissive hosts reduced phage growth rate and thus prolonged the time before phage density exceeded cell density, so transfers could be performed less often. This lower fitness led to a slightly later optimal lysis time (which ST-1K-12 still did not attain). The mutagen N-methyl-N'-nitro-N-nitrosoguanidine was added at a concentration of 0.5 ug/ul to one flask at 17.5 hr of the ST-1 adaptation to promote adaptation by generating new mutations.
Viral fitness was measured under conditions essentially identical to those used for adaptation. The rate of increase was determined by growth at low phage/cell ratios (not exceeding 0.1 by the end of transfer) across five consecutive transfers, using phage titers measured from the end of the second flask to the end of the fifth [see 17 for details]. Fitness is reported as doublings per hour, calculated as [log 2(Nt/N0)]/t, where Nt is the total number of phage in the flask at time t hours, corrected for dilutions over multiple transfers. The intrinsic rate of increase per minute, r, is simply ln(2)/60 times the fitness in db/hr, but db/hr is easier to comprehend.
Eclipse time, lysis interval and burst size
To measure eclipse, E, 5 × 107 phages were added to 10 ml of LB with 5 mM CaCl2 and cells previously grown to 108/ml in an orbital water bath, then incubated for three minutes. Cultures were diluted 103 fold in a new (preheated) flask to halt further adsorption. Samples were removed at a number of time points and diluted 10-fold in HFB  containing 2 mg/ml lysozyme saturated with chloroform, which lyses cells infected with ΦX174-like phages . Eclipse time was estimated by plating these samples, then fitting the resultant curve to an empirical least squares model .
A similar protocol was used to determine lysis interval and burst size of phage infections. Cultures were diluted 104-fold and 105-fold to new flasks after 3 min. The total number of free phage plus infected cells was estimated by direct plating of appropriate dilutions; plating after chloroform treatment but prior to 7.5 min (before progeny production) gave the number of free phage. This allowed calculation of infected cell density. Burst size was calculated by comparing this number with the final phage count after phage production had reached a plateau (in chloroform-treated samples, ~30 min after infection).
The data from direct plating at a number of time points allowed determination of lysis time and variance of lysis by an empirical least squares model . The asynchrony of lysis is reported as the standard deviation (σ) of the time infected cells took to lyse, and is larger when individual cells in the population lyse at very different times due to phenotypic noise.
Plaques were resuspended in LB within one day of plating and used for attachment assays within 3 days to ensure freshness. Between 105 - 106 phages were added to cells grown to 108 cells/ml at 37°C (as described above). After 5 min, part of the sample was spun down to remove attached phages. Spun and unspun cultures were plated to calculate free (N free ) and total (N total ) phage densities. Adsorption rate α was calculated from N free = N total × e-5Cα, where C is the cell density in cells/ml.
Phage genome sequences used dideoxy chain termination reactions with ABI Big Dye Terminator and were analyzed with an ABI3100 automated sequencer. Sequence files were checked and compared with DNA Star Lasergene Seqman software (version 5.05). All sequencing was from polymerase chain reaction products, except the wild-type ST-1, which was sequenced from genomic DNA. The entire phage genomes from the end of ST-1 phase I adaptation, ST-1K-12, and ΦX174C were sequenced from two isolates apiece, and one of the sequenced clones from each was selected randomly for phenotypic assays. ΦX174pif-, one of the initial phages for the ΦX174 adaptation, was sequenced over regions that differed between the ΦX174+ phage and ΦX174C in order to determine which changes arose during the final adaptation. All ΦX174 DNA nucleotide changes are reported using the reference ΦX174 sequence [GenBank: NC_001422].
The full genome of ST-1 was sequenced and is reported as part of this study [GenBank: GQ149088]; homologues to all ΦX174 genes were found. ST-1's genome is most similar to the ΦK and α3 genomes, which are ΦX174-like but members of a distinct clade . The genome includes a region encoding a number of putative genes that are homologous to some other ΦX174-like phages, such as α3 and WA13, but not ΦX174 itself .
The simple optimality model for lysis presented in the background can be extended to account for the effects of lysis time variance, which is important for ΦX174 and may play a role in other phages as well . If
(where is the lysis time at equilibrium and E the invariant eclipse time), and the lysis interval (L - E) is assumed to vary in a gamma distribution with variance σ2 between individual infections, the following equation
is obtained . This can be solved numerically for . This equation is used throughout this paper to predict optima.
Another equation allows calculation of expected phage fitness (r), which is useful when trying to determine how much higher the fitness would be if lysis interval changed.
when lysis interval has a gamma distribution with mean μ and variance σ2 . R is the constant rate at which phage accumulate within a cell between eclipse and lysis, c is the product of cell density and adsorption.
Estimating fitness benefit of optimal lysis
The expected fitness of phages was calculated based on empirical parameters of eclipse time, lysis time, lysis variance, intracellular growth rate, and adsorption by solving equation (4) numerically. This parameter-based fitness estimation, which usually differed from direct measurements of fitness, was used to approximate the fitness effect of setting lysis time equal to the optimum.
To test whether L - E - X = 0 (that is, whether the lysis time was at an optimum in an environment), we used empirical data to simulate a normal distribution for eclipse time, post-infection lysis time, lysis variance, and r (growth rate), parameterized with the means and variances calculated from the data. Numbers were randomly drawn from these distributions, with all numbers less than zero rejected as biologically implausible. Lysis variance and r from the distributions were used to solve for X by equation (3). This was repeated 5000 times, recording the proportion of times L - E - X had the opposite sign from the mean deviation from the optimum. The result multiplied by two is the 2-tailed probability that the means were identical. The significance of all other differences was determined by t-tests assuming unequal variances [58–63].
Lewontin RC: A natural selection. Nature. 1989, 339: 107-10.1038/339107a0.
Gould SJ, Lewontin RC: The spandrels of San Marco and the Panglossian paradigm: a critique of the adaptationist programme. Proc R Soc Lond B Biol Sci. 1979, 205: 581-598. 10.1098/rspb.1979.0086.
Pribil S, Searcy WA: Experimental confirmation of the polygyny threshold model for red-winged blackbirds. Proc Biol Sci. 2001, 268: 1643-1646. 10.1098/rspb.2001.1720.
Frankino WA, Zwaan BJ, Stern DL, Brakefield PM: Natural selection and developmental constraints in the evolution of allometries. Science. 2005, 307: 718-720. 10.1126/science.1105409.
Jensen KH, Little T, Skorping A, Ebert D: Empirical Support for Optimal Virulence in a Castrating Parasite. PLoS Biol. 2006, 4: e197-10.1371/journal.pbio.0040197.
Roff DA, Heibo E, Vollestad LA: The importance of growth and mortality costs in the evolution of the optimal life history. J Evol Biol. 2006, 19: 1920-1930. 10.1111/j.1420-9101.2006.01155.x.
Noblin X, Mahadevan L, Coomaraswamy IA, Weitz DA, Holbrook NM, Zwieniecki MA: Optimal vein density in artificial and real leaves. Proc Natl Acad Sci USA. 2008, 105: 9140-9144. 10.1073/pnas.0709194105.
Herre EA: Sex Ratio Adjustment in Fig Wasps. Science. 1985, 228: 896-898. 10.1126/science.228.4701.896.
Carvalho AB, Sampaio MC, Varandas FR, Klaczko LB: An experimental demonstration of Fisher's principle: evolution of sexual proportion by natural selection. Genetics. 1998, 148: 719-731.
West SA: Sex Allocation (Monographs in Population Biology). 2009, Princeton, NJ: Princeton University Press
Hanifin CT, Brodie ED: Phenotypic mismatches reveal escape from arms-race coevolution. PLoS Biol. 2008, 6: e60-10.1371/journal.pbio.0060060.
Miller SP, Lunzer M, Dean AM: Direct demonstration of an adaptive constraint. Science. 2006, 314: 458-461. 10.1126/science.1133479.
Dekel E, Alon U: Optimality and evolutionary tuning of the expression level of a protein. Nature. 2005, 436: 588-592. 10.1038/nature03842.
Poelwijk FJ, Heyning PD, de Vos MG, Kiviet DJ, Tans SJ: Optimality and evolution of transcriptionally regulated gene expression. BMC Syst Biol. 2011, 5: 128-10.1186/1752-0509-5-128.
Wang IN: Lysis timing and bacteriophage fitness. Genetics. 2006, 172: 17-26.
Shao Y, Wang IN: Bacteriophage adsorption rate and optimal lysis time. Genetics. 2008, 180: 471-482. 10.1534/genetics.108.090100.
Heineman RH, Bull JJ: Testing optimality with experimental evolution: lysis time in a bacteriophage. Evolution Int J Org Evolution. 2007, 61: 1695-1709. 10.1111/j.1558-5646.2007.00132.x.
Hutchinson CA, Sinsheimer RL: The process of infection with bacteriophage phiX174. J Mol Biol. 1966, 18: 429-447. 10.1016/S0022-2836(66)80035-9.
Josslin R: The lysis mechanism of phage T4: mutants affecting lysis. Virology. 1970, 40: 719-726. 10.1016/0042-6822(70)90216-3.
Heineman RH, Molineux IJ, Bull JJ: Evolutionary robustness of an optimal phenotype: Re-evolution of lysis in a bacteriophage deleted for its lysin gene. J Mol Evol. 2005, 61: 181-191. 10.1007/s00239-004-0304-4.
Miyazaki JI, Ryo Y, Fujisawa H, Minagawa T: Mutation in bacteriophage T3 affecting host cell lysis. Virology. 1978, 89: 327-329. 10.1016/0042-6822(78)90067-3.
Wang IN, Dykhuizen DE, Slobodkin LB: The evolution of phage lysis timing. Evol Ecol. 1996, 10: 545-558. 10.1007/BF01237884.
Bull JJ: Optimality models of phage life history and parallels in disease evolution. J Theor Biol. 2006, 241: 928-938. 10.1016/j.jtbi.2006.01.027.
Charnov EL: Optimal foraging, the marginal value theorem. Theor Popul Biol. 1976, 9: 129-136. 10.1016/0040-5809(76)90040-X.
Abedon ST, Hyman P, Thomas C: Experimental examination of bacteriophage latent-period evolution as a response to bacterial availability. Appl Environ Microbiol. 2003, 69: 7499-7506. 10.1128/AEM.69.12.7499-7506.2003.
Zheng Y, Struck DK, Dankenbring CA, Young R: Evolutionary dominance of holin lysis systems derives from superior genetic malleability. Microbiology. 2008, 154: 1710-1718. 10.1099/mic.0.2008/016956-0.
Gillam S, Astell CR, Jahnke P, Hutchison CA: Smith M: Construction and properties of a ribosome-binding site mutation in gene E of phi X174 bacteriophage. J Virol. 1984, 52: 892-896.
Rokyta DR, Burch CL, Caudle SB, Wichman HA: Horizontal gene transfer and the evolution of microvirid coliphage genomes. J Bacteriol. 2006, 188: 1134-1142. 10.1128/JB.188.3.1134-1142.2006.
Maratea D, Young K, Young R: Deletion and fusion analysis of the phage phi X174 lysis gene E. Gene. 1985, 40: 39-46. 10.1016/0378-1119(85)90022-8.
Young KD, Young R: Lytic action of cloned phi X174 gene E. J Virol. 1982, 44: 993-1002.
Dokland T, McKenna R, Ilag LL, Bowman BR, Incardona NL, Fane BA, Rossmann MG: Structure of a viral procapsid with molecular scaffolding. Nature. 1997, 389: 308-313. 10.1038/38537.
Young R, Wang I, Roof WD: Phages will out: strategies of host cell lysis. Trends Microbiol. 2000, 8: 120-128. 10.1016/S0966-842X(00)01705-4.
Bernhardt TG, Struck DK, Young R: The lysis protein E of phi X174 is a specific inhibitor of the MraY-catalyzed step in peptidoglycan synthesis. J Biol Chem. 2001, 276: 6093-6097. 10.1074/jbc.M007638200.
Wang IN, Smith DL, Young R: Holins: the protein clocks of bacteriophage infections. Annu Rev Microbiol. 2000, 54: 799-825. 10.1146/annurev.micro.54.1.799.
Bernhardt TG, Roof WD, Young R: Genetic evidence that the bacteriophage phi X174 lysis protein inhibits cell wall synthesis. Proc Natl Acad Sci USA. 2000, 97: 4297-4302. 10.1073/pnas.97.8.4297.
Zheng Y, Struck DK, Young R: Purification and functional characterization of phiX174 lysis protein E. Biochemistry. 2009, 48: 4999-5006. 10.1021/bi900469g.
Boyle DS, Donachie WD: mraY is an essential gene for cell growth in Escherichia coli. J Bacteriol. 1998, 180: 6429-6432.
Roof WD, Young R: Phi X174 lysis requires slyD, a host gene which is related to the FKBP family of peptidyl-prolyl cis-trans isomerases. FEMS Microbiol Rev. 1995, 17: 213-218.
Roof WD, Horne SM, Young KD, Young R: slyD, a host gene required for phi X174 lysis, is related to the FK506-binding protein family of peptidyl-prolyl cis-trans-isomerases. J Biol Chem. 1994, 269: 2902-2910.
Bernhardt TG, Roof WD, Young R: The Escherichia coli FKBP-type PPIase SlyD is required for the stabilization of the E lysis protein of bacteriophage phi X174. Mol Microbiol. 2002, 45: 99-108. 10.1046/j.1365-2958.2002.02984.x.
Bull JJ, Heineman RH, Wilke CO: The phenotype-fitness map in experimental evolution of phages. PLoS One. 2011, 6: e27796-10.1371/journal.pone.0027796.
Bull JJ, Badgett MR, Wichman HA: Big-benefit mutations in a bacteriophage inhibited with heat. Mol Biol Evol. 2000, 17: 942-950.
Bull JJ, Badgett MR, Springman R, Molineux IJ: Genome properties and the limits of adaptation in bacteriophages. Evolution Int J Org Evolution. 2004, 58: 692-701.
Guyader S, Burch CL: Optimal foraging predicts the ecology but not the evolution of host specialization in bacteriophages. PLoS One. 2008, 3: e1946-10.1371/journal.pone.0001946.
Heineman RH, Springman R, Bull JJ: Optimal Foraging by Bacteriophages through Host Avoidance. Am Nat. 2008, 171: E149-E157. 10.1086/528962.
Duffy S, Turner PE, Burch CL: Pleiotropic costs of niche expansion in the RNA bacteriophage deleted for its lysin gene. Genetics. 2006, 172: 751-757.
Fisher RA: The Genetical Theory of Natural Selection. 1930, Oxford: Oxford University Press
Orr HA: The evolutionary genetics of adaptation: a simulation study. Genet Res. 1999, 74: 207-214. 10.1017/S0016672399004164.
Heineman RH, Bull JJ, Molineux IJ: Layers of evolvability in a bacteriophage life history trait. Mol Biol Evol. 2009, 26: 1289-1298. 10.1093/molbev/msp037.
Young R, Wang I: Phage Lysis. The Bacteriophage. Edited by: Calender R. 2005, Oxford: Oxford University Press, 104-125. 2
Bowes JM, Dowell CE: Purification and some properties of bacteriophage ST-1. J Virol. 1974, 13: 53-61.
Bradley DE: A comparative study of some properties of the phi-X174 type bacteriophages. Can J Microbiol. 1970, 16: 965-971. 10.1139/m70-165.
Baba T, Ara T, Hasegawa M, Takai Y, Okumura Y, Baba M, Datsenko KA, Tomita M, Wanner BL, Mori H: Construction of Escherichia coli K-12 in-frame, single-gene knockout mutants: the Keio collection. Mol Syst Biol. 2006, 2: 2006.0008
Scott JF, Eisenberg S, Bertsch LL, Kornberg A: A mechanism of duplex DNA replication revealed by enzymatic studies of phage phi X174: catalytic strand separation in advance of replication. Proc Natl Acad Sci USA. 1977, 74: 193-197. 10.1073/pnas.74.1.193.
Fane BA, Hayashi M: Second-site suppressors of a cold-sensitive prohead accessory protein of bacteriophage phi X174. Genetics. 1991, 128: 663-671.
Uchiyama A, Heiman P, Fane BA: N-terminal deletions of the phiX174 external scaffolding protein affect the timing and fidelity of assembly. Virology. 2009, 386: 303-309. 10.1016/j.virol.2009.01.030.
Dennehy JJ, Wang IN: Factors influencing lysis time stochasticity in bacteriophage lambda. BMC Microbiol. 2011, 11: 174-10.1186/1471-2180-11-174.
Fane BA, Brentlinger KL, Burch AD, Chen M, Hafenstein S, Moore E, Novak CR, Uchiyama A: phiX174 et al., the Microviridae. The Bacteriophages. Edited by: Abedon ST. 2006, New York: Oxford University Press (US), 129-145.
Hayashi M, Aoyama A, Richardson DL, Hayashi MN: Biology of the Bacteriophage phiX174. The Bacteriophages. Edited by: Calender R. 1988, New York: Plenum Press, 2:
Ewald PW: Evolution of Infectious Disease. 1996, Oxford: Oxford University Press
Gandon S, Mackinnon M, Nee S, Read A: Imperfect vaccination: some epidemiological and evolutionary consequences. Proc Biol Sci. 2003, 270: 1129-1136. 10.1098/rspb.2003.2370.
Gandon S, Mackinnon MJ, Nee S, Read AF: Imperfect vaccines and the evolution of pathogen virulence. Nature. 2001, 414: 751-756. 10.1038/414751a.
Mackinnon MJ, Gandon S, Read AF: Virulence evolution in response to vaccination: the case of malaria. Vaccine. 2008, 26 (Suppl 3): C42-C52.
We thank J.J. Bull for vital assistance and H.A. Wichman for ΦX174 background. I.N. Wang provided us with early drafts of a relevant paper. T.E. Keller provided coding assistance. W.R. Harcombe, R. Springman, and M.L. Wayne gave edits across multiple drafts. I. Molineux provided us with ST-1 bacteriophage. The National BioResource Project (NIG, Japan): E. coli provided the Keio collection lines used for ST-1 adaptation. This work was funded by NIH Grant Number 57756 to J.J. Bull, HHMI Award #52005907 (Undergraduate Science Education) to the Freshman Research Initiative, and Grant Number T32DE007288 from the National Institute of Dental and Craniofacial Research. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institute of Dental & Craniofacial Research, the National Institutes of Health, or the Howard Hughes Medical Institute.
The authors declare that they have no competing interests.
LC performed the adaptations, sequencing, the majority of the phenotypic assays, and helped draft the manuscript. RH performed calculations of predicted optima, phage release assays with preadsorption, adsorption assays, and drafted the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Release curves for ST-1+ (solid line) and ST-1K-12 (dotted line). (A) Samples treated with lysis solution containing 2 mg/ml lysozyme and chloroform, used to estimate eclipse. (B) Untreated samples, used to estimate lysis and lysis interval. (PDF 55 KB)
About this article
Cite this article
Chantranupong, L., Heineman, R.H. A common, non-optimal phenotypic endpoint in experimental adaptations of bacteriophage lysis time. BMC Evol Biol 12, 37 (2012). https://doi.org/10.1186/1471-2148-12-37
- Experimental evolution
- Life history evolution
- Genetic constraint
- Molecular evolution
- Phenotype prediction