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.     

The isolation and characterization of two Stenotrophomonas maltophilia bacteriophages capable of cross-taxonomic order infectivity

Background

A rapid worldwide increase in the number of human infections caused by the extremely antibiotic resistant bacterium Stenotrophomonas maltophilia is prompting alarm. One potential treatment solution to the current antibiotic resistance dilemma is “phage therapy”, the clinical application of bacteriophages to selectively kill bacteria.

Results

Towards that end, phages DLP1 and DLP2 (vB_SmaS-DLP_1 and vB_SmaS-DLP_2, respectively) were isolated against S. maltophilia strain D1585. Host range analysis for each phage was conducted using 27 clinical S. maltophilia isolates and 11 Pseudomonas aeruginosa strains. Both phages exhibit unusually broad host ranges capable of infecting bacteria across taxonomic orders. Transmission electron microscopy of the phage DLP1 and DLP2 morphology reveals that they belong to the Siphoviridae family of bacteriophages. Restriction fragment length polymorphism analysis and complete genome sequencing and analysis indicates that phages DLP1 and DLP2 are closely related but different phages, sharing 96.7 % identity over 97.2 % of their genomes. These two phages are also related to P. aeruginosa phages vB_Pae-Kakheti_25 (PA25), PA73, and vB_PaeS_SCH_Ab26 (Ab26) and more distantly related to Burkholderia cepacia complex phage KL1, which together make up a taxonomic sub-family. Phages DLP1 and DLP2 exhibited significant differences in host ranges and growth kinetics.

Conclusions

The isolation and characterization of phages able to infect two completely different species of bacteria is an exciting discovery, as phages typically can only infect related bacterial species, and rarely infect bacteria across taxonomic families, let alone across taxonomic orders.

Electronic supplementary material

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

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.     

Methyltransferases acquired by lactococcal 936-type phage provide protection against restriction endonuclease activity

Background

So-called 936-type phages are among the most frequently isolated phages in dairy facilities utilising Lactococcus lactis starter cultures. Despite extensive efforts to control phage proliferation and decades of research, these phages continue to negatively impact cheese production in terms of the final product quality and consequently, monetary return.

Results

Whole genome sequencing and in silico analysis of three 936-type phage genomes identified several putative (orphan) methyltransferase (MTase)-encoding genes located within the packaging and replication regions of the genome. Utilising SMRT sequencing, methylome analysis was performed on all three phages, allowing the identification of adenine modifications consistent with N-6 methyladenine sequence methylation, which in some cases could be attributed to these phage-encoded MTases. Heterologous gene expression revealed that M.Phi145I/M.Phi93I and M.Phi93DAM, encoded by genes located within the packaging module, provide protection against the restriction enzymes HphI and DpnII, respectively, representing the first functional MTases identified in members of 936-type phages.

Conclusions

SMRT sequencing technology enabled the identification of the target motifs of MTases encoded by the genomes of three lytic 936-type phages and these MTases represent the first functional MTases identified in this species of phage. The presence of these MTase-encoding genes on 936-type phage genomes is assumed to represent an adaptive response to circumvent host encoded restriction-modification systems thereby increasing the fitness of the phages in a dynamic dairy environment.

Electronic supplementary material

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

Campylobacter group II phage CP21 is the prototype of a new subgroup revealing a distinct modular genome organization and host specificity

Background

The application of phages is a promising tool to reduce the number of Campylobacter along the food chain. Besides the efficacy against a broad range of strains, phages have to be safe in terms of their genomes. Thus far, no genes with pathogenic potential (e.g., genes encoding virulence factors) have been detected in Campylobacter phages. However, preliminary studies suggested that the genomes of group II phages may be diverse and prone to genomic rearrangements.

Results

