- Open Access
Genome sequence and analysis of a broad-host range lytic bacteriophage that infects the Bacillus cereus group
Virology Journal volume 10, Article number: 48 (2013)
Comparatively little information is available on members of the Myoviridae infecting low G+C content, Gram-positive host bacteria of the family Firmicutes. While numerous Bacillus phages have been isolated up till now only very few Bacillus cereus phages have been characterized in detail.
Here we present data on the large, virulent, broad-host-range B. cereus phage vB_BceM_Bc431v3 (Bc431v3). Bc431v3 features a 158,618 bp dsDNA genome, encompassing 239 putative open reading frames (ORFs) and, 20 tRNA genes encoding 17 different amino acids. Since pulsed-field gel electrophoresis indicated that the genome of this phage has a mass of 155-158 kb Bc431v3 DNA appears not to contain long terminal repeats that are found in the genome of Bacillus phage SPO1.
Bc431v3 displays significant sequence similarity, at the protein level, to B. cereus phage BCP78, Listeria phage A511 and Enterococcus phage ØEF24C and other morphologically related phages infecting Firmicutes such as Staphylococcus phage K and Lactobacillus phage LP65. Based on these data we suggest that Bc431v3 should be included as a member of the Spounavirinae; however, because of all the diverse taxonomical information has been addressed recently, it is difficult to determine the genus. The Bc431v3 phage contains some highly unusual genes such as gp143 encoding putative tRNAHis guanylyltransferase. In addition, it carries some genes that appear to be related to the host sporulation regulators. These are: gp098, which encodes a putative segregation protein related to FstK/SpoIIIE DNA transporters; gp105, a putative segregation protein; gp108, RNA polymerase sigma factor F/B; and, gp109 encoding RNA polymerase sigma factor G.
The genus Bacillus can be subdivided into two groups: the Bacillus subtilis group and the Bacillus cereus group. Members of the latter subdivision include B. anthracis, the causative agent of the fatal human and animal disease, anthrax, B. cereus, B. weihenstephanensis, B. thuringiensis, B. mycoides and B. pseudomycoides. B. cereus and B. weihenstephanensis are implicated in foodborne illnesses as well as food spoilages. Their psychrotrophic properties enable them to cause problems in foods stored at low temperatures especially dairy products[3–6]. B. thuringiensis, B. mycoides and B. pseudomycoides have also been noted to cause food spoilage. In addition, there is some evidence to indicate they potentially can cause foodborne illnesses[1, 8–11].
With the exception of phages AP50, AP50-04, AP50-11, AP50-23, AP50-26, AP50-27, and Bam35 which belong to the family of Tectiviridae[12, 13], all other Bacillus phages belong to the three tailed phage families: Myoviridae, Siphoviridae and Podoviridae. They range from very large viruses such as the G, to very small phages such as ϕ29 and ϕ15. All of the Bacillus phages possess dsDNA, and some of them are quite distinct as they contain unusual bases in their DNA. In the case of phage SPO1, thymine residues are replaced with 5-(hydroxymethyl) uracil (HmUra); while in phage PBS2 DNA the thymine are totally replaced by uracil. Among the myoviruses, the subfamily Spounavirinae (SPO plus “una”, Latin for “one”) is comprised of two genera, “SPO1-like viruses” and “Twort-like viruses”. Members of Spounavirinae feature isometric heads 87–94 nm in diameter and tails 140–219 nm long[15–17]. The dsDNA molecules of this subfamily are quite large (127–142 kb) and nonpermuted with 3.1 – 20 kb terminal redundancies. The genus “Spo1likevirus” contains at least ten Bacillus phages, Lactobacillus phage 222a as well as Enterococcus faecalis phage ØEF24C with only the genomes of SPO1 and ØEF24C phages having been sequenced[16, 18–21]. The distinct criterion of their genomes (i.e., “Spo1likevirus”) is the substitution of thymine with HmUra, and the presence dUMP hydroxymethylase activity. The “Twortlikevirus”, on the other hand, contains a group of virulent phages of Staphylococcus including Twort, G1 and K phages as well as Listeria phages A511 and P100[22–25]. As far as we know the DNA of these phages is not modified.
It is worthwhile mentioning that because of the close relatedness between members of B. cereus group; most of the phages isolated against B. anthracis were also able to infect B. cereus and vice-versa. The first phage reported against B. anthracis was phage W, which was originally induced from a soil strain of B. cereus. A lytic variant of W, phage γ, was isolated shortly after and it is now used as a standard protocol for B. anthracis identification by the Centers for Disease Control and Prevention and other public health laboratories in the United States[27, 28].
Very few phages of B. cereus have been isolated and characterized in detail. Some temperate phages were characterized and utilized in the transductional analyses of B. cereus strains[25, 29–33]. In addition, a number of lytic phages have been isolated and used in phage typing schemes for epidemiological studies. It is noteworthy that phage typing could be considered as an efficient tool especially because of the close relatedness between B. cereus, B. thuringiensis and B. anthracis, providing a cheap, convenient and fairly accurate tool to identify such closely related species. Lytic phages have also been used in the control of B. cereus in mashed potatoes[25, 35, 36].
In the present study, we isolated a lytic phage, vB_BceM_Bc431v3, that showed a strong potential to control B. cereus. Prior to its application in food preservation, we needed to characterize it and confirm that it will not affect the pathogenicity of this bacterium. Therefore, the objectives of this study were to determine the genomic sequence of phage Bc431v3, annotate its genes and appraise its structural functions in conjunction with existing phages present in the scientific literature and in the GenBank.
Isolation and morphology of Bc431v3
Phage Bc431v3 was isolated from sewage and produced small (1.8 mm) clear plaques with turbid borders on B. cereus strain LJH431. This virus was negatively stained with 2% uranyl acetate and examined by transmission electron microscopy. The micrographs show that this phage has isometric heads 85.4±3 nm in diameter with individual capsomers visible. The virus possesses a long contractile tail 180±3 nm in length by 12±4 nm in width. The base plate has a cluster of projections and what appears to be a central tail fibre. Collectively these features indicate that this virus belongs to the family Myoviridae (Figure1).
In a host range study, Bc431v3 was able to infect all tested B. cereus, B. anthracis, B. licheniformis and B. weihenstephanensis strains with different degrees of lysis (Table1; Figure2). Phage Bc431v3 also infected B. thuringiensis, B. psychrosaccharolyticus and B. megaterium, but was not able to lyse a strain of B. subtilis.
One-step growth curve
The growth characteristics of this virus in B. cereus strain LJH431 revealed a latent period of 85 ± 5 min and a burst size of 318 ± 5 (Figure3).
Genome of phage Bc431v3
The phage genome was sequenced to 19-fold coverage using 454 (pyrosequencing) technology. The genome is 158,618 bp long and possesses a G+C content of 39%. BLASTX analysis against the nonredundant NCBI database failed to reveal any frameshifts. Since pulsed-field gel electrophoresis indicated that the phage genome is a single linear DNA molecule of 155–158 kb it would suggest that this phage does not possess long terminally redundant ends such as are present on Bacillus phage SPO1.
Identification and analysis of open reading frames (ORFs)
The genome was initially analyzed using AutoFACT, with the ORF calls verified using Kodon and BLAST. A total of 239 putative ORFs were identified in the genome (Figure4; Additional file1: Table S1). A total of 143,817 bp nucleotides (90.7% of the genome) are involved in coding for putative proteins. Of these, 38 ORFs (15.9%) had identified functions, while 76 ORFs (31.8%) were identified as possessing homologs with proteins in the nonredundant NCBI database; however, their functions are undefined (conserved hypothetical proteins). A large percentage (125 ORFs; 52.3%) were considered to be hypothetical proteins unique to this phage. Using the homology approach introduced by Abbasifar et. al., no toxins were identified among the proteins specified by this phage. Three different start codons utilized in the genome; ATG, GTG and TTG were used at frequencies 88.7%, 5.9% and 5.5%, respectively.
The genome of this phage was exposed to DraI and SspI, EcoRI and NdeI which resulted in poor digestion, while HindIII and AccI failed to digest (data not shown). In silico restriction endonucleases digestion of phage Bc431v3 DNA revealed that the DNA possessed 200 AccI sites, 28 NdeI sites and 23 HindIII sites, suggesting that the DNA of this virus may be modified. Two C-5 cytosine-specific DNA methyltransferases were identified in the genome encoded by genes 172 and 173, but their impact on the digestibility of the DNA is unknown.
In the following sections the properties of some of the genes will be discussed.
Nucleotide metabolism and DNA replication
Sequence-based predictions identified many genes as being involved in Bc431v3 nucleotide metabolism and DNA synthesis. The former included thymidylate synthase (gp021), ribonucleotide reductase alpha (gp185) and beta subunits (gp183), dihydrofolate reductase (gp024), deoxyuridine 5′-triphosphate nucleotidohydrolase (gp190), an exonuclease I homolog (gp194), and a putative exonuclease II (gp195). The proteins involved in DNA replication included a DNA polymerase (gp170), DNA primase (gp192), and two DNA helicases (gp197 and gp200).
DNA packaging and morphogenesis – genomic and proteomic analysis
Several genes encoding proteins directly involved in DNA packaging and morphogenesis were identified in Bc431v3 genome. Terminase small and large subunits homologs were determined to be the products of genes 010 and 011, respectively. Genes encoding capsid morphogenesis were determined to be 232, 235 and 236 for major head, prohead protease, and portal vertex proteins, respectively. Tail proteins were expressed by a set of genes, which are 224 (tail core protein), 225 (tail sheath protein), 203 and 215 (putative tail proteins), 205 and 206 (baseplate proteins); and, 204 for putative tail fiber adapter protein. It is highly unusual for the terminase genes to be separated from the capsid-tail gene complex.
A cluster of genes starting from orf203 to orf230 carried the genes that are involved in morphogenesis, and their products were structurally similar to B. cereus phage Bcp78 and Listeria phage A511[23, 25]. Additionally, it was found that the product of orf218 is a putative tail protein possessing endo-beta-N-acetylglucosaminidase activity that also contains three transmembrane domains. BLASTP homology showed that the gp218 has the motif of C-terminus similar to the putative tail lysin of the Enterococcus phage ØEF24C as well as the pfam lysozyme sub-family 2. Additionally, gp218 also showed a homology with Listeria phage A511 gp97, a putative peptidase.
To achieve reliable protein identification, the purified phage were digested by two different enzymes (trypsin and chymotrypsin) and analyzed by high resolution UPLC and LTQ-FT mass spectrometry in duplicates. The combined dataset from Mascot database search identified 45 phage proteins, in which the top 22 proteins were present at high scores and sequence coverages (Additional file2: Table S2). As expected, this short list included the highly abundant structural proteins of phage Bc431v3: major capsid membrane protein (orf232; sequence coverage 50%), tail protein (orf203; 13%), tail sheath protein (orf225; 63%), tail fiber protein (orf215; 14%), tail protein possessing endo-beta-N-acetylglucosaminidase activity (orf218; 35%), baseplate protein (orf205; 23%), and related structural protein (orf166; 39%). Other phage interacting proteins and enzymes were also detected such as N-acetylmuramoyl-L alanine amidase (orf012; coverage 87%) and prohead protease (orf235, 12%). Prohead protease is a proteolytic enzyme essential for phage capsid morphogenesis, and N-acetylmuramoyl-L alanine amidase plays an important role in DNA synthesis through cleavage of amide bonds between N-acetylmuramoyl and L–amino acids in the bacteriophage cell walls.
We know relatively little about the regulation of transcription in large Bacillus phages. Using a combination of two techniques, visual scanning for sequences similar to SigA promoters (Sigma43; TTGACA(−35)-N14-tgnTATAAT(−10);) and MEME analysis we identified 27 putative promoters, 12 of which contained the extended −10 region (Additional file3: Table S3). All of these sequences possessed typical AT-rich UTR sequences, but atypically the −35 region has the consensus TTGTTGAC, not TTGACA as would be expected for promoters recognized by SigA. The majority of putative promoters were located upstream of genes involved in nucleotide metabolism (e.g. tRNAHis guanylyltransferase, dihydrofolate reductase) or DNA synthesis (DNA polymerase, primase, polymerase) which, based upon analysis of other phages, correspond to middle or delayed early genes. Because of the lack of homologs we cannot define the early genes of this phage. One of the middle genes (orf163) specifies a protein which is related to proteins defined as sigma factors for Enterococcus phage phiEF24C (YP_001504175), Staphylococcus phage K (YP_024522), Bacillus phages SP10 (BAK53012) and SPO1 (YP_002300430). In the latter virus, gp34 is known to facilitate, along with gp33 late transcription. Based upon similarity to CGTTAGA(N17-19)GATATT and allowing for 2 bp mismatch we identified four putative late promoters upstream of genes 056, 162, 190 and 198. Examination product of Bc431v3 orf164 using HHPred, failed to demonstrate homology to SPO1 gp33.
Interestingly, the Bc431v3 genome, like that of temperate phage Bcp1, was found to carry two other genes encoding for bacterial-type sigma factors. These are the products of genes 108 and 109, which appear to be members of the RNA polymerase sigma F/B and sigma factor G subfamilies, respectively. MEME analysis using a variety of Bacillus, Acyclobacillus, Geobacillus and Paenibacillus sigma factors revealed that gp108 shared the following motif Wx(9)Ix(2)Lx(3)Ex(2)Ix(6)KDx(2)QSx (2)Ax(2)LGx(5)V with SigmaF-like proteins, while gp109 shared DDxFQxGxIGLx(3)Ix(2)FDx(6)FSTYAV with related SigG proteins in the Bacillaceae.
The final stage of the phage lytic cycle is degradation of the bacterial cell wall and release of progeny phages. The lysis of the cell wall is typically induced by two phage-encoded proteins, holin and endolysin[41, 42]. Holin forms a hole in the cell membrane, and endolysin passes through the hole and destroys the peptidoglycan structure[41, 42]. Two gene products were determined to be the possible lysin, an N-acetylmuramoyl-L-alanine amidases encoded by gene 012 and/or L-alanoyl-D-glutamate peptidase encoded by gene 222. The former contains a MurNAc-LAA [cd02696], N-acetylmuramoyl-L-alanine amidase domain while gp222 contains a COG3584 (uncharacterized protein conserved in bacteria) domain. Both share homology to defined phage lysins. As for the holin, while ten proteins possessed two transmembrane domains, and two three domains, none possessed the small size generally associated with holins.
Codon usage and tRNAs
The Bc431v3 genome was found to contain 20 tRNA genes for 17 amino acids (Additional file4: Table S4). By comparing the codon usage pattern of the phage with that of its host (Bacillus cereus) we were able to see that in ten cases the phage-encoded tRNAs may significantly enhance translation of phage mRNAs (Additional file5: Table S5). The products of two genes 143 and 177 appear to encode proteins, a tRNAHis guanylyltransferase and predicted a nucleotidyltransferase, respectively, which may modify the phage or host tRNAs. As is quite common, this virus also includes a putative tRNAMet, which TFAM identified as a selenocysteinyl tRNA.
The phage’s Bc431v3 genome was found to contain several genes that have rarely been detected in other phages. In addition to sigma factors which are related to similar host proteins, we identified a putative segregation protein (gp098) related to FtsK-SpoIIIE, a DNA-binding protein (gp174) related to integration host factor (IHF), and a putative segregation protein gp105.
Using BLASTN the sequence of Bc431v3 showed significant sequence similarity to the terminases of B. cereus phage Bcp1, and B. anthracis phages PlyM19 and 1102ϕ3-1. At the protein level TBLASTX analysis (Figure4) revealed similarity between Bc431v3 and Listeria phage A511 and Enterococcus phage ØEF24C localized to regions: 5′-end to 20 kb, 60–120 kb, and 130 kb - 3′-end. The 20-60 kb regions encode a wide range of hypothetical and conserved hypothetical proteins, plus the FtsK-SpoIIIE, and two sporulation-specific sigma factors, while the 120-130 kb region species a number of minor structural proteins. More detailed analysis using CoreGenes[47, 48] showed that this phage shared almost 36.4% sequence homology with Listeria phage A511 and Enterococcus phage ØEF24C while it only showed 24.3% sequence homology with Bacillus phage SPO1. In a recent reclassification of the Myoviridae, the subfamily Spounavirinae was proposed with two genera the “Spo1likevirus” and the “Twortlikevirus.” Recently a considerable number of new complete genomes of Bacillus and Staphylococcus genomes have been deposited in GenBank. Rather than verifying the ICTV taxonomy proposal (2009.009a-pB.A.v3) or the paper by Klumpp et al. the new data, as shown through CoreGenes proteomic analyses and phylogeny (Figure5) strongly suggests that the Spounavirinae are far more diverse than originally realized. It would, we think be premature to assign Bc431v3 to any genus in the Spounavirinae.
In this study we studied the biological, genomic and proteomic characteristics of Bc431v3, a strong lytic phage that is capable infecting a broad range of Bacillus species in general and members of B. cereus group in particular. This virus lacked repressor, site-specific integrase, virulence and antibiotic resistance determinants, which increase its potential application to the biocontrol of members of the B. cereus group.
One of the unique features of this virus is the presence of several rare or unique genes. Gene 174 encodes a putative DNA-binding protein related to bacteria integration host factor (IHF). The latter protein is nucleoid-associated and is implicated in a number of chromosomal functions including DNA compaction. Bc431v3 also carries several genes directly related to host sporulation regulators, such as gene 105 encoding a putative segregation protein related to the FstK/SpoIIIE family of DNA transporters (DNA translocase). DNA transporters are responsible for the translocation of chromosomes after the completion of genome division. In Bacillus, the DNA translocase SpoIIIE is involved in transporting miss-segregated DNA during vegetative growth[50, 51] and also plays an essential role in chromosome segregation during the process of sporulation. Additionally, genes 108 and 109 encode sigma factor F/B σF/B and sigma factor σG, respectively. After verification of protein motifs with MEME software, we found that gp108 displays >70% sequence similarity to Bacillus sigma factors F and B while gp109 is > 80% related to Bacillus sigma factor G. These proteins are both forespore transcription factors functioning respectively early and late in sporulation. Occurrence of sporulation-related transcriptional factors on the phage Bc431v3 genome raises intriguing questions about their role in phage transcription. Schuch and Fischetti noted similar sigma factors F and G genes on the temperate phage Bcp1 genome that infects B. anthracis.
Materials and methods
Bacterial strains and culture media
Bacillus cereus strain LJH431 used in this study was obtained from the Canadian Research Institute for Food Safety culture collection (CRIFS, University of Guelph, Guelph, ON, Canada). Another fifteen strains of B. cereus were provided generously by T. Abee (Laboratory of Food Microbiology, Wageningen University, Netherlands). Strains A2 and D2 of B. cereus were kindly provided by A. Hudson (Christchurch Science Centre, Institute of Environmental Science and Research, New Zealand). Strains of B. weihenstephanensis were kindly provided by S. Scherer (Department of Biosciences, Wissenschaftszentrum Weihenstephan, Germany). Finally B. anthracis ΔSterne was obtained from K. Amoako (Canadian Food Inspection Agency, Ottawa, ON).
B. cereus strain LJH431 was grown in half-strength media of Tryptic Soy Broth (TSB; Difco Laboratories, Detroit, MI), half-strength Tryptic Soy Agar (TSA: TSB plus 1.5% agar) and half-strength Tryptic Soy Soft Agar (TSB + 0.5% agar) at 30°C for 24 h.
Transmission electron microscopy
The morphology of the phage Bc431v3 was examined by transmission electron microscopy as described. Briefly, 1 ml of each phage (~ 9 log10 PFU ml-1) was centrifuged at 16,000 × g for 1 h at 4°C (Beckman J-20 centrifuge, Beckman Coulter Inc., Mississauga, ON), pellets were washed once using CM buffer (2.5 g/L MgSO4.7H2O; 0.05 g/L gelatin; 6 ml/L 1 M Tris buffer; 0.735 g/L CaCl2; pH 7.5) and then resuspended in 300 μl of CM buffer. Five microliters of the phage suspension were applied onto 300-mesh copper grids coated with formvar and allowed to stand for 1 min. The excess liquid was removed by filter paper and grids were stained with 2% uranyl acetate for 30 s then carefully blotted with filter paper to remove excess stain solution. Negatively stained phages were examined with a LEO 912AB transmission electron microscope (Energy filtered TEM [EFTEM], LEO 912AB model at 100kv, Zeiss, Germany).
Phage isolation and purification and, extraction of DNA
Phage vB_BceM_Bc431v3 was isolated from sludge samples collected from a local wastewater management plant, and plaque purified. Lysates were prepared using the solid propagation method as described elsewhere, clarified from bacterial cell debris by centrifugation (5,400× g, 15 min) and filtered through 0.45 μm mixed cellulose ester (MCE) Fisher brand syringe filters (Fisher Scientific Company, Ottawa, ON, Canada). DNase I and RNase A (Sigma-Aldrich, Oakville, ON, Canada) were added to 10 μg/ml and held at 37°C for 30 min. Subsequent DNA purification steps were carried on according to the handbook of QIAGEN® Lambda midi Kit (Qiagen Inc., Mississauga, ON, Canada). DNA concentration was determined by absorbance at 260 nm using a Nanodrop™ 1000 spectrophotometer (Thermo Scientific, Ottawa, ON, Canada).
Determination of host specificity using Bioscreen C technology
The host range of the phage Bc431v3 on the selected strains of Bacillus spp. was determined using the spot test technique and also by measuring the optical density (OD) in liquid medium of the tested bacterium in the presence of phage using the Bioscreen C Microbiology Plate Reader (Labsystems, Helsinki, Finland) as described elsewhere. The following settings were used: wide band (wb) wavelength; 25°C incubation temperature; 10 min preheating time; kinetic measurement; measurement time 24 h at time intervals of 30 min with medium intensity shaking for 10 s before and after each measurement. Fifty microliters of each phage lysate were transferred to each of the 100 wells of the sterilized honeycomb plates of the Bioscreen C reader (Fischer Scientific), and then each of the wells was inoculated with 125 μl of an overnight bacterial culture at final concentrations of 3 log10 CFU ml-1. The multiplicity of infection (M.O.I) used to determine the host range of the phage was around 103. In this experimental design, three types of control samples were used: phage only in broth, bacteria only and sterile medium. The absorbance data were analyzed using the Bioscreen C data processing software version 5.26 (Labsystems) to determine the detection time (time required for each test well to increase by 0.3 OD units). Detection times (h:min) were converted to decimal values, averaged and the mean control detection time was subtracted from all test data for each isolate tested and expressed as detection time difference (DT diff.). Based on the host specificity, phages were divided into three groups: (++), when phages completely inhibited the growth of the host bacterium; (+), indicating that phages delayed the growth of the host bacterium; and (−), indicating that the phage had no effect on bacterial growth.
One-step growth curve
Burst sizes and latent periods of the selected phages were determined by a one-step growth experiment as described by Anany. Phages were added to its host bacterium at an MOI of around 0.1 and incubated in a water bath at 30°C for 5 min. One ml was removed and added to 100 μl of chloroform and mixed well. One hundred microliters of this mixture were added to 100 μl of an overnight culture of the host bacterium and mixed with 4 ml of overlay media and poured onto TSA agar plates to determine the degree of adsorption of the phage to bacterial cells. After an additional 30 s at 30°C, 100 μl were transferred to a tube containing 9.9 ml of fresh TSB and then diluted 10 times in fresh TSB (0.1 log10 ml-1). One ml of the 10-1 dilution tube was additionally diluted 10 times in fresh TSB (0.01 log10 ml-1) and all three tubes (the original plus 10-1 and 10-2 dilutions) were incubated in a water bath at 30°C. After 6 min, samples were collected every 5 min for 3 h and phages were titrated in each respective sample as previously described. The relative burst size was determined according to the equation:
Relative burst size = [(Final titre – Initial titre)/ Initial titre]
The relative burst size at different times was plotted against time to determine the latent period.
Determination of phage genome size using PFGE
Phage genome size was determined by pulsed-field gel electrophoresis (PFGE) as described elsewhere. Phage particles were embedded in 1% Seakem Gold agarose (Mandel Scientific, Guelph, ON) and subjected to electrophoresis in 0.5X TBE buffer (5X: 20 ml of 0.5 M EDTA [pH 8.0], 53 g/l Tris base, and 27.5 g/l boric acid) at 14°C for 18 h, using a CHEF DR-III Mapper electrophoresis system (Bio-Rad, Mississauga, ON) with pulse times of 2.2-54.2 s, at 6 V/cm. Low range DNA marker and phage lambda DNA concatemers (New England Biolabs) were used as size standards. The gels stained with ethidium bromide and DNA bands were visualized under UV transillumination. PFGE results were analyzed using BioNumerics software (Applied Maths Inc., Austin, TX).
Genome sequencing and annotation
The sequencing of phage Bc431v3 DNA was carried out at the McGill University Genome Quebec Innovation Centre (Montreal, QC, Canada) using pyrosequencing (454 technology). AutoFACT automated annotation software was initially used for genome annotation and then all open reading frames (ORFs) were confirmed using Kodon version 2.0 (Applied Maths). The individual proteins were analyzed using BLASTP against the protein databases at NCBI (http://www.ncbi.nlm.nih.gov). Protein motifs structures were identified using Pfam 24 (http://pfam.sanger.ac.uk/). In the case of the phage sigma factors a motif analysis was carried out using MEME v.4.5.0. Phage-encoded tRNA genes were identified with Aragorn and tRNAscan-SE v.1.21 using the default parameters. DNAMAN (Lynnon Corp., Vaudreuil-Dorion, QC, Canada) was used to determine the codon usage information of both Bc431v3 and its bacterial host Bacillus cereus, which also provided information on the GC content and presence of direct repeats in Bc431v3 genome. Putative promoters were identified by visual inspection for sequence similarity to TTGACA-N15-18-TATAAT, and by MEME analysis on the 5′ sequences extracted using extractUpStreamDNA (http://lfz.corefacility.ca/extractUpStreamDNA/). Rho-independent terminators were determined by ARNold and TransTerm software[60, 61] and verified by examining the secondary structure of the DNA adjacent to polyT sequences using Mfold. The physicochemical parameters of the gene products were determined using Molecule Weight and Isoelectric Point Finder (http://greengene.uml.edu/programs/FindMW.html). Transmembrane domains were predicted using TMHMM v2.0 and Phobius or Split 4.0[63–65]. Sequences of bacterial sigma factors and sigma-factor binding sites were identified through DBTBS database (http://dbtbs.hgc.jp).
Genomic comparisons at the proteomic level were made using CoreGenes[47, 48]. For alignments of multiple genomes and defining sequence homology percentage with related phages, progressive Mauve was used. For genomic map visualization and annotation pipelines, the CGview software was used (http://www.wishart.biology.ualberta.ca/cgview/).
GenBank accession number
The sequence of genome of this phage has been deposited with GenBank under accession number JX094431.
Phage Bc431v3 purified through CsCl gradients, was reduced with 10 mM dithiothreitol (56°C, 1 hr) and alkylated by 55 mM iodoacetamide (room temperature, dark, 1 hr), and then dialyzed against 10 mM NH4HCO3, and dried by SpeedVac concentrator (Savant, Fisher scientific, Nepean, Ontario). Enzymatic digestions were performed on ~20 μg of the purified protein using either sequencing grade trypsin or chymotrypsin (100 ng, Roche Diagnostics GmbH, Indianapolis, IN) for 4 hours. The digests were subsequently diluted by 0.2% formic acid and analyzed by online nanoAcquity ultra-performance liquid chromatography (UPLC, Waters, Milford, MA ) coupled with linear ion-trap Fourier transform ion cyclotron resonance (LTQ-FT ICR, Thermo Fisher, San Jose, CA) mass spectrometry. Peptides were trapped by a RP Symmetry C18 column (180 μm i.d. × 20 mm length, 5 μm) at 5 μl/min, and subsequently separated on a C18 analytical column (100 μm i.d. × 100 mm, 1.7 μm, BEH 130) at 400 nl/min. Peptide elution was achieved using mobile phases consisting of solvent A (0.1% FA) and solvent B (acetonitrile/0.1% FA) at a linear gradient from 5% to 30%, and then 85% of solvent B (65 min run). FT-MS scans were acquired with high resolution (100,000) at the mass range of m/z 300 to 2000, and low resolution MS/MS measurements in linear ion-trap mode were obtained by data-dependent scans of the top eight most intense precursor ions at multiply charged states of 2+, 3+, and 4+. Dynamic exclusion was enabled for a period of 180 S.
Protein identification was performed using an in-house Mascot Server (version 2.3.0, Matrix Science, London, UK), and the raw data were searched against the Bc431v3 protein database. The parameter settings allowed specific trypsin digestion for maximum 2 missed cleavage sites, and non-specific digestion of chymotrypsin. Cystein carbamidomethylation was designated as a fixed modification of peptides, and deamidation of asparagine and glutamine, methionine oxidation, pyro-Glu of Gln conversion at the N-terminus were considered as variable modifications. Mass tolerances were set up to 10 ppm for the FT MS ions and 1 Da for ion trap MS/MS fragment ions. Peptide assignments were filtered by an ion score cut off of 20, and the significance threshold was adjusted to 0.001 to achieve a false discovery rate (FDR) of less than 3%.
TFEA is currently a post-doctoral fellow in Canadian Research Institute for Food Safety, Food Science Department, University of Guelph, Ontario, Canada.
Basic Local Alignment Search Tool
Double-stranded deoxyribonucleic acid
Multiplicity of infection (ratio of phage to bacteria)
National Center for Biotechnology Information
Open reading frames
Pulsed-field gel electrophoresis
Plaque forming unit
Tris borate EDTA
Tryptic Soy Broth
TransMembrane prediction using Hidden Markov Models
Griffiths MW: Bacillus cereus and Other Bacillus spp. In Pathogens and Toxins in Foods: Challenges and Interventions. Edited by: Juneja VK, Sofos JN. Washignton, DC: ASM Press; 2010.
Garber KB: Anthrax (Bacillus anthracis). KGaA: Molecular Biology of. Wiley-VCH Verlag GmbH & Co; 2006.
Bandara N, Jo J, Ryu S, Kim KP: Bacteriophages BCP1-1 and BCP8-2 require divalent cations for efficient control of Bacillus cereus in fermented foods. Food Microbiol 2012,31(1):9-16. 10.1016/j.fm.2012.02.003
Bottone EJ: Bacillus cereus, a volatile human pathogen. Clin Microbiol Rev 2010, 23: 382-398. 10.1128/CMR.00073-09
Shin H, Bandara N, Shin E, Ryu S, Kim KP: Prevalence of Bacillus cereus bacteriophages in fermented foods and characterization of phage JBP901. Res Microbiol 2011,162(8):791-797. 10.1016/j.resmic.2011.07.001
Stenfors L, Mayr R, Scherer S, Granum P: Pathogenic potential of fifty Bacillus weihenstephanensis strains. FEMS Microbiol Lett 2002, 215: 47-51. 10.1111/j.1574-6968.2002.tb11368.x
Barfod K, Poulsen S, Hammer M, Larsen S: Sub-chronic lung inflammation after airway exposures to Bacillus thuringiensis biopesticides in mice. BMC Microbiol 2010, 10: 233. 10.1186/1471-2180-10-233
Ankolekar C, Rahmati T, Labbé RG: Detection of toxigenic Bacillus cereus and Bacillus thuringiensis spores in U.S. rice. Int. J. Food Microbiol 2009, 128: 460-466. 10.1016/j.ijfoodmicro.2008.10.006
McIntyre L, Bernard K, Beniac D, Isaac-Renton JL, Naseby DC: Identification of Bacillus cereus group species associated with food poisoning outbreaks in British Columbia. Canada Appl Environ Microbiol 2008, 74: 7451-7453. 10.1128/AEM.01284-08
Zhou G, Liu H, He J, Yuan Y, Yuan Z: The occurrence of Bacillus cereus, B. thuringiensis and B. mycoides in Chinese pasteurized full fat milk. Int J Food Microbiol 2008, 121: 195-200. 10.1016/j.ijfoodmicro.2007.11.028
Zhou G, Zheng D, Dou L, Cai Q, Yuan Z: Occurrence of psychrotolerant Bacillus cereus group strains in ice creams. Int J Food Microbiol 2010, 137: 143-146. 10.1016/j.ijfoodmicro.2009.12.005
Ackermann HW, Abedon ST 2000.http://www.mansfield.ohio-state.edu/~sabedon/names/2000/
Ackermann HW, Azizbekyan RR, Emadi Konjin HP, Lecadet MM, Seldin L, Yu MX: New Bacillus bacteriophage species. Arch Virol 1994, 135: 333-344. 10.1007/BF01310018
Ackermann HW: Phage Classification and Characterization. Methods Mol Biol 2009, 501: 127-140. 10.1007/978-1-60327-164-6_13
Allan BJ, Davies P, Carstens EB, Kropinski AM: Characterization of the genome of Pseudomonas aeruginosa bacteriophage φPLS27 with particular reference to the ends of the DNA. J Virol 1989, 63: 1587-1594.
Lavigne R, Darius Pl, Summer E, Seto D, Mahadevan P, Nilsson A, Ackermann H, Kropinski A: Classification of Myoviridae bacteriophages using protein sequence similarity. BMC Microbiol 2009, 9: 224. 10.1186/1471-2180-9-224
Perkus ME, Shub DA: Mapping the genes in the terminal redundancy of bacteriophage SPO1 with restriction endonucleases. J Virol 1985, 56: 40-48.
Klumpp J, Lavigne R, Loessner M, Ackermann HW: The SPO1-related bacteriophages. Arch Virol 2010, 155: 1547-1561. 10.1007/s00705-010-0783-0
Stewart CR, Casjens SR, Cresawn SG, Houtz JM, Smith AL, Ford ME, Peebles CL, Hatfull GF, Hendrix RW, Huang WM, Pedulla ML: The Genome of Bacillus subtilis bacteriophage SPO1. J Mol Biol 2009, 388: 48-70. 10.1016/j.jmb.2009.03.009
Uchiyama J, Rashel M, Maeda Y, Takemura I, Sugihara S, Akechi K, Muraoka A, Wakiguchi H, Matsuzaki S: Isolation and characterization of a novel Enterococcus faecalis bacteriophage φEF24C as a therapeutic candidate. FEMS Microbiol Lett 2008, 278: 200-206. 10.1111/j.1574-6968.2007.00996.x
Uchiyama J, Rashel M, Takemura I, Wakiguchi H, Matsuzaki S: In silico and in vivo evaluation of bacteriophage φEF24C, a candidate for treatment of Enterococcus faecalis infections. Appl Environ Microbiol 2008, 74: 4149-4163. 10.1128/AEM.02371-07
Carlton RM, Noordman WH, Biswas B, de Meester ED, Loessner MJ: Bacteriophage P100 for control of Listeria monocytogenes in foods: Genome sequence, bioinformatic analyses, oral toxicity study, and application. Regul Toxicol Pharmacol 2005, 43: 301-312. 10.1016/j.yrtph.2005.08.005
Klumpp J, Dorscht J, Lurz R, Bielmann R, Wieland M, Zimmer M, Calendar R, Loessner MJ: The terminally redundant, nonpermuted genome of Listeria bacteriophage A511: a model for the SPO1-like myoviruses of Gram-positive bacteria. J Bacteriol 2008, 190: 5753-5765. 10.1128/JB.00461-08
Kwan T, Liu J, DuBow M, Gros P, Pelletier J: The complete genomes and proteomes of 27 Staphylococcus aureus bacteriophages. Proc Natl Acad Sci USA 2005, 102: 5174-5179. 10.1073/pnas.0501140102
Lee J, Shin H, Son B, Ryu S: Complete genome sequence of Bacillus cereus bacteriophage BCP78. J Virol 2012, 86: 637-638. 10.1128/JVI.06520-11
Ruhfel RE, Robillard NJ, Thorne CB: Interspecies transduction of plasmids among Bacillus anthracis, B. cereus, and B. thuringiensis. J Bacteriol 1984, 157: 708-711.
Mccloy EW: Lysogenicity and immunity to Bacillus phage. W J Gen Microbiol 1958, 18: 198-220.
Schuch R, Pelzek AJ, Kan S, Fischetti VA: Prevalence of Bacillus anthracis-like organisms and bacteriophages in the intestinal tract of the earthworm Eisenia fetida. Appl Environ Microbiol 2010, 76: 2286-2294. 10.1128/AEM.02518-09
Kong M, Kim M, Ryu S: Complete genome sequence of Bacillus cereus bacteriophage PBC1. J Virol 2012,86(11):6379-6380. 10.1128/JVI.00706-12
Schuch R, Fischetti VA: The secret life of the anthrax agent Bacillus anthracis: bacteriophage-mediated ecological adaptations. PLoS One 2009, 4: e6532. 10.1371/journal.pone.0006532
Thorne CB: Transduction in Bacillus cereus and Bacillus anthracis. Microbiol Mol Biol Rev 1968, 32: 358-361.
Thorne CB, Holt SC: Cold lability of Bacillus cereus bacteriophage CP-51. J Virol 1974, 14: 1008-1012.
Yelton DB, Thorne CB: Transduction in Bacillus cereus by each of two bacteriophages. J Bacteriol 1970, 102: 573-579.
Ahmed R, Sankar-Mistry P, Jackson S, Ackermann HW, Kasatiya SS: Bacillus cereus phage typing as an epidemiological tool in outbreaks of food poisoning. J Clin Microbiol 1995, 33: 636-640.
Lee W, Billington C, Hudson JA, Heinemann JA: Isolation and characterization of phages infecting Bacillus cereus. Microbiol: Lett. Appl; 2011.
Loessner MJ, Maier SK, Daubek-Puza H, Wendlinger G, Scherer S: Three Bacillus cereus bacteriophage endolysins are unrelated but reveal high homology to cell wall hydrolases from different bacilli. J Bacteriol 1997, 179: 2845-2851.
Abbasifar RAM, Kropinski PM, Sabour HWA, Lingohr EJ, Griffiths MW: Complete Genome Sequence of Cronobacter sakazakii Bacteriophage vB_CsaM_GAP161. J Virol 2012, 86: 13806-13807. 10.1128/JVI.02546-12
Liu J, Mushegian A: Displacements of prohead protease genes in the late operons of double-stranded-DNA bacteriophages. J Bacteriol 2004, 186: 4369-4375. 10.1128/JB.186.13.4369-4375.2004
Ishii T, Yoshida K, Terai G, Fujita Y, Nakai K: DBTBS: a database of Bacillus subtilis promoters and transcription factors. Nucl Acids Res 2001, 29: 278-280. 10.1093/nar/29.1.278
Söding J, Biegert A, Lupas AN: The HHpred interactive server for protein homology detection and structure prediction. Nucl Acids Res 33: W244-W248.
Hanlon GW: Bacteriophages: an appraisal of their role in the treatment of bacterial infections. Int J Antimicrob Agents 2007, 30: 118-128.
Kutter E, Sulakvelidze A: Bacteriophages: biology and applications. Boca Raton, FL: CRC Press; 2005.
Hyde SJ, Eckenroth BE, Smith BA, Eberley WA, Heintz NH, Jackman JE, Doublié S: tRNAHis guanylyltransferase (THG1), a unique 3′-5′ nucleotidyl transferase, shares unexpected structural homology with canonical 5′-3′ DNA polymerases. Proc Natl Acad Sci USA 2010, 107: 20305-20310. 10.1073/pnas.1010436107
Tåquist H, Cui Y, Ardell DH: TFAM 1.0: an online tRNA function classifier. Nucl. Acids Res 2007, 35: W350-W353. 10.1093/nar/gkm393
Sarkar T, Petrov AS, Vitko JR, Santai CT, Harvey SC, Mukerji I, Hud NV: Integration host factor (IHF) dictates the structure of polyamine-DNA condensates: implications for the role of IHF in the compaction of bacterial Chromatin. Biochemistry 2009, 48: 667-675. 10.1021/bi8019965
Schmitz JE, Schuch R, Fischetti VA: Identifying active phage lysins through functional viral metagenomics. Appl Environ Microbiol 2010, 76: 7181-7187. 10.1128/AEM.00732-10
Kropinski AM, Borodovsky M, Carver TJ, Cerdeño-Tárraga AM, Darling A, Lomsadze A, Mahadevan P, Stothard P, Seto D, Domselaar G, Wishart DS: In silico identification of genes in bacteriophage DNA. Methods Mol Biol 2009, 502: 57-89. 10.1007/978-1-60327-565-1_6
Zafar N, Mazumder R, Seto D: CoreGenes: A computational tool for identifying and cataloging “core” genes in a set of small genomes. BMC Bioinforma 2002, 3: 12. 10.1186/1471-2105-3-12
Łobocka M, Hejnowicz MS, Dąbrowski K, Gozdek A, Kosakowski J, Witkowska M, Ulatowska MI, Weber-Dąbrowska B, Kwiatek M, Parasion S, Gawor J, Kosowska H, Głowacka A: Genomics of staphylococcal Twort-like phages–potential therapeutics of the post-antibiotic era. Adv Virus Res 2012, 83: 143-216.
Camp AH, Losick R: A feeding tube model for activation of a cell-specific transcription factor during sporulation in Bacillus subtilis. Genes Dev 2009, 23: 1014-1024. 10.1101/gad.1781709
Sharpe ME, Errington J: Postseptational chromosome partitioning in bacteria. Proc Natl Acad Sci USA 1995, 92: 8630-8634. 10.1073/pnas.92.19.8630
Wu LJ, Errington J: Use of asymmetric cell division and spoIIIE mutants to probe chromosome orientation and organization in Bacillus subtilis. Molec Microbiol 1998, 27: 777-786. 10.1046/j.1365-2958.1998.00724.x
Anany HE: Biocontrol of FoodBorne Bacterial Pathogens Using Immobilized Bacteriophages. University of Guelph Libraries, Ph.D. University of Guelph; 2010.
Abedon S, Hyman P, Thomas C: Experimental examination of bacteriophage latent-period evolution as a response to bacterial availability. Appl Environ Microbiol 2003, 69: 7499-7506. 10.1128/AEM.69.12.7499-7506.2003
Lingohr E, Frost S, Johnson RP: Determination of bacteriophage genome size by pulsed-field gel electrophoresis. Methods Mol Biol 2009, 502: 19-25. 10.1007/978-1-60327-565-1_3
Koski L, Gray M, Lang BF, Burger G: AutoFACT: An automatic functional annotation and classification tool. BMC Bioinforma 2005, 6: 151. 10.1186/1471-2105-6-151
Finn RD, Mistry J, Tate J, Coggill P, Heger A, Pollington JE, Gavin OL, Gunasekaran P, Ceric G, Forslund K, Holm L, Sonnhammer ELL, Eddy SR, Bateman A: The Pfam protein families database. Nucl Acids Res 2010, 38: D211-D222. 10.1093/nar/gkp985
Laslett D, Canback B: ARAGORN, a program to detect tRNA genes and tmRNA genes in nucleotide sequences. Nucl Acids Res 2004, 32: 11-16. 10.1093/nar/gkh152
Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucl Acids Res 1997, 25: 955-964.
Ermolaeva MD, Khalak HG, White O, Smith HO, Salzberg SL: Prediction of transcription terminators in bacterial genomes. J Mol Biol 2000, 301: 27-33. 10.1006/jmbi.2000.3836
Gautheret D, Lambert A: Direct RNA motif definition and identification from multiple sequence alignments using secondary structure profiles. J Mol Biol 2001, 313: 1003-1011. 10.1006/jmbi.2001.5102
Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucl Acids Res 2003, 31: 3406-3415. 10.1093/nar/gkg595
Juretić D, Zoranić L, Zucić D: Basic charge clusters and predictions of membrane protein topology. J Chem Inf Comput Sci 2002, 42: 620-632. 10.1021/ci010263s
Käll L, Krogh A, Sonnhammer ELL: Advantages of combined transmembrane topology and signal peptide prediction—the Phobius web server. Nucl Acids 2007, 35: 429-432. 10.1093/nar/gkm256
Krogh A, Larsson B, von Heijne G, Sonnhammer ELL: Predicting transmembrane protein topology with a hidden markov model: application to complete genomes. J Mol Biol 2001, 305: 567-580. 10.1006/jmbi.2000.4315
Darling AE, Mau B, Perna NT: progressiveMauve: multiple genome alignment with gene gain, loss and rearrangement. PLoS One 2010, 5: e11147. 10.1371/journal.pone.0011147
Stothard P, Wishart DS: Circular genome visualization and exploration using CGView. Bioinformatics 2005, 21: 537-539. 10.1093/bioinformatics/bti054
During this research A.M.K. was supported by a Discovery Grant from the Natural Sciences and Engineering Research Council of Canada. The funding body did not have any role in the design of the experiments, in the collection, analysis, and interpretation of data; in the writing of the manuscript; or in the decision to submit the manuscript for publication.
The authors have no competing interests to disclose.
The bulk of the research and writing was conducted by TFEA, under the supervision of MWG. YMS carried out the proteomic studies. The initial annotation was performed by AV, and then extended by AMK and TFEA. EJL carried out the PFGE analyses. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Table S1: General features of putative ORFs of Bc431v3 and homology to proteins in the database. In addition, protein motifs including transmembrane domains are included. (DOC 590 KB)
Additional file 3: Table S3: Potential promoters and rho-independent terminators in vB_BceM_Bc431v3 with WebLogo of the consensus. (DOC 218 KB)
Additional file 5: Table S5: Comparison of codon usage in Bc431v3 genes with that of Bacillus cereus Q1. Those codons which are 50% more common in the phage are coloured orange. (DOC 132 KB)
About this article
Cite this article
El-Arabi, T.F., Griffiths, M.W., She, YM. et al. Genome sequence and analysis of a broad-host range lytic bacteriophage that infects the Bacillus cereus group. Virol J 10, 48 (2013). https://doi.org/10.1186/1743-422X-10-48
- Bacillus cereus
- Host range