Phage Genome Annotation Using the RAST Pipeline

Isolation, characterization, and application of a lytic bacteriophage SSP49 to control Staphylococcus aureus contamination on baby spinach leaves

Staphylococcus aureus, a major foodborne pathogen, is frequently detected in fresh produce. It often causes food poisoning accompanied by abdominal pain, diarrhea, and vomiting. Additionally, the abuse of antibiotics to control S. aureus has resulted in the emergence of antibiotics-resistant bacteria, such as methicillin resistant S. aureus. Therefore, bacteriophage, a natural antimicrobial agent, has been suggested as an alternative to antibiotics. In this study, a lytic phage SSP49 that specifically infects S. aureus was isolated from a sewage sample, and its morphological, biological, and genetic characteristics were determined. We found that phage SSP49 belongs to the Straboviridae family (Caudoviricetes class) and maintained host growth inhibition for 30 h in vitro. In addition, it showed high host specificity and a broad host range against various S. aureus strains. Receptor analysis revealed that phage SSP49 utilized cell wall teichoic acid as a host receptor. Whole genome sequencing revealed that the genome size of SSP49 was 137,283 bp and it contained 191 open reading frames. The genome of phage SSP49 did not contain genes related to lysogen formation, bacterial toxicity, and antibiotic resistance, suggesting its safety in food application. The activity of phage SSP49 was considerably stable under various high temperature and pH conditions. Furthermore, phage SSP49 effectively inhibited S. aureus growth on baby spinach leaves both at 4 °C and 25 °C while maintaining the numbers of active phage during treatments (reductions of 1.2 and 2.1 log CFU/cm2, respectively). Thus, this study demonstrated the potential of phage SSP49 as an alternative natural biocontrol agent against S. aureus contamination in fresh produce.

Physiochemical characterization of a potential Klebsiella phage MKP-1 and analysis of its application in reducing biofilm formation

The common intestinal pathogen Klebsiella pneumoniae (K. pneumoniae) is one of the leading causes of fatal superbug infections that can resist the effects of commonly prescribed medicines. The uncontrolled use or misuse of antibiotics has increased the prevalence of drug-resistant K. pneumoniae strains in the environment. In the quest to search for alternative therapeutics for treating these drug-resistant infections, bacteriophages (bacterial viruses) emerged as potential candidates for in phage therapy against Klebsiella. The effective formulation of phage therapy against drug-resistant Klebsiella infections demands thorough characterization and screening of many bacteriophages. To contribute effectively to the formulation of successful phage therapy against superbug infections by K. pneumoniae, this study includes the isolation and characterization of a novel lytic bacteriophage MKP-1 to consider its potential to be used as therapeutics in treating drug-resistant Klebsiella infections. Morphologically, having a capsid attached to a long non-contractile tail, it was found to be a siphovirus that belongs to the class Caudoviricetes and showed infectivity against different strains of the target host bacterium. Comparatively, this double-stranded DNA phage has a large burst size and is quite stable in various physiological conditions. More interestingly, it has the potential to degrade the tough biofilms formed by K. pneumoniae (Klebsiella pneumoniae subsp. pneumoniae (Schroeter) Trevisan [ATCC 15380]) significantly. Thus, the following study would contribute effectively to considering phage MKP-1 as a potential candidate for phage therapy against Klebsiella infection.

Isolation and characterization of two Staphylococcus aureus lytic bacteriophages "Huma" and "Simurgh"

Nowadays finding the new antimicrobials is necessary due to the emerging of multidrug resistant strains. The present study aimed to isolate and characterize bacteriophages against S. aureus. Strains Huma and Simurgh were the two podovirus morphology phages which isolated and then characterized. Huma and Simurgh had a genome size of 16,853 and 17,245 bp, respectively and both were Rosenblumvirus with G + C content of 29%. No lysogeny-related genes, nor virulence genes were identified in their genomes. They were lytic only against two out of four S. aureus strains. They also were able to inhibit S. aureus for 8 h in-vitro. Both showed a rapid adsorption. Huma and Simurgh had the latent period of 80 and 60 m and the burst sizes of 45 and 40 PFU/ml and also, they showed very low cell toxicity of 1.23%-1.79% on HT-29 cells, respectively. Thus, they can be considered potential candidates for biocontrol applications.

Isolation and Characterization of Six Novel Fusobacterium necrophorum Phages

Introduction

Fusobacterium necrophorum, a human and animal pathogen, is the primary etiologic agent of bovine liver abscesses and a driving factor for prophylactic antibiotic use in the fed cattle industry. Considering calls to reduce agricultural antibiotic use, we isolated phages capable of killing F. necrophorum as an alternative or complementary biocontrol strategy.

Methods

Six novel phages (φFN37, φRTG5, φKSUM, φHugo, φPaco, and φBB) were isolated from rumen fluid or ruminal F. necrophorum isolates and subjected to host range testing on both F. necrophorum subspecies. Four F. necrophorum subspecies, necrophorum phages, were tested for cross-resistance and host growth inhibition individually and in pairs. Additionally, genomic sequencing, annotation, and analysis were performed.s.

Results

Four of six isolated phages were able to form lysogens, although all six contained lysogeny-related genes. φKSUM and φBB, did not form lysogens and were able to infect both subspecies. Four phages could infect F. necrophorum 8L1 (a liver abscess model challenge strain) in vitro. Genomic analysis showed that these phages belong to class Caudoviricetes with genome sizes ranging from 35 kbp to 111 kbp and GC values ranging from 26% to 36% and have extremely limited similarity to other deposited phage genomes infecting Fusobacterium or other genera.

Conclusions

Although all phages isolated contained sequences bearing similarities to genes implicated in lysogeny, the four selected for use in cocktails showed potential in inhibiting host growth, with several demonstrating promising attributes for biocontrol and therapeutic applications. Phage cocktails that may offer enhanced antibacterial activity were also identified, indicating the potential of some lysogenic phages to be adapted for biocontrol or therapeutic purposes when lytic phages are difficult to obtain.

Phage resistance in Klebsiella pneumoniae and bidirectional effects impacting antibiotic susceptibility

OBJECTIVES

Bacteriophage (phage) therapy is a promising anti-infective option to combat antimicrobial resistance. However, the clinical utilization of phage therapy has been severely compromised by the potential emergence of phage resistance. Although certain phage resistance mechanisms can restore bacterial susceptibility to certain antibiotics, a lack of knowledge of phage resistance mechanisms hinders optimal use of phages and their combination with antibiotics.

METHODS

Genome-wide transposon screening was performed with a mutant library of Klebsiella pneumoniae MKP103 to identify phage pKMKP103_1-resistant mutants. Phage-resistant phenotypes were evaluated by time-kill kinetics and efficiency of plating assays. Phage resistance mechanisms were investigated with adsorption, one-step growth, and mutation frequency assays. Antibiotic susceptibility was determined with broth microdilution and population analysis profiles.

RESULTS

We observed a repertoire of phage resistance mechanisms in K pneumoniae, such as disruption of phage binding (fhuA::Tn and tonB::Tn), extension of the phage latent period (mnmE::Tn and rpoN::Tn), and increased mutation frequency (mutS::Tn and mutL::Tn). Notably, in contrast to the prevailing view that phage resistance re-sensitizes antibiotic-resistant bacteria, we observed a bidirectional steering effect on bacterial antibiotic susceptibility. Specifically, rpoN::Tn increased susceptibility to colistin while mutS::Tn and mutL::Tn increased resistance to rifampicin and colistin.

DISCUSSION

