Genome sequence and analysis of a broad-host range lytic bacteriophage that infects the Bacillus cereus group

Isolation of Bacteriophages from Soil Samples in a Poorly Equipped Field Laboratory in Kruger National Park

Bacteriophages are viruses that infect bacteria. Bacteriophages are ubiquitous and are the most abundant organisms on the planet. Despite this, very little is known about the influence and effect of bacteriophages within terrestrial environments. Additionally, the natural soil microbiome profiles remain largely unexplored. Here we describe protocols that can be used, in field or rural laboratories containing only basic equipment, to make bacteriophage isolation more accessible and to facilitate such research.

Phage-Encoded Sigma Factors Alter Bacterial Dormancy

By entering a reversible state of reduced metabolic activity, dormant microorganisms are able to tolerate suboptimal conditions that would otherwise reduce their fitness. Dormancy may also benefit bacteria by serving as a refuge from parasitic infections. Here, we focus on dormancy in the Bacillota, where endospore development is transcriptionally regulated by the expression of sigma factors. A disruption of this process could influence the survivorship or reproduction of phages that infect spore-forming hosts with implications for coevolutionary dynamics. We characterized the distribution of sigma factors in over 4,000 genomes of diverse phages capable of infecting hosts that span the bacterial domain. From this, we identified homologs of sporulation-specific sigma factors in phages that infect spore-forming hosts. Unlike sigma factors required for phage reproduction, we provide evidence that sporulation-like sigma factors are nonessential for lytic infection. However, when expressed in the spore-forming Bacillus subtilis, some of these phage-derived sigma factors can activate the bacterial sporulation gene network and lead to a reduction in spore yield. Our findings suggest that the acquisition of host-like transcriptional regulators may allow phages to manipulate a complex and ancient trait in one of the most abundant cell types on Earth. IMPORTANCE As obligate parasites, phages exert strong top-down pressure on host populations with eco-evolutionary implications for community dynamics and ecosystem functioning. The process of phage infection, however, is constrained by bottom-up processes that influence the energetic and nutritional status of susceptible hosts. Many phages have acquired auxiliary genes from bacteria, which can be used to exploit host metabolism with consequences for phage fitness. In this study, we demonstrate that phages infecting spore-forming bacteria carry homologs of sigma factors, which their hosts use to orchestrate gene expression during spore development. By tapping into regulatory gene networks, phages may manipulate the physiology and survival strategies of nongrowing bacteria in ways that influence host-parasite coevolution.

Phenotypic and genotypic characterization of the new Bacillus cereus phage SWEP1

A new Bacillus cereus phage, SWEP1, was isolated from black soil. The host lysis activity of phage SWEP1 has a relatively short latent time (20 min) and a small burst size of 83 PFU. The genome of SWEP1 consists of 162,461 bp with 37.77% G+C content. The phage encodes 278 predicted proteins, 103 of which were assigned functionally. No tRNA genes were found. Comparative genomics analysis indicated that SWEP1 is related to Bacillus phage B4 (86.91% identity, 90% query coverage). Phenotypic and genotypic characterization suggested that SWEP1 is a new member of a new species in the genus Bequatrovirus, family Herelleviridae.

Isolation, Characterization, and Genome Sequence Analysis of a Novel Lytic Phage, Xoo-sp15 Infecting Xanthomonas oryzae pv. oryzae

Xanthomonas oryzae pv. oryzae (X. oryzae) is a bacterial pathovar of rice diseases all over the world. Owing to emerging antibacterial resistance, phage therapies have gained significant attention to treat various bacterial infections. Nevertheless, comprehensive research is needed for their use as a safe biocontrol agent. In this study, isolation and characterization of a novel phage Xoo-sp15, that infects X. oryzae was ascertained through experimental and bioinformatics analyses to determine its virulent potency and reliability. High throughput sequencing demonstrated that Xoo-sp15 has a dsDNA genome with a total size of 157,091 bp and 39.9% GC content lower than its host (63.6%). Morphological and phylogenetic analyses characterized it as a new member of the Bastille-like group within the family Herelleviridae. In silico analysis revealed that it contains 229 open reading frames and 16 tRNAs. Additionally, this novel phage does not contain any resistant determinants and can infect nine X. oryzae strains. Therefore, Xoo-sp15 has the potential to serve as a novel candidate for phage therapy.

Expanding Diversity of Firmicutes Single-Strand Annealing Proteins: A Putative Role of Bacteriophage-Host Arms Race

