A quantifiable phenotype of viral propagation

The dynamics of phage predation on a microcolony

Phage predation is an important factor for controlling the bacterial biomass. At face value, dense microbial habitats are expected to be vulnerable to phage epidemics due to the abundance of fresh hosts immediately next to any infected bacteria. Despite this, the bacterial microcolony is a common habitat for bacteria in nature. Here, we experimentally quantify the fate of microcolonies of Escherichia coli exposed to virulent phage T4. It has been proposed that the outer bacterial layers of the colony will shield the inner layers from the phage invasion and thereby constrain the phage to the colony's surface. We develop a dynamical model that incorporates this shielding mechanism and fit the results with experimental measurements to extract important phage-bacteria interaction parameters. The analysis suggests that, while the shielding mechanism delays phage attack, T4 phage are able to diffuse so deep into the dense bacterial environment that colony-level survival of the bacterial community is challenged.

Computational multigene interactions in virus growth and infection spread

Viruses persist in nature owing to their extreme genetic heterogeneity and large population sizes, which enable them to evade host immune defenses, escape antiviral drugs, and adapt to new hosts. The persistence of viruses is challenging to study because mutations affect multiple virus genes, interactions among genes in their impacts on virus growth are seldom known, and measures of viral fitness are yet to be standardized. To address these challenges, we employed a data-driven computational model of cell infection by a virus. The infection model accounted for the kinetics of viral gene expression, functional gene-gene interactions, genome replication, and allocation of host cellular resources to produce progeny of vesicular stomatitis virus, a prototype RNA virus. We used this model to computationally probe how interactions among genes carrying up to eleven deleterious mutations affect different measures of virus fitness: single-cycle growth yields and multicycle rates of infection spread. Individual mutations were implemented by perturbing biophysical parameters associated with individual gene functions of the wild-type model. Our analysis revealed synergistic epistasis among deleterious mutations in their effects on virus yield; so adverse effects of single deleterious mutations were amplified by interaction. For the same mutations, multicycle infection spread indicated weak or negligible epistasis, where single mutations act alone in their effects on infection spread. These results were robust to simulation in high- and low-host resource environments. Our work highlights how different types and magnitudes of epistasis can arise for genetically identical virus variants, depending on the fitness measure. More broadly, gene-gene interactions can differently affect how viruses grow and spread.

Inhibition of host cell division by T5 protein 008 (Hdi)

IMPORTANCE

We have identified a novel phage-encoded inhibitor of the major cytoskeletal protein in bacterial division, FtsZ. The inhibition is shown to confer T5 bacteriophage with a growth advantage in dividing hosts. Our studies demonstrate a strategy in bacteriophages to maximize their progeny number by inhibiting escape of one of the daughter cells of an infected bacterium. They further emphasize that FtsZ is a natural target for bacterial growth inhibition.

Ecology and Evolutionary Biology of Hindering Phage Therapy: The Phage Tolerance vs. Phage Resistance of Bacterial Biofilms

As with antibiotics, we can differentiate various acquired mechanisms of bacteria-mediated inhibition of the action of bacterial viruses (phages or bacteriophages) into ones of tolerance vs. resistance. These also, respectively, may be distinguished as physiological insensitivities (or protections) vs. resistance mutations, phenotypic resistance vs. genotypic resistance, temporary vs. more permanent mechanisms, and ecologically vs. also near-term evolutionarily motivated functions. These phenomena can result from multiple distinct molecular mechanisms, many of which for bacterial tolerance of phages are associated with bacterial biofilms (as is also the case for the bacterial tolerance of antibiotics). The resulting inhibitions are relevant from an applied perspective because of their potential to thwart phage-based treatments of bacterial infections, i.e., phage therapies, as well as their potential to interfere more generally with approaches to the phage-based biological control of bacterial biofilms. In other words, given the generally low toxicity of properly chosen therapeutic phages, it is a combination of phage tolerance and phage resistance, as displayed by targeted bacteria, that seems to represent the greatest impediments to phage therapy's success. Here I explore general concepts of bacterial tolerance of vs. bacterial resistance to phages, particularly as they may be considered in association with bacterial biofilms.

Phage co-transport with hyphal-riding bacteria fuels bacterial invasion in a water-unsaturated microbial model system