Our findings demonstrate that K pneumoniae employs multiple strategies to overcome phage infection, which may result in enhanced or reduced antibiotic susceptibility. Mechanism-guided phage steering should be incorporated into phage therapy to better inform clinical decisions on phage-antibiotic combinations.

Complete genome sequences of cluster F1 and cluster B1 Mycobacterium smegmatis phages Karhdo and Basato

We present the complete genome sequences of Mycobacterium smegmatis phages Karhdo and Basato, isolated in Clark County, Nevada. The phages were isolated and annotated by students enrolled in undergraduate research courses over two semesters at the University of Nevada, Las Vegas.

A novel lytic phage exhibiting a remarkable in vivo therapeutic potential and higher antibiofilm activity against Pseudomonas aeruginosa

BACKGROUND

Pseudomonas aeruginosa is a nosocomial bacterium responsible for variety of infections. Inappropriate use of antibiotics could lead to emergence of multidrug-resistant (MDR) P. aeruginosa strains. Herein, a virulent phage; vB_PaeM_PS3 was isolated and tested for its application as alternative to antibiotics for controlling P. aeruginosa infections.

METHODS

Phage morphology was observed using transmission electron microscopy (TEM). The phage host range and efficiency of plating (EOP) in addition to phage stability were analyzed. One-step growth curve was performed to detect phage growth kinetics. The impact of isolated phage on planktonic cells and biofilms was assessed. The phage genome was sequenced. Finally, the therapeutic potential of vB_PaeM_PS3 was determined in vivo.

RESULTS

Isolated phage has an icosahedral head and a contractile tail and was assigned to the family Myoviridae. The phage vB_PaeM_PS3 displayed a broad host range, strong bacteriolytic ability, and higher environmental stability. Isolated phage showed a short latent period and large burst size. Importantly, the phage vB_PaeM_PS3 effectively eradicated bacterial biofilms. The genome of vB_PaeM_PS3 consists of 93,922 bp of dsDNA with 49.39% G + C content. It contains 171 predicted open reading frames (ORFs) and 14 genes as tRNA. Interestingly, the phage vB_PaeM_PS3 significantly attenuated P. aeruginosa virulence in host where the survival of bacteria-infected mice was markedly enhanced following phage treatment. Moreover, the colonizing capability of P. aeruginosa was markedly impaired in phage-treated mice as compared to untreated infected mice.

CONCLUSION

Based on these findings, isolated phage vB_PaeM_PS3 could be potentially considered for treating of P. aeruginosa infections.

Isolation and Characterization of a Weizmannia coagulans Bacteriophage Youna2 and Its Endolysin PlyYouna2

Weizmannia coagulans (formerly Bacillus coagulans) is Gram-positive, and spore-forming bacteria causing food spoilage, especially in acidic canned food products. To control W. coagulans, we isolated a bacteriophage Youna2 from a sewage sludge sample. Morphological analysis revealed that phage Youna2 belongs to the Siphoviridae family with a non-contractile and flexible tail. Youna2 has 52,903 bp double-stranded DNA containing 61 open reading frames. There are no lysogeny-related genes, suggesting that Youna2 is a virulent phage. plyYouna2, a putative endolysin gene was identified in the genome of Youna2 and predicted to be composed of a N-acetylmuramoyl-L-alanine amidase domain (PF01520) at the N-terminus and unknown function DUF5776 domain (PF19087) at the C-terminus. While phage Youna2 has a narrow host range, infecting only certain strains of W. coagulans, PlyYouna2 exhibited a broad antimicrobial spectrum beyond the Bacillus genus. Interestingly, PlyYouna2 can lyse Gram-negative bacteria such as Escherichia coli, Yersinia enterocolitica, Pseudomonas putida and Cronobacter sakazakii without other additives to destabilize bacterial outer membrane. To the best of our knowledge, Youna2 is the first W. coagulans-infecting phage and we speculate its endolysin PlyYouna2 can provide the basis for the development of a novel biocontrol agent against various foodborne pathogens.

Characteristics and Whole-Genome Analysis of Limosilactobacillus fermentum Phage LFP02

Limosilactobacillus fermentum is a bacterium widely used in food production, medicine, and industrial fermentation. However, fermentation could fail due to phage contamination. L. fermentum bacteriophage LFP02 can be induced from L. fermentum IMAU 32579 using mitomycin C. To better understand the characteristics of this phage, its physiological and genomic characteristics were evaluated. The results showed that its optimal multiplicity of infection was 0.01, and the burst size was 148.03 ± 2.65 pfu/infective center. Compared to temperature, pH had a more obvious influence on phage viability, although its adsorption capacity was not affected by the divalent cations (Ca²⁺ and Mg²⁺) or chloramphenicol. Its genome size was 43,789 bp and the GC content was 46.06%, including 53 functional proteins. Compared to other L. fermentum phages, phage LFP02 had chromosome deletion, insertion, and inversion, which demonstrated that it was a novel phage. This study could expand the knowledge of the biological characteristics of L. fermentum bacteriophages and provide some theoretical basis for bacteriophage prevention during fermentation.

Epitopes in the capsular polysaccharide and the porin OmpK36 receptors are required for bacteriophage infection of Klebsiella pneumoniae

To kill bacteria, bacteriophages (phages) must first bind to a receptor, triggering the release of the phage DNA into the bacterial cell. Many bacteria secrete polysaccharides that had been thought to shield bacterial cells from phage attack. We use a comprehensive genetic screen to distinguish that the capsule is not a shield but is instead a primary receptor enabling phage predation. Screening of a transposon library to select phage-resistant Klebsiella shows that the first receptor-binding event docks to saccharide epitopes in the capsule. We discover a second step of receptor binding, dictated by specific epitopes in an outer membrane protein. This additional and necessary event precedes phage DNA release to establish a productive infection. That such discrete epitopes dictate two essential binding events for phages has profound implications for understanding the evolution of phage resistance and what dictates host range, two issues critically important to translating knowledge of phage biology into phage therapies.

Transcriptomics-Driven Characterization of LUZ100, a T7-like Pseudomonas Phage with Temperate Features

Autographiviridae is a diverse yet distinct family of bacterial viruses marked by a strictly lytic lifestyle and a generally conserved genome organization. Here, we characterized Pseudomonas aeruginosa phage LUZ100, a distant relative of type phage T7. LUZ100 is a podovirus with a limited host range which likely uses lipopolysaccharide (LPS) as a phage receptor. Interestingly, infection dynamics of LUZ100 indicated moderate adsorption rates and low virulence, hinting at temperate characteristics. This hypothesis was supported by genomic analysis, which showed that LUZ100 shares the conventional T7-like genome organization yet carries key genes associated with a temperate lifestyle. To unravel the peculiar characteristics of LUZ100, ONT-cappable-seq transcriptomics analysis was performed. These data provided a bird's-eye view of the LUZ100 transcriptome and enabled the discovery of key regulatory elements, antisense RNA, and transcriptional unit structures. The transcriptional map of LUZ100 also allowed us to identify new RNA polymerase (RNAP)-promoter pairs that can form the basis for biotechnological parts and tools for new synthetic transcription regulation circuitry. The ONT-cappable-seq data revealed that the LUZ100 integrase and a MarR-like regulator (proposed to be involved in the lytic/lysogeny decision) are actively cotranscribed in an operon. In addition, the presence of a phage-specific promoter transcribing the phage-encoded RNA polymerase raises questions on the regulation of this polymerase and suggests that it is interwoven with the MarR-based regulation. This transcriptomics-driven characterization of LUZ100 supports recent evidence that T7-like phages should not automatically be assumed to have a strictly lytic life cycle. IMPORTANCE Bacteriophage T7, considered the "model phage" of the Autographiviridae family, is marked by a strictly lytic life cycle and conserved genome organization. Recently, novel phages within this clade have emerged which display characteristics associated with a temperate life cycle. Screening for temperate behavior is of utmost importance in fields like phage therapy, where strictly lytic phages are generally required for therapeutic applications. In this study, we applied an omics-driven approach to characterize the T7-like Pseudomonas aeruginosa phage LUZ100. These results led to the identification of actively transcribed lysogeny-associated genes in the phage genome, pointing out that temperate T7-like phages are emerging more frequent than initially thought. In short, the combination of genomics and transcriptomics allowed us to obtain a better understanding of the biology of nonmodel Autographiviridae phages, which can be used to optimize the implementation of phages and their regulatory elements in phage therapy and biotechnological applications, respectively.