Bacteriophage-encoded single strand annealing proteins (SSAPs) are recombinases which can substitute the classical, bacterial RecA and manage the DNA metabolism at different steps of phage propagation. SSAPs have been shown to efficiently promote recombination between short and rather divergent DNA sequences and were exploited for in vivo genetic engineering mainly in Gram-negative bacteria. In opposition to the conserved and almost universal bacterial RecA protein, SSAPs display great sequence diversity. The importance for SSAPs in phage biology and phage-bacteria evolution is underlined by their role as key players in events of horizontal gene transfer (HGT). All of the above provoke a constant interest for the identification and study of new phage recombinase proteins in vivo, in vitro as well as in silico. Despite this, a huge body of putative ssap genes escapes conventional classification, as they are not properly annotated. In this work, we performed a wide-scale identification, classification and analysis of SSAPs encoded by the Firmicutes bacteria and their phages. By using sequence similarity network and gene context analyses, we created a new high quality dataset of phage-related SSAPs, substantially increasing the number of annotated SSAPs. We classified the identified SSAPs into seven distinct families, namely RecA, Gp2.5, RecT/Redβ, Erf, Rad52/22, Sak3, and Sak4, organized into three superfamilies. Analysis of the relationships between the revealed protein clusters led us to recognize Sak3-like proteins as a new distinct SSAP family. Our analysis showed an irregular phylogenetic distribution of ssap genes among different bacterial phyla and specific phages, which can be explained by the high rates of ssap HGT. We propose that the evolution of phage recombinases could be tightly linked to the dissemination of bacterial phage-resistance mechanisms (e.g., abortive infection and CRISPR/Cas systems) targeting ssap genes and be a part of the constant phage-bacteria arms race.

Bacteria have numerous distinctive groups of phage-plasmids with conserved phage and variable plasmid gene repertoires

Plasmids and temperate phages are key contributors to bacterial evolution. They are usually regarded as very distinct. However, some elements, termed phage–plasmids, are known to be both plasmids and phages, e.g. P1, N15 or SSU5. The number, distribution, relatedness and characteristics of these phage–plasmids are poorly known. Here, we screened for these elements among ca. 2500 phages and 12000 plasmids and identified 780 phage–plasmids across very diverse bacterial phyla. We grouped 92% of them by similarity of gene repertoires to eight defined groups and 18 other broader communities of elements. The existence of these large groups suggests that phage–plasmids are ancient. Their gene repertoires are large, the average element is larger than an average phage or plasmid, and they include slightly more homologs to phages than to plasmids. We analyzed the pangenomes and the genetic organization of each group of phage–plasmids and found the key phage genes to be conserved and co-localized within distinct groups, whereas genes with homologs in plasmids are much more variable and include most accessory genes. Phage–plasmids are a sizeable fraction of the sequenced plasmids (∼7%) and phages (∼5%), and could have key roles in bridging the genetic divide between phages and other mobile genetic elements.

Isolation and Characterization of Bacillus cereus Phage vB_BceP-DLc1 Reveals the Largest Member of the Φ29-Like Phages

Bacillus phage φ29 and its relatives have been considered as one of the most important model organisms for DNA replication, transcription, morphogenesis, DNA packaging studies, and nanotechnology applications. Here, we isolated and characterized a new member of the φ29-like phage, named Bacillus cereus phage vB_BceP-DLc1. This phage, with a unique inserted gene cluster, has the largest genome among known φ29-like phages. DLc1 can use the surface carbohydrate structures of the host cell as receptors and only infects the most related B. cereus strains, showing high host-specificity. The adsorption rate constant and life cycle of DLc1 under experimental conditions were also determined. Not only stable under temperatures below 55 °C and pH range from 5 to 11, the new phage also showed tolerance to high concentrations of NaCl, 75% ethanol, chloroform, and mechanical vortex, which is preferable for practical use in the food and pharmaceutical industries.

Genomic analyses of two novel biofilm-degrading methicillin-resistant Staphylococcus aureus phages

Background

Methicillin-resistant Staphylococcus aureus (MRSA) biofilm producers represent an important etiological agent of many chronic human infections. Antibiotics and host immune responses are largely ineffective against bacteria within biofilms. Alternative actions and novel antimicrobials should be considered. In this context, the use of phages to destroy MRSA biofilms presents an innovative alternative mechanism.

Results

Twenty-five MRSA biofilm producers were used as substrates to isolate MRSA-specific phages. Despite the difficulties in obtaining an isolate of this phage, two phages (UPMK_1 and UPMK_2) were isolated. Both phages varied in their ability to produce halos around their plaques, host infectivity, one-step growth curves, and electron microscopy features. Furthermore, both phages demonstrated antagonistic infectivity on planktonic cultures. This was validated in an in vitro static biofilm assay (in microtiter-plates), followed by the visualization of the biofilm architecture in situ via confocal laser scanning microscopy before and after phage infection, and further supported by phages genome analysis. The UPMK_1 genome comprised 152,788 bp coding for 155 putative open reading frames (ORFs), and its genome characteristics were between the Myoviridae and Siphoviridae family, though the morphological features confined it more to the Siphoviridae family. The UPMK_2 has 40,955 bp with 62 putative ORFs; morphologically, it presented the features of the Podoviridae though its genome did not show similarity with any of the S. aureus in the Podoviridae family. Both phages possess lytic enzymes that were associated with a high ability to degrade biofilms as shown in the microtiter plate and CLSM analyses.

Conclusions

The present work addressed the possibility of using phages as potential biocontrol agents for biofilm-producing MRSA.

Electronic supplementary material

The online version of this article (10.1186/s12866-019-1484-9) contains supplementary material, which is available to authorized users.

Global In-Silico Scenario of tRNA Genes and Their Organization in Virus Genomes

Viruses are known to be highly dependent on the host translation machinery for their protein synthesis. However, tRNA genes are occasionally identified in such organisms, and in addition, few of them harbor tRNA gene clusters comprising dozens of genes. Recently, tRNA gene clusters have been shown to occur among the three domains of life. In such a scenario, the viruses could play a role in the dispersion of such structures among these organisms. Thus, in order to reveal the prevalence of tRNA genes as well as tRNA gene clusters in viruses, we performed an unbiased large-scale genome survey. Interestingly, tRNA genes were predicted in ssDNA (single-stranded DNA) and ssRNA (single-stranded RNA) viruses as well in many other dsDNA viruses of families from Caudovirales order. In the latter group, tRNA gene clusters composed of 15 to 37 tRNA genes were characterized, mainly in bacteriophages, enlarging the occurrence of such structures within viruses. These bacteriophages were from hosts that encompass five phyla and 34 genera. This in-silico study presents the current global scenario of tRNA genes and their organization in virus genomes, contributing and opening questions to be explored in further studies concerning the role of the translation apparatus in these organisms.

Characterization and comparative genomic analysis of virulent and temperate Bacillus megaterium bacteriophages

Next-Generation Sequencing (NGS) technologies provide unique possibilities for the comprehensive assessment of the environmental diversity of bacteriophages. Several Bacillus bacteriophages have been isolated, but very few Bacillus megaterium bacteriophages have been characterized. In this study, we describe the biological characteristics, whole genome sequences, and annotations for two new isolates of the B. megaterium bacteriophages (BM5 and BM10), which were isolated from Egyptian soil samples. Growth analyses indicated that the phages BM5 and BM10 have a shorter latent period (25 and 30 min, respectively) and a smaller burst size (103 and 117 PFU, respectively), in comparison to what is typical for Bacillus phages. The genome sizes of the phages BM5 and BM10 were 165,031 bp and 165,213 bp, respectively, with modular organization. Bioinformatic analyses of these genomes enabled the assignment of putative functions to 97 and 65 putative ORFs, respectively. Comparative analysis of the BM5 and BM10 genome structures, in conjunction with other B. megaterium bacteriophages, revealed relatively high levels of sequence and organizational identity. Both genomic comparisons and phylogenetic analyses support the conclusion that the sequenced phages (BM5 and BM10) belong to different sub-clusters (L5 and L7, respectively), within the L-cluster, and display different lifestyles (lysogenic and lytic, respectively). Moreover, sequenced phages encode proteins associated with Bacillus pathogenesis. In addition, BM5 does not contain any tRNA sequences, whereas BM10 genome codes for 17 tRNAs.

Development and Validation of a Microtiter Plate-Based Assay for Determination of Bacteriophage Host Range and Virulence