Nonmotile microorganisms often enter new habitats by co-transport with motile microorganisms. Here, we report that also lytic phages can co-transport with hyphal-riding bacteria and facilitate bacterial colonization of a new habitat. This is comparable to the concept of biological invasions in macroecology. In analogy to invasion frameworks in plant and animal ecology, we tailored spatially organized, water-unsaturated model microcosms using hyphae of Pythium ultimum as invasion paths and flagellated soil-bacterium Pseudomonas putida KT2440 as carrier for co-transport of Escherichia virus T4. P. putida KT2440 efficiently dispersed along P. ultimum to new habitats and dispatched T4 phages across air gaps transporting ≈0.6 phages bacteria⁻¹. No T4 displacement along hyphae was observed in the absence of carrier bacteria. If E. coli occupied the new habitat, T4 co-transport fueled the fitness of invading P. putida KT2440, while the absence of phage co-transport led to poor colonization followed by extinction. Our data emphasize the importance of hyphal transport of bacteria and associated phages in regulating fitness and composition of microbial populations in water-unsaturated systems. As such co-transport seems analogous to macroecological invasion processes, hyphosphere systems with motile bacteria and co-transported phages could be useful models for testing hypotheses in invasion ecology.

Sustainability of spatially distributed bacteria-phage systems

Virulent phages can expose their bacterial hosts to devastating epidemics, in principle leading to complete elimination of their hosts. Although experiments indeed confirm a large reduction of susceptible bacteria, there are no reports of complete extinctions. We here address this phenomenon from the perspective of spatial organization of bacteria and how this can influence the final survival of them. By modelling the transient dynamics of bacteria and phages when they are introduced into an environment with finite resources, we quantify how time delayed lysis, the spatial separation of initial bacterial positions, and the self-protection of bacteria growing in spherical colonies favour bacterial survival. Our results suggest that spatial structures on the millimetre and submillimetre scale play an important role in maintaining microbial diversity.

Phage "delay" towards enhancing bacterial escape from biofilms: a more comprehensive way of viewing resistance to bacteriophages

In exploring bacterial resistance to bacteriophages, emphasis typically is placed on those mechanisms which completely prevent phage replication. Such resistance can be detected as extensive reductions in phage ability to form plaques, that is, reduced efficiency of plating. Mechanisms include restriction-modification systems, CRISPR/Cas systems, and abortive infection systems. Alternatively, phages may be reduced in their “vigor” when infecting certain bacterial hosts, that is, with phages displaying smaller burst sizes or extended latent periods rather than being outright inactivated. It is well known, as well, that most phages poorly infect bacteria that are less metabolically active. Extracellular polymers such as biofilm matrix material also may at least slow phage penetration to bacterial surfaces. Here I suggest that such “less-robust” mechanisms of resistance to bacteriophages could serve bacteria by slowing phage propagation within bacterial biofilms, that is, delaying phage impact on multiple bacteria rather than necessarily outright preventing such impact. Related bacteria, ones that are relatively near to infected bacteria, e.g., roughly 10+ µm away, consequently may be able to escape from biofilms with greater likelihood via standard dissemination-initiating mechanisms including erosion from biofilm surfaces or seeding dispersal/central hollowing. That is, given localized areas of phage infection, so long as phage spread can be reduced in rate from initial points of contact with susceptible bacteria, then bacterial survival may be enhanced due to bacteria metaphorically “running away” to more phage-free locations. Delay mechanisms—to the extent that they are less specific in terms of what phages are targeted—collectively could represent broader bacterial strategies of phage resistance versus outright phage killing, the latter especially as require specific, evolved molecular recognition of phage presence. The potential for phage delay should be taken into account when developing protocols of phage-mediated biocontrol of biofilm bacteria, e.g., as during phage therapy of chronic bacterial infections.

Bacteriophage T4 Infection of Stationary Phase E. coli: Life after Log from a Phage Perspective