Phage Encounters Recorded in CRISPR Arrays in the Genus Oenococcus

The Oenococcus genus comprises four recognized species, and members have been found in different types of beverages, including wine, kefir, cider and kombucha. In this work, we implemented two complementary strategies to assess whether oenococcal hosts of different species and habitats were connected through their bacteriophages. First, we investigated the diversity of CRISPR-Cas systems using a genome-mining approach, and CRISPR-endowed strains were identified in three species. A census of the spacers from the four identified CRISPR-Cas loci showed that each spacer space was mostly dominated by species-specific sequences. Yet, we characterized a limited records of potentially recent and also ancient infections between O. kitaharae and O. sicerae and phages of O. oeni, suggesting that some related phages have interacted in diverse ways with their Oenococcus hosts over evolutionary time. Second, phage-host interaction analyses were performed experimentally with a diversified panel of phages and strains. None of the tested phages could infect strains across the species barrier. Yet, some infections occurred between phages and hosts from distinct beverages in the O. oeni species.

Preliminary Characterization of Phage-Like Particles from the Male-Killing Mollicute Spiroplasma poulsonii (an Endosymbiont of Drosophila)

Bacteriophages are vastly abundant, diverse, and influential, but with few exceptions (e.g. the Proteobacteria genera Wolbachia and Hamiltonella), the role of phages in heritable bacteria-arthropod interactions, which are ubiquitous and diverse, remains largely unexplored. Despite prior studies documenting phage-like particles in the mollicute Spiroplasma associated with Drosophila flies, genomic sequences of such phage are lacking, and their effects on the Spiroplasma-Drosophila interaction have not been comprehensively characterized. We used a density step gradient to isolate phage-like particles from the male-killing bacterium Spiroplasma poulsonii (strains NSRO and MSRO-Br) harbored by Drosophila melanogaster. Isolated particles were subjected to DNA sequencing, assembly, and annotation. Several lines of evidence suggest that we recovered phage-like particles of similar features (shape, size, DNA content) to those previously reported in Drosophila-associated Spiroplasma strains. We recovered three ~ 19 kb phage-like contigs (two in NSRO and one in MSRO-Br) containing 21-24 open reading frames, a read-alignment pattern consistent with circular permutation, and terminal redundancy (at least in NSRO). Although our results do not allow us to distinguish whether these phage-like contigs represent infective phage-like particles capable of transmitting their DNA to new hosts, their encoding of several typical phage genes suggests that they are at least remnants of functional phage. We also recovered two smaller non-phage-like contigs encoding a known Spiroplasma toxin (Ribosome Inactivating Protein; RIP), and an insertion element, suggesting that they are packaged into particles. Substantial homology of our particle-derived contigs was found in the genome assemblies of members of the Spiroplasma poulsonii clade.

Isolation, Characterization, and Genome Analysis of a Novel Bacteriophage, Escherichia Phage vB_EcoM-4HA13, Representing a New Phage Genus in the Novel Phage Family Chaseviridae

Shiga toxin-producing Escherichia coli (STEC) is one of the leading causes of foodborne illnesses in North America and can lead to severe symptoms, with increased fatality risk for young children. While E. coli O157:H7 remains the dominant STEC serotype associated with foodborne outbreaks, there has been an increasing number of non-O157 STEC outbreaks in recent years. For the food industry, lytic bacteriophages offer an organic, self-limiting alternative to pathogen reduction-one that could replace or reduce the use of chemical and physical food processing methods. From EHEC-enriched sewage, we isolated a novel bacteriophage, vB_EcoM-4HA13 (4HA13). Phenotypic characterizations revealed 4HA13 to possess a myoviral morphotype, with a high specificity to non-motile O111 serotype, and a long latent period (90 min). Through genomic analyses, this 52,401-bp dsDNA phage was found to contain 81 CDS, but no detectable presence of antibiotic resistance, integrase, or virulence genes. A BLASTn search for each of the identified 81 CDS yielded homologues with low levels of similarity. Comparison of RNA polymerase and terminase large subunit amino acid sequences led to the proposal and acceptance of a new bacteriophage family, Chaseviridae, with 4HA13 representing a new species and genus. The discovery of this phage has broadened our current knowledge of bacteriophage diversity.

Isolation, characterization, and application of a novel polyvalent lytic phage STWB21 against typhoidal and nontyphoidal Salmonella spp

Salmonella is one of the common causal agents of bacterial gastroenteritis-related morbidity and mortality among children below 5 years and the elderly populations. Salmonellosis in humans is caused mainly by consuming contaminated food originating from animals. The genus Salmonella has several serovars, and many of them are recently reported to be resistant to multiple drugs. Therefore, isolation of lytic Salmonella bacteriophages in search of bactericidal activity has received importance. In this study, a Salmonella phage STWB21 was isolated from a lake water sample and found to be a novel lytic phage with promising potential against the host bacteria Salmonella typhi. However, some polyvalence was observed in their broad host range. In addition to S. typhi, the phage STWB21 was able to infect S. paratyphi, S. typhimurium, S. enteritidis, and a few other bacterial species such as Sh. flexneri 2a, Sh. flexneri 3a, and ETEC. The newly isolated phage STWB21 belongs to the Siphoviridae family with an icosahedral head and a long flexible non-contractile tail. Phage STWB21 is relatively stable under a wide range of pH (4–11) and temperatures (4°C–50°C) for different Salmonella serovars. The latent period and burst size of phage STWB21 against S. typhi were 25 min and 161 plaque-forming units per cell. Since Salmonella is a foodborne pathogen, the phage STWB21 was applied to treat a 24 h biofilm formed in onion and milk under laboratory conditions. A significant reduction was observed in the bacterial population of S. typhi biofilm in both cases. Phage STWB21 contained a dsDNA of 112,834 bp in length, and the GC content was 40.37%. Also, genomic analysis confirmed the presence of lytic genes and the absence of any lysogeny or toxin genes. Overall, the present study reveals phage STWB21 has a promising ability to be used as a biocontrol agent of Salmonella spp. and proposes its application in food industries.

Genomic features of a new head-tail halovirus VOLN27B infecting a Halorubrum strain

Relatively few viruses infecting haloarchaea (haloviruses) have been reported. In this study, the genome sequence of VOLN27B, a recently described archaeal tailed virus (arTV) with a myovirus morphotype was described, along with the sequence of its host, Halorubrum spp. LN27. Halovirus VOLN27B contains a linear, dsDNA genome of 76,891 bp which is predicted to encode 109 proteins and four tRNAs (tRNA^Thr, tRNA^Arg, tRNA^Gly and tRNA^Asn). The DNA G+C content of VOLN27B genome is 56.1 mol%, nearly 10% lower than that of its host strain. A 315 bp LTR (long terminal repeat) was detected in the genome. The genome of its host strain LN27 was 3,301,211 bp (chromosome and 1 plasmid) with a DNA G+C content of 68.3 mol% and 3,142 annotated protein coding genes. At least two hypothetical proviruses were detected in the genome. It lacked a CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats) locus. Sequence similarity and phylogenetic tree reconstructions placed it within the genus Halorubrum as a potential new species. VOLN27B exhibits a distinct difference in the frequency of codon usage against its host strain Halorubrum sp. LN27. The organization of VOLN27B genome shows remarkable synteny and amino acid sequence similarity to the genomes and predicted proteins of HF1-like haloviruses (genus Haloferacalesvirus) and a provirus in the genome of Halorubrum depositum Y78. VOLN27B and its host Halorubrum sp. LN27 comprise a new virus-host system from a hypersaline ecosystem and can be used to further understand the novel biology at extreme salt concentration.