Bacteriophages, which are the natural predators of bacteria, have re-emerged as an attractive alternative to combat antibiotic resistant bacteria. Phages are highly specific at the species and strain level and measurement of the phage host range plays an important role in utilizing the phage as antimicrobials. The most common method for phage host range determination has been to spot phage lysates on soft agar overlays and observe plaque formation. In this study, a liquid culture-based assay was developed in a 96-well microtiter plate format to measure the phage host range and virulence for a collection of 15 Salmonella phages against a panel of 20 Salmonella strains representing 11 serovars. This method was compared to a traditional spot method. The majority of the host range results from two methods were in agreement including in cases where a bacterial strain was insensitive to the phage. Each method produced a false-negative result in 19/300 (6%) of the measured phage-host combinations when compared to the other method. The spot method tended to indicate greater phage sensitivity than the microtiter assay even though direct comparisons of the response magnitude between the two methods is difficult since they operate on different mechanisms. The microtiter plate assay was able to provide data on both the phage host range and virulence in greater resolution in a high-throughput format.

Isolation and Characterization of Bacillus cereus Bacteriophages from Foods and Soil

The aim of this study was to isolate and characterize Bacillus cereus bacteriophages of various origins. Twenty-seven bacteriophages against B. cereus were isolated from various Korean traditional fermented foods and soils. Plaque size, transmission electron microscopy, virulence profile, and in vitro lytic activity of bacteriophage isolates were examined. Transmission electron microscopy confirmed B. cereus bacteriophages belonging to the family Siphoviridae. Among B. cereus bacteriophages with broad host range, 18 isolates (66.7%) did not harbor any B. cereus virulence factors. Among them, bacteriophage strain CAU150036, CAU150038, CAU150058, CAU150064, CAU150065, and CAU150066 effectively inhibited B. cereus in vitro within 1 h. Therefore, they are considered potential candidates for controlling the contamination of B. cereus in food or other applications.

Complete Genome Sequences of Bacillus Phages DirtyBetty and Kida

We report here the genome sequences of two bacteriophages of the Bacillus cereus group, DirtyBetty and Kida. These bacteriophages are double-stranded DNA-containing Myoviridae isolated from soil samples using Bacillus thuringiensis subsp. kurstaki as their host bacteria.

Safety assessment of Staphylococcus phages of the family Myoviridae based on complete genome sequences

Staphylococcus phages of the Myoviridae family have a wide host range and potential applications in phage therapy. In this report, safety assessments of these phages were conducted based on their complete genome sequences. The complete genomes of Staphylococcus phages of the Myoviridae family were analyzed, and the Open Reading Frame (ORFs) were compared with a pool of virulence and antibiotic resistance genes using the BLAST algorithm. In addition, the lifestyle of the phages (virulent or temperate) was also confirmed using PHACTS. The results showed that all phages were lytic and did not contain resistance or virulence genes based on bioinformatic analyses, excluding the possibility that they could be vectors for the dissemination of these undesirable genes. These findings suggest that the phages are safe at the genome level. The SceD-like transglycosylase, which is a biomarker for vancomycin-intermediate strains, was widely distributed in the phage genomes. Approximately 70% of the ORFs encoded in the phage genomes have unknown functions; therefore, their roles in the antibiotic resistance and virulence of Staphylococcus aureus are still unknown and require consideration before use in phage therapy.

Complete Genome Sequence of Bacteriophage Deep-Blue Infecting Emetic Bacillus cereus

The Bacillus cereus emetic pathotype is responsible for important food-borne intoxications. Here, we describe the complete genome sequence of bacteriophage Deep-Blue, which is able to infect emetic strains of B. cereus Deep-Blue is a 159-kb myophage of the Bastille-like group within the Spounavirinae.

Complete genome sequence of the cold-active bacteriophage VMY22 from Bacillus cereus

The cold-active bacteriophage VMY22, belonging to the Podoviridae family, was isolated from Mingyong Glacier in China. Sequence analysis revealed that the genome is 18,609 bp long, with an overall G + C content of 36.4 mol%, and 25 open reading frames (ORFs). The sequence contains 46 potential promoters, 6 transcription terminators, and no tRNAs. Most of the ORFs show a high degree of similarity to B103 (NC_004165). Two noteworthy findings were made. First, one of the predicted proteins, ORF 19, shows high sequence similarity to the bacteriocin biosynthesis protein from Bacillus cereus. From this information, we propose that the VMY22 phage is at an intermediate phase in its coevolution with its bacterial host. Second, seven of the hypothetical proteins appear to be unique to this cold-active B. cereus phage (i.e., not found in temperate-active B. cereus phages). These observations add to our current knowledge about the coevolution of bacteriophages and their hosts. The identification of a novel group of gene and protein structures and functions will lead to a better understanding of cold-adaptation mechanisms in bacteria and their bacteriophages.

Putative type 1 thymidylate synthase and dihydrofolate reductase as signature genes of a novel Bastille-like group of phages in the subfamily Spounavirinae

Background

Spounavirinae viruses have received an increasing interest as tools for the control of harmful bacteria due to their relatively broad host range and strictly virulent phenotype.

Results

In this study, we collected and analyzed the complete genome sequences of 61 published phages, either ICTV-classified or candidate members of the Spounavirinae subfamily of the Myoviridae. A set of comparative analyses identified a distinct, recently proposed Bastille-like phage group within the Spounavirinae. More importantly, type 1 thymidylate synthase (TS1) and dihydrofolate reductase (DHFR) genes were shown to be unique for the members of the proposed Bastille-like phage group, and are suitable as molecular markers. We also show that the members of this group encode beta-lactamase and/or sporulation-related SpoIIIE homologs, possibly questioning their suitability as biocontrol agents.

Conclusions

We confirm the creation of a new genus—the “Bastille-like group”—in Spounavirinae, and propose that the presence of TS1- and DHFR-encoding genes could serve as signatures for the new Bastille-like group. In addition, the presence of metallo-beta-lactamase and/or SpoIIIE homologs in all members of Bastille-like group phages makes questionable their suitability for use in biocontrol.

Electronic supplementary material

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

Complete genome sequence analysis and identification of putative metallo-beta-lactamase and SpoIIIE homologs in Bacillus cereus group phage BCP8-2, a new member of the proposed Bastille-like group

Bacillus cereus group-specific bacteriophage BCP8-2 exhibits a broad lysis spectrum among food and human isolates (330/364) of B. cereus while not infecting B. subtilis (50) or B. licheniformis (12) strains. Its genome is 159,071 bp long with 220 open reading frames, including genes for putative methyltransferases, metallo-beta-lactamase, and a sporulation-related SpoIIIE homolog, as wells as 18 tRNAs. Comparative genome analysis showed that BCP8-2 is related to the recently proposed Bastille-like phages, but not with either SPO1-like or Twort-like phages of the subfamily Spounavirinae.

Bacteriophages and bacteriophage-derived endolysins as potential therapeutics to combat Gram-positive spore forming bacteria

Since their discovery in 1915, bacteriophages have been routinely used within Eastern Europe to treat a variety of bacterial infections. Although initially ignored by the West due to the success of antibiotics, increasing levels and diversity of antibiotic resistance is driving a renaissance for bacteriophage-derived therapy, which is in part due to the highly specific nature of bacteriophages as well as their relative abundance. This review focuses on the bacteriophages and derived lysins of relevant Gram-positive spore formers within the Bacillus cereus group and Clostridium genus that could have applications within the medical, food and environmental sectors.

Complete genome sequence and phylogenetic position of the Bacillus cereus group phage JBP901

Bacteriophage JBP901, isolated from fermented food, is specific for Bacillus cereus group species and exhibits a broad host spectrum among a large number of B. cereus isolates. Genome sequence analysis revealed a linear 159,492-bp genome with overall G+C content of 39.7 mol%, and 201 ORFs. The presence of a putative methylase, the large number of tRNAs, and the large number of nucleotide-metabolism- and replication-related genes in JBP901 reflects its broad lytic capacity. Most of the ORFs showed a high degree of similarity to Bcp1, Bc431v3 and BCP78, and various comparative genomics analyses also consistently clustered JBP901 with orphan (unclassified) Bacillus phages in the subfamily Spounavirinae of the family Myoviridae, supporting the presence of a distinguishable group in the subfamily.

Isolation and characterization of glacier VMY22, a novel lytic cold-active bacteriophage of Bacillus cereus

As a unique ecological system with low temperature and low nutrient levels, glaciers are considered a "living fossil" for the research of evolution. In this work, a lytic cold-active bacteriophage designated VMY22 against Bacillus cereus MYB41-22 was isolated from Mingyong Glacier in China, and its characteristics were studied. Electron microscopy revealed that VMY22 has an icosahedral head (59.2 nm in length, 31.9 nm in width) and a tail (43.2 nm in length). Bacteriophage VMY22 was classified as a Podoviridae with an approximate genome size of 18 to 20 kb. A one-step growth curve revealed that the latent and the burst periods were 70 and 70 min, respectively, with an average burst size of 78 bacteriophage particles per infected cell. The pH and thermal stability of bacteriophage VMY22 were also investigated. The maximum stability of the bacteriophage was observed to be at pH 8.0 and it was comparatively stable at pH 5.0-9.0. As VMY22 is a cold-active bacteriophage with low production temperature, its characterization and the relationship between MYB41-22 and Bacillus cereus bacteriophage deserve further study.

