Phage cluster relationships identified through single gene analysis

Isolation, characterization, and application of bacteriophage on Vibrio parahaemolyticus biofilm to control seafood contamination

OBJECTIVE

This study intended to isolate a Vibrio-particular phage from the natural environment, analyse its characteristics and genome sequence, and investigate its reduction effect on V. parahaemolyticus biofilm as a biocontrol agent in squid and mackerel.

METHODS

Among 21 phages, phage CAU_VPP01, isolated from beach mud, was chosen for further experiments based on host range and EOP tests. When examining the reduction effect of phage CAU_VPP01 against Vibrio parahaemolyticus biofilms on surfaces (stainless steel [SS] and polyethylene terephthalate [PET]) and food surfaces (squid and mackerel).

RESULTS

The phage showed the most excellent reduction effect at a multiplicity-of-infection (MOI) 10. Three-dimensional images acquired with confocal laser scanning microscopy (CLSM) analysis were quantified using COMSTAT, which showed that biomass, average thickness, and roughness coefficient decreased when treated with the phage. Colour and texture analysis confirmed that the quality of squid and mackerel was maintained after the phage treatment. Finally, a comparison of gene expression levels determined by qRT-PCR analysis showed that the phage treatment induced a decrease in the gene expression of flaA, vp0962, andluxS, as examples.

CONCLUSION

This study indicated that Vibrio-specific phage CAU_VPP01 effectively controlled V. parahaemolyticus biofilms under various conditions and confirmed that the isolated phage could possibly be used as an effective biocontrol weapon in the seafood manufacturing industry.

Insights into the genomic features and lifestyle of B1 subcluster mycobacteriophages

Bacteriophages infecting Mycobacterium smegmatis mc2155 are numerous and, hence, are classified into clusters based on nucleotide sequence similarity. Analyzing phages belonging to clusters/subclusters can help gain deeper insights into their biological features and potential therapeutic applications. In this study, for genomic characterization of B1 subcluster mycobacteriophages, a framework of online tools was developed, which enabled functional annotation of about 55% of the previously deemed hypothetical proteins in B1 phages. We also studied the phenotype, lysogeny status, and antimycobacterial activity of 10 B1 phages against biofilm and an antibiotic-resistant M. smegmatis strain (4XR1). All 10 phages belonged to the Siphoviridae family, appeared temperate based on their spontaneous release from the putative lysogens and showed antibiofilm activity. The highest inhibitory and disruptive effects on biofilm were 64% and 46%, respectively. This systematic characterization using a combination of genomic and experimental tools is a promising approach to furthering our understanding of viral dark matter.

Diversity and conservation of the genome architecture of phages infecting the Alphaproteobacteria

IMPORTANCE

This study reports the results of the largest analysis of genome sequences from phages that infect the Alphaproteobacteria class of bacterial hosts. We analyzed over 100 whole genome sequences of phages to construct dotplots, categorize them into genetically distinct clusters, generate a bootstrapped phylogenetic tree, compute protein orthologs, and predict packaging strategies. We determined that the phage sequences primarily cluster by the bacterial host family, phage morphotype, and genome size. We expect that the findings reported in this seminal study will facilitate future analyses that will improve our knowledge of the phages that infect these hosts.

Genomes of a Novel Group of Phages That Use Alternative Genetic Code Found in Human Gut Viromes

Metagenomics provides detection of phage genome sequences in various microbial communities. However, the use of alternative genetic codes by some phages precludes the correct analysis of their genomes. In this study, the unusual phage genome (phAss-1, 135,976 bp) was found after the de novo assembly of the human gut virome. Genome analysis revealed the presence of the TAG stop codons in 41 ORFs, including characteristic phage ORFs, and three genes of suppressor tRNA. Comparative analysis indicated that no phages with similar genomes were described. However, two phage genomes (BK046881_ctckW2 and BK025033_ct6IQ4) with substantial similarity to phAss-1 were extracted from the human gut metagenome data. These two complete genomes demonstrated 82.7% and 86.4% of nucleotide identity, respectively, similar genome synteny to phAss-1, the presence of suppressor tRNA genes and suppressor TAG stop codons in many characteristic phage ORFs. These data indicated that phAss-1, BK046881_ctckW2, and BK025033_ct6IQ4 are distinct species within the proposed Phassvirus genus. Moreover, a monophyletic group of divergent phage genomes containing the proposed Phassvirus genus was found among metagenome data. Several phage genomes from the group also contain ORFs with suppressor TAG stop codons, indicating the need to use various translation tables when depositing phage genomes in GenBank.

Identification, Characterization, and Genome Analysis of Two Novel Temperate Pseudomonas protegens Phages PseuP_222 and PseuP_224

Two novel P. protegens bacteriophages PseuP_222 and Pseu_224 and their host P. protegens CEMTC 4060 were isolated from the same sample (Inya river, Siberia). Both phages have siphovirus morphology and belong to lambdoid phages. Comparative genome analysis revealed a low nucleotide and amino acid sequence similarity of PseuP_222 and PseuP_224 between themselves, and between them and other lambdoid phages. Bioinformatics analysis indicated that PseuP_222 and PseuP_224 are members of a genetically diverse group of phages of environmental Pseudomonas spp.; this group is distant from a large group of P. aeruginosa phages. In phylogenetic trees, the positioning of the terminase large subunits, major capsid proteins, tail tape measure proteins, and CI-like repressors of PseuP_222 and PseuP_224 were remote and changed relative to those of the Escherichia lambda phage and lambdoid phages of Pseudomonas spp. However, the nucleoid-associated protein NdpA/YejK and P5-like structural protein from both phages showed high similarity and were not found in lambda phage and other lambdoid phages of Pseudomonas spp. Substantial divergences of the PseuP_222 and PseuP_224 genomes and proteomes indicated that the evolutionary history of these phages was mostly independent and they probably began to use one host only recently.

Characterization and genomic study of EJP2, a novel jumbo phage targeting antimicrobial resistant Escherichia coli

The emergence of antimicrobial resistance (AMR) Escherichia coli has noticeably increased in recent years worldwide and causes serious public health concerns. As alternatives to antibiotics, bacteriophages are regarded as promising antimicrobial agents. In this study, we isolated and characterized a novel jumbo phage EJP2 that specifically targets AMR E. coli strains. EJP2 belonged to the Myoviridae family with an icosahedral head (120.9 ± 2.9 nm) and a non-contractile tail (111.1 ± 0.6 nm), and contained 349,185 bp double-stranded DNA genome with 540 putative ORFs, suggesting that EJP2 could be classified as jumbo phage. The functions of genes identified in EJP2 genome were mainly related to nucleotide metabolism, DNA replication, and recombination. Comparative genomic analysis revealed that EJP2 was categorized in the group of Rak2-related virus and presented low sequence similarity at the nucleotide and amino acid level compared to other E. coli jumbo phages. EJP2 had a broad host spectrum against AMR E. coli as well as pathogenic E. coli and recognized LPS as a receptor for infection. Moreover, EJP2 treatment could remove over 80% of AMR E. coli biofilms on 96-well polystyrene, and exhibit synergistic antimicrobial activity with cefotaxime against AMR E. coli. These results suggest that jumbo phage EJP2 could be used as a potential biocontrol agent to combat the AMR issue in food processing and clinical environments.

Whole genome sequencing and annotation of a lysogenic phage vB_EcoP_DE5 isolated from donkey-derived Escherichia coli

