The Effect of Mutation and Selection on Codon Adaptation in Escherichia coli Bacteriophage

Studying phage codon adaptation is important not only for understanding the process of translation elongation, but also for reengineering phages for medical and industrial purposes. To evaluate the effect of mutation and selection on phage codon usage, we developed an index to measure selection imposed by host translation machinery, based on the difference in codon usage between all host genes and highly expressed host genes. We developed linear and nonlinear models to estimate the C→T mutation bias in different phage lineages and to evaluate the relative effect of mutation and host selection on phage codon usage. C→T-biased mutations occur more frequently in single-stranded DNA (ssDNA) phages than in double-stranded DNA (dsDNA) phages and affect not only synonymous codon usage, but also nonsynonymous substitutions at second codon positions, especially in ssDNA phages. The host translation machinery affects codon adaptation in both dsDNA and ssDNA phages, with a stronger effect on dsDNA phages than on ssDNA phages. Strand asymmetry with the associated local variation in mutation bias can significantly interfere with codon adaptation in both dsDNA and ssDNA phages.     

Horizontal gene transfer and the evolution of microvirid coliphage genomes

Bacteriophage genomic evolution has been largely characterized by rampant, promiscuous horizontal gene transfer involving both homologous and nonhomologous source DNA. This pattern has emerged through study of the tailed double-stranded DNA (dsDNA) phages and is based upon a sparse sampling of the enormous diversity of these phages. The single-stranded DNA phages of the family Microviridae, including phiX174, appear to evolve through qualitatively different mechanisms, possibly as result of their strictly lytic lifestyle and small genome size. However, this apparent difference could reflect merely a dearth of relevant data. We sought to characterize the forces that contributed to the molecular evolution of the Microviridae and to examine the genetic structure of this single family of bacteriophage by sequencing the genomes of microvirid phage isolated on a single bacterial host. Microvirids comprised 3.5% of the detectable phage in our environmental samples, and sequencing yielded 42 new microvirid genomes. Phylogenetic analysis of the genes contained in these and five previously described microvirid phages identified three distinct clades and revealed at least two horizontal transfer events between clades. All members of one clade have a block of five putative genes that are not present in any member of the other two clades. Our data indicate that horizontal transfer does contribute to the evolution of the microvirids but is both quantitatively and qualitatively different from what has been observed for the dsDNA phages.     

Erratum: Causes and Consequences of Bacteriophage Diversification via Genetic Exchanges across Lifestyles and Bacterial Taxa

Bacteriophages (phages) evolve rapidly by acquiring genes from other phages. This results in mosaic genomes. Here, we identify numerous genetic transfers between distantly related phages and aim at understanding their frequency, consequences, and the conditions favoring them. Gene flow tends to occur between phages that are enriched for recombinases, transposases, and nonhomologous end joining, suggesting that both homologous and illegitimate recombination contribute to gene flow. Phage family and host phyla are strong barriers to gene exchange, but phage lifestyle is not. Even if we observe four times more recent transfers between temperate phages than between other pairs, there is extensive gene flow between temperate and virulent phages, and between the latter. These predominantly involve virulent phages with large genomes previously classed as low gene flux, and lead to the preferential transfer of genes encoding functions involved in cell energetics, nucleotide metabolism, DNA packaging and injection, and virion assembly. Such exchanges may contribute to the observed twice larger genomes of virulent phages. We used genetic transfers, which occur upon coinfection of a host, to compare phage host range. We found that virulent phages have broader host ranges and can mediate genetic exchanges between narrow host range temperate phages infecting distant bacterial hosts, thus contributing to gene flow between virulent phages, as well as between temperate phages. This gene flow drastically expands the gene repertoires available for phage and bacterial evolution, including the transfer of functional innovations across taxa.     

Motif depletion in bacteriophages infecting hosts with CRISPR systems

Background

CRISPR is a microbial immune system likely to be involved in host-parasite coevolution. It functions using target sequences encoded by the bacterial genome, which interfere with invading nucleic acids using a homology-dependent system. The system also requires protospacer associated motifs (PAMs), short motifs close to the target sequence that are required for interference in CRISPR types I and II. Here, we investigate whether PAMs are depleted in phage genomes due to selection pressure to escape recognition.