Computational Tools for the Analysis of Uncultivated Phage Genomes

Over a century of bacteriophage research has uncovered a plethora of fundamental aspects of their biology, ecology, and evolution. Furthermore, the introduction of community-level studies through metagenomics has revealed unprecedented insights on the impact that phages have on a range of ecological and physiological processes. It was not until the introduction of viral metagenomics that we began to grasp the astonishing breadth of genetic diversity encompassed by phage genomes. Novel phage genomes have been reported from a diverse range of biomes at an increasing rate, which has prompted the development of computational tools that support the multilevel characterization of these novel phages based solely on their genome sequences. The impact of these technologies has been so large that, together with MAGs (Metagenomic Assembled Genomes), we now have UViGs (Uncultivated Viral Genomes), which are now officially recognized by the International Committee for the Taxonomy of Viruses (ICTV), and new taxonomic groups can now be created based exclusively on genomic sequence information. Even though the available tools have immensely contributed to our knowledge of phage diversity and ecology, the ongoing surge in software programs makes it challenging to keep up with them and the purpose each one is designed for. Therefore, in this review, we describe a comprehensive set of currently available computational tools designed for the characterization of phage genome sequences, focusing on five specific analyses: (i) assembly and identification of phage and prophage sequences, (ii) phage genome annotation, (iii) phage taxonomic classification, (iv) phage-host interaction analysis, and (v) phage microdiversity.

SGP-C: A Broad Host Range Temperate Bacteriophage; Against Salmonella gallinarum

Salmonella gallinarum is a poultry restricted-pathogen causing fowl-typhoid disease in adult birds with mortality rates up-to 80% and exhibit resistance against commonly used antibiotics. In this current study, a temperate broad host range bacteriophage SGP-C was isolated against S. gallinarum from poultry digesta. It showed infection ability in all the 15 tested field strains of S. gallinarum. The SGP-C phage produced circular, turbid plaques with alternate rings. Its optimum activity was observed at pH 7.0 and 37-42°C, with a latent period of 45 min and burst size of 187 virions/bacterial cell. The SGP-C lysogens, SGPC-L5 and SGPC-L6 exhibited super-infection immunity against the same phage, an already reported feature of lysogens. A virulence index of 0.5 and 0.001 as MV50 of SGP-C suggests its moderate virulence. The genome of SGP-C found circular double stranded DNA of 42 Kbp with 50.04% GC content, which encodes 63 ORFs. The presence of repressor gene at ORF49, and absence of tRNA sequence in SGP-C genome indicates its lysogenic nature. Furthermore, from NGS analysis of lysogens we propose that SGP-C genome might exist either as an episome, or both as integrated and temporary episome in the host cell and warrants further studies. Phylogenetic analysis revealed its similarity with Salmonella temperate phages belonging to family Siphoviridae. The encoded proteins by SGP-C genome have not showed homology with any known toxin and virulence factor. Although plenty of lytic bacteriophages against this pathogen are already reported, to our knowledge SGP-C is the first lysogenic phage against S. gallinarum reported so far.

Phage Annotation Guide: Guidelines for Assembly and High-Quality Annotation

All sequencing projects of bacteriophages (phages) should seek to report an accurate and comprehensive annotation of their genomes. This article defines 14 questions for those new to phage genomics that should be addressed before submitting a genome sequence to the International Nucleotide Sequence Database Collaboration or writing a publication.

Novel PhoH-encoding vibriophages with lytic activity against environmental Vibrio strains

Cholera is a devastating diarrheal disease that accounts for more than 10% of children's lives worldwide, but its treatment is hampered by a rise in antibiotic resistance. One promising alternative to antibiotic therapy is the use of bacteriophages to treat antibiotic-resistant cholera infections, and control Vibrio cholera in clinical cases and in the environment, respectively. Here, we report four novel, closely related environmental myoviruses, VP4, VP6, VP18, and VP24, which we isolated from two environmental toxigenic Vibrio cholerae strains from river Kuja and Usenge beach in Kenya. High-throughput sequencing followed by bioinformatics analysis indicated that the genomes of the four bacteriophages have closely related sequences, with sizes of 148,180 bp, 148,181 bp, 148,179 bp, and 148,179 bp, and a G + C content of 36.4%. The four genomes carry the phoH gene, which is overrepresented in marine cyanophages. The isolated phages displayed a lytic activity against 15 environmental, as well as one clinical, Vibrio cholerae strains. Thus, these novel lytic vibriophages represent potential biocontrol candidates for water decontamination against pathogenic Vibrio cholerae and ought to be considered for future studies of phage therapy.

Isolation and characterization of a novel bacteriophage WO from Allonemobius socius crickets in Missouri

Wolbachia are endosymbionts of numerous arthropod and some nematode species, are important for their development and if present can cause distinct phenotypes of their hosts. Prophage DNA has been frequently detected in Wolbachia, but particles of Wolbachia bacteriophages (phage WO) have been only occasionally isolated. Here, we report the characterization and isolation of a phage WO of the southern ground cricket, Allonemobius socius, and provided the first whole-genome sequence of phage WO from this arthropod family outside of Asia. We screened A. socius abdomen DNA extracts from a cricket population in eastern Missouri by quantitative PCR for Wolbachia surface protein and phage WO capsid protein and found a prevalence of 55% and 50%, respectively, with many crickets positive for both. Immunohistochemistry using antibodies against Wolbachia surface protein showed many Wolbachia clusters in the reproductive system of female crickets. Whole-genome sequencing using Oxford Nanopore MinION and Illumina technology allowed for the assembly of a high-quality, 55 kb phage genome containing 63 open reading frames (ORF) encoding for phage WO structural proteins and host lysis and transcriptional manipulation. Taxonomically important regions of the assembled phage genome were validated by Sanger sequencing of PCR amplicons. Analysis of the nucleotides sequences of the ORFs encoding the large terminase subunit (ORF2) and minor capsid (ORF7) frequently used for phage WO phylogenetics showed highest homology to phage WOAu of Drosophila simulans (94.46% identity) and WOCin2USA1 of the cherry fruit fly, Rhagoletis cingulata (99.33% identity), respectively. Transmission electron microscopy examination of cricket ovaries showed a high density of phage particles within Wolbachia cells. Isolation of phage WO revealed particles characterized by 40–62 nm diameter heads and up to 190 nm long tails. This study provides the first detailed description and genomic characterization of phage WO from North America that is easily accessible in a widely distributed cricket species.

A Tail Fiber Protein and a Receptor-Binding Protein Mediate ICP2 Bacteriophage Interactions with Vibrio cholerae OmpU