Bacteriophage PBC1 and its endolysin as an antimicrobial agent against Bacillus cereus

Bacillus cereus is an opportunistic human pathogen responsible for food poisoning and other, nongastrointestinal infections. Due to the emergence of multidrug-resistant B. cereus strains, the demand for alternative therapeutic options is increasing. To address these problems, we isolated and characterized a Siphoviridae virulent phage, PBC1, and its lytic enzymes. PBC1 showed a very narrow host range, infecting only 1 of 22 B. cereus strains. Phylogenetic analysis based on the major capsid protein revealed that PBC1 is more closely related to the Bacillus clarkii phage BCJA1c and phages of lactic acid bacteria than to the phages infecting B. cereus. Whole-genome comparison showed that the late-gene region, including the terminase gene, structural genes, and holin gene of PBC1, is similar to that from B. cereus temperate phage 250, whereas their endolysins are different. Compared to the extreme host specificity of PBC1, its endolysin, LysPBC1, showed a much broader lytic spectrum, albeit limited to the genus Bacillus. The catalytic domain of LysPBC1 when expressed alone also showed Bacillus-specific lytic activity, which was lower against the B. cereus group but higher against the Bacillus subtilis group than the full-length protein. Taken together, these results suggest that the virulent phage PBC1 is a useful component of a phage cocktail to control B. cereus, even with its exceptionally narrow host range, as it can kill a strain of B. cereus that is not killed by other phages, and that LysPBC1 is an alternative biocontrol agent against B. cereus.

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.

Genome analysis of Cronobacter phage vB_CsaP_Ss1 reveals an endolysin with potential for biocontrol of Gram-negative bacterial pathogens

Bacteriophages and their derivatives are continuously gaining impetus as viable alternative therapeutic agents to control harmful multidrug-resistant bacterial pathogens, particularly in the food industry. The reduced efficacy of conventional antibiotics has resulted in a quest to find novel alternatives in the war against infectious disease. This study describes the full-genome sequence of Cronobacter phage vB_CsaP_Ss1, with subsequent cloning and expression of its endolysin, capable of hydrolysing Gram-negative peptidoglycan. Cronobacter phage vB_CsaP_Ss1 is composed of 42 205 bp of dsDNA with a G+C content of 46.1 mol%. A total of 57 ORFs were identified of which 18 could be assigned a putative function based on similarity to characterized proteins. The genome of Cronobacter phage vB_CsaP_Ss1 showed little similarity to any other bacteriophage genomes available in the database and thus was considered unique. In addition, functional analysis of the predicted endolysin (LysSs1) was also investigated. Zymographic experiments demonstrated the hydrolytic activity of LysSs1 against Gram-negative peptidoglycan, and this endolysin thus represents a novel candidate with potential for use against Gram-negative pathogens.

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.

Possible use of bacteriophages active against Bacillus anthracis and other B. cereus group members in the face of a bioterrorism threat

Anthrax is an infectious fatal disease with epidemic potential. Nowadays, bioterrorism using Bacillus anthracis is a real possibility, and thus society needs an effective weapon to neutralize this threat. The pathogen may be easily transmitted to human populations. It is easy to store, transport, and disseminate and may survive for many decades. Recent data strongly support the effectiveness of bacteriophage in treating bacterial diseases. Moreover, it is clear that bacteriophages should be considered a potential incapacitative agent against bioterrorism using bacteria belonging to B. cereus group, especially B. anthracis. Therefore, we have reviewed the possibility of using bacteriophages active against Bacillus anthracis and other species of the B. cereus group in the face of a bioterrorism threat.

The genomes, proteomes, and structures of three novel phages that infect the Bacillus cereus group and carry putative virulence factors