We determined and analysed the genomic sequence (182,761 bp) of group II phage CP21 that is closely related to the already characterized group II phages CP220 and CPt10. The genomes of these phages are comprised of four modules separated by very similar repeat regions, some of which harbouring open reading frames (ORFs). Though, the arrangement of the modules and the location of some ORFs on the genomes are different in CP21 and in CP220/CPt10. In this work, a PCR system was established to study the modular genome organization of other group II phages demonstrating that they belong to different subgroups of the CP220-like virus genus, the prototypes of which are CP21 and CP220. The subgroups revealed different restriction patterns and, interestingly enough, also distinct host specificities, tail fiber proteins and tRNA genes. We additionally analysed the genome of group II phage vB_CcoM-IBB_35 (IBB_35) for which to date only five individual contigs could be determined. We show that the contigs represent modules linked by long repeat regions enclosing some yet not identified ORFs (e.g., for a head completion protein). The data suggest that IBB_35 is a member of the CP220 subgroup.

Conclusion

Campylobacter group II phages are diverse regarding their genome organization. Since all hitherto characterized group II phages contain numerous genes for transposases and homing endonucleases as well as similar repeat regions, it cannot be excluded that these phages are genetically unstable. To answer this question, further experiments and sequencing of more group II phages should be performed.

Electronic supplementary material

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

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.     

The Human-Bacterial Pathogen Protein Interaction Networks of Bacillus anthracis,Francisella tularensis,and Yersinia pestis

Background

Bacillus anthracis, Francisella tularensis, and Yersinia pestis are bacterial pathogens that can cause anthrax, lethal acute pneumonic disease, and bubonic plague, respectively, and are listed as NIAID Category A priority pathogens for possible use as biological weapons. However, the interactions between human proteins and proteins in these bacteria remain poorly characterized leading to an incomplete understanding of their pathogenesis and mechanisms of immune evasion.

Methodology

In this study, we used a high-throughput yeast two-hybrid assay to identify physical interactions between human proteins and proteins from each of these three pathogens. From more than 250,000 screens performed, we identified 3,073 human-B. anthracis, 1,383 human-F. tularensis, and 4,059 human-Y. pestis protein-protein interactions including interactions involving 304 B. anthracis, 52 F. tularensis, and 330 Y. pestis proteins that are uncharacterized. Computational analysis revealed that pathogen proteins preferentially interact with human proteins that are hubs and bottlenecks in the human PPI network. In addition, we computed modules of human-pathogen PPIs that are conserved amongst the three networks. Functionally, such conserved modules reveal commonalities between how the different pathogens interact with crucial host pathways involved in inflammation and immunity.

Significance

These data constitute the first extensive protein interaction networks constructed for bacterial pathogens and their human hosts. This study provides novel insights into host-pathogen interactions.     

Intraspecies comparison of Streptomyces pratensis genomes reveals high levels of recombination and gene conservation between strains of disparate geographic origin

Background

Streptomyces are widespread bacteria that contribute to the terrestrial carbon cycle and produce the majority of clinically useful antibiotics. While interspecific genomic diversity has been investigated among Streptomyces, information is lacking on intraspecific genomic diversity. Streptomyces pratensis has high rates of homologous recombination but the impact of such gene exchange on genome evolution and the evolution of natural product gene clusters remains uncharacterized.

Results

We report draft genome sequences of four S. pratensis strains and compare to the complete genome of Streptomyces flavogriseus IAF-45-CD (=ATCC 33331), a strain recently reclassified to S. pratensis. Despite disparate geographic origins, the genomes are highly similar with 85.9% of genes present in the core genome and conservation of all natural product gene clusters. Natural products include a novel combination of carbapenem and beta-lactamase inhibitor gene clusters. While high intraspecies recombination rates abolish the phylogenetic signal across the genome, intraspecies recombination is suppressed in two genomic regions. The first region is centered on an insertion/deletion polymorphism and the second on a hybrid NRPS-PKS gene. Finally, two gene families accounted for over 25% of the divergent genes in the core genome. The first includes homologs of bldB (required for spore development and antibiotic production) while the second includes homologs of an uncharacterized protein with a helix-turn-helix motif (hpb). Genes from these families co-occur with fifteen pairs spread across the genome. These genes have evidence for co-evolution of co-localized pairs, supporting previous assertions that these genes may function akin to a toxin-antitoxin system.