ICP2 is a virulent bacteriophage (phage) that preys on Vibrio cholerae. ICP2 was first isolated from cholera patient stool samples. Some of these stools also contained ICP2-resistant isogenic V. cholerae strains harboring missense mutations in the trimeric outer membrane porin protein OmpU, identifying it as the ICP2 receptor. In this study, we identify the ICP2 proteins that mediate interactions with OmpU by selecting for ICP2 host range mutants within infant rabbits infected with a mixture of wild-type and OmpU mutant strains. ICP2 host range mutants that can now infect OmpU mutant strains have missense mutations in the putative tail fiber gene gp25 and the putative adhesin gene gp23. Using site-specific mutagenesis, we show that single or double mutations in gp25 are sufficient to generate the host range mutant phenotype. However, at least one additional mutation in gp23 is required for robust plaque formation on specific OmpU mutants. Mutations in gp23 alone were insufficient to produce a host range mutant phenotype. All ICP2 host range mutants retained the ability to form plaques on wild-type V. cholerae cells. The strength of binding of host range mutants to V. cholerae correlated with plaque morphology, indicating that the selected mutations in gp25 and gp23 restore molecular interactions with the receptor. We propose that ICP2 host range mutants evolve by a two-step process. First, gp25 mutations are selected for their broad host range, albeit accompanied by low-level phage adsorption. Subsequent selection occurs for gp23 mutations that further increase productive binding to specific OmpU alleles, allowing for near-wild-type efficiencies of adsorption and subsequent phage multiplication. IMPORTANCE Concern over multidrug-resistant bacterial pathogens, including Vibrio cholerae, has led to renewed interest in phage biology and the potential for phage therapy. ICP2 is a genetically unique virulent phage isolated from cholera patient stool samples. It is also one of three phages in a prophylactic cocktail that have been shown to be effective in animal models of infection and the only one of the three that requires a protein receptor (OmpU). This study identifies an ICP2 tail fiber and a receptor binding protein and examines how ICP2 responds to the selective pressures of phage-resistant OmpU mutants. We found that this particular coevolutionary arms race presents fitness costs to both ICP2 and V. cholerae.

How Broad Is Enough: The Host Range of Bacteriophages and Its Impact on the Agri-Food Sector

Novel bacteriophages (phages) possessing a broad host range are consistently and routinely reported, yet there is presently no consensus on the definition of "broad host range." As phages are increasingly being used in the development of methods for the detection and biocontrol of human pathogens, it is important to address the limitations associated with the host range. For instance, unanticipated host range breadth may result in the detection of nonpathogenic targets, thereby increasing the false-positive rate. Moreover, a broad host range is generally favored in biocontrol applications despite the risk of undesirable ancillary effects against nontarget species. Here, we discuss the research progress, applications, and implications of broad host range phages with a focus on tailed broad host range phages infecting human pathogens of concern in the Agri-Food sector.

Component Parts of Bacteriophage Virions Accurately Defined by a Machine-Learning Approach Built on Evolutionary Features

Antimicrobial resistance (AMR) continues to evolve as a major threat to human health, and new strategies are required for the treatment of AMR infections. Bacteriophages (phages) that kill bacterial pathogens are being identified for use in phage therapies, with the intention to apply these bactericidal viruses directly into the infection sites in bespoke phage cocktails. Despite the great unsampled phage diversity for this purpose, an issue hampering the roll out of phage therapy is the poor quality annotation of many of the phage genomes, particularly for those from infrequently sampled environmental sources. We developed a computational tool called STEP³ to use the “evolutionary features” that can be recognized in genome sequences of diverse phages. These features, when integrated into an ensemble framework, achieved a stable and robust prediction performance when benchmarked against other prediction tools using phages from diverse sources. Validation of the prediction accuracy of STEP³ was conducted with high-resolution mass spectrometry analysis of two novel phages, isolated from a watercourse in the Southern Hemisphere. STEP³ provides a robust computational approach to distinguish specific and universal features in phages to improve the quality of phage cocktails and is available for use at http://step3.erc.monash.edu/.

IMPORTANCE In response to the global problem of antimicrobial resistance, there are moves to use bacteriophages (phages) as therapeutic agents. Selecting which phages will be effective therapeutics relies on interpreting features contributing to shelf-life and applicability to diagnosed infections. However, the protein components of the phage virions that dictate these properties vary so much in sequence that best estimates suggest failure to recognize up to 90% of them. We have utilized this diversity in evolutionary features as an advantage, to apply machine learning for prediction accuracy for diverse components in phage virions. We benchmark this new tool showing the accurate recognition and evaluation of phage component parts using genome sequence data of phages from undersampled environments, where the richest diversity of phage still lies.

Bacteriophages: from isolation to application

Bacteriophages are considered as a potential alternative to fight pathogenic bacteria during the antibiotic resistance era. With their high specificity, they are being widely used in various applications: medicine, food industry, agriculture, animal farms, biotechnology, diagnosis, etc. Many techniques have been designed by different researchers for phage isolation, purification, and amplification, each of which has strengths and weaknesses. However, all aim at having a reasonably pure phage sample that can be further characterized. Phages can be characterized based on their physiological, morphological or inactivation tests. Microscopy, in particular, has opened a wide gate not only for visualizing phage morphological structure, but also for monitoring biochemistry and behavior. Meanwhile, computational analysis of phage genomes provides more details about phage history, lifestyle, and potential for toxigenic or lysogenic conversion, which translate to safety in biocontrol and phage therapy applications. This review summarizes phage application pipelines at different levels and addresses specific restrictions and knowledge gaps in the field. Recently developed computational approaches, which are used in phage genome analysis, are critically assessed. We hope that this assessment provides researchers with useful insights for selection of suitable approaches for Phage-related research aims and applications.

Distribution of Prophages in the Oenococcus oeni Species

Oenococcus oeni is the most exploited lactic acid bacterium in the wine industry and drives the malolactic fermentation of wines. Although prophage-like sequences have been identified in the species, many are not characterized, and a global view of their integration and distribution amongst strains is currently lacking. In this work, we analyzed the complete genomes of 231 strains for the occurrence of prophages, and analyzed their size and positions of insertion. Our data show the limited variation in the number of prophages in O. oeni genomes, and that six sites of insertion within the bacterial genome are being used for site-specific recombination. Prophage diversity patterns varied significantly for different host lineages, and environmental niches. Overall, the findings highlight the pervasive presence of prophages in the O. oeni species, their role as a major source of within-species bacterial diversity and drivers of horizontal gene transfer. Our data also have implications for enhanced understanding of the prophage recombination events which occurred during evolution of O. oeni, as well as the potential of prophages in influencing the fitness of these bacteria in their distinct niches.

Characterization and Genomic Analysis of PALS2, a Novel Staphylococcus Jumbo Bacteriophage

Staphylococcus aureus is an important human pathogen that can be frequently encountered in clinical and food-processing surroundings. Among the various countermeasures, bacteriophages have been considered to be promising alternatives to antibiotics. In this study, the bacteriophage PALS2 was isolated from bird feces, and the genomic and biological characteristics of this phage were investigated. PALS2 was determined to belong to the Myoviridae family and exhibited extended host inhibition that persisted for up to 24 h with repeated bursts of 12 plaque-forming units/cell. The complete genome of PALS2 measured 268,746 base pairs (bp), indicating that PALS2 could be classified as a jumbo phage. The PALS2 genome contained 279 ORFs and 1 tRNA covering asparagine, and the majority of predicted PALS2 genes encoded hypothetical proteins. Additional genes involved in DNA replication and repair, nucleotide metabolism, and genes encoding multisubunit RNA polymerase were identified in the PALS2 genome, which is a common feature of typical jumbo phages. Comparative genomic analysis indicated that PALS2 is a phiKZ-related virus and is more similar to typical jumbo phages than to staphylococcal phages. Additionally, the effective antimicrobial activities of phage PALS2 suggest its possible use as a biocontrol agent in various clinical and food processing environments.