A lysogenic phage vB_EcoP_DE5 (hereafter designated DE5) was isolated from donkey-derived Escherichia coli. The bacteriophage was examined by transmission electron microscopy, and the result showed that DE5 belonged to the genus Kuravirus. DE5 was sensitive to changes in temperature and pH, and it could maintain its activity at pH 7 and below 60 ℃. The whole genome sequencing revealed that DE5 had a double-stranded DNA genome of 77, 305 bp with 42.09% G+C content. A total of 126 open reading frames (ORFs) were identified, including functional genes related to phage integration, DNA replication and modification, transcriptional regulation, structural and packaging proteins, and host cell lysis. One phage integrase gene, one autotransporter adhesin gene, and one tRNA gene were predicted in the whole genome, and no genes associated with drug resistance were identified. The phage DE5 integrase contained 187 amino acids and belonged to the small serine recombinase family. BLASTn analysis revealed that phage DE5 had a high-sequence identity (96%) with E. coli phage SU10. Phylogenetic analysis showed that phage DE5 was a member of the genus Kuravirus. The whole genome sequencing of lysogenic phage DE5 enhanced our understanding of lysogenic phages and their therapeutic applications.

PCR Assay for Rapid Taxonomic Differentiation of Virulent Staphylococcus aureus and Klebsiella pneumoniae Bacteriophages

Phage therapy is now seen as a promising way to overcome the current global crisis in the spread of multidrug-resistant bacteria. However, phages are highly strain-specific, and in most cases one will have to isolate a new phage or search for a phage suitable for a therapeutic application in existing libraries. At an early stage of the isolation process, rapid screening techniques are needed to identify and type potential virulent phages. Here, we propose a simple PCR approach to differentiate between two families of virulent Staphylococcus phages (Herelleviridae and Rountreeviridae) and eleven genera of virulent Klebsiella phages (Przondovirus, Taipeivirus, Drulisvirus, Webervirus, Jiaodavirus, Sugarlandvirus, Slopekvirus, Jedunavirus, Marfavirus, Mydovirus and Yonseivirus). This assay includes a thorough search of a dataset comprising S. aureus (n = 269) and K. pneumoniae (n = 480) phage genomes available in the NCBI RefSeq/GenBank database for specific genes that are highly conserved at the taxonomic group level. The selected primers showed high sensitivity and specificity for both isolated DNA and crude phage lysates, which permits circumventing DNA purification protocols. Our approach can be extended and applied to any group of phages, given the large number of available genomes in the databases.

Pseudomonas Phage MD8: Genetic Mosaicism and Challenges of Taxonomic Classification of Lambdoid Bacteriophages

Pseudomonas phage MD8 is a temperate phage isolated from the freshwater lake Baikal. The organisation of the MD8 genome resembles the genomes of lambdoid bacteriophages. However, MD8 gene and protein sequences have little in common with classified representatives of lambda-like phages. Analysis of phage genomes revealed a group of other Pseudomonas phages related to phage MD8 and the genomic layout of MD8-like phages indicated extensive gene exchange involving even the most conservative proteins and leading to a high degree of genomic mosaicism. Multiple horizontal transfers and mosaicism of the genome of MD8, related phages and other λ-like phages raise questions about the principles of taxonomic classification of the representatives of this voluminous phage group. Comparison and analysis of various bioinformatic approaches applied to λ-like phage genomes demonstrated different efficiency and contradictory results in the estimation of genomic similarity and relatedness. However, we were able to make suggestions for the possible origin of the MD8 genome and the basic principles for the taxonomic classification of lambdoid phages. The group comprising 26 MD8-related phages was proposed to classify as two close genera belonging to a big family of λ-like phages.

Isolation, Genomic Analysis, and Preliminary Application of a Bovine Klebsiella pneumoniae Bacteriophage vB_Kpn_B01

Klebsiella pneumoniae is an important pathogen that can infect both humans and cattle. The widespread K. pneumoniae and its high drug resistance make it difficult to treat Klebsiella infections/diseases. In this study, a lytic K. pneumoniae bacteriophage vB_Kpn_B01 was isolated from a dairy farm trough in Sichuan Province, and its biological properties were studied, and the entire genome of vB_Kpn_B01 was sequenced. The therapeutic effects of the phage on disease-causing mice were preliminarily tested. Phages found in this study are double-stranded DNA bacterial viruses belonging to the family Siphoviridae, Sugarlandvirus. The results suggest that vB_Kpn_B01 has strong specificity and low adaptability to different adverse conditions. Meanwhile, the predicted gene products of phage vB_Kpn_B01 comprised 149 coding sequences (CDS) and 25 tRNAs, of which 34 CDS had known functions. Of course, vB_Kpn_B01 did not contain any known antibiotic-resistant or virulent genes. The pathological sections of the liver and lungs of mice showed that the inflammatory scores of the treatment group were lower than in the bacterial group. Phage vB_Kpn_B01 alleviated the inflammatory response in the organs of the infected mice, and the organ tissue bacterial load of the treatment group was significantly lower than that of the bacterial group. Therefore, vB_Kpn_B01 can inhibit the proliferation of K. pneumoniae 18 in vivo and can alleviate the inflammation of target organs caused by infectious bacteria, which preliminarily indicates that vB_Kpn_B01 has a certain therapeutic effect on laboratory-infected mice.

Novel Acinetobacter baumannii Bacteriophage Aristophanes Encoding Structural Polysaccharide Deacetylase

Acinetobacter baumannii appears to be one of the most crucial nosocomial pathogens. A possible component of antimicrobial therapy for infections caused by extremely drug-resistant A. baumannii strains may be specific lytic bacteriophages or phage-derived enzymes. In the present study, we observe the biological features, genomic organization, and phage–host interaction strategy of novel virulent bacteriophage Aristophanes isolated on A. baumannii strain having K26 capsular polysaccharide structure. According to phylogenetic analysis phage Aristophanes can be classified as a representative of a new distinct genus of the subfamily Beijerinckvirinae of the family Autographiviridae. This is the first reported A. baumannii phage carrying tailspike deacetylase, which caused O-acetylation of one of the K26 sugar residues.

Tailed Lytic Bacteriophages of Soft Rot Pectobacteriaceae

The study of the ecological and evolutionary traits of Soft Rot Pectobacteriaceae (SRP) comprising genera Pectobacterium and Dickeya often involves bacterial viruses (bacteriophages). Bacteriophages are considered to be a prospective tool for the ecologically safe and highly specific protection of plants and harvests from bacterial diseases. Information concerning bacteriophages has been growing rapidly in recent years, and this has included new genomics-based principles of taxonomic distribution. In this review, we summarise the data on phages infecting Pectobacterium and Dickeya that are available in publications and genomic databases. The analysis highlights not only major genomic properties that assign phages to taxonomic families and genera, but also the features that make them potentially suitable for phage control applications. Specifically, there is a discussion of the molecular mechanisms of receptor recognition by the phages and problems concerning the evolution of phage-resistant mutants.

Bi- and Multi-directional Gene Transfer in the Natural Populations of Polyvalent Bacteriophages, and Their Host Species Spectrum Representing Foodborne Versus Other Human and/or Animal Pathogens