This article reports the results of studying three novel bacteriophages, JL, Shanette, and Basilisk, which infect the pathogen Bacillus cereus and carry genes that may contribute to its pathogenesis. We analyzed host range and superinfection ability, mapped their genomes, and characterized phage structure by mass spectrometry and transmission electron microscopy (TEM). The JL and Shanette genomes were 96% similar and contained 217 open reading frames (ORFs) and 220 ORFs, respectively, while Basilisk has an unrelated genome containing 138 ORFs. Mass spectrometry revealed 23 phage particle proteins for JL and 15 for Basilisk, while only 11 and 4, respectively, were predicted to be present by sequence analysis. Structural protein homology to well-characterized phages suggested that JL and Shanette were members of the family Myoviridae, which was confirmed by TEM. The third phage, Basilisk, was similar only to uncharacterized phages and is an unrelated siphovirus. Cryogenic electron microscopy of this novel phage revealed a T=9 icosahedral capsid structure with the major capsid protein (MCP) likely having the same fold as bacteriophage HK97 MCP despite the lack of sequence similarity. Several putative virulence factors were encoded by these phage genomes, including TerC and TerD involved in tellurium resistance. Host range analysis of all three phages supports genetic transfer of such factors within the B. cereus group, including B. cereus, B. anthracis, and B. thuringiensis. This study provides a basis for understanding these three phages and other related phages as well as their contributions to the pathogenicity of B. cereus group bacteria. Importance: The Bacillus cereus group of bacteria contains several human and plant pathogens, including B. cereus, B. anthracis, and B. thuringiensis. Phages are intimately linked to the evolution of their bacterial hosts and often provide virulence factors, making the study of B. cereus phages important to understanding the evolution of pathogenic strains. Herein we provide the results of detailed study of three novel B. cereus phages, two highly related myoviruses (JL and Shanette) and an unrelated siphovirus (Basilisk). The detailed characterization of host range and superinfection, together with results of genomic, proteomic, and structural analyses, reveal several putative virulence factors as well as the ability of these phages to infect different pathogenic species.

Phages preying on Bacillus anthracis, Bacillus cereus, and Bacillus thuringiensis: past, present and future

Many bacteriophages (phages) have been widely studied due to their major role in virulence evolution of bacterial pathogens. However, less attention has been paid to phages preying on bacteria from the Bacillus cereus group and their contribution to the bacterial genetic pool has been disregarded. Therefore, this review brings together the main information for the B. cereus group phages, from their discovery to their modern biotechnological applications. A special focus is given to phages infecting Bacillus anthracis, B. cereus and Bacillus thuringiensis. These phages belong to the Myoviridae, Siphoviridae, Podoviridae and Tectiviridae families. For the sake of clarity, several phage categories have been made according to significant characteristics such as lifestyles and lysogenic states. The main categories comprise the transducing phages, phages with a chromosomal or plasmidial prophage state, γ-like phages and jumbo-phages. The current genomic characterization of some of these phages is also addressed throughout this work and some promising applications are discussed here.

Clostridium difficile phages: still difficult?

Phages that infect Clostridium difficile were first isolated for typing purposes in the 1980s, but their use was short lived. However, the rise of C. difficile epidemics over the last decade has triggered a resurgence of interest in using phages to combat this pathogen. Phage therapy is an attractive treatment option for C. difficile infection, however, developing suitable phages is challenging. In this review we summarize the difficulties faced by researchers in this field, and we discuss the solutions and strategies used for the development of C. difficile phages for use as novel therapeutics. Epidemiological data has highlighted the diversity and distribution of C. difficile, and shown that novel strains continue to emerge in clinical settings. In parallel with epidemiological studies, advances in molecular biology have bolstered our understanding of C. difficile biology, and our knowledge of phage–host interactions in other bacterial species. These three fields of biology have therefore paved the way for future work on C. difficile phages to progress and develop. Benefits of using C. difficile phages as therapeutic agents include the fact that they have highly specific interactions with their bacterial hosts. Studies also show that they can reduce bacterial numbers in both in vitro and in vivo systems. Genetic analysis has revealed the genomic diversity among these phages and provided an insight into their taxonomy and evolution. No strictly virulent C. difficile phages have been reported and this contributes to the difficulties with their therapeutic exploitation. Although treatment approaches using the phage-encoded endolysin protein have been explored, the benefits of using “whole-phages” are such that they remain a major research focus. Whilst we don’t envisage working with C. difficile phages will be problem-free, sufficient study should inform future strategies to facilitate their development to combat this problematic pathogen.

Characterization and genome analysis of the Bacillus cereus-infecting bacteriophages BPS10C and BPS13