Isolation and Characterization of Klebsiella Phages for Phage Therapy

Introduction: Klebsiella is a clinically important pathogen causing a variety of antimicrobial resistant infections in both community and nosocomial settings, particularly pneumonia, urinary tract infection, and sepsis. Bacteriophage (phage) therapy is being considered a primary option for the treatment of drug-resistant infections of these types. Methods: We report the successful isolation and characterization of 30 novel, genetically diverse Klebsiella phages. Results: The isolated phages span six different phage families and nine genera, representing both lysogenic and lytic lifestyles. Individual Klebsiella phage isolates infected up to 11 of the 18 Klebsiella capsule types tested, and all 18 capsule-types were infected by at least one of the phages. Conclusions: Of the Klebsiella-infecting phages presented in this study, the lytic phages are most suitable for phage therapy, based on their broad host range, high virulence, short lysis period and given that they encode no known toxin or antimicrobial resistance genes. Phage isolates belonging to the Sugarlandvirus and Slopekvirus genera were deemed most suitable for phage therapy based on our characterization. Importantly, when applied alone, none of the characterized phages were able to suppress the growth of Klebsiella for more than 12 h, likely due to the inherent ease of Klebsiella to generate spontaneous phage-resistant mutants. This indicates that for successful phage therapy, a cocktail of multiple phages would be necessary to treat Klebsiella infections.

Galaxy and Apollo as a biologist-friendly interface for high-quality cooperative phage genome annotation

In the modern genomic era, scientists without extensive bioinformatic training need to apply high-power computational analyses to critical tasks like phage genome annotation. At the Center for Phage Technology (CPT), we developed a suite of phage-oriented tools housed in open, user-friendly web-based interfaces. A Galaxy platform conducts computationally intensive analyses and Apollo, a collaborative genome annotation editor, visualizes the results of these analyses. The collection includes open source applications such as the BLAST+ suite, InterProScan, and several gene callers, as well as unique tools developed at the CPT that allow maximum user flexibility. We describe in detail programs for finding Shine-Dalgarno sequences, resources used for confident identification of lysis genes such as spanins, and methods used for identifying interrupted genes that contain frameshifts or introns. At the CPT, genome annotation is separated into two robust segments that are facilitated through the automated execution of many tools chained together in an operation called a workflow. First, the structural annotation workflow results in gene and other feature calls. This is followed by a functional annotation workflow that combines sequence comparisons and conserved domain searching, which is contextualized to allow integrated evidence assessment in functional prediction. Finally, we describe a workflow used for comparative genomics. Using this multi-purpose platform enables researchers to easily and accurately annotate an entire phage genome. The portal can be accessed at https://cpt.tamu.edu/galaxy-pub with accompanying user training material.

New insights into the biodiversity of coliphages in the intestine of poultry

Despite phages' ubiquitous presence and great importance in shaping microbial communities, little is known about the diversity of specific phages in different ecological niches. Here, we isolated, sequenced, and characterized 38 Escherichia coli-infecting phages (coliphages) from poultry faeces to gain a better understanding of the coliphage diversity in the poultry intestine. All phages belonged to either the Siphoviridae or Myoviridae family and their genomes ranged between 44,324 and 173,384 bp, with a G+C content between 35.5 and 46.4%. Phylogenetic analysis was performed based on single "marker" genes; the terminase large subunit, portal protein, and exonucleases, as well as the full draft genomes. Single gene analysis resulted in six distinct clusters. Only minor differences were observed between the different phylogenetic analyses, including branch lengths and additional duplicate or triplicate subclustering. Cluster formation was according to genome size, G+C content and phage subfamily. Phylogenetic analysis based on the full genomes supported these clusters. Moreover, several of our Siphoviridae phages might represent a novel unclassified phage genus. This study allowed for identification of several novel coliphages and provides new insights to the coliphage diversity in the intestine of poultry. Great diversity was observed amongst the phages, while they were isolated from an otherwise similar ecosystem.

Experimental approaches to tracking mobile genetic elements in microbial communities

Horizontal gene transfer is an important mechanism of microbial evolution and is often driven by the movement of mobile genetic elements between cells. Due to the fact that microbes live within communities, various mechanisms of horizontal gene transfer and types of mobile elements can co-occur. However, the ways in which horizontal gene transfer impacts and is impacted by communities containing diverse mobile elements has been challenging to address. Thus, the field would benefit from incorporating community-level information and novel approaches alongside existing methods. Emerging technologies for tracking mobile elements and assigning them to host organisms provide promise for understanding the web of potential DNA transfers in diverse microbial communities more comprehensively. Compared to existing experimental approaches, chromosome conformation capture and methylome analyses have the potential to simultaneously study various types of mobile elements and their associated hosts. We also briefly discuss how fermented food microbiomes, given their experimental tractability and moderate species complexity, make ideal models to which to apply the techniques discussed herein and how they can be used to address outstanding questions in the field of horizontal gene transfer in microbial communities.

Development of a broad-spectrum Salmonella phage cocktail containing Viunalike and Jerseylike viruses isolated from Thailand

Salmonella is one of the most common agents of foodborne disease worldwide. As natural alternatives to traditional antimicrobial agents, bacteriophages (phages) are emerging as highly effective biocontrol agents against Salmonella and other foodborne bacteria. Due to the high diversity within the Salmonella genus and emergence of drug resistant strains, improved efforts are necessary to find broad range and strictly lytic Salmonella phages for use in food biocontrol. Here, we describe the isolation and characterization of two Salmonella phages: ST-W77 isolated on S. Typhimurium and SE-W109 isolated on S. Enteritidis with extraordinary Salmonella specificity. Whole genome sequencing identified ST-W77 as a Myovirus within the Viunalikevirus genus and SE-W109 as a Siphovirus within the Jerseylikevirus genus. Infectivity studies using a panel of S. Typhimurium cell wall mutants revealed both phages require the lipopolysaccharide O-antigen, with SE-W109 also recognizing the flagella, during infection of Salmonella. A combination of both phages was capable of prolonged (one-week) antibacterial activity when added to milk or chicken meat contaminated with Salmonella. Due to their broad host ranges, strictly lytic lifestyles and lack of lysogeny-related genes or virulence genes in their genomes, ST-W77 and SE-W109 are ideal phages for further development as Salmonella biocontrol agents for food production.

Integrative omics analysis of Pseudomonas aeruginosa virus PA5oct highlights the molecular complexity of jumbo phages

Pseudomonas virus vB_PaeM_PA5oct is proposed as a model jumbo bacteriophage to investigate phage-bacteria interactions and is a candidate for phage therapy applications. Combining hybrid sequencing, RNA-Seq and mass spectrometry allowed us to accurately annotate its 286,783 bp genome with 461 coding regions including four non-coding RNAs (ncRNAs) and 93 virion-associated proteins. PA5oct relies on the host RNA polymerase for the infection cycle and RNA-Seq revealed a gradual take-over of the total cell transcriptome from 21% in early infection to 93% in late infection. PA5oct is not organized into strictly contiguous regions of temporal transcription, but some genomic regions transcribed in early, middle and late phases of infection can be discriminated. Interestingly, we observe regions showing limited transcription activity throughout the infection cycle. We show that PA5oct upregulates specific bacterial operons during infection including operons pncA-pncB1-nadE involved in NAD biosynthesis, psl for exopolysaccharide biosynthesis and nap for periplasmic nitrate reductase production. We also observe a downregulation of T4P gene products suggesting mechanisms of superinfection exclusion. We used the proteome of PA5oct to position our isolate amongst other phages using a gene-sharing network. This integrative omics study illustrates the molecular diversity of jumbo viruses and raises new questions towards cellular regulation and phage-encoded hijacking mechanisms.

