Next Article in Journal
Isoconazole and Clemizole Hydrochloride Partially Reverse the Xeroderma Pigmentosum C Phenotype
Next Article in Special Issue
Hypervirulent Klebsiella pneumoniae Sequence Type 420 with a Chromosomally Inserted Virulence Plasmid
Previous Article in Journal
Microbiota Signals during the Neonatal Period Forge Life-Long Immune Responses
Previous Article in Special Issue
Mammalian Neuropeptides as Modulators of Microbial Infections: Their Dual Role in Defense versus Virulence and Pathogenesis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bacillus pumilus 15.1, a Strain Active against Ceratitis capitata, Contains a Novel Phage and a Phage-Related Particle with Bacteriocin Activity

1
Institute of Biotechnology, Faculty of Sciences, University of Granada, 18071 Granada, Spain
2
Department of Biochemistry and Molecular Biology I, Faculty of Sciences, University of Granada, 18071 Granada, Spain
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(15), 8164; https://doi.org/10.3390/ijms22158164
Submission received: 18 June 2021 / Revised: 19 July 2021 / Accepted: 26 July 2021 / Published: 29 July 2021
(This article belongs to the Collection Microbial Virulence Factors)

Abstract

:
A 98.1 Kb genomic region from B. pumilus 15.1, a strain isolated as an entomopathogen toward C. capitata, the Mediterranean fruit fly, has been characterised in search of potential virulence factors. The 98.1 Kb region shows a high number of phage-related protein-coding ORFs. Two regions with different phylogenetic origins, one with 28.7 Kb in size, highly conserved in Bacillus strains, and one with 60.2 Kb in size, scarcely found in Bacillus genomes are differentiated. The content of each region is thoroughly characterised using comparative studies. This study demonstrates that these two regions are responsible for the production, after mitomycin induction, of a phage-like particle that packages DNA from the host bacterium and a novel phage for B. pumilus, respectively. Both the phage-like particles and the novel phage are observed and characterised by TEM, and some of their structural proteins are identified by protein fingerprinting. In addition, it is found that the phage-like particle shows bacteriocin activity toward other B. pumilus strains. The effect of the phage-like particles and the phage in the toxicity of the strain toward C. capitata is also evaluated.

1. Introduction

Bacillus pumilus 15.1, a natural bacterium isolated from a partially decomposed common reed plant as a strain active against larvae of the Mediterranean fruit fly (Ceratitis capitata) [1], has been extensively characterised in our research group for the last 10 years. B. pumilus 15.1 is a motile bacterium, shows the high capability of metabolising a wide range of compounds, and can form biofilms [2]. The strain contains a 7.8 Kb cryptic plasmid, with 33 copies per chromosome [3], and at least one uncharacterised megaplasmid. The strain genome has been previously described [4] and is publicly available (GenBank accession number LBDK00000000.1). B. pumilus 15.1 produces, during the sporulation phase, a parasporal crystal [2] resembling the one produced by Bacillus thuringiensis, a well-known entomopathogen. The crystal is of proteinaceous nature, is solubilised at low temperature and is mainly composed of the enzyme oxalate decarboxylase. The enzyme can catalyse the transformation of oxalate to formate, a previously recognised neurotoxic compound. It has been proposed that oxalate decarboxylase could be a potential virulence factor [5], although the mechanism of toxicity is still under study.
Many phage and phage-related structures can be associated with virulence factors in bacteria, such as the cholera toxin from Corynebacterium diphtheriae [6], the Shiga toxin from the enterohaemorrhagic Escherichia coli [7] or the botulism toxin from Clostridium botulinum [8,9].
Bacteriophages or phages are viruses that use bacteria as a host, either replicating, causing the death of the host (lytic cycle), or inserting their genome into the bacterial chromosome (prophage) and linking the viral multiplication to the bacterial growth cycle (lysogenic cycle). The ‘decision’ between lysis and lysogeny depends on the phage nature (not all phages are capable of lysogeny), host physiological state, genetic compatibility and host and phage abundance, among other factors [10].
Bacteriophages are extremely common in nature. Early studies showed that most of the bacteria (60–70%) whose genomes have been sequenced have at least one prophage in their genomes [11]. However, recent studies show that this percentage is around 30% [12]. Lysogens can reach up to 90% in some species, such as Pectobacterium spp. or Dickeya spp. [13]. In the case of Bacillus species, such as B. thuringiensis, all of the strains investigated for prophage prediction contained putative prophage regions, and 81% of the strains (50 out of 61) have complete prophage regions with an average of 2.21 prophages/genome [14].
Bacteriophages display high levels of mosaicism in their genomes. This level of mosaicism is particularly relevant in Lambdoid phages. It is explained, in part, by the large number of horizontal gene transfer events that take place thanks to the action of recombinases, often encoded by phages themselves [15]. Prophages and the so-called cryptic (or defective) prophages represent the largest reservoir of functional genes for temperate phages [11]. Virulent and defective prophages can exchange genes recovering the activity in the case of defective prophages or modifying the activity of the virulent ones [16].
Some prophages also contain what is known as lysogenic conversion genes [17], independent transcriptional units with non-essential functions for phage replication, but which can change bacterial fitness, enhancing the survival under certain conditions of the lysogenic bacteria over the non-lysogenic one. These lysogenic conversion genes can range from antibiotic resistance genes [18], to the production of toxins [9], or genes that provide increased adhesion properties to the host bacteria [19], among others.
It is widely accepted that once phage DNA is inserted in a bacterial genome, the prophage might be subjected to mutations or partial deletions that could lead to ‘phage domestication’ by losing genes or through gene inactivation, becoming the so-called cryptic or defective phages [20]. This process of domestication had led to the evolution of multiple molecular systems, which are morphologically similar to phage structures, and provide significant advantages to the bacteria that bear them. Among such systems, Gene Transfer Agents (GTAs) [21], Phage-Inducible Chromosomal Islands (PICIs) [22], bacterial type VI secretion systems (T6SS) [23], some bacteriocins, such as pyocins [24], the Serratia anti-feeding prophage (AFP) [25] or the Photorhabdus virulence cassette (PVC) [26] can be mentioned.
GTAs are phage-like particles with the ability to pack small fragments of host DNA instead of the phage-like particle genome and delivering them into a recipient bacterial cell. The DNA fragments packaged in GTAs particles are smaller than their coding region, and therefore, are incapable of self-replication by horizontal transfer. GTA production is highly controlled by host systems, including quorum-sensing [27] and is restricted to ~3% of the population [28], so the vast majority of the population can be a receptor of the DNA encapsulated in these particles.
PICIs are the Phage Inducible Chromosomal Island family of mobile elements. The best characterised are the Staphylococcus aureus Pathogenic Islands (SAPIs), a 15 Kb mobile element that encodes virulence factors and uses helper phages for transmission [29].
Bacteriocins are ribosomally-synthesised peptides produced by bacteria with a narrow bactericide activity against related strains. Among them, the phage-tail-like Pseudomonas aeruginosa pyocins are a kind of bacteriocins similar to bacteriophage structures. R-type pyocins are syringe-like contractile particles, similar to phage tails from the Myoviridae family, which insert themselves into the bacterial membrane and kill the bacteria by forming pores in the membrane and disrupting the membrane potential [30]. In addition, F-type pyocins also resemble phage tails, but in this case to the non-contractile form of siphoviridae tails [24].
Although all of the previously described systems are very different in their function (injecting DNA from bacteria to bacteria in the case of GTAs, injecting virulence factors from bacteria to eukaryotic cells in the case of PVCs and AFPs, injecting lytic enzymes from bacteria to bacteria in the case of pyocins, or secreting toxic effectors in the case of T6SS) all have the common characteristic that they are evolutionarily-related phage structures designed to efficiently transport macromolecules from bacteria to other cells.
The PBSX particle, a phage-like particle produced by Bacillus subtilis, is another phage-related molecular system that shows some interesting features [31]. This phage-like particle is a mitomycin-inducible structure (dependent on the SOS system) that shows bactericidal activity against B. subtilis related bacterial strains, similar to bacteriocins [32]. The PBSX particle is encoded by an approximately 33 Kb region identified in the B. subtilis 168 chromosome (Accession number NC_000964.3) [33]. The most remarkable feature of PBSX particles is that they package random 13 Kb DNA fragments from the host cell [34] inside the head when induced with mitomycin.
The interesting functions that phages and phage-related systems represent for bacteria is an ongoing research area in Microbiology and Molecular Biology, especially in the case of pathogenic bacteria, as a high number of the remarkable features that pathogens show are related to them.
With the aim of searching for potential virulence factors in B. pumilus 15.1 strain, bioinformatics analysis was performed, and a region of 113 kb in the genome of the bacteria (Contig 59), containing a high number of phage proteins coding ORFs was found. Contig 59 contained a 98.1 Kb in silico-predicted prophage present in its sequence. We found that part of the 98.1 Kb region showed a genomic organisation and structure similar to B. subtilis PBSX particles. The PBSX-like particle, called Bp15.1PLP, was observed under TEM in mitomycin-induced cultures and was concentrated by polyethylene glycol precipitation. The Bp15.1PLP packages host bacterial DNA in its capsid, with an average size of 9 Kb, instead of the 13 Kb DNA fragments that B. subtilis PBSX packages. In addition, B. pumilus PBSX-like particles show bactericidal activity against all the B. pumilus strains assayed. The major structural proteins of the particle were analysed by SDS-PAGE and identified by MALDI-TOF, and the possible effect of the Bp15.1PLP particle on bacterial toxicity on C. capitata larvae was investigated. Apart from Bp15.1PLP, a novel B. pumilus phage, named Bp15.1Hope, was also found in the mitomycin-induced B. pumilus 15.1 lysates. It was proven that Bp15.1Hope genome is also contained in the studied 98.1 Kb region (Contig 59).

2. Results and Discussion

2.1. Prophage Prediction in the B. pumilus 15.1 Genome

During characterisation of the strain B. pumilus 15.1 [2], we observed that under some culture conditions (such as elevated temperatures (45–50 °C) or oxygen depletion), a phage-like lysis process occurred, as the optical density of the cultures suddenly dropped after a normal initial growth phase (unpublished work). This observation reminded us of the behaviour of the induction of the lytic cycle of a temperate phage and led us to hypothesise the presence of a prophage in the genome of B. pumilus 15.1.
The draft genome sequence of B. pumilus 15.1, previously published by our group [4] and available under the accession number LBDK00000000.1, was analysed for the detection of prophage regions using PHASTER [35]. Only one region in the bacterial genome, with a potential phage organisation, was identified in contig 59 (NZ_LBDK01000059).
The in silico analysis of the phage-predicted region, with a length of 98.1 Kb (Figure 1), showed that it contained 127 predicted open reading frames (ORFs), with an average of 44.4% in GC content, a GC content higher than the average reported by B. pumilus strains (41.6% GC). That was a striking fact, as previous studies have established that bacteriophage usually has lower GC content than their hosts [36,37]. One tRNA for Alanine was also identified with the tRNA-scanSE algorithm with a score of 98.8.
A preliminary comparative analysis of the sequence identified by PHASTER at nucleotide level showed two distinct regions, with different prevalence levels in databases. A region of 28.695 Kb, widely spread in the genome of bacteria belonging to the B. pumilus group, was detected at the 5’ end of the 98.1 Kb region (Figure 1; shadowed in green), while another region of 60.213 Kb, at the 3´ end of the 98.1 Kb region (shadowed in purple), contained nucleotide sequences rarely found in sequenced genomes from bacteria belonging to the B. pumilus group (only five matches in the NCBI database; specified below). This region was flanked by a repeat sequence of 45 bp that comprised part of the tRNA-Ala identified by tRNAScan-se.
The two identified regions showed a significantly different GC content (Figure 1); while the 28.7 Kb region showed a GC content (40.5%) similar to the host (41.6%), the 60.2 Kb region had a GC content of 47.7%. The significant difference in the GC% content, together with the fact that this region is flanked by short, direct repeats, led us to hypothesise that this DNA region had a heterologous origin and is present in the genome of B. pumilus 15.1, due to insertion or transposition events. Based on these results, we decided to study these two regions separately.

