Dependence of epistasis on environment and mutation severity as revealed by in silico mutagenesis of phage t7

Understanding how interactions among deleterious mutations affect fitness may shed light on a variety of fundamental biological phenomena, including the evolution of sex, the buffering of genetic variations, and the topography of fitness landscapes. It remains an open question under what conditions and to what extent such interactions may be synergistic or antagonistic. To address this question, we employed a computer model for the intracellular growth of bacteriophage T7. We created in silico 90,000 mutants of phage T7, each carrying from 1 to 30 mutations, and evaluated the fitness of each by simulating its growth cycle. The simulations sought to account for the severity of single deleterious mutations on T7 growth, as well as the effect of the resource environment on our fitness measures. We found that mildly deleterious mutations interacted synergistically in poor-resource environments but antagonistically in rich-resource environments. However, severely deleterious mutations always interacted antagonistically, irrespective of environment. These results suggest that synergistic epistasis may be difficult to experimentally distinguish from nonepistasis because its effects appear to be most pronounced when the effects of mutations on fitness are most challenging to measure. Our approach demonstrates how computer simulations of developmental processes can be used to quantitatively study genetic interactions at the population level.     

Effects of Escherichia coli physiology on growth of phage T7 in vivo and in silico

Phage development depends not only upon phage functions but also on the physiological state of the host, characterized by levels and activities of host cellular functions. We established Escherichia coli at different physiological states by continuous culture under different dilution rates and then measured its production of phage T7 during a single cycle of infection. We found that the intracellular eclipse time decreased and the rise rate increased as the growth rate of the host increased. To develop mechanistic insight, we extended a computer simulation for the growth of phage T7 to account for the physiology of its host. Literature data were used to establish mathematical correlations between host resources and the host growth rate; host resources included the amount of genomic DNA, pool sizes and elongation rates of RNA polymerases and ribosomes, pool sizes of amino acids and nucleoside triphosphates, and the cell volume. The in silico (simulated) dependence of the phage intracellular rise rate on the host growth rate gave quantitatively good agreement with our in vivo results, increasing fivefold for a 2.4-fold increase in host doublings per hour, and the simulated dependence of eclipse time on growth rate agreed qualitatively, deviating by a fixed delay. When the simulation was used to numerically uncouple host resources from the host growth rate, phage growth was found to be most sensitive to the host translation machinery, specifically, the level and elongation rate of the ribosomes. Finally, the simulation was used to follow how bottlenecks to phage growth shift in response to variations in host or phage functions.     

In Vitro Packaging of UV Radiation-Damaged DNA from Bacteriophage T7

When DNA from bacteriophage T7 is irradiated with UV light, the efficiency with which this DNA can be packaged in vitro to form viable phage particles is reduced. A comparison between irradiated DNA packaged in vitro and irradiated intact phage particles shows almost identical survival as a function of UV dose when Escherichia coli wild type or polA or uvrA mutants are used as the host. Although uvrA mutants perform less host cell reactivation, the polA strains are identical with wild type in their ability to support the growth of irradiated T7 phage or irradiated T7 DNA packaged in vitro into complete phage. An examination of in vitro repair performed by extracts of T7-infected E.coli suggests that T7 DNA polymerase may substitute for E. coli DNA polymerase I in the resynthesis step of excision repair. Also tested was the ability of a similar in vitro repair system that used extracts from uninfected cells to restore biological activity of irradiated DNA. When T7 DNA damaged by UV irradiation was treated with an endonuclease from Micrococcus luteus that is specific for pyrimidine dimers and then was incubated with an extract of uninfected E. coli capable of removing pyrimidine dimers and restoring the DNA of its original (whole genome size) molecular weight, this DNA showed a higher packaging efficiency than untreated DNA, thus demonstrating that the in vitro repair system partially restored the biological activity of UV-damaged DNA.     

Compensatory evolution in response to a novel RNA polymerase: orthologous replacement of a central network gene

