Gene order constrains adaptation in bacteriophage T7

Genome-scale top-down strategy to generate viable genome-reduced phages

Efforts have been made to reduce the genomes of living cells, but phage genome reduction remains challenging. It is of great interest to investigate whether genome reduction can make phages obtain new infectious properties. We developed a CRISPR/Cas9-based iterative phage genome reduction (CiPGr) approach and applied this to four distinct phages, thereby obtaining heterogeneous genome-reduced mutants. We isolated and sequenced 200 mutants with loss of up to 8-23% (3.3-35 kbp) of the original sequences. This allowed the identification of non-essential genes for phage propagation, although loss of these genes is mostly detrimental to phage fitness to various degrees. Notwithstanding this, mutants with higher infectious efficiency than their parental strains were characterized, indicating a trade-off between genome reduction and infectious fitness for phages. In conclusion, this study provides a foundation for future work to leverage the information generated by CiPGr in phage synthetic biology research.

Generating dynamic gene expression patterns without the need for regulatory circuits

Synthetic biology has successfully advanced our ability to design and implement complex, time-varying genetic circuits to control the expression of recombinant proteins. However, these circuits typically require the production of regulatory genes whose only purpose is to coordinate expression of other genes. When designing very small genetic constructs, such as viral genomes, we may want to avoid introducing such auxiliary gene products while nevertheless encoding complex expression dynamics. To this end, here we demonstrate that varying only the placement and strengths of promoters, terminators, and RNase cleavage sites in a computational model of a bacteriophage genome is sufficient to achieve solutions to a variety of basic gene expression patterns. We discover these genetic solutions by computationally evolving genomes to reproduce desired gene expression time-course data. Our approach shows that non-trivial patterns can be evolved, including patterns where the relative ordering of genes by abundance changes over time. We find that some patterns are easier to evolve than others, and comparable expression patterns can be achieved via different genetic architectures. Our work opens up a novel avenue to genome engineering via fine-tuning the balance of gene expression and gene degradation rates.

Genome Modularization Reveals Overlapped Gene Topology Is Necessary for Efficient Viral Reproduction

Sequence overlap between two genes is common across all genomes, with viruses having high proportions of these gene overlaps. Genome modularization and refactoring is the process of disrupting natural gene overlaps to separate coding sequences to enable their individual manipulation. The biological function and fitness effects of gene overlaps are not fully understood, and their effects on gene cluster and genome-level refactoring are unknown. The bacteriophage φX174 genome has ∼26% of nucleotides involved in encoding more than one gene. In this study we use an engineered φX174 phage containing a genome with all gene overlaps removed to show that gene overlap is critical to maintaining optimal viral fecundity. Through detailed phenotypic measurements we reveal that genome modularization in φX174 causes virion replication, stability, and attachment deficiencies. Quantitation of the complete phage proteome across an infection cycle reveals 30% of proteins display abnormal expression patterns. Taken together, we have for the first time comprehensively demonstrated that gene modularization severely perturbs the coordinated functioning of a bacteriophage replication cycle. This work highlights the biological importance of gene overlap in natural genomes and that reducing gene overlap disruption should be an integral part of future genome engineering projects.

Transcript degradation and codon usage regulate gene expression in a lytic phage

Many viral genomes are small, containing only single- or double-digit numbers of genes and relatively few regulatory elements. Yet viruses successfully execute complex regulatory programs as they take over their host cells. Here, we propose that some viruses regulate gene expression via a carefully balanced interplay between transcription, translation, and transcript degradation. As our model system, we employ bacteriophage T7, whose genome of approximately sixty genes is well annotated and for which there is a long history of computational models of gene regulation. We expand upon prior modeling work by implementing a stochastic gene expression simulator that tracks individual transcripts, polymerases, ribosomes, and ribonucleases participating in the transcription, translation, and transcript-degradation processes occurring during a T7 infection. By combining this detailed mechanistic modeling of a phage infection with high-throughput gene expression measurements of several strains of bacteriophage T7, evolved and engineered, we can show that both the dynamic interplay between transcription and transcript degradation, and between these two processes and translation, appear to be critical components of T7 gene regulation. Our results point to targeted degradation as a generic gene regulation strategy that may have evolved in many other viruses. Further, our results suggest that detailed mechanistic modeling may uncover the biological mechanisms at work in both evolved and engineered virus variants.

