Novel Giant Siphovirus from Bacillus anthracis Features Unusual Genome Characteristics

Here we present vB_BanS-Tsamsa, a novel temperate phage isolated from Bacillus anthracis, the agent responsible for anthrax infections in wildlife, livestock and humans. Tsamsa phage is a giant siphovirus (order Caudovirales), featuring a long, flexible and non-contractile tail of 440 nm (not including baseplate structure) and an isometric head of 82 nm in diameter. We induced Tsamsa phage in samples from two different carcass sites in Etosha National Park, Namibia. The Tsamsa phage genome is the largest sequenced Bacillus siphovirus, containing 168,876 bp and 272 ORFs. The genome features an integrase/recombinase enzyme, indicative of a temperate lifestyle. Among bacterial strains tested, the phage infected only certain members of the Bacillus cereus sensu lato group (B. anthracis, B. cereus and B. thuringiensis) and exhibited moderate specificity for B. anthracis. Tsamsa lysed seven out of 25 B. cereus strains, two out of five B. thuringiensis strains and six out of seven B. anthracis strains tested. It did not lyse B. anthracis PAK-1, an atypical strain that is also resistant to both gamma phage and cherry phage. The Tsamsa endolysin features a broader lytic spectrum than the phage host range, indicating possible use of the enzyme in Bacillus biocontrol.     

Erratum to: 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.     

Complete genome sequence of Bacillus cereus bacteriophage PBC1

Bacillus cereus is a ubiquitous, spore-forming bacterium associated with food poisoning cases. To develop an efficient biocontrol agent against B. cereus, we isolated lytic phage PBC1 and sequenced its genome. PBC1 showed a very low degree of homology to previously reported phages, implying that it is novel. Here we report the complete genome sequence of PBC1 and describe major findings from our analysis.     

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.     

A mutant bacteriophage evolved to infect resistant bacteria gained a broader host range

Bacteriophages (phages) are the most abundant entities in nature, yet little is known about their capacity to acquire new hosts and invade new niches. By exploiting the Gram‐positive soil bacterium Bacillus subtilis (B. subtilis) and its lytic phage SPO1 as a model, we followed the coevolution of bacteria and phages. After infection, phage‐resistant bacteria were readily isolated. These bacteria were defective in production of glycosylated wall teichoic acid (WTA) polymers that served as SPO1 receptor. Subsequently, a SPO1 mutant phage that could infect the resistant bacteria evolved. The emerging phage contained mutations in two genes, encoding the baseplate and fibers required for host attachment. Remarkably, the mutant phage gained the capacity to infect non‐host Bacillus species that are not infected by the wild‐type phage. We provide evidence that the evolved phage lost its dependency on the species‐specific glycosylation pattern of WTA polymers. Instead, the mutant phage gained the capacity to directly adhere to the WTA backbone, conserved among different species, thereby crossing the species barrier.     

Pervasive prophage recombination occurs during evolution of spore-forming Bacilli

Phages are the main source of within-species bacterial diversity and drivers of horizontal gene transfer, but we know little about the mechanisms that drive genetic diversity of these mobile genetic elements (MGEs). Recently, we showed that a sporulation selection regime promotes evolutionary changes within SPβ prophage of Bacillus subtilis, leading to direct antagonistic interactions within the population. Herein, we reveal that under a sporulation selection regime, SPβ recombines with low copy number phi3Ts phage DNA present within the B. subtilis population. Recombination results in a new prophage occupying a different integration site, as well as the spontaneous release of virulent phage hybrids. Analysis of Bacillus sp. strains suggests that SPβ and phi3T belong to a distinct cluster of unusually large phages inserted into sporulation-related genes that are equipped with a spore-related genetic arsenal. Comparison of Bacillus sp. genomes indicates that similar diversification of SPβ-like phages takes place in nature. Our work is a stepping stone toward empirical studies on phage evolution, and understanding the eco-evolutionary relationships between bacteria and their phages. By capturing the first steps of new phage evolution, we reveal striking relationship between survival strategy of bacteria and evolution of their phages.Subject terms: Bacterial genetics, Evolution     

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.     

Characterization of Two Virulent Phages of Lactobacillus plantarum