Results

To this end, we analyzed two data sets. Phages infecting all bacterial hosts were analyzed first, followed by a detailed analysis of phages infecting the genus Streptococcus, where PAMs are best understood. We use two different measures of motif underrepresentation that control for codon bias and the frequency of submotifs. We compare phages infecting species with a particular CRISPR type to those infecting species without that type. Since only known PAMs were investigated, the analysis is restricted to CRISPR types I-C and I-E and in Streptococcus to types I-C and II. We found evidence for PAM depletion in Streptococcus phages infecting hosts with CRISPR type I-C, in Vibrio phages infecting hosts with CRISPR type I-E and in Streptococcus thermopilus phages infecting hosts with type II-A, known as CRISPR3.

Conclusions

The observed motif depletion in phages with hosts having CRISPR can be attributed to selection rather than to mutational bias, as mutational bias should affect the phages of all hosts. This observation implies that the CRISPR system has been efficient in the groups discussed here.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-663) contains supplementary material, which is available to authorized users.     

Computational Prediction of Genomic Functional Cores Specific to Different Microbes

Computational and experimental attempts tried to characterize a universial core of genes representing the minimal set of functional needs for an organism. Based on the increasing number of available complete genomes, comparative genomics has concluded that the universal core contains < 50 genes. In contrast, experiments suggest a much larger set of essential genes (certainly more than several hundreds, even under the most restrictive hypotheses) that is dependent on the biological complexity and environmental specificity of the organism. Highly biased genes, which are generally also the most expressed in translationally biased organisms, tend to be over represented in the class of genes deemed to be essential for any given bacterial species. This association is far from perfect; nevertheless, it allows us to propose a new computational method to detect, to a certain extent, ubiquitous genes, nonorthologous genes, environment-specific genes, genes involved in the stress response, and genes with no identified function but highly likely to be essential for the cell. Most of these groups of genes cannot be identified with previously attempted computational and experimental approaches. The large variety of life-styles and the unusually detectable functional signals characterizing translationally biased organisms suggest using them as reference organisms to infer essentiality in other microbial species. The case of small parasitic genomes is discussed. Data issued by the analysis are compared with previous computational and experimental studies. Results are discussed both on methodological and biological grounds. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Martin Kreitman]     

Evaluation of codon bias perspectives in phage therapy of Mycobacterium tuberculosis by multivariate analysis

To reveal the relative synonymous codon usage and base composition variation in bacteriophages, six mycobacteriophages were used as a model system here and both parameters in these phages and their host bacteria, Mycobacterium tuberculosis, have been determined and compared. As expected for GC-rich genomes, there are predominantly G and C ending codons in all 6 phages. Both N_{c} plot and correspondence analysis on relative synonymous codon usage indicate that mutation bias and translation selection influences codon usage variation in the 6 phages. Further analysis indicates that among 6 Mycobacterium phages Che9c, Bxz1 and TM4 may be extremely virulent in nature as most of their genes have high translation efficiency. Based on our data we suggest that the genes of above three phages are expressed rapidly by host's translation machinery. The information might be used to select the extremely virulent Mycobacterium tuberculosis phages suitable for phage therapy.     

Viral adaption of staphylococcal phage: A genome-based analysis of the selective preference based on codon usage Bias

Given the high therapeutic value of the staphylococcal phage, the genome co-evolution of the phage and the host has gained great attention. Though the genome-wide AT richness in staphylococcal phages has been well-studied with nucleotide usage bias, here we proved that host factor, lifestyle and taxonomy are also important factors in understanding the phage nucleotide usages bias using information entropy formula. Such correlation is especially prominent when it comes to the synonymous codon usages of staphylococcal phages, despite the overall scattered codon usage pattern represented by principal component analysis. This strong relationship is explained by nucleotide skew which testified that the usage biases of nucleotide at different codon positions are acting on synonymous codons. Therefore, our study reveals a hidden relationship of genome evolution with host limitation and phagic phenotype, providing new insight into phage genome evolution at genetic level.     

