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

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

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

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.

Characterization of G-type Clostridium perfringens bacteriophages and their disinfection effect on chicken meat

OBJECTIVE

Clostridium perfringens is one of most important bacterial pathogens in the poultry industry and mainly causes necrotizing enteritis (NE). This pathogen and its toxins can cause foodborne diseases in humans through the food chain. In China, with the rise of antibiotic resistance and the banning of antibiotic growth promoters (AGPs) in poultry farming, food contamination and NE are becoming more prevalent. Bacteriophages are a viable technique to control C. perfringens as an alternative to antibiotics. We isolated Clostridium phage from the environment, providing a new method for the prevention of NE and C. perfringens contamination in meat.

METHODS

In this study, we selected C. perfringens strains from various regions and animal sources in China for phage isolation. The biological characteristics of Clostridium phage were studied in terms of host range, MOI, one-step curve, temperature and pH stability. We sequenced and annotated the genome of the Clostridium phage and performed phylogenetic and pangenomic analyses. Finally, we studied its antibacterial activity against bacterial culture and its disinfection effect against C. perfringens in meat.

RESULTS

A Clostridium phage, named ZWPH-P21 (P21), was isolated from chicken farm sewage in Jiangsu, China. P21 has been shown to specifically lyse C. perfringens type G. Further analysis of basic biological characteristics showed that P21 was stable under the conditions of pH 4-11 and temperature 4-60 °C, and the optimal multiple severity of infection (MOI) was 0.1. In addition, P21 could form a "halo" on agar plates, suggesting that the phage may encode depolymerase. Genome sequence analysis showed that P21 was the most closely related to Clostridium phage CPAS-15 belonging to the Myoviridae family, with a recognition rate of 97.24% and a query coverage rate of 98%. No virulence factors or drug resistance genes were found in P21. P21 showed promising antibacterial activity in vitro and in chicken disinfection experiments. In conclusion, P21 has the potential to be used for preventing and controlling C. perfringens in chicken food production.

Genome, biology and stability of the Thurquoise phage – A new virus from the Bastillevirinae subfamily

Bacteriophages from the Bastillevirinae subfamily (Herelleviridae family) have proven to be effective against bacteria from the Bacillus genus including organisms from the B. cereus group, which cause food poisoning and persistent contamination of industrial installations. However, successful application of these phages in biocontrol depends on understanding of their biology and stability in different environments. In this study, we isolated a novel virus from garden soil in Wrocław (Poland) and named it ‘Thurquoise’. The genome of that phage was sequenced and assembled into a single continuous contig with 226 predicted protein-coding genes and 18 tRNAs. The cryo-electron microscopy revealed that Thurquoise has complex virion structure typical for the Bastillevirinae family. Confirmed hosts include selected bacteria from the Bacillus cereus group–specifically B. thuringiensis (isolation host) and B. mycoides, but susceptible strains display different efficiency of plating (EOP). The eclipse and latent periods of Thurquoise in the isolation host last ~ 50 min and ~ 70 min, respectively. The phage remains viable for more than 8 weeks in variants of the SM buffer with magnesium, calcium, caesium, manganese or potassium and can withstand numerous freeze–thaw cycles if protected by the addition of 15% glycerol or, to a lesser extent, 2% gelatine. Thus, with proper buffer formulation, this virus can be safely stored in common freezers and refrigerators for a considerable time. The Thurquoise phage is the exemplar of a new candidate species within the Caeruleovirus genus in the Bastillevirinae subfamily of the Herelleviridae family with a genome, morphology and biology typical for these taxa.

More's the Same-Multiple Hosts Do Not Select for Broader Host Range Phages

Bacteriophage host range is a result of the interactions between phages and their hosts. For phage therapy, phages with a broader host range are desired so that a phage can infect and kill the broadest range of pathogen strains or related species possible. A common, but not well-tested, belief is that using multiple hosts during the phage isolation will make the isolation of broader host range phage more likely. Using a Bacillus cereus group system, we compared the host ranges of phages isolated on one or four hosts and found that there was no difference in the breadth of host ranges of the isolated phages. Both narrow and broader host range phage were also equally likely to be isolated from either isolation procedure. While there are methods that reliably isolate broader host range phages, such as sequential host isolation, and there are other reasons to use multiple hosts during isolation, multiple hosts are not a consistent way to obtain broader host range phages.

Genomic diversity and comprehensive taxonomical classification of 61 Bacillus subtilis group member infecting bacteriophages, and the identification of ortholog taxonomic signature genes

BACKGROUND

Despite the applications of Bacillus subtilis group species in various sectors, limited information is available regarding their phages. Here, 61 B. subtilis group species-infecting phages (BSPs) were studied for their taxonomic classification considering the genome-size, genomic diversity, and the host, followed by the identification of orthologs taxonomic signature genes.

RESULTS

BSPs have widely ranging genome sizes that can be bunched into groups to demonstrate correlations to family and subfamily classifications. Comparative analysis re-confirmed the existing, BSPs-containing 14 genera and 21 species and displayed inter-genera similarities within existing subfamilies. Importantly, it also revealed the need for the creation of new taxonomic classifications, including 28 species, nine genera, and two subfamilies (New subfamily1 and New subfamily2) to accommodate inter-genera relatedness. Following pangenome analysis, no ortholog shared by all BSPs was identified, while orthologs, namely, the tail fibers/spike proteins and poly-gamma-glutamate hydrolase, that are shared by more than two-thirds of the BSPs were identified. More importantly, major capsid protein (MCP) type I, MCP type II, MCP type III and peptidoglycan binding proteins that are distinctive orthologs for Herelleviridae, Salasmaviridae, New subfamily1, and New subfamily2, respectively, were identified and analyzed which could serve as signatures to distinguish BSP members of the respective taxon.

CONCLUSIONS

In this study, we show the genomic diversity and propose a comprehensive classification of 61 BSPs, including the proposition for the creation of two new subfamilies, followed by the identification of orthologs taxonomic signature genes, potentially contributing to phage taxonomy.

The State of the Art in Biodefense Related Bacterial Pathogen Detection Using Bacteriophages: How It Started and How It's Going

Accurate pathogen detection and diagnosis is paramount in clinical success of treating patients. There are two general paradigms in pathogen detection: molecular and immuno-based, and phage-based detection is a third emerging paradigm due to its sensitivity and selectivity. Molecular detection methods look for genetic material specific for a given pathogen in a sample usually by polymerase chain reaction (PCR). Immuno-methods look at the pathogen components (antigens) by antibodies raised against that pathogen specific antigens. There are different variations and products based on these two paradigms with advantages and disadvantages. The third paradigm at least for bacterial pathogen detection entails bacteriophages specific for a given bacterium. Sensitivity and specificity are the two key parameters in any pathogen detection system. By their very nature, bacteriophages afford the best sensitivity for bacterial detection. Bacteria and bacteriophages form the predator-prey pair in the evolutionary arms race and has coevolved over time to acquire the exquisite specificity of the pair, in some instances at the strain level. This specificity has been exploited for diagnostic purposes of various pathogens of concern in clinical and other settings. Many recent reviews focus on phage-based detection and sensor technologies. In this review, we focus on a very special group of pathogens that are of concern in biodefense because of their potential misuse in bioterrorism and their extremely virulent nature and as such fall under the Centers for Disease and Prevention (CDC) Category A pathogen list. We describe the currently available phage methods that are based on the usual modalities of detection from culture, to molecular and immuno- and fluorescent methods. We further highlight the gaps and the needs for more modern technologies and sensors drawing from technologies existing for detection and surveillance of other pathogens of clinical relevance.

Classifying the Unclassified: A Phage Classification Method

This work reports the method ClassiPhage to classify phage genomes using sequence derived taxonomic features. ClassiPhage uses a set of phage specific Hidden Markov Models (HMMs) generated from clusters of related proteins. The method was validated on all publicly available genomes of phages that are known to infect Vibrionaceae. The phages belong to the well-described phage families of Myoviridae, Podoviridae, Siphoviridae, and Inoviridae. The achieved classification is consistent with the assignments of the International Committee on Taxonomy of Viruses (ICTV), all tested phages were assigned to the corresponding group of the ICTV-database. In addition, 44 out of 58 genomes of Vibrio phages not yet classified could be assigned to a phage family. The remaining 14 genomes may represent phages of new families or subfamilies. Comparative genomics indicates that the ability of the approach to identify and classify phages is correlated to the conserved genomic organization. ClassiPhage classifies phages exclusively based on genome sequence data and can be applied on distinct phage genomes as well as on prophage regions within host genomes. Possible applications include (a) classifying phages from assembled metagenomes; and (b) the identification and classification of integrated prophages and the splitting of phage families into subfamilies.

Isolation and Characterization of a Bacteriophage Preying an Antifungal Bacterium

Several Bacillus species were isolated from rice field soils, and 16S rRNA gene sequence analysis showed that Bacillus cereus was the most abundant. A strain named BC1 showed antifungal activity against Rhizoctonia solani. Bacteriophages infecting strain BC1 were isolated from the same soil sample. The isolated phage PK16 had an icosahedral head of 100 ± 5 nm and tail of 200 ± 5 nm, indicating that it belonged to the family Myoviridae. Analysis of the complete linear dsDNA genome revealed a 158,127-bp genome with G + C content of 39.9% comprising 235 open reading frames as well as 19 tRNA genes (including 1 pseudogene). Blastp analysis showed that the proteins encoded by the PK16 genome had the closest hits to proteins of seven different bacteriophages. A neighbor-joining phylogenetic tree based on the major capsid protein showed a robust clustering of phage PK16 with phage JBP901 and BCP8-2 isolated from Korean fermented food.

Genome Sequence of Bacillus cereus Phage vB_BceS-MY192