Conclusions

S. pratensis genomes are highly similar with exceptional levels of recombination which erase phylogenetic signal among strains of the species. This species has a large core genome and variable terminal regions that are smaller than those found in interspecies comparisons. There is no geographic differentiation between these strains, but there is evidence for local linkage disequilibrium affecting two genomic regions. We have also shown further observational evidence that the DUF397-HTH (bldB and hpb) are a novel toxin-antitoxin pair.     

Bacillus cereus Spores Release Alanine that Synergizes with Inosine to Promote Germination

Background

The first step of the bacterial lifecycle is the germination of bacterial spores into their vegetative form, which requires the presence of specific nutrients. In contrast to closely related Bacillus anthracis spores, Bacillus cereus spores germinate in the presence of a single germinant, inosine, yet with a significant lag period.

Methods and Findings

We found that the initial lag period of inosine-treated germination of B. cereus spores disappeared in the presence of supernatants derived from already germinated spores. The lag period also dissipated when inosine was supplemented with the co-germinator alanine. In fact, HPLC-based analysis revealed the presence of amino acids in the supernatant of germinated B. cereus spores. The released amino acids included alanine in concentrations sufficient to promote rapid germination of inosine-treated spores. The alanine racemase inhibitor D-cycloserine enhanced germination of B. cereus spores, presumably by increasing the L-alanine concentration in the supernatant. Moreover, we found that B. cereus spores lacking the germination receptors gerI and gerQ did not germinate and release amino acids in the presence of inosine. These mutant spores, however, germinated efficiently when inosine was supplemented with alanine. Finally, removal of released amino acids in a washout experiment abrogated inosine-mediated germination of B. cereus spores.

Conclusions

We found that the single germinant inosine is able to trigger a two-tier mechanism for inosine-mediated germination of B. cereus spores: Inosine mediates the release of alanine, an essential step to complete the germination process. Therefore, B. cereus spores appear to have developed a unique quorum-sensing feedback mechanism to monitor spore density and to coordinate germination.     

Receptor Diversity and Host Interaction of Bacteriophages Infecting Salmonella enterica Serovar Typhimurium

Background

Salmonella enterica subspecies enterica serovar Typhimurium is a Gram-negative pathogen causing salmonellosis. Salmonella Typhimurium-targeting bacteriophages have been proposed as an alternative biocontrol agent to antibiotics. To further understand infection and interaction mechanisms between the host strains and the bacteriophages, the receptor diversity of these phages needs to be elucidated.

Methodology/Principal Findings

Twenty-five Salmonella phages were isolated and their receptors were identified by screening a Tn5 random mutant library of S. Typhimurium SL1344. Among them, three types of receptors were identified flagella (11 phages), vitamin B₁₂ uptake outer membrane protein, BtuB (7 phages) and lipopolysaccharide-related O-antigen (7 phages). TEM observation revealed that the phages using flagella (group F) or BtuB (group B) as a receptor belong to Siphoviridae family, and the phages using O-antigen of LPS as a receptor (group L) belong to Podoviridae family. Interestingly, while some of group F phages (F-I) target FliC host receptor, others (F-II) target both FliC and FljB receptors, suggesting that two subgroups are present in group F phages. Cross-resistance assay of group B and L revealed that group L phages could not infect group B phage-resistant strains and reversely group B phages could not infect group L SPN9TCW-resistant strain.

Conclusions/Significance

In this report, three receptor groups of 25 newly isolated S. Typhimurium-targeting phages were determined. Among them, two subgroups of group F phages interact with their host receptors in different manner. In addition, the host receptors of group B or group L SPN9TCW phages hinder other group phage infection, probably due to interaction between receptors of their groups. This study provides novel insights into phage-host receptor interaction for Salmonella phages and will inform development of optimal phage therapy for protection against Salmonella.     

Peptides panned from a phage-displayed random peptide library are useful for the detection of Bacillus anthracis surrogates B. cereus 4342 and B. anthracis Sterne

Recent use of Bacillus anthracis as a bioweapon has highlighted the need for a sensitive monitoring system. Current bacterial detection tests use antibodies as bio-molecular recognition elements which have limitations with regard to time, specificity and sensitivity, creating the need for new and improved cost-effective high-affinity detection probes. In this study, we screened a commercially available bacteriophage-displayed random peptide library using Bacillus cereus 4342 cells as bait to identify peptides that could be used for detection of Bacillus. The method enabled us to identify two 12-amino acid consensus peptide sequences that specifically bind to B. cereus 4342 and B. anthracis Sterne, the nonpathogenic surrogates of B. anthracis strain. The two Bacillus-binding peptides (named BBP-1 and BBP-2) were synthesized with biotin tag to confirm their binding by four independent detection assays. Dot-blot analysis revealed that the peptides bind specifically to B. cereus 4342 and B. anthracis Sterne. Quantitative analysis of this interaction by ELISA and fluorometry demonstrated a detection sensitivity of 10² colony forming U/ml (CFU/ml) by both assays. When the peptides were used in combination with Qdots, the sensitivity was enhanced further by enabling detection of even a single bacterium by fluorescence microscopy. Immunoblot analysis and protein sequencing showed that BBP-1 and BBP-2 bound to the S-layer protein of B. anthracis Sterne. Overall, our findings validate the usefulness of synthetic versions of phage-derived peptides in combination with Qdot-liquid nanocrystals as high sensitivity bioprobes for various microbial detection platforms.     

Molecular Characterization of Korean Bacillus anthracis Isolates by Amplified Fragment Length Polymorphism Analysis and Multilocus Variable-Number Tandem Repeat Analysis

We analyzed the genetic relationships and molecular characteristics of 34 Bacillus anthracis isolates from soil and clinical samples in various regions of Korea and 17 related Bacillus species, using the amplified fragment length polymorphism (AFLP) and multilocus variable-number tandem repeat (MLVA) approaches. Triplicate AFLP profiles of these strains showed high reproducibility and identified 376 polymorphisms. AFLP phylogenetic analysis of B. anthracis isolates showed a high level of similarity, 0.93, and this monomorphic fragment profile proved to be useful to differentiate B. anthracis strains from other Bacillus species. The B. cereus group was separated from other Bacillus species at a level of similarity of 0.68. Among them, some B. cereus strains showed genetic interspersion with B. thuringiensis strains. The evolutionary pattern of nucleotide differences among B. anthracis strains with the eight MLVA markers showed nine MLVA types. Three MLVA types, M1 to M3, were pathogenic B. anthracis isolates and were assigned as new genotypes belonging to the A4 and B3 clusters, compared with 89 genotypes deduced from previous data. This indicates that differences in cluster prevalence and distribution may be influenced more by MLVA markers on two plasmids loci and human activity. Consequently, we suggest that the novel MLVA type may represent significant evidence for historic adaptation to environmental conditions of the Asian continent, particularly Korea. Therefore, MLVA techniques may be available for molecular monitoring on anthrax-release-related bioterrorism and further study is required for the continuous epidemiological study of variable anthrax collections.     

Bacillus anthracis Multiplication, Persistence, and Genetic Exchange in the Rhizosphere of Grass Plants