Core and accessory genome architecture in a group of Pseudomonas aeruginosa Mu-like phages

Background

Bacteriophages that infect the opportunistic pathogen Pseudomonas aeruginosa have been classified into several groups. One of them, which includes temperate phage particles with icosahedral heads and long flexible tails, bears genomes whose architecture and replication mechanism, but not their nucleotide sequences, are like those of coliphage Mu. By comparing the genomic sequences of this group of P. aeruginosa phages one could draw conclusions about their ontogeny and evolution.

Results

Two newly isolated Mu-like phages of P. aeruginosa are described and their genomes sequenced and compared with those available in the public data banks. The genome sequences of the two phages are similar to each other and to those of a group of P. aeruginosa transposable phages. Comparing twelve of these genomes revealed a common genomic architecture in the group. Each phage genome had numerous genes with homologues in all the other genomes and a set of variable genes specific for each genome. The first group, which comprised most of the genes with assigned functions, was named “core genome”, and the second group, containing mostly short ORFs without assigned functions was called “accessory genome”. Like in other phage groups, variable genes are confined to specific regions in the genome.

Conclusion

Based on the known and inferred functions for some of the variable genes of the phages analyzed here, they appear to confer selective advantages for the phage survival under particular host conditions. We speculate that phages have developed a mechanism for horizontally acquiring genes to incorporate them at specific loci in the genome that help phage adaptation to the selective pressures imposed by the host.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1146) contains supplementary material, which is available to authorized users.     

Codon engineering for improved antibody expression in mammalian cells

While well established in bacterial hosts, the effect of coding sequence variation on protein expression in mammalian systems is poorly characterized outside of viral proteins or proteins from distant phylogenetic families. The potential impact is substantial given the extensive use of mammalian expression systems in research and manufacturing of protein biotherapeutics. We are studying the effect of codon engineering on expression of recombinant antibodies with an emphasis on developing manufacturing cell lines. CNTO 888, a human mAb specific for human MCP-1, was obtained by antibody phage display in collaboration with MorphoSys AG. The isolated DNA sequence of the antibody was biased towards bacterial codons, reflecting the engineering of the Fab library for phage display expression in Escherichia coli. We compared the expression of CNTO 888 containing the parental V-region sequences with two engineered coding variants. In the native codon exchanged (NCE) variant, the V-region codons were replaced with those used in naturally derived human antibody genes. In the human codon optimized (HCO) variant the V-region codons were those used at the highest frequency based on a human codon usage table. The antibody expression levels from stable transfections in mammalian host cells were measured. The HCO codon variant of CNTO 888 yielded the highest expressing cell lines and the highest average expression for the screened populations. This had a significant positive effect on the process to generate a CNTO 888 production cell line and indicates the potential to improve antibody expression in mammalian expression systems by codon engineering.     

Evolutionarily conserved orthologous families in phages are relatively rare in their prokaryotic hosts

We have identified conserved orthologs in completely sequenced genomes of double-strand DNA phages and arranged them into evolutionary families (phage orthologous groups [POGs]). Using this resource to analyze the collection of known phage genomes, we find that most orthologs are unique in their genomes (having no diverged duplicates [paralogs]), and while many proteins contain multiple domains, the evolutionary recombination of these domains does not appear to be a major factor in evolution of these orthologous families. The number of POGs has been rapidly increasing over the past decade, the percentage of genes in phage genomes that have orthologs in other phages has also been increasing, and the percentage of unknown "ORFans" is decreasing as more proteins find homologs and establish a family. Other properties of phage genomes have remained relatively stable over time, most notably the high fraction of genes that are never or only rarely observed in their cellular hosts. This suggests that despite the renowned ability of phages to transduce cellular genes, these cellular "hitchhiker" genes do not dominate the phage genomic landscape, and a large fraction of the genes in phage genomes maintain an evolutionary trajectory that is distinct from that of the host genes.     

Marine phage genomics

Marine phages are the most abundant biological entities in the oceans. They play important roles in carbon cycling through marine food webs, gene transfer by transduction and conversion of hosts by lysogeny. The handful of marine phage genomes that have been sequenced to date, along with prophages in marine bacterial genomes, and partial sequencing of uncultivated phages are yielding glimpses of the tremendous diversity and physiological potential of the marine phage community. Common gene modules in diverse phages are providing the information necessary to make evolutionary comparisons. Finally, deciphering phage genomes is providing clues about the adaptive response of phages and their hosts to environmental cues.     

Revelation of Influencing Factors in Overall Codon Usage Bias of Equine Influenza Viruses

Equine influenza viruses (EIVs) of H3N8 subtype are culprits of severe acute respiratory infections in horses, and are still responsible for significant outbreaks worldwide. Adaptability of influenza viruses to a particular host is significantly influenced by their codon usage preference, due to an absolute dependence on the host cellular machinery for their replication. In the present study, we analyzed genome-wide codon usage patterns in 92 EIV strains, including both H3N8 and H7N7 subtypes by computing several codon usage indices and applying multivariate statistical methods. Relative synonymous codon usage (RSCU) analysis disclosed bias of preferred synonymous codons towards A/U-ended codons. The overall codon usage bias in EIVs was slightly lower, and mainly affected by the nucleotide compositional constraints as inferred from the RSCU and effective number of codon (ENc) analysis. Our data suggested that codon usage pattern in EIVs is governed by the interplay of mutation pressure, natural selection from its hosts and undefined factors. The H7N7 subtype was found less fit to its host (horse) in comparison to H3N8, by possessing higher codon bias, lower mutation pressure and much less adaptation to tRNA pool of equine cells. To the best of our knowledge, this is the first report describing the codon usage analysis of the complete genomes of EIVs. The outcome of our study is likely to enhance our understanding of factors involved in viral adaptation, evolution, and fitness towards their hosts.     

Local biotic environment shapes the spatial scale of bacteriophage adaptation to bacteria

The ecological, epidemiological, and evolutionary consequences of host-parasite interactions are critically shaped by the spatial scale at which parasites adapt to hosts. The scale of interaction between hyperparasites and their parasites is likely to be influenced by the host of the parasite and potentially likely to differ among within-host environments. Here we examine the scale at which bacteriophages adapt to their host bacteria by studying natural isolates from the surface or interior of horse chestnut leaves. We find that phages are more infective to bacteria from the same tree relative to those from other trees but do not differ in infectivity to bacteria from different leaves within the same tree. The results suggest that phages target common bacterial species, including an important plant pathogen, within plant host tissues; this result has important implications for therapeutic phage epidemiology. Furthermore, we show that phages from the leaf interior are more infective to their local hosts than phages from the leaf surface are to theirs, suggesting either increased resistance of bacteria on the leaf surface or increased phage adaptation within the leaf. These results highlight that biotic environment can play a key role in shaping the spatial scale of parasite adaptation and influencing the outcome of coevolutionary interactions.     

Reticulate representation of evolutionary and functional relationships between phage genomes

Bacteriophage genomes show pervasive mosaicism, indicating the importance of horizontal gene exchange in their evolution. Phage genomes represent unique combinations of modules, each of them with a different phylogenetic history. The traditional classification, based on a variety of criteria such as nucleic acid type (single/double-stranded DNA/RNA), morphology, and host range, appeared inconsistent with sequence analyses. With the genomic era, an ever increasing number of sequenced phages cannot be classified, in part due to a lack of morphological information and in part to the intrinsic incapability of tree-based methods to efficiently deal with mosaicism. This problem led some virologists to call for a moratorium on the creation of additional taxa in the order Caudovirales, in order to let virologists discuss classification schemes that might better suit phage evolution. In this context, we propose a framework for a reticulate classification of phages based on gene content. Starting from gene families, we built a weighted graph, where nodes represent phages and edges represent phage-phage similarities in terms of shared genes. We then apply various measures of graph topology to analyze the resulting graph. Most double-stranded DNA phages are found in a single component. The values of the clustering coefficient and closeness distinguish temperate from virulent phages, whereas chimeric phages are characterized by a high betweenness coefficient. We apply a 2-step clustering method to this graph to generate a reticulate classification of phages: Each phage is associated with a membership vector, which quantitatively characterizes its membership to the set of clusters. Furthermore, we cluster genes based on their "phylogenetic profiles" to define "evolutionary cohesive modules." In virulent phages, evolutionary modules span several functional categories, whereas in temperate phages they correspond better to functional modules. Moreover, despite the fact that modules only cover a fraction of all phage genes, phage groups can be distinguished by their different combination of modules, serving the bases for a higher level reticulate classification. These 2 classification schemes provide an automatic and dynamic way of representing the relationships within the phage population and can be extended to include newly sequenced phage genomes, as well as other types of genetic elements.     