A bacteriophage genome was forced to evolve a new system of regulation by replacing its RNA polymerase (RNAP) gene, a central component of the phage developmental pathway, with that of a relative. The experiment used the obligate lytic phage T7 and the RNAP gene of phage T3. T7 RNAP uses 17 phage promoters, which are responsible for all middle and late gene expression, DNA replication, and progeny maturation, but the enzyme has known physical contacts with only 2 other phage proteins. T3 RNAP was supplied in trans by the bacterial host to a T7 genome lacking its own RNAP gene and the phage population was continually propagated on naive bacteria throughout the adaptation. Evolution of the T3 RNAP gene was thereby prevented, and selection was for the evolution of regulatory signals throughout the phage genome. T3 RNAP transcribes from T7 promoters only at low levels, but a single mutation in the promoter confers high expression, providing a ready mechanism for reevolution of gene expression in this system. When selected for rapid growth, fitness of the engineered phage evolved from a low of 5 doublings/h to 33 doublings/h, close to the expected maximum of 37 doublings/h. However, the experiment was terminated before it could be determined accurately that fitness had reached an obvious plateau, and it is not known whether further adaptation could have resulted in complete recovery of fitness. More than 30 mutations were observed in the evolved genome, but changes were found in only 9 of the 16 promoters, and several coding changes occurred in genes with no known contacts with the RNAP. Surprisingly, the T7 genome adapted to T3 RNAP also maintained high fitness when using T7 RNAP, suggesting that the extreme incompatibility of T7 elements with T3 RNAP is not an invariant property of divergence in these expression systems.     

Inadequate inhibition of host RNA polymerase restricts T7 bacteriophage growth on hosts overexpressing udk

Overexpression of udk, an Escherichia coli gene encoding a uridine/cytidine kinase, interferes with T7 bacteriophage growth. We show here that inhibition of T7 phage growth by udk overexpression can be overcome by inhibition of host RNA polymerase. Overexpression of gene 2, whose product inhibits host RNA polymerase, restores T7 phage growth on hosts overexpressing udk. In addition, rifampicin, an inhibitor of host RNA polymerase, restores the burst size of T7 phage on udk-overexpressing hosts to normal. In agreement with these findings, suppressor mutants that overcome the inhibition arising from udk overexpression gain the ability to grow on hosts that are resistant to inhibition of RNA polymerase by gene 2 protein, and suppressor mutants that overcome a lack of gene 2 protein gain the ability to grow on hosts that overexpress udk. Mutations that eliminate or weaken strong promoters for host RNA polymerase in T7 DNA, and mutations in T7 gene 3.5 that affect its interaction with T7 RNA polymerase, also reduce the interference with T7 growth by host RNA polymerase. We propose a general model for the requirement of host RNA polymerase inhibition.     

Energy-efficient growth of phage Q Beta in Escherichia coli

The role of natural selection in the optimal design of organisms is controversial. Optimal forms, functions, or behaviors of organisms have long been claimed without knowledge of how genotype contributes to phenotype, delineation of design constraints, or reference to alternative designs. Moreover, arguments for optimal designs have been often based on models that were difficult, if not impossible, to test. Here, we begin to address these issues by developing and probing a kinetic model for the intracellular growth of bacteriophage Q beta in Escherichia coli. The model accounts for the energetic costs of all template-dependent polymerization reactions, in ATP equivalents, including RNA-dependent RNA elongation by the phage replicase and synthesis of all phage proteins by the translation machinery of the E. coli host cell. We found that translation dominated phage growth, requiring 85% of the total energy expenditure. Only 10% of the total energy was applied to activities other than the direct synthesis of progeny phage components, reflecting primarily the cost of making the negative-strand RNA template that is needed for replication of phage genomic RNA. Further, we defined an energy efficiency of phage growth and showed its direct relationship to the yield of phage progeny. Finally, we performed a sensitivity analysis and found that the growth of wild-type phage was optimized for progeny yield or energy efficiency, suggesting that phage Q beta has evolved to optimally utilize the finite resources of its host cells.     

Predicting evolution from genomics: experimental evolution of bacteriophage T7