Unraveling the trends of phage-host versus phage-phage coevolution is critical for avoiding possible undesirable outcomes from the use of phage preparations intended for therapeutic, food safety or environmental safety purposes. We aimed to investigate a phenomenon of intergeneric recombination and its trajectories across the natural populations of phages predominantly linked to foodborne pathogens. The results from the recombination analyses, using a large array of the recombination detection algorithms imbedded in SplitsTree, RDP4, and Simplot software packages, provided strong evidence (fit: 100; P  ≤ 0.014) for both bi- and multi-directional intergeneric recombination of the genetic loci involved collectively in phage morphogenesis, host specificity, virulence, replication, and persistence. Intergeneric recombination was determined to occur not only among conspecifics of the virulent versus temperate phages but also between the phages with these different lifestyles. The recombining polyvalent phages were suggested to interact with fairly large host species networks, including sometimes genetically very distinct species, such as e.g., Salmonella enterica and/or Escherichia coli versus Staphylococcus aureus or Yersinia pestis. Further studies are needed to understand whether phage-driven intergeneric recombination can lead to undesirable changes of intestinal and other microbiota in humans and animals.

Origin and Evolution of Studiervirinae Bacteriophages Infecting Pectobacterium: Horizontal Transfer Assists Adaptation to New Niches

Black leg and soft rot are devastating diseases causing up to 50% loss of potential potato yield. The search for, and characterization of, bacterial viruses (bacteriophages) suitable for the control of these diseases is currently a sought-after task for agricultural microbiology. Isolated lytic Pectobacterium bacteriophages Q19, PP47 and PP81 possess a similar broad host range but differ in their genomic properties. The genomic features of characterized phages have been described and compared to other Studiervirinae bacteriophages. Thorough phylogenetic analysis has clarified the taxonomy of the phages and their positioning relative to other genera of the Autographiviridae family. Pectobacterium phage Q19 seems to represent a new genus not described previously. The genomes of the phages are generally similar to the genome of phage T7 of the Teseptimavirus genus but possess a number of specific features. Examination of the structure of the genes and proteins of the phages, including the tail spike protein, underlines the important role of horizontal gene exchange in the evolution of these phages, assisting their adaptation to Pectobacterium hosts. The results provide the basis for the development of bacteriophage-based biocontrol of potato soft rot as an alternative to the use of antibiotics.

Pseudomonas Phage PaBG-A Jumbo Member of an Old Parasite Family

Bacteriophage PaBG is a jumbo Myoviridae phage isolated from water of Lake Baikal. This phage has limited diffusion ability and thermal stability and infects a narrow range of Pseudomonas aeruginosa strains. Therefore, it is hardly suitable for phage therapy applications. However, the analysis of the genome of PaBG presents a number of insights into the evolutionary history of this phage and jumbo phages in general. We suggest that PaBG represents an ancient group distantly related to all known classified families of phages.

A Unique Isolation of a Lytic Bacteriophage Infected Bacillus anthracis Isolate from Pafuri, South Africa

Bacillus anthracis is a soil-borne, Gram-positive endospore-forming bacterium and the causative agent of anthrax. It is enzootic in Pafuri, Kruger National Park in South Africa. The bacterium is amplified in a wild ungulate host, which then becomes a source of infection to the next host upon its death. The exact mechanisms involving the onset (index case) and termination of an outbreak are poorly understood, in part due to a paucity of information about the soil-based component of the bacterium’s lifecycle. In this study, we present the unique isolation of a dsDNA bacteriophage from a wildebeest carcass site suspected of having succumbed to anthrax. The aggressively lytic bacteriophage hampered the initial isolation of B. anthracis from samples collected at the carcass site. Classic bacteriologic methods were used to test the isolated phage on B. anthracis under different conditions to simulate deteriorating carcass conditions. Whole genome sequencing was employed to determine the relationship between the bacterium isolated on site and the bacteriophage-dubbed Bacillus phage Crookii. The 154,012 bp phage belongs to Myoviridae and groups closely with another African anthrax carcass-associated Bacillus phage WPh. Bacillus phage Crookii was lytic against B. cereus sensu lato group members but demonstrated a greater affinity for encapsulated B. anthracis at lower concentrations (<1 × 10⁸ pfu) of bacteriophage. The unusual isolation of this bacteriophage demonstrates the phage’s role in decreasing the inoculum in the environment and impact on the life cycle of B. anthracis at a carcass site.

Functional characterization of the endolysins derived from mycobacteriophage PDRPxv

Bacteriophage-derived endolysin enzymes play a critical role in disintegration of the host bacterial cell wall and hence have gained considerable attention as possible therapeutics for the treatment of drug-resistant infections. Endolysins can target both dividing and non-dividing cells and given the vital role peptidoglycan plays in bacterial survival, bacteria are less likely to modify it even if continuously exposed to lysins. Hence, probability of bacteria developing resistance to lysins appear bleak. Endolysins from mycobacteriophages offer great potential as alternative therapeutics for the drug-resistant TB. However, considering that a large number of mycobacteriophages have been discovered so far, the information on endolysins come from only a few mycobacteriophages. In this study, we report the structural and functional characterization of endolysins (LysinA and LysinB) encoded by mycobacteriophage PDRPxv which belongs to B1 sub cluster. On in silico analysis, we found LysinA to be a modular protein having peptidase domain at the N-terminal (104 aa), a central amidase domain (174 aa) and the peptidoglycan binding domain (62 aa) at the C-terminal. Additionally, 'H-X-H', which is a conserved motif and characteristic of peptidase domains, and the conserved residues His-His-Asp, which are characteristic of amidase domain were also observed. In LysinB enzyme, a single α/β hydrolase domain having a catalytic triad (Ser-Asp-His) and G-X-S-X-G motif, which are characteristic of the serine esterase enzymes were predicted to be present. Both the enzymes were purified as recombinant proteins and their antimycobacterial activity against M. smegmatis was demonstrated through turbidimetric experiments and biochemical assay. Interesting observation in this study is the secretory nature of LysinA evident by its periplasmic expression in E.coli, which might explain the ability of PDRPxv to lyse the bacterial host in the absence of transmembrane Holin protein.

Weirdo19ES is a novel singleton mycobacteriophage that selects for glycolipid deficient phage-resistant M. smegmatis mutants

The sequencing and bioinformatics analysis of bacteriophages infecting mycobacteria has yielded a large amount of information on their evolution, including that on their environmental propagation on other genera such as Gordonia, closely related to Mycobacterium. However, little is known on mycobacteriophages cell biology such as the nature of their receptor(s) or their replication cycle. As part of our on-going screening for novel mycobacteriophages, we herein report the isolation and genome bioinformatics analysis of Weirdo19ES, a singleton Siphoviridae temperate mycobacteriophage with a 70.19% GC content. Nucleotide and protein sequence comparison to actinobacteriophage databases revealed that Weirdo19ES shows low homology to Gordonia phage Ruthy and mycobacteriophages falling in clusters Q and G and to singleton DS6A.Weirdo19ES also displays uncommon features such as a very short Lysin A gene (with only one enzymatic domain) and two putative HNH endonucleases. Mycobacterium smegmatis mutants resistant to Weirdo19ES are cross- resistant to I3. In agreement with that phenotype, analysis of cell envelope of those mutants showed that Weirdo19ES shares receptors with the transducing mycobacteriophage I3.This singleton mycobacteriophage adds up to the uncommonness of local mycobacteriophages previously isolated by our group and helps understanding the nature of mycobacteriophage receptors.

Characterization and genome analysis of B1 sub-cluster mycobacteriophage PDRPxv

Mycobacteriophages are viruses specific to mycobacteria that have gained attention as alternative therapeutic strategies for treating antibiotic-resistant infections. Mycobacteriophages are highly diverse and have been grouped into 29 clusters, 71 sub-clusters and 10 singletons based on the genome sequence. Here, we annotate the genome of PDRPxv, a lytic mycobacteriophage isolated from New Delhi; it belongs to the Siphoviridae family as determined by transmission electron microscopy. This phage survives at higher temperatures (up to 55 °C) and in alkaline conditions (up to pH11). PDRPxv phage genome is 69,171 bp in length with 66.35 % GC content and encodes 107 putative open reading frames and belongs to the B1 sub-cluster. Genome annotation indicated that genes for DNA encapsidation, structural proteins, replication/transcription and lysis of the host are present in functional clusters. Structural proteins encoded by Gp10-Gp12, Gp18, Gp25 and Gp28-Gp33 were identified by mass spectrometry. Interestingly, no gene encoding a holin function was found. Single-step growth curve revealed that PDRPxv has an adsorption time of 45 min, a latency time of 135 min and an average burst size of 99 phage particles per infected cell. The short latency period and the large burst size mark the lytic nature of the PDRPxv phage, which could therefore be a promising therapeutic candidate against pathogenic Mycobacterium species.

Phage therapy as a renewed therapeutic approach to mycobacterial infections: a comprehensive review

Mycobacterial infections are considered to a serious challenge of medicine, and the emergence of MDR and XDR tuberculosis is a serious public health problem. Tuberculosis can cause high morbidity and mortality around the world, particularly in developing countries. The emergence of drug-resistant Mycobacterium infection following limited therapeutic technologies coupled with the serious worldwide tuberculosis epidemic has adversely affected control programs, thus necessitating the study of the role bacteriophages in the treatment of mycobacterial infection. Bacteriophages are viruses that are isolated from several ecological specimens and do not exert adverse effects on patients. Phage therapy can be considered as a significant alternative to antibiotics for treating MDR and XDR mycobacterial infections. The useful ability of bacteriophages to kill Mycobacterium spp has been explored by numerous research studies that have attempted to investigate the phage therapy as a novel therapeutic/diagnosis approach to mycobacterial infections. However, there are restricted data about phage therapy for treating mycobacterial infections. This review presents comprehensive data about phage therapy in the treatment of mycobacterial infection, specifically tuberculosis disease.

A Novel, Highly Related Jumbo Family of Bacteriophages That Were Isolated Against Erwinia

Erwinia amylovora is a plant pathogen from the Erwiniaceae family and a causative agent of the devastating agricultural disease fire blight. Here we characterize eight lytic bacteriophages of E. amylovora that we isolated from the Wasatch front (Utah, United States) that are highly similar to vB_EamM_Ea35-70 which was isolated in Ontario, Canada. With the genome size ranging from 271 to 275 kb, this is a novel jumbo family of bacteriophages. These jumbo bacteriophages were further characterized through genomic and proteomic comparison, mass spectrometry, host range and burst size. Their proteomes are highly unstudied, with over 200 putative proteins with no known homologs. The production of 27 of these putative proteins was confirmed by mass spectrometry analysis. These bacteriophages appear to be most similar to bacteriophages that infect Pseudomonas and Ralstonia rather than Enterobacteriales bacteria by protein similarity, however, we were only able to detect infection of Erwinia and the closely related strains of Pantoea.

Genomic comparison of 60 completely sequenced bacteriophages that infect Erwinia and/or Pantoea bacteria

Erwinia and Pantoea are closely related bacterial plant pathogens in the Gram negative Enterobacteriales order. Sixty tailed bacteriophages capable of infecting these pathogens have been completely sequenced by investigators around the world and are in the current databases, 30 of which were sequenced by our lab. These 60 were compared to 991 other Enterobacteriales bacteriophage genomes and found to be, on average, just over twice the overall average length. These Erwinia and Pantoea phages comprise 20 clusters based on nucleotide and protein sequences. Five clusters contain only phages that infect the Erwinia and Pantoea genera, the other 15 clusters are closely related to bacteriophages that infect other Enterobacteriales; however, within these clusters the Erwinia and Pantoea phages tend to be distinct, suggesting ecological niche may play a diversification role. The failure of many of their encoded proteins to have predicted functions highlights the need for further study of these phages.

A major-capsid-protein-based multiplex PCR assay for rapid identification of selected virulent bacteriophage types

Bacteriophages represent a promising alternative for controlling pathogenic bacteria. They are ubiquitous in the environment, and their isolation is usually simple and fast. However, not every phage is suitable for biocontrol applications. It must be virulent (i.e., strictly lytic), non-transducing, and safe. We have developed a method for identifying selected types of virulent phages at an early stage of the isolation process to simplify the search for suitable candidates. Using the major capsid protein (MCP) as a phylogenetic marker, we designed degenerate primers for the identification of Felix O1-, GJ1-, N4-, SP6-, T4-, T7-, and Vi1-like phages in multiplex PCR setups with single phage plaques as templates. Performance of the MCP PCR assay was evaluated with a set of 26 well-characterized phages. Neither false-positive nor false-negative results were obtained. In addition, 154 phages from enrichment cultures from various environmental samples were subjected to MCP PCR analysis. Eight of them, specific for Salmonella enterica, Escherichia coli, or Erwinia amylovora, belonged to one of the selected phage types. Their PCR-based identification was successfully confirmed by pulsed-field gel electrophoresis of the phage genomes, electron microscopy, and sequencing of the amplified mcp gene fragment. The MCP PCR assay was shown to be a simple method for preliminary assignment of new phages to a certain group and thus to identify candidates for biocontrol immediately after their isolation. Given that sufficient sequence data are available, this method can be extended to any phage group of interest.

Mycobacteriophage CRB2 defines a new subcluster in mycobacteriophage classification

Mycobacteriophages are viruses -mostly temperates- that infect Mycobacterium smegmatis and sometimes Mycobacterium tuberculosis. Mycobacteriophages are grouped in clusters on the basis of the overall nucleotide sequence homology, being further divided in subclusters as more mycobacteriophage genomes are sequenced and annotated. As part of our on-going screening for novel isolates, we herein report the bioinformatics analysis of CRB2, a mycobacteriophage belonging into the Siphoviridae family that propagates at 30°C. CRB2 has a 72,217 bp genome with a 69.78% GC content that belongs to Cluster B; nucleotide comparison with other B cluster members positions CRB2 as the sole member of a new subcluster, B9, being mycobacteriophage Saguaro (belonging into subcluster B7) its closest relative. Sequencing and annotation of 14 mycobacteriophages isolated by our group has yielded six cluster A members, a singleton, four of the five members of subcluster B6, one of the three reported members of subcluster G4, and CRB2 which defines subcluster B9. Considering the massive mycobacteriophage search performed in USA and the relatively rarity of our phages, we propose that factors other than size of the sampling determine the variability of mycobacteriophage distribution, and thus a world-wide concerted mining would most likely bring extremely rare and yet undiscovered mycobacteriophages.

Genomic characterization of three novel Basilisk-like phages infecting Bacillus anthracis

Background

In the present study, we sequenced the complete genomes of three novel bacteriophages v_B-Bak1, v_B-Bak6, v_B-Bak10 previously isolated from historical anthrax burial sites in the South Caucasus country of Georgia. We report here major trends in the molecular evolution of these phages, which we designate as “Basilisk-Like-Phages” (BLPs), and illustrate patterns in their evolution, genomic plasticity and core genome architecture.

Results

Comparative whole genome sequence analysis revealed a close evolutionary relationship between our phages and two unclassified Bacillus cereus group phages, phage Basilisk, a broad host range phage (Grose JH et al., J Vir. 2014;88(20):11846-11860) and phage PBC4, a highly host-restricted phage and close relative of Basilisk (Na H. et al. FEMS Microbiol. letters. 2016;363(12)). Genome comparisons of phages v_B-Bak1, v_B-Bak6, and v_B-Bak10 revealed significant similarity in sequence, gene content, and synteny with both Basilisk and PBC4. Transmission electron microscopy (TEM) confirmed the three phages belong to the Siphoviridae family. In contrast to the broad host range of phage Basilisk and the single-strain specificity of PBC4, our three phages displayed host specificity for Bacillus anthracis. Bacillus species including Bacillus cereus, Bacillus subtilis, Bacillus anthracoides, and Bacillus megaterium were refractory to infection.

Conclusions

Data reported here provide further insight into the shared genomic architecture, host range specificity, and molecular evolution of these rare B. cereus group phages. To date, the three phages represent the only known close relatives of the Basilisk and PBC4 phages and their shared genetic attributes and unique host specificity for B. anthracis provides additional insight into candidate host range determinants.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5056-4) contains supplementary material, which is available to authorized users.

Characterization of a lytic vibriophage VP06 of Vibrio parahaemolyticus

Vibrio parahaemolyticus is a human enteropathogenic bacterium and is also pathogenic to shrimp and finfish. In a search for a biocontrol agent for V. parahaemolyticus and other pathogenic Vibrio species, a lytic phage VP06 was isolated from oyster using V. parahaemolyticus as the host. VP06 is a Siphoviridae phage with a polyhedral head and a long tail. The genome sequence of VP06 was 75,893 nucleotides in length and the G + C content was 49%; a total of 101 CDSs were identified in VP06, of which 39 exhibited functional domains/motifs. The genomic sequence of VP06 is similar to those of a lytic Vibrio vulnificus phage SSP002 and a temperate V. parahaemolyticus phage vB_VpaS_MAR10, although VP06 has distinct features in the CDS arrangement and 14 unique CDSs. Phylogenetic analysis revealed that VP06, SSP002 and vB_VpaS_MAR10 belong to a novel genus cluster of Siphoviridae phages. This phage lysed 28.1% of various Vibrio strains, and the efficiency of plating method revealed that VP06 was highly effective in lysing strains of Vibrio alginolyticus, Vibrio azureus, Vibrio harveyi and V. parahaemolyticus. The properties of VP06, including its broad range of hosts and resistance to environmental stresses, indicate that it may be a candidate biocontrol agent.

Isolation and characterization of bacteriophages from India, with lytic activity against Mycobacterium tuberculosis

Bacteriophages are being considered as a promising natural resource for the development of alternative strategies against mycobacterial diseases, especially in the context of the wide-spread occurrence of drug resistance among the clinical isolates of Mycobacterium tuberculosis. However, there is not much information documented on mycobacteriophages from India. Here, we report the isolation of 17 mycobacteriophages using Mycobacterium smegmatis as the bacterial host, where 9 phages also lyse M. tuberculosis H37Rv. We present detailed analysis of one of these mycobacteriophages - PDRPv. Transmission electron microscopy and polymerase chain reaction analysis (of a conserved region within the TMP gene) show PDRPv to belong to the Siphoviridae family and B1 subcluster, respectively. The genome (69 110 bp) of PDRPv is circularly permuted double-stranded DNA with ∼66% GC content and has 106 open reading frames (ORFs). On the basis of sequence similarity and conserved domains, we have assigned function to 28 ORFs and have broadly categorized them into 6 groups that are related to replication and genome maintenance, DNA packaging, virion release, structural proteins, lysogeny-related genes and endolysins. The present study reports the occurrence of novel antimycobacterial phages in India and highlights their potential to contribute to our understanding of these phages and their gene products as potential antimicrobial agents.

Recovery of mycobacteriophages from archival stocks stored for approximately 50 years in Japan

Mycobacteriophage archival stocks have been kept for ca. 20-50 years in Japan. In this study, we attempted to recover mycobacteriophages from 50 archival stocks and briefly analyzed the recovered phages. The phages were recovered from 72.2% (13/18) of the lyophilized stocks that had been stored for 47-56 years. Moreover, the analysis of 12 representative recovered phages led to their classification as belonging to the family Siphoviridae, and seven of them were typed by polymerase chain reaction (PCR) targeting the gene that encodes the tape measure protein. Considering these results, lyophilization seems to be suitable for phage archival storage.

A Lytic Providencia rettgeri Virus of Potential Therapeutic Value Is a Deep-Branching Member of the T5virus Genus

Providencia rettgeri is emerging as a new opportunistic pathogen with high antibiotic resistance. The need to find alternative methods to control antibiotic-resistant bacteria and the recent advances in phage therapy motivate the search for new phages able to infect Providencia spp. This study describes the isolation and characterization of an obligatory lytic phage, vB_PreS_PR1 (PR1), with therapeutic potential against drug-resistant P. rettgeri PR1 is a siphovirus. Its virion DNA size (118,537 bp), transcriptional organization, terminal repeats (10,461 bp), and nicks in the 3'-to-5' strand are similar to those of phage T5. However, sequence similarities of PR1 to phages of the T5virus genus at the DNA and protein levels are limited, suggesting that it belongs to a new species within the Siphoviridae family. PR1 exhibits the ability to kill P. rettgeri antibiotic-resistant strains, is highly specific to the species, and did not present known genomic markers indicating a temperate lifestyle. The lack of homologies between its proteins and proteins of the only other sequenced Providencia prophage, Redjac, suggests that these two phages evolved separately and may target different host proteins.IMPORTANCE The alarming increase in the number of bacteria resistant to antibiotics has been observed worldwide. This is particularly true for Gram-negative bacteria. For certain of their strains, no effective antibiotics are available. Providencia sp. has been a neglected pathogen but is emerging as a multidrug-resistant bacterium. This has revived interest in bacteriophages as alternative therapeutic agents against this bacterium. We describe the morphological, physiological, and genomic characterization of a novel lytic virus, PR1, which is able to kill drug-resistant P. rettgeri clinical isolates. Genomic and phylogenetic analyses indicate that PR1 is a distant relative of T5virus genus representatives. The lack of known virulence- or temperate lifestyle-associated genes in the genome of PR1 makes this phage a potential candidate for therapeutic use. Analysis of its genome also improves our knowledge of the ecology and diversity of T5-like siphoviruses, providing a new link for evolutionary studies of this phage group.

Genome Sequences of Mycobacteriophages Findley, Hurricane, and TBond007

We report here the genome sequences of three newly isolated phages that infect Mycobacterium smegmatis mc²155. Phages Findley, Hurricane, and TBond007 were discovered in geographically distinct locations and are related to cluster K mycobacteriophages, with Findley being similar to subcluster K2 phages and Hurricane and TBond007 being similar to subcluster K3 phages.

Comparing Phylogenetic Trees by Matching Nodes Using the Transfer Distance Between Partitions