Due to the emergence of antibiotic-resistant strains, bacteriophages are considered to be an alternative approach for the control of pathogens. In this study, the bacteriophages BPS10C and BPS13 were isolated and characterized to investigate their ability to control food-borne pathogenic Bacillus cereus. Phage BPS13 exhibited slightly higher host lysis activity compared with phage BPS10C. In addition, phage BPS13 exhibited greater stability under various pH and temperature conditions. To extend our knowledge of the lysis of B. cereus by these phages, their genomes were completely sequenced and analyzed, revealing that these phage genomes encode endolysin and two tail lysins, which are likely involved in host lysis and invasion mechanisms, respectively. These lysis-related proteins may increase the bactericidal activities of these phages, suggesting that they may be good candidates for the potential control of B. cereus.

Genomic analysis of a phage and prophage from a Bacillus thuringiensis strain

Bacteriophages have been found to be the most abundant and also potentially most diverse biological entities on Earth. In the present study, Bacillus phages were isolated rapidly and shown to have a high degree of diversity. The genomes of a newly isolated phage, phiCM3, and a prophage, proCM3, from the Bacillus thuringiensis strain YM-03 were sequenced and characterized. Comparative genome analysis showed that the phiCM3 genome is highly similar to the genomes of eight other Bacillus phages and seven of these phages were classified as the Wβ group of phages. Analysis of the differential evolution of the genes in the Wβ-group phages indicated that the genes encoding the antirepressor and tail fibre protein were more highly conserved than those encoding the major capsid protein, DNA replication protein, and RNA polymerase σ factor, which might have diverged to acquire mechanisms suitable for survival in different microbial hosts. Genome analysis of proCM3 revealed that proCM3 might be a defective phage because of mutations in the minor structural protein, and it was not inducible by mitomycin C treatment. The proCM3 genome was similar to those of two lytic Bacillus phages in sequence, but had a different genomic structure, composed of three regions in a different order. These data suggest that the three phages might have had a common ancestor and that genome rearrangement might have occurred during evolution. The findings of this study enrich our current knowledge of Bacillus phage diversity and evolution, especially for the Wβ-group and TP21-L-like phages, and may help the development of practical applications of Bacillus phages.

Characterization and comparative genomic analysis of bacteriophages infecting members of the Bacillus cereus group

The Bacillus cereus group phages infecting B. cereus, B. anthracis, and B. thuringiensis (Bt) have been studied at the molecular level and, recently, at the genomic level to control the pathogens B. cereus and B. anthracis and to prevent phage contamination of the natural insect pesticide Bt. A comparative phylogenetic analysis has revealed three different major phage groups with different morphologies (Myoviridae for group I, Siphoviridae for group II, and Tectiviridae for group III), genome size (group I > group II > group III), and lifestyle (virulent for group I and temperate for group II and III). A subsequent phage genome comparison using a dot plot analysis showed that phages in each group are highly homologous, substantiating the grouping of B. cereus phages. Endolysin is a host lysis protein that contains two conserved domains: a cell-wall-binding domain (CBD) and an enzymatic activity domain (EAD). In B. cereus sensu lato phage group I, four different endolysin groups have been detected, according to combinations of two types of CBD and four types of EAD. Group I phages have two copies of tail lysins and one copy of endolysin, but the functions of the tail lysins are still unknown. In the B. cereus sensu lato phage group II, the B. anthracis phages have been studied and applied for typing and rapid detection of pathogenic host strains. In the B. cereus sensu lato phage group III, the B. thuringiensis phages Bam35 and GIL01 have been studied to understand phage entry and lytic switch regulation mechanisms. In this review, we suggest that further study of the B. cereus group phages would be useful for various phage applications, such as biocontrol, typing, and rapid detection of the pathogens B. cereus and B. anthracis and for the prevention of phage contamination of the natural insect pesticide Bt.

Characterization and complete genome sequence of a virulent bacteriophage B4 infecting food-borne pathogenic Bacillus cereus

Bacillus cereus causes food poisoning, resulting in vomiting and diarrhea, due to production of enterotoxins. As a means of controlling this food-borne pathogen, the virulent bacteriophage B4 was isolated and characterized. Bacterial challenge assays showed that phage B4 effectively inhibited growth of members of the B. cereus group as well as B. subtilis, and growth inhibition persisted for over 20 h. One-step growth analysis also revealed the host lysis activity of phage B4, with relatively short eclipse/latent times (10/15 min) and a large burst size (>200 PFU). The complete genome of phage B4, containing a 162-kb DNA with 277 ORFs, was analyzed. The endolysin encoded by the phage B4 genome accounts for the cell lysis activity of this phage. These results suggest that phage B4 has potential as a biological agent to control B. cereus propagation.