Viral attenuation by engineered protein fragmentation

A possible but untested method of viral attenuation is protein fragmentation, engineering wild-type proteins as two or more peptides that self-assemble after translation. Here, the bacteriophage T7 was engineered to encode its essential RNA polymerase as two peptides. Initial fitness was profoundly suppressed. Subjecting the engineered virus to over 100 generations of adaptation by serial transfer resulted in a large fitness increase, still remaining below that of evolved wild-type. The fitness increase was accompanied by three substitutions in the fragmented peptides as well as six mutations in other parts of the genome, but the fragmentation was retained. This study thereby demonstrates the feasibility of using gene fragmentation as a possibly permanent method of attenuation, but the initial fitness of the engineered genome may be a poor measure of its fitness on extended adaptation.

Reduced Protein Expression in a Virus Attenuated by Codon Deoptimization

A general means of viral attenuation involves the extensive recoding of synonymous codons in the viral genome. The mechanistic underpinnings of this approach remain unclear, however. Using quantitative proteomics and RNA sequencing, we explore the molecular basis of attenuation in a strain of bacteriophage T7 whose major capsid gene was engineered to carry 182 suboptimal codons. We do not detect transcriptional effects from recoding. Proteomic observations reveal that translation is halved for the recoded major capsid gene, and a more modest reduction applies to several coexpressed downstream genes. We observe no changes in protein abundances of other coexpressed genes that are encoded upstream. Viral burst size, like capsid protein abundance, is also decreased by half. Together, these observations suggest that, in this virus, reduced translation of an essential polycistronic transcript and diminished virion assembly form the molecular basis of attenuation.

Multiple Barriers to the Evolution of Alternative Gene Orders in a Positive-Strand RNA Virus

The order in which genes are organized within a genome is generally not conserved between distantly related species. However, within virus orders and families, strong conservation of gene order is observed. The factors that constrain or promote gene-order diversity are largely unknown, although the regulation of gene expression is one important constraint for viruses. Here we investigate why gene order is conserved for a positive-strand RNA virus encoding a single polyprotein in the context of its authentic multicellular host. Initially, we identified the most plausible trajectory by which alternative gene orders could evolve. Subsequently, we studied the accessibility of key steps along this evolutionary trajectory by constructing two virus intermediates: (1) duplication of a gene followed by (2) loss of the ancestral gene. We identified five barriers to the evolution of alternative gene orders. First, the number of viable positions for reordering is limited. Second, the within-host fitness of viruses with gene duplications is low compared to the wild-type virus. Third, after duplication, the ancestral gene copy is always maintained and never the duplicated one. Fourth, viruses with an alternative gene order have even lower fitness than viruses with gene duplications. Fifth, after more than half a year of evolution in isolation, viruses with an alternative gene order are still vastly inferior to the wild-type virus. Our results show that all steps along plausible evolutionary trajectories to alternative gene orders are highly unlikely. Hence, the inaccessibility of these trajectories probably contributes to the conservation of gene order in present-day viruses.

Analysis of whole genome sequencing for the Escherichia coli O157:H7 typing phages

Background

Shiga toxin producing Escherichia coli O157 can cause severe bloody diarrhea and haemolytic uraemic syndrome. Phage typing of E. coli O157 facilitates public health surveillance and outbreak investigations, certain phage types are more likely to occupy specific niches and are associated with specific age groups and disease severity. The aim of this study was to analyse the genome sequences of 16 (fourteen T4 and two T7) E. coli O157 typing phages and to determine the genes responsible for the subtle differences in phage type profiles.

Results

The typing phages were sequenced using paired-end Illumina sequencing at The Genome Analysis Centre and the Animal Health and Veterinary Laboratories Agency and bioinformatics programs including Velvet, Brig and Easyfig were used to analyse them. A two-way Euclidian cluster analysis highlighted the associations between groups of phage types and typing phages. The analysis showed that the T7 typing phages (9 and 10) differed by only three genes and that the T4 typing phages formed three distinct groups of similar genomic sequences: Group 1 (1, 8, 11, 12 and 15, 16), Group 2 (3, 6, 7 and 13) and Group 3 (2, 4, 5 and 14). The E. coli O157 phage typing scheme exhibited a significantly modular network linked to the genetic similarity of each group showing that these groups are specialised to infect a subset of phage types.

Conclusion

Sequencing the typing phage has enabled us to identify the variable genes within each group and to determine how this corresponds to changes in phage type.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1470-z) contains supplementary material, which is available to authorized users.

Evolutionary reversion of live viral vaccines: Can genetic engineering subdue it?

Attenuated, live viral vaccines have been extraordinarily successful in protecting against many diseases. The main drawbacks in their development and use have been reliance on an unpredictable method of attenuation and the potential for evolutionary reversion to high virulence. Methods of genetic engineering now provide many safer alternatives to live vaccines, so if live vaccines are to compete with these alternatives in the future, they must either have superior immunogenicity or they must be able to overcome these former disadvantages. Several live vaccine designs that were historically inaccessible are now feasible because of advances in genome synthesis. Some of those methods are addressed here, with an emphasis on whether they enable predictable levels of attenuation and whether they are stable against evolutionary reversion. These new designs overcome many of the former drawbacks and position live vaccines to be competitive with alternatives. Not only do new methods appear to retard evolutionary reversion enough to prevent vaccine-derived epidemics, but it may even be possible to permanently attenuate live vaccines that are transmissible but cannot evolve to higher virulence under prolonged adaptation.

Phylogenomic network and comparative genomics reveal a diverged member of the ΦKZ-related group, marine vibrio phage ΦJM-2012

Bacteriophages are the largest reservoir of genetic diversity. Here we describe the novel phage ΦJM-2012. This natural isolate from marine Vibrio cyclitrophicus possesses very few gene contents relevant to other well-studied marine Vibrio phages. To better understand its evolutionary history, we built a mathematical model of pairwise relationships among 1,221 phage genomes, in which the genomes (nodes) are linked by edges representing the normalized number of shared orthologous protein families. This weighted network revealed that ΦJM-2012 was connected to only five members of the Pseudomonas ΦKZ-like phage family in an isolated network, strongly indicating that it belongs to this phage group. However, comparative genomic analyses highlighted an almost complete loss of colinearity with the ΦKZ-related genomes and little conservation of gene order, probably reflecting the action of distinct evolutionary forces on the genome of ΦJM-2012. In this phage, typical conserved core genes, including six RNA polymerase genes, were frequently displaced and the hyperplastic regions were rich in both unique genes and predicted unidirectional promoters with highly correlated orientations. Further, analysis of the ΦJM-2012 genome showed that segments of the conserved N-terminal parts of ΦKZ tail fiber paralogs exhibited evidence of combinatorial assortment, having switched transcriptional orientation, and there was recruitment and/or structural changes among phage endolysins and tail spike protein. Thus, this naturally occurring phage appears to have branched from a common ancestor of the ΦKZ-related groups, showing a distinct genomic architecture and unique genes that most likely reflect adaptation to its chosen host and environment.

Evolutionarily stable attenuation by genome rearrangement in a virus

Live, attenuated viruses provide many of the most effective vaccines. For the better part of a century, the standard method of attenuation has been viral growth in novel environments, whereby the virus adapts to the new environment but incurs a reduced ability to grow in the original host. The downsides of this approach were that it produced haphazard results, and even when it achieved sufficient attenuation for vaccine production, the attenuated virus was prone to evolve back to high virulence. Using bacteriophage T7, we apply a synthetic biology approach for creating attenuated genomes and specifically study their evolutionary stability. Three different genome rearrangements are used, and although some initial fitness recovery occurs, all exhibit greatly impaired abilities to recover wild-type fitness over a hundred or more generations. Different degrees of stable attenuation appear to be attainable by different rearrangements. Efforts to predict fitness recovery using the extensive background of T7 genetics and biochemistry were only sometimes successful. The use of genome rearrangement thus offers a practical mechanism of evolutionary stable viral attenuation, with some progress toward prediction.

Multiple genetic pathways to similar fitness limits during viral adaptation to a new host

The gain in fitness during adaptation depends on the supply of beneficial mutations. Despite a good theoretical understanding of how evolution proceeds for a defined set of mutations, there is little understanding of constraints on net fitness-whether fitness will reach a limit despite ongoing selection and mutation, and if there is a limit, what determines it. Here, the dsDNA bacteriophage SP6, a virus of Salmonella, was adapted to Escherichia coli K-12. From an isolate capable of modest growth on E. coli, four lines were adapted for rapid growth by protocols differing in use of mutagen, propagation method, and duration, but using the same media, temperature, and a continual excess of the novel host. Nucleotide changes underlying those adaptations differed greatly in number and identity, but the four lines achieved similar absolute fitness at the end, an increase of more than 4000-fold phage descendants per hour. Thus, the fitness landscape allows multiple genetic paths to the same approximate fitness limit. The existence and causes of fitness limits have ramifications to genome engineering, vaccine design, and "lethal mutagenesis" treatments to cure viral infections.

Multiplicity of infection does not accelerate infectivity evolution of viral parasites in laboratory microcosms

Coinfection with multiple parasite genotypes [multiplicity of infection (MOI)] creates within-host competition and opportunities for parasite recombination and is therefore predicted to be important for both parasite and host evolution. We tested for a difference in the infectivity of viral parasites (lytic phage Φ2) and resistance of their bacterial hosts (Pseudomonas fluorescens SBW25) under both high and low MOI during coevolution in laboratory microcosms. Results show that MOI has no effect on infectivity and resistance evolution during coevolution over ∼80 generations of host growth, and this is true when the experiment is initiated with wild-type viruses and hosts, or with viruses and hosts that have already been coevolving for ∼330 generations. This suggests that MOI does not have a net effect of accelerating parasite adaptation to hosts through recombination, or slowing adaptation to hosts through between-parasite conflict in this system.

Experimental evolution, loss-of-function mutations, and "the first rule of adaptive evolution"

Adaptive evolution can cause a species to gain, lose, or modify a function; therefore, it is of basic interest to determine whether any of these modes dominates the evolutionary process under particular circumstances. Because mutation occurs at the molecular level, it is necessary to examine the molecular changes produced by the underlying mutation in order to assess whether a given adaptation is best considered as a gain, loss, or modification of function. Although that was once impossible, the advance of molecular biology in the past half century has made it feasible. In this paper, I review molecular changes underlying some adaptations, with a particular emphasis on evolutionary experiments with microbes conducted over the past four decades. I show that by far the most common adaptive changes seen in those examples are due to the loss or modification of a pre-existing molecular function, and I discuss the possible reasons for the prominence of such mutations.

Layers of evolvability in a bacteriophage life history trait

Functional redundancy in genomes arises from genes with overlapping functions, allowing phenotypes to persist after gene knockouts. Evolutionary redundancy or evolvability of a genome is one step removed, in that functional redundancy is absent but the genome has the potential to evolve to restore a lost phenotype. Exploring the extent to which this recovery alters gene networks can illuminate how functional gene interactions change through time. Here, the evolvability of lysis was studied in bacteriophage T7, revealing hidden functional interactions. Lysis is the destruction of host cell wall and membranes that releases progeny and is therefore essential for phage propagation. In most phages, lysis is mediated by a two-component genetic module: a muralytic enzyme that degrades the bacterial cell wall (endolysin) and a holin that permeabilizes the inner membrane to allow the endolysin access to the cell wall. T7 carries one known holin, one endolysin, and a second muralytic enzyme that plays little role in lysis by wild-type phage. If the primary endolysin is deleted, the second muralytic enzyme evolves to restore normal lysis after selection for faster growth. Here, a second level of evolutionary redundancy was revealed. When the second muralytic enzyme was prevented from adapting in a genome lacking the primary endolysin, the phage reevolved lysis de novo in the absence of any known muralytic enzymes by changes in multiple genes outside the original lysis module. This second level of redundancy proved to be evolutionarily inferior to the first, and both result in a lower fitness and slower lysis than wild-type T7. Deletion of the holin gene delayed lysis time modestly; fitness was restored by compensatory substitutions in genes that lack known roles in lysis of the wild type.

Exploring tree-building methods and distinct molecular data to recover a known asymmetric phage phylogeny