Virtually all studies of phage infections investigate bacteria growing exponentially in rich media. In nature, however, phages largely encounter non-growing cells. Bacteria entering stationary phase often activate well-studied stress defense mechanisms that drastically alter the cell, facilitating its long-term survival. An understanding of phage-host interactions in such conditions is of major importance from both an ecological and therapeutic standpoint. Here, we show that bacteriophage T4 can efficiently bind to, infect and kill E. coli in stationary phase, both in the presence and absence of a functional stationary-phase sigma factor, and explore the response of T4-infected stationary phase cells to the addition of fresh nutrients 5 or 24 h after that infection. An unexpected new mode of response has been identified. "Hibernation" mode is a persistent but reversible dormant state in which the infected cells make at least some phage enzymes, but halt phage development until appropriate nutrients become available before producing phage particles. Our evidence indicates that the block in hibernation mode occurs after the middle-mode stage of phage development; host DNA breakdown and the incorporation of the released nucleotides into phage DNA indicate that the enzymes of the nucleotide synthesizing complex, under middle-mode control, have been made and assembled into a functional state. Once fresh glucose and amino acids become available, the standard lytic infection process rapidly resumes and concentrations of up to 10(11) progeny phage (an average of about 40 phage per initially present cell) are produced. All evidence is consistent with the hibernation-mode control point lying between middle mode and late mode T4 gene expression. We have also observed a "scavenger" response, where the infecting phage takes advantage of whatever few nutrients are available to produce small quantities of progeny within 2 to 5 h after infection. The scavenger response seems able to produce no more than an average of one phage per originally available cell, and few if any further progeny are produced by cells in this mode even if fresh nutrients are made available later.

Population Dynamics of Phage and Bacteria in Spatially Structured Habitats Using Phage λ and Escherichia coli

Bacteria living in physically structured habitats are exposed heterogeneously to both resources and different types of phages. While there have been numerous experimental approaches to examine spatially distributed bacteria exposed to phages, there is little theory to guide the design of these experiments, interpret their results, or expand the inferences drawn to a broader ecological and evolutionary context. Plaque formation provides a window into understanding phage-bacterium interactions in physically structured populations, including surfaces, semisolids, and biofilms. We develop models to address the plaque dynamics for a temperate phage and its virulent mutants. The models are compared with phage λ-Escherichia coli system to quantify their applicability. We found that temperate phages gave an increasing number of gradually smaller colonies as the distance increased from the plaque center. For low-lysogen frequency this resulted in plaques with most of the visible colonies at an intermediate distance between the center and periphery. Using spot inoculation, where phages in excess of bacteria were inoculated in a circular area, we measured the frequency and spatial distribution of lysogens. The spot morphology of cII-negative (cII(-)) and cIII(-) mutants of phage λ displays concentric rings of high-density lysogenic colonies. The simplest of these ring morphologies was reproduced by including multiplicity of infection (MOI) sensitivity in lysis-lysogeny decisions, but its failure to explain the occasional observation of multiple rings in cIII(-) mutant phages highlights unknown features of this phage. Our findings demonstrated advantages of temperate phages over virulent phages in exploiting limited resources in spatially distributed microbial populations.

IMPORTANCE

Phages are the most abundant organisms on earth, and yet little is known about how phages and bacterial hosts are influencing each other in density and evolution. Phages can be either virulent or temperate, a difference that is highlighted when a spatially structured bacterial population is infected. Phage λ is a temperate phage, with a capacity for dormancy that can be modified by single gene knockouts. The stochastic bias in the lysis-lysogeny decision's probability is reflected in plaque morphologies on bacterial lawns. We present a model for plaque morphology of both virulent and temperate phages, taking into account the underlying survival of bacterial microcolonies. It reproduces known plaque morphologies and speaks to advantages of temperate phages in a spatially structured environment.

Virulence evolution at the front line of spreading epidemics

Understanding and predicting the spatial spread of emerging pathogens is a major challenge for the public health management of infectious diseases. Theoretical epidemiology shows that the speed of an epidemic is governed by the life-history characteristics of the pathogen and its ability to disperse. Rapid evolution of these traits during the invasion may thus affect the speed of epidemics. Here we study the influence of virulence evolution on the spatial spread of an epidemic. At the edge of the invasion front, we show that more virulent and transmissible genotypes are expected to win the competition with other pathogens. Behind the front line, however, more prudent exploitation strategies outcompete virulent pathogens. Crucially, even when the presence of the virulent mutant is limited to the edge of the front, the invasion speed can be dramatically altered by pathogen evolution. We support our analysis with individual-based simulations and we discuss the additional effects of demographic stochasticity taking place at the front line on virulence evolution. We confirm that an increase of virulence can occur at the front, but only if the carrying capacity of the invading pathogen is large enough. These results are discussed in the light of recent empirical studies examining virulence evolution at the edge of spreading epidemics.

Phage therapy: eco-physiological pharmacology