A wealth of molecular biology has been exploited in designing and interpreting experimental evolution studies with bacteriophage T7. The modest size of its genome (40 kb dsDNA) and the ease of making genetic constructs, combined with the many genetic resources for its host (Escherichia coli), have enabled comprehensive and detailed studies of experimental adaptations. In several studies, the genome was specifically altered (gene knockouts, gene replacements, reordering of genetic elements) such that a priori knowledge of genetics and biochemistry of the phage could be used to predict the pathways of compensatory evolution when the modified phage is adapted to recover fitness. In other work, the phage has been adapted to specific environmental conditions chosen to select phenotypic outcomes with a quantitative basis, and the molecular bases of that evolution have been explored. Predicting the outcomes of these adaptations has been challenging. In hindsight, one-third to one-half of the compensatory nucleotide changes observed during the adaptation can be rationalized based on T7 biology. This rationalization usually only applies at the genetic level-a gene product may be known to be involved in the affected pathway, but it usually remains unknown how the observed change affects activity. The progress is encouraging, but the prediction of experimental evolution pathways remains far from complete, and is still sometimes confounded by observation when an adaptation yields a completely unexpected outcome.     

Virulence reaction norms across a food gradient

Host-parasite interactions involve competition for nutritional resources between hosts and the parasites growing within them. Consuming part of a host's resources is one cause of a parasite's virulence, i.e. part of the fitness cost imposed on the host by the parasite. The influence of a host's nutritional conditions on the virulence of a parasite was experimentally tested using the mosquito Aedes aegypti and the microsporidian parasite Vavraia culicis. A condition-dependent expression of virulence was found and a positive relation between virulence and transmissibility was established. Spore production was positively influenced by host food availability, indicating that the parasite's within-host growth is limited by host condition. We also investigated how the fitness of each partner varied across the nutritional gradient and demonstrated that the sign of the correlation between host fitness and parasite fitness depended on the amount of nutritional resources available to the host.     

Evolutionary principles and synthetic biology: avoiding a molecular tragedy of the commons with an engineered phage

ABSTRACT: BACKGROUND: In prior work, adding a gene to phage T7 that degraded the host K1 capsule facilitated growth when plated on capsulated hosts. However, the transgenic protein (an endosialidase) is expressed as an exoenzyme, released from the cell at lysis but unattached to the phage particle. There is thus the possibility that the gene will be subject to a tragedy of the commons and be selected against, if the enzyme benefits other genomes. Results: This evolutionary perspective was supported in short term experiments. The genome carrying the endosialidase gene was favored on a capsulated host if grown in physical isolation of control genomes (lacking the gene) but was selected against otherwise. Conclusions: These results challenge efforts to engineer phages with exoenzymes that degrade biofilm polymers. If biofilms do not facilitate spatially structured phage growth, the transgenic enzymes may be rapidly eliminated from the phage population after release in the environment, even if the transgene benefits overall phage growth on the biofilm.     

Phenotypic Correlation of Juvenile Growth Rate between Different Consecutive foraging Environments in a Salmonid fish: A Field Experiment

Many organisms occupy considerably different environments during individual's lifespan. We are interested in how the phenotypic characteristics that are favourable in the earlier environment predict fitness in the later environment. High predictability of fitness between the two consecutive environments suggests that the either the same traits are favored in both environments, or that the favourable traits are genetically correlated. In this study, we ask how similarity of consecutive foraging environments affects the phenotypic correlation of juvenile brown trout growth rate. More specifically, we used a genetically narrow stock of hatchery-bred fish to contrast individual growth rates between high and low density hatchery environments, and thereafter between hatchery and natural lake environment. As expected, growth rate was highly dependent on the environment, and the fish showed considerable phenotypic plasticity. Furthermore, we found a strong positive correlation in growth rate between similar foraging environments, for example, between high and low density hatchery stocks, and between hatchery and a lake with small fish as main prey. However, hatchery growth did not predict growth rate in lakes where fish had to forage on bottom-dwelling invertebrates. Our results suggest that when the consecutive environments differed dramatically with respect to traits that fish use for foraging, relative performance of individual fish changed, earlier performance not being an accurate predictor of performance in the new environment. In this case, fitness of the fish was determined by an environment-specific set of traits that were not the same between the two consecutive environments. The result indicates that assessment of individual performance may be highly environment specific in trout. This revised version was published online in July 2006 with corrections to the Cover Date.     