An experimental phylogeny was constructed using bacteriophage T7 and a propagation protocol, in the presence of the mutagen N-methyl-N'-nitro-N'-nitrosoguanidine, based on Hillis et al. [Hillis, D.M., Bull, J.J., White, M.E., Badgett, M.R., Molineux, I.J., 1992. Experimental phylogenetics, generation of a known phylogeny. Science 255, 589-592]. The topology presented in this study has a considerable variation in branch lengths and is less symmetric than the one presented by Hillis et al. [Hillis, D.M., Bull, J.J., White, M.E., Badgett, M.R., Molineux, I.J., 1992. Experimental phylogenetics, generation of a known phylogeny. Science 255, 589-592]. These features are known to present additional difficulties to phylogenetic inference methods. The performance of several phylogenetic methods (conventional and less conventional) was tested using restriction site and nucleotide data. Only methods that encompassed a molecular clock or those based on sequence signatures recovered the true phylogeny. Nevertheless a likelihood ratio test rejected the hypothesis of the existence of a molecular clock when the whole sequence data set was considered. This fact or the particular substitution pattern (mainly G-->A and C-->T) may be related to the unexpected performance of distance methods based on sequence signatures. To test if the results could have been predicted by simulation studies we estimated the evolution parameters from the real phylogeny and used them to simulate evolution along the same tree (parametric bootstrap). We found that simulation could predict most but not all of the problems encountered by phylogenetic inference methods in the real phylogeny. Short interior branches may be more prone to error than predicted by theoretical studies.

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.

Complete genomic sequence of bacteriophage phiEcoM-GJ1, a novel phage that has myovirus morphology and a podovirus-like RNA polymerase

The complete genome of phiEcoM-GJ1, a lytic phage that attacks porcine enterotoxigenic Escherichia coli of serotype O149:H10:F4, was sequenced and analyzed. The morphology of the phage and the identity of the structural proteins were also determined. The genome consisted of 52,975 bp with a G+C content of 44% and was terminally redundant and circularly permuted. Seventy-five potential open reading frames (ORFs) were identified and annotated, but only 29 possessed homologs. The proteins of five ORFs showed homology with proteins of phages of the family Myoviridae, nine with proteins of phages of the family Podoviridae, and six with proteins of phages of the family Siphoviridae. ORF 1 encoded a T7-like single-subunit RNA polymerase and was preceded by a putative E. coli sigma(70)-like promoter. Nine putative phage promoters were detected throughout the genome. The genome included a tRNA gene of 95 bp that had a putative 18-bp intron. The phage morphology was typical of phages of the family Myoviridae, with an icosahedral head, a neck, and a long contractile tail with tail fibers. The analysis shows that phiEcoM-GJ1 is unique, having the morphology of the Myoviridae, a gene for RNA polymerase, which is characteristic of phages of the T7 group of the Podoviridae, and several genes that encode proteins with homology to proteins of phages of the family Siphoviridae.

Comparative genomics of Bacillus thuringiensis phage 0305phi8-36: defining patterns of descent in a novel ancient phage lineage

BACKGROUND

The recently sequenced 218 kb genome of morphologically atypical Bacillus thuringiensis phage 0305phi8-36 exhibited only limited detectable homology to known bacteriophages. The only known relative of this phage is a string of phage-like genes called BtI1 in the chromosome of B. thuringiensis israelensis. The high degree of divergence and novelty of phage genomes pose challenges in how to describe the phage from its genomic sequences.

RESULTS

Phage 0305phi8-36 and BtI1 are estimated to have diverged 2.0 - 2.5 billion years ago. Positionally biased Blast searches aligned 30 homologous structure or morphogenesis genes between 0305phi8-36 and BtI1 that have maintained the same gene order. Functional clustering of the genes helped identify additional gene functions. A conserved long tape measure gene indicates that a long tail is an evolutionarily stable property of this phage lineage. An unusual form of the tail chaperonin system split to two genes was characterized, as was a hyperplastic homologue of the T4gp27 hub gene. Within this region some segments were best described as encoding a conservative array of structure domains fused with a variable component of exchangeable domains. Other segments were best described as multigene units engaged in modular horizontal exchange. The non-structure genes of 0305phi8-36 appear to include the remnants of two replicative systems leading to the hypothesis that the genome plan was created by fusion of two ancestral viruses. The case for a member of the RNAi RNA-directed RNA polymerase family residing in 0305phi8-36 was strengthened by extending the hidden Markov model of this family. Finally, it was noted that prospective transcriptional promoters were distributed in a gradient of small to large transcripts starting from a fixed end of the genome.

CONCLUSION

Genomic organization at a level higher than individual gene sequence comparison can be analyzed to aid in understanding large phage genomes. Methods of analysis include 1) applying a time scale, 2) augmenting blast scores with positional information, 3) categorizing genomic rearrangements into one of several processes with characteristic rates and outcomes, and 4) correlating apparent transcript sizes with genomic position, gene content, and promoter motifs.