Bacterial virus use as antibacterial agents, in the guise of what is commonly known as phage therapy, is an inherently physiological, ecological, and also pharmacological process. Physiologically we can consider metabolic properties of phage infections of bacteria and variation in those properties as a function of preexisting bacterial states. In addition, there are patient responses to pathogenesis, patient responses to phage infections of pathogens, and also patient responses to phage virions alone. Ecologically, we can consider phage propagation, densities, distribution (within bodies), impact on body-associated microbiota (as ecological communities), and modification of the functioning of body "ecosystems" more generally. These ecological and physiological components in many ways represent different perspectives on otherwise equivalent phenomena. Comparable to drugs, one also can view phages during phage therapy in pharmacological terms. The relatively unique status of phages within the context of phage therapy as essentially replicating antimicrobials can therefore result in a confluence of perspectives, many of which can be useful towards gaining a better mechanistic appreciation of phage therapy, as I consider here. Pharmacology more generally may be viewed as a discipline that lies at an interface between organism-associated phenomena, as considered by physiology, and environmental interactions as considered by ecology.

Imaging the propagation of viruses

The propagation of viruses in a growing plaque has been measured using a digital image acquisition and analysis system. Plaques of phage T7 incubated at 37 degrees C and illuminated against a dark field emerged as dark growing spots against a background of host bacteria. Images of the growth were acquired using a charge-coupled device (CCD) camera at 1-h intervals over 24 h. The first 10 h of plaque development coincided with rapid growth of the agar-immobilized Escherichia coli host, measured as a reduction in gray value. Following this period, the average radial velocity of plaque growth remained constant at 0.059 mm/h while the standard deviation about this velocity increased. These results suggest the suitability of the system for spatially resolving the dynamics of viral evolution during plaque growth. (c) 1996 John Wiley & Sons, Inc.

Complex spatial dynamics of oncolytic viruses in vitro: mathematical and experimental approaches

Oncolytic viruses replicate selectively in tumor cells and can serve as targeted treatment agents. While promising results have been observed in clinical trials, consistent success of therapy remains elusive. The dynamics of virus spread through tumor cell populations has been studied both experimentally and computationally. However, a basic understanding of the principles underlying virus spread in spatially structured target cell populations has yet to be obtained. This paper studies such dynamics, using a newly constructed recombinant adenovirus type-5 (Ad5) that expresses enhanced jellyfish green fluorescent protein (EGFP), AdEGFPuci, and grows on human 293 embryonic kidney epithelial cells, allowing us to track cell numbers and spatial patterns over time. The cells are arranged in a two-dimensional setting and allow virus spread to occur only to target cells within the local neighborhood. Despite the simplicity of the setup, complex dynamics are observed. Experiments gave rise to three spatial patterns that we call "hollow ring structure", "filled ring structure", and "disperse pattern". An agent-based, stochastic computational model is used to simulate and interpret the experiments. The model can reproduce the experimentally observed patterns, and identifies key parameters that determine which pattern of virus growth arises. The model is further used to study the long-term outcome of the dynamics for the different growth patterns, and to investigate conditions under which the virus population eliminates the target cells. We find that both the filled ring structure and disperse pattern of initial expansion are indicative of treatment failure, where target cells persist in the long run. The hollow ring structure is associated with either target cell extinction or low-level persistence, both of which can be viewed as treatment success. Interestingly, it is found that equilibrium properties of ordinary differential equations describing the dynamics in local neighborhoods in the agent-based model can predict the outcome of the spatial virus-cell dynamics, which has important practical implications. This analysis provides a first step towards understanding spatial oncolytic virus dynamics, upon which more detailed investigations and further complexity can be built.

Effects of bacteriophage traits on plaque formation

Background

The appearance of plaques on a bacterial lawn is one of the enduring imageries in modern day biology. The seeming simplicity of a plaque has invited many hypotheses and models in trying to describe and explain the details of its formation. However, until now, there has been no systematic experimental exploration on how different bacteriophage (phage) traits may influence the formation of a plaque. In this study, we constructed a series of isogenic λ phages that differ in their adsorption rate, lysis timing, or morphology so that we can determine the effects if these changes on three plaque properties: size, progeny productivity, and phage concentration within plaques.

Results

We found that the adsorption rate has a diminishing, but negative impact on all three plaque measurements. Interestingly, there exists a concave relationship between the lysis time and plaque size, resulting in an apparent optimal lysis time that maximizes the plaque size. Although suggestive in appearance, we did not detect a significant effect of lysis time on plaque productivity. Nonetheless, the combined effects of plaque size and productivity resulted in an apparent convex relationship between the lysis time and phage concentration within plaques. Lastly, we found that virion morphology also affected plaque size. We compared our results to the available models on plaque size and productivity. For the models in their current forms, a few of them can capture the qualitative aspects of our results, but not consistently in both plaque properties.

Conclusions

By using a collection of isogenic phage strains, we were able to investigate the effects of individual phage traits on plaque size, plaque productivity, and average phage concentration in a plaque while holding all other traits constant. The controlled nature of our study allowed us to test several model predictions on plaque size and plaque productivity. It seems that a more realistic theoretical approach to plaque formation is needed in order to capture the complex interaction between phage and its bacterium host in a spatially restricted environment.

Dynamics of virus spread in the presence of fluid flow

The dynamics of viral infection spread, whether in laboratory cultures or in naturally infected hosts, reflects a coupling of biological and physical processes that remain to be fully elucidated. Biological processes include the kinetics of virus growth in infected cells while physical processes include transport of virus progeny from infected cells, where they are produced, to susceptible cells, where they initiate new infections. Mechanistic models of infection spread have been widely developed for systems where virus growth is coupled with transport of virus particles by diffusion, but they have yet to be developed for systems where viruses move under the influence of fluid flows. Recent experimental observations of flow-enhanced infection spread in laboratory cultures motivate here the development of initial continuum and discrete virus-particle models of infection spread. The magnitude of a dimensionless group, the Damköhler number, shows how parameters that characterize particle adsorption to cells, strain rates that reflect flow profiles, and diffusivities of virus particles combine to influence the spatial pattern of infection spread.

Infection on a chip: a microscale platform for simple and sensitive cell-based virus assays

The plaque assay has long served as the "gold standard" to measure virus infectivity and test antiviral drugs, but the assay is labor-intensive, lacks sensitivity, uses excessive reagents, and is hard to automate. Recent modification of the assay to exploit flow-enhanced virus spread with quantitative imaging has increased its sensitivity. Here we performed flow-enhanced infection assays in microscale channels, employing passive fluid pumping to inoculate cell monolayers with virus and drive infection spread. Our test of an antiviral drug (5-fluorouracil) against vesicular stomatitis virus infections of BHK cell monolayers yielded a two-fold improvement in sensitivity, relative to the standard assay based on plaque counting. The reduction in scale, simplified fluid handling, image-based quantification, and higher assay sensitivity will enable infection measurements for high-throughput drug screening, sero-conversion testing, and patient-specific diagnosis of viral infections.

Bacteriophage plaques: theory and analysis

Laboratory characterization of bacteriophage growth traditionally is done either in broth cultures or in semisolid agar media. These two environments may be distinguished in terms of their spatial structure, i.e., the degree to which they limit diffusion, motility, and environmental mixing. Well-mixed broth, for example, represents the microbiological ideal of a non-spatially structured environment. Agar, by contrast, imposes significant limitations on phage and bacterial movement and therefore gives rise to spatial structure. The study of phage growth within spatially structured environments, such as that seen during phage plaque formation, is important for three reasons. First, a large fraction of environmental bacteria live within spatially structured environments such as within biofilms, within soil, or when growing in or on the tissues of plants and animals. Second, phage growth as plaques is a central technique to phage studies, yet appears to be under appreciated by phage workers in terms of its complexity. Third, selective pressures acting on phage during plaque growth differ from those seen during broth growth. In this chapter we will discuss just what a plaque is, how one forms, and what can affect plaque size. We will describe methods, both experimental and theoretical, that have been employed to study plaque growth. As caveats we will discuss why plaque formation failure is not necessarily equivalent to virion inviability (Note 1). We also will consider problems with inferring phage broth growth fitness as a function of plaque characteristics (Note 2).

Bacteriophage evolution given spatial constraint

Spatial structure can impede mixing, diffusion, and motility. In microbiology laboratories, spatial structure is commonly achieved via formation of agar gels, within which bacteriophage (phage) replication results in localized clearings called plaques. Developing a better understanding of phage plaque formation is relevant because of the ubiquity of phage plaquing in the laboratory; because plaque size has been employed as a measure of phage fitness; because many bacteria exist within environments that display significant spatial structure (e.g., biofilms, soils, sediments, and in or on plant or animal tissues); and because spatial structure could impede phage exploitation of bacterial communities. There is, however, a relative dearth of experimentation and analysis considering phage plaque formation from the perspective of selection acting on individual phage growth parameters-latent period, burst size, and adsorption rate. Here we consider the impact of these parameters on rates of plaque wavefront velocity (rates of radial plaque enlargement), especially as functions of existing phage and environmental properties. We do so based on analyses of published equations which predict plaque enlargement rates. These indicate that greater wavefront velocities should be associated with (i) latent period reductions, (ii) larger burst sizes, or (iii) faster virion binding to bacteria. We suggest, however, that deviations could occur, respectively, (i) if virion adsorption is "slow" or if burst sizes are large, (ii) if burst sizes are already large, or (iii) if virion binding rates are already fast, bacterial densities are especially high, or burst sizes are large. Higher initial lawn bacterial densities could also contribute to faster plaque expansion, but only if adsorption is otherwise slow or burst sizes are large. By contrast, faster virion diffusion is always expected to result in greater plaque wavefront velocities. Overall, we provide a snapshot of how phage populations may respond evolutionarily to selection for more-rapid propagation during spatially constrained growth.

Propagation of fluorescent viruses in growing plaques

To study virus propagation, we have developed a method by which the propagation of the Lambda bacteriophage can be observed and quantified. This is done by creating a fusion protein of the capsid protein gpD and the enhanced yellow fluorescent protein (EYFP). We show that this fusion allows capsid formation and that the modified viruses propagate on a surface covered with host bacteria thus forming fluorescent plaques. The intensity of fluorescence in a growing plaque determines the distribution of phages. This provides a new tool to study the propagation of infection at the microscopic level.

Spatial patterns of protein expression in focal infections of human cytomegalovirus

Human cytomegalovirus (HCMV) is a medically significant human pathogen that infects a wide range of cell and tissue types. During infection, HCMV activates a variety of signal transduction pathways that induce profound changes in cellular processes and dramatically affect cellular gene expression patterns. To better define how these virus-host interactions affect the local microenvironment and influence the spatial and temporal spread of HCMV, we initiated HCMV focal infections on normal human dermal fibroblast monolayers and monitored viral gene expression patterns and infection spread over 45 days. To establish baseline temporal measurements of HCMV infection and spread in cell monolayers, we characterized the influence of three experimental variables on viral gene expression: cell plating density, the presence of serum, and neutralization of cellular antiviral responses with an antibody against interferon-beta. We found that high cell plating density or the inclusion of serum correlated with enhanced HCMV infection spread. Dramatic differences in the expression pattern of the viral immediate early 2 (IE2) gene were observed under these conditions as compared to low plating density or the absence of serum. In the latter case round, uniform foci were observed with a clear wave of IE2 expression visible in advance of a late stage viral protein, envelope glycoprotein B. By contrast, larger irregular foci with arms of IE2 expression were observed in the presence of serum. Addition of the antibody had little effect on the rate of spread, which is consistent with the knowledge that HCMV represses antiviral responses during infection. This experimental system provides a useful means to visualize and quantify complex virus-host interactions.

Time-delayed spread of viruses in growing plaques

The spread of viruses in growing plaques predicted by classical models is greater than that measured experimentally. There is a widespread belief that this discrepancy is due to biological factors. Here we show that the observed speeds can be satisfactorily predicted by a purely physical model that takes into account the delay time due to virus reproduction inside infected cells. No free or adjustable parameters are used.

Dispersal probability distributions and the wave-front speed problem

The speed and width of front solutions to reaction-dispersal models are analyzed both analytically and numerically. We perform our analysis for Laplace and Gaussian distribution kernels, both for delayed and nondelayed models. The results are discussed in terms of the characteristic parameters of the models.

Independent contrasts succeed where ancestor reconstruction fails in a known bacteriophage phylogeny

Methods of ancestor reconstruction are important tools for evolutionary inference that are difficult to test empirically because ancestral states are rarely known with certainty. We evaluated reconstruction methods for continuous phenotypic characters using taxa from an experimentally generated bacteriophage phylogeny. Except for one slowly evolving character, the estimated ancestral states of continuous phenotypic characters were highly inaccurate and biased, even when including a known ancestor at the root. This error was caused by a directional trend in character evolution and by rapid rates of character evolution. Computer simulations confirmed that such factors affect reconstruction of continuous characters in general. We also used phenotypic viral characters to evaluate two methods that attempt to estimate the correlation between characters during evolution. Whereas a nonphylogenetic regression was relatively inaccurate and biased, independent contrasts accurately estimated the correlation between characters with little bias.

Amplification and spread of viruses in a growing plaque

The two-dimensional propagation of viruses through a "lawn" of receptive hosts, commonly called plaque growth, reflects the dynamics of interactions between viruses and host cells. Here we treat the amplification of viruses during plaque growth as a reaction-diffusion system, where interactions among the virus, uninfected host cells, and virus-producing host-virus complexes are accounted for using rates of viral adsorption to and desorption from the host-cell surface, rates of reproduction and release of progeny viruses by lysis of the host, and by the coupling of these reactions with diffusion of free virus within the agar support. Numerical solution of the system shows the development of a traveling wave of reproducing viruses, where the velocity of the wave is governed by the kinetic and diffusion parameters. The model has been applied to predict the propagation velocity of a bacteriophage plaque. Different mechanisms may account for the dependence of this velocity on the host density during early stages of a growing plaque. The model provides a means to explore how changes in the virus-host interactions may be manifest in a growing plaque.

Detection of evolving viruses

The spread of viruses on a homogeneous lawn of receptive hosts provides an opportunity to detect the dynamics of their evolution. We have previously found that when repeated virus passages are confined to the expanding perimeter of a growing plaque, the appearance and outgrowth of genetically diverse strains (all descended from the same parent strain) can be traced along different radii of the plaque. As a plaque grows, the random mutation and selection of new fast-growing strains reduce the roundness or circularity of the growing plaque. Here we have quantified such changes in growing plaques of bacteriophage T7 using a digital imaging system. We find that T7 populations not adapted for fast growth exhibit a broader diversity of growth rates than populations adapted for fast growth. These results provide a foundation for understanding how viruses exploit mutation and selection processes to persist in nature.

Whole-virus vaccine development by continuous culture on a complementing host

We have evaluated an adaptive strategy for generating whole-virus vaccines using a bacteriophage model. Wildtype phage T7 was cultivated in a two-stage continuous stirred-tank reactor (CSTR) utilizing a recombinant E. coli host that constitutively expressed T7 RNA polymerase, an essential enzyme of the early viral metabolism. Over the course of 180 generations a diversity of phage variants emerged, outgrew the wildtype, and were subsequently eclipsed by yet fitter variants, based on host-ranges, restriction patterns, and one-step growth responses of isolated clones. The fittest variant, which required complementation by the recombinant host in order to grow, deleted at least 12 percent of its genome and replicated twice as fast as the wildtype. Moreover, this variant was immunogenically indistinguishable from the wildtype, based on cross-reactivities of antisera raised against both. These results suggest the feasibility of the proposed strategy for the development of safe whole-virus vaccines.

Evolution of bacteriophage T7 in a growing plaque

The emergence of mutants during the 10(9)-fold amplification of a bacteriophage was spatially resolved in a growing plaque. When wild-type phage T7 was grown on an Escherichia coli host which expressed an essential early enzyme of the phage infection cycle, the T7 RNA polymerase, mutant phage relying on this enzyme appeared in 10(8) phage replications and outgrew the wild type. Spatial resolution of the selection process was achieved by analyzing stab samples taken along a plaque radius. Different mutants were selected at different rates along different radii of the plaque, based on host range and restriction patterns of the isolates. The mutants deleted up to 11% of their genomes, including the gene for their own RNA polymerase. They gained an advantage over the wild type by replicating more efficiently, as determined by one-step growth cultures.

Replication of viruses in a growing plaque: a reaction-diffusion model

An understanding of the viral replication process commonly referred to as "plaque growth" is developed in the context of a reaction-diffusion model. The interactions among three components: the virus, the healthy host, and the infected host are represented using rates of viral adsorption and desorption to the cell surface, replication and release by host lysis, and diffusion. The solution to the full model reveals a maximum in the dependence of the velocity of viral propagation on its equilibrium adsorption constant, suggesting that conditions can be chosen where viruses which adsorb poorly to their hosts will replicate faster in plaques than those which adsorb well. Analytic expressions for the propagation velocity as a function of the kinetic and diffusion parameters are presented for the limiting cases of equilibrated adsorption, slow adsorption, fast adsorption, and large virus yields. Hindered diffusion at high host concentrations must be included for quantitative agreement with experimental data.