2.2. Genome Organisation of the 28.7 Kb Region

All the ORFs identified in the 28.7 Kb region were analysed with BLASTp to predict conserved protein domains and to determine the potential function of the putative proteins. The results of this analysis are detailed in Table 1. The predicted function for each ORF was determined with BLASTp unless otherwise specified. DNA and amino acid sequences of the identified ORFs are included in Table S1 of Supplementary Materials. The 28.7 Kb region contained 40 predicted ORFs ranging from 141 bp (ORF 34) to 4272 bp (ORF 21). All of the ORFs, except for one (ORF 1), were encoded in the direct strand.
The 28.7 Kb region showed a similar genome organisation to PBSX, a well-characterised defective prophage from Bacillus subtilis 168 [33,34] that produce phage-like particles. Even the different cassettes present in PBSX defective phage (early, middle and late cassettes) could be easily identified in the 28.7 Kb region from B. pumilus 15.1 (PBSX-like region from now on) (Table 1). We predicted that the early cassette could be constituted by ORF 1, a backward-transcribed ORF, with a predicted Helix-Turn-Helix (HTH) conserved domain, similar to the Xenobiotic Response Element (XRE) transcriptional regulators (42.86% identity), the PBSX repressor responsible for maintaining the lysogenic state of the phage-like particle by repressing the transcription of the structural genes [38]. In addition, we predicted that the region from ORF 2 to ORF 6 could constitute the middle cassette, as putative proteins encoded by this region are related to DNA transcription. Although ORF 2, ORF 3 and ORF 4 showed no homology in the PBSX particle, ORF 6 showed a sigma-like factor conserved domain and 49% identity with Xpf, the RNA polymerase sigma factor of the PBSX particle. This protein is known to regulate the expression of the late cassette, and therefore, the expression of the structural genes of the PBSX particle [39]. ORF 2 showed an XRE family transcriptional regulator domain, and ORF 4 was homologous to proteins identified as Holliday junction resolvases, enzymes that both resolve the Holliday structures resulting from the recombination process, and are responsible for (i) debranching DNA structures in phages prior to the viral genome packaging into head particles and for (ii) degrading host DNA [40].
The late cassette, comprising all the structural genes that form the phage particle and genes responsible for lysis, was identified between ORF 7 and ORF 37. This region included the capsid protein (ORF 11), the tail sheath protein (ORF 17), the tail tube protein (ORF 18) and the baseplate J-like protein (ORF 26), among several others. In addition, a lysis module, with 3 ORFs (ORF 38-ORF 40) and with the same organisation as PBSX defective phage, was predicted at the 3´ end of the late cassette. ORF 38 and ORF 39 were identified as putative holins, proteins that accumulate and form aggregates to hydrolyse the host plasmatic membrane, and facilitate access to the endolysin (ORF 40) for the peptidoglycan degradation at the lysis stage [41,42].

2.3. The PBSX-like Region Is Highly Conserved in Bacteria Belonging to the B. pumilus Group

A search for homologous sequences to the PBSX-like region in NCBI database using Megablast, identified 38 complete genomes belonging to the B. pumilus group (Table 2) and two genomes from strains classified as Bacillus sp. These sequences were compared for Average Nucleotide Identity (ANI) with the JSpeciesWS algorithm, and their nucleotide synteny was determined with a dotplot (Figure 2).
The selected sequences were highly similar and showed extremely conserved synteny (Figure 2 and Figure 3), although, as expected, the level of synteny decreased in the more distant strains. Nucleotide similarity ranged between 87.08 to 96.15%, with an average similarity of 90.76%.
When the PBSX-like region from B. pumilus 15.1 was compared to B. subtilis defective phages PBSX and PBSZ, using Discontiguous Megablast algorithm, identity and cover sequence were 68.37% and 41% for PBSX and 68.72% and 41% for PBSZ. Using JSpeciesWS, the ANI observed was 66.13% for PBSX and 65.08% for PBSZ.
As shown in Figure 2, the synteny of all genes in PBSX-like regions, compared by dotplot, was highly conserved, except for two regions (marked with arrows); one at the middle (black arrow) and one at the 3´-end of the 28.7 Kb region (light grey arrow). Both variable regions were close to genes related to cell lysis. The first region (inside of the B. pumilus 15.1 ORF 21) was similar to the xkdO gene from PBSX, encoding a putative tape measure protein (TMP), a protein present in the tail tube, in fully functional and defective phages, that determines the tail length. TMP is known to be ejected after the irreversible binding of the phage to the host receptors and before the DNA injection in T4 phage (E. coli phage of Myoviridae family) [44]. TMP contains a conserved transglycosylase domain involved in peptidoglycan degradation. Some TMPs proteins have been proposed to participate in phage infection by hydrolysing the beta 1–4 glycosidic bond in the peptidoglycan, allowing the phage particle to inject the DNA into the host [45]. The variable region was found in the middle of ORF21, right outside of the transglycosylase domain, which is highly conserved and normally localised at the C-terminal end of the predicted protein. We could speculate that ORF 21 encodes a protein with transglycosylase activity at the C-terminal end and a variable domain to recognise specific molecules. In addition, the next ORF (ORF 22) encodes a protein with a LysM conserved domain and a peptidoglycan-binding domain, which is a homologue of XkdP (72% similarity) in PBSX.
The second variable region ranged from ORF 32 to ORF 37, a region that encodes putative proteins of unknown function but, given their position in the genome they could be related to proteins in the tail fibres [46]. Two of them, the ORF 31 (63 aa) and ORF 34 (46 aa), showed an Xkdx conserved domain of unknown function. These ORFs are normally found near holins and endolysins, so it is tempting to speculate that they could be part of the lysis module. In fact, ORF 38 encodes a holin, which surprisingly has more identity with Bhla holin from B. subtilis prophage SPBeta (100% cover and 77% similarity) than with the PBSX holin Xhla (40% cover and 34% similarity), an indication of the genomic plasticity that phages and phage-related particles have.
The ORF 32 to ORF 37 region did not show sequence homology in PBSX and constitute the most variable region in all the prophage genomes, with modifications and gene reorganisations that are not observed in other regions. This observation is in agreement with the idea of the region encoding tail fibres, as they are responsible for the recognition of the bacterial receptors that determine the specificity of the particle [47] (Figure 3).

2.4. Genome Organisation of the 60.2 kb Region

A comparative analysis of the 60.2 Kb region was also performed, to detect conserved protein domains and determine the potential function of the putative proteins. The result of this analysis is shown in Table 3. The predicted function for each ORF was determined by BLASTp, unless otherwise specified. The 60.2 Kb region contained 86 ORFs ranging from 141 bp (ORF 77) to 2805 bp (ORF 71). Most of the ORFs were encoded on the direct strand except for 10 ORFs (ORF 12, 13, 42, 47, 48, 56, 82, 83, 84 and 86). Six of these 10 reverse-transcribed genes are related to transcriptional regulation (ORF 12, 13, 82 and 84) or DNA integration/recombination (ORF 56 and ORF 86).
A close analysis of the 60.2 Kb region shows that the predicted ORFs can be grouped into two subregions according to their putative function. This separation does not imply that they have different origin, and it was only made for description purposes. The 5´end of the 60.2 Kb region (ORF 1 to 48; Table 3 shadowed in light grey; named subregion 1) of approximately 31.5 Kb, contained around 20 ORFs, with known functions related to DNA replication, DNA regulation and DNA metabolism; specifically, a DNA polymerase I (ORF 21), a helicase (ORF 10), a DNA primase (ORF 11), a Holiday junction resolvase (ORF 25), and several genes involved in nucleotide metabolism (ORF 31 to 33, ORF 36, ORF 38) and transcriptional regulation (ORF 2, ORF 12, ORF 13 and ORF 46).
The 3´ end subregion (ORF 49 to ORF 86; Table 3, not shadowed; named subregion 2), with a size of almost 28.7 Kb, contained phage-related proteins and ORFs encoding putative proteins with integrase/recombinase conserved domains (ORF 56 and ORF 86, both of which are backward-transcribed). Subregion 2 contains two ORFs encoding large and small subunits of the terminase (ORF 57 and ORF 58), a portal protein (ORF 59), capsid proteins (ORF 60 and 62), tail proteins (ORF 66, ORF 71 to 73), and a baseplate protein (ORF 75). In addition, putative proteins related to the host lysis can be observed, such as an N-acetylmuramoyl-L-alanine amidase (lysin) (ORF 79) and a holin (ORF 80). Finally, five ORFs related to DNA metabolism/regulation, four of which are transcribed backwards (ORF 82, ORF 83, ORF 84 and ORF 86), can be found at the 3′ terminal end of subregion 2. ORF 82 showed a conserved YoID domain with unknown function, but related to the UmuD subunit of Polimerase V in Gram-negative bacteria. ORF 84 and 85 (transcribed backward and forward, respectively) are two putative transcriptional regulators showing HTH conserved domains, and ORF 86 showed a conserved domain in proteins related to site-specific recombinases/integrases (XerD superfamily).
Both subregions 1 and 2 showed a similar GC content (47.84% for subregion 1 and 47.54% for subregion 2), so no differences in the phylogenetic origin are predicted. Determining the origin of the 60.2 Kb fragment could be a very useful piece of information, so it is planned to be addressed in the near future.
The 60.2 Kb region, found just downstream of the PBSX-like region, comprise 86 ORFs with no similar genome organisation to any other known phage, although PHASTER predicted this region to be part of a prophage, given that several phage-related putative proteins are encoded by this region, especially at the 3´end (subregion 2). The 60.2 Kb region was less common in bacterial genomes compared to the PBSX-like region, and only five strong coincidences were found in NCBI database, all present in genomes of strains belonging to the B. pumilus group. These strains were B. pumilus PDSLzg-1 (Accession number CP016784), B. cellulasensis GLB197 (Accession number CP018574), B. pumilus NCTC10337 (Accession number LT906438), B. altitudinis P-10 (Accession number CP024204) and B. aerophilus strain 232 (Accession number CP026008.1). We compared the 60.2 Kb region with the five homolog sequences and found they were highly similar (82–89% identity) with an ANI of 84.17–87.6%. In addition, the 60.2 Kb region showed a conserved synteny in all the strains with no variable regions (Figure 4 and Figure 5).
The analysis of the 60.2 Kb region in the five Bacillus genomes showed that the region seems to be a DNA insertion, as it was predicted for B. pumilus 15.1, as (i) the CG content of the region is different from the host and (ii) this region is flanked by a direct repeat sequence that is part of a tRNA gene.
Comparative studies of all the strains containing the 60.2 Kb insert showed two different insertions sites: (i) The tRNA gene for Alanine, with a 45 bp direct repeat sequence (in B. pumilus 15.1 and B. pumilus NCTC10337) and (ii) the tRNA gene for Arginine, with a 48 bp direct repeat sequence in the other four strains (B. altitudinis P10, B. aerophilus 232, B. cellulasensis GLB197 and B. pumilus PDSLzg-1). The insertions, however, seem to reconstitute the tRNA sequence, so the function of these genes is presumably maintained in the six strains analysed.
The relative position between the 60.2 Kb region and the PBSX-like region in the six genomes was also studied (Figure 6). B. pumilus 15.1 and B. pumilus NCTC10337 showed the same organisation and insertion site, in the tRNA-Ala gene, next to PBSX-like region. Both of them shared the same short 45 bp direct repeat, which includes an important part of the tRNA sequence (data not shown); however, we observed differences downstream, where the region showed no homology between these two strains.
In B. aerophilus 232, B. cellulasensis GLB197, B. altitudinis P-10 and B. pumilus PDSLzg-1 the insertion was observed in a tRNA-Arg located ~300 Kb upstream the PBSX-like region. All of the strains showed the same insertion site and genomic organisation of the flanking areas.