Draft genome sequence data of a T7like phage 3A_8767 isolated from wastewater of a butcher house near Palar river

Phage 3A_8767 is a newly isolates phage from river water sample against Salmonella typhi 8767 (MTCC). The genome of the phage is linear, double stranded and 38,821 bp long in size. A total 49 functional ORF (open reading frame) were annotated and no tRNA was predicted. Phage 3A_8767 has icosahedral shaped head with stubby tail which comes under family Podoviridae, and genera T7 like virus.

Isolation, characterization, and comparative genomic analysis of a phage infecting high-level aminoglycoside-resistant (HLAR) Enterococcus faecalis

Enterococcus is a genus of Gram-positive bacteria that are commensal to the gastrointestinal tracts of humans but some species have been increasingly implicated as agents of nosocomial infections. The increase in infections and the spread of antibiotic-resistant strains have contributed to renewed interest in the discovery of Enterococcus phages. The aims of this study were (1) the isolation, characterization, and genome sequencing of a phage capable of infecting an antibiotic-resistant E. faecalis strain, and (2) the comparative genomic analysis of publicly-available Enterococcus phages. For this purpose, multiple phages were isolated from wastewater treatment plant (WWTP) influent using a high-level aminoglycoside-resistant (HLAR) E. faecalis strain as the host. One phage, phiNASRA1, demonstrated a high lytic efficiency (∼97.52%). Transmission electron microscopy (TEM) and whole-genome sequencing (WGS) showed that phiNASRA1 belongs to the Siphoviridae family of double-stranded DNA viruses. The phage was approximately 250 nm in length and its complete genome (40,139 bp, 34.7% GC) contained 62 open reading frames (ORFs). Phylogenetic comparisons of phiNASRA1 and 31 publicly-available Enterococcus phages, based on the large subunit terminase and portal proteins, grouped phage by provenance, size, and GC content. In particular, both phylogenies grouped phages larger than 100 kbp into distinct clades. A phylogeny based on a pangenome analysis of the same 32 phages also grouped phages by provenance, size, and GC content although agreement between the two single-locus phylogenies was higher. Per the pangenome phylogeny, phiNASRA1 was most closely related to phage LY0322 that was similar in size, GC content, and number of ORFs (40,139 and 40,934 bp, 34.77 and 34.80%, and 60 and 64 ORFs, respectively). The pangenome analysis did illustrate the high degree of sequence diversity and genome plasticity as no coding sequence was homologous across all 32 phages, and even 'conserved' structural proteins (e.g., the large subunit terminase and portal proteins) were homologous in no more than half of the 32 phage genomes. These findings contribute to a growing body of literature devoted to understanding phage biology and diversity. We propose that this high degree of diversity limited the value of the single-locus and pangenome phylogenies. By contrast, the high degree of homology between phages larger than 100 kbp suggests that pangenome analyses of more similar phages is a viable method for assessing subclade diversity. Future work is focused on validating phiNASRA1 as a potential therapeutic agent to eradicate antibiotic-resistant E. faecalis infections in an animal model.

The PATRIC Bioinformatics Resource Center: expanding data and analysis capabilities