Bacillus cereus is an opportunistic foodborne pathogen. The phage vB_BceS-MY192 was isolated from B. cereus 192 in a cooked rice sample. The temperate phage belongs to the Siphoviridae family, Caudovirales order. Here we announce the phage genome sequence and its annotation, which may expand the understanding of B. cereus siphophages.

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.

Whole genome comparison of a large collection of mycobacteriophages reveals a continuum of phage genetic diversity

The bacteriophage population is large, dynamic, ancient, and genetically diverse. Limited genomic information shows that phage genomes are mosaic, and the genetic architecture of phage populations remains ill-defined. To understand the population structure of phages infecting a single host strain, we isolated, sequenced, and compared 627 phages of Mycobacterium smegmatis. Their genetic diversity is considerable, and there are 28 distinct genomic types (clusters) with related nucleotide sequences. However, amino acid sequence comparisons show pervasive genomic mosaicism, and quantification of inter-cluster and intra-cluster relatedness reveals a continuum of genetic diversity, albeit with uneven representation of different phages. Furthermore, rarefaction analysis shows that the mycobacteriophage population is not closed, and there is a constant influx of genes from other sources. Phage isolation and analysis was performed by a large consortium of academic institutions, illustrating the substantial benefits of a disseminated, structured program involving large numbers of freshman undergraduates in scientific discovery.

DOI:http://dx.doi.org/10.7554/eLife.06416.001

Viruses are unable to replicate independently. To generate copies of itself, a virus must instead invade a target cell and commandeer that cell's replication machinery. Different viruses are able to invade different types of cell, and a group of viruses known as bacteriophages (or phages for short) replicate within bacteria. The enormous number and diversity of phages in the world means that they play an important role in virtually every ecosystem.

Despite their importance, relatively little is known about how different phage populations are related to each other and how they evolved. Many phages contain their genetic information in the form of strands of DNA. Using genetic sequencing to find out where and how different genes are encoded in the DNA can reveal information about how different viruses are related to each other. These relationships are particularly complicated in phages, as they can exchange genes with other viruses and microbes.

Previous studies comparing the genomes—the complete DNA sequence—of reasonably small numbers of phages that infect the Mycobacterium group of bacteria have found that the phages can be sorted into ‘clusters’ based on similarities in their genes and where these are encoded in their DNA. However, the number of phages investigated so far has been too small to conclude how different clusters are related. Are the clusters separate, or do they form a ‘continuum’ with different genes and DNA sequences shared between different clusters?

Here, Pope, Bowman, Russell et al. compare the individual genomes of 627 bacteriophages that infect the bacterial species Mycobacterium smegmatis. This is by far the largest number of phage genomes analyzed from a single host species. The large number of genomes analyzed allowed a much clearer understanding of the complexity and diversity of these phages to be obtained. The isolation, sequencing and analysis of the hundreds of M. smegmatis bacteriophage genomes was performed by an integrated research and education program, called the Science Education Alliance Phage Hunters Advancing Genomics and Evolutionary Science (SEA-PHAGES) program. This enabled thousands of undergraduate students from different institutions to contribute to the phage discovery and sequencing project, and co-author the report. SEA-PHAGES therefore shows that it is possible to successfully incorporate genuine scientific research into an undergraduate course, and that doing so can benefit both the students and researchers involved.

The results show that while the genomes could be categorized into 28 clusters, the genomes are not completely unrelated. Instead, a spread of diversity is seen, as genes and groups of genes are shared between different clusters. Pope, Bowman, Russell et al. further reveal that the phage population is in a constant state of change, and continuously acquires genes from other microorganisms and viruses.

DOI:http://dx.doi.org/10.7554/eLife.06416.002

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.

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.

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.

Prevalence, genetic diversity, and host range of tectiviruses among members of the Bacillus cereus group

GIL01, Bam35, GIL16, AP50, and Wip1 are tectiviruses preying on the Bacillus cereus group. Despite the significant contributions of phages in different biological processes, little is known about the dealings taking place between tectiviruses and their Gram-positive bacterial hosts. Therefore, this work focuses on characterizing the interactions between tectiviruses and the B. cereus group by assessing their occurrence and genetic diversity and evaluating their host range. To study the occurrence of tectiviruses in the B. cereus group, 2,000 isolates were evaluated using primers designed to be specific to two variable regions detected in previously described elements. PCR and propagation tests revealed that tectivirus-like elements occurred in less than 3% of the isolates. Regardless of this limited distribution, several novel tectiviruses were found, and partial DNA sequencing indicated that a greater diversity exists within the family Tectiviridae. Analyses of the selected variable regions, along with their host range, showed that tectiviruses in the B. cereus group can be clustered mainly into two different groups: the ones infecting B. anthracis and those isolated from other B. cereus group members. In order to address the host range of some novel tectiviruses, 120 strains were tested for sensitivity. The results showed that all the tested tectiviruses produced lysis in at least one B. cereus sensu lato strain. Moreover, no simple relationship between the infection patterns of the tectiviruses and their diversity was found.