Synonymous codon preferences in bacteriophage T4: a distinctive use of transfer RNAs from T4 and from its host Escherichia coli.

Codon usage data of bacteriophage T4 genes were compiled and synonymous codon preferences were investigated in comparison with tRNA availabilities in an infected cell. Since the genome of T4 is highly AT rich and its codon usage pattern is significantly different from that of its host Escherichia coli, certain codons of T4 genes need to be translated by appropriate host transfer RNAs present in minor amounts. To avoid this predicament, T4 phage seems to direct the synthesis of its own tRNA molecules and these phage tRNAs are suggested to supplement the host tRNA population with isoacceptors that are normally present in minor amounts. A positive correlation was found in that the frequency of E. coli optimal codons in T4 genes increases as the number of protein monomers per phage particle increases. A negative correlation was also found between the number of protein monomers per phage and the frequency of "T4 optimal codons", which are defined as those codons that are efficiently recognized by T4 tRNAs. From these observations it was proposed that tRNAs from the host are predominantly used for translation of highly expressed T4 genes while tRNAs from T4 tend to be used for translation of weakly expressed T4 genes. This distinctive tRNA-usage in T4 may be an optimization of translational efficiency, and an adjustment of T4-encoded tRNAs to the synonymous codon preferences, which are largely influenced by the high genomic AT-content, would have occurred during evolution.     

Within-host competition determines reproductive success of temperate bacteriophages

Within-host competition between parasites is frequently invoked as a major force for parasite evolution, yet quantitative studies on its extent in an organismal group are lacking. Temperate bacteriophages are diverse and abundant parasites of bacteria, distinguished by their ability to enter a facultative dormant state in their host. Bacteria can accumulate multiple phages that may eventually abandon dormancy in response to host stress. Host resources are then converted into phage particles, whose release requires cell death. To study within-host competition between phages, I used the bacterium Escherichia coli and 11 lambdoid phages to construct single and double lysogens. Lysogenic bacterial cultures were then induced and time to host cell lysis and productivity of phages was measured. In double lysogens, this revealed strong competitive interactions as in all cases productivity of at least one phage declined. The outcome of within-host competition was often asymmetrical, and phages were found to vary hierarchically in within-host competitive ability. In double infections, the phage with the shorter lysis time determined the timing of cell lysis, which was associated with a competitive advantage when time differences were large. The results emphasize that within-host competition greatly affects phage fitness and that multiple infections should be considered an integral part of bacteriophage ecology.     

Inhibition of gene expression of T7-related phages by prophage P1

Summary The gene expression of nine phages of the T7 group was compared after infection of Escherichia coli B(P1). With the exception of phage 13a which grew normally, all of them infected E. coli B(P1) abortively. Differences were found in the efficiency of host killing which ranged from 100% for phage 13a to 37% for phage A1122. Infection by T7 prevented colony formation by about 70% of the cells but they showed filamentous growth until about 2h after infection. It was shown by SDS-polyacrylamide gel electrophoresis and autoradiography of [³⁵S]methionine-labelled phage-coded proteins that all phages except for 13a showed measurable expression only of the early genes. No correlation was observed between killing capacity and the pattern of gene expression, and the ability to hydrolyse S-adenosyl-methionine (SAM, a cofactor for the P1 restriction endonuclease) by means of a phage-coded SAMase. Mixed infection of E. coli B(P1) with 13a and T7 yielded mixed progeny indistinguishable from that observed after mixed infection of the normal host E. coli B. Genetic crosses with amber mutants of 13a and T7 showed that the 13a marker opo ⁺ (overcomes P one), required for growth on B(P1), is located in the early region, to the left of gene 1 (RNA polymerase gene).     

重组T4噬菌体gp37基因可扩大其在沙门氏菌中的宿主范围

[目的]将T4噬菌体WG01宿主决定区的gp37基因片段,与另一株T4噬菌体QL01的相应基因进行同源重组,从而获得嵌合噬菌体并进行宿主谱分析,为阐明T4噬菌体的宿主谱形成机制以及快速筛选针对特定病原菌的噬菌体奠定了基础。[方法]通过同源重组的方法将WG01 gp37上的8个基因片段分别替换给QL01,用沙门氏菌作为宿主菌筛选嵌合噬菌体,并对嵌合噬菌体进行宿主谱、最佳感染复数、一步生长曲线和遗传稳定性测定。[结果]本研究共获得了5株嵌合噬菌体(QWA、QWC、QWF、QWG、QWFG)。宿主谱试验结果表明,与噬菌体QL01相比,嵌合噬菌体对21株沙门宿主菌分别可以多裂解7、8、4、10和9株菌,即嵌合噬菌体都获得了相对较宽的宿主谱,其中QWG的沙门氏菌宿主菌拓宽最多。生物学特性试验结果表明,嵌合噬菌体QWG生物学特性稳定。嵌合噬菌体QWG经连续传代培养20代,测序分析第1代和第20代嵌合噬菌体尾丝蛋白基因在传代过程中的稳定性,测序结果表明,嵌合噬菌体改造部分的基因能稳定遗传。[结论]用基因改造的方法可以产生宿主谱拓宽且能稳定遗传的嵌合噬菌体,为快速筛选针对特定病原菌的噬菌体提供了可能。     

The role of ATP and membrane potential in the penetration of phage T17 DNA into the cell during infection

The role of ATP and membrane potential in phage T7 DNA injection into E. coli during infection has been studied. Entrance of phage T7 genes of class II and III was shown to be prevented by arsenate, indicating the requirement for phosphorylated macroergs in the phage DNA injection. The injection process was also inhibited by exposition of the cells to the uncoupler of oxidative phosphorylation. Dependence of the injection efficiency on the membrane-potential value has been shown to be sigmoidal, which suggests a regulatory role of the membrane potential in phage T7 DNA injection from the virion into the host cell.     

Biological consequences of infection of Escherichia coli B by alkylated T7 bacteriophage.

Alkylation of T7 bacteriophage considerably delayed phage development and reduced the phage's killing action on host cells. Only a small fraction of infected cells produced phage. For these phages, the latent period was markedly prolonged but the burst was equivalent to or only slightly lower than that of untreated phage. In the progeny of alkylated phage, there was an increase in the fraction of defective particles as well as a change in their morphology. These data show that infection with alkylated T7 bacteriophage is to a large degree abortive; hence, biological consequences of this infection are very different from those characteristic of a normal virus infection.     

Ongoing phenotypic and genomic changes in experimental coevolution of RNA bacteriophage Qβ and Escherichia coli

According to the Red Queen hypothesis or arms race dynamics, coevolution drives continuous adaptation and counter-adaptation. Experimental models under simplified environments consisting of bacteria and bacteriophages have been used to analyze the ongoing process of coevolution, but the analysis of both parasites and their hosts in ongoing adaptation and counter-adaptation remained to be performed at the levels of population dynamics and molecular evolution to understand how the phenotypes and genotypes of coevolving parasite-host pairs change through the arms race. Copropagation experiments with Escherichia coli and the lytic RNA bacteriophage Qβ in a spatially unstructured environment revealed coexistence for 54 days (equivalent to 163-165 replication generations of Qβ) and fitness analysis indicated that they were in an arms race. E. coli first adapted by developing partial resistance to infection and later increasing specific growth rate. The phage counter-adapted by improving release efficiency with a change in host specificity and decrease in virulence. Whole-genome analysis indicated that the phage accumulated 7.5 mutations, mainly in the A2 gene, 3.4-fold faster than in Qβ propagated alone. E. coli showed fixation of two mutations (in traQ and csdA) faster than in sole E. coli experimental evolution. These observations suggest that the virus and its host can coexist in an evolutionary arms race, despite a difference in genome mutability (i.e., mutations per genome per replication) of approximately one to three orders of magnitude.