The PathoSystems Resource Integration Center (PATRIC) is the bacterial Bioinformatics Resource Center funded by the National Institute of Allergy and Infectious Diseases (https://www.patricbrc.org). PATRIC supports bioinformatic analyses of all bacteria with a special emphasis on pathogens, offering a rich comparative analysis environment that provides users with access to over 250 000 uniformly annotated and publicly available genomes with curated metadata. PATRIC offers web-based visualization and comparative analysis tools, a private workspace in which users can analyze their own data in the context of the public collections, services that streamline complex bioinformatic workflows and command-line tools for bulk data analysis. Over the past several years, as genomic and other omics-related experiments have become more cost-effective and widespread, we have observed considerable growth in the usage of and demand for easy-to-use, publicly available bioinformatic tools and services. Here we report the recent updates to the PATRIC resource, including new web-based comparative analysis tools, eight new services and the release of a command-line interface to access, query and analyze data.

vB_PaeM_MIJ3, a Novel Jumbo Phage Infecting Pseudomonas aeruginosa, Possesses Unusual Genomic Features

Phages are the most abundant biological entity on Earth. There are many variants in phage virion sizes, morphology, and genome sizes. Large virion sized phages, with genome sizes greater than 200 kbp have been identified and termed as Jumbo phages. These phages exhibit certain characteristics that have not been reported in phages with smaller genomes. In this work, a jumbo phage named MIJ3 (vB_PaeM_MIJ3) that infects Pseudomonas aeruginosa PAO1 was isolated from an equine livery yard in Leicestershire, United Kingdom. The genome and biological characteristics of this phage have been investigated. MIJ3 is a Myovirus with multiple long tail fibers. Assessment of the host range of MIJ3 revealed that it has the ability to infect many clinical isolates of P. aeruginosa. Bioinformatics analysis of the phage genome indicated that MIJ3 is closely related to the Pseudomonas phage, PA5oct. MIJ3 possesses several unusual features that are either rarely present in other phages or have not yet been reported. In particular, MIJ3 encodes a FtsH-like protein, and a putative lysidine synthase, TilS. These two proteins have not been reported in phages. MIJ3 also possesses a split DNA polymerase B with a novel intein. Of particular interest, unlike other jumbo phages infecting Pseudomonas spp., MIJ3 lacks the genetic elements required for the formation of the phage nucleus, which was believed to be conserved across jumbo Pseudomonas phages.

Constructing and Characterizing Bacteriophage Libraries for Phage Therapy of Human Infections

Phage therapy requires libraries of well-characterized phages. Here we describe the generation of phage libraries for three target species: Escherichia coli, Pseudomonas aeruginosa, and Enterobacter cloacae. The basic phage characteristics on the isolation host, sequence analysis, growth properties, and host range and virulence on a number of contemporary clinical isolates are presented. This information is required before phages can be added to a phage library for potential human use or sharing between laboratories for use in compassionate use protocols in humans under eIND (emergency investigational new drug). Clinical scenarios in which these phages can potentially be used are discussed. The phages presented here are currently being characterized in animal models and are available for eINDs.

PHANOTATE: a novel approach to gene identification in phage genomes

MOTIVATION

Currently there are no tools specifically designed for annotating genes in phages. Several tools are available that have been adapted to run on phage genomes, but due to their underlying design, they are unable to capture the full complexity of phage genomes. Phages have adapted their genomes to be extremely compact, having adjacent genes that overlap and genes completely inside of other longer genes. This non-delineated genome structure makes it difficult for gene prediction using the currently available gene annotators. Here we present PHANOTATE, a novel method for gene calling specifically designed for phage genomes. Although the compact nature of genes in phages is a problem for current gene annotators, we exploit this property by treating a phage genome as a network of paths: where open reading frames are favorable, and overlaps and gaps are less favorable, but still possible. We represent this network of connections as a weighted graph, and use dynamic programing to find the optimal path.

RESULTS

We compare PHANOTATE to other gene callers by annotating a set of 2133 complete phage genomes from GenBank, using PHANOTATE and the three most popular gene callers. We found that the four programs agree on 82% of the total predicted genes, with PHANOTATE predicting more genes than the other three. We searched for these extra genes in both GenBank's non-redundant protein database and all of the metagenomes in the sequence read archive, and found that they are present at levels that suggest that these are functional protein-coding genes.

AVAILABILITY AND IMPLEMENTATION

https://github.com/deprekate/PHANOTATE.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Long-read based de novo assembly of low-complexity metagenome samples results in finished genomes and reveals insights into strain diversity and an active phage system

BACKGROUND

Complete and contiguous genome assemblies greatly improve the quality of subsequent systems-wide functional profiling studies and the ability to gain novel biological insights. While a de novo genome assembly of an isolated bacterial strain is in most cases straightforward, more informative data about co-existing bacteria as well as synergistic and antagonistic effects can be obtained from a direct analysis of microbial communities. However, the complexity of metagenomic samples represents a major challenge. While third generation sequencing technologies have been suggested to enable finished metagenome-assembled genomes, to our knowledge, the complete genome assembly of all dominant strains in a microbiome sample has not been demonstrated. Natural whey starter cultures (NWCs) are used in cheese production and represent low-complexity microbiomes. Previous studies of Swiss Gruyère and selected Italian hard cheeses, mostly based on amplicon metagenomics, concurred that three species generally pre-dominate: Streptococcus thermophilus, Lactobacillus helveticus and Lactobacillus delbrueckii.

RESULTS

Two NWCs from Swiss Gruyère producers were subjected to whole metagenome shotgun sequencing using the Pacific Biosciences Sequel and Illumina MiSeq platforms. In addition, longer Oxford Nanopore Technologies MinION reads had to be generated for one to resolve repeat regions. Thereby, we achieved the complete assembly of all dominant bacterial genomes from these low-complexity NWCs, which was corroborated by a 16S rRNA amplicon survey. Moreover, two distinct L. helveticus strains were successfully co-assembled from the same sample. Besides bacterial chromosomes, we could also assemble several bacterial plasmids and phages and a corresponding prophage. Biologically relevant insights were uncovered by linking the plasmids and phages to their respective host genomes using DNA methylation motifs on the plasmids and by matching prokaryotic CRISPR spacers with the corresponding protospacers on the phages. These results could only be achieved by employing long-read sequencing data able to span intragenomic as well as intergenomic repeats.

CONCLUSIONS

Here, we demonstrate the feasibility of complete de novo genome assembly of all dominant strains from low-complexity NWCs based on whole metagenomics shotgun sequencing data. This allowed to gain novel biological insights and is a fundamental basis for subsequent systems-wide omics analyses, functional profiling and phenotype to genotype analysis of specific microbial communities.

virMine: automated detection of viral sequences from complex metagenomic samples

Metagenomics has enabled sequencing of viral communities from a myriad of different environments. Viral metagenomic studies routinely uncover sequences with no recognizable homology to known coding regions or genomes. Nevertheless, complete viral genomes have been constructed directly from complex community metagenomes, often through tedious manual curation. To address this, we developed the software tool virMine to identify viral genomes from raw reads representative of viral or mixed (viral and bacterial) communities. virMine automates sequence read quality control, assembly, and annotation. Researchers can easily refine their search for a specific study system and/or feature(s) of interest. In contrast to other viral genome detection tools that often rely on the recognition of viral signature sequences, virMine is not restricted by the insufficient representation of viral diversity in public data repositories. Rather, viral genomes are identified through an iterative approach, first omitting non-viral sequences. Thus, both relatives of previously characterized viruses and novel species can be detected, including both eukaryotic viruses and bacteriophages. Here we present virMine and its analysis of synthetic communities as well as metagenomic data sets from three distinctly different environments: the gut microbiota, the urinary microbiota, and freshwater viromes. Several new viral genomes were identified and annotated, thus contributing to our understanding of viral genetic diversity in these three environments.

Phages for Phage Therapy: Isolation, Characterization, and Host Range Breadth

For a bacteriophage to be useful for phage therapy it must be both isolated from the environment and shown to have certain characteristics beyond just killing strains of the target bacterial pathogen. These include desirable characteristics such as a relatively broad host range and a lack of other characteristics such as carrying toxin genes and the ability to form a lysogen. While phages are commonly isolated first and subsequently characterized, it is possible to alter isolation procedures to bias the isolation toward phages with desirable characteristics. Some of these variations are regularly used by some groups while others have only been shown in a few publications. In this review I will describe (1) isolation procedures and variations that are designed to isolate phages with broader host ranges, (2) characterization procedures used to show that a phage may have utility in phage therapy, including some of the limits of such characterization, and (3) results of a survey and discussion with phage researchers in industry and academia on the practice of characterization of phages.

Minimum Information about an Uncultivated Virus Genome (MIUViG)

We present an extension of the Minimum Information about any (x) Sequence (MIxS) standard for reporting sequences of uncultivated virus genomes. Minimum Information about an Uncultivated Virus Genome (MIUViG) standards were developed within the Genomic Standards Consortium framework and include virus origin, genome quality, genome annotation, taxonomic classification, biogeographic distribution and in silico host prediction. Community-wide adoption of MIUViG standards, which complement the Minimum Information about a Single Amplified Genome (MISAG) and Metagenome-Assembled Genome (MIMAG) standards for uncultivated bacteria and archaea, will improve the reporting of uncultivated virus genomes in public databases. In turn, this should enable more robust comparative studies and a systematic exploration of the global virosphere.

Characterization and Genomic Analyses of Pseudomonas aeruginosa Podovirus TC6: Establishment of Genus Pa11virus

Phages have attracted a renewed interest as alternative to chemical antibiotics. Although the number of phages is 10-fold higher than that of bacteria, the number of genomically characterized phages is far less than that of bacteria. In this study, phage TC6, a novel lytic virus of Pseudomonas aeruginosa, was isolated and characterized. TC6 consists of an icosahedral head with a diameter of approximately 54 nm and a short tail with a length of about 17 nm, which are characteristics of the family Podoviridae. TC6 can lyse 86 out of 233 clinically isolated P. aeruginosa strains, thus showing application potentials for phage therapy. The linear double-stranded genomic DNA of TC6 consisted of 49796 base pairs and was predicted to contain 71 protein-coding genes. A total of 11 TC6 structural proteins were identified by mass spectrometry. Comparative analysis revealed that the P. aeruginosa phages TC6, O4, PA11, and IME180 shared high similarity at DNA sequence and proteome levels, among which PA11 was the first phage discovered and published. Meanwhile, these phages contain 54 core genes and have very close phylogenetic relationships, which distinguish them from other known phage genera. We therefore proposed that these four phages can be classified as Pa11virus, comprising a new phage genus of Podoviridae that infects Pseudomonas spp. The results of this work promoted our understanding of phage biology, classification, and diversity.

Conserved Genome Organization and Core Transcriptome of the Lactobacillus acidophilus Complex

The Lactobacillus genus encompasses a genetically and functionally diverse group of species, and contains many strains widely formulated in the human food supply chain as probiotics and starter cultures. Within this genetically expansive group, there are several distinct clades that have high levels of homology, one of which is the Lactobacillus acidophilus group. Of the uniting features, small genomes, low GC content, adaptation to dairy environments, and fastidious growth requirements, are some of the most defining characteristics of this group. To better understand what truly links and defines this clade, we sought to characterize the genomic organization and content of the genomes of several members of this group. Through core genome analysis we explored the synteny and intrinsic genetic underpinnings of the L. acidophilus clade, and observed key features related to the evolution and adaptation of these organisms. While genetic content is able to provide a large map of the potential of each organism, it does not always reflect their functionality. Through transcriptomic data we inferred the core transcriptome of the L. acidophilus complex to better define the true metabolic capabilities that unite this clade. Using this approach we have identified seven small ORFs that are both highly conserved and transcribed in diverse members of this clade and could be potential novel small peptide or untranslated RNA regulators. Overall, our results reveal the core features of the L. acidophilus complex and open new avenues for the enhancement and formulation and of next generation probiotics and starter cultures.