Multi-species host range of staphylococcal phages isolated from wastewater

The host range of bacteriophages defines their impact on bacterial communities and genome diversity. Here, we characterize 94 novel staphylococcal phages from wastewater and establish their host range on a diversified panel of 117 staphylococci from 29 species. Using this high-resolution phage-bacteria interaction matrix, we unveil a multi-species host range as a dominant trait of the isolated staphylococcal phages. Phage genome sequencing shows this pattern to prevail irrespective of taxonomy. Network analysis between phage-infected bacteria reveals that hosts from multiple species, ecosystems, and drug-resistance phenotypes share numerous phages. Lastly, we show that phages throughout this network can package foreign genetic material enclosing an antibiotic resistance marker at various frequencies. Our findings indicate a weak host specialism of the tested phages, and therefore their potential to promote horizontal gene transfer in this environment.

The Bacillus phage SPβ and its relatives: A temperate phage model system reveals new strains, species, prophage integration loci, conserved proteins and lysogeny management components

The Bacillus phage SPβ has been known for about 50 years, but only a few strains are avalible. We isolated four new wild type strains of the SPbeta species. Phage vB_BsuS-Goe14 introduces its prophage into the spoVK locus, previously not observed to be used by SPβ-like phages. We could also reveal the SPβ-like phage genome replication strategy, the genome packaging mode, and the phage genome opening point. We extracted 55 SPβ-like prophages from public Bacillus genomes, thereby discovering three more integration loci and one additional type of integrase. The identified prophages resembled four new species clusters and three species orphans in the genus Spbetavirus. The determined core proteome of all SPβ-like prophages consists of 38 proteins. The integration cassette proved to be not conserved even though present in all strains. It consists of distinct integrases. Analysis of SPβ transcriptomes revealed three conserved genes, yopQ, yopR, and yokI, to be transcribed from a dormant prophage. While yopQ and yokI could be deleted from the prophage without activating the prophage, damaging of yopR led to a clear-plaque phenotype. Under the applied laboratory conditions, the yokI mutant showed an elevated virion release implying the YokI protein being a component of the arbitrium system.

Categorization of prophage genes in Bacillus subtilis 168 and assessing their relative importance through RNA-seq gene expression analysis

Bacteriophage evolves to control the population of fast-growing bacterial cells, without which explosion in bacterial population may induce unimaginable harm to diverse ecosystems. But, bacteriophage also hide in bacterial genomes when nutritional and environmental circumstances are unfavourable. This involves the integration of phage genome into the host genome at appropriate genomic loci in a process known as lysogeny. This work sought to delineate the prophages present in the annotated genome of Bacillus subtilis 168, and assess their relative importance through RNA-seq expression analysis. Firstly, examination of the annotated genome of the model Gram-positive bacterium revealed five distinct prophage regions: SPBeta, prophage 6, PBSX, prophage 3 region, and prophage 1 region. All prophage regions contain host genes, which suggests that host transposase activity have swapped in host genes for phage genes in the prophage genome. Given the significant number of phage genes that have been swapped into each of the prophage genome, all prophage regions are deemed to be defective. BLAST analysis further highlighted that many of the prophages in B. subtilis are extinct given that they do not have ancestral or daughter brethren. However, RNA-seq transcriptome analysis of B. subtilis turned out an interesting paradox indicative of the important role that host transposase have in swapping in host promoters for prophage genes. Specifically, a significant number of prophage genes are highly expressed, which is implausible given that phage genes should be transcriptionally silent. The result and phenomenon further suggests the relative facile nature in which host promoters could be swapped in for phage genes, which is indicative of presence of genomic motifs in prophage genome recognizable by host transposase. Existence of such sequence motifs is thus indicative of possible co-evolution of transposase and phages where transposases were originally a part of the phage genome, which latter jumped out into the host genome to aid the swapping in of host genes into the prophage genome for augmenting prophage genetic repertoire in the face of changing environmental conditions. Overall, it is not uncommon for bacterial species to harbour multiple prophages. But, lysogeny may not be a viable option for long-term preservation of prophage genetic repertoire given that host transposase would inevitable swap in host genes at random locations in the prophage genome.

DeePVP: Identification and classification of phage virion protein using deep learning

The poor annotation of phage virion protein (PVP) is the bottleneck of many areas of viral research, such as viral phylogenetic analysis, viral host identification and antibacterial drug design. Because of the high diversity of the PVP sequences, the PVP annotation remains a great challenging bioinformatic task. Based on deep learning, we present DeePVP that contains a main module and an extended module. The main module aims to identify the PVPs from non-PVP over a phage genome, while the extended module can further classify the predicted PVP into one of the ten major classes of PVP. Compared with the state-of-the-art tools that can distinguish PVP from non-PVP, DeePVP's main module performs much better, with an F1-score 9.05% higher in the PVP identification task. Compared with PhANNs, a tool that can further classify the predicted PVP into a specific class, the overall accuracy of DeePVP's extended module is approximately 3.72% higher in the PVP classification task. Two application cases on the genome of mycobacteriophage PDRPxv and Escherichia phage HP3 show that the predictions of DeePVP are much more reliable and can better reveal the compact PVP-enriched region, which may be conserved during the viral evolution process, over the phage genome.

The Characterization of a Novel Phage, pPa_SNUABM_DT01, Infecting Pseudomonas aeruginosa