Ability to quantify dissimilarity of different phylogenetic trees describing the relationship between the same group of taxa is required in various types of phylogenetic studies. For example, such metrics are used to assess the quality of phylogeny construction methods, to define optimization criteria in supertree building algorithms, or to find horizontal gene transfer (HGT) events. Among the set of metrics described so far in the literature, the most commonly used seems to be the Robinson-Foulds distance. In this article, we define a new metric for rooted trees-the Matching Pair (MP) distance. The MP metric uses the concept of the minimum-weight perfect matching in a complete bipartite graph constructed from partitions of all pairs of leaves of the compared phylogenetic trees. We analyze the properties of the MP metric and present computational experiments showing its potential applicability in tasks related to finding the HGT events.

Functional and comparative genome analysis of novel virulent actinophages belonging to Streptomyces flavovirens

Background

Next Generation Sequencing (NGS) technologies provide exciting possibilities for whole genome sequencing of a plethora of organisms including bacterial strains and phages, with many possible applications in research and diagnostics. No Streptomyces flavovirens phages have been sequenced to date; there is therefore a lack in available information about S. flavovirens phage genomics. We report biological and physiochemical features and use NGS to provide the complete annotated genomes for two new strains (Sf1 and Sf3) of the virulent phage Streptomyces flavovirens, isolated from Egyptian soil samples.

Results

The S. flavovirens phages (Sf1 and Sf3) examined in this study show higher adsorption rates (82 and 85%, respectively) than other actinophages, indicating a strong specificity to their host, and latent periods (15 and 30 min.), followed by rise periods of 45 and 30 min. As expected for actinophages, their burst sizes were 1.95 and 2.49 virions per mL. Both phages were stable and, as reported in previous experiments, showed a significant increase in their activity after sodium chloride (NaCl) and magnesium chloride (MgCl₂.6H₂O) treatments, whereas after zinc chloride (ZnCl2) application both phages showed a significant decrease in infection.

The sequenced phage genomes are parts of a singleton cluster with sizes of 43,150 bp and 60,934 bp, respectively. Bioinformatics analyses and functional characterizations enabled the assignment of possible functions to 19 and 28 putative identified ORFs, which included phage structural proteins, lysis components and metabolic proteins.

Thirty phams were identified in both phages, 10 (33.3%) of them with known function, which can be used in cluster prediction. Comparative genomic analysis revealed significant homology between the two phages, showing the highest hits among Sf1, Sf3 and the closest Streptomyces phage (VWB phages) in a specific 13Kb region. However, the phylogenetic analysis using the Major Capsid Protein (MCP) sequences highlighted that the isolated phages belong to the BG Streptomyces phage group but are clearly separated, representing a novel sub-cluster.

Conclusion

The results of this study provide the first physiological and genomic information for S. flavovirens phages and will be useful for pharmaceutical industries based on S. flavovirens and future phage evolution studies.

Characterization of Five Novel Brevibacillus Bacteriophages and Genomic Comparison of Brevibacillus Phages

Brevibacillus laterosporus is a spore-forming bacterium that causes a secondary infection in beehives following European Foulbrood disease. To better understand the contributions of Brevibacillus bacteriophages to the evolution of their hosts, five novel phages (Jenst, Osiris, Powder, SecTim467, and Sundance) were isolated and characterized. When compared with the five Brevibacillus phages currently in NCBI, these phages were assigned to clusters based on whole genome and proteome synteny. Powder and Osiris, both myoviruses, were assigned to the previously described Jimmer-like cluster. SecTim467 and Jenst, both siphoviruses, formed a novel phage cluster. Sundance, a siphovirus, was assigned as a singleton phage along with the previously isolated singleton, Emery. In addition to characterizing the basic relationships between these phages, several genomic features were observed. A motif repeated throughout phages Jenst and SecTim467 was frequently upstream of genes predicted to function in DNA replication, nucleotide metabolism, and transcription, suggesting transcriptional co-regulation. In addition, paralogous gene pairs that encode a putative transcriptional regulator were identified in four Brevibacillus phages. These paralogs likely evolved to bind different DNA sequences due to variation at amino acid residues predicted to bind specific nucleotides. Finally, a putative transposable element was identified in SecTim467 and Sundance that carries genes homologous to those found in Brevibacillus chromosomes. Remnants of this transposable element were also identified in phage Jenst. These discoveries provide a greater understanding of the diversity of phages, their behavior, and their evolutionary relationships to one another and to their host. In addition, they provide a foundation with which further Brevibacillus phages can be compared.

Genome characterization of a novel vibriophage VpKK5 (Siphoviridae) specific to fish pathogenic strain of Vibrio parahaemolyticus

Vibrio parahaemolyticus has long been known pathogenic to shrimp but only recently it is also reported pathogenic to tropical cultured marine finfish. Traditionally, bacterial diseases in aquaculture are often treated using synthetic antibiotics but concern due to side effects of these chemicals is elevating hence, new control strategies which are both environmental and consumer friendly, are urgently needed. One promising control strategy is the bacteriophage therapy. In this study, we report the isolation and characterization of a novel vibriophage (VpKK5), belonging to the family Siphoviridae that was specific and capable of complete lysing the fish pathogenic strain of V. parahaemolyticus. The VpKK5 exhibited short eclipse and latent periods of 24 and 36 min, respectively, but with a large burst size of 180 pfu/cell. The genome analysis revealed that the VpKK5 is a novel bacteriophage with the estimated genome size of 56,637 bp and has 53.1% G + C content. The vibriophage has about 80 predicted open reading frames consisted of 37 complete coding sequences which did not match to any protein databases. The analysis also found no lysogeny and virulence genes in the genome of VpKK5. With such genome features, we suspected the vibriophage is novel and could be explored for phage therapy against fish pathogenic strains of V. parahaemolyticus in the near future.

Genome sequence of a cluster A13 mycobacteriophage detected in Mycobacterium phlei over a half century ago

A phage infecting Mycobacterium phlei was isolated in 1958 from a soil sample in Hungary. Some physicochemical and biological properties of the virus were described in independent studies over the years. Here, we report the genome sequence of this early mycobacteriophage isolate. The Phlei phage genome measured 50,418 bp, had a GC content of 60.1 % and was predicted to encode 81 proteins and three tRNAs. Phylogeny of the tape measure protein revealed genetic relatedness to other early isolates of mycobacteriophages within subcluster A2. The genomic organization and genetic relationships to other strains showed that the Phlei phage belongs to a novel genetic cluster, designated A13.

Understanding the enormous diversity of bacteriophages: the tailed phages that infect the bacterial family Enterobacteriaceae

Bacteriophages are the predominant biological entity on the planet. The recent explosion of sequence information has made estimates of their diversity possible. We describe the genomic comparison of 337 fully sequenced tailed phages isolated on 18 genera and 31 species of bacteria in the Enterobacteriaceae. These phages were largely unambiguously grouped into 56 diverse clusters (32 lytic and 24 temperate) that have syntenic similarity over >50% of the genomes within each cluster, but substantially less sequence similarity between clusters. Most clusters naturally break into sets of more closely related subclusters, 78% of which are correlated with their host genera. The largest groups of related phages are superclusters united by genome synteny to lambda (81 phages) and T7 (51 phages). This study forms a robust framework for understanding diversity and evolutionary relationships of existing tailed phages, for relating newly discovered phages and for determining host/phage relationships.

Correction: genomic comparison of 93 Bacillus phages reveals 12 clusters, 14 singletons and remarkable diversity

BACKGROUND

The Bacillus genus of Firmicutes bacteria is ubiquitous in nature and includes one of the best characterized model organisms, B. subtilis, as well as medically significant human pathogens, the most notorious being B. anthracis and B. cereus. As the most abundant living entities on the planet, bacteriophages are known to heavily influence the ecology and evolution of their hosts, including providing virulence factors. Thus, the identification and analysis of Bacillus phages is critical to understanding the evolution of Bacillus species, including pathogenic strains.

RESULTS

Whole genome nucleotide and proteome comparison of the 83 extant, fully sequenced Bacillus phages revealed 10 distinct clusters, 24 subclusters and 15 singleton phages. Host analysis of these clusters supports host boundaries at the subcluster level and suggests phages as vectors for genetic transfer within the Bacillus cereus group, with B. anthracis as a distant member. Analysis of the proteins conserved among these phages reveals enormous diversity and the uncharacterized nature of these phages, with a total of 4,442 protein families (phams) of which only 894 (20%) had a predicted function. In addition, 2,583 (58%) of phams were orphams (phams containing a single member). The most populated phams were those encoding proteins involved in DNA metabolism, virion structure and assembly, cell lysis, or host function. These included several genes that may contribute to the pathogenicity of Bacillus strains.

CONCLUSIONS

This analysis provides a basis for understanding and characterizing Bacillus and other related phages as well as their contributions to the evolution and pathogenicity of Bacillus cereus group bacteria. The presence of sparsely populated clusters, the high ratio of singletons to clusters, and the large number of uncharacterized, conserved proteins confirms the need for more Bacillus phage isolation in order to understand the full extent of their diversity as well as their impact on host evolution.

Automated classification of tailed bacteriophages according to their neck organization

Background

The genetic diversity observed among bacteriophages remains a major obstacle for the identification of homologs and the comparison of their functional modules. In the structural module, although several classes of homologous proteins contributing to the head and tail structure can be detected, proteins of the head-to-tail connection (or neck) are generally more divergent. Yet, molecular analyses of a few tailed phages belonging to different morphological classes suggested that only a limited number of structural solutions are used in order to produce a functional virion. To challenge this hypothesis and analyze proteins diversity at the virion neck, we developed a specific computational strategy to cope with sequence divergence in phage proteins. We searched for homologs of a set of proteins encoded in the structural module using a phage learning database.

Results

We show that using a combination of iterative profile-profile comparison and gene context analyses, we can identify a set of head, neck and tail proteins in most tailed bacteriophages of our database. Classification of phages based on neck protein sequences delineates 4 Types corresponding to known morphological subfamilies. Further analysis of the most abundant Type 1 yields 10 Clusters characterized by consistent sets of head, neck and tail proteins. We developed Virfam, a webserver that automatically identifies proteins of the phage head-neck-tail module and assign phages to the most closely related cluster of phages. This server was tested against 624 new phages from the NCBI database. 93% of the tailed and unclassified phages could be assigned to our head-neck-tail based categories, thus highlighting the large representativeness of the identified virion architectures. Types and Clusters delineate consistent subgroups of Caudovirales, which correlate with several virion properties.

Conclusions

Our method and webserver have the capacity to automatically classify most tailed phages, detect their structural module, assign a function to a set of their head, neck and tail genes, provide their morphologic subtype and localize these phages within a “head-neck-tail” based classification. It should enable analysis of large sets of phage genomes. In particular, it should contribute to the classification of the abundant unknown viruses found on assembled contigs of metagenomic samples.

Electronic supplementary material

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

PHIRE and TWiV: Experiences in Bringing Virology to New Audiences

Virology encompasses a broad spectrum of topics touching upon many aspects of our everyday lives. However, appreciation of this impact is too often restricted to those who have specialized training and participate in virology research. The Phage Hunters Integrating Research and Education (PHIRE) program and the This Week in Virology (TWiV) podcast seek to bring virology to new audiences through two different approaches—direct involvement of undergraduates in discovering and genomically characterizing bacteriophages (PHIRE) and clear, accessible, and free discussions among experts of all topics in virology (TWiV). Here we discuss these two high-impact programs, the audiences that they serve, their broader impacts, and their future potential.

Isolation and characterization of vB_ArS-ArV2 - first Arthrobacter sp. infecting bacteriophage with completely sequenced genome

This is the first report on a complete genome sequence and biological characterization of the phage that infects Arthrobacter. A novel virus vB_ArS-ArV2 (ArV2) was isolated from soil using Arthrobacter sp. 68b strain for phage propagation. Based on transmission electron microscopy, ArV2 belongs to the family Siphoviridae and has an isometric head (∼63 nm in diameter) with a non-contractile flexible tail (∼194×10 nm) and six short tail fibers. ArV2 possesses a linear, double-stranded DNA genome (37,372 bp) with a G+C content of 62.73%. The genome contains 68 ORFs yet encodes no tRNA genes. A total of 28 ArV2 ORFs have no known functions and lack any reliable database matches. Proteomic analysis led to the experimental identification of 14 virion proteins, including 9 that were predicted by bioinformatics approaches. Comparative phylogenetic analysis, based on the amino acid sequence alignment of conserved proteins, set ArV2 apart from other siphoviruses. The data presented here will help to advance our understanding of Arthrobacter phage population and will extend our knowledge about the interaction between this particular host and its phages.

Genomic comparison of 93 Bacillus phages reveals 12 clusters, 14 singletons and remarkable diversity

Background

The Bacillus genus of Firmicutes bacteria is ubiquitous in nature and includes one of the best characterized model organisms, B. subtilis, as well as medically significant human pathogens, the most notorious being B. anthracis and B. cereus. As the most abundant living entities on the planet, bacteriophages are known to heavily influence the ecology and evolution of their hosts, including providing virulence factors. Thus, the identification and analysis of Bacillus phages is critical to understanding the evolution of Bacillus species, including pathogenic strains.

Results

Whole genome nucleotide and proteome comparison of the 93 extant Bacillus phages revealed 12 distinct clusters, 28 subclusters and 14 singleton phages. Host analysis of these clusters supports host boundaries at the subcluster level and suggests phages as vectors for genetic transfer within the Bacillus cereus group, with B. anthracis as a distant member of the group. Analysis of the proteins conserved among these phages reveals enormous diversity and the uncharacterized nature of these phages, with a total of 4,922 protein families (phams) of which only 951 (19%) had a predicted function. In addition, 3,058 (62%) of phams were orphams (phams containing a gene product from a single phage). The most populated phams were those encoding proteins involved in DNA metabolism, virion structure and assembly, cell lysis, or host function. These included several genes that may contribute to the pathogenicity of Bacillus strains.

Conclusions

This analysis provides a basis for understanding and characterizing Bacillus phages and other related phages as well as their contributions to the evolution and pathogenicity of Bacillus cereus group bacteria. The presence of sparsely populated clusters, the high ratio of singletons to clusters, and the large number of uncharacterized, conserved proteins confirms the need for more Bacillus phage isolation in order to understand the full extent of their diversity as well as their impact on host evolution.

Mycobacteriophage-drived diversification of Mycobacterium abscessus

Background

Mycobacterium abscessus is an emerging opportunistic pathogen which diversity was acknowledged by the recent description of two subspecies accommodating M. abscessus, Mycobacterium bolletii and Mycobacterium massiliense isolates.

Results