We characterized two Lactobacillus plantarum virulent siphophages, ATCC 8014-B1 (B1) and ATCC 8014-B2 (B2), previously isolated from corn silage and anaerobic sewage sludge, respectively. Phage B2 infected two of the eight L. plantarum strains tested, while phage B1 infected three. Phage adsorption was highly variable depending on the strain used. Phage defense systems were found in at least two L. plantarum strains, LMG9211 and WCSF1. The linear double-stranded DNA genome of the pac-type phage B1 had 38,002 bp, a G+C content of 47.6%, and 60 open reading frames (ORFs). Surprisingly, the phage B1 genome has 97% identity with that of Pediococcus damnosus phage clP1 and 77% identity with that of L. plantarum phage JL-1; these phages were isolated from sewage and cucumber fermentation, respectively. The double-stranded DNA (dsDNA) genome of the cos-type phage B2 had 80,618 bp, a G+C content of 36.9%, and 127 ORFs with similarities to those of Bacillus and Lactobacillus strains as well as phages. Some phage B2 genes were similar to ORFs from L. plantarum phage LP65 of the Myoviridae family. Additionally, 6 tRNAs were found in the phage B2 genome. Protein analysis revealed 13 (phage B1) and 9 (phage B2) structural proteins. To our knowledge, this is the first report describing such high identity between phage genomes infecting different genera of lactic acid bacteria.     

A Virulent Phage Infecting Lactococcus garvieae,with Homology to Lactococcus lactis Phages

A new virulent phage belonging to the Siphoviridae family and able to infect Lactococcus garvieae strains was isolated from compost soil. Phage GE1 has a prolate capsid (56 by 38 nm) and a long noncontractile tail (123 nm). It had a burst size of 139 and a latent period of 31 min. Its host range was limited to only two L. garvieae strains out of 73 tested. Phage GE1 has a double-stranded DNA genome of 24,847 bp containing 48 predicted open reading frames (ORFs). Putative functions could be assigned to only 14 ORFs, and significant matches in public databases were found for only 17 ORFs, indicating that GE1 is a novel phage and its genome contains several new viral genes and encodes several new viral proteins. Of these 17 ORFs, 16 were homologous to deduced proteins of virulent phages infecting the dairy bacterium Lactococcus lactis, including previously characterized prolate-headed phages. Comparative genome analysis confirmed the relatedness of L. garvieae phage GE1 to L. lactis phages c2 (22,172 bp) and Q54 (26,537 bp), although its genome organization was closer to that of phage c2. Phage GE1 did not infect any of the 58 L. lactis strains tested. This study suggests that phages infecting different lactococcal species may have a common ancestor.     

Comparative whole genome analysis of six diagnostic brucellaphages

Whole genome sequencing of six diagnostic brucellaphages, Tbilisi (Tb), Firenze (Fz), Weybridge (Wb), S708, Berkeley (Bk) and R/C, was followed with genomic comparisons including recently described genomes of the Tb phage from Mexico (Tb_M) and Pr phage to elucidate genomic diversity and candidate host range determinants. Comparative whole genome analysis revealed high sequence homogeneity among these brucellaphage genomes and resolved three genetic groups consistent with defined host range phenotypes. Group I was composed of Tb and Fz phages that are predominantly lytic for Brucella abortus and Brucella neotomae; Group II included Bk, R/C, and Pr phages that are lytic mainly for B. abortus, Brucella melitensis and Brucella suis; Group III was composed of Wb and S708 phages that are lytic for B. suis, B. abortus and B. neotomae. We found that the putative phage collar protein is a variable locus with features that may be contributing to the host specificities exhibited by different brucellaphage groups. The presence of several candidate host range determinants is illustrated herein for future dissection of the differential host specificity observed among these phages.     

Development and application of CRISPR-based genetic tools in Bacillus species and Bacillus phages

Recently, the clustered regularly interspaced short palindromic repeats (CRISPR) system has been developed into a precise and efficient genome editing tool. Since its discovery as an adaptive immune system in prokaryotes, it has been applied in many different research fields including biotechnology and medical sciences. The high demand for rapid, highly efficient and versatile genetic tools to thrive in bacteria-based cell factories accelerates this process. This review mainly focuses on significant advancements of the CRISPR system in Bacillus subtilis, including the achievements in gene editing, and on problems still remaining. Next, we comprehensively summarize this genetic tool's up-to-date development and utilization in other Bacillus species, including B. licheniformis, B. methanolicus, B. anthracis, B. cereus, B. smithii and B. thuringiensis. Furthermore, we describe the current application of CRISPR tools in phages to increase Bacillus hosts' resistance to virulent phages and phage genetic modification. Finally, we suggest potential strategies to further improve this advanced technique and provide insights into future directions of CRISPR technologies for rendering Bacillus species cell factories more effective and more powerful.     