Bacillus anthracis, the causative agent of anthrax, is known for its rapid proliferation and dissemination in mammalian hosts. In contrast, little information exists regarding the lifestyle of this important pathogen outside of the host. Considering that Bacillus species, including close relatives of B. anthracis, are saprophytic soil organisms, we investigated the capacity of B. anthracis spores to germinate in the rhizosphere and to establish populations of vegetative cells that could support horizontal gene transfer in the soil. Using a simple grass plant-soil model system, we show that B. anthracis strains germinate on and around roots, growing in characteristic long filaments. From 2 to 4 days postinoculation, approximately one-half of the B. anthracis CFU recovered from soil containing grass seedlings arose from heat-sensitive organisms, while B. anthracis CFU retrieved from soil without plants consisted of primarily heat-resistant spores. Coinoculation of the plant-soil system with spores of a fertile B. anthracis strain carrying the tetracycline resistance plasmid pBC16 and a selectable B. anthracis recipient strain resulted in transfer of pBC16 from the donor to the recipient as early as 3 days postinoculation. Our findings demonstrate that B. anthracis can survive as a saprophyte outside of the host. The data suggest that horizontal gene transfer in the rhizosphere of grass plants may play a role in the evolution of the Bacillus cereus group species.     

Compositional studies and Biological activities of Perovskia abrotanoides Kar. oils

Background

Current study has been designed to evaluate the chemical composition of essential and fixed oils from stem and leaves of Perovskia abrotanoides and antioxidant and antimicrobial activities of these oils.

Results

GC-MS analysis of essential oil identified 19 compounds with (E)-9-dodecenal being the major component in stem and hexadecanoic acid in leaves. In contrast, GC-MS analysis of fixed oil showed 40 constituents with α-amyrin the major component in stem and α-copaene in leaves. The antioxidant activity showed the highest value of 76.7% in essential oil from leaves in comparison with fixed oil from stem (45.9%) through inhibition of peroxidation in linoleic acid system. The antimicrobial assay tested on different microorganisms (e.g. E. coli, S. aureus, B. cereus, Nitrospira, S. epidermis, A. niger, A. flavus and C. albicans) showed the higher inhibition zone at essential oil from leaves (15.2 mm on B. cereus) as compared to fixed oil from stem (8.34 mm on S. aureus) and leaves (11.2 mm on S. aureus).

Conclusions

The present study revealed the fact that essential oil analyzed from Perovskia abrotanoides stem and leaves could be a promising source of natural products with potential antioxidant and antimicrobial activities, as compared to fixed oil.     

Genomic and transcriptomic analysis of the endophytic fungus Pestalotiopsis fici reveals its lifestyle and high potential for synthesis of natural products

Background

In recent years, the genus Pestalotiopsis is receiving increasing attention, not only because of its economic impact as a plant pathogen but also as a commonly isolated endophyte which is an important source of bioactive natural products. Pestalotiopsis fici Steyaert W106-1/CGMCC3.15140 as an endophyte of tea produces numerous novel secondary metabolites, including chloropupukeananin, a derivative of chlorinated pupukeanane that is first discovered in fungi. Some of them might be important as the drug leads for future pharmaceutics.

Results

Here, we report the genome sequence of the endophytic fungus of tea Pestalotiopsis fici W106-1/CGMCC3.15140. The abundant carbohydrate-active enzymes especially significantly expanding pectinases allow the fungus to utilize the limited intercellular nutrients within the host plants, suggesting adaptation of the fungus to endophytic lifestyle. The P. fici genome encodes a rich set of secondary metabolite synthesis genes, including 27 polyketide synthases (PKSs), 12 non-ribosomal peptide synthases (NRPSs), five dimethylallyl tryptophan synthases, four putative PKS-like enzymes, 15 putative NRPS-like enzymes, 15 terpenoid synthases, seven terpenoid cyclases, seven fatty-acid synthases, and five hybrids of PKS-NRPS. The majority of these core enzymes distributed into 74 secondary metabolite clusters. The putative Diels-Alderase genes have undergone expansion.

Conclusion

The significant expansion of pectinase encoding genes provides essential insight in the life strategy of endophytes, and richness of gene clusters for secondary metabolites reveals high potential of natural products of endophytic fungi.

Electronic supplementary material

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

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.