The bacterial genus Pseudomonas is a common causative agent of infections in veterinary medicine. In this study, we focused on Pseudomonas aeruginosa canine otitis externa isolates. Due to prolonged antibiotic treatment of otitis externa, antibiotic resistance is common and has become a major complication. Many alternatives to antibiotics have been studied, with bacteriophages emerging as the most promising alternatives. Here, we isolated and characterized a novel phage, pPa_SNUABM_DT01, by investigating its morphology, growth, lysis kinetics, and genomic characteristics. Phages have a vigorous capacity to eliminate bacterial cells through bacterial lysis. This capacity is dependent on the multiplicity of infection (MOI), but even at low MOIs, the phage successfully inhibited bacterial regrowth. The phage genome was 265,520 bp in size and comprised 312 putative open reading frames (ORFs). Comparative genome analysis demonstrated that the phage is a novel species in Myoviridae. The nucleotide similarity was moderately high compared with the Pseudomonas virus, Noxifer. However, a phylogenetic analysis and a dot plot indicated that pPa_SNUABM_DT01 is not closely related to the Phikzvirus or Noxifervirus genus but, instead, belongs to a novel one. The genome comparisons also indicate that the phage, pPa_SNUABM_DT01, could be a novel genus.

Engineering a dynamic, controllable infectivity switch in bacteriophage T7

Transcriptional repressors play an important role in regulating phage genomes. Here, we examined how synthetic regulation based on repressors can be to create a dynamic, controllable infectivity switch in bacteriophage T7. We engineered T7 by replacing a large region of the early phage genome with combinations of ligand-responsive promoters and ribosome binding sites (RBS) designed to control the phage RNA polymerase. Phages with the engineered switch showed virulence comparable to wildtype when not repressed, indicating the phage can be engineered without a loss of fitness. When repressed, the most effective switch used a TetR promoter and a weak RBS, resulting in a two-fold increase in latent period (time to lyse host) and change in phage titer. Further, phage activity could be tuned by varying inducer concentrations. Our study provides a proof of concept for a simple circuit for user control over phage infectivity.

Characterization of the first Pseudomonas grimontii bacteriophage, PMBT3

The complete genome sequence of the virulent bacteriophage PMBT3, isolated on the proteolytic Pseudomonas grimontii strain MBTL2-21, showed no significant similarity to other known phage genome sequences, making this phage the first reported to infect a strain of P. grimontii. Electron microscopy revealed PMBT3 to be a member of the family Siphoviridae, with notably long and flexible whiskers. The linear, double-stranded genome of 87,196 bp has a mol% G+C content of 60.4 and contains 116 predicted protein-encoding genes. A putative tellurite resistance (terB) gene, originally reported to occur in the genome of a bacterium, was detected in the genome of phage PMBT3.

Mining bacterial NGS data vastly expands the complete genomes of temperate phages

Temperate phages (active prophages induced from bacteria) help control pathogenicity, modulate community structure, and maintain gut homeostasis. Complete phage genome sequences are indispensable for understanding phage biology. Traditional plaque techniques are inapplicable to temperate phages due to the lysogenicity of these phages, which curb the identification and characterization of temperate phages. Existing in silico tools for prophage prediction usually fail to detect accurate and complete temperate phage genomes. In this study, by a novel computational method mining both the integrated active prophages and their spontaneously induced forms (temperate phages), we obtained 192,326 complete temperate phage genomes from bacterial next-generation sequencing (NGS) data, hence expanded the existing number of complete temperate phage genomes by more than 100-fold. The reliability of our method was validated by wet-lab experiments. The experiments demonstrated that our method can accurately determine the complete genome sequences of the temperate phages, with exact flanking sites (attP and attB sites), outperforming other state-of-the-art prophage prediction methods. Our analysis indicates that temperate phages are likely to function in the evolution of microbes by 1) cross-infecting different bacterial host species; 2) transferring antibiotic resistance and virulence genes; and 3) interacting with hosts through restriction-modification and CRISPR/anti-CRISPR systems. This work provides a comprehensive complete temperate phage genome database and relevant information, which can serve as a valuable resource for phage research.

Genome-scale top-down reduction of phages to generate viable minimal phage genomes

Reduction of tailed-phage genomes to generate viable minimal genome phages is important for expanding our understanding of phage biology, providing insights for phage synthetic biology. Many efforts have been made to minimize living cells, but such work remains a challenge for phages due to the extraordinary genomic diversity and lack of genome-scale editing techniques. Here, we developed a CRISPR/Cas9-based iterative phage genome reduction (CiPGr) approach to detect the nonessential gene set of phages and minimize phage genomes. By CiPGr, inactivated genes accumulated on the phage genome, and mutant progeny with robust growth gradually arose, eventually becoming predominant in the populations. CiPGr was applied to four distinct tailed phages (model phages T7 and T4; wild-type phages seszw and selz), resulting in mutants of these phages with deletion of 8–20% (3.3–33 kbp) sequences, and leading to minimal genomes. Metagenomic sequencing of the mutant phage populations generated showed that 46.7 to 65.4% of genes of these phages were removed. Loss of some genes (39.6%-50%) in the removable gene sets was likely severely detrimental to phage growth. This made the corresponding mutant progenies recede in the populations, leading to the failure of detection of these genes in the genomes of the isolated mutants. In summary, our results for these four distinct tailed phages demonstrated that CiPGr is a generic yet effective approach suitable for use in novel phages without prior knowledge.