Analysis of six prophages in Lactococcus lactis IL1403: different genetic structure of temperate and virulent phage populations

We report the genetic organisation of six prophages present in the genome of Lactococcus lactis IL1403. The three larger prophages (36–42 kb), belong to the already described P335 group of temperate phages, whereas the three smaller ones (13–15 kb) are most probably satellites relying on helper phage(s) for multiplication. These data give a new insight into the genetic structure of lactococcal phage populations. P335 temperate phages have variable genomes, sharing homology over only 10–33% of their length. In contrast, virulent phages have highly similar genomes sharing homology over >90% of their length. Further analysis of genetic structure in all known groups of phages active on other bacterial hosts such as Escherichia coli, Bacillus subtilis, Mycobacterium and Streptococcus thermophilus confirmed the existence of two types of genetic structure related to the phage way of life. This might reflect different intensities of horizontal DNA exchange: low among purely virulent phages and high among temperate phages and their lytic homologues. We suggest that the constraints on genetic exchange among purely virulent phages reflect their optimal genetic organisation, adapted to a more specialised and extreme form of parasitism than temperate/lytic phages.     

Novel Virulent and Broad-Host-Range Erwinia amylovora Bacteriophages Reveal a High Degree of Mosaicism and a Relationship to Enterobacteriaceae Phages

A diverse set of 24 novel phages infecting the fire blight pathogen Erwinia amylovora was isolated from fruit production environments in Switzerland. Based on initial screening, four phages (L1, M7, S6, and Y2) with broad host ranges were selected for detailed characterization and genome sequencing. Phage L1 is a member of the Podoviridae, with a 39.3-kbp genome featuring invariable genome ends with direct terminal repeats. Phage S6, another podovirus, was also found to possess direct terminal repeats but has a larger genome (74.7 kbp), and the virus particle exhibits a complex tail fiber structure. Phages M7 and Y2 both belong to the Myoviridae family and feature long, contractile tails and genomes of 84.7 kbp (M7) and 56.6 kbp (Y2), respectively, with direct terminal repeats. The architecture of all four phage genomes is typical for tailed phages, i.e., organized into function-specific gene clusters. All four phages completely lack genes or functions associated with lysogeny control, which correlates well with their broad host ranges and indicates strictly lytic (virulent) lifestyles without the possibility for host lysogenization. Comparative genomics revealed that M7 is similar to E. amylovora virus ΦEa21-4, whereas L1, S6, and Y2 are unrelated to any other E. amylovora phage. Instead, they feature similarities to enterobacterial viruses T7, N4, and ΦEcoM-GJ1. In a series of laboratory experiments, we provide proof of concept that specific two-phage cocktails offer the potential for biocontrol of the pathogen.     

Novel Phage Group Infecting Lactobacillus delbrueckii subsp. lactis,as Revealed by Genomic and Proteomic Analysis of Bacteriophage Ldl1

Ldl1 is a virulent phage infecting the dairy starter Lactobacillus delbrueckii subsp. lactis LdlS. Electron microscopy analysis revealed that this phage exhibits a large head and a long tail and bears little resemblance to other characterized phages infecting Lactobacillus delbrueckii. In vitro propagation of this phage revealed a latent period of 30 to 40 min and a burst size of 59.9 ± 1.9 phage particles. Comparative genomic and proteomic analyses showed remarkable similarity between the genome of Ldl1 and that of Lactobacillus plantarum phage ATCC 8014-B2. The genomic and proteomic characteristics of Ldl1 demonstrate that this phage does not belong to any of the four previously recognized L. delbrueckii phage groups, necessitating the creation of a new group, called group e, thus adding to the knowledge on the diversity of phages targeting strains of this industrially important lactic acid bacterial species.     

Identification of Bacillus cereus Group Species Associated with Food Poisoning Outbreaks in British Columbia, Canada

Food poisoning laboratories identify Bacillus cereus using routine methods that may not differentiate all Bacillus cereus group species. We recharacterized Bacillus food-poisoning strains from 39 outbreaks and identified B. cereus in 23 outbreaks, B. thuringiensis in 4, B. mycoides in 1, and mixed strains of Bacillus in 11 outbreaks.     

Isolation and characterization of phages infecting Bacillus cereus