Here, genome analysis found 1–8 prophage regions in 47/48 M. abscessus genomes ranging from small prophage-like elements to complete prophages. A total of 20,304 viral and phage proteins clustered into 853 orthologous groups. Phylogenomic and phylogenetic analyses based on prophage region homology found three main clusters corresponding to M. abscessus, M. bolletii and M. massiliense. Analysing 135 annotated Tape Measure Proteins found thirteen clusters and four singletons, suggesting that at least 17 mycobacteriophages had infected M. abscessus during its evolution. The evolutionary history of phages differed from that of their mycobacterial hosts. In particular, 33 phage-related proteins have been horizontally transferred within M. abscessus genomes. They comprise of an integrase, specific mycobacteriophage proteins, hypothetical proteins and DNA replication and metabolism proteins. Gene exchanges, loss and gains which occurred in M. abscessus genomes have been driven by several mycobacteriophages.

Conclusions

This analysis of phage-mycobacterium co-evolution suggests that mycobacteriophages are playing a key-role in the on-going diversification of M. abscessus.

Reviewers

This article was reviewed by Eric Bapteste, Patrick Forterre and Eugene Koonin.

Characterization of Paenibacillus larvae bacteriophages and their genomic relationships to firmicute bacteriophages

Background

Paenibacillus larvae is a Firmicute bacterium that causes American Foulbrood, a lethal disease in honeybees and is a major source of global agricultural losses. Although P. larvae phages were isolated prior to 2013, no full genome sequences of P. larvae bacteriophages were published or analyzed. This report includes an in-depth analysis of the structure, genomes, and relatedness of P. larvae myoviruses Abouo, Davis, Emery, Jimmer1, Jimmer2, and siphovirus phiIBB_Pl23 to each other and to other known phages.

Results

P. larvae phages Abouo, Davies, Emery, Jimmer1, and Jimmer2 are myoviruses with ~50 kbp genomes. The six P. larvae phages form three distinct groups by dotplot analysis. An annotated linear genome map of these six phages displays important identifiable genes and demonstrates the relationship between phages. Sixty phage assembly or structural protein genes and 133 regulatory or other non-structural protein genes were identifiable among the six P. larvae phages. Jimmer1, Jimmer2, and Davies formed stable lysogens resistant to superinfection by genetically similar phages. The correlation between tape measure protein gene length and phage tail length allowed identification of co-isolated phages Emery and Abouo in electron micrographs. A Phamerator database was assembled with the P. larvae phage genomes and 107 genomes of Firmicute-infecting phages, including 71 Bacillus phages. Phamerator identified conserved domains in 1,501 of 6,181 phamilies (only 24.3%) encoded by genes in the database and revealed that P. larvae phage genomes shared at least one phamily with 72 of the 107 other phages. The phamily relationship of large terminase proteins was used to indicate putative DNA packaging strategies. Analyses from CoreGenes, Phamerator, and electron micrograph measurements indicated Jimmer1, Jimmer2, Abouo and Davies were related to phages phiC2, EJ-1, KC5a, and AQ113, which are small-genome myoviruses that infect Streptococcus, Lactobacillus, and Clostridium, respectively.

Conclusions

This paper represents the first comparison of phage genomes in the Paenibacillus genus and the first organization of P. larvae phages based on sequence and structure. This analysis provides an important contribution to the field of bacteriophage genomics by serving as a foundation on which to build an understanding of the natural predators of P. larvae.

Electronic supplementary material

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

A web-based restriction endonuclease tool for mycobacteriophage cluster prediction

A recent explosion in the amount of genomic data has revealed a large genetic diversity in the bacteriophages that infect Mycobacterium smegmatis. In an effort to assess the novelty of newly described mycobacteriophage isolates and provide a preliminary determination of their probable cluster assignment prior to full genome sequencing, we have developed a systematic approach that relies on restriction endonuclease analysis. We demonstrate that a web-based tool, the Phage Enzyme Tool (or PET), is capable of rapidly facilitating this analysis and exhibits reliability in the putative placement of mycobacteriophages into specific clusters of previously sequenced phages. We propose that this tool represents a useful analytical step in the initial study of phage genomes and that this tool will increase the efficiency of phage genome characterization and enhance the educational activities involving mycobacteriophage discovery.

A broadly implementable research course in phage discovery and genomics for first-year undergraduate students

Engaging large numbers of undergraduates in authentic scientific discovery is desirable but difficult to achieve. We have developed a general model in which faculty and teaching assistants from diverse academic institutions are trained to teach a research course for first-year undergraduate students focused on bacteriophage discovery and genomics. The course is situated within a broader scientific context aimed at understanding viral diversity, such that faculty and students are collaborators with established researchers in the field. The Howard Hughes Medical Institute (HHMI) Science Education Alliance Phage Hunters Advancing Genomics and Evolutionary Science (SEA-PHAGES) course has been widely implemented and has been taken by over 4,800 students at 73 institutions. We show here that this alliance-sourced model not only substantially advances the field of phage genomics but also stimulates students' interest in science, positively influences academic achievement, and enhances persistence in science, technology, engineering, and mathematics (STEM) disciplines. Broad application of this model by integrating other research areas with large numbers of early-career undergraduate students has the potential to be transformative in science education and research training.

IMPORTANCE

Engagement of undergraduate students in scientific research at early stages in their careers presents an opportunity to excite students about science, technology, engineering, and mathematics (STEM) disciplines and promote continued interests in these areas. Many excellent course-based undergraduate research experiences have been developed, but scaling these to a broader impact with larger numbers of students is challenging. The Howard Hughes Medical Institute (HHMI) Science Education Alliance Phage Hunting Advancing Genomics and Evolutionary Science (SEA-PHAGES) program takes advantage of the huge size and diversity of the bacteriophage population to engage students in discovery of new viruses, genome annotation, and comparative genomics, with strong impacts on bacteriophage research, increased persistence in STEM fields, and student self-identification with learning gains, motivation, attitude, and career aspirations.

A broadly implementable research course in phage discovery and genomics for first-year undergraduate students

Engaging large numbers of undergraduates in authentic scientific discovery is desirable but difficult to achieve. We have developed a general model in which faculty and teaching assistants from diverse academic institutions are trained to teach a research course for first-year undergraduate students focused on bacteriophage discovery and genomics. The course is situated within a broader scientific context aimed at understanding viral diversity, such that faculty and students are collaborators with established researchers in the field. The Howard Hughes Medical Institute (HHMI) Science Education Alliance Phage Hunters Advancing Genomics and Evolutionary Science (SEA-PHAGES) course has been widely implemented and has been taken by over 4,800 students at 73 institutions. We show here that this alliance-sourced model not only substantially advances the field of phage genomics but also stimulates students' interest in science, positively influences academic achievement, and enhances persistence in science, technology, engineering, and mathematics (STEM) disciplines. Broad application of this model by integrating other research areas with large numbers of early-career undergraduate students has the potential to be transformative in science education and research training.

IMPORTANCE

Engagement of undergraduate students in scientific research at early stages in their careers presents an opportunity to excite students about science, technology, engineering, and mathematics (STEM) disciplines and promote continued interests in these areas. Many excellent course-based undergraduate research experiences have been developed, but scaling these to a broader impact with larger numbers of students is challenging. The Howard Hughes Medical Institute (HHMI) Science Education Alliance Phage Hunting Advancing Genomics and Evolutionary Science (SEA-PHAGES) program takes advantage of the huge size and diversity of the bacteriophage population to engage students in discovery of new viruses, genome annotation, and comparative genomics, with strong impacts on bacteriophage research, increased persistence in STEM fields, and student self-identification with learning gains, motivation, attitude, and career aspirations.