Plasticity of the gene functions for DNA replication in the T4-like phages

We have completely sequenced and annotated the genomes of several relatives of the bacteriophage T4, including three coliphages (RB43, RB49 and RB69), three Aeromonas salmonicida phages (44RR2.8t, 25 and 31) and one Aeromonas hydrophila phage (Aeh1). In addition, we have partially sequenced and annotated the T4-like genomes of coliphage RB16 (a close relative of RB43), A. salmonicida phage 65, Acinetobacter johnsonii phage 133 and Vibrio natriegens phage nt-1. Each of these phage genomes exhibited a unique sequence that distinguished it from its relatives, although there were examples of genomes that are very similar to each other. As a group the phages compared here diverge from one another by several criteria, including (a) host range, (b) genome size in the range between approximately 160 kb and approximately 250 kb, (c) content and genetic organization of their T4-like genes for DNA metabolism, (d) mutational drift of the predicted T4-like gene products and their regulatory sites and (e) content of open-reading frames that have no counterparts in T4 or other known organisms (novel ORFs). We have observed a number of DNA rearrangements of the T4 genome type, some exhibiting proximity to putative homing endonuclease genes. Also, we cite and discuss examples of sequence divergence in the predicted sites for protein-protein and protein-nucleic acid interactions of homologues of the T4 DNA replication proteins, with emphasis on the diversity in sequence, molecular form and regulation of the phage-encoded DNA polymerase, gp43. Five of the sequenced phage genomes are predicted to encode split forms of this polymerase. Our studies suggest that the modular construction and plasticity of the T4 genome type and several of its replication proteins may offer resilience to mutation, including DNA rearrangements, and facilitate the adaptation of T4-like phages to different bacterial hosts in nature.     

Schematic for efficient computation of GC, GC3, and AT3 bias spectra of genome

Selection of synonymous codons for an amino acid is biased in protein translation process. This biased selection causes repetition of synonymous codons in structural parts of genome that stands for high N/3 peaks in DNA spectrum. Period-3 spectral property is utilized here to produce a 3-phase network model based on polyphase filterbank concepts for derivation of codon bias spectra (CBS). Modification of parameters in this model can produce GC, GC3, and AT3 bias spectra. Complete schematic in LabVIEW platform is presented here for efficient and parallel computation of GC, GC3, and AT3 bias spectra of genomes alongwith results of CBS patterns. We have performed the correlation coefficient analysis of GC, GC3, and AT3 bias spectra with codon bias patterns of CBS for biological and statistical significance of this model.     

Rule-based simulation of temperate bacteriophage infection: Restriction–modification as a limiter to infection in bacterial populations

An individual-based model (IbM) for bacterial adaptation and evolution, COSMIC-Rules, has been employed to simulate interactions of virtual temperate bacteriophages (phages) and their bacterial hosts. Outcomes of infection mimic those of a phage such as lambda, which can enter either the lytic or lysogenic cycle, depending on the nutritional status of the host. Infection of different hosts possessing differing restriction and modification systems is also simulated. Phages restricted upon infection of one restricting host can be adapted (by host-controlled modification of the phage genome) and subsequently propagate with full efficiency on this host. However, such ability is lost if the progeny phages are passaged through a new host with a different restriction and modification system before attempted re-infection of the original restrictive host. The simulations show that adaptation and re-adaptation to a particular host-controlled restriction and modification system result in lower efficiency and delayed lysis of bacterial cells compared with infection of non-restricting host bacteria.