Aim: To isolate and characterize bacteriophages (phages) that infect the foodborne pathogen Bacillus cereus. Methods and Results: Two phages were isolated from soil based on their ability to form plaques on four indicator hosts including Bacillus thuringiensis subsp. israelensis, and three isolates of B. cereus. The purified phages were characterized by morphology, host range, single‐step growth curves and restriction enzyme digestion profiles. The phages appeared to be of the Myoviridae family based on their structure in electron micrographs. The phages lysed bacteria of several species, produced average burst sizes of 322 and 300 phages per infected cell, and both had genomes over 90 kb. The phages were chloroform‐resistant and stable at 4°C. They reduced the concentration of B. cereus in mashed potatoes by >6 log₁₀ CFU ml^?1 within 24 h at room temperature, when applied at a high concentration. Conclusions: The relatively narrow host range within B. cereus might mean that these phages need to be used as part of a ‘cocktail’ of phages for biocontrol, but their efficacy for the control of their host in food was demonstrated. Significance and Impact of the Study: This is the first report of biocontrol by phages of B. cereus in food.     

Alteration of Host Specificity to Lytic Bacteriophages in Streptococcus cremoris

A mutant of Streptococcus cremoris strain ML1 was isolated based on its resistance to acriflavine. The mutant strain showed resistance to the growth of virulent bacteriophages to which the parental strain was sensitive whereas it became sensitive to a number of other virulent phages to which the parental strain was resistant. At the same time, infection of the mutant strain by another bacteriophage sc607 resulted in killing of cells without production of progeny phages. The phage adsorption appeared normal, suggesting that the killing was a postadsorption event. Such killing of bacterial cells was prevented by chloramphenicol treatment, indicating that involvement of some protein either synthesized by phage or phage-induced cellular protein. Synthesis of ribonucleic acid was abruptly terminated after infection of the mutant strain by phage sc607 but not of the parental strain. The alteration of host specificity in the mutant to different lytic bacteriophages and especially abortive infection by phage sc607 resembles the prophage-mediated interference observed in other bacteria.     

Romulus and Remus,Two Phage Isolates Representing a Distinct Clade within the Twortlikevirus Genus,Display Suitable Properties for Phage Therapy Applications

The renewed interest in controlling Staphylococcus aureus infections using their natural enemies, bacteriophages, has led to the isolation of a limited number of virulent phages so far. These phages are all members of the Twortlikevirus, displaying little variance. We present two novel closely related (95.9% DNA homology) lytic myoviruses, Romulus and Remus, with double-stranded DNA (dsDNA) genomes of 131,333 bp and 134,643 bp, respectively. Despite their relatedness to Staphylococcus phages K, G1, ISP, and Twort and Listeria phages A511 and P100, Romulus and Remus can be proposed as isolates of a new species within the Twortlikevirus genus. A distinguishing feature for these phage genomes is the unique distribution of group I introns compared to that in other staphylococcal myoviruses. In addition, a hedgehog/intein domain was found within their DNA polymerase genes, and an insertion sequence-encoded transposase exhibits splicing behavior and produces a functional portal protein. From a phage therapy application perspective, Romulus and Remus infected approximately 70% of the tested S. aureus isolates and displayed promising lytic activity against these isolates. Furthermore, both phages showed a rapid initial adsorption and demonstrated biofilm-degrading capacity in a proof-of-concept experiment.     

Beyond arbitrium: identification of a second communication system in Bacillus phage phi3T that may regulate host defense mechanisms

The evolutionary stability of temperate bacteriophages at low abundance of susceptible bacterial hosts lies in the trade-off between the maximization of phage replication, performed by the host-destructive lytic cycle, and the protection of the phage-host collective, enacted by lysogeny. Upon Bacillus infection, Bacillus phages phi3T rely on the “arbitrium” quorum sensing (QS) system to communicate on their population density in order to orchestrate the lysis-to-lysogeny transition. At high phage densities, where there may be limited host cells to infect, lysogeny is induced to preserve chances of phage survival. Here, we report the presence of an additional, host-derived QS system in the phi3T genome, making it the first known virus with two communication systems. Specifically, this additional system, coined “Rapφ-Phrφ”, is predicted to downregulate host defense mechanisms during the viral infection, but only upon stress or high abundance of Bacillus cells and at low density of population of the phi3T phages. Post-lysogenization, Rapφ-Phrφ is also predicted to provide the lysogenized bacteria with an immediate fitness advantage: delaying the costly production of public goods while nonetheless benefiting from the public goods produced by other non-lysogenized Bacillus bacteria. The discovered “Rapφ-Phrφ” QS system hence provides novel mechanistic insights into how phage communication systems could contribute to the phage-host evolutionary stability.Subject terms: Bacteriophages, Viral genetics