2.5. Phylogenetic Analysis of the 98.1 Kb Region by Whole-Genome and Single-Gene Comparison

The phylogenetic relationships between bacteriophages are currently analysed by whole genome comparison [48] and single gene analysis. Different phage genes have been used for this last purpose; the major capsid protein [49], the integrase [50], or the large subunit of the terminase [51], among others. To determine the relationship between the prophage regions identified by PHASTER and previously known phages or phage-related structures, these two approaches were used.
For the whole-genome approach, all of the complete sequenced phages using B. pumilus as host and present in the Virus-host database (https://www.genome.jp/virushostdb/view/; last access 1 December 2020) were compared for nucleotide similarity by Average Nucleotide Identity and dotplot matrix comparison with the 28.7 and 60.2 Kb regions. No sequence similarity was found between them (data not shown). For that reason, a different phylogenetic study was conducted using the large subunit of the terminase, a protein that is commonly used for phage classification. Depending on the type of terminase, phages can be classified by their DNA packaging strategy, and the type of phage DNA ends [52] can be predicted. An analysis of the 98.1 Kb region identified by PHASTER rendered two sequences encoding a large subunit of terminases (ORF 8 in the PBSX-like region and ORF 57 in the 60.2 Kb region). Both sequences were compared to the 109 large terminase subunit sequences found in databases, some of which are from phages with well-known DNA packaging mechanisms (marked by * in Figure 7). As shown, most of the terminases from phages seemed to group together in four large clades, although not all of them were grouped (Listeria phage B054, Paenibacillus phage Emery, Bacillus phage SPBeta, Bacillus phage 0305 phi8-36, Mycobacterium phage Nigel, Bacillus phage phBC6A52, or Campylobacter phage Cpt1). Terminases from the seven well-defined DNA packaging strategies [52,53] grouped together in eight clades.
The headful packaging strategy was present in two clusters, the T4 like packaging system (Yellow), and the P22-like headful packaging system (Blue). Surprisingly, terminases from Streptococcus phage O1205 and sfi11 did not group in the main headful packaging P22-like clade. We observed differences in phylogenetic relationship between the two terminases found in B. pumilus 15.1 (marked with red square boxes). While the terminase found in the PBSX-like region (ORF8) clearly clustered into the P22-like HeadFul packaging group, where it formed a defined cluster with PBSX-like terminases, the terminase found in the 60.2 Kb region (tRNA-Ala insertion) (ORF 57) did not convincingly cluster in any group of well-characterised terminases and remained unclassified. These results suggested that both terminases are responsible for different DNA packaging strategies and support the hypothesis that they have an independent origin.

2.6. Isolation and Morphological Characterisation of Phage-Like Particles from B. pumilus 15.1

To demonstrate the functional capacity of the PBSX-like ORFs found in the genome of B. pumilus 15.1 to form phage-like particles, we performed an induction assay with mitomycin C. Seven hours after the addition of mitomycin C, an important decrease in the optical density of the B. pumilus 15.1 culture (final OD ≤ 0.2) was observed. Following induction, the lysate was examined by TEM, and phage-like particles with a typical Myoviridae morphology were observed (Figure 8A). The particles had very distinctive heads (39.53 ± 2.85 nm wide (n ≥ 3)), and a tail (227.34 nm ± 3.99 nm in length and 22.66 ± 1.81 nm in wide (n ≥ 3)) that seemed to be contractile, as two different forms of the particles, an extended form (white arrows) and a contracted form (black arrows), were observed. The contracted forms showed axial striations, whereas cross-striations were observed in the extended forms. The tails seemed to be comprised of an internal structure (11.53 ± 1.01 nm wide, black arrow head) and an external sheath (24.19 ± 2.76 nm wide, white arrow head) (Figure 8A,B). In addition, fibres were observed at the distal part of the tail (empty arrows). Some of the particles did not show the head, and only the tail was observed (empty arrow heads). The observed Phage-Like Particles obtained from a mitomycin C-induced B. pumilus 15.1 culture were called Bp15.1PLPs.
As Bp15.1PLP particles were extremely similar to the PBSX defective phage [32], it was tempting to speculate that those particles were encoded by the 28.7 Kb region found by PHASTER.

2.7. DNA Content in Bp15.1PLP Particles

After obtaining PEG-precipitated Bp15.1PLP particles from a mitomycin C-induced culture, the DNA present inside their capsids was isolated and visualised in an agarose gel. One single band around 9 Kb was observed, a size very different than the predicted for the PBSX-like region by in silico studies (28.695 Kb). The digestion of the 9 Kb DNA fragment with five different restriction enzymes (BamHI, HindIII, SalI, EcoRI and BglII) (Figure 9A) rendered no defined bands in most of the cases; instead, a DNA smear was obtained, which was especially patent in HindIII and EcoRI digestions, and seems to indicate the heterogeneity of the extracted DNA.
To identify the packed DNA in Bp15.1PLP, the EcoRI and BamHI digested DNA was introduced in EcoRI and BamHI digested pUC19 vectors, ligated and used to transform E. coli DH5α cells. Plasmid content from nine white colonies was sequenced using M13 F (5′-GTTTTCCCAGTCACGAC 3′) and M13 R (5′-CAGGAAACAGCTATGAC 3) primers from the vector.
Sequencing results showed that the cloned DNA belonged to different regions of the bacterial genome of B. pumilus15.1, randomly packaged into the Bp15.1PLP particle heads, without any relation to the prophage genome identified by PHASTER and with no homology with phage-like particle proteins. To rule out the possibility of bacterial DNA contamination in the purified particles, the experiment was repeated three times, with increasing concentration of DNAse after bacterial lysis and previous to Bp15.1PLP PEG-precipitation. The results were the same, so we concluded that the particle did not pack its encoding DNA, but random DNA from the host. In the bibliography, we further found that this phenomenon was also described in the defective phage PBSX [34].

2.8. Protein Characterisation of Bp15.1PLPs

To identify the major structural proteins present in Bp15.1PLP particle, a sucrose-gradient purified Bp15.1PLP suspension was subjected to 12% SDS-PAGE electrophoresis and stained with Coomassie brilliant blue. The Bp15.1PLP particle contained at least 15 proteins ranging from 110 to 17 KDa (Figure 9B). The most abundant bands were excised (shown with arrows) and identified by MALDI-TOF (Table 4). The results showed that all the Bp15.1PLP proteins identified by MALDI-TOF were encoded by the 98.1 Kb region under study. Proteins bands 1, 3 and 4 were encoded by ORFs present in the PBSX-like region, specifically ORF 17 (Phage tail sheath protein), ORF 11 (Phage major capsid protein) and ORF 18 (Phage tail tube protein). However, protein in band 2 was encoded by the ORF 62, identified as a major capsid protein by PHASTER (and as a Hypothetical Protein by Blastp) at the 60.2 Kb region (subregion 2), 41 Kb downstream the end of the PBSX-like region.
The surprising fact that Bp15.1PLP could be formed with proteins encoded from distant genes in the chromosome is similar to that which happens in other phage-like elements, the GTAs. The components of the most studied GTAs, RcGTA from Rhodobacter capsulatus and VSH-1 from Brachyspira hyodysenteriae, have a multilocus genome organisation and are encoded in five and two different loci, respectively, separated by hundreds of kilobases, and some of them with different evolutionary origins [54]. However, although this level of mosaicism is frequent in genomes, a simpler explanation could not be ruled out. Given that the PBSX-like region is highly distributed among the Bacillus genus does not make much sense that its production is dependent on the expression of a gene present in a DNA fragment scarcely found in bacterial genomes. Hence, the hypothesis of B. pumilus 15.1 having a second phage or phage-like particle, encoded by the 60.2 Kb region and also induced by mitomycin, was made. This hypothesis was supported by the fact that two different, non-phylogenetically related terminases were present in the 98.1 Kb fragment, one in the PBSX-like region and the other in the 60.2 Kb region.

2.9. B. pumilus 15.1 Produces a Phage That Is Co-Purified with Bp15.1PLP Particles

To test the hypothesis of a second phage-related particle produced by B. pumilus 15.1, TEM micrographs obtained from sucrose purified samples were closely inspected again, in search for less abundant phages or phages-like particles that could have been overlooked in the first analysis. The result of this close inspection was the identification of another phage structure, much less abundant than the Bp15.1PLP particle, which could be in accordance with our hypothesis (Figure 10). The relative abundance of this phage is also in agreement with the amount of protein found in SDS-PAGE analysis (band 2), identified as a major capsid protein by PHASTER, but being the least abundant of the selected proteins analysed by finger printing. The newly discovered particle showed a typical phage structure, with a 73.41 ± 0.15 (n = 3) nm hexagonal shaped head, highly electron-dense in the centre and clear at the edges, and a tail 155.03 ± 4.6 nm long and 12.01 ± 0.4 nm wide. The tail was slightly curved, non-contractile, showed cross-striation and seamed to end in a wider area (Figure 10C). Phage morphology was consistent with bacteriophages from the Siphoviridae family of the Caudovirales order (small heads and long flexible and non-contractile tails).
Given the fact that one of the proteins obtained in the sucrose purified lysates is produced by ORF 62 at the 60.2 kb region and that this region also contains two ORFs encoding large and small subunits of a terminase (ORF 57 and ORF 58), a portal protein (ORF 59), other capsid protein apart from the one codified by ORF 62 (ORF 60), tail proteins (ORF 66, ORF 71 to 73), a baseplate protein (ORF 75), putative proteins related to host lysis, such as an N-acetylmuramoyl-L-alanine amidase (lysin) (ORF 79) and a holin (ORF 80), make us draw the conclusion that this region could be responsible for the synthesis of a phage, as many of the required elements for phage assembly and function are present in this region.
After a bibliographic search, no publications reporting the production of phages in any of the five strains containing the 60.2 Kb homolog region in B. pumilus 15.1 (B. pumilus PDSLzg-1, B. cellulasensis GLB197, B. pumilus NCTC10337, B. altitudinis P-10 and B. aerophilus strain 232) were found (May 2021). To search for similar phages that the one observed in B. pumilus 15.1 lysates, an extensive search was conducted using the Viral Genome Database at NCBI (https://www.ncbi.nlm.nih.gov/genome/viruses/; last access 1 May 2021). This database compiles all the known and currently available viral genomes. In this database, there are 143 completely sequenced phage genomes (plus 2 incomplete) that have been either found in a Bacillus strain or which can use a Bacillus strain as a host. The Bacillus species with the highest number of genome-sequenced phages were B. thuringiensis (44 phages), B. cereus (31 phages), B. subtilis (21 phages), B. megaterium (17 phages), B. pumilus (14 phages) and B. anthracis (3 phages) (Table 5). Other Bacillus species were hosts for only one phage, and 7 phages used unclassified Bacillus species (Bacillus sp.) as a host. The size of Bacillus phages ranged from 18,379 nt (phage BsuP-Goel) from B. subtilis to 497,513 nt (phage G) from B. megaterium (data not shown). Phages found in B. pumilus (Table 6) showed a genome size ranging from 18,466 nt (phage Bpu_PumA1) to 160,281 nt (phage Bobb from B. pumilus SAFR32). The size of the 60.2 Kb inserted fragment found in B. pumilus 15.1 lies into this size range.
A whole sequence comparison of the 60.2 Kb region found in B. pumilus 15.1 with the 14 phage genome sequences found in B. pumilus was performed using tBlastx and plotted by Easyfig. No relevant similarity (nor genomic synteny) of the 60.2 Kb region with any of the previously described B. pumilus phages was found after the analysis. This result supports the hypothesis that the 60.2 Kb region encodes for a novel B. pumilus phage, named Bp15.1Hope. The characterisation of the novel phage and the search for a suitable host has been already undertaken in our lab.

2.10. Bacteriocin-Like Activity of PEG-Precipitated Phage like Particles

It has been reported that the B. subtilis PBSX particle has a bacteriocin-like activity against host-related strains, which seems to play an important role for the strain in competing with other bacteria in the same environment [32]. To determine whether the Bp15.1PLP particle has similar bacteriocin-like activity, fifteen strains from our laboratory collection were tested as described in the Materials and Methods section. Several strains were susceptible to Bp15.1PLP particles (Table 7), showing a clear halo when a drop of a suspension of Bp15.1PLP (and Bp15.1Hope) was placed over the strain in a spot test. All of the susceptible strains were B. pumilus strains, including B. pumilus 15.1C, an extra-chromosomal DNA cured strain obtained from B. pumilus 15.1 and previously described by our group [5]. Bp15.1Hope was not able to form plaques in any of the strains tested in lawn assays (data not shown), so the halos observed seem to indicate that Bp15.1PLP has a strong bacteriocin-like activity. This activity must be considered very narrow, as only strains from the same species were susceptible. This specificity agrees with the previously described diversity in the tail fibre sequences found in Bp15.1PLP, as fibres are responsible for recognising bacterial receptors.
Moreover, it was observed that Bp15.1PLP was not active against B. pumilus 15.1 strain, the strain that produces the particles. This is a common feature among phages, as lysogenic hosts are not susceptible to the lytic action of the phage that they produce, and apparently, Bp15.1PLP has the same behaviour. Surprisingly, the closely related B. pumilus 15.1C strain, previously obtained by treatment with acridine orange for plasmid curing, is susceptible to the action of Bp15.1PLP particles. This phenomenon will be further investigated in the near future as it is not straight-forward to understand how curing plasmids in a strain can make them susceptible to the action of their own producing PLP.

2.11. Bioassaying Ceratitis Capitata Larvae with Phage like Particles

Phage-related systems, such as the Photorabdus virulence cassette (PVC) and the anti-feeding prophage (AFP) from Serratia, have been previously described to be the delivery system for virulence factors in entomopathogenic bacteria [25,26]. To study the contribution of Bp15.1PLP and Bp15.1Hope in the toxicity of B. pumilus 15.1 toward C. capitata, diet contamination bioassays using first-instar larvae were performed. A Bp15.1PLP/Bp15.1Hope preparation consisting of a mitomycin-induced lysate culture was bioassayed in combination or not with a B. pumilus 15.1 sporulated culture (bioassayed in sub-lethal conditions). The rationale behind this experiment was that if phage particles are somehow involved in C. capitata toxicity, the increase in the number of particles and the induction of the genes related to the particles would have a synergistic effect on the toxicity of a B. pumilus 15.1 sporulated culture. Particle preparation showed no toxicity against C. capitata (even when it was 200 times concentrated) (Table 8), and no synergistic effect between the sporulated culture and the phage particles was observed. These results seem to suggest that neither Bp15.1PLP nor Bp15.1Hope are related to B. pumilus 15.1 toxicity against C. capitata.

2.12. Concluding Remarks

The 98.1 kb DNA region, found in Contig 59 of the entomopathogenic strain B. pumilus 15.1, encodes a PBSX-like particle, called Bp15.1PLP, and an undescribed B. pumilus phage named Bp15.1Hope. The Bp15.1PLP particle is encoded by a 28.7 kb region, and Bp15.1Hope seems to be a genomic insertion in a tRNA gene. The Bp15.1PLP particle packages bacterial host DNA in its head, and shows bacteriocin activity toward all of the B. pumilus strains tested, but not to other related bacteria (such as B. subtilis). The Bp15.1Hope phage is not similar to any of the currently described B. pumilus phages, and it can be classified as a siphovirus according to its morphology. Both particles are expressed in liquid cultures when the strain B. pumilus 15.1 is induced by mitomycin, but at a very different rate, as Bp15.1PLP is much more abundant than Bp15.1Hope. None of these particles is related to the entomopathogenic activity that B. pumilus 15.1 shows toward C. capitata. The separation of these two particles and DNA sequencing of Bp15.1Hope have become our next objectives to definitely prove that Bp15.1Hope particles are produced by the 60.2 kb region.

3. Materials and Methods

3.1. Bacterial Strains and Growth Conditions

The following bacterial strains were used in this study: Bacillus pumilus 15.1 [1], B. pumilus15.1_C [2] (extra-chromosomal-cured B. pumilus 15.1 strain), B. pumilus M1, B. pumilus M2, (both kindly provided by Dr C. Calvo) [55], B. pumilus 8A3 (ATCC 7061), (kindly provided by Dr C. Berry, Cardiff University), B. thuringiensis 78/11, Bacillus subtilis 168C, B. cereus 169, Micrococcus sp., Arthrobacter sp., Lactococcus garvieae, Lactococcus NZ9000 and Listeria innocua (CECT 4030) (kindly provided by Dr Martinez Bueno, Granada University).
Luria-Bertani (LB) medium broth or agar was routinely used for growing bacteria at 30 °C. Brain-Heart Infusion (BHI) broth was used for growing Lactococcus, and Listeria strains at 30 °C. When sporulation was required, T3 medium [56] was used, and the culture was incubated at 30 °C for 72h.

3.2. Induction, Concentration and Purification of Phage-Like Particles

Phage-like particles were directly obtained from B. pumilus 15.1 strain by adding mitomycin C to the bacterial culture. Fresh LB broth was inoculated with an overnight culture (1:100 dilution) and incubated at 30 °C until the culture reached an optical density at 600 nm (OD600) of 0.4. Then, mitomycin C (0.5 µg/mL final concentration) was added, and the culture was incubated until lysis occurred (OD600 <0.2). Two percent (v/v) chloroform was added to the lysate and incubated for 20 min under agitation.
For the rapid concentration of phage-like particles, when no extra purification was necessary, the lysate was centrifuged at 4000× g for 15 min to remove cell debris. The supernatant was then filtered through a 0.2 µm nitrocellulose filter and treated with 10 µg/mL DNAse and RNAse for 30 min. Finally, the supernatant was centrifuged at 100,000× g for 90 min to pellet phage-like particles. The pellet was resuspended in 1/200 of the initial volume of the culture in SM medium (NaCl 0.1 M; MgSO4 0.1 M; Tris 0.05 M; Gelatine 0.01%; pH 7.5).
Phage-like particles were also concentrated and partially purified by a polyethylene glycol (PEG) precipitation method according to [57] with minor modifications. Briefly, lysates were centrifuged twice at 4000× g for 30 min to remove cell debris. Lysates were treated with 1 µg/mL DNAse and RNAse at 37 °C for 30 min, and then, solid NaCl was added to a final concentration of 1 M. Lysates were incubated for 1 h on ice and centrifuged at 12,000× g for 30 min. To pellet phage-like particles, 10% w/v PEG 8000 was added to the supernatant, incubated on ice for at least 1 h, and then centrifuged at 12,000× g for 30 min. Pellets were resuspended in 1/400 of the initial culture volume in SM medium, vortexed with the same volume of chloroform and centrifuged for 10 min at 16,000× g. The upper aqueous phases, containing the phage-like particles, were collected and stored at 4 °C until use.
When necessary, particle suspension was placed on a preformed discontinuous sucrose gradient (60-67-75 and 83%) and ultra-centrifuged (53,000× g, 4 °C) in a Beckman JSW24.38 rotor for 16 h. Fractions were collected from the top of the gradient by pipetting, washed three times in SM buffer and analysed by SDS-PAGE and Transmission Electronic Microscopy (TEM) to observe proteins and phage-like particles, respectively.

3.3. Transmission Electron Microscopy (TEM)

A B. pumilus 15.1 culture was induced with mitomycin C and centrifuged (4000× g 15 min 4 °C) when lysis was observed. The supernatant was filtered through a 0.2 µm nitrocellulose filter and sent to the Microscopy service to be fixed and negatively stained with 1% uranyl acetate. Samples were examined under a ZEISS EM 902 transmission electron microscope (CIC, Universidad de Granada, Granada, Spain).
High-titre phage-like samples obtained by sucrose gradient, ultracentrifugation (100,000× g for 1 h) or PEG 8000 precipitation were diluted and also analysed by TEM to check the structural integrity of phage-like particles.

3.4. Bacteriocin-like Activity

The bacteriocin-like activity of phage-like particles was carried out in a spot test, with some modifications [58]. Overnight cultures (12–16 h) of the test bacteria were streaked on an LB plate with the help of a sterile toothpick; once dried, a drop of phage-like particle suspension, after PEG-precipitation, was placed on one end of the streak, incubated O/N at 37 °C and observed for clear zones on the following day.

3.5. Bp15.1PLP DNA Isolation

DNA extraction from sucrose-purified phage-like particles was performed using the DNeasy Blood and Tissue DNA kit (Qiagen) according to the manufacturer’s instructions from a maximum of 200 µL of a phage-like particle suspension.

3.6. Phage-Like Particle DNA Mapping and Sequencing

The DNA contained in the phage-like particles was isolated as described above and digested with restriction enzymes (NEB) according to the supplier’s instructions. Digestions were analysed in a 0.8% w/v agarose gel. HindIII and EcoRI fragments were ligated to HindIII or EcoRI-digested pUC19 vectors and introduced in E. coli DH5α competent cells by heat shock [59].
The resulting transformants were screened by toothpick assay [60]. Positive clones were isolated, plasmids extracted and sequenced with M13 Forward and Reverse primers at the Genomic Service of Instituto de Parasitología y Biología Molecular Lopez Neyra, Granada, Spain.

3.7. Phage-Like Particle Analysis by SDS-PAGE and Protein Identification

Ten microlitres of a suspension of sucrose-purified phage-like particles were analysed in a 12% wt/v SDS-PAGE gel and stained with Coomassie brilliant blue according to standard procedures. PAGE Ruler Plus Pre-Stained Protein Ladder (Thermo Scientific) was used as size protein markers.
Protein bands were excised and digested with trypsin before identification by MALDI-TOF. The data obtained were analysed by MASCOT at the Centro de Biología Molecular Severo Ochoa, Madrid, Spain.

3.8. Ceratitis Capitata Larval Bioassay

C. capitata bioassays were performed as described previously by our group [1] with some modifications. A B. pumilus 15.1 sporulated culture was 40-fold concentrated by lyophilisation, and the lysate obtained after mitomycin C induction was used directly or after PEG concentration. To determine whether the phage-like particles had some synergistic effect on the toxicity of the bacteria, both phage-like particles and bacteria were bioassayed together. Bacterial or phage-like samples were mixed with the artificial C. capitata larva diet in a diet contamination bioassay. Deionised water and SM buffer were used as a negative control in the bioassays. All bioassays were performed at least twice from separate cultures under the conditions described previously by our group [1].Mortality was recorded 14 days after the beginning of the bioassay. Results were statistically analysed using a one-way ANOVA method.

3.9. Bioinformatics

B. pumilus 15.1 draft genome sequence (Accession number LBDK00000000.1) was analysed with PHASTER algorithm [35] (PHAge Search Tool Enhanced Release http://phaster.ca, University of Alberta, last access date 1 September 2019) for potential prophage detection. Predicted phage sequences were analysed by comparative analysis with BLASTp [61] against the total NCBI protein database. tRNA identification was performed using tRNAscan-SE [62]. DotPlot comparison was performed with Gepard 1.41 [43], using a word size of 10, and the average nucleotide identity based on Blast+ (ANIb) analysis was performed with JspeciesWS software [63]. The neighbour-joining tree was constructed using ClustalW [64] and MEGAX [65], following the Poisson model with a bootstrapping of 1000 trials.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ijms22158164/s1.

Author Contributions

A.F.-F. carried out all the experiments and wrote the first draft of the manuscript. S.V. conceived the idea, obtained the funding, designed the experiments and supported in some experimental work. A.O. helped supervise the project and the interpretation of the results. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by the Andalucian Operative Program, grant number B-BIO-081-UGR18. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data presented in this study are insert in the article or in the Supplementary Materials.

Acknowledgments

We thank María José Guerrero, from the Scientific Instrumentation Centre at the University of Granada for her technical support with TEM. We thank the Economy and Knowledge Council from the Andalusia Government for funding this project.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Molina, C.; Caña-Roca, J.F.; Osuna, A.; Vilchez, S. Selection of a Bacillus pumilus strain highly active against Ceratitis capitata (wiedemann) larvae. Appl. Environ. Microbiol. 2010, 76, 1320–1327. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Garcia-Ramon, D.C.; Molina, C.A.; Osuna, A.; Vílchez, S. An in-depth characterization of the entomopathogenic strain Bacillus pumilus 15.1 reveals that it produces inclusion bodies similar to the parasporal crystals of Bacillus thuringiensis. Appl. Microbiol. Biotechnol. 2016, 100, 3637–3654. [Google Scholar] [CrossRef]
  3. Garcia-Ramon, D.C.; Luque-Navas, M.J.; Molina, C.A.; del Val, C.; Osuna, A.; Vilchez, S. Identification, sequencing and comparative analysis of pBp15.S plasmid from the newly described entomopathogen Bacillus pumilus 15.1. Plasmid 2015, 82, 17–27. [Google Scholar] [CrossRef] [PubMed]
  4. García-Ramón, D.C.; Palma, L.; Osuna, A.; Berry, C.; Vílchez, S. Draft genome of the entomopathogenic bacterium Bacillus pumilus 15.1, a strain highly toxic towards Ceratitis capitata, the Mediterranean fruit fly. Genome Annoucement 2015, 3, 1–2. [Google Scholar] [CrossRef]
  5. Garcia-Ramon, D.C.; Berry, C.; Tse, C.; Fernández-Fernández, A.; Osuna, A.; Vílchez, S. The parasporal crystals of Bacillus pumilus strain 15.1: A potential virulence factor? Microb. Biotechnol. 2018, 11, 302–316. [Google Scholar] [CrossRef]
  6. Waldor, M.K.; Mekalanos, J.J. Lysogenic conversion by a filamentous phage encoding cholera toxin. Science 1996, 272, 1910–1913. [Google Scholar] [CrossRef] [Green Version]
  7. O’Brien, A.D.; Newland, J.W.; Miller, S.F.; Holmes, R.K.; Smith, H.W.; Formal, S.B. Shiga-like toxin-converting phages from Escherichia coli strains that cause hemorrhagic colitis or infantile diarrhea. Science 1984, 226, 694–696. [Google Scholar] [CrossRef]
  8. Freeman, V.J. Studies on the virulence of bacteriophage-infected strains of Corynebacterium diphtheriae. J. Bacteriol. 1951, 61, 675–688. [Google Scholar] [CrossRef] [Green Version]
  9. Eklund, M.W.; Poysky, F.T.; Reed, S.M.; Smith, C.A. Bacteriophage and the toxigenicity of Clostridium botulinum type C. Science 1971, 172, 480–482. [Google Scholar] [CrossRef]
  10. Howard-Varona, C.; Hargreaves, K.R.; Abedon, S.T.; Sullivan, M.B. Lysogeny in nature: Mechanisms, impact and ecology of temperate phages. ISME J. 2017, 11, 1511–1520. [Google Scholar] [CrossRef] [Green Version]
  11. Casjens, S. Prophages and bacterial genomics: What have we learned so far? Mol. Microbiol. 2003, 49, 277–300. [Google Scholar] [CrossRef] [PubMed]
  12. Roux, S.; Hallam, S.J.; Woyke, T.; Sullivan, M.B. Viral dark matter and virus–host interactions resolved from publicly available microbial genomes: For each bacteria or archaea phylum (or phylum-level group), corresponding viruses in RefSeq (gray) and VirSorter curated data set (red) are indicated with c. Elife 2015, 4, e08490. [Google Scholar] [CrossRef]
  13. Czajkowski, R. May the phage be with you? Prophage-like elements in the genomes of soft rot Pectobacteriaceae: Pectobacterium spp. and Dickeya spp. Front. Microbiol. 2019, 10, 138. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Fu, Y.; Wu, Y.; Yuan, Y.; Gao, M. Prevalence and diversity analysis of candidate prophages to provide an understanding on their roles in Bacillus thuringiensis. Viruses 2019, 11, 388. [Google Scholar] [CrossRef] [Green Version]
  15. Bobay, L.M.; Touchon, M.; Rocha, E.P.C. Manipulating or superseding host recombination functions: A dilemma that shapes phage evolvability. PLoS Genet. 2013, 9, e1003825. [Google Scholar] [CrossRef] [Green Version]
  16. De Paepe, M.; Hutinet, G.; Son, O.; Amarir-Bouhram, J.; Schbath, S.; Petit, M.A. Temperate phages acquire DNA from defective prophages by relaxed homologous recombination: The role of Rad52-like recombinases. PLoS Genet. 2014, 10, e1004181. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Brüssow, H.; Canchaya, C.; Hardt, W.; Bru, H. Phages and the evolution of bacterial pathogens: From genomic rearrangements to lysogenic conversion. Microbiol. Mol. Biol. Rev. 2004, 68, 560–602. [Google Scholar] [CrossRef] [Green Version]
  18. Colomer-Lluch, M.; Jofre, J.; Muniesa, M. Antibiotic resistance genes in the bacteriophage DNA fraction of environmental samples. PLoS ONE 2011, 6, e17549. [Google Scholar] [CrossRef] [Green Version]
  19. Robinson, C.M.; Sinclair, J.F.; Smith, M.J.; O’Brien, A.D. Shiga toxin of enterohemorrhagic Escherichia coli type 0157:H7 promotes intestinal colonization. Proc. Natl. Acad. Sci. USA 2006, 103, 9667–9672. [Google Scholar] [CrossRef] [Green Version]
  20. Canchaya, C.; Fournous, G.; Brüssow, H. The impact of prophages on bacterial chromosomes. Mol. Microbiol. 2004, 53, 9–18. [Google Scholar] [CrossRef]
  21. Lang, A.S.; Beatty, J.T. Genetic analysis of a bacterial genetic exchange element: The gene transfer agent of Rhodobacter capsulatus. Proc. Natl. Acad. Sci. USA 2000, 97, 859–864. [Google Scholar] [CrossRef] [Green Version]
  22. Penadés, J.R.; Christie, G.E. The Phage-inducible chromosomal islands: A family of highly evolved molecular parasites. Annu. Rev. Virol. 2015, 2, 181–201. [Google Scholar] [CrossRef]
  23. Leiman, P.G.; Basler, M.; Ramagopal, U.A.; Bonanno, J.B.; Sauder, J.M.; Pukatzki, S.; Burley, S.K.; Almo, S.C.; Mekalanos, J.J. Type VI secretion apparatus and phage tail-associated protein complexes share a common evolutionary origin. Proc. Natl. Acad. Sci. USA 2009, 106, 4154–4159. [Google Scholar] [CrossRef] [Green Version]
  24. Michel-Briand, Y.; Baysse, C. The pyocins of Pseudomonas aeruginosa. Biochimie 2002, 84, 499–510. [Google Scholar] [CrossRef]
  25. Hurst, M.R.H.; Beard, S.S.; Jackson, T.A.; Jones, S.M. Isolation and characterization of the Serratia entomophila antifeeding prophage. FEMS Microbiol. Lett. 2007, 270, 42–48. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Yang, G.; Dowling, A.J.; Gerike, U.; Ffrench-Constant, R.H.; Waterfield, N.R. Photorhabdus virulence cassettes confer injectable insecticidal activity against the wax moth. J. Bacteriol. 2006, 188, 2254–2261. [Google Scholar] [CrossRef] [Green Version]
  27. Westbye, A.B.; Leung, M.M.; Florizone, S.M.; Taylor, T.A.; Johnson, J.A.; Fogg, P.C.; Beatty, J.T. Phosphate concentration and the putative sensor kinase protein CckA modulate cell lysis and release of the Rhodobacter capsulatus gene transfer agent. J. Bacteriol. 2013, 195, 5025–5040. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Hynes, A.P.; Mercer, R.G.; Watton, D.E.; Buckley, C.B.; Lang, A.S. DNA packaging bias and differential expression of gene transfer agent genes within a population during production and release of the Rhodobacter capsulatus gene transfer agent, RcGTA. Mol. Microbiol. 2012, 85, 314–325. [Google Scholar] [CrossRef] [PubMed]
  29. Penadés, J.R.; Chen, J.; Quiles-Puchalt, N.; Carpena, N.; Novick, R.P. Bacteriophage-mediated spread of bacterial virulence genes. Curr. Opin. Microbiol. 2015, 23, 171–178. [Google Scholar] [CrossRef] [PubMed]
  30. Uratani, Y.; Hoshino, T. Pyocin R1 inhibits active transport in Pseudomonas aeruginosa and depolarizes membrane potential. J. Bacteriol. 1984, 157, 632–636. [Google Scholar] [CrossRef] [Green Version]
  31. Seaman, E.; Tarmy, E.; Marmur, J. Inducible phages of Bacillus subtilis. Biochemistry 1964, 3, 607–613. [Google Scholar] [CrossRef]
  32. Okamoto, K.; Mudd, J.A.; Mangan, J.; Huang, W.M.; Subbaiah, T.V.; Marmur, J. Properties of the defective phage of Bacillus subtilis. J. Mol. Biol. 1968, 34, 413–428. [Google Scholar] [CrossRef]
  33. Wood, H.E.; Dawson, M.T.; Devine, K.M.; McConnell, D.J. Characterization of PBSX, a defective prophage of Bacillus subtilis. J. Bacteriol. 1990, 172, 2667–2674. [Google Scholar] [CrossRef] [Green Version]
  34. Andersont, L.M.; Bott, K.F. DNA packaging by the Bacillus subtilis defective bacteriophage PBSX. J. Virol. 1985, 54, 773–780. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Arndt, D.; Grant, J.R.; Marcu, A.; Sajed, T.; Pon, A.; Liang, Y.; Wishart, D.S. PHASTER: A better, faster version of the PHAST phage search tool. Nucleic Acids Res. 2016, 44, W16–W21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Rocha, E.P.C.; Danchin, A. Base composition bias might result from competition for metabolic resources. Trends Genet. 2002, 18, 291–294. [Google Scholar] [CrossRef]
  37. Almpanis, A.; Swain, M.; Gatherer, D.; McEwan, N. Correlation between bacterial G+C content, genome size and the G+C content of associated plasmids and bacteriophages. Microb. Genomics 2018, 4, 168. [Google Scholar] [CrossRef] [PubMed]
  38. Wood, H.E.; Devine, K.M.; McConnell, D.J. Characterisation of a repressor gene (xre) and a temperature-sensitive allele from the Bacillus subtilis prophage, PBSX. Gene 1990, 96, 83–88. [Google Scholar] [CrossRef]
  39. Mcdonnell, G.E.; Wood, H.; Devine, K.M.; Mcconnell, D.J. Genetic control of bacterial suicide: Regulation of the induction of PBSX in Bacillus subtilis. J. Bacteriol. 1994, 176, 5820–5830. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Wyatt, H.D.M.; West, S.C. Holliday junction resolvases. Cold Spring Harb. Perspect. Biol. 2014, 6, a023192. [Google Scholar] [CrossRef] [Green Version]
  41. Krogh, S.; Jørgensen, S.T.; Devine, K.M. Lysis genes of the Bacillus subtilis defective prophage PBSX. J. Bacteriol. 1998, 180, 2110–2117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Longchamp, P.F.; Mauel, C.; Karamata, D. Lytic enzymes associated with defective prophages of Bacillus subtilis: Sequencing and characterization of the region comprising the N-acetylmuramoyl-L-alanine amidase gene of prophage PBSX. Microbiology 1994, 140, 1855–1867. [Google Scholar] [CrossRef] [Green Version]
  43. Krumsiek, J.; Arnold, R.; Rattei, T. Gepard: A rapid and sensitive tool for creating dotplots on genome scale. Bioinformatics 2007, 23, 1026–1028. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Hu, B.; Margolin, W.; Molineux, I.J.; Liu, J. Structural remodeling of bacteriophage T4 and host membranes during infection initiation. Proc. Natl. Acad. Sci. USA 2015, 112, E4919–E4928. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Rodríguez-Rubio, L.; Gutiérrez, D.; Martínez, B.; Rodríguez, A.; Götz, F.; García, P. The tape measure protein of the Staphylococcus aureus bacteriophage vB_SauS-phiIPLA35 has an active muramidase domain. Appl. Environ. Microbiol. 2012, 78, 6369–6371. [Google Scholar] [CrossRef] [Green Version]
  46. Li, X.; Koç, C.; Kühner, P.; Stierhof, Y.D.; Krismer, B.; Enright, M.C.; Penadés, J.R.; Wolz, C.; Stehle, T.; Cambillau, C.; et al. An essential role for the baseplate protein Gp45 in phage adsorption to Staphylococcus aureus. Sci. Rep. 2016, 6, 26455. [Google Scholar] [CrossRef] [Green Version]
  47. Dowah, A.S.A.; Clokie, M.R.J. Review of the nature, diversity and structure of bacteriophage receptor binding proteins that target Gram-positive bacteria. Biophys. Rev. 2018, 10, 535–542. [Google Scholar] [CrossRef] [Green Version]
  48. Pope, W.H.; Bowman, C.A.; Russell, D.A.; Jacobs-Sera, D.; Asai, D.J.; Cresawn, S.G.; Jacobs, W.R.; Hendrix, R.W.; Lawrence, J.G.; Hatfull, G.F. Whole genome comparison of a large collection of mycobacteriophages reveals a continuum of phage genetic diversity. Elife 2015, 4, e06416. [Google Scholar] [CrossRef]
  49. Zhan, Y.; Huang, S.; Voget, S.; Simon, M.; Chen, F. A novel roseobacter phage possesses features of podoviruses, siphoviruses, prophages and gene transfer agents. Sci. Rep. 2016, 6, 4–11. [Google Scholar] [CrossRef]
  50. Bobay, L.M.; Rocha, E.P.C.; Touchon, M. The adaptation of temperate bacteriophages to their host genomes. Mol. Biol. Evol. 2013, 30, 737–751. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Merrill, B.D.; Grose, J.H.; Breakwell, D.P.; Burnett, S.H. Characterization of Paenibacillus larvae bacteriophages and their genomic relationships to firmicute bacteriophages. BMC Genomics 2014, 15, 745. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Casjens, S.R.; Gilcrease, E.B.; Winn-Stapley, D.A.; Schicklmaier, P.; Schmieger, H.; Pedulla, M.L.; Ford, M.E.; Houtz, J.M.; Hatfull, G.F.; Hendrix, R.W. The generalized transducing Salmonella bacteriophage ES18: Complete genome sequence and DNA packaging strategy. J. Bacteriol. 2005, 187, 1091–1104. [Google Scholar] [CrossRef] [Green Version]
  53. Grose, J.H.; Belnap, D.M.; Jensen, J.D.; Mathis, A.D.; Prince, J.T.; Merrill, B.; Burnett, S.H.; Breakwell, D.P. The genomes, proteomes and structure of three novel phages that infect the Bacillus cereus group and carry putative virulence factors. J. Virol. 2014, 88, 11846–11860. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Hynes, A.P.; Shakya, M.; Mercer, R.G.; Grüll, M.P.; Bown, L.; Davidson, F.; Steffen, E.; Matchem, H.; Peach, M.E.; Berger, T.; et al. Functional and evolutionary characterization of a gene transfer agent’s multilocus “genome”. Mol. Biol. Evol. 2016, 33, 2530–2543. [Google Scholar] [CrossRef] [Green Version]
  55. Uad, I.; Toledo, F.L.; Pozo, C.; Silva, G.A.; Gonzalez-Lopez, J.; Calvo, C. Caracterización de bacterias productoras de bioemulgentes aisladas de ambientes marinos. In Proceedings of the XXI Congreso Nacional de Microbiología Sociedad Española de Microbiología, Seville, Spain, 1 February 2021; p. 292. [Google Scholar]
  56. Travers, R.S.; Martin, P.A.; Reichelderfer, C.F. Selective process for efficient isolation of soil Bacillus spp. Appl. Environ. Microbiol. 1987, 53, 1263–1266. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Yamamoto, K.R.; Alberts, B.M.; Benzinger, R.; Lawhorne, L.; Treiber, G. Rapid bacteriophage sedimentation in the presence of polyethylene glycol and its application to large-scale virus purification. Virology 1970, 40, 734–744. [Google Scholar] [CrossRef]
  58. Bradley, D.E. Ultrastructure of bacteriophage and bacteriocins. Bacteriol. Rev. 1967, 31, 230–314. [Google Scholar] [CrossRef]
  59. Mandel, M.; Higa, A. Calcium-dependent bacteriophage DNA infection. J. Mol. Biol. 1970, 53, 159–162. [Google Scholar] [CrossRef]
  60. Barnes, W. Plasmid detection and sizing in single colony lysates. Science 1977, 195, 393–394. [Google Scholar] [CrossRef]
  61. Altschul, S.F.; Wootton, J.C.; Gertz, E.M.; Agarwala, R.; Morgulis, A.; Schaffer, A.A.; Yu, Y.-K. Protein database searches using compositionally adjusted substitution matrices. FEBS J. 2005, 272, 5101–5109. [Google Scholar] [CrossRef]
  62. Lowe, T.M.; Eddy, S.R. tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25, 955–964. [Google Scholar] [CrossRef] [PubMed]
  63. Richter, M.; Rosselló-Móra, R.; Oliver Glöckner, F.; Peplies, J. JSpeciesWS: A web server for prokaryotic species circumscription based on pairwise genome comparison. Bioinformatics 2016, 32, 929–931. [Google Scholar] [CrossRef] [PubMed]
  64. Madeira, F.; Park, Y.M.; Lee, J.; Buso, N.; Gur, T.; Madhusoodanan, N.; Basutkar, P.; Tivey, A.R.N.; Potter, S.C.; Finn, R.D.; et al. The EMBL-EBI search and sequence analysis tools APIs in 2019. Nucleis Acid Res. 2019, 47, W636–W641. [Google Scholar] [CrossRef] [Green Version]
  65. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Schematic representation of the genomic area detected by PHASTER (coloured) and flanking regions (black). The PBSX-like particle is highlighted in green, while the tRNA-Ala is marked in purple. The subregions found in the tRNA-Ala insertion are also indicated. The GC content is showed at the top of the figure. The blue line indicates the mean GC content of the 98 Kb fragment, and the yellow line shows the average GC content for B. pumilus. The potential function of each ORF was inferred by Blastp.
Figure 1. Schematic representation of the genomic area detected by PHASTER (coloured) and flanking regions (black). The PBSX-like particle is highlighted in green, while the tRNA-Ala is marked in purple. The subregions found in the tRNA-Ala insertion are also indicated. The GC content is showed at the top of the figure. The blue line indicates the mean GC content of the 98 Kb fragment, and the yellow line shows the average GC content for B. pumilus. The potential function of each ORF was inferred by Blastp.
Ijms 22 08164 g001
Figure 2. Whole-genome dot-plot comparison of some PBSX-like particles sequences from bacteria belonging to the B. pumilus group (shown in bold in Table 2). A comparison was made with Gepard 1.40 [43] with a word length of 10 nucleotides.
Figure 2. Whole-genome dot-plot comparison of some PBSX-like particles sequences from bacteria belonging to the B. pumilus group (shown in bold in Table 2). A comparison was made with Gepard 1.40 [43] with a word length of 10 nucleotides.
Ijms 22 08164 g002
Figure 3. Comparative analysis of the PBSX-like regions among different B. pumilus related species. The potential function of each ORF was inferred by Blastp and shown as coloured arrows. tRNAs are indicated by asterisks. Homology studies were made by tBlastx and plotted with Easyfig with a minimal length of 50 bp.
Figure 3. Comparative analysis of the PBSX-like regions among different B. pumilus related species. The potential function of each ORF was inferred by Blastp and shown as coloured arrows. tRNAs are indicated by asterisks. Homology studies were made by tBlastx and plotted with Easyfig with a minimal length of 50 bp.
Ijms 22 08164 g003
Figure 4. Comparative analysis of the 60.2 Kb region present in B. pumilus related species by dotplot. In B. pumilus 15.1 and B. pumilus NCTC10337 the 60.2 Kb was found inserted in a tRNA-Ala gene, while in B. cellulasensis GLB197, B. altitudinis P-10, B. aerophilus 232 and B. pumilus PDSLzg-1 the insertion was observed in a tRNA-Arg gene.
Figure 4. Comparative analysis of the 60.2 Kb region present in B. pumilus related species by dotplot. In B. pumilus 15.1 and B. pumilus NCTC10337 the 60.2 Kb was found inserted in a tRNA-Ala gene, while in B. cellulasensis GLB197, B. altitudinis P-10, B. aerophilus 232 and B. pumilus PDSLzg-1 the insertion was observed in a tRNA-Arg gene.
Ijms 22 08164 g004
Figure 5. Comparative analysis of the 60.2 Kb insertion fragment among different B. pumilus related species. The potential function of each ORF was inferred by Blastp and shown as coloured arrows. tRNAs position are indicated by asterisks. Homology studies were made by tBlastx and plotted with Easyfig with a minimal length of 50 bp.
Figure 5. Comparative analysis of the 60.2 Kb insertion fragment among different B. pumilus related species. The potential function of each ORF was inferred by Blastp and shown as coloured arrows. tRNAs position are indicated by asterisks. Homology studies were made by tBlastx and plotted with Easyfig with a minimal length of 50 bp.
Ijms 22 08164 g005
Figure 6. Genomic organisation comparison of the 60.2 Kb insertion in B. pumilus related strains. The relative position of the tRNA inserts and PBSX-like particles are shown in the six strains and compared to the model strain B. pumilus SAFR 032 (panel 1), which lacks tRNA inserts. tRNA-Arg gene is shown in blue, and tRNA-Ala is shown in pink. In most Bacillus species (panel 1), the PBSX-like region is adjacent to the tRNA-Ala. In B. pumilus 15.1 and NCTC 10337 (panel 2), tRNA insert is adjacent to PBSX-like region, in the tRNA-Ala gene. However, in B. cellulasensis GLB197, B. pumilus PDSLzg-1, B. altitudinis P-10 and B. aerophilus 232, inserts are present in the tRNA-Arg gene, circa 300 bp upstream PBSX-like genome (panel 3).
Figure 6. Genomic organisation comparison of the 60.2 Kb insertion in B. pumilus related strains. The relative position of the tRNA inserts and PBSX-like particles are shown in the six strains and compared to the model strain B. pumilus SAFR 032 (panel 1), which lacks tRNA inserts. tRNA-Arg gene is shown in blue, and tRNA-Ala is shown in pink. In most Bacillus species (panel 1), the PBSX-like region is adjacent to the tRNA-Ala. In B. pumilus 15.1 and NCTC 10337 (panel 2), tRNA insert is adjacent to PBSX-like region, in the tRNA-Ala gene. However, in B. cellulasensis GLB197, B. pumilus PDSLzg-1, B. altitudinis P-10 and B. aerophilus 232, inserts are present in the tRNA-Arg gene, circa 300 bp upstream PBSX-like genome (panel 3).
Ijms 22 08164 g006
Figure 7. Neighbor-joining tree of the large terminase subunit in different phage and phage related structures. The amino acid sequences of large subunit terminases from 99 phages were aligned with ClustalW, and a neighbor-joining tree was prepared with MEGA X, with bootstrapping set to 1000. The phylogenetic relationship was inferred by the Poisson model. The numbers near bifurcations are bootstrap percentages for 1000 trials, and all the branches with a bootstrapping level of less than 30% were collapsed. Phages whose DNA termini was reported are marked with *, and clusters of similar packaging methods are highlighted in different colours marking previously clustered terminases in [51,52], unclustered terminases are not highlighted. B. pumilus phages are marked by black-vertical lines, and the terminases from B. pumilus 15.1 sequence are highlighted by red boxes. Phages are named according to the EBI virus database and virus hostdb.
Figure 7. Neighbor-joining tree of the large terminase subunit in different phage and phage related structures. The amino acid sequences of large subunit terminases from 99 phages were aligned with ClustalW, and a neighbor-joining tree was prepared with MEGA X, with bootstrapping set to 1000. The phylogenetic relationship was inferred by the Poisson model. The numbers near bifurcations are bootstrap percentages for 1000 trials, and all the branches with a bootstrapping level of less than 30% were collapsed. Phages whose DNA termini was reported are marked with *, and clusters of similar packaging methods are highlighted in different colours marking previously clustered terminases in [51,52], unclustered terminases are not highlighted. B. pumilus phages are marked by black-vertical lines, and the terminases from B. pumilus 15.1 sequence are highlighted by red boxes. Phages are named according to the EBI virus database and virus hostdb.
Ijms 22 08164 g007
Figure 8. (A) Transmission Electron Microscopy (TEM) of partially purified Bp15.1PLPs by sucrose gradient. Two different forms were observed, an extended particle (white arrows) and a contracted form (black arrows). In the contracted form, a tail tube (black arrowhead), a tail sheath (white arrowhead) and fibres at the end of the tails (empty arrows) were distinctly observed. Some of the particles lacked a head (empty arrowhead). (B) Schematic representation of the observed phage-like particles with the average size of each part (obtained from at least three measurements).
Figure 8. (A) Transmission Electron Microscopy (TEM) of partially purified Bp15.1PLPs by sucrose gradient. Two different forms were observed, an extended particle (white arrows) and a contracted form (black arrows). In the contracted form, a tail tube (black arrowhead), a tail sheath (white arrowhead) and fibres at the end of the tails (empty arrows) were distinctly observed. Some of the particles lacked a head (empty arrowhead). (B) Schematic representation of the observed phage-like particles with the average size of each part (obtained from at least three measurements).
Ijms 22 08164 g008
Figure 9. (A) Agarose gel electrophoresis of the DNA contained in Bp15.1PLP before and after restriction enzyme digestion. The nucleic acids were extracted from partially purified Bp15.1PLPs after mitomycin C induction and DNAse treatment. The DNA was digested with BamHI B), HinDIII (H), SalI (S), EcoRI (E) and BglII (Bg) and compared to the control without enzymatic digestion (C). HyperLadder™ 1 kb (BioLine) (M) was used as a weight marker. (B) SDS-PAGE analysis of partially purified Bp15.1PLPs particles along with the PageRuler Plus Prestained Protein Ladder (ThermoFisher) (M). The most abundant proteins (arrows 1 to 4) were excised and identified by MALDI-TOF.
Figure 9. (A) Agarose gel electrophoresis of the DNA contained in Bp15.1PLP before and after restriction enzyme digestion. The nucleic acids were extracted from partially purified Bp15.1PLPs after mitomycin C induction and DNAse treatment. The DNA was digested with BamHI B), HinDIII (H), SalI (S), EcoRI (E) and BglII (Bg) and compared to the control without enzymatic digestion (C). HyperLadder™ 1 kb (BioLine) (M) was used as a weight marker. (B) SDS-PAGE analysis of partially purified Bp15.1PLPs particles along with the PageRuler Plus Prestained Protein Ladder (ThermoFisher) (M). The most abundant proteins (arrows 1 to 4) were excised and identified by MALDI-TOF.
Ijms 22 08164 g009
Figure 10. Transmission Electron Microscopy of partially purified (sucrose gradient) bacteriophage-like particles. Apart from Bp15.1PLP, another less abundant phage-like particle was observed. Heads (H) and tails (T) of this particle are shown with arrowheads.
Figure 10. Transmission Electron Microscopy of partially purified (sucrose gradient) bacteriophage-like particles. Apart from Bp15.1PLP, another less abundant phage-like particle was observed. Heads (H) and tails (T) of this particle are shown with arrowheads.
Ijms 22 08164 g010
Table 1. Predicted ORFs in the 28.7 Kb region identified by PHASTER.
Table 1. Predicted ORFs in the 28.7 Kb region identified by PHASTER.
ORF DirectionPredicted Nº of aaPredicted FunctionModule
1-125XRE family transcriptional regulatorEarly
2+87XRE family transcriptional regulator
3+75Hypothetical protein
4+133Holliday junction resolvaseMiddle
5+68Hypothetical protein XtrA Domain
6+172Positive control sigma-like factor
7+214Small subunit of the terminase*
8+437Phage terminase large subunit
9+485portal protein superfamily
10+362Putative phage serine proteaseLate
11+307Phage major capsid protein
12+127Head-tail connector protein
13+118Tail shaft protein
14+165Putative tail-component
15+151Hypothetical protein
16+70Hypothetical protein
17+448Phage tail sheath protein
18+147Phage tail tube protein
19+73Hypothetical protein
20+146Tail assembly chaperone protein
21+1423TMP with transglycosylase domain
22+222LysM peptidoglycan binding domain
23+336Hypothetical protein
24+88DUF2577 Protein of unknown function
25+140DUF 2634Protein of unknown function
26+349baseplate J like protein
27+307DUF 2313 Family of phage tail proteins
28+127DUF3751 Phage tail-collar fibre protein
29+290Hypothetical protein
30+109XkdW protein
31+63XkdX protein
32+371PHA01818 superfamily
33+106Hypothetical protein
34+46XkdX protein
35+470Hypothetical protein
36+115Putative short tail fibre *
37+48XkdX protein
38+70Holin BhlA
39+87Holin
40+267N-acetylmuramoyl-L-alanineamidase
* Function predicted by PHASTER. +/− Indicate Forward/Backward direction of transcription, respectively. Proteins found in fingerprinting analysis are shown in bold.
Table 2. Strains and accession number of the highest hits were obtained to compare the 28.7 kb region from B. pumilus 15.1 with sequences at the Genbank. The coverage and identity of the found regions are also indicated. Genomes selected for whole-genome dot-plot comparison are shown in bold.
Table 2. Strains and accession number of the highest hits were obtained to compare the 28.7 kb region from B. pumilus 15.1 with sequences at the Genbank. The coverage and identity of the found regions are also indicated. Genomes selected for whole-genome dot-plot comparison are shown in bold.
Organism CoverIdentity Accession Number
1Bacillus sp. WP80.8698.01%CP010075.1
2Bacillus safensis strain U17-10.8597.14%CP015611.1
3Bacillus safensis strain U410.8597.14%CP015610.1
4Bacillus safensis strain KCTC 12796BP0.8595.98%CP018197.1
5Bacillus safensis strain U14-50.7995.95%CP015607.1
6Bacillus safensis strain FO-36b0.8495.93%CP010405.1
7Bacillus safensis strain BRM10.8395.58%CP018100.1
8Bacillus safensis strain sami0.8595.52%CP032830.1
9Bacillus safensis strain pgKB200.8995.52%CP043404.1
10Bacillus safensis strain IDN10.8794.79%AP021906.1
11Bacillus altitudinis strain Cr-10.8887.44%CP031774.1
12Bacillus pumilus strain C40.8787.30%CP011109.1
13Bacillus pumilus strain SH-B110.9486.86%CP010997.1
14Bacillus pumilus strain PDSLzg-10.8192.32%CP016784.1
15Bacillus pumilus strain ZB2017010.8392.30%CP029464.1
16Bacillus pumilus strain MTCC B60330.8987.92%CP007436.1
17Bacillus altitudinis strain P-100.7987.91%CP024204.1
18Bacillus altitudinis strain GQYP1010.8587.84%CP040514.1
19Bacillus pumilus strain NCTC103370.8391.72%LT906438.1
20Bacillus altitudinis strain W30.8687.80%CP011150.1
21Bacillus phage PBP1800.8188.05%KC847113.1
22Bacillus pumilus strain 1450.8291.54%CP027116.1
23Bacillus cellulasensis strain ku-bf10.7787.68%CP014165.1
24Bacillus cellulasensis strain NJ-V0.8287.63%CP012330.1
25Bacillus cellulasensis strain NJ-V20.8287.63%CP012482.1
26Bacillus cellulasensis strain NJ-M20.8287.63%CP012329.1
27Bacillus cellulasensis strain GLB1970.8387.61%CP018574.1
28Bacillus pumilus strain TUAT10.8387.61%AP014928.1
29Bacillus aerophilus strain 2320.8487.59%CP026008.1
30Bacillus altitudinis strain SGAir00310.8687.58%CP022319.2
31Bacillus altitudinis strain GR-80.8387.67%CP009108.1
32Bacillus altitudinis strain HQ-51-Ba0.7787.56%CP040747.1
33Bacillus pumilus strain SH-B90.7891.26%CP011007.1
34Bacillus altitudinis strain CHB190.7787.42%CP043559.1
35Bacillus pumilus strain SF-40.7891.01%CP047089.1
36Bacillus pumilus strain SAFR-0320.890.94%CP000813.4
37Bacillus pumilus strain 150a0.7590.93%CP027034.1
38Bacillus xiamenensis strain VV30.7788.58%CP017786.1
39Bacillus altitudinis strain FD480.887.77%CP025643.1
40Bacillus sp. ms-22 0.8288.48%CP046653.1
Table 3. Predicted ORFs in the 60.213 Kb region identified by PHASTER.
Table 3. Predicted ORFs in the 60.213 Kb region identified by PHASTER.
ORF DirectionPredicted Nº of aaPredicted Function
1+150Hypothetical protein
2+138Transcriptional regulator Helix-turn-helix domain
3+121Hypothetical protein
4+75Hypothetical protein
5+125Hypothetical protein
6+122Hypothetical protein
7+269Hypothetical protein
8+253DNA replication protein with ATP binding domain
9+216Hypothetical protein
10+463Phage DNA helicase
11+333DNA primase
12104XRE family transcriptional regulator
1368XRE family transcriptional regulator
14+388Hypothetical protein
15+176Hypothetical protein
16+167Hypothetical protein
17+257Hypothetical protein
18+83Hypothetical protein
19+144Hypothetical protein
20+237Hypothetical protein
21+759DNA Polimerase I
22+349Hypothetical protein
23+84Hypothetical protein
24+65Hypothetical protein
25+184Holiday junction resolvase
26+164Hypothetical protein
27+87Hypothetical protein
28+121Hypothetical protein
29+67Hypothetical protein
30+64Hypothetical protein
31+118Class 1b ribonucleoside-diphosphate reductase assembly flavoprotein NrdI
32+695Class 1b ribonucleoside-diphosphate reductase subunit alpha
33+324Class 1b ribonucleoside-diphosphate reductase subunit beta
34+64Hypothetical protein
35+70Hypothetical protein
36+108Deoxyuridine 5’-triphosphate nucleotide hydrolase
37+166Hypothetical protein
38+264FAD-dependent thymidylate synthase
39+164This protein contain the critical active site aspartate of mltA-like lytic transglycosylases
40+283Hypothetical protein
41+186Dephospho-CoA kinase
42124Hypothetical protein
43+124DUF134 domain
44+68Hypothetical protein
45+162Hypothetical protein
46+278MarR family transcriptional regulator
4761Hypothetical protein
48127Hypothetical protein
49+70XRE family transcriptional regulator
50+93Hypothetical protein
51+93Hypothetical protein
52+95Hypothetical protein
53+73Hypothetical protein
54+94Hypothetical protein
55+100Hypothetical protein
56181Phage integrase family
57+588Large subunit terminase
58+141Hypothetical protein with HTH TNP domain
59+533Phage portal protein
60+276Minor capsid protein *
61+227Scaffold protein *
62+363Major capsid protein *
63+73Hypothetical protein
64+128Hypothetical protein
65+112Hypothetical protein
66+136Bacteriophage putative tail-component
67+124Hypothetical protein
68+180Hypothetical protein
69+126Hypothetical protein
70+71Hypothetical protein
71+934Phage tail tape measure protein
72+475Phage tail protein
73+463Prophage endopeptidase tail
74+626right-handed parallel beta-helix repeat-containing protein
75+483Phage baseplate upper protein
76+95Hypothetical protein
77+46XkdX family protein
78+103Hypothetical protein
79+277N-acetylmuramoyl-L-alanineamidase
80+79Phage holin
81+183Hypothetical protein
82114YOID-like family protein
8376Hypothetical protein
8471HTH transcriptional regulator
85+263phage protein HTH domain containing
86353Site specificrecombinase/integraseXerD
* Function predicted by PHASTER. +/− Indicate Forward/Backward direction of transcription, respectively. The 31.5 Kb subregion (shadowed) and 28.7 Kb subregion (not shadowed) are indicated in the table. The protein found in fingerprinting analysis is shown in bold.
Table 4. Proteins identified by MALDI-TOF.
Table 4. Proteins identified by MALDI-TOF.
Protein Function
BandMass (kDa)MASCOT IdentificationPhaster ORFScoreNº MatchesBy PhasterBy Blastp
149,06Phage portal protein [Bacillus pumilus]PBSX-like ORF 17869XkdK-like tail sheath proteinXkdK-like tail sheath protein
240,12Hypothetical protein [Bacillus sp.]tRNA insert ORF 62919Major capsid proteinHypothetical protein
334,14Phage Major Capsid Protein [Bacillus pumilus]PBSX-like ORF 119510Capsid proteinPhage Major Capsid Protein
416,45Phage Portal Protein [Bacillus sp.]PBSX-like ORF 1811512phage portal proteinPhage tail tube protein
Table 5. Number of phages, with the available genome at NCBI *, described using a Bacillus strain as host (in alphabetical order).
Table 5. Number of phages, with the available genome at NCBI *, described using a Bacillus strain as host (in alphabetical order).
Bacillus SpecieNo. PhagesBacillus SpecieNo. Phages
B. alcalophilus1B. pumilus14
B. anthracis3B. safensis1
B. bogoriensis1B. subtilis21
B. cereus31B. thuringiensis44
B. clarkii1B. velezensis1
B. halmapalus1B. weihenstephanensis1
B. licheniformis1Bacillus sp. 7
B. megaterium17
* NCBI: https://www.ncbi.nlm.nih.gov/genome/viruses/; last access 1 May 2021.
Table 6. B. pumilus phages found in viral genome database at NCBI *. The name of the phage, the strain used for multiplication, the access number and genome size are detailed.
Table 6. B. pumilus phages found in viral genome database at NCBI *. The name of the phage, the strain used for multiplication, the access number and genome size are detailed.
Phage NameB. pumilus(B. p) StrainAcc. NumberGenome Size
BobbB. p SAFR32NC_024792160,281
Bp8p-CB. p GR8NC_029121151,417
Bp8p-TB. p GR8NC_047744151,419
phiGATEB. p GL1NC_020081149,844
Bpu_PumAB. pNC_04997118,466
Bpu_PumBB. pNC_04997218,932
AndromedaB. p BL8NC_02047849,259
BlastoidB. p BL8NC_02277350,354
CurlyB. p BL8NC_02047949,425
EoghanB. p BL8NC_02047749,458
FinnB. p BL8NC_02048050,161
GlitteringB. p BL8NC_02276649,246
RiggiB. p BL8NC_02276549,836
TaylorB. p BL8NC_04185849,492
* NCBI: https://www.ncbi.nlm.nih.gov/genome/viruses/; last access 1 May 2021.
Table 7. Activity of Bp15.1PLP on bacterial strains from a laboratory collection.
Table 7. Activity of Bp15.1PLP on bacterial strains from a laboratory collection.
Strain/GroupLysis ActivityInhibition Halo Morphology
B. pumilus 15.1
B. pumilus 15.1C+++Translucent
B. pumilus 8A3+++Translucent
B. pumilus M1+++Translucent
B. pumilus M2+++Translucent
B. subtilis sp.
B. subtilis 168C
Bacillus thuringiensis 78/11
Bacillus cereus 169
Latococcus garvaeae
Lactococcus NZ9000
Escherichia coli DH5α
Arthrobacter sp.
Listeria innocua CECT 4030
Micrococcus sp.
(−) indicates no lysis activity; (+++) indicates strong lysis activity.
Table 8. Mortality results, obtained after 10 days in C. capitata larvae bioassays using sporulated cultures of B. pumilus 15.1 and different concentrations of Bp15.1PLPs. The increase in toxicity compared to the negative control (water) was also calculated (Fold increase column).
Table 8. Mortality results, obtained after 10 days in C. capitata larvae bioassays using sporulated cultures of B. pumilus 15.1 and different concentrations of Bp15.1PLPs. The increase in toxicity compared to the negative control (water) was also calculated (Fold increase column).
Bioassay% MortalityFold Increase
H2O61.0
SM buffer122.0
Sporulated culture203.3
Bp15.1PLP (Lysate)00.0
Sporulated culture + Bp15.1PLP (Lysate)142.3
Bp15.1PLP (200x concentrated)122.0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fernández-Fernández, A.; Osuna, A.; Vilchez, S. Bacillus pumilus 15.1, a Strain Active against Ceratitis capitata, Contains a Novel Phage and a Phage-Related Particle with Bacteriocin Activity. Int. J. Mol. Sci. 2021, 22, 8164. https://doi.org/10.3390/ijms22158164

AMA Style

Fernández-Fernández A, Osuna A, Vilchez S. Bacillus pumilus 15.1, a Strain Active against Ceratitis capitata, Contains a Novel Phage and a Phage-Related Particle with Bacteriocin Activity. International Journal of Molecular Sciences. 2021; 22(15):8164. https://doi.org/10.3390/ijms22158164

Chicago/Turabian Style

Fernández-Fernández, Alberto, Antonio Osuna, and Susana Vilchez. 2021. "Bacillus pumilus 15.1, a Strain Active against Ceratitis capitata, Contains a Novel Phage and a Phage-Related Particle with Bacteriocin Activity" International Journal of Molecular Sciences 22, no. 15: 8164. https://doi.org/10.3390/ijms22158164

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop