Evolution of bacteriophage T7 in a growing plaque

Virus-like Particles: Measures and Biological Functions

Virus-like particles resemble infectious virus particles in size, shape, and molecular composition; however, they fail to productively infect host cells. Historically, the presence of virus-like particles has been inferred from total particle counts by microscopy, and infectious particle counts or plaque-forming-units (PFUs) by plaque assay; the resulting ratio of particles-to-PFUs is often greater than one, easily 10 or 100, indicating that most particles are non-infectious. Despite their inability to hijack cells for their reproduction, virus-like particles and the defective genomes they carry can exhibit a broad range of behaviors: interference with normal virus growth during co-infections, cell killing, and activation or inhibition of innate immune signaling. In addition, some virus-like particles become productive as their multiplicities of infection increase, a sign of cooperation between particles. Here, we review established and emerging methods to count virus-like particles and characterize their biological functions. We take a critical look at evidence for defective interfering virus genomes in natural and clinical isolates, and we review their potential as antiviral therapeutics. In short, we highlight an urgent need to better understand how virus-like genomes and particles interact with intact functional viruses during co-infection of their hosts, and their impacts on the transmission, severity, and persistence of virus-associated diseases.

Historical Perspective on the Discovery of the Quasispecies Concept

Viral quasispecies are dynamic distributions of nonidentical but closely related mutant and recombinant viral genomes subjected to a continuous process of genetic variation, competition, and selection that may act as a unit of selection. The quasispecies concept owes its theoretical origins to a model for the origin of life as a collection of mutant RNA replicators. Independently, experimental evidence for the quasispecies concept was obtained from sampling of bacteriophage clones, which revealed that the viral populations consisted of many mutant genomes whose frequency varied with time of replication. Similar findings were made in animal and plant RNA viruses. Quasispecies became a theoretical framework to understand viral population dynamics and adaptability. The evidence came at a time when mutations were considered rare events in genetics, a perception that was to change dramatically in subsequent decades. Indeed, viral quasispecies was the conceptual forefront of a remarkable degree of biological diversity, now evident for cell populations and organisms, not only for viruses. Quasispecies dynamics unveiled complexities in the behavior of viral populations,with consequences for disease mechanisms and control strategies. This review addresses the origin of the quasispecies concept, its major implications on both viral evolution and antiviral strategies, and current and future prospects.

Biological Characterization and Evolution of Bacteriophage T7-△holin During the Serial Passage Process

Bacteriophage T7 gene 17.5 coding for the only known holin is one of the components of its lysis system, but the holin activity in T7 is more complex than a single gene, and evidence points to the existence of additional T7 genes with holin activity. In this study, a T7 phage with a gene 17.5 deletion (T7-△holin) was rescued and its biological characteristics and effect on cell lysis were determined. Furthermore, the genomic evolution of mutant phage T7-△holin during serial passage was assessed by whole-genome sequencing analysis. It was observed that deletion of gene 17.5 from phage T7 delays lysis time and enlarges the phage burst size; however, this biological characteristic recovered to normal lysis levels during serial passage. Scanning electron microscopy showed that the two opposite ends of E. coli BL21 cells swell post-T7-△holin infection rather than drilling holes on cell membrane when compared with T7 wild-type infection. No visible progeny phage particle accumulation was observed inside the E. coli BL21 cells by transmission electron microscopy. Following serial passage of T7-△holin from the 1st to 20th generations, the mRNA levels of gene 3.5 and gene 19.5 were upregulated and several mutation sites were discovered, especially two missense mutations in gene 19.5, which indicate a potential contribution to the evolution of the T7-△holin. Although the burst size of T7-△holin increased, high titer cultivation of T7-△holin was not achieved by optimizing the culture process. Accordingly, these results suggest that gene 19.5 is a potential lysis-related component that needs to be studied further and that the T7-△holin strain with its gene 17.5 deletion is not adequate to establish the high-titer phage cultivation process.

Patterns of virus growth across the diversity of life

Although viruses in their natural habitats add up to less than 10% of the biomass, they contribute more than 90% of the genome sequences [1]. These viral sequences or 'viromes' encode viruses that populate the Earth's oceans [2, 3] and terrestrial environments [4, 5], where their infections impact life across diverse ecological niches and scales [6, 7], including humans [8-10]. Most viruses have yet to be isolated and cultured [11-13], and surprisingly few efforts have explored what analysis of available data might reveal about their nature. Here, we compiled and analyzed seven decades of one-step growth and other data for viruses from six major families, including their infections of archaeal, bacterial and eukaryotic hosts [14-191]. We found that the use of host cell biomass for virus production was highest for archaea at 10%, followed by bacteria at 1% and eukarya at 0.01%, highlighting the degree to which viruses of archaea and bacteria exploit their host cells. For individual host cells, the yield of virus progeny spanned a relatively narrow range (10-1000 infectious particles per cell) compared with the million-fold difference in size between the smallest and largest cells. Furthermore, healthy and infected host cells were remarkably similar in the time they needed to multiply themselves or their virus progeny. Specifically, the doubling time of healthy cells and the delay time for virus release from infected cells were not only correlated (r = 0.71, p < 10-10, n = 101); they also spanned the same range from tens of minutes to about a week. These results have implications for better understanding the growth, spread and persistence of viruses in complex natural habitats that abound with diverse hosts, including humans and their associated microbes.

The RNA-binding protein Hfq assembles into foci-like structures in nitrogen starved Escherichia coli

The initial adaptive responses to nutrient depletion in bacteria often occur at the level of gene expression. Hfq is an RNA-binding protein present in diverse bacterial lineages that contributes to many different aspects of RNA metabolism during gene expression. Using photoactivated localization microscopy and single-molecule tracking, we demonstrate that Hfq forms a distinct and reversible focus-like structure in Escherichia coli specifically experiencing long-term nitrogen starvation. Using the ability of T7 phage to replicate in nitrogen-starved bacteria as a biological probe of E. coli cell function during nitrogen starvation, we demonstrate that Hfq foci have a role in the adaptive response of E. coli to long-term nitrogen starvation. We further show that Hfq foci formation does not depend on gene expression once nitrogen starvation has set in and occurs indepen-dently of the transcription factor N-regulatory protein C, which activates the initial adaptive response to N starvation in E. coli These results serve as a paradigm to demonstrate that bacterial adaptation to long-term nutrient starvation can be spatiotemporally coordinated and can occur independently of de novo gene expression during starvation.

Gut virome of mammals and birds reveals high genetic diversity of the family Microviridae

Nineteen families of phages infecting bacteria or archaea are currently recognized by the International Committee on Taxonomy of Viruses (ICTV). Of these, only two have single-stranded DNA genomes, namely Inoviridae and Microviridae. The distribution, genetic characteristics, and ecological roles of Microviridae remain largely under explored. Here, using viral metagenomics, we investigate the intestinal virome from human and twenty-four species of animals, as well as freshwater samples, containing abundant sequence reads showing similarity to the Microviridae. Eight hundred and sixty complete or near complete Microviridae-related genomes were generated, showing high levels of co-infections and sequence divergence. Sequence comparison and phylogenetic analysis showed that the Microviridae subfamily Gokushovirinae was highly prevalent and that some strains may qualify as new subfamilies. This study significantly augments our knowledge of the genetic diversity, genome evolution, and distribution in animal species of members of the family Microviridae.

In plaque-mass spectrometry imaging of a bloom-forming alga during viral infection reveals a metabolic shift towards odd-chain fatty acid lipids

Tapping into the metabolic cross-talk between a host and its virus can reveal unique strategies employed during infection. Viral infection is a dynamic process that generates an evolving metabolic landscape. Gaining a continuous view into the infection process is highly challenging and is limited by current metabolomics approaches, which typically measure the average of the entire population at various stages of infection. Here, we took an innovative approach to study the metabolic basis of host-virus interactions between the bloom-forming alga Emiliania huxleyi and its specific virus. We combined a classical method in virology, the plaque assay, with advanced mass spectrometry imaging (MSI), an approach we termed ‘in plaque-MSI’. Taking advantage of the spatial characteristics of the plaque, we mapped the metabolic landscape induced during infection in a high spatiotemporal resolution, unfolding the infection process in a continuous manner. Further unsupervised spatially-aware clustering, combined with known lipid biomarkers, revealed a systematic metabolic shift during infection towards lipids containing the odd-chain fatty acid pentadecanoic acid (C15:0). Applying ‘in plaque-MSI’ might facilitate the discovery of bioactive compounds that mediate the chemical arms race of host-virus interactions in diverse model systems.

Kinetic Modeling of Virus Growth in Cells

When a virus infects a host cell, it hijacks the biosynthetic capacity of the cell to produce virus progeny, a process that may take less than an hour or more than a week. The overall time required for a virus to reproduce depends collectively on the rates of multiple steps in the infection process, including initial binding of the virus particle to the surface of the cell, virus internalization and release of the viral genome within the cell, decoding of the genome to make viral proteins, replication of the genome, assembly of progeny virus particles, and release of these particles into the extracellular environment. For a large number of virus types, much has been learned about the molecular mechanisms and rates of the various steps. However, in only relatively few cases during the last 50 years has an attempt been made-using mathematical modeling-to account for how the different steps contribute to the overall timing and productivity of the infection cycle in a cell. Here we review the initial case studies, which include studies of the one-step growth behavior of viruses that infect bacteria (Qβ, T7, and M13), human immunodeficiency virus, influenza A virus, poliovirus, vesicular stomatitis virus, baculovirus, hepatitis B and C viruses, and herpes simplex virus. Further, we consider how such models enable one to explore how cellular resources are utilized and how antiviral strategies might be designed to resist escape. Finally, we highlight challenges and opportunities at the frontiers of cell-level modeling of virus infections.

T7 phage factor required for managing RpoS in Escherichia coli

T7 development in Escherichia coli requires the inhibition of the housekeeping form of the bacterial RNA polymerase (RNAP), Eσ⁷⁰, by two T7 proteins: Gp2 and Gp5.7. Although the biological role of Gp2 is well understood, that of Gp5.7 remains to be fully deciphered. Here, we present results from functional and structural analyses to reveal that Gp5.7 primarily serves to inhibit Eσ^S, the predominant form of the RNAP in the stationary phase of growth, which accumulates in exponentially growing E. coli as a consequence of the buildup of guanosine pentaphosphate [(p)ppGpp] during T7 development. We further demonstrate a requirement of Gp5.7 for T7 development in E. coli cells in the stationary phase of growth. Our finding represents a paradigm for how some lytic phages have evolved distinct mechanisms to inhibit the bacterial transcription machinery to facilitate phage development in bacteria in the exponential and stationary phases of growth.

Spatial-Temporal Patterns of Viral Amplification and Interference Initiated by a Single Infected Cell

When viruses infect their host cells, they can make defective virus-like particles along with intact virus. Cells coinfected with virus and defective particles often exhibit interference with virus growth caused by the competition for resources by defective genomes. Recent reports of the coexistence and cotransmission of such defective interfering particles (DIPs) in vivo, across epidemiological length and time scales, suggest a role in viral pathogenesis, but it is not known how DIPs impact infection spread, even under controlled culture conditions. Using fluorescence microscopy, we quantified coinfections of vesicular stomatitis virus (VSV) expressing a fluorescent reporter protein and its DIPs on BHK-21 host cell monolayers. We found that viral gene expression was more delayed, infections spread more slowly, and patterns of spread became more “patchy” with higher DIP inputs to the initial cell. To examine how infection spread might depend on the behavior of the initial coinfected cell, we built a computational model, adapting a cellular automaton (CA) approach to incorporate kinetic data on virus growth for the first time. Specifically, changes in observed patterns of infection spread could be directly linked to previous high-throughput single-cell measures of virus-DIP coinfection. The CA model also provided testable hypotheses on the spatial-temporal distribution of the DIPs, which remain governed by their predator-prey interaction. More generally, this work offers a data-driven computational modeling approach for better understanding of how single infected cells impact the multiround spread of virus infections across cell populations.

IMPORTANCE Defective interfering particles (DIPs) compete with intact virus, depleting host cell resources that are essential for virus growth and infection spread. However, it is not known how such competition, strong or weak, ultimately affects the way in which infections spread and cause disease. In this study, we address this unmet need by developing an integrated experimental-computational approach, which sheds new light on how infections spread. We anticipate that our approach will also be useful in the development of DIPs as therapeutic agents to manage the spread of viral infections.

Front propagation speeds of T7 virus mutants

We propose a new reaction-diffusion model with an eclipse time to study the spread of viruses on bacterial populations. This new model is both biologically and physically sound, unlike previous ones. We determine important parameter values from experimental data, such as the one-step growth. We verify the proposed model by comparing theoretical and experimental data of the front propagation speed for several T7 virus strains.

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.

RNA editing by T7 RNA polymerase bypasses InDel mutations causing unexpected phenotypic changes

DNA-dependent T7 RNA polymerase (T7 RNAP) is the most powerful tool for both gene expression and in vitro transcription. By using a Next Generation Sequencing (NGS) approach we have analyzed the polymorphism of a T7 RNAP-generated mRNA pool of the mboIIM2 gene. We find that the enzyme displays a relatively high level of template-dependent transcriptional infidelity. The nucleotide misincorporations and multiple insertions in A/T-rich tracts of homopolymers in mRNA (0.20 and 0.089%, respectively) cause epigenetic effects with significant impact on gene expression that is disproportionally high to their frequency of appearance. The sequence-dependent rescue of single and even double InDel frameshifting mutants and wild-type phenotype recovery is observed as a result. As a consequence, a heterogeneous pool of functional and non-functional proteins of almost the same molecular mass is produced where the proteins are indistinguishable from each other upon ordinary analysis. We suggest that transcriptional infidelity as a general feature of the most effective RNAPs may serve to repair and/or modify a protein function, thus increasing the repertoire of phenotypic variants, which in turn has a high evolutionary potential.

Pseudolysogeny

Pseudolysogeny can be defined as the stage of stalled development of a bacteriophage in a host cell without either multiplication of the phage genome (as in lytic development) or its replication synchronized with the cell cycle and stable maintenance in the cell line (as in lysogenization), which proceeds with no viral genome degradation, thus allowing the subsequent restart of virus development. This phenomenon is usually caused by unfavorable growth conditions for the host cell (such as starvation) and is terminated with initiation of either true lysogenization or lytic growth when growth conditions improve. Pseudolysogeny has been known for tens of years; however, its role has often been underestimated. Currently, it is being considered more often as an important aspect of phage-host interactions. The reason for this is mostly an increased interest in phage-host interactions in the natural environment. Pseudolysogeny seems to play an important role in phage survival, as bacteria in a natural environment are starved or their growth is very slow. This phenomenon can be an important aspect of phage-dependent bacterial mortality and may influence the virulence of some bacterial strains.

Genetic hitchhiking in spatially extended populations

When a mutation with selective advantage s spreads through a panmictic population, it may cause two lineages at a linked locus to coalesce; the probability of coalescence is exp(-2rT), where T∼log(2Ns)/s is the time to fixation, N is the number of haploid individuals, and r is the recombination rate. Population structure delays fixation, and so weakens the effect of a selective sweep. However, favourable alleles spread through a spatially continuous population behind a narrow wavefront; ancestral lineages are confined at the tip of this front, and so coalesce rapidly. In extremely dense populations, coalescence is dominated by rare fluctuations ahead of the front. However, we show that for moderate densities, a simple quasi-deterministic approximation applies: the rate of coalescence within the front is λ∼2g(η)/(ρℓ), where ρ is the population density and ℓ=σ2/s is the characteristic scale of the wavefront; g(η) depends only on the strength of random drift, η=ρσs/2. The net effect of a sweep on coalescence also depends crucially on whether two lineages are ever both within the wavefront at the same time: even in the extreme case when coalescence within the front is instantaneous, the net rate of coalescence may be lower than in a single panmictic population. Sweeps can also have a substantial impact on the rate of gene flow. A single lineage will jump to a new location when it is hit by a sweep, with mean square displacement σeff(2)/σ(2)=(8/3)(L/ℓ)(Λ/R); this can be substantial if the species' range, L, is large, even if the species-wide rate of sweeps per map length, Λ/R, is small. This effect is half as strong in two dimensions. In contrast, the rate of coalescence between lineages, at random locations in space and on the genetic map, is proportional to (c/L)(Λ/R), where c is the wavespeed: thus, on average, one-dimensional structure is likely to reduce coalescence due to sweeps, relative to panmixis. In two dimensions, genes must move along the front before they can coalesce; this process is rapid, being dominated by rare fluctuations. This leads to a dramatically higher rate of coalescence within the wavefront than if lineages simply diffused along the front. Nevertheless, the net rate of coalescence due to a sweep through a two-dimensional population is likely to be lower than it would be with panmixis.

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.

Viral quasispecies evolution

Evolution of RNA viruses occurs through disequilibria of collections of closely related mutant spectra or mutant clouds termed viral quasispecies. Here we review the origin of the quasispecies concept and some biological implications of quasispecies dynamics. Two main aspects are addressed: (i) mutant clouds as reservoirs of phenotypic variants for virus adaptability and (ii) the internal interactions that are established within mutant spectra that render a virus ensemble the unit of selection. The understanding of viruses as quasispecies has led to new antiviral designs, such as lethal mutagenesis, whose aim is to drive viruses toward low fitness values with limited chances of fitness recovery. The impact of quasispecies for three salient human pathogens, human immunodeficiency virus and the hepatitis B and C viruses, is reviewed, with emphasis on antiviral treatment strategies. Finally, extensions of quasispecies to nonviral systems are briefly mentioned to emphasize the broad applicability of quasispecies theory.

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.

Intracellular kinetics of a growing virus: a genetically structured simulation for bacteriophage T7

Viruses have evolved to efficiently direct the resources of their hosts toward their own reproduction. A quantitative understanding of viral growth will help researchers develop antiviral strategies, design metabolic pathways, construct vectors for gene therapy, and engineer molecular systems that self-assemble. As a model system we examine here the growth of bacteriophage T7 in Escherichia coli using a chemical-kinetic framework. Data published over the last three decades on the genetics, physiology, and biophysics of phage T7 are incorporated into a genetically structured simulation that accounts for entry of the T7 genome into its host, expression of T7 genes, replication of T7 DNA, assembly of T7 procapsids, and packaging of T7 DNA to finally produce intact T7 progeny. Good agreement is found between the simulated behavior and experimental observations for the shift in transcription capacity from the host to the phage, the initiation times of phage protein synthesis, and the intracellular assembly of both wild-type phage and a fast-growing deletion mutant. The simulation is utilized to predict the effect of antisense molecules targeted to different T7 mRNA. Further, a postulated mechanism for the down regulation of T7 transcription in vivo is quantitatively examined and shown to agree with available data. The simulation is found to be a useful tool for exploring and understanding the dynamics of virus growth at the molecular level. (c) 1997 John Wiley & Sons, Inc. Biotechnol Bioeng 55: 375-389, 1997.

Metal recognition by in-vitro selection

We propose an in vitro selection strategy to identify bacteriophage variants that recognize metal ions in solution. In 6 M urea, phage T7 loses 99.9% of its activity in less than 5 min. Inactivation is accelerated by gold, but slowed by zinc and magnesium. Selection of phage over five generations in the presence of gold, zinc, and magnesium increases phage half-lives 4-, 10-, and 70-fold, respectively. As selections are repeated, phage become increasingly dependent on the specific metal employed in the selection, indicating the suitability of the strategy for optimization of metal-ion recognition.

Memory effects in a nonequilibrium growth model

We study memory effects in a kinetic roughening model. For d=1, a different dynamic scaling is uncovered in the memory dominated phases; the Kardar-Parisi-Zhang scaling is restored in the absence of noise. dc=2 represents the critical dimension where memory is shown to smoothen the roughening front (alpha<or=0). Studies on a discrete atomistic model in the same universality class reconfirm the analytical results in the large time limit, while a different scaling behavior shows up for t<tau, with tau being the memory characteristic of the atomistic model. Results can be generalized for other nonconservative systems.

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.

Spatial invasion by a mutant pathogen

Imagine a pathogen that is spreading radially as a circular wave through a population of susceptible hosts. In the interior of this circular region, the infection dies out due to a subcritical density of susceptibles. If a mutant pathogen, having some advantage over wild-type pathogens, arises in this region it is likely to die out without leaving a noticeable trace. Mutants that arise closer to the infection wavefront have access to more susceptible hosts and thus are more likely to become established and perhaps (locally) out-compete the original pathogen. Among the factors (position, transmission rate, pathogen-induced death rate) that influence the fate of a mutant, which are most important? What does this tell us about the types of mutants that are likely to invade and become established? How do such tendencies serve to steer the evolution of pathogens in a spatial setting? Do different types of models of the same phenomena lead to similar conclusions? We address these issues from the point of view of an individual-based stochastic spatial model of host-pathogen interactions. We consider the probability of a successful invasion by a single mutant as a function of the transmissibility and virulence strengths and the mutant position in the wavefront. Next, for a version of the model in which mutations arise spontaneously, we obtain analytical and simulation results on the mean time to a successful invasion. We also use our model predictions to gain insight into experimental data on bacteriophage plaques. Finally, we compare our results to those based on ordinary and partial differential equations to better understand how different models might influence our predictions on the fate of a mutant pathogen.

Quasispecies dynamics and RNA virus extinction

The extinction of foot-and-mouth disease virus (FMDV) is strongly influenced by mutation rates, types of mutations, relative viral fitness and virus population regimens during infection. Here we review experimental results and theoretical models that describe a contrast between the effective extinction of FMDV subjected to increased mutagenesis, and the remarkable resistance to extinction of the same and related FMDV clones subjected to serial bottleneck events. The results suggest procedures to master key parameters to develop effective antiviral strategies based on virus entry into error catastrophe.

Assessing the likelihood of recurrence during RNA evolution in vitro

Recurrence is the possibility of resulting in the same endpoint multiple times when a living system is allowed to evolve repeatedly starting from a given initial point. This concept is of concern to both evolutionary theoreticians and molecular biologists who use nucleic acid selection techniques to mimic biotic and computorial processes in the test tube. Using the continuous in vitro evolution methodology, many replicate experimental evolutionary lineages with populations of catalytic RNA were performed to gain insight into the parameters that could affect recurrence. The likelihood that the same genotype will result in parallel trials of an evolution experiment in vitro depends on several factors, including the phenotype under selection, the size and composition of the initial diverse pool of nucleic acids used in the experiment, the degree of mutation possible during the experiment, the shape of the fitness landscape through which the population evolves, and the strategies used to invoke selection and to search the landscape, among others. By considering these factors, it can be predicted that recurrence is more likely when a small, wild-type-based starting pool is used with efficient selection and search strategies involving little online mutagenesis within a rugged adaptive landscape with a strong local optimum. The recurrence experiments performed here on the 150-nucleotide ligase ribozyme demonstrate that it repeatedly jumps from one peak in a fitness landscape to another, apparently hurdling a deep fitness valley. These predictions can and should be tested by additional multiple replicates of actual evolution experiments in the laboratory.

Divergent evolution during an experimental adaptive radiation

How repeatable a process is evolution? Comparative studies of multicellular eukaryotes and experimental studies with unicellular prokaryotes document the repeated evolution of adaptive phenotypes during similar adaptive radiations, suggesting that the outcome of adaptive radiation is broadly reproducible. The goal of this study was to test this hypothesis by using phenotypic traits to infer the genetic basis of adaptation to simple carbon-limited environments in an extensive adaptive radiation. We used a clone of the bacterium Pseudomonas fluorescens to found two sets of experimental lines. The first set of lines was allowed to adapt to one of 23 novel environments for 1100 generations while the second set of lines was allowed to accumulate mutations by drift for 2000 generations. All lines were then assayed in the 95 environments provided by Biolog microplates to determine the phenotypic consequences of selection and drift. Replicate selection lines propagated in a common environment evolved similar adaptive components of their phenotype but showed extensive variation in non-adaptive phenotypic traits. This variation in non-adaptive phenotypic traits primarily resulted from the ascendance of different beneficial mutations in different lines. We argue that these results reconcile experimental and comparative approaches to studying adaptation by demonstrating that the convergent phenotypic evolution that occurs during adaptive radiation may be associated with radically different sets of beneficial mutations.

Propagation of viruses on micropatterned host cells

We have developed a technique to characterize the in vitro propagation of viruses. Microcontact printing was used to generate linear arrays of alkanethiols on gold surfaces, which served as substrates for the patterned culture of baby hamster kidney (BHK-21) cells. Vesicular stomatitis virus (VSV) was added to unpatterned cell reservoirs adjacent to the patterned cells and incubated, setting in motion a continuously advancing viral infection into the patterned cells. At different incubation times, multiple arrays were chemically fixed to stop the viral propagation. Viral propagation distances into the patterned cells were determined by indirect immunofluorescent labeling and visualization of the VSV surface glycoprotein (G). The infection spread at approximately 50 microm/h in the 140-microm lines. Moreover, different temporal stages of the infection process were simultaneously visualized along individual lines. These stages included initiation of infection, based on G protein expression; cell-cell fusion, based on virus-induced clustering of cell nuclei; and cytoskeletal degradation, based on localized release of cells from the surface. This work sets a foundation for parallel, high-throughput characterization of viral and cellular processes.

Macroevolution simulated with autonomously replicating computer programs

The process of adaptation occurs on two timescales. In the short term, natural selection merely sorts the variation already present in a population, whereas in the longer term genotypes quite different from any that were initially present evolve through the cumulation of new mutations. The first process is described by the mathematical theory of population genetics. However, this theory begins by defining a fixed set of genotypes and cannot provide a satisfactory analysis of the second process because it does not permit any genuinely new type to arise. The evolutionary outcome of selection acting on novel variation arising over long periods is therefore difficult to predict. The classical problem of this kind is whether 'replaying the tape of life' would invariably lead to the familiar organisms of the modern biota. Here we study the long-term behaviour of populations of autonomously replicating computer programs and find that the same type, introduced into the same simple environment, evolves on any given occasion along a unique trajectory towards one of many well-adapted end points.

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.

Growth pulsations in symmetric dendritic crystallization in thin polymer blend films

The crystallization of polymeric and metallic materials normally occurs under conditions far from equilibrium, leading to patterns that grow as propagating waves into the surrounding unstable fluid medium. The Mullins-Sekerka instability causes these wave fronts to break up into dendritic arms, and we anticipate that the normal modes of the dendrite tips have a significant influence on pattern growth. To check this possibility, we focus on the dendritic growth of polyethylene oxide in a thin-film geometry. This crystalline polymer is mixed with an amorphous polymer (polymethyl-methacrylate) to "tune" the morphology and clay was added to nucleate the crystallization. The tips of the main dendrite trunks pulsate during growth and the sidebranches, which grow orthogonally to the trunk, pulsate out of phase so that the tip dynamics is governed by a limit cycle. The pulsation period P increases sharply with decreasing film thickness L and then vanishes below a critical value L(c) approximately 80 nm. A change of dendrite morphology accompanies this transition.

An error-prone T7 RNA polymerase mutant generated by directed evolution

Viruses replicate their genomes at exceptionally high mutation rates. Their offspring evolve rapidly and therefore, are able to evade common immunological and chemical antiviral agents. In parallel, virus genomes cannot tolerate a further increase in mutation rate: Experimental evidence exists that even few additional mutations are sufficient for the extinction of a viral population. A future antiviral strategy might therefore aim at increasing the error-producing capacity of viral replication enzymes. We employed the principles of directed evolution and developed a scheme for the stringent positive selection of error-prone polymerase activity. A mutant T7 RNA polymerase with a nucleotide substitution error rate at least 20-fold greater than that of the wild-type was selected. This enzyme synthesized highly heterogeneous RNA products in vitro or in vivo and also decreased the replication efficiency of wild-type bacteriophage T7 during infection.

Bacteriophage K1-5 encodes two different tail fiber proteins, allowing it to infect and replicate on both K1 and K5 strains of Escherichia coli

A virulent double-stranded DNA bacteriophage, Phi K1-5, has been isolated and found to be capable of infecting Escherichia coli strains that possess either the K1 or the K5 polysaccharide capsule. Electron micrographs show that the virion consists of a small icosohedral head with short tail spikes, similar to members of the Podoviridae family. DNA sequence analysis of the region encoding the tail fiber protein showed two open reading frames encoding previously characterized hydrolytic phage tail fiber proteins. The first is the K5 lyase protein gene of Phi K5, which allows this phage to specifically infect K5 E. coli strains. A second open reading frame encodes a protein almost identical in amino acid sequence to the N-acetylneuraminidase (endosialidase) protein of Phi K1E, which allows this phage to specifically infect K1 strains of E. coli. We provide experimental evidence that mature phage particles contain both tail fiber proteins, and mutational analysis indicates that each protein can be independently inactivated. A comparison of the tail gene regions of Phi K5, Phi K1E, and Phi K1-5 shows that the genes are arranged in a modular or cassette configuration and suggests that this family of phages can broaden host range by horizontal gene transfer.

Big-benefit mutations in a bacteriophage inhibited with heat

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

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.

Exceptional convergent evolution in a virus

Replicate lineages of the bacteriophage phiX 174 adapted to growth at high temperature on either of two hosts exhibited high rates of identical, independent substitutions. Typically, a dozen or more substitutions accumulated in the 5.4-kilobase genome during propagation. Across the entire data set of nine lineages, 119 independent substitutions occurred at 68 nucleotide sites. Over half of these substitutions, accounting for one third of the sites, were identical with substitutions in other lineages. Some convergent substitutions were specific to the host used for phage propagation, but others occurred across both hosts. Continued adaptation of an evolved phage at high temperature, but on the other host, led to additional changes that included reversions of previous substitutions. Phylogenetic reconstruction using the complete genome sequence not only failed to recover the correct evolutionary history because of these convergent changes, but the true history was rejected as being a significantly inferior fit to the data. Replicate lineages subjected to similar environmental challenges showed similar rates of substitution and similar rates of fitness improvement across corresponding times of adaptation. Substitution rates and fitness improvements were higher during the initial period of adaptation than during a later period, except when the host was changed.

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.

Evidence for multiple adaptive peaks from populations of bacteria evolving in a structured habitat

Natural selection tends to promote the divergence of populations living in different environments. Even in identical environments, however, replicate populations may diverge if they find alternative adaptive solutions. We describe the evolution of 18 bacterial populations (Comamonas sp.) founded from a single progenitor genotype and propagated separately for 1000 generations in two distinct environments, one physically unstructured (mass-action liquid) and the other structured (agar surfaces). Phenotypic diversity, as reflected in colony morphology, was greater in the structured habitat than in the unstructured habitat. More importantly, the trajectories for mean fitness, as measured by competition against the common ancestor, were more divergent for populations in the structured habitat than those in the unstructured habitat. Structured environments may accelerate evolutionary diversification by promoting genetic polymorphisms within populations, thereby increasing the complexity of genetic constraints that allow divergence among replicate populations.