Layers of evolvability in a bacteriophage life history trait

Implication of a gene deletion on a Salmonella Enteritidis phage growth parameters

Synthetic biology has been applied countless times for the modification and improvement of bacterial strains and for the synthesis of products that do not exist in nature. Phages are natural predators of bacteria controlling their population levels; however, their genomes carry several genes with unknown functions. In this work, Bacteriophage Recombineering of Electroporated DNA was used to assess the influence of deletion of a single gene with unknown function in the overall replication parameters of Salmonella phage PVP-SE2. Deletion of ORF_01, transcribed immediately after infection, reduced both the latent and rise periods by 5 min in PVP-SE2ΔORF_01 compared to the wild-type phage. A direct consequence of the deletion led to a smaller progeny release per infected cell by the mutant compared to the wild-type phage. Despite the difference in growth characteristics, the mutant phage remained infective towards exponentially growing cells. The mutation engineered endured for at least ten passages, showing that there is no reversion back to the wild-type sequence. This study provides proof of concept that methodologies used for phage engineering should always be complemented by phage growth characterization to assess whether a mutation can trigger undesirable characteristics.

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.

Efficient Genome Engineering of a Virulent Klebsiella Bacteriophage Using CRISPR-Cas9

Klebsiella pneumoniae is one of the most common nosocomial opportunistic pathogens and usually exhibits multiple-drug resistance. Phage therapy, a potential therapeutic to replace or supplement antibiotics, has attracted much attention. However, very few Klebsiella phages have been well characterized because of the lack of efficient genome-editing tools. Here, Cas9 from Streptococcus pyogenes and a single guide RNA (sgRNA) were used to modify a virulent Klebsiella bacteriophage, phiKpS2. We first evaluated the distribution of sgRNA activity in phages and proved that it is largely inconsistent with the predicted activity from current models trained on eukaryotic cell data sets. A simple CRISPR-based phage genome-editing procedure was developed based on the discovery that homologous arms as short as 30 to 60 bp were sufficient to introduce point mutation, gene deletion, and swap. We also demonstrated that weak sgRNAs could be used for precise phage genome editing but failed to select random recombinants, possibly because inefficient cleavage can be tolerated through continuous repair by homologous recombination with the uncut genomes. Small frameshift deletion was proved to be an efficient way to evaluate the essentiality of phage genes. By using the abovementioned strategies, a putative promoter and nine genes of phiKpS2 were successfully deleted. Interestingly, the holin gene can be deleted with little effect on phiKpS2 infection, but the reason is not yet clear. This study established an efficient, time-saving, and cost-effective procedure for phage genome editing, which is expected to significantly promote the development of bacteriophage therapy.IMPORTANCE In the present study, we have addressed efficient, time-saving, and cost-effective CRISPR-based phage genome editing of Klebsiella phage, which has the potential to significantly expand our knowledge of phage-host interactions and to promote applications of phage therapy. The distribution of sgRNA activity was first evaluated in phages. Short homologous arms were proven to be enough to introduce point mutation, small frameshift deletion, gene deletion, and swap into phages, and weak sgRNAs were proven useful for precise phage genome editing but failed to select random recombinants, all of which makes the CRISPR-based phage genome-editing method easier to use.

Experimental evolution of UV resistance in a phage

The dsDNA bacteriophage T7 was subjected to 30 cycles of lethal ultraviolet light (UV) exposure to select increased resistance to UV. The exposure effected a 0.9999 kill of the ancestral population, and survival of the ending population was nearly 50-fold improved. At the end point, a 2.1 kb deletion of early genes and three substitutions in structural-genes were the only changes observed at high frequency throughout the 40 kb genome; no changes were observed in genes affecting DNA metabolism. The deletion accounted for only a two-fold improvement in survival. One possible explanation of its benefit is that it represents an error catastrophe, whereby the genome experiences a reduced mutation rate. The mechanism of benefit provided by the three structural-gene mutations remains unknown. The results offer some hope of artificially evolving greater protection against sunlight damage in applications of phage therapy to plants, but the response of T7 is weak compared to that observed in bacteria selected to resist ionizing radiation. Because of the weak response, mathematical analysis of the selection process was performed to determine how the protocol might have been modified to achieve a greater response, but the greatest protection may well come from evolving phages to bind materials that block the UV.

Experimental evolution reveals hidden diversity in evolutionary pathways

Replicate populations of natural and experimental organisms often show evidence of parallel genetic evolution, but the causes are unclear. The wrinkly spreader morph of Pseudomonas fluorescens arises repeatedly during experimental evolution. The mutational causes reside exclusively within three pathways. By eliminating these, 13 new mutational pathways were discovered with the newly arising WS types having fitnesses similar to those arising from the commonly passaged routes. Our findings show that parallel genetic evolution is strongly biased by constraints and we reveal the genetic bases. From such knowledge, and in instances where new phenotypes arise via gene activation, we suggest a set of principles: evolution proceeds firstly via pathways subject to negative regulation, then via promoter mutations and gene fusions, and finally via activation by intragenic gain-of-function mutations. These principles inform evolutionary forecasting and have relevance to interpreting the diverse array of mutations associated with clinically identical instances of disease in humans.

DOI:http://dx.doi.org/10.7554/eLife.07074.001

Different living things often develop similar strategies to adapt to the environments in which they live. Sometimes two species that share a common ancestor independently evolve the same trait by changing the exact same genes. This is called ‘parallel evolution’, and it has led some scientists to ask: are there certain traits that can only evolve in a limited number of ways? Or are there other ways to evolve the same trait that, for some reason, are not explored?

Experimentally, investigating these questions is challenging, but parallel evolution occurs in the laboratory as well as in the wild. Many commonly studied organisms—such as fruit flies or bacteria—can be used in relevant studies, because they can be grown in large numbers and then exposed to identical environments. However, if this method fails to find a new way that a trait can evolve, it doesn't mean that alternative mechanisms do not exist.

Lind et al. used a different approach that instead relies on removing all of the known pathways that can be mutated to produce a given trait and then seeing if that trait can still evolve via mutations elsewhere. The experiments involved a bacterium called Pseudomonas fluorescens that can evolve to grow flattened and wrinkled colonies (instead of smooth, round ones) when it has to compete for access to oxygen.

Previous experiments had shown that the evolution of the so-called ‘wrinkly spreader’ form can be caused by mutations in one of three biological pathways. But P. fluorescens can survive unharmed without these pathways, which enabled Lind et al. to ask if there might be other ways that this trait could evolve. Bacteria without these three pathways were engineered and then grown under oxygen-deprived conditions. This experiment produced 91 new mutants that each had the wrinkly spreader phenotype. Further experiments revealed that together these mutants represented 13 previously unrecognized ways that the ‘wrinkly spreader’ phenotype can evolve.

The new rare mutants had similar fitness as the previously known, common ones—so this cannot explain why they hadn't been seen before. Lind et al. instead suggest a set of principles to explain why these newly discovered pathways are rarely mutated and how genetic constraints can bias the outcome of evolution. Further work could investigate whether these principles can help us to predict the course of evolution in other biological contexts, such as in the evolution of antibiotic resistance.

DOI:http://dx.doi.org/10.7554/eLife.07074.002

Bacteriophages use an expanded genetic code on evolutionary paths to higher fitness

Bioengineering advances have made it possible to fundamentally alter the genetic codes of organisms. However, the evolutionary consequences of expanding an organism's genetic code with a non-canonical amino acid are poorly understood. Here we show that bacteriophages evolved on a host that incorporates 3-iodotyrosine at the amber stop codon acquired neutral and beneficial mutations to this new amino acid in their proteins, demonstrating that an expanded genetic code increases evolvability.

Lysis delay and burst shrinkage of coliphage T7 by deletion of terminator Tφ reversed by deletion of early genes

Bacteriophage T7 terminator Tϕ is a class I intrinsic terminator coding for an RNA hairpin structure immediately followed by oligo(U), which has been extensively studied in terms of its transcription termination mechanism, but little is known about its physiological or regulatory functions. In this study, using a T7 mutant phage, where a 31-bp segment of Tϕ was deleted from the genome, we discovered that deletion of Tϕ from T7 reduces the phage burst size but delays lysis timing, both of which are disadvantageous for the phage. The burst downsizing could directly result from Tϕ deletion-caused upregulation of gene 17.5, coding for holin, among other Tϕ downstream genes, because infection of gp17.5-overproducing Escherichia coli by wild-type T7 phage showed similar burst downsizing. However, the lysis delay was not associated with cellular levels of holin or lysozyme or with rates of phage adsorption. Instead, when allowed to evolve spontaneously in five independent adaptation experiments, the Tϕ-lacking mutant phage, after 27 or 29 passages, recovered both burst size and lysis time reproducibly by deleting early genes 0.5, 0.6, and 0.7 of class I, among other mutations. Deletion of genes 0.5 to 0.7 from the Tϕ-lacking mutant phage decreased expression of several Tϕ downstream genes to levels similar to that of the wild-type phage. Accordingly, phage T7 lysis timing is associated with cellular levels of Tϕ downstream gene products. This suggests the involvement of unknown factor(s) besides the known lysis proteins, lysozyme and holin, and that Tϕ plays a role of optimizing burst size and lysis time during T7 infection. IMPORTANCE Bacteriophages are bacterium-infecting viruses. After producing numerous progenies inside bacteria, phages lyse bacteria using their lysis protein(s) to get out and start a new infection cycle. Normally, lysis is tightly controlled to ensure phage progenies are maximally produced and released at an optimal time. Here, we have discovered that phage T7, besides employing its known lysis proteins, additionally uses its transcription terminator Tϕ to guarantee the optimal lysis of the E. coli host. Tϕ, positioned in the middle of the T7 genome, must be inactivated at least partially to allow for transcription-driven translocation of T7 DNA into hosts and expression of Tϕ downstream but promoter-lacking genes. What role is played by Tϕ before inactivation? Without Tϕ, not only was lysis time delayed but also the number of progenies was reduced in this study. Furthermore, T7 can overcome Tϕ deletion by further deleting some genes, highlighting that a phage has multiple strategies for optimizing lysis.

A common, non-optimal phenotypic endpoint in experimental adaptations of bacteriophage lysis time

Background

Optimality models of evolution, which ignore genetic details and focus on natural selection, are widely used but sometimes criticized as oversimplifications. Their utility for quantitatively predicting phenotypic evolution can be tested experimentally. One such model predicts optimal bacteriophage lysis interval, how long a virus should produce progeny before lysing its host bacterium to release them. The genetic basis of this life history trait is well studied in many easily propagated phages, making it possible to test the model across a variety of environments and taxa.

Results

We adapted two related small single-stranded DNA phages, ΦX174 and ST-1, to various conditions. The model predicted the evolution of the lysis interval in response to host density and other environmental factors. In all cases the initial phages lysed later than predicted. The ΦX174 lysis interval did not evolve detectably when the phage was adapted to normal hosts, indicating complete failure of optimality predictions. ΦX174 grown on slyD-defective hosts which initially entirely prevented lysis readily recovered to a lysis interval similar to that attained on normal hosts. Finally, the lysis interval still evolved to the same endpoint when the environment was altered to delay optimal lysis interval. ST-1 lysis interval evolved to be ~2 min shorter, qualitatively in accord with predictions. However, there were no changes in the single known lysis gene. Part of ST-1's total lysis time evolution consisted of an earlier start to progeny production, an unpredicted phenotypic response outside the boundaries of the optimality model.

Conclusions

The consistent failure of the optimality model suggests that constraint and genetic details affect quantitative and even qualitative success of optimality predictions. Several features of ST-1 adaptation show that lysis time is best understood as an output of multiple traits, rather than in isolation.

Functional analysis of the holin-like proteins of mycobacteriophage Ms6

The mycobacteriophage Ms6 is a temperate double-stranded DNA (dsDNA) bacteriophage which, in addition to the predicted endolysin (LysA)-holin (Gp4) lysis system, encodes three additional proteins within its lysis module: Gp1, LysB, and Gp5. Ms6 Gp4 was previously described as a class II holin-like protein. By analysis of the amino acid sequence of Gp4, an N-terminal signal-arrest-release (SAR) domain was identified, followed by a typical transmembrane domain (TMD), features which have previously been observed for pinholins. A second putative holin gene (gp5) encoding a protein with a predicted single TMD at the N-terminal region was identified at the end of the Ms6 lytic operon. Neither the putative class II holin nor the single TMD polypeptide could trigger lysis in pairwise combinations with the endolysin LysA in Escherichia coli. One-step growth curves and single-burst-size experiments of different Ms6 derivatives with deletions in different regions of the lysis operon demonstrated that the gene products of gp4 and gp5, although nonessential for phage viability, appear to play a role in controlling the timing of lysis: an Ms6 mutant with a deletion of gp4 (Ms6(Δgp4)) caused slightly accelerated lysis, whereas an Ms6(Δgp5) deletion mutant delayed lysis, which is consistent with holin function. Additionally, cross-linking experiments showed that Ms6 Gp4 and Gp5 oligomerize and that both proteins interact. Our results suggest that in Ms6 infection, the correct and programmed timing of lysis is achieved by the combined action of Gp4 and Gp5.

Optimality models in the age of experimental evolution and genomics

Optimality models have been used to predict evolution of many properties of organisms. They typically neglect genetic details, whether by necessity or design. This omission is a common source of criticism, and although this limitation of optimality is widely acknowledged, it has mostly been defended rather than evaluated for its impact. Experimental adaptation of model organisms provides a new arena for testing optimality models and for simultaneously integrating genetics. First, an experimental context with a well-researched organism allows dissection of the evolutionary process to identify causes of model failure--whether the model is wrong about genetics or selection. Second, optimality models provide a meaningful context for the process and mechanics of evolution, and thus may be used to elicit realistic genetic bases of adaptation--an especially useful augmentation to well-researched genetic systems. A few studies of microbes have begun to pioneer this new direction. Incompatibility between the assumed and actual genetics has been demonstrated to be the cause of model failure in some cases. More interestingly, evolution at the phenotypic level has sometimes matched prediction even though the adaptive mutations defy mechanisms established by decades of classic genetic studies. Integration of experimental evolutionary tests with genetics heralds a new wave for optimality models and their extensions that does not merely emphasize the forces driving evolution.

Analysis of high-throughput sequencing and annotation strategies for phage genomes

Background

Bacterial viruses (phages) play a critical role in shaping microbial populations as they influence both host mortality and horizontal gene transfer. As such, they have a significant impact on local and global ecosystem function and human health. Despite their importance, little is known about the genomic diversity harbored in phages, as methods to capture complete phage genomes have been hampered by the lack of knowledge about the target genomes, and difficulties in generating sufficient quantities of genomic DNA for sequencing. Of the approximately 550 phage genomes currently available in the public domain, fewer than 5% are marine phage.

Methodology/Principal Findings

To advance the study of phage biology through comparative genomic approaches we used marine cyanophage as a model system. We compared DNA preparation methodologies (DNA extraction directly from either phage lysates or CsCl purified phage particles), and sequencing strategies that utilize either Sanger sequencing of a linker amplification shotgun library (LASL) or of a whole genome shotgun library (WGSL), or 454 pyrosequencing methods. We demonstrate that genomic DNA sample preparation directly from a phage lysate, combined with 454 pyrosequencing, is best suited for phage genome sequencing at scale, as this method is capable of capturing complete continuous genomes with high accuracy. In addition, we describe an automated annotation informatics pipeline that delivers high-quality annotation and yields few false positives and negatives in ORF calling.

Conclusions/Significance

These DNA preparation, sequencing and annotation strategies enable a high-throughput approach to the burgeoning field of phage genomics.

Tuning a genetic switch: experimental evolution and natural variation of prophage induction

Genetic switches allow organisms to modulate their phenotype in response to environmental changes. Understanding the evolutionary processes by which switches are tuned is central to understanding how phenotypic variation is realized. Prophage induction by phage lambda is the classic example of a genetic switch and allows lambda to move between two different modes of transmission: as a lysogen it reproduces vertically as a component of the host genome; as a free phage it reproduces horizontally by infectious epidemic spread. We show that the lambda switch can respond rapidly to selection for alteration in sensitivity and threshold. Sequencing of candidate genes in the genetic circuitry underlying the switch revealed mutations of likely adaptive significance in some, but not all candidates, suggesting that the core genetic circuitry plays a limited role in the fine-tuning of the switch in vivo. The relative ease with which the switch could be tuned by selection was further indicated by extensive variation in sensitivity and threshold of its response function among wild lambdoid phages. Together, our findings emphasize the adaptive significance of a finely tuned switch and draw attention to the selective factors shaping prophage induction in natural phage populations.