# The Revisited Genome of Bacillus subtilis Bacteriophage SPP1

Lia M. Godinho, Mehdi El Sadek Fadel, Céline Monniot, Lina Jakutytè, Isabelle Auzat, Audrey Labarde, Karima Djacem, Leonor Oliveira, Rut Carballido-Lopez, Silvia Ayora, et al.

► To cite this version:

Lia M. Godinho, Mehdi El Sadek Fadel, Céline Monniot, Lina Jakutytè, Isabelle Auzat, et al. The Revisited Genome of Bacillus subtilis Bacteriophage SPP1. Viruses, 2018, 10 (12), pp.705. 10.3390/v10120705. hal-02178891

HAL Id: hal-02178891

https://hal.science/hal-02178891

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# The Revisited Genome of *Bacillus subtilis*
## Bacteriophage SPP1

Lia M. Godinho<fnref n="1" />, Mehdi El Sadek Fadel<fnref n="1" />, Céline Monniot<fnref n="2" />, Lina Jakutyte<fnref n="3" />, Isabelle Auzat<fnref n="1" />, Audrey Labarde<fnref n="1" />, Karima Djacem <fnref n="1" />,†, Leonor Oliveira<fnref n="1" />, Rut Carballido-Lopez<fnref n="2" />, Silvia Ayora <fnref n="4" /> and Paulo Tavares <fnref n="1" />,*

Received: 23 October 2018; Accepted: 6 December 2018; Published: 11 December 2018

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** *Bacillus subtilis* bacteriophage SPP1 is a lytic siphovirus first described 50 years ago. Its complete DNA sequence was reported in 1997. Here we present an updated annotation of the 44,016 bp SPP1 genome and its correlation to different steps of the viral multiplication process. Five early polycistronic transcriptional units encode phage DNA replication proteins and lysis functions together with less characterized, mostly non-essential, functions. Late transcription drives synthesis of proteins necessary for SPP1 viral particles assembly and for cell lysis, together with a short set of proteins of unknown function. The extensive genetic, biochemical and structural biology studies on the molecular mechanisms of SPP1 DNA replication and phage particle assembly rendered it a model system for tailed phages research. We propose SPP1 as the reference species for a new SPP1-like viruses genus of the *Siphoviridae* family.

**Keywords:** Bacteriophage; *Siphoviridae*; SPP1; *Bacillus subtilis*; genome organization; DNA replication; virus assembly; virus DNA packaging; virus evolution

## 1. Introduction

The lytic bacteriophage SPP1 (Subtilis Phage Pavia 1) that infects the soil bacterium *Bacillus subtilis* was isolated in the Botanical Garden of Pavia in Italy [1]. The capacity of its purified DNA to transfect *B. subtilis* competent cells and the difference of density between its two DNA strands, rendering easy their separation following denaturation, were quickly recognized. They attracted initial interest to this phage system for studying DNA uptake into bacteria [2,3], mismatch repair [2,3] and recombination [4–6]. Establishment of the genetic [7] and physical [8,9] maps of the SPP1 genome then paved the way for studies of phage gene expression, DNA replication and assembly of the viral particle. The intensive research that followed on the molecular mechanisms supporting SPP1 infection rendered it one of the best characterized phages of Gram-positive bacteria. SPP1 generalized transduction was

<footnote n="1">¹ Institut de Biologie Intégrative de la Cellule (I2BC), French Alternative Energies and Atomic Energy Commission (CEA), Centre National de la Recherche Scientifique (CNRS), Univ Paris-Sud, Université Paris-Saclay, 91190 Gif-sur-Yvette, France; lia.marques-godinho@i2bc.paris-saclay.fr (L.M.G.); mehdi.elsadekfadel@i2bc.paris-saclay.fr (M.E.S.F.); Isabelle.AUZAT@i2bc.paris-saclay.fr (I.A.); audrey.labarde@i2bc.paris-saclay.fr (A.L.); karima.djacem@gmail.com (K.D.); Leonor.OLIVEIRA@i2bc.paris-saclay.fr (L.O.)</footnote>

<footnote n="2">² MICALIS, Institut National de la Recherche Agronomique (INRA), AgroParisTech, Université Paris-Saclay, 78350 Jouy-en-Josas, France; celine.monniot@jouy.inra.fr (C.M.); rut.carballido-lopez@inra.fr (R.C.-L.)</footnote>

<footnote n="3">³ Unité de Virologie Moléculaire et Structurale (VMS), CNRS, 91198 Gif-sur-Yvette, France; lina.jakutyte@gmail.com</footnote>

<footnote n="4">⁴ Centro Nacional de Biotecnología (CNB-CSIC), 28049 Madrid, Spain; sayora@cnb.csic.es</footnote>

<footnote>* Correspondence: paulo.tavares@i2bc.paris-saclay.fr</footnote>

<footnote>† Current address: Folium Science, Unit DX, St Philips Central, Albert road, Bristol BS2 0XJ, UK.</footnote>
also instrumental for fine genetic mapping of the *B. subtilis* chromosome, strain construction, and plasmid horizontal transfer [5,10–14].

SPP1 belongs to the family *Siphoviridae*. Its viral particle has an isometric icosahedral capsid with a diameter of 61 nm [15] and a 190 nm-long non-contractile tail [16]. The linear double-stranded DNA (dsDNA) molecule contained in the viral capsid has a length of ~45.9 kbp with a variation in size of ~2.5 kbp [17,18]. The phage genome size is 44,016 bp. The DNA molecules in a phage population are terminally redundant and partially circularly permuted, resulting from a headful packaging mechanism [19,20].

SPP1 infection is initiated by reversible adsorption of the viral particle to glycosylated teichoic acids [21] followed by irreversible binding of SPP1 to YueB. This integral membrane protein is a component of a type VII secretion system [22] that crosses the *B. subtilis* cell wall to be exposed at the bacterial surface [23–25]. The interaction of the SPP1 tail fiber with YueB triggers ejection of SPP1 DNA through its tail tube, committing the phage to infection. Phage DNA circularizes in the bacterial cytoplasm most probably by recombination between its redundant ends [26–30]. DNA replication then ensues in a discrete position of the bacterial cytoplasm [24]. DNA synthesis was proposed to initiate by theta replication of circular molecules followed by a switch to rolling circle (sigma) replication that generates concatemers of the SPP1 genome [30,31]. These are the substrate for phage DNA encapsidation into a preformed procapsid structure. Packaging is initiated by specific recognition and cleavage of a *pac* sequence within the SPP1 genome, followed by its translocation into the procapsid interior through a specialized portal vertex. Encapsidation is terminated by an imprecise [18], sequence-independent, endonucleolytic cleavage of the substrate concatemer when a threshold amount of DNA is reached inside the capsid (headful packaging mechanism) [19,20]. Subsequent encapsidation cycles follow processively along the concatemer. After packaging, DNA is retained inside the capsid by proteins that close the portal vertex and build the interface for attachment of the phage tail. Binding of the tail, which is assembled in an independent pathway, yields the infectious particle (virion). The SPP1 tail is formed by the adsorption apparatus, responsible for interaction with the host cell surface, which is connected by a 160-nm long helical tail tube to the tail tapered end that binds to the capsid portal [16,32]. Virions accumulate in the infected bacterium until lysis. Lysis is promoted by membrane proteins (holins) that concentrate in the cytoplasmic membrane leading to its disruption and release of an endolysin that digests the *B. subtilis* cell wall [33].

Here we revisit the sequence and organization of the bacteriophage SPP1 genome. A detailed updated annotation of its genes combining available experimental data and bioinformatics completes this comprehensive study of SPP1 genetics and biology.

## 2. Materials and Methods

The SPP1 genome was fully sequenced using Illumina sequencing of genome libraries at the I2BC NGS platform facility, as described [34,35]. Sanger sequencing of PCR fragments covering the SPP1 genome was carried out at GATC Biotech (Germany).

Open reading frames (ORFs) identification and visualization of the SPP1 chromosome organization was performed using the SnapGene Viewer software (GSL Biotech, Chicago, IL) and Fgenesb annotator (http://www.softberry.com/berry.phtml?topic=fgenesb&group=programs&subgroup=gfindb) [36]. A cut-off of ORFs initiated at AUG, GUG or UUG codons that code for putative proteins of at least 45 amino acids-long was used as criterion to identify the ORFs for annotation. Among those, there are 20 ORFs without ribosome binding site (RBS) that are embedded in longer ORFs. Such internal ORFs were eliminated from our downstream analysis because they are most likely not translated. The strong distribution bias of ORFs coded by the SPP1 DNA heavy strand together with transcriptional and functional data (compiled in [26]) support that this strand is used, possibly exclusively, as the coding DNA template. Therefore, the analysis described here is limited to 80 ORFs, defined according to the above criteria, which are encoded by the heavy strand. The assignments made in [26] and a survey of SPP1 functional studies were used as the primary information for ORF
annotation. Visual inspection of the ORF 5' region was used to confirm the initiation codon and to assess if it is preceded by an RBS whose sequence is complementary to the *B. subtilis* 16 S rRNA 3' sequence UCUUUCCUCCACUAG [37]. A spacing of 8-14 nucleotides was allowed between the center of the RBS (complementary to the U underlined) and the nucleotide preceding the initiation codon. Sequence and structural homology searches of the ORFs encoded proteins were carried out with BLASTp [38,39] and HHPred [40], respectively. Protein properties were analyzed with the Expasy ProtParam tool (https://web.expasy.org/protparam/). Prediction of transmembrane helices was carried out with the TMHMM server v.2.0 (http://www.cbs.dtu.dk/services/TMHMM/).

SPP1 putative early promoters were predicted using BPROM. BPROM recognizes promoter sequences bound by bacterial sigma factors of the $\sigma^{70}$ family with about 80% accuracy and specificity (http://www.softberry.com/berry.phtml?topic=bprom&group=programs&subgroup=gfindb) [36]. These include promoters recognized by the primary sigma factor $\sigma^A$ of *B. subtilis*. Promoter sequence prediction was carried out within the regions where early promoters were previously mapped [26,41,42]. The -10 and -35 sequences of late promoter PL1 were used for word scanning of the SPP1 genome to search for putative late promoters without success. Rho-independent transcriptional terminators were determined using ARNold (http://rna.igmors.u-psud.fr/toolbox/arnold/index.php) [43-46], performing a whole sequence search analysis on the coding strand using two complementary programs, Erpin and RNAmotif. The free energy ($\Delta G^\circ$) of the predicted terminator stem-loop structure was computed with ARNold [44].

Codon usage bias of the 80 SPP1 ORFs was calculated using the Codon Usage Calculator from Biologics International Corp (https://www.biologicscorp.com/tools/CodonUsageCalculator/). The *B. subtilis* genome codon usage frequencies were obtained from [47,48].

The revised version of the complete SPP1 genome sequence and its annotation are in the process of submission to Genbank with accession code "X97918.3". Table 1 summarizes the changes identified relative to the previous sequence.

<table><caption>**Table 1.** Revision of the SPP1 sequence. Changes in the GenBank X97918.2 sequence relative to the GenBank X97918.3 revised version are listed. The position coordinates presented are the ones found in the GenBank X97918.3 sequence. The ORF concerned, the nucleotide sequence change(s), and the effect(s) on the gene product amino acid sequence are displayed.</caption><thead><tr><th>Position (nt)</th><th>ORF</th><th>ORF Sequence Change</th><th>ORF Product Amino Acid Change</th></tr></thead><tbody><tr><td>230</td><td>1</td><td>insertion of a G; frameshift</td><td>gp1 Cter sequence changed and shortened</td></tr><tr><td>766</td><td>2</td><td>T→C</td><td>none (silent mutation)</td></tr><tr><td>1997</td><td>4</td><td>insertion of an A; frameshift</td><td>gp4 Cter changed and lengthened</td></tr><tr><td>5001</td><td>intergenic region ORFs 7 and 8</td><td>insertion of a T</td><td>none</td></tr><tr><td>23,794</td><td>30</td><td>deletion of a C; frameshift</td><td>gp30 Cter changed and lengthened</td></tr><tr><td>34,450</td><td>37.1</td><td>T→G</td><td>gp37.1 C<sub>18</sub>→G<sub>18</sub></td></tr><tr><td>37,304 and 37,306</td><td>40</td><td>GGC→CGG</td><td>gp40 R<sub>194</sub>R<sub>195</sub>→P<sub>194</sub>G<sub>195</sub></td></tr><tr><td>38,371 to 38,374</td><td>41</td><td>GTGT→TGTG</td><td>gp41 K<sub>101</sub>C<sub>102</sub>→N<sub>101</sub>V<sub>102</sub></td></tr><tr><td>40,240</td><td>44 and former 45</td><td>insertion of a G; frameshift; ORFs fusion</td><td>gp44 Cter lengthened by former gp45 sequence</td></tr><tr><td>41,787</td><td>48</td><td>insertion of an A; frameshift</td><td>gp48 changed and shortened</td></tr><tr><td>42,750</td><td>51</td><td>insertion of an A; frameshift</td><td>gp51 changed and lengthened</td></tr><tr><td>42,778 to 42,781</td><td>51</td><td>CGCG→GCGC</td><td>gp51 R<sub>56</sub>E<sub>57</sub>→A<sub>56</sub>Q<sub>57</sub></td></tr><tr><td>42,819</td><td>51</td><td>insertion of a G; frameshift</td><td>gp51 Cter lengthened</td></tr></tbody></table>
# 3. Results and Discussion

## 3.1. Properties of the SPP1 DNA Molecule

The genome of SPP1 was completely sequenced in the context of research projects on phage DNA replication and viral particle assembly, as well as by the necessity for an accurate genome annotation for downstream SPP1 "omics" studies. Very few changes were identified relative to the deposited sequence (Genbank accession code X97918) confirming the high quality of the original Sanger sequencing work [26] and its 2006 revision (X97918.2; [16]). Changes in the nucleotide sequence are listed in Table 1. The 44,016 bp genome has a GC content of 43.7% which is similar to the one found in the host *B. subtilis* genome (43.5%) [49]. The distribution of purines and pyrimidines in the two DNA strands is very different. This asymmetry provided the physical basis for their separation after denaturation by density using isopycnic centrifugation [1]. The heavy strand, whose purine (dA+dG) content is 58.4%, has a density (ρ) of 1.725 g/cm³. It is the SPP1 genome coding strand. The light chain density is 1.713 g/cm³ [1,3]. The strongest sequence bias are dA+dT rich-islands found in intergenic regions involved in transcriptional regulation (see Section 3.2) and in the two origins of SPP1 DNA replication *oriR* and *oriL* (see Section 3.5) (Figure 1).

<figure><img src="image_4.png" /></figure>

**Figure 1.** Organization of the SPP1 genome. The continuous bar represents the 44,016 bp-long genome where coordinate 1 is the main *pac* cleavage position [35]. The two origins of replication *oriR* and *oriL* (magenta), the DNA packaging signal *pac* (black) and non-essential regions of the SPP1 genome defined by deletions (dashed) are highlighted in the bar. The sequence inverted in SPP1 *invmir* is displayed by a pink line underneath the genome bar. The position of promoters (Figure 2) and potential Rho-independent transcriptional terminators that form stem loops (red) in mRNA (Table 2) is displayed on top of the bar. Transcription is from left to right. DNA packaging initiated at *pac* occurs in the same direction (arrow on the top left). The set of SPP1 genes and ORFs, identified as described in Materials and Methods (see Section 2), are presented above the genome bar and colored according to their function assignment shown on the bottom legend.
## 3.2. Transcription and Translation of the SPP1 Genome

Transcription of the SPP1 genome was reported to occur from the DNA heavy strand [1,50–52]. This is not an absolute requirement, as phage SPP1 *invmir* carries an inversion within its genome that leads to transcription from the light chain of a large genome segment which includes the DNA replication genes (Figure 1) [53]. The temporal gene expression program is defined by early and late transcription carried out by the *B. subtilis* RNA polymerase [54]. Five promoters with a canonical sequence recognized by *B. subtilis* σA direct transcription of early ORFs/genes (we name gene the ORFs for which a function was demonstrated experimentally or supported by substantial bioinformatics data) (Figures 1 and 2) [33,41,55,56]. These polycistronic transcriptional units are bracketed by the promoter and by putative Rho-independent transcriptional terminators coding for potential RNA stem-loop-forming sequences (Figure 1 and Table 2). The two strongest promoters (PE3 and PE2) drive transcription of gene sets that include essential DNA replication functions (Figure 1) [41,52,56]. PE1, also a strong promoter, and the weaker PE4 and PE5 promoters control expression of less characterized ORFs (see Section 3.4). Early genes/ORFs are encoded contiguously in the SPP1 genome. Transcription of the late genes segment requires translation of (an) early, yet unidentified, SPP1 factor(s) [55,57]. The only late promoters experimentally characterized are the adjacent PL1 and PL2 (Figure 1). PL1, which accounts for more than 95% of the transcriptional activity of the region upstream from gene 1, has a canonical –10 sequence and an atypical –35 region ([57]; Figure 2). Several clusters of late genes are delimited by potential stem-loop transcription terminator sequences ([26]; Figure 1), mostly found in intergenic regions, but no consensual late promoter sequence related to PL1 was identified.

<table><thead><tr><th></th><th>-35</th><th>-10</th><th>References</th></tr></thead><tbody><tr><td>PE1 (40523) TAAAAAGTGTTCACAACGTAACCCTAAGATGAFATTATATGATTACAACGTAACCACGAGAGGGGAAATTAAAAATG (40598)</td><td></td><td></td><td>[55]</td></tr><tr><td>PE2 (35313) AAAAAACTATAGACAAGCTTAGACATTTGAGTTAAGATATCTTTAGATTCAGTTGGACATACTTAGACATGCAGGAGGAAAAAAGGAAATG (35402)</td><td>+1</td><td></td><td>[56]</td></tr><tr><td>PE3 (30301) AAATAGTGGTCCCTTCTATGTTTTCTATGTTTAAAGAATCATAGAGAGGGGGGACAACATG (30366)</td><td></td><td></td><td>[42,55]</td></tr><tr><td>PE4 (27975) ATTTTAGGGGTTACGTGTGAATCATTTGTAGTCATACTGTGTCAAGGGGGAGGTG GAAAACCCTTG (28041)</td><td></td><td></td><td>[26]</td></tr><tr><td>PE5 (20757) ATTTCCACGCTCTCAACTGTTGTAGTAGGTGCGTATACITFACAATTGTATGGAACCTTGTAAAGAACTAATGAAGGAGAGTGAAAGAATATG (20847)</td><td></td><td></td><td>[33]</td></tr><tr><td>PL1 (43785) GCACCTATTTGGGTGCTTTTTTGTTGTATAATTAGGTTTATATAAGGTTTTATCATTCAGATATGAGGTTCAAATATAGTTTTAAGGAGGTTTTTTCTATG (43886)</td><td>+1</td><td></td><td>(PL1 is the stronger of the two promoters, accounting for 95% of transcripts of gene 1 operon) [26,57]</td></tr><tr><td>PL2 (43797) GGTGCTTTTTTGTTGTATAATTAGGTTTATATAAGGTTTTATCATTCAGATATGAGGTTCAAATATAGTTTTAAGGAGGTTTTTTCTATG (43886)</td><td>+1</td><td></td><td></td></tr></tbody></table>

**Figure 2.** SPP1 promoters. The sequence of SPP1 promoters and the initiation codon of their downstream gene (double underline) are displayed. The –35 and –10 promoter regions are shaded in grey for promoters whose transcription start position (+1) was determined experimentally. A dashed box denotes the atypical –35 sequence of PL1. SPP1 putative early promoters identified by sequence similarity to *B. subtilis* vegetative promoters are highlighted in black with sequence characters in white. Note that the approximate position of transcription initiation and promoter strength was determined for all early promoters by electron microscopy of DNA-RNA polymerase complexes [41,42].
<table><caption>**Table 2.** Rho-independent transcriptional terminators of SPP1. The terminators were identified as described in Material and Methods (see Section 2). Putative stem (underlined) and loop (bold) regions are highlighted. The coordinates of the sequence presented are shown (nt—nucleotide). The ΔG° of the RNA stem loop is shown on the right column.</caption><thead><tr><th>Terminator Sequence (Loop, Stem)</th><th>Start (nt)</th><th>End (nt)</th><th>ΔG° (kcal.mol<sup>-1</sup>)</th></tr></thead><tbody><tr><td>AAAAGCGAGTTAACCGACGTAAAAATGCGTCGGTTTTTTTCGTGCTTC</td><td>4718</td><td>4765</td><td>-12.20</td></tr><tr><td>GTGTCAGTGCGCGGGTTCAATTCCCGGAGTCCGTTTTTACCGCCC</td><td>6780</td><td>6824</td><td>-5.50</td></tr><tr><td>TAAAAAGAGAAGAGGGGCTAACGCCTCTCTTTTTTGAAAG</td><td>8775</td><td>8815</td><td>-14.00</td></tr><tr><td>ATTGCCAGCAGAGAGCACGGGTTAATTCCCGCGCTCTTTTTTGTATTCA</td><td>11,255</td><td>11,304</td><td>-9.20</td></tr><tr><td>TGAATTAGACAGGGCCGCGCAAGTGGCTCTTTTTAATAGGT</td><td>12,009</td><td>12,049</td><td>-15.20</td></tr><tr><td>TGAAAAGACGCGCGCCGCTAACGGCGTCCGTTTTGAACATGA</td><td>16,711</td><td>16,753</td><td>-12.10</td></tr><tr><td>GATTGACGAAGTAAAGGGCCATGTGCCCTTTATTTTTTGCAAA</td><td>21,897</td><td>21,940</td><td>-10.70</td></tr><tr><td>ACCTTCGCTTGCCGCCCGGCTGATGGGCGGTTTTTATTTTT</td><td>24,095</td><td>24,136</td><td>-13.20</td></tr><tr><td>TGTCCTAAAATCGGCCCCGTTCCCAGTCGGGCCACTTTTTTATTTTA</td><td>27,934</td><td>27,980</td><td>-13.60</td></tr><tr><td>GACTTTGAAAGGAACCGTTCTCTAACGGTTCTTTTTTATTTTC</td><td>30,205</td><td>30,247</td><td>-9.90</td></tr><tr><td>TATTTTGATTGAGTCGGGGAAACCCGGCTTTTTTATTTGGG</td><td>35,232</td><td>35,273</td><td>-14.90</td></tr><tr><td>ATCAAAGTTATGGTGGGAGTAATCCCGCCTTTTTCTATTTT</td><td>40,466</td><td>40,506</td><td>-14.40</td></tr><tr><td>AACACAGAGAGGCACCCTATTTGGGTGCTTTTTTGTTGTA</td><td>43,774</td><td>43,813</td><td>-11.90</td></tr></tbody></table>

Genome-wide sequence analysis and subsequent individual inspection identified 80 ORFs encoded by the SPP1 wild type DNA heavy chain with a length longer than 45 codons (135 bp), as described in Material and Methods (see Section 2) (Table 3; Figure 1). A stringent criterion was used to annotate only ORFs that are most likely translated, an assignment further supported by the subsequent search for their 5' RBS (see below). The ORFs defined cover 94% of the SPP1 sequence. Non-coding segments longer than 45 bp were characterized in most cases by the presence of transcriptional promoters and/or terminators sequences, but were also found in five cases between adjacent genes within a transcriptional unit (genes 24–24.1, 26.1–27, 31–31.1, 46.1–47, 48–49). Only one ORF longer than 15 codons was identified in these intergenic regions. It has 33 codons starting by AUG and preceded by an RBS in the segment between ORFs 46.1 and 47 (not shown). It cannot be excluded at present that other regions of the genome code polypeptides smaller than 45 amino acids-long, as found for several phages [58,59]. We are pursuing proteomic and genetic studies for an accurate annotation of short SPP1 ORFs. An RBS whose sequence is complementary to the *B. subtilis* 16 S rRNA 3' was found in the 5' region of 77 from the 80 ORFs investigated here (Table 3). No RBS was identified for genes 2 and 16 that encode well-characterized proteins essential for phage particle assembly (Table 3). Note that in case of gene 16 the initiation codon was confirmed experimentally by amino-terminus sequencing of its encoded protein [60]. The initiation codon of gene 2 overlaps the stop codon of its 5' gene 1 (AUGA; the initiation codon is shown in bold and the termination codon in double underlined) while the initiation codon of gene 16 is spaced by one nucleotide from the gene 15 termination codon (UAAGAUG). This organization possibly allows coupling initiation of translation of the genes lacking an RBS with translational termination of their upstream gene [61]. It might also provide a strategy to control the ratio between proteins encoded by adjacent genes. Those proteins interact directly during phage assembly: gene product 1 (gp1 (note that the designations gpX (gene product X) and GXP (gene X product) are synonymous in the SPP1 literature)) and gp2 form the DNA packaging terminase while gp15 and gp16 bind sequentially to the capsid portal after viral DNA packaging (see Section 3.8). ORF 42.2 also lacks a canonical RBS. Its initiation codon overlaps the termination codon of ORF 42.1 (AUGA) and is preceded by a GGGG sequence that could act as a weak RBS (Table 3), probably ensuring gp42.2 synthesis. It is possible that other SPP1 ORFs without RBS, a situation found for ~10% of the host *B. subtilis* ORFs [62], are expressed during SPP1 infection. These can be either smaller than 45 codons or embedded within SPP1 annotated ORFs (see Section 2).
<table><caption>**Table 3.** SPP1 ORFs. The 80 ORFs of the SPP1 transcribed heavy chain were assigned as described in Material and Methods (see Section 2) and numbered following the original nomenclature of Alonso et al. [26]. The RBS are sequences complementary to the *B. subtilis* 16 S rRNA 3' sequence (mismatches are shown in small case; the position after which spacing is calculated between the RBS and the nucleotide preceding the initiation codon is underlined; n.d.—not determined). The ORFs coding regions coordinates are listed. Genes are defined essential when their inactivation in conditional lethal mutants prevents phage multiplication. The length, molecular mass (MM) and presence of putative transmembrane segments in proteins translated from the ORFs coding frame are listed on the right side of the Table. The coordinates of X-ray crystallography, NMR or cryo-electron microscopy structures available for individual proteins or their complexes are also provided. The function of individual proteins is based on experimental data while putative function is deduced from bioinformatics analysis. Structural components of the SPP1 viral particle (st.) are assigned based on biochemical, structural and/or robust bioinformatics data. Proteins are grouped according to function following the color code used in Figure 1.</caption><thead><tr><th rowspan="2">ORF</th><th rowspan="2">RBS (mRNA)</th><th colspan="7">ORF</th><th colspan="4">Protein</th><th rowspan="2">Protein Function</th><th rowspan="2">References</th></tr><tr><th>Spacing</th><th>Start</th><th>nt</th><th>Stop</th><th>nt</th><th>Essential Gene</th><th>Length (aa)</th><th>MM (kDa)</th><th>Predicted TMM Segments</th><th>3D Structure (PDB or EMD)</th></tr></thead><tbody><tr><td>1</td><td>AAGGAGGU</td><td>10</td><td>AUG</td><td>43,884</td><td>UGA</td><td>311</td><td>yes</td><td>147</td><td>16.3</td><td>no</td><td>3ZQQ<sup>a</sup> (Xtal<sup>b</sup>)</td><td>small terminase subunit (TerS)</td><td>[57,63–69]</td></tr><tr><td>2</td><td>n.d.<sup>c</sup></td><td></td><td>AUG<sup>d</sup></td><td>308<sup>d</sup></td><td>UAG</td><td>1576</td><td>yes</td><td>422<sup>d</sup></td><td>48.8<sup>d</sup></td><td>1<sup>e</sup></td><td>2WBN<sup>f</sup>; 2WC9<sup>f</sup> (Xtal<sup>b</sup>)</td><td>large terminase subunit (TerL)</td><td>[57,67,70–74]</td></tr><tr><td>3</td><td>AAAGGAGG</td><td>11</td><td>AUG</td><td>1567</td><td>UAA</td><td>1782</td><td>no</td><td>71</td><td>8.5</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>4</td><td>GGgGGU</td><td>10</td><td>AUG</td><td>1782</td><td>UAA</td><td>2072</td><td>no</td><td>96</td><td>11.4</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>5</td><td>AAAGAGG</td><td>11</td><td>AUG</td><td>2065</td><td>UGA</td><td>2334</td><td>no</td><td>89</td><td>10.3</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>6</td><td>AGGAGGU</td><td>11</td><td>AUG</td><td>2336</td><td>UGA</td><td>3847</td><td>yes</td><td>503</td><td>57.3</td><td>no</td><td>2JES (Xtal<sup>b</sup>); 5A20, 5A21 (cryoEM<sup>g</sup>)</td><td>st.; portal protein</td><td>[20,60,67,73–82]</td></tr><tr><td>7</td><td>AGGAGG</td><td>12</td><td>AUG</td><td>3804</td><td>UAA</td><td>4730</td><td>no</td><td>308</td><td>35.1</td><td>no</td><td>n.d.</td><td>st.; initiation of infection; binds to portal</td><td>[83,84]</td></tr><tr><td>8</td><td>AAAGGAG</td><td>12</td><td>AUG</td><td>5067</td><td>UGA</td><td>6215</td><td>n.d.</td><td>382</td><td>43.7</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>9</td><td>AAcGGAGG</td><td>9</td><td>AUG</td><td>6217</td><td>UAA</td><td>6555</td><td>n.d.</td><td>112</td><td>12.6</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>10</td><td>GGuGGUG<sup>h</sup></td><td>12<sup>h</sup></td><td>AUG</td><td>6583</td><td>UAG</td><td>6750</td><td>n.d.</td><td>55</td><td>6.2</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>11</td><td>AGGAG</td><td>9</td><td>AUG</td><td>6917</td><td>UAA</td><td>7561</td><td>yes</td><td>214</td><td>23.4</td><td>no</td><td>n.d.</td><td>procapsid scaffolding protein</td><td>[85–87]</td></tr><tr><td>12</td><td>AAGG<sub>g</sub>GG</td><td>11</td><td>AUG</td><td>7576</td><td>UAA</td><td>7770</td><td>no</td><td>64</td><td>6.6</td><td>no</td><td>n.d.</td><td>st.; capsid accessory protein with collagen-like fold</td><td>[15,88]</td></tr><tr><td>13</td><td>AAAGGAG</td><td>9</td><td>AUG</td><td>7803</td><td>UAA</td><td>8777</td><td>yes</td><td>324</td><td>35.3</td><td>no</td><td>4AN5 (cryoEM<sup>g</sup>)</td><td>st.; major capsid protein (MCP)</td><td>[15,85,86]</td></tr><tr><td>14</td><td>AAAGGAG</td><td>10</td><td>AUG</td><td>8828</td><td>UGA</td><td>9004</td><td>no</td><td>58</td><td>6.7</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>15</td><td>AAuGAGG</td><td>10</td><td>AUG</td><td>9015</td><td>UAA</td><td>9323</td><td>yes</td><td>102</td><td>11.6</td><td>no</td><td>2KBZ (NMR<sup>i</sup>); 5A20, 5A21 (cryoEM<sup>g</sup>)</td><td>st.; connector adaptor protein</td><td>[60,78,81,89]</td></tr><tr><td>16</td><td>n.d.<sup>c</sup></td><td></td><td>AUG<sup>j</sup></td><td>9325<sup>j</sup></td><td>UAG</td><td>9654</td><td>yes</td><td>109<sup>j</sup></td><td>12.5<sup>j</sup></td><td>no</td><td>2KCA (NMR<sup>i</sup>); 5A20, 5A21 (cryoEM<sup>g</sup>)</td><td>st.; connector stopper protein</td><td>[60,78,81,89]</td></tr><tr><td>16.1</td><td>AAaGAGG</td><td>11</td><td>AUG</td><td>9644</td><td>UGA</td><td>10,069</td><td>n.d.</td><td>141</td><td>15.9</td><td>no</td><td>n.d.</td><td>putative tail protein</td><td></td></tr><tr><td>17</td><td>AGGAGGU</td><td>10</td><td>AUG</td><td>10,066</td><td>UGA</td><td>10,470</td><td>yes</td><td>134</td><td>15</td><td>no</td><td>2LFP (NMR<sup>i</sup>)</td><td>st.; tail-to-head joining protein (THJP)</td><td>[32,90]</td></tr><tr><td>17.1</td><td>AGGAGG</td><td>10</td><td>AUG</td><td>10,484<sup>k</sup></td><td>UAA</td><td>11,017</td><td>yes</td><td>177</td><td>19.2</td><td>no</td><td>n.d.</td><td>st.; tail tube protein (TTP)</td><td>[16,91,92]</td></tr></tbody></table>
<table><caption>Table 3. *Cont.*</caption><thead><tr><th colspan="5">ORF</th><th rowspan="2">nt</th><th rowspan="2">Stop</th><th rowspan="2">nt</th><th rowspan="2">Essential Gene</th><th rowspan="2">Length (aa)</th><th rowspan="2">MM (kDa)</th><th rowspan="2">Predicted TMM Segments</th><th colspan="2">Protein</th><th rowspan="2">Protein Function</th><th rowspan="2">References</th></tr><tr><th>ORF</th><th>RBS (mRNA)</th><th>Spacing</th><th>Start</th><th>nt</th><th>3D Structure (PDB or EMD)</th><th>3D Structure (PDB or EMD)</th></tr></thead><tbody><tr><td>17.1*</td><td>AGGAGG</td><td>10</td><td>AUG</td><td>10,484 k</td><td>UAA</td><td>11,279</td><td>no</td><td></td><td>264</td><td>28.2</td><td>no</td><td>n.d.</td><td>st.; tail tube protein; Cter FN3 motif</td><td>[91]</td></tr><tr><td>17.5</td><td>GAGG</td><td>12</td><td>AUG</td><td>11,363 l</td><td>UAA</td><td>11,884</td><td>n.d. m</td><td></td><td>173</td><td>20.2</td><td>no</td><td>n.d.</td><td>tail chaperone protein</td><td>[91]</td></tr><tr><td>17.5*</td><td>GAGG</td><td>12</td><td>AUG</td><td>11,363 l</td><td>UAG</td><td>12,255</td><td>n.d. m</td><td></td><td>297</td><td>34</td><td>no</td><td>n.d.</td><td>tail chaperone protein</td><td>[91]</td></tr><tr><td>18</td><td>AGGAGG</td><td>9</td><td>AUG</td><td>12,267</td><td>UGA</td><td>15,365</td><td>n.d. m</td><td></td><td>1032</td><td>110.9</td><td>4 n</td><td>n.d.</td><td>st.; tape measure protein (TMP)</td><td>[16]</td></tr><tr><td>19.1</td><td>GAGG</td><td>10</td><td>AUG</td><td>15,362</td><td>UAA</td><td>16,123</td><td>n.d. m</td><td></td><td>253</td><td>28.6</td><td>no</td><td>2X8K (Xtal<sup>b</sup>)</td><td>st.; distal tail protein (Dit)</td><td>[93,94]</td></tr><tr><td>21</td><td>AAGA<underline>GGUGA</underline></td><td>10</td><td>UUG</td><td>16,137</td><td>UAA</td><td>19,463</td><td>n.d. m</td><td></td><td>1108</td><td>123.6</td><td>no</td><td>n.d.</td><td>st.; tail tip protein; Tal; anti-receptor protein</td><td>[93,95]</td></tr><tr><td>22</td><td>AAGGAGG</td><td>9</td><td>AUG</td><td>19,476</td><td>UAA</td><td>19,916</td><td>n.d.</td><td></td><td>146</td><td>16.7</td><td>no</td><td>2XC8 (Xtal<sup>b</sup>)</td><td>putative tail protein</td><td>[96]</td></tr><tr><td>23</td><td>AGGAGGU</td><td>10</td><td>AUG</td><td>19,932</td><td>UGA</td><td>20,096</td><td>n.d.</td><td></td><td>54</td><td>6.1</td><td>no</td><td>n.d.</td><td>n.d.</td><td>n.d.</td></tr><tr><td>23.1</td><td>GGAG</td><td>9</td><td>AUG</td><td>20,089</td><td>UAA</td><td>20,244</td><td>n.d.</td><td></td><td>51</td><td>5.8</td><td>no</td><td>2XF7 (Xtal<sup>b</sup>)</td><td>n.d.</td><td>[97]</td></tr><tr><td>24</td><td>GgGGUG h</td><td>10 h</td><td>AUG</td><td>20,237</td><td>UAG</td><td>20,467</td><td>n.d.</td><td></td><td>76</td><td>8.4</td><td>no</td><td>n.d.</td><td>n.d.</td><td></td></tr><tr><td>24.1</td><td>AAAGGgGG</td><td>11</td><td>AUG</td><td>20,547</td><td>UAA</td><td>20,825</td><td>n.d.</td><td></td><td>92</td><td>10.6</td><td>1</td><td>n.d.</td><td>component of holin; cell lysis</td><td>[33,98]</td></tr><tr><td>24.1*</td><td>AGGAGGU</td><td>10</td><td>AUG</td><td>20,574</td><td>UAA</td><td>20,825</td><td>n.d.</td><td></td><td>83</td><td>9.5</td><td>1</td><td>n.d.</td><td>component of holin; cell lysis</td><td>[33]</td></tr><tr><td>25</td><td>AAGGAG</td><td>12</td><td>AUG</td><td>20,845</td><td>UAA</td><td>21,660</td><td>n.d.</td><td></td><td>271</td><td>29.9</td><td>no</td><td>n.d.</td><td>endolysin; cell lysis</td><td>[33,98]</td></tr><tr><td>26</td><td>AAAGGAG</td><td>8</td><td>AUG</td><td>21,662</td><td>UAA</td><td>21,910</td><td>n.d.</td><td></td><td>82</td><td>9.4</td><td>2</td><td>n.d.</td><td>component of holin; cell lysis</td><td>[33,98]</td></tr><tr><td>26*</td><td>AAAGGAG o</td><td>14 o</td><td>AUG</td><td>21,668</td><td>UAA</td><td>21,910</td><td>n.d.</td><td></td><td>80</td><td>9.1</td><td>2</td><td>n.d.</td><td>component of holin; cell lysis</td><td>[33]</td></tr><tr><td>26.1</td><td>AAGGgGG o</td><td>10 o</td><td>AUG</td><td>22,009</td><td>UAG</td><td>22,152</td><td>n.d.</td><td></td><td>47</td><td>5.8</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>27</td><td>AAGGAGG</td><td>12</td><td>UUG</td><td>22,277</td><td>UAA</td><td>22,831</td><td>n.d.</td><td></td><td>184</td><td>20.8</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>28</td><td>GGAGG</td><td>9</td><td>AUG</td><td>22,834</td><td>UGA</td><td>23,121</td><td>n.d.</td><td></td><td>95</td><td>10.8</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>29</td><td>AGGAGG</td><td>13</td><td>AUG</td><td>23,069</td><td>UGA</td><td>23,371</td><td>n.d.</td><td></td><td>100</td><td>12</td><td>no</td><td>n.d.</td><td>putative DNA binding protein</td><td></td></tr><tr><td>29.1</td><td>AGGgGG</td><td>9</td><td>GUG</td><td>23,358</td><td>UGA</td><td>23,675</td><td>n.d.</td><td></td><td>105</td><td>12.3</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>30</td><td>AGGgGG</td><td>10</td><td>AUG</td><td>23,675</td><td>UAA</td><td>23,854</td><td>n.d.</td><td></td><td>59</td><td>7.2</td><td>1</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>30.1</td><td>AGGgGG</td><td>9</td><td>AUG</td><td>23,859</td><td>UGA</td><td>24,029</td><td>n.d.</td><td></td><td>56</td><td>6.4</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>31</td><td>AacGGAGGU</td><td>12</td><td>AUG</td><td>24,209</td><td>UAA</td><td>24,493</td><td>n.d.</td><td></td><td>94</td><td>11</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>31.1</td><td>GAGG</td><td>12</td><td>AUG</td><td>24,589</td><td>UAA</td><td>24,951</td><td>n.d.</td><td></td><td>120</td><td>12.9</td><td>3</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>31.2</td><td>GgGGUG h</td><td>10 h</td><td>AUG</td><td>24,964</td><td>UGA</td><td>25,281</td><td>no P</td><td></td><td>105</td><td>11.5</td><td>2</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>32</td><td>GGAGGUG</td><td>8</td><td>AUG</td><td>25,278</td><td>UAA</td><td>27,788</td><td>no P</td><td></td><td>836</td><td>96.3</td><td>no</td><td>n.d.</td><td>putative ATP-binding protein</td><td></td></tr><tr><td>32.5</td><td>GGAGGUG</td><td>11</td><td>UUG</td><td>28,039</td><td>UAA</td><td>28,209</td><td>n.d.</td><td></td><td>56</td><td>6.7</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>33</td><td>AAAaGgGGU</td><td>11</td><td>AUG</td><td>28,226</td><td>UAA</td><td>29,995</td><td>no</td><td></td><td>589</td><td>64.9</td><td>no</td><td>n.d.</td><td>putative bacteria surface binding protein</td><td></td></tr><tr><td>33.1</td><td>AacGGAGG</td><td>9</td><td>AUG</td><td>30,011</td><td>UAA</td><td>30,229</td><td>n.d.</td><td></td><td>72</td><td>8.4</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>34</td><td>AGGgGG</td><td>9</td><td>AUG</td><td>30,364</td><td>UAG</td><td>30,522</td><td>n.d.</td><td></td><td>52</td><td>6.3</td><td>no</td><td>n.d.</td><td>putative transcriptional repressor</td><td></td></tr><tr><td>34.1</td><td>AGGAGG</td><td>9</td><td>AUG</td><td>30,534</td><td>UGA</td><td>31,469</td><td>no</td><td></td><td>311</td><td>35.9</td><td>no</td><td>n.d.</td><td>5'-3' exonuclease</td><td>[13,99]</td></tr><tr><td>34.2</td><td>GGAGG</td><td>12</td><td>AUG</td><td>31,466</td><td>UAA</td><td>31,639</td><td>n.d.</td><td></td><td>57</td><td>6.7</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>34.3</td><td>AAGGAGG</td><td>11</td><td>AUG</td><td>31,641</td><td>UAA</td><td>31,895</td><td>n.d.</td><td></td><td>84</td><td>9.8</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>34.4</td><td>GGAGG</td><td>11</td><td>AUG</td><td>31,897</td><td>UAA</td><td>32,187</td><td>n.d.</td><td></td><td>96</td><td>11.1</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>35</td><td>GGAGGU</td><td>11</td><td>AUG</td><td>32,177</td><td>UAG</td><td>33,040</td><td>yes</td><td></td><td>287</td><td>32</td><td>no</td><td>n.d.</td><td>recT-like recombinase</td><td>[13,100]</td></tr><tr><td>36</td><td>AGGAGGgGA</td><td>10</td><td>AUG</td><td>33,033</td><td>UAA</td><td>33,512</td><td>no</td><td></td><td>159</td><td>17.1</td><td>no</td><td>n.d.</td><td>SSB</td><td>[101]</td></tr><tr><td>36.1</td><td>AAaGgGGUGA h</td><td>10 h</td><td>AUG</td><td>33,537</td><td>UGA</td><td>34,028</td><td>n.d.</td><td></td><td>163</td><td>18.9</td><td>no</td><td>n.d.</td><td>putative HNH endonuclease</td><td></td></tr><tr><td>37</td><td>GGAGG</td><td>10</td><td>AUG</td><td>34,032</td><td>UGA</td><td>34,406</td><td>n.d.</td><td></td><td>124</td><td>14.3</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr></tbody></table>
<table><caption>Table 3. *Cont.*</caption><thead><tr><th colspan="5">ORF</th><th colspan="5">Protein</th></tr><tr><th>ORF</th><th>RBS (mRNA)</th><th>Spacing</th><th>Start</th><th>nt</th><th>Stop</th><th>nt</th><th>Essential Gene</th><th>Length (aa)</th><th>MM (kDa)</th><th>Predicted TMM Segments</th><th>3D Structure (PDB or EMD)</th><th>Protein Function</th><th>References</th></tr></thead><tbody><tr><td>37.1</td><td>GGAGG</td><td>11</td><td>AUG</td><td>34,399</td><td>UGA</td><td>34,992</td><td>n.d.</td><td>197</td><td>22.3</td><td>no</td><td>n.d.</td><td>putative poly-gamma-glutamate hydrolase</td><td></td></tr><tr><td>37.2</td><td>GAGG</td><td>14</td><td>AUG</td><td>34,989</td><td>UGA</td><td>35,243</td><td>n.d.</td><td>84</td><td>9.7</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>37.3</td><td>AGGAGG</td><td>11</td><td>AUG</td><td>35,400</td><td>UAG</td><td>35,573</td><td>n.d.</td><td>57</td><td>6.7</td><td>no</td><td>n.d.</td><td>putative DNA binding protein</td><td></td></tr><tr><td>38</td><td>AAGGAGG</td><td>13</td><td>AUG</td><td>35,580</td><td>UGA</td><td>36,350</td><td>yes</td><td>256</td><td>30</td><td>no</td><td>n.d.</td><td>SPP1 origin binding protein and replication re-start (PriA-like)</td><td>[56,102]</td></tr><tr><td>39</td><td>AGGAGG</td><td>9</td><td>AUG</td><td>36,347</td><td>UGA</td><td>36,727</td><td>yes</td><td>126</td><td>14.6</td><td>no</td><td>1NO1 (Xtal<sup>b</sup>)</td><td>gp40 helicase loader</td><td>[56,103,104]</td></tr><tr><td>40</td><td>GAGG</td><td>11</td><td>AUG</td><td>36,724</td><td>UAA</td><td>38,052</td><td>yes</td><td>442</td><td>49.7</td><td>no</td><td>3BGW (Xtal<sup>b</sup>)</td><td>replicative DNA helicase; binds host DnaG and DnaX</td><td>[56,103,105–109]</td></tr><tr><td>41</td><td>AAAGGgGG</td><td>10</td><td>AUG</td><td>38,069</td><td>UGA</td><td>38,569</td><td>n.d.</td><td>166</td><td>19.1</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>42</td><td>AAAGGAG</td><td>11</td><td>AUG</td><td>38,566</td><td>UAA</td><td>38,961</td><td>no</td><td>131</td><td>16</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>42.1</td><td>GGAGG</td><td>9</td><td>AUG</td><td>38,964</td><td>UGA</td><td>39,134</td><td>no</td><td>56</td><td>6.5</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>42.2</td><td>GgGG<sup>c, q</sup></td><td>12<sup>c, q</sup></td><td>AUG</td><td>39,131<sup>q</sup></td><td>UAA</td><td>39,427</td><td>no</td><td>98<sup>q</sup></td><td>10.7<sup>q</sup></td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>43</td><td>GGAGG</td><td>10</td><td>GUG</td><td>39,431</td><td>UGA</td><td>39,784</td><td>no</td><td>117</td><td>14.2</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>44</td><td>AAGGAG</td><td>11</td><td>AUG</td><td>39,777</td><td>UAA</td><td>40,487</td><td>no</td><td>236</td><td>27.5</td><td>no</td><td>n.d.</td><td>Holliday junction resolvase</td><td>[13,31]</td></tr><tr><td>46</td><td>GAGG</td><td>12</td><td>AUG</td><td>40,596</td><td>UAA</td><td>40,898</td><td>no</td><td>100</td><td>11.5</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>46.1</td><td>AGGAGG</td><td>9</td><td>AUG</td><td>40,898</td><td>UAA</td><td>41,209</td><td>no</td><td>103</td><td>11.7</td><td>3</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>47</td><td>AGGgGG</td><td>9</td><td>AUG</td><td>41,304</td><td>UGA</td><td>41,663</td><td>no</td><td>119</td><td>13.7</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>48</td><td>GGAG</td><td>13</td><td>AUG</td><td>41,645</td><td>UGA</td><td>41,995</td><td>no</td><td>116</td><td>13.2</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>49</td><td>GGAGG</td><td>9</td><td>GUG</td><td>42,075</td><td>UGA</td><td>42,248</td><td>no</td><td>57</td><td>6.5</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>50</td><td>AAGGAGG</td><td>9</td><td>GUG</td><td>42,245</td><td>UAA</td><td>42,418</td><td>no</td><td>57</td><td>6.8</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>50.1</td><td>AAAGGAGG</td><td>9</td><td>GUG</td><td>42,434</td><td>UGA</td><td>42,616</td><td>no</td><td>60</td><td>6.7</td><td>2</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>51</td><td>AAGGAGG</td><td>9</td><td>AUG</td><td>42,613</td><td>UAA</td><td>43,014</td><td>no</td><td>133</td><td>14.7</td><td>1</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>51.1</td><td>AAGGAG</td><td>9</td><td>AUG</td><td>43,027</td><td>UGA</td><td>43,182</td><td>no</td><td>51</td><td>6.1</td><td>1</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>52</td><td>AAAGGAG</td><td>10</td><td>AUG</td><td>43,179</td><td>UGA</td><td>43,421</td><td>no</td><td>80</td><td>9.5</td><td>no</td><td>n.d.</td><td>unknown</td><td></td></tr><tr><td>53</td><td>AAAGGAG</td><td>10</td><td>AUG</td><td>43,405</td><td>UGA</td><td>43,611</td><td>no</td><td>68</td><td>7.4</td><td>1</td><td>n.d.</td><td>unknown</td><td></td></tr></tbody></table>

<footnote>a structure of gp1 from the SPP1-related phage SF6 (83% amino acid sequence identity with SPP1 gp1); b X-ray crystallography structure PDB access code; c No RBS identified according to the criteria defined in Material and Methods (see Section 2). In case of ORF 42.2 a G-rich sequence is identified as a potential site for ribosome binding; d the gene 2 beginning and the resulting length of gp2 is based exclusively on the position of the initiation codon assigned during annotation as no RBS was identified (see Section 3.2 for details) e gp2 is not a membrane protein according to presently available biochemical data; f Structures of the gp2 nuclease domain (residues 232 to 422 of the gp2 amino acid sequence); g cryo-electron microscopy structure EMD access code; h The RBS sequence is compatible with different spacings relative to the ORF initiation codon. i NMR structure PDB access code; j the gene 16 beginning and the resulting length of gp16 is based on the position of the initiation codon assigned during annotation, as no RBS was identified, and on the amino terminus sequencing of gp16 (see Section 3.2 for details); k genes 17.1 and 17.1* have the same 5' sequence because the product of 17.1* results from a +1 frameshift at the end of their common reading frame [91] (see Section 3.2 for details); l genes 17.5 and 17.5* have the same 5' sequence as the product of 17.5* results from a putative -1 frameshift within their common reading frame [91] (see Section 3.2 for details); m the essential nature of SPP1 genes coding for the tail chaperones, TMP, Dit and Tal proteins was not demonstrated experimentally but their functional homologs are essential in the phage systems presently characterized; n the SPP1 TMP gp18 features four predicted transmembrane segments but it is an anticipated component of the phage particle occupying the tail tube internal space [16] (see Section 3.8 for details); o the same RBS is used for translation at the initiation codons of genes 26 and 26* [33] (see Section 3.6 for details); p the sequence inversion in SPP1invmir (Figure 1) disrupts ORF 31.2 and renders gene 32 promoterless indicating that they are non-essential; q poor RBS sequence which might not ensure putative ORF 42.2 translation.</footnote>
AUG is the translation initiation codon of 72 ORFs (90%) of the SPP1 genome. It is found in all genes coding proteins of known function with the exception of the tail spike gene 21 that starts with UUG (Table 3). ORFs 27 and 32.5 also initiate with UUG while five other ORFs start with GUG (6%). These percentages differ somehow from the ones found for *B. subtilis* ORFs but AUG remains the initiation codon most frequently used (78%) by the host bacterium [110]. A significantly higher frequency of UGA termination codons (39%) is found in the SPP1 genome when compared to *B. subtilis* (24%). This is accounted by a reduction of UAA codon usage from 62% in *B. subtilis* to 50% in SPP1 (Table 4). UAG remains the rarest codon in the two genomes. The overrepresentation of UGA in SPP1 is surprising because there is an abnormally high 6% read-through of this codon in *B. subtilis* that extends the length of polypeptide chains [111]. However, the frequency of UGA is reduced to 29% when considering the SPP1 genes of known function that play key functions in DNA replication, viral particle assembly, and cell lysis (the "core" genome). In case of gene 6, which encodes the essential capsid portal protein, correct polypeptide chain termination is ensured by an arrangement of stop codons in tandem (UGAUAА). Furthermore, all genes coding proteins required in high amounts for phage particle assembly feature a UAA stop codon (hundreds of copies of the (pro)capsid proteins gp11, gp12, gp13 as well as of the major tail tube protein (TTP (note that TTPs are also designated major tail proteins (MTPs)) gp17.1 are used for assembly of one viral particle (see Section 3.8)). The global codon usage frequency shows also some differences between SPP1 and *B. subtilis* (frequency differences above 10% are highlighted in bold and rare codons (<10% frequency) are underlined in Table 4).

There are two documented cases of programmed translational frameshifts in SPP1. These recoding events result from slippage of ribosomes into a different coding frame during translation of mRNA. The frequency of the frameshift dictates a constant ratio of the two proteins synthesized which have an identical amino terminus but a different carboxyl region sequence. The two TTPs gp17.1 and gp17.1* share an identical sequence, but the gp17.1* carboxyl terminus is extended by 87 additional amino acids [91]. The +1 translational frameshift results of ribosomes pausing at CCC rare proline codons and their shift to the overlapping frequent CCU proline codon in the CCCUAA sequence which codes also for the gene 17.1 termination codon [91]. The amount of gp17.1/gp17.1* synthesized is compatible with their 3:1 ratio found in SPP1 tails. The other programmed frameshift was identified by bioinformatics in genes 17.5 and 17.5* which encode functional analogs of phage lambda gpG/gpGT [91,112]. These are chaperones of tail tube assembly [113,114]. In SPP1, the putative -1 frameshift occurs at a UUUUUUC heptanucleotide slippage sequence within gene 17.5 that leads some ribosomes to change coding frame yielding gp17.5* [91]. Gp17.5* has the first 112 amino acid sequence identical to gp17.5.

Genes with two initiation codons preceded by canonical RBSs for the same coding frame and spaced by a few codons were identified in the two holin genes (genes 24.1/24.1* and 26/26*) [33]; see Section 3.6).
<table><caption>**Table 4.** Codon usage bias in SPP1 and *B. subtilis*. The total number of codons (No) used in the complete set of ORFs of SPP1 (this work) or *B. subtilis* [47,48] and the fraction of each codon used to code a specific amino acid are listed. When the fraction value differs by more than 0.1 this variation of codon usage frequency is highlighted in bold. Rare codons (fraction of usage below 0.1) are underlined.</caption><thead><tr><th>Amino Acid</th><th>Codon</th><th>Fraction</th><th>No</th><th>Amino Acid</th><th>Codon</th><th>Fraction</th><th>No</th><th>Amino Acid</th><th>Codon</th><th>Fraction</th><th>No</th></tr></thead><tbody><tr><td rowspan="2">Ala</td><td>GCG</td><td>SPP1<br>B.s</td><td>0.26<br>24,574</td><td rowspan="2">Gly</td><td>GGA</td><td>SPP1<br>B.s</td><td>0.31<br>26,381</td><td rowspan="2">Pro</td><td>CCU</td><td>SPP1<br>B.s</td><td>0.22<br>12,824</td></tr><tr><td>GCA</td><td>SPP1<br>B.s</td><td>0.27<br>272</td><td>GGU</td><td>SPP1<br>B.s</td><td>0.24<br>245</td><td>Pro</td><td>CCC</td><td>SPP1<br>B.s</td><td>0.08<br>42</td></tr><tr><td rowspan="2">Ala</td><td>GCU</td><td>SPP1<br>B.s</td><td>0.27<br>279</td><td rowspan="2">Gly</td><td>GGC</td><td>SPP1<br>B.s</td><td>0.25<br>262</td><td rowspan="2">Ser</td><td>AGU</td><td>SPP1<br>B.s</td><td>0.16<br>123</td></tr><tr><td>GCC</td><td>SPP1<br>B.s</td><td>0.20<br>201</td><td>CAU</td><td>SPP1<br>B.s</td><td>0.34<br>131</td><td>AGC</td><td>SPP1<br>B.s</td><td>0.11<br>8096</td></tr><tr><td rowspan="2">Arg</td><td>AGG</td><td>SPP1<br>B.s</td><td>0.18<br>130</td><td rowspan="2">His</td><td>CAC</td><td>SPP1<br>B.s</td><td>0.47<br>119</td><td rowspan="2">Ser</td><td>UCG</td><td>SPP1<br>B.s</td><td>0.10<br>79</td></tr><tr><td>AGA</td><td>SPP1<br>B.s</td><td>0.10<br>4788</td><td>AUA</td><td>SPP1<br>B.s</td><td>0.27<br>261</td><td>UCG</td><td>SPP1<br>B.s</td><td>0.10<br>7717</td></tr><tr><td rowspan="2">Arg</td><td>CGG</td><td>SPP1<br>B.s</td><td>0.12<br>83</td><td rowspan="2">Ile</td><td>AUU</td><td>SPP1<br>B.s</td><td>0.13<br>335</td><td rowspan="2">Ser</td><td>UCU</td><td>SPP1<br>B.s</td><td>0.17<br>132</td></tr><tr><td>CGA</td><td>SPP1<br>B.s</td><td>0.15<br>7329</td><td>CUA</td><td>SPP1<br>B.s</td><td>0.24<br>11,517</td><td>UCA</td><td>SPP1<br>B.s</td><td>0.20<br>15,615</td></tr><tr><td rowspan="2">Arg</td><td>CGU</td><td>SPP1<br>B.s</td><td>0.18<br>125</td><td rowspan="2">Ile</td><td>AUC</td><td>SPP1<br>B.s</td><td>0.38<br>365</td><td rowspan="2">Ser</td><td>UCC</td><td>SPP1<br>B.s</td><td>0.11<br>83</td></tr><tr><td>CGC</td><td>SPP1<br>B.s</td><td>0.18<br>8755</td><td>CUG</td><td>SPP1<br>B.s</td><td>0.37<br>32,872</td><td>ACG</td><td>SPP1<br>B.s</td><td>0.13<br>9757</td></tr><tr><td rowspan="2">Arg</td><td>CGU</td><td>SPP1<br>B.s</td><td>0.18<br>125</td><td rowspan="2">Leu (s)</td><td>UUG</td><td>SPP1<br>B.s</td><td>0.24<br>246</td><td rowspan="2">Thr</td><td>ACG</td><td>SPP1<br>B.s</td><td>0.26<br>211</td></tr><tr><td>CGC</td><td>SPP1<br>B.s</td><td>0.14<br>99</td><td>CUA</td><td>SPP1<br>B.s</td><td>0.16<br>18,745</td><td>ACA</td><td>SPP1<br>B.s</td><td>0.27<br>17,693</td></tr><tr><td rowspan="2">Asn</td><td>AAU</td><td>SPP1<br>B.s</td><td>0.45<br>344</td><td rowspan="2">Leu</td><td>CUU</td><td>SPP1<br>B.s</td><td>0.24<br>251</td><td rowspan="2">Thr</td><td>ACU</td><td>SPP1<br>B.s</td><td>0.47<br>379</td></tr><tr><td>AAC</td><td>SPP1<br>B.s</td><td>0.55<br>425</td><td>CUA</td><td>SPP1<br>B.s</td><td>0.13<br>139</td><td>ACC</td><td>SPP1<br>B.s</td><td>0.16<br>10,620</td></tr><tr><td rowspan="2">Asn</td><td>GAU</td><td>SPP1<br>B.s</td><td>0.52<br>442</td><td rowspan="2">Leu</td><td>CUU</td><td>SPP1<br>B.s</td><td>0.20<br>203</td><td rowspan="2">Thr</td><td>UGG</td><td>SPP1<br>B.s</td><td>0.41<br>27,117</td></tr><tr><td>GAC</td><td>SPP1<br>B.s</td><td>0.48<br>415</td><td>CUC</td><td>SPP1<br>B.s</td><td>0.24<br>28,226</td><td>ACU</td><td>SPP1<br>B.s</td><td>0.15<br>120</td></tr><tr><td rowspan="2">Asp</td><td>GAU</td><td>SPP1<br>B.s</td><td>0.52<br>442</td><td rowspan="2">Leu</td><td>CUU</td><td>SPP1<br>B.s</td><td>0.20<br>203</td><td rowspan="2">Trp</td><td>UGG</td><td>SPP1<br>B.s</td><td>1.00<br>190</td></tr><tr><td>GAC</td><td>SPP1<br>B.s</td><td>0.64<br>40,291</td><td>CUC</td><td>SPP1<br>B.s</td><td>0.11<br>13,232</td><td>UAC</td><td>SPP1<br>B.s</td><td>0.16<br>10,497</td></tr><tr><td rowspan="2">Cys</td><td>UGU</td><td>SPP1<br>B.s</td><td>0.54<br>50</td><td rowspan="2">Lys</td><td>AAG</td><td>SPP1<br>B.s</td><td>0.38<br>460</td><td rowspan="2">Tyr</td><td>UAU</td><td>SPP1<br>B.s</td><td>0.50<br>276</td></tr><tr><td>UGC</td><td>SPP1<br>B.s</td><td>0.45<br>4429</td><td>UAC</td><td>SPP1<br>B.s</td><td>0.30<br>25,647</td><td>UAC</td><td>SPP1<br>B.s</td><td>0.50<br>278</td></tr><tr><td rowspan="2">Cys</td><td>UGC</td><td>SPP1<br>B.s</td><td>0.46<br>42</td><td rowspan="2">Lys</td><td>AAA</td><td>SPP1<br>B.s</td><td>0.62<br>760</td><td rowspan="2">Val (s)</td><td>GUG</td><td>SPP1<br>B.s</td><td>0.23<br>211</td></tr><tr><td>UGC</td><td>SPP1<br>B.s</td><td>0.55<br>5322</td><td>CUC</td><td>SPP1<br>B.s</td><td>0.70<br>60,072</td><td>UAC</td><td>SPP1<br>B.s</td><td>0.26<br>21,585</td></tr></tbody></table>
<table><caption>Table 4. *Cont.*</caption><thead><tr><th>Amino Acid</th><th>Codon</th><th>Fraction</th><th>No</th><th>Amino Acid</th><th>Codon</th><th>Fraction</th><th>No</th><th>Amino Acid</th><th>Codon</th><th>Fraction</th><th>No</th></tr></thead><tbody><tr><td rowspan="2">Gln</td><td rowspan="2">CAG</td><td>SPP1</td><td>0.38</td><td rowspan="2">Met (s)</td><td rowspan="2">AUG</td><td>SPP1</td><td>1.00</td><td>412</td><td rowspan="2">Val</td><td rowspan="2">GUA</td><td>SPP1</td><td>0.26</td><td>240</td></tr><tr><td>B.s</td><td>0.46</td><td>B.s</td><td>1.00</td><td>32,918</td><td>B.s</td><td>0.20</td><td>16,296</td></tr><tr><td rowspan="2">Gln</td><td rowspan="2">CAA</td><td>SPP1</td><td>0.62</td><td rowspan="2"><b>Phe</b></td><td rowspan="2"><b>UUU</b></td><td>SPP1</td><td><b>0.48</b></td><td>270</td><td rowspan="2">Val</td><td rowspan="2">GUU</td><td>SPP1</td><td>0.32</td><td>292</td></tr><tr><td>B.s</td><td>0.54</td><td>B.s</td><td><b>0.68</b></td><td>37,445</td><td>B.s</td><td>0.28</td><td>23,440</td></tr><tr><td rowspan="2">Glu</td><td rowspan="2">GAG</td><td>SPP1</td><td>0.39</td><td rowspan="2"><b>Phe</b></td><td rowspan="2"><b>UUC</b></td><td>SPP1</td><td><b>0.52</b></td><td>294</td><td rowspan="2">Val</td><td rowspan="2">GUC</td><td>SPP1</td><td>0.18</td><td>168</td></tr><tr><td>B.s</td><td>0.32</td><td>B.s</td><td><b>0.32</b></td><td>17,253</td><td>B.s</td><td>0.26</td><td>21,143</td></tr><tr><td rowspan="2">Glu</td><td rowspan="2">GAA</td><td>SPP1</td><td>0.61</td><td rowspan="2">Pro</td><td rowspan="2">CCG</td><td>SPP1</td><td>0.44</td><td>219</td><td rowspan="2"><b>End</b></td><td rowspan="2"><b>UGA</b></td><td>SPP1</td><td><b>0.39</b></td><td>31</td></tr><tr><td>B.s</td><td>0.68</td><td>B.s</td><td>0.43</td><td>19,421</td><td>B.s</td><td><b>0.24</b></td><td>965</td></tr><tr><td rowspan="2">Gly</td><td rowspan="2">GGG</td><td>SPP1</td><td>0.21</td><td rowspan="2">Pro</td><td rowspan="2">CCA</td><td>SPP1</td><td>0.26</td><td>127</td><td rowspan="2">End</td><td rowspan="2">UAG</td><td>SPP1</td><td>0.11</td><td>9</td></tr><tr><td>B.s</td><td>0.16</td><td>B.s</td><td>0.19</td><td>8541</td><td>B.s</td><td>0.14</td><td>591</td></tr><tr><td colspan="5"></td><td colspan="2"><b>End</b></td><td><b>UAA</b></td><td><b>SPP1</b></td><td><b>0.50</b></td><td>40</td></tr><tr><td colspan="5"></td><td colspan="2"></td><td><b>B.s</b></td><td><b>0.62</b></td><td>2542</td><td></td></tr></tbody></table>
### 3.3. *Organization and Function of the SPP1 Genes*

The 80 ORFs annotated here have a compact arrangement leaving only 6% non-coding sequences in the SPP1 genome. The longest intergenic regions carry transcriptional regulation sequences and, one of them, oriR (Figure 1).

The core SPP1 genome is presently composed of 28 genes: (i) 13 essential genes for which conditional lethal mutations (suppressor sensitive mutations (*sus*) or temperature sensitive mutations (*ts*)) were obtained by chemical mutagenesis of the overall SPP1 genome [7] (genes 1, 2, 6, 11, 13, 15, 16, 17, 17.1, 35, 38, 39, 40; see Sections 3.5 and 3.8); (ii) 11 genes coding proteins that are anticipated to act in phage tail assembly (17.5, 17.5*, 18, 19.1, 21, 22; see Section 3.8) or in host cell lysis (24.1, 24.1*, 25, 26 and 26*; see Section 3.6); and (iii) 4 non-essential genes whose inactivation is detrimental for SPP1 multiplication (7, 34.1, 36, and 44; see Sections 3.5 and 3.8). Genes 12 and 17.1* are excluded from this group because their inactivation has no detectable effect in SPP1 fitness in laboratory conditions (see Section 3.8). The core genes identified include most of the minimal genetic set necessary for lytic phage multiplication ensuring viral DNA replication, viral particle assembly, and host lysis. A present omission in SPP1 research is dissection of the genetic circuitry that controls genome transcription and, in particular, late genes expression. Several putative DNA-binding proteins encoded by early genes are good candidates to be transcriptional regulators (Table 3; see Section 3.4). Bioinformatics allowed us to assign putative function and/or biochemical activity to genes 16.1, 29, 32, 33, 34, 36.1, 37.1, and 37.3 (Figure 1; Table 3; see Sections 3.4 and 3.8). In total, 38 genes were functionally annotated (48% of SPP1 ORFs). They cluster in segments within four early transcriptional units, coding DNA replication and early lysis proteins together with uncharacterized polypeptides, and in one large genome region transcribed late. In this late region the synteny of genes encoding structural proteins of the SPP1 virion and late lysis genes is conserved when compared to the genomes of other *Siphoviridae* (Figure 1).

The function of other SPP1 ORF products is poorly understood. Genome deletions or inversions showed that ORFs 3–5, 7, 12, 14, 31.2, 32, 33 and 42–53 are non-essential (Figure 1; [26,35,53,55,83,88]; this work; P.T., unpublished) which corresponds to 19% of the SPP1 genome. Among these, a function was assigned experimentally only to the products of genes 7, 12 (see Section 3.8) and 44 (see Section 3.5) while bioinformatics revealed a putative role for gp32 and gp33 (see Section 3.4). Within the SPP1 genome regions not tested by deletion analysis, the (putative) function of the products from ORFs 8, 9, 10, 23, 23.1, 24, 26.1, 27, 28, 29.1, 30, 30.1, 31, 31.1, 32.5, 33.1, 34.2, 34.3, 34.4, 37, 37.2, and 41 remains unknown. Most of these 22 ORFs are probably non-essential for SPP1 multiplication. In total, 42 ORFs (52% of the total SPP1 ORFs) have no assigned function. Protein sequence (BLASTp) and structural (HHPred) homology searches revealed 18 orphan ORFs among those (8, 10, 14, 23, 23.1, 24, 28, 30.1, 32.5, 33.1, 34.3, 34.4, 37, 37.2, 42.1, 46.1, 51.1, and 52) whose products are unrelated to protein sequences in databases. The other SPP1 ORF protein products are homologous in their vast majority to proteins with no known function from *Bacillus* spp. or from their phages, suggesting that they belong to a mobile gene pool of bacilli and their viruses. Most of the less characterized SPP1 ORFs are likely additions to the lytic phage core genome resulting from insertions of DNA that add “more on to it” (“morons” [115–117]). Morons designated originally complete transcriptional units added to phage genomes [115] but include presently also gene insertions that are unique or found in a limited set of genomes. They frequently encode beneficial features for phage adaptation to its host and environment. We consider SPP1 morons the ORFs 3–5, 8–10, 12, 14, 23–24, 26.1–28, 29.1–31.2, 32.5, 33.1, 34.2–34.4, 37, 37.2, 41–43, and 46–53.

The SPP1 ORFs length is highly variable, encoding proteins with an average length of 179 ± 197 amino acids and an average molecular mass of 20.3 ± 21.9 kDa. The tail tape measure protein and in some cases the tail tip protein, both with more than 1000 residues in SPP1, is (are) the longest protein(s) encoded by phages with long tails. They are landmarks to map the phage tail genes which, when combined with the synteny of genes coding virion assembly proteins, provide a first approximation to the genomic organization of late genes. The SPP1 core genes (most coding proteins >100 amino acids-long) tend to be longer than moron genes (most coding proteins <100 amino acids-long).
A BLASTn of the complete SPP1 genome showed extensive nucleotide sequence homology only to phages of the SPP1 group rho15, SF6, 41c [9,118], and the recently identified Lurz phage series [35]; P.T. unpublished) defining the genetic basis for a SPP1-like *genus*. Strong DNA sequence homology hits to other phages and to *Bacillus* spp. genomes (probably to prophages, defective phages or other elements of the bacilli mobile genome) were limited to a few individual genes (5, 9, 10, 17.1*, 25, 33, 48, 49, 50.1, and 51). An interesting case is *B. subtilis* (natto) phage PM1 [119] that has four blocks of DNA homology to four SPP1 ORFs/genes (9, 25, 33 and 51) which are separated by unrelated sequences. Multiple horizontal gene transfer events thus occur within the genetic mobile landscape of *Bacillus* spp. and its phages.

DNA nucleotide homology provides a sensitive criterion to detect recent genetic exchanges, while protein amino acid homology, found for a much larger set of SPP1 gene products, assesses far-reaching functional relationships with proteins from other phages. This is particularly well illustrated by bacteriophage GBK2 that infects the termophilic bacterium *Geobacillus kaustophilus*. Its genome shows no nucleotide sequence homology to SPP1 which infects a mesophilic bacillus. However, a large region encodes proteins homologous to SPP1 gp26 (holin), gp29, gp31, gp31.1, gp32, known DNA replication proteins (gp34.1, gp35, gp36, gp39, gp40), gp42, gp42.2 and gp43 [120] (see Sections 3.4 and 3.5). Their genes order is conserved in the two phage genomes, being spaced by genes of unrelated proteins. In contrast, phage particle assembly proteins of the two phages lack significant homology, apart from SPP1 proteins gp16.1 and gp17 which have GBK2 homologs, while the ensemble of their coding genes conserves the synteny found in siphoviruses [120] (this study). Thus, the genome of both phages is assembled from a common ancient genome module of early genes, traceable by protein homology, and of an evolutionarily distinct module coding the viral particle assembly proteins. The first module diversified by acquisition and eventual loss of different morons (single genes and also transcriptional units like the one controlled by PE4 in SPP1; Figure 1) while the viral particle assembly module acquired the genes 16.1–17 cluster, most probably from horizontal exchange into a conserved position of the module. The genomes of SPP1 and GBK2 are a compelling case of tailed phages evolution by combination of modules, horizontal gene transfer, and diversification by morons acquisition/loss [115–117,121,122].

In the following Sections we describe current knowledge on SPP1 genes and ORFs with particular attention to those that are less characterized. In-depth reviews on SPP1 biology [123], DNA replication [30], and viral particle assembly [124–126] are available.

### 3.4. The SPP1 Genes Set. I. Uncharacterized Early Genes

Putative moron genes are spread throughout the five early SPP1 transcriptional units (Figure 1).

The dispensable set of short ORFs 46 to 53, under the control of PE1, has no known or putative function deduced from bioinformatics. ORFs 46.1, 53, and the sequential ORFs 50.1, 51, 51.1 have predicted transmembrane segments (Table 3) indicating that these polypeptides insert most likely in the cytoplasmic membrane to achieve early roles in the host cell. Segments of DNA homologous to different combinations of ORFs 48 through 51 are found in different SPP1 phages or *Bacillus* spp. strains, suggesting common functions.

Protein sequence analysis and predicted structural homology led to functional assignments to several genes in transcriptional units controlled by promoters PE4 and PE5. PE4 drives expression of ORFs 32.5 and 33.1, whose function is unknown, together with the dispensable gene 33. Gp33 is a 589-long protein that shares 88% amino acid sequence identity with a protein of *B. subtilis* (natto) phage PM1. It is also highly homologous to a large number of proteins from *Bacillus* spp. annotated to have cell wall hydrolysis activity or as proteins of *Bacillus* phages with a predicted right-handed parallel β-helix repeat fold. HHPred extends this structural homology to adhesins binding to the cell wall surface as well as to phage tail spikes of *B. subtilis* phage phi29 and of enterobacteria phages like sf6, HK620, and LKA1. Gp33 could be the trace of a tail fiber used by an ancestor of SPP1, like the Ur-lambda tail fibers lost during laboratory evolution [127]. However, this does not explain why SPP1
maintained gene 33 expressed early during infection in an apparently functional form. We privilege the hypothesis that gp33 accumulates in the cytoplasm and its release upon lysis acts on the wall of the infected bacterium to facilitate lysis and/or on the envelope of other *B. subtilis* cells to facilitate their subsequent infection.

PE5 controls expression of 14 ORFs including genes 25, 26 and 26* involved in cell lysis (see Section 3.6). Most of the other ORFs encode short polypeptides of unknown function. Gp29 and its homologous hypothetical protein of *Geobacillus* phage GBK2 give strong hits in HHPred to DNA-binding proteins with a winged helix fold like transcriptional regulators and excionases. ORFs 30, 31.1, and 31.2 products have predicted transmembrane segments suggesting an early role in the cytoplasmic membrane during infection. Gp31.1 and gp31.2, which are encoded by contiguous ORFs, are homologous to numerous proteins of *Bacillus* spp. and their phages, likely sharing a widespread function. The non-essential gp32 (Table 3) is a basic 836-long protein with 75% amino acid sequence identity to a protein of unknown function from *Geobacillus* phage GBK2. Gp32 has also extensive similarity to ATP-binding proteins from *Bacillus* spp. and its phages. HHPred reveals structural homology of its carboxyl terminus, with high confidence scores, to the conjugation protein TrwB, to VirB4 of type IV secretion systems and to proteins of the FtsK family. These machines translocate DNA or proteins across the bacterial membrane. Since gp32 has no predicted transmembrane segments, it could act together with gp31.1 and/or gp31.2 to build a trans-membrane translocon of macromolecules. The synteny of genes coding for proteins homologous to gp31, gp31.1 and gp32 is found also in phage GBK2 [120] suggesting a conserved activity. The exact protein composition, functionality and role of such potential machine remain to be established.

The strongest early promoters, PE3 and PE2 [41], drive expression of DNA replication genes (see 3.5) that alternate with genes coding uncharacterized proteins (Figure 1). Bioinformatics allowed us to assign putative activities to four proteins of the latter group. Gp34 is a putative 6.3 kDa basic polypeptide with predicted structural homology to DNA binding proteins, whose strong hits are transcriptional repressors. Its small size advises, nevertheless, some caution on this functional assignment. Gp36.1 is highly homologous to HNH endonucleases of Gram-positive bacteria, and HHPred uncovered a relationship to the structure of *B. subtilis* phage SPO1 HNH homing endonuclease I-HmuI (PDB accession number 1U3E) [128]. Proteins containing the HNH motif carry out intron homing in phages T4, SPO1 and SP82 [129,130], and are essential for DNA packaging in a group of cos-phages [131,132]. In other phages they play dispensable functions like in *Staphylococcus aureus* phage 80α [133]. Gp37.3 is a small basic polypeptide whose amino acid sequence and structure prediction show relatedness to numerous bacterial and phage DNA binding proteins. SPP1 gp34, gp36.1, gp37.3 together with gp29, and possibly other early polypeptides that did not deliver robust bioinformatics hits yet, are strong candidates to participate in SPP1 DNA metabolism and in gene expression regulation processes calling for further research. Gene 37.1 that is found within this set of genes codes for an enzyme with a distinct function. Gp37.1 is highly homologous to poly-γ-glutamate (γ-PGA) hydrolases from phages of Gram-positive bacteria, mainly infecting *Bacillus* spp. like vB BsuM-Goe3, BSP10, Grass, BSNPO1, PM1, phiNIT1, PBS1, AR9 and others. HHPred provided a strong hit to the structure of the γ-PGA hydrolase PghP from phage phiNIT1 (PDB accession number 3A9L). The supposed role of the enzyme in phages is to degrade the γ-PGA polymer layer. This layer forms a shield at the surface of several microorganisms, mainly *Bacillus* spp., to protect them from environmental attacks as diverse as phagocytosis or phage infection [134,135]. Release of γ-PGA hydrolase from phage-infected bacteria would thus open the way for infection of new host bacteria protected by a γ-PGA layer. Experimental demonstration of γ-PGA hydrolase activity from *B. subtilis* phages phiNIT1 (protein PghP; [136]), BSP10 [137], as well as from prophages SPβ and prophage-like element 5 (YokZ and YndL, respectively) [135] support this hypothesis.

The large number of early ORFs not conserved among SPP1 and other tailed phages, whose majority has an unknown function, is a rich and diverse patrimony. We anticipate that this genetic pool is mostly dedicated to strategies engaged at the beginning of infection to take over the host cell [138]
(and references therein), possibly to mediate superinfection exclusion of the infected bacterium by other phages [139,140] (and references therein), and/or to accumulate molecules that will be released upon lysis to support subsequent infection of new host bacteria [134–137]. Most of the phage effectors involved are dispensable (e.g., the proteins encoded by transcriptional units under the control of PE1 and PE4 in SPP1) but their combination might provide a determinant fitness advantage, globally or in specific environmental settings, to SPP1-like phages.

### 3.5. *The SPP1 Genes Set. II. DNA Replication Early Genes*

The bacteriophage SPP1 DNA replication is a well-characterized process. SPP1 was reported to have two origins of replication, *oriR* and *oriL* (Figure 1), localized ~13 kbp apart in the circularized genome. The two sequences have a similar organization composed of direct repeats and an AT-rich region that acts as a DNA unwinding element (DUE) [102]. The SPP1 origin binding protein gp38 binds to the direct repeats of *oriL*, positioned within gene 38 (Figure 1), and of *oriR* with nanomolar affinity [102]. Replication initiation requires gp38, the helicase loader gp39 [103,104] (PDB accession number 1NO1), and the replicative helicase gp40 that belongs to the DnaB family [103,105,106] (PDB accession number 3BGW). After melting of the origin sequence and unwinding, the gp40 hub interacts with the host DnaG [107,108] and with the DnaX subunit of the clamp loader [109]. These interactions recruit the host replisome. The theta replication reaction was reconstituted in vitro with a supercoiled plasmid bearing *oriL* and purified SPP1 gp38, gp39 and gp40 together with *B. subtilis* DNA polymerases PolC and DnaE, the processivity clamp DnaN, the hetero-pentameric clamp loader complex (3xDnaX+YqeN+HolB), the primase DnaG, DNA gyrase, and a single-stranded DNA binding protein (SSB) (bacterial SsbA, or SPP1 gp36) [141]. Cellular SsbA can replace the non-essential SPP1-encoded SSB gp36 in SPP1 DNA replication. In contrast, gp36 does not support *B. subtilis* DNA replication and acts in vitro as an inhibitor of the cellular process [101].

After a few rounds of circular molecules DNA synthesis, SPP1 DNA replication switches to generate linear head-to-tail concatemers [28,29] which are the substrate for DNA encapsidation into viral particles. The switch was proposed to be triggered by stalling of the replication fork, probably when it collides with gp38 firmly bound to *oriR* [30,142]. Four SPP1 proteins are likely involved in the process of resuming DNA replication after stalling: the Holliday junction resolvase gp44 [31], the 5'→3' exonuclease gp34.1 [99], the ATP-independent single-strand annealing recombinase gp35 [100], and the SSB gp36 [101]. The current model is that they would act coordinately to process the stalled replication fork and make a double-strand break, followed by homology-directed recombination generating the template for DNA replication by a rolling circle mode (sigma-type) [30,31,142]. A DNA substrate assembled in vitro which mimics the sigma template was used to show that subsequent rolling circle replication is achieved by the same set of proteins that carry out theta-type DNA replication, with the exception of DNA gyrase that is not necessary [101]. The assays also highlighted that gp38 may act like the bacterial PriA protein, because no specific DNA region (*ori*) is needed for this reaction and DNA replication can start at any site after fork pausing [101]. Gp35 recombinase is the only SPP1 essential protein among the switch putative effectors. The lack of gp34.1 and gp36 is detrimental while deletion of the gp44-encoding gene has only a marginal effect in SPP1 viability [13]. Functionally related cellular proteins most likely (partially) compensate for the roles of those viral effectors. Homologs of the SPP1 DNA replication proteins are found in many phages, suggesting that such recombination-driven strategy is a common mechanism to produce concatemers.

The SPP1 replication proteins are encoded by two transcriptional units under the control of promoters PE2 and PE3 (Figure 1), as in phage GBK2 that shares a similar operon organization and gene synteny [120]. Genes 38, 39, and 40 whose products initiate DNA replication form a cluster, but the other DNA replication genes are spaced by ORFs coding proteins with putative DNA binding, host takeover, or unknown functions (see Section 3.4). The selective pressure to cluster genes encoding proteins that closely interact [143] appears stronger for the viral particle assembly genes module (see Section 3.8) than for the DNA replication genes module.
### 3.6. The SPP1 Genes Set. III. Lysis Early and Late Genes

Tailed phage lysis cassettes code typically for a holin that inserts in the cytoplasmic membrane to create holes and for an endolysin that hydrolyses the cell wall leading to efficient bacterial disruption [33,98,144]. Interestingly, The SPP1 genome codes two holin proteins that localize in the cytoplasmic membrane (gp24.1 and gp26) and an endolysin (gp25 or LysSPP1). They act together to lyse *B. subtilis* at the end of the infectious cycle allowing viral particles to escape from the host [33]. Strong nucleotide homology was only detected between gene 25 and the endolysin gene of phage PM1 [119] but proteins homologous to gp24.1, gp25 and gp26 are encoded by numerous phages of Gram-positive bacteria.

A lysis system with two holins (XhIA and XhIB) combined with an endolysin (XlyA) was also described for the defective *B. subtilis* phage PBSX [145]. XhIA and XhIB are homologous to SPP1 gp24.1 and gp26, respectively. In SPP1 their coding genes flank the endolysin gene 25, while the two holin genes precede the endolysin gene in the PBSX genome. Gp24.1 and XhIA feature a predicted transmembrane helix in their carboxyl terminus while gp26 and XhIB have the canonical organization of phage holins with two transmembrane segments. Each individual holin of these phages does not appear detrimental to *B. subtilis*, while their co-production leads to immediate cell death followed by lysis, showing that the two holins cooperate for disruption of the bacterial membrane [33,145].

Both SPP1 genes 24.1 and 26 have one in-frame internal initiation codon AUG preceded by a correctly spaced RBS (Table 3). Translation started at these AUG codes gp24.1* and gp26* which lack the first 9 and 2 amino terminus amino acids when compared to gp24.1 and gp26, respectively. Such dual start motifs leading to synthesis of two highly related proteins might be involved in lysis regulation as found for the phage lambda S protein whose longer polypeptide has an antiholin function that counteracts the shorter polypeptide holin activity [146]. In contrast to most phages, SPP1 genes 25 and 26 are transcribed early during infection from promoter PE5 localized within the gene 24.1 coding sequence [33]. It is likely that the production of gp24.1, which results from late transcription, defines the tempo for the three SPP1 lysis proteins to concur for disruption of the infected bacterium.

### 3.7. The SPP1 Genes Set. IV. Uncharacterized Late Genes

The SPP1 late genome region has an order of genes coding proteins involved in assembly of the viral particle that is conserved among numerous tailed phages. This arrangement is interrupted by a genome segment, bracketed between transcriptional signals, which codes for the origin of replication *oriR* and three moron genes (Figure 1). Gp8 and gp10 are unrelated to known proteins, while gp9 is highly homologous to a protein from *B. subtilis* (natto) phage PM1 [119] and to numerous hypothetical proteins encoded by *Bacillus* spp. The function of these SPP1 proteins and if they are necessary for phage amplification remain to be established.

The genome segment between the tail fiber gene 21 and the lysis proteins encodes also four proteins whose precise role is unknown. Gp22 has homology to a group of *Bacillus* spp. putative proteins and its structure (PDB access code 2XC8) reveals a fold similar to a domain of the tail receptor binding protein from lactococcal phage p2 [96] (see Section 3.8). The structure of the gp23.1 hexamer was also determined (PDB access code 2XF7) but provided no conclusive insight on protein function [97]. Gp23, gp23.1 and gp24 gave no strong homology hits to proteins in the databank.

The role of the dispensable short ORFs 3, 4, 5 and 14 that separate essential genes for SPP1 capsid assembly is not known. The products of ORFs 3 and 4 show homology to a hypothetical protein of *Bacillus* phages vB_BsuM-Goe3 and PM1, respectively, while gp5 has similarity to a large number of bacterial and phage uncharacterized proteins. ORF 14 is an orphan.

### 3.8. The SPP1 Genes Set. V. Viral Particle Assembly Late Genes

SPP1 devotes ~40% of its genome information to assembly of the SPP1 viral particle, a process that was extensively studied. DNA-filled capsids (nucleocapsids) and tails are built in two independent
assembly pathways. The two structures then join in a final reaction to yield the infectious particle (virion), like in all studied phages with a long tail.

Construction of the SPP1 nucleocapsid follows the similar assembly pathway used by viruses of the tailed phages-herpesviruses lineage [147,148] (and references therein). A spherically shaped icosahedral DNA-free procapsid with a diameter of ~55 nm is formed first. It is built by polymerization of the major capsid protein (MCP) gp13 that is chaperoned by the internal scaffolding protein gp11 [85–87]. One of the 12 vertexes of the procapsid icosahedron is a specialized structure defined by presence of a cyclical dodecamer of the portal protein gp6 ([20,60,75–81]; PDB access code 2JES). Assembly of the procapsid initiates most likely at the portal vertex by co-interaction between gp6, gp11 and gp13. In absence of gp6, gp11 and gp13 assemble procapsids of normal and of a smaller size, revealing that the portal ensures correct size determination of the procapsid [86]. Co-production of gp6, gp11 and gp13 is sufficient and necessary for assembly of procapsids functional for DNA packaging, showing that they are the minimal set of essential components to build these structures [86]. Gp6 makes a strong interaction with gp7 that targets this protein to the procapsid interior in a few number of copies [84]. This strategy was likely exploited in numerous lysogenic phages for targeting to their capsid interior a toxin fused to the carboxyl terminus of gp7-like proteins (also designated Mufs) [149]. The toxin was proposed to be subsequently delivered to the host at the beginning of infection [149,150]. SPP1 gp7 is a dispensable component of the virion. Phage assembly is not affected in its absence, but only ~25% of virions lacking gp7 are infectious, showing that it supports initiation of SPP1 infection [83].

Head-to-tail concatemers of the SPP1 genome synthesized during phage DNA replication (see Section 3.5) are the substrate for encapsidation into procapsids. The *pac* sequence in the SPP1 genome [151] (Figure 1) is specifically recognized and cleaved by the SPP1 terminase complex to initiate DNA packaging [70,152]. The terminase is composed of the small subunit (TerS) gp1 that binds to *pac* [57,63–66] (PDB access code 3ZQQ of gp1 from the SPP1-related phage SF6) and of the large subunit (TerL) gp2, a two-domain protein with ATPase and nuclease activities [57,67,70–72] (PDB access codes 2WBN and 2WC9 of the gp2 nuclease domain). Re-analysis of gene 1 sequence showed that gp1 is shorter (Table 1) and more similar to the TerS from other SPP1-related phages [118] than initially reported. The *pac* sequence overlaps the PL1/PL2 promoters (see Section 3.2), the RBS, and the sequence of gene 1 encoding the gp1 DNA-binding domain (Figure 1) [63]. Gp1 binds to the *pacL* and *pacR* regions of *pac* [63,68] and recruits gp2 for making a double-strand cleavage within *pacC*, a sequence that is flanked by *pacL* and *pacR* [70]. Gp2 has a non-specific nuclease activity [57,67,70–72] that is controlled by its appropriate positioning in the gp1-gp2-DNA complex [35]. Cleavage at *pac* is auto-regulated [69] leaving most *pac* sequences uncut in the concatemer used for packaging. This allows encapsidation of DNA molecules longer than the unit-genome length during processive headful packaging (see Section 1) [17,19]. The terminal redundancy of SPP1 encapsidated DNA molecules is essential for their re-circularization at the beginning of infection. The auto-regulated cleavage of *pac* was reproduced in a plasmid minimal system bearing genes 1, 2 and *pac*, showing that gp1 and gp2 are necessary and sufficient to carry out this reaction in *B. subtilis* [72].

The gp1-gp2-SPP1 DNA complex docks at the procapsid portal vertex to assemble the DNA packaging motor. DNA translocation powered by the ATPase activity of gp2 involves an intricate cross-talk between gp1, gp2, and gp6 [67,73,74,82]. The scaffolding protein gp11 leaves the capsid interior and the capsid lattice undergoes a major conformational change leading to its expansion, maximizing the space for phage DNA packing [85,86]. This process uncovers binding sites at the capsid surface for attachment of gp12 trimers that have a central collagen-like fold [88]. The DNA packaging reaction was characterized in vitro, either using extracts of infected cells [153] or purified linear DNA, gp1, gp2 and procapsids [154]. When a threshold amount of DNA is packaged inside the capsid, the portal protein senses the level of DNA headfilling and triggers the sequence-independent headful cleavage of the DNA concatemer that is most likely achieved by the gp2 nuclease domain [20,77,125,155]. Disassembly of the DNA packaging motor is coordinated with sequential binding of the head completion protein gp15 (PDB access code 2KBZ) to gp6, extending the portal channel, and of gp16 (PDB access code
2KCA) that closes the channel, preventing release of the packaged DNA [60,78,89]. The gp6-gp15-gp16 complex is named connector.

The SPP1 long non-contractile tail has a building plan similar to the long tails of *Siphoviridae* and *Myoviridae*. These structures are structurally related to a variety of cellular, tube-like, delivery devices used in bacterial warfare like phage tail-like bacteriocins (PTLBs), phage tail-like complexes that confer toxicity against eukaryotic cells, or type VI secretion systems [156,157].

The SPP1 tail features a host adsorption apparatus that binds selectively to the *B. subtilis* surface and promotes DNA transfer across the bacterial envelope. This ~31 nm-long structure binds selectively to glycosylated teichoic acids [21], which facilitates its subsequent strong interaction with the membrane protein YueB [23,158], committing the phage particle to infection. The known components of the adsorption apparatus are the tail fiber (or Tal) gp21 [93], which binds YueB [95], and the distal tail protein (Dit) gp19.1 (PDB access code 2X8K) [94]. The gp19.1 hexamer forms a complex with the trimeric gp21 amino terminus that closes the tail tube in the virion and opens at the beginning of infection [93]. Gp19.1 defines one end of the ~160 nm-long tail tube built by a helical array of the TTPs gp17.1/gp17.1* that are found at a ratio of ~3:1 in the tube [91]. Similar tubes can be formed exclusively by gp17.1 in vivo [91] and in vitro [92]. The TTPs form hexameric rings organized around the tape measure protein (TMP) gp18 [16]. The tail tube of sipho and myoviruses is tapered by (a) tail completion protein(s), which provide(s) the interface for tail attachment to the connector. In SPP1 it is the tail-to-head joining protein (THJP) gp17 (PDB access code 2LFP) that interacts with the capsid connector in the viral particle [32,90]. The connector-tail completion proteins structure is named head-to-tail interface [81].

SPP1 tail assembly is anticipated to follow the pathway of phages with long tails typified by T4 [159] and lambda [160]. In such case, the tail adsorption apparatus is formed first providing a platform for initiation of gp17.1/gp17.1* helical polymerization around the TMP gp18. During this reaction gp18 is probably pre-shielded by the chaperones gp17.5 and gp17.5*, as proposed for phage lambda [114]. When the tail tube reaches a defined length, determined by gp18, polymerization stops and tail completion proteins bind to the tube end. Gp17 that is exposed at this end joins the tail to the capsid, a reaction that was characterized in vitro [32].

The function of proteins engaged in assembly of the SPP1 nucleocapsid was defined. They have a large number of homologous proteins from phages or prophages of Gram-positive bacteria that play similar roles in viral particle assembly. The exception is the non-essential capsid accessory protein gp12, whose collagen-like motif is found only in some capsid associated proteins [88] and in tail fibers [161]. The order of genes coding proteins necessary for nucleocapsid assembly follows the conserved genome organization of siphoviruses [162-165] with the particularity that ORFs 3–5, 8–10 and 14 interrupt the usual gene order. The DNA packaging/portal module (coding gp1, gp2, (-, -,-), gp6, (gp7) spaced by ORFs 8–10, is followed by the icosahedral capsid building module (coding gp11, (gp12), gp13) which is separated by ORF 14 from the head completion proteins module (coding gp15, gp16). Note that in this protein list “-” identifies a protein encoded by a moron gene while proteins that have a known function but are non-essential for SPP1 multiplication are shown within brackets.

Gene 16 is followed by 16.1 and a set of genes whose synteny and protein products are conserved in the tail module of siphoviruses [16,91,162–167]. The function of SPP1 gp16.1 and its homologous proteins that are widespread among phages with long tails is yet unknown. Their coding gene is consistently found between the capsid connector and the tail proteins coding genes [166]. The tail module codes for gp17, gp17.1, (gp17.1*), gp17.5, gp17.5*, gp18, gp19.1, and gp21 (Figure 1) [16]. They all have extensive homology to proteins from Gram-positive bacteria and their phages. Gp17.1* and gp17.5*, which result from programmed translational frameshifts (see Section 3.2), feature two domains defined by similarity to different proteins. The amino terminus domains boundaries are roughly delimited by their region of identity with gp17.1 and gp17.5, respectively. The gp17.1* amino terminus shows robust similarity to TTP annotated proteins, while its non-essential carboxyl terminus is homologous to FN3 motifs of proteins from *Bacillus* spp. that are involved in binding to cell
surfaces [91]. Their exposure in phage structures suggests a role in adhesion to bacteria [91,168]. The amino and carboxyl domains of gp17.5* have similarity to two distinct sets of proteins with unknown function from Gram-positive bacteria or from their phages. The full-length gp21 Tal has numerous homologs. In addition, a subset of proteins has similarity only to its amino terminus (~400 residues), which is the conserved region of Tal proteins that closes the tail tube [93]. The gp21 long carboxyl terminus interacts with the SPP1 bacterial receptor [95]. This organization is consistent with the modular organization of Tal proteins whose carboxyl termini differs both in length and sequence. Such variation results, namely, from whether they carry or not a domain for interaction with the bacterium and of the selective pressure on this domain to generate diverse strategies to target a specific host [169,170] (and references therein). It remains to be established if proteins gp22, gp23, gp23.1 and/or gp24 (see Section 3.7) which are encoded by genes downstream of gene 21 play a role in tail assembly that is specific for the SPP1 system. Note that the three-dimensional structure of the gp22 monomer shows structural similarity to the shoulder of the tail receptor binding protein of lactococcal phage p2 [96] supporting that gp22 is a SPP1 tail component.

**4. Conclusions**

This reannotation of the SPP1 genome provides an actualized view of our understanding of this phage genetic patrimony and how it supports its multiplication. Phage DNA replication and assembly of the viral particle are particularly well studied, rendering SPP1 one of the forefront systems to understand the molecular basis of these processes in tailed bacteriophages. The available knowledge and tools are excellent assets to pursue research on both themes. Recent work provided also insight on the SPP1 lysis mechanism. The genes involved in these essential steps of the virus cycle occupy ~60% of the SPP1 genome. Much less is known about the function of other genes in spite of 50 years of SPP1 research. One or several of them code yet unidentified proteins that regulate SPP1 gene expression circuitry, a less studied aspect of SPP1 molecular biology. Apart from such regulators, a majority of uncharacterized genes likely code effectors that participate in host cell hijacking to optimize viral multiplication, that provide immunity to super-infection by other phages, or that support phage dissemination in the natural environment. These functions are frequently non-essential for phage survival, having subtle effects that render their study difficult. Transcriptomics, proteomics, and metabolomics combined with systematic knock-outs and sensitive phenotyping assays in different infection settings appears as a promising approach to deliver a complete functional map of SPP1. Those studies will provide insights on the wild side of SPP1 that remains to be explored.

Current knowledge on this phage system and its genetic landscape, that is distinct from other model phages, clearly recommend SPP1 as the reference virus for a new SPP1-like virus genus of the *Siphoviridae* family. The genus includes SPP1-related phages rho15, SF6, 41c [9,118], and the recently identified Lurz phage series [35] (P.T. unpublished) whose DNA sequence and genome organization are similar to SPP1. Phages GBK2 and PM1 have modules of the genome that are evolutionarily linked to SPP1 but other modules of essential genes code proteins without detectable similarity to SPP1 gene products. This mosaic genome organization brings to debate if such type of relatedness [122] has enough taxonomic value to include GBK2 and PM1 in the same taxon as SPP1. Therefore, we limit at present the SPP1-like genus proposal to its close genetic neighborhood until more robust phylogenetic and biological evidence is obtained to expand the genus to other phages.

**Author Contributions:** Conceptualization, L.M.G., M.E.S.F., S.A. and P.T.; Investigation, L.M.G., M.E.S.F., C.M., L.J., A.L., K.D. and S.A.; Supervision, S.A., R.C.-L. and P.T.; Validation, L.M.G., M.E.S.F., I.A., K.D., L.O., R.C.-L., S.A. and P.T.; Visualization, L.M.G., M.E.S.F. and P.T.; Writing—original draft P.T.; Writing—review and editing L.M.G., M.E.S.F., I.A., L.O., K.D., A.L. and S.A.

**Funding:** Research in our teams was funded by the CNRS, INRA, MINECO grant BFU2015-67065-P (to S.A.), ANR grants 09-BLAN-0149-0 (to P.T.), ANR-12-Blanc-BSV3-0021 (to P.T. and R.C.L.), ANR-15-CE11-0010-01 (to P.T. and R.C.L.), and an "Equipe de la Fondation Médicale (FRM)" grant (to P.T.).
**Acknowledgments:** We thank past and present researchers working on bacteriophage SPP1 for their contributions to our current understanding of this exciting biological system.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Riva, S.; Polsinelli, M.; Falaschi, A. A new phage of *Bacillus subtilis* with infectious DNA having separable strands. *J. Mol. Biol.* **1968**, *35*, 347–356. [CrossRef]

2. Riva, S.; Polsinelli, M. Relationship between competence for transfection and for transformation. *J. Virol.* **1968**, *2*, 587–593. [PubMed]

3. Spatz, H.C.; Trautner, T.A. One way to do experiments on gene conversion? *Mol. Gen. Genet.* **1970**, *109*, 84–106. [CrossRef] [PubMed]

4. Spatz, H.C.; Trautner, T.A. The role of recombination in transfection of *B. subtilis*. *Mol. Gen. Genet.* **1971**, *113*, 174–190. [CrossRef] [PubMed]

5. Deichelbohrer, I.; Alonso, J.C.; Lüder, G.; Trautner, T.A. Plasmid transduction by *Bacillus subtilis* bacteriophage SPP1: Effects of DNA homology between plasmid and bacteriophage. *J. Bacteriol.* **1985**, *162*, 1238–1243.

6. Alonso, J.C.; Lüder, G.; Trautner, T.A. Intramolecular homologous recombination in *Bacillus subtilis* **168**. *Mol. Gen. Genet.* **1992**, *236*, 60–64. [CrossRef]

7. Behrens, B.; Lüder, G.; Behncke, M.; Trautner, T.A.; Ganesan, A.T. The genome of *B. subtilis* phage SPP1: Physical arrangement of phage genes. *Mol. Gen. Genet.* **1979**, *175*, 351–357. [CrossRef]

8. Ratcliff, S.W.; Luh, J.; Ganesan, A.T.; Behrens, B.; Thompson, R.; Montenegro, M.A.; Morelli, G.; Trautner, T.A. The genome of *Bacillus subtilis* phage SPP1: The arrangement of restriction endonuclease generated fragments. *Mol. Gen. Genet.* **1979**, *168*, 165–172. [CrossRef]

9. Santos, M.A.; Almeida, J.; de Lencastre, H.; Morelli, G.; Kamke, M.; Trautner, T.A. Genomic organization of the related *Bacillus subtilis* bacteriophages SPP1, 41c, rho 15, and SF6. *J. Virol.* **1986**, *60*, 702–707.

10. Yasbin, R.E.; Frank, E. Young Transduction in *Bacillus subtilis* by bacteriophage SPP1. *J. Virol.* **1974**, *14*, 1343–1348.

11. Ferrari, E.; Canosi, U.; Galizzi, A.; Mazza, G. Studies on Transduction Process by SPP1 Phage. *J. Gen. Virol.* **1978**, *41*, 563–572. [CrossRef] [PubMed]

12. Canosi, U.; Lüder, G.; Trautner, T.A. SPP1-mediated plasmid transduction. *J. Virol.* **1982**, *44*, 431–436. [PubMed]

13. Valero-Rello, A.; López-Sanz, M.; Quevedo-Olmos, A.; Sorokin, A.; Ayora, S. Molecular mechanisms that contribute to horizontal transfer of plasmids by the bacteriophage SPP1. *Front. Microbiol.* **2017**, *8*, 1–13. [CrossRef] [PubMed]

14. Alonso, J.C.; Lüder, G.; Trautner, T.A. Requirements for the formation of plasmid-transducing particles of *Bacillus subtilis* bacteriophage SPP1. *EMBO J.* **1986**, *5*, 3723–3728. [CrossRef] [PubMed]

15. White, H.E.; Sherman, M.B.; Brasilès, S.; Jacquet, E.; Seavers, P.; Tavares, P.; Orlova, E.V. Capsid structure and its stability at the late stages of bacteriophage SPP1 assembly. *J. Virol.* **2012**, *86*, 6768–6777. [CrossRef] [PubMed]

16. Plisson, C.; White, H.E.; Auzat, I.; Zafarani, A.; São-José, C.; Lhuillier, S.; Tavares, P.; Orlova, E.V. Structure of bacteriophage SPP1 tail reveals trigger for DNA ejection. *EMBO J.* **2007**, *26*, 3720–3728. [CrossRef] [PubMed]

17. Tavares, P.; Lurz, R.; Stiege, A.; Rückert, B.; Trautner, T.A. Sequential headful packaging and fate of the cleaved DNA ends in bacteriophage SPP1. *J. Mol. Biol.* **1996**, *264*, 954–967. [CrossRef] [PubMed]

18. Humphreys, G.O.; Trautner, T.A. Maturation of bacteriophage SPP1 DNA: Limited precision in the sizing of mature bacteriophage genomes. *J. Virol.* **1981**, *37*, 832–835.

19. Morelli, G.; Fisseau, C.; Behrens, B.; Trautner, T.A.; Luh, J.; Ratcliff, S.W.; Allison, D.P.; Ganesan, A.T. The genome of *B. subtilis* phage SPP1: The topology of DNA molecules. *Mol. Gen. Genet.* **1979**, *168*, 153–164. [CrossRef]

20. Tavares, P.; Santos, M.A.; Lurz, R.; Morelli, G.; De Lencastre, H.; Trautner, T.A. Identification of a gene in *Bacillus subtilis* bacteriophage SPP1 determining the amount of packaged DNA. *J. Mol. Biol.* **1992**, *225*, 81–92. [CrossRef]
21. Baptista, C.; Santos, M.A.; São-José, C. Phage SPP1 reversible adsorption to *Bacillus subtilis* cell wall teichoic acids accelerates virus recognition of membrane receptor YueB. *J. Bacteriol.* **2008**, *190*, 4989–4996. [CrossRef] [PubMed]

22. Baptista, C.; Barreto, H.C.; São-José, C. High levels of DegU-P activate an Esat-6-Like secretion system in *Bacillus subtilis*. *PLoS ONE* **2013**, *8*, e67840. [CrossRef] [PubMed]

23. São-José, C.; Lhuillier, S.; Lurz, R.; Melki, R.; Lepault, J.; Santos, M.A.; Tavares, P. The ectodomain of the viral receptor YueB forms a fiber that triggers ejection of bacteriophage SPP1 DNA. *J. Biol. Chem.* **2006**, *281*, 11464–11470. [CrossRef] [PubMed]

24. Jakutyte, L.; Baptista, C.; São-José, C.; Daugelavičius, R.; Carballido-López, R.; Tavares, P. Bacteriophage infection in rod-shaped gram-positive bacteria: Evidence for a preferential polar route for phage SPP1 entry in *Bacillus subtilis*. *J. Bacteriol.* **2011**, *193*, 4893–4903. [CrossRef] [PubMed]

25. Jakutyte, L.; Lurz, R.; Baptista, C.; Carballido-López, R.; São-José, C.; Tavares, P.; Daugelavičius, R. First steps of bacteriophage SPP1 entry into *Bacillus subtilis*. *Virology* **2012**, *422*, 425–434. [CrossRef] [PubMed]

26. Alonso, J.C.; Lüder, G.; Stiege, A.C.; Chai, S.; Weise, F.; Trautner, T.A. The complete nucleotide sequence and functional organization of *Bacillus subtilis* bacteriophage SPP1. *Gene* **1997**, *204*, 201–212. [CrossRef]

27. Burger, K. Biochemische und genetische Untersuchungen zur DNA-Replikation des *B. subtilis*-Bakteriophagen SPP1. Ph.D. Thesis, Freien Universität, Berlin, Germany, 1978.

28. Burger, K.J.; Trautner, T.A. Specific labelling of replicating SPP1 DNA: Analysis of viral DNA synthesis and identification of phage DNA-genes. *Mol. Gen. Genet.* **1978**, *166*, 277–285. [CrossRef] [PubMed]

29. Mastromei, G.; Riva, S.; Fietta, A.; Pagani, L. SPP1 DNA replicative forms: Growth of phage SPP1 in *Bacillus subtilis* mutants temperature-sensitive in DNA synthesis. *Mol. Gen. Genet.* **1978**, *167*, 157–164. [CrossRef]

30. Lo Piano, A.; Martínez-Jiménez, M.I.; Zecchi, L.; Ayora, S. Recombination-dependent concatemeric viral DNA replication. *Virus Res.* **2011**, *160*, 1–14. [CrossRef]

31. Zecchi, L.; Lo Piano, A.; Suzuki, Y.; Cañas, C.; Takeyasu, K.; Ayora, S. Characterization of the Holliday Junction resolving enzyme encoded by the *Bacillus subtilis* bacteriophage SPP1. *PLoS ONE* **2012**, *7*, e48440. [CrossRef]

32. Auzat, I.; Petitpas, I.; Lurz, R.; Weise, F.; Tavares, P. A touch of glue to complete bacteriophage assembly: The tail-to-head joining protein (THJP) family. *Mol. Microbiol.* **2014**, *91*, 1164–1178. [CrossRef] [PubMed]

33. Fernandes, S.; São-José, C. Probing the function of the two holin-like proteins of bacteriophage SPP1. *Virology* **2017**, *500*, 184–189. [CrossRef] [PubMed]

34. Van Dijk, E.L.; Auger, H.; Jaszczyszyn, Y.; Thermes, C. Ten years of next-generation sequencing technology. *Trends Genet.* **2014**, *30*, 418–426. [CrossRef] [PubMed]

35. Djacem, K.; Tavares, P.; Oliveira, L. Bacteriophage SPP1 pac cleavage: A precise cut without sequence specificity requirement. *J. Mol. Biol.* **2017**, *429*, 1381–1395. [CrossRef] [PubMed]

36. Solovyev, V.; Salamov, A. Automatic annotation of microbial genomes and metagenomic sequences. In *Metagenomics and Its Applications in Agriculture, Biomedicine, and Environmental Studies*; Li, R.W., Ed.; Nova Science Publisher's: Hauppauge, NY, USA, 2011; pp. 61–78, ISBN 978-1-61668-682-6.

37. Moran, C.P.; Lang, N.; LeGrice, S.F.J.; Lee, G.; Stephens, M.; Sonenshein, A.L.; Pero, J.; Losick, R. Nucleotide sequences that signal the initiation of transcription and translation in *Bacillus subtilis*. *Mol. Gen. Genet.* **1982**, *186*, 339–346. [CrossRef] [PubMed]

38. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. *J. Mol. Biol.* **1990**, *215*, 403–410. [CrossRef]

39. Madden, T. The BLAST Sequence Analysis Tool. Available online: https://www.ncbi.nlm.nih.gov/books/NBK153387/ (accessed on 22 August 2018).

40. Zimmermann, L.; Stephens, A.; Nam, S.-Z.; Rau, D.; Kübler, J.; Lozajic, M.; Gabler, F.; Söding, J.; Lupas, A.N.; Alva, V. A completely reimplemented MPI Bioinformatics toolkit with a new HHpred server at its core. *J. Mol. Biol.* **2018**, *430*, 2237–2243. [CrossRef] [PubMed]

41. Stüber, D.; Morelli, G.; Bujard, H.; Montenegro, M.A.; Trautner, T.A. Promoter sites in the genome of *B. subtilis* phage SPP1. *Mol. Gen. Genet.* **1981**, *181*, 518–521. [CrossRef] [PubMed]

42. Tailor, R.; Bensi, G.; Morelli, G.; Canosi, U.; Trautner, T.A. The genome of *Bacillus subtilis* phage SPP1: structure of an early promoter. *J. Gen. Microbiol.* **1985**, *131*, 1259–1262. [CrossRef]

43. Gautherset, D.; Lambert, A. Direct RNA motif definition and identification from multiple sequence alignments using secondary structure profiles. *J. Mol. Biol.* **2001**, *313*, 1003–1011. [CrossRef] [PubMed]
44. Hofacker, I.L.; Fontana, W.; Stadler, P.F.; Bonhoeffer, L.S.; Tacker, M.; Schuster, P. Fast folding and comparison of RNA secondary structures. *Monatshefte für Chemie* **1994**, *125*, 167–188. [CrossRef]

45. Lesnik, E.A.; Sampath, R.; Levene, H.B.; Henderson, T.J.; McNeil, J.A.; Ecker, D.J. Prediction of rho-independent transcriptional terminators in *Escherichia coli*. *Nucleic Acids Res.* **2001**, *29*, 3583–3594. [CrossRef] [PubMed]

46. Macke, T.J.; Ecker, D.J.; Gutell, R.R.; Gautheret, D.; Case, D.A.; Sampath, R. RNAMotif, an RNA secondary structure definition and search algorithm. *Nucleic Acids Res.* **2001**, *29*, 4724–4735. [CrossRef] [PubMed]

47. Moszer, I.; Rocha, E.P.C.; Danchin, A. Codon usage and lateral gene transfer in *Bacillus subtilis*. *Curr. Opin. Microbiol.* **1999**, *2*, 524–528. [CrossRef]

48. Gupta, S.K.; Ghosh, T.C. CUCG: A non-redundant codon usage database from complete genomes. *Curr. Sci.* **2000**, *78*, 28–29.

49. Kunst, F.; Ogasawara, N.; Moszer, I.; Albertini, A.M.; Alloni, G.; Azevedo, V.; Bertero, M.G.; Bessières, P.; Bolotin, A.; Borchert, S.; et al. The complete genome sequence of the gram-positive bacterium *Bacillus subtilis*. *Nature* **1997**, *390*, 249–256. [CrossRef] [PubMed]

50. Riva, S.C. Asymmetric transcription of *B. subtilis* phage SPP1 DNA in vitro. *Biochem. Biophys. Res. Commun.* **1969**, *34*, 824–830. [CrossRef]

51. Chenciner, N.; Milanesi, G. Restriction fragment analysis of bacteriophage SPP1 in vitro transcription by host RNA polymerase. *J. Virol.* **1978**, *28*, 95–105. [PubMed]

52. Montenegro, M.A.; Trautner, T.A. In vivo transcription of *Bacillus subtilis* bacteriophage SPP1. *Mol. Gen. Genet.* **1981**, *181*, 512–517. [CrossRef]

53. Desmyter, A.; Reeve, J.N.; Morelli, G.; Trautner, T.A. Inversion and deletion mutants in *Bacillus subtilis* bacteriophage SPP1 as a consequence of cloning. *Mol. Gen. Genet.* **1985**, *198*, 537–539. [CrossRef] [PubMed]

54. Milanesi, G.; Cassani, G. Transcription after bacteriophage SPP1 infection in *Bacillus subtilis*. *J. Virol.* **1972**, *10*, 187–192. [PubMed]

55. Chai, S.; Szepan, U.; Liider, G.; Trautner, T.A.; Alonso, J.C. Sequence analysis of the left end of the *Bacillus subtilis* bacteriophage SPP1 genome. *Gene* **1993**, *129*, 41–49. [CrossRef]

56. Pedré, X.; Weise, F.; Chai, S.; Lüder, G.; Alonso, J.C. Analysis of cis and trans acting elements required for the initiation of DNA replication in the *Bacillus subtilis* bacteriophage SPP1. *J. Mol. Biol.* **1994**, *236*, 1324–1340. [CrossRef]

57. Chai, S.; Bravo, A.; Lüder, G.; Nedlin, A.; Trautner, T.A.; Alonso, J.C. Molecular analysis of the *Bacillus subtilis* bacteriophage SPP1 region encompassing genes 1 to 6. The products of gene 1 and gene 2 are required for pac cleavage. *J. Mol. Biol.* **1992**, *224*, 87–102. [CrossRef]

58. Pope, W.H.; Jacobs-Sera, D.; Russell, D.A.; Rubin, D.H.F.; Kajee, A.; Msibi, Z.N.P.; Larsen, M.H.; Jacobs, W.R.; Lawrence, J.G.; Hendrix, R.W.; et al. Genomics and proteomics of mycobacteriophage patience, an accidental tourist in the *Mycobacterium* neighborhood. *MBio* **2014**, *5*, e02145. [CrossRef] [PubMed]

59. Erez, Z.; Steinberger-Levy, I.; Shamir, M.; Doron, S.; Stokar-Avihail, A.; Peleg, Y.; Melamed, S.; Leavitt, A.; Savidor, A.; Albeck, S.; et al. Communication between viruses guides lysis-lysogeny decisions. *Nature* **2017**, *541*, 488–493. [CrossRef] [PubMed]

60. Lurz, R.; Orlova, E.V.; Günther, D.; Dube, P.; Dröge, A.; Weise, F.; van Heel, M.; Tavares, P. Structural organisation of the head-to-tail interface of a bacterial virus. *J. Mol. Biol.* **2001**, *310*, 1027–1037. [CrossRef] [PubMed]

61. Sprengel, R.; Reiss, B.; Schaller, H. Translationally coupled initiation of protein synthesis in *Bacillus subtilis*. *Nucleic Acids Res.* **1985**, *13*, 893–909. [CrossRef] [PubMed]

62. Ma, J.; Campbell, A.; Karlin, S. Correlations between Shine-Dalgarno sequences and gene features such as predicted expression levels and operon structures. *J. Bacteriol.* **2002**, *184*, 5733–5745. [CrossRef]

63. Chai, S.; Lurz, R.; Alonso, J.C. The small subunit of the terminase enzyme of *Bacillus subtilis* bacteriophage SPP1 forms a specialized nucleoprotein complex with the packaging initiation region. *J. Mol. Biol.* **1995**, *252*, 386–398. [CrossRef]

64. Gual, A.; Alonso, J.C. Characterization of the small subunit of the terminase enzyme of the *Bacillus subtilis* bacteriophage SPP1. *Virology* **1998**, *242*, 279–287. [CrossRef] [PubMed]

65. Büttner, C.R.; Chechik, M.; Ortiz-Lombardía, M.; Smits, C.; Ebong, I.-O.; Chechik, V.; Jeschke, G.; Dykeman, E.; Benini, S.; Robinson, C.V.; et al. Structural basis for DNA recognition and loading into a viral packaging motor. *Proc. Natl. Acad. Sci. USA* **2012**, *109*, 811–816. [CrossRef] [PubMed]
66. Greive, S.J.; Fung, H.K.H.; Chechik, M.; Jenkins, H.T.; Weitzel, S.E.; Aguiar, P.M.; Brentnall, A.S.; Glousieau, M.; Gladyshev, G.V.; Potts, J.R.; et al. DNA recognition for virus assembly through multiple sequence-independent interactions with a helix-turn-helix motif. *Nucleic Acids Res.* **2016**, *44*, 776–789.
[CrossRef] [PubMed]

67. Camacho, A.G.; Gual, A.; Lurz, R.; Tavares, P.; Alonso, J.C. *Bacillus subtilis* bacteriophage SPP1 DNA packaging motor requires terminase and portal proteins. *J. Biol. Chem.* **2003**, *278*, 23251–23259. [CrossRef]

68. Chai, S.; Alonso, J.C. Distamycin-induced inhibition of formation of a nucleoprotein complex between the terminase small subunit G1P and the non-encapsidated end (pacL site) of *Bacillus subtilis* bacteriophage SPP1. *Nucleic Acids Res.* **1996**, *24*, 282–288. [CrossRef] [PubMed]

69. Chai, S.; Szepan, U.; Alonso, J.C. *Bacillus subtilis* bacteriophage SPP1 terminase has a dual activity: It is required for the packaging initiation and represses its own synthesis. *Gene* **1997**, *184*, 251–256. [CrossRef]

70. Gual, A.; Camacho, A.G.; Alonso, J.C. Functional analysis of the terminase large subunit, G2P, of *Bacillus subtilis* Bacteriophage SPP1. *J. Biol. Chem.* **2000**, *275*, 35311–35319. [CrossRef]

71. Smits, C.; Chechik, M.; Kovalevskiy, O.V.; Shevtsov, M.B.; Foster, A.W.; Alonso, J.C.; Antson, A.A. Structural basis for the nuclease activity of a bacteriophage large terminase. *EMBO Rep.* **2009**, *10*, 592–598. [CrossRef]

72. Cornilleau, C.; Atmane, N.; Jacquet, E.; Smits, C.; Alonso, J.C.; Tavares, P.; Oliveira, L. The nuclease domain of the SPP1 packaging motor coordinates DNA cleavage and encapsidation. *Nucleic Acids Res.* **2013**, *41*, 340–354. [CrossRef]

73. Oliveira, L.; Henriques, A.O.; Tavares, P. Modulation of the viral ATPase activity by the portal protein correlates with DNA packaging efficiency. *J. Biol. Chem.* **2006**, *281*, 21914–21923. [CrossRef]

74. Oliveira, L.; Cuervo, A.; Tavares, P. Direct interaction of the bacteriophage SPP1 packaging ATPase with the portal protein. *J. Biol. Chem.* **2010**, *285*, 7366–7373. [CrossRef] [PubMed]

75. Dube, P.; Tavares, P.; Lurz, R.; Van Heel, M. The portal protein of bacteriophage SPP1: a DNA pump with 13-fold symmetry. *EMBO J.* **1993**, *1*, 1303–1309. [CrossRef]

76. Jekow, P.; Behlke, J.; Tichelaar, W.; Lurz, R.; Regalla, M.; Hinrichs, W.; Tavares, P. Effect of the ionic environment on the molecular structure of bacteriophage SPP1 portal protein. *Eur. J. Biochem.* **1999**, *264*, 724–735. [CrossRef] [PubMed]

77. Orlova, E.V.; Dube, P.; Beckmann, E.; Zemlin, F.; Lurz, R.; Trautner, T.A.; Tavares, P.; Van Heel, M. Structure of the 13-fold symmetric portal protein of bacteriophage SPP1. *Nat. Struct. Biol.* **1999**, *6*, 842–846. [CrossRef] [PubMed]

78. Orlova, E.V.; Gowen, B.; Dröge, A.; Stiege, A.; Weise, F.; Lurz, R.; van Heel, M.; Tavares, P. Structure of a viral DNA gatekeeper at 10 Å resolution by cryo-electron microscopy. *EMBO J.* **2003**, *22*, 1255–1262. [CrossRef] [PubMed]

79. Isidro, A.; Henriques, A.O.; Tavares, P. The portal protein plays essential roles at different steps of the SPP1 DNA packaging process. *Virology* **2004**, *322*, 253–263. [CrossRef]

80. Isidro, A.; Santos, M.A.; Henriques, A.O.; Tavares, P. The high-resolution functional map of bacteriophage SPP1 portal protein. *Mol. Microbiol.* **2004**, *51*, 949–962. [CrossRef]

81. Chaban, Y.; Lurz, R.; Brasilès, S.; Cornilleau, C.; Karreman, M.; Zinn-Justin, S.; Tavares, P.; Orlova, E.V. Structural rearrangements in the phage head-to-tail interface during assembly and infection. *Proc. Natl. Acad. Sci. USA* **2015**, *112*, 7009–7014. [CrossRef]

82. Cuervo, A.; Vaney, M.-C.; Antson, A.A.; Tavares, P.; Oliveira, L. Structural rearrangements between portal protein subunits are essential for viral DNA translocation. *J. Biol. Chem.* **2007**, *282*, 18907–18913. [CrossRef]

83. Vinga, I.; Dröge, A.; Stiege, A.C.; Lurz, R.; Santos, M.A.; Daugelavičius, R.; Tavares, P. The minor capsid protein gp7 of bacteriophage SPP1 is required for efficient infection of *Bacillus subtilis*. *Mol. Microbiol.* **2006**, *61*, 1609–1621. [CrossRef]

84. Stiege, A.C.; Isidro, A.; Dröge, A.; Tavares, P. Specific targeting of a DNA-binding protein to the SPP1 procapsid by interaction with the portal oligomer. *Mol. Microbiol.* **2003**, *49*, 1201–1212. [CrossRef] [PubMed]

85. Becker, B.; de la Fuente, N.; Gassel, M.; Günther, D.; Tavares, P.; Lurz, R.; Trautner, T.A.; Alonso, J.C. Head morphogenesis genes of the *Bacillus subtilis* bacteriophage SPP1. *J. Mol. Biol.* **1997**, *268*, 822–839. [CrossRef] [PubMed]

86. Dröge, A.; Santos, M.A.; Stiege, A.C.; Alonso, J.C.; Lurz, R.; Trautner, T.A.; Tavares, P. Shape and DNA packaging activity of bacteriophage SPP1 procapsid: protein components and interactions during assembly. *J. Mol. Biol.* **2000**, *296*, 117–132. [CrossRef] [PubMed]
87. Poh, S.L.; el Khadali, F.; Berrier, C.; Lurz, R.; Melki, R.; Tavares, P. Oligomerization of the SPP1 scaffolding Protein. *J. Mol. Biol.* **2008**, *378*, 551–564. [CrossRef] [PubMed]

88. Zairi, M.; Stiege, A.C.; Nhiri, N.; Jacquet, E.; Tavares, P. The collagen-like protein gp12 is a temperature-dependent reversible binder of SPP1 viral capsids. *J. Biol. Chem.* **2014**, *289*, 27169–27181. [CrossRef] [PubMed]

89. Lhuillier, S.; Gallopin, M.; Gilquin, B.; Brasilès, S.; Lancelot, N.; Letellier, G.; Gilles, M.; Dethan, G.; Orlova, E.V.; Couprie, J.; et al. Structure of bacteriophage SPP1 head-to-tail connection reveals mechanism for viral DNA gating. *Proc. Natl. Acad. Sci. USA* **2009**, *106*, 8507–8512. [CrossRef] [PubMed]

90. Chagot, B.; Auzat, I.; Gallopin, M.; Petitpas, I.; Gilquin, B.; Tavares, P.; Zinn-Justin, S. Solution structure of gp17 from the Siphoviridae bacteriophage SPP1: Insights into its role in virion assembly. *Proteins Struct. Funct. Bioinform.* **2012**, *80*, 319–326. [CrossRef]

91. Auzat, I.; Dröge, A.; Weise, F.; Lurz, R.; Tavares, P. Origin and function of the two major tail proteins of bacteriophage SPP1. *Mol. Microbiol.* **2008**, *70*, 557–569. [CrossRef]

92. Langlois, C.; Ramboarina, S.; Cukkemane, A.; Auzat, I.; Chagot, B.; Gilquin, B.; Ignatiou, A.; Petitpas, I.; Kasotakis, E.; Paternostre, M.; et al. Bacteriophage SPP1 tail tube protein self-assembles into β-structure-rich tubes. *J. Biol. Chem.* **2015**, *290*, 3836–3849. [CrossRef]

93. Goulet, A.; Lai-Kee-Him, J.; Veesler, D.; Auzat, I.; Robin, G.; Shepherd, D.A.; Ashcroft, A.E.; Richard, E.; Lichière, J.; Tavares, P.; et al. The opening of the SPP1 bacteriophage tail, a prevalent mechanism in Gram-positive-infecting siphophages. *J. Biol. Chem.* **2011**, *286*, 25397–25405. [CrossRef]

94. Veesler, D.; Robin, G.; Lichière, J.; Auzat, I.; Tavares, P.; Bron, P.; Campanacci, V.; Cambillau, C. Crystal structure of bacteriophage SPP1 distal tail protein (gp19.1): A baseplate hub paradigm in gram-positive infecting phages. *J. Biol. Chem.* **2010**, *285*, 36666–36673. [CrossRef] [PubMed]

95. Vinga, I.; Baptista, C.; Auzat, I.; Petipas, I.; Lurz, R.; Tavares, P.; Santos, M.A.; São-José, C. Role of bacteriophage SPP1 tail spike protein gp21 on host cell receptor binding and trigger of phage DNA ejection. *Mol. Microbiol.* **2012**, *83*, 289–303. [CrossRef] [PubMed]

96. Veesler, D.; Blangy, S.; Spinelli, S.; Tavares, P.; Campanacci, V.; Cambillau, C. Crystal structure of *Bacillus subtilis* SPP1 phage gp22 shares fold similarity with a domain of lactococcal phage p2 RBP. *Protein Sci.* **2010**, *19*, 1439–1443. [CrossRef] [PubMed]

97. Veesler, D.; Blangy, S.; Lichière, J.; Ortiz-Lombardía, M.; Tavares, P.; Campanacci, V.; Cambillau, C. Crystal structure of *Bacillus subtilis* SPP1 phage gp23.1, a putative chaperone. *Protein Sci.* **2010**, *19*, 1812–1816. [CrossRef] [PubMed]

98. Fernandes, S.; São-José, C. More than a hole: The holin lethal function may be required to fully sensitize bacteria to the lytic action of canonical endolysins. *Mol. Microbiol.* **2016**, *102*, 92–106. [CrossRef] [PubMed]

99. Martínez-Jiménez, M.I.; Alonso, J.C.; Ayora, S. *Bacillus subtilis* bacteriophage SPP1-encoded gene 34.1 product is a recombination-dependent DNA replication protein. *J. Mol. Biol.* **2005**, *351*, 1007–1019. [CrossRef]

100. Ayora, S.; Missich, R.; Mesa, P.; Lurz, R.; Yang, S.; Egelman, E.H.; Alonso, J.C. Homologous-pairing activity of the *Bacillus subtilis* bacteriophage SPP1 replication protein G35P. *J. Biol. Chem.* **2002**, *277*, 35969–35979. [CrossRef]

101. Seco, E.M.; Zinder, J.C.; Manhart, C.M.; Lo Piano, A.; McHenry, C.S.; Ayora, S. Bacteriophage SPP1 DNA replication strategies promote viral and disable host replication in vitro. *Nucleic Acids Res.* **2013**, *41*, 1711–1721. [CrossRef]

102. Missich, R.; Weise, F.; Chai, S.; Lurz, R.; Pedré, X.; Alonso, J.C. The replisome organizer (G38P) of *Bacillus subtilis* bacteriophage SPP1 forms specialized nucleoprotein complexes with two discrete distant regions of the SPP1 genome. *J. Mol. Biol.* **1997**, *270*, 50–64. [CrossRef]

103. Ayora, S.; Stasiak, A.; Alonso, J.C. The *Bacillus subtilis* bacteriophage SPP1 G39P delivers and activates the G40P DNA helicase upon interacting with the G38P-bound replication origin. *J. Mol. Biol.* **1999**, *288*, 71–85. [CrossRef]

104. Bailey, S.; Sedelnikova, S.E.; Mesa, P.; Ayora, S.; Waltho, J.P.; Ashcroft, A.E.; Baron, A.J.; Alonso, J.C.; Rafferty, J.B. Structural analysis of *Bacillus subtilis* SPP1 phage helicase loader protein G39P. *J. Biol. Chem.* **2003**, *278*, 15304–15312. [CrossRef] [PubMed]

105. Ayora, S.; Weise, F.; Mesa, P.; Stasiak, A.; Alonso, J.C. *Bacillus subtilis* bacteriophage SPP1 hexameric DNA helicase, G40P, interacts with forked DNA. *Nucleic Acids Res.* **2002**, *30*, 2280–2289. [CrossRef] [PubMed]
106. Mesa, P.; Alonso, J.C.; Ayora, S. *Bacillus subtilis* Bacteriophage SPP1 G40P helicase lacking the N-terminal domain unwinds DNA bidirectionally. *J. Mol. Biol.* **2006**, *357*, 1077–1088. [CrossRef] [PubMed]

107. Ayora, S.; Langer, U.; Alonso, J.C. *Bacillus subtilis* DnaG primase stabilises the bacteriophage SPP1 G40P helicase-ssDNA complex. *FEBS Lett.* **1998**, *439*, 59–62. [CrossRef]

108. Wang, G.; Klein, M.G.; Tokonzaba, E.; Zhang, Y.; Holden, L.G.; Chen, X.S. The structure of a DnaB-family replicative helicase and its interactions with primase. *Nat. Struct. Mol. Biol.* **2008**, *15*, 94–100. [CrossRef] [PubMed]

109. Martínez-Jiménez, M.I.; Mesa, P.; Alonso, J.C. *Bacillus subtilis* tau subunit of DNA polymerase III interacts with bacteriophage SPP1 replicative DNA helicase G40P. *Nucleic Acids Res.* **2002**, *30*, 5056–5064. [CrossRef] [PubMed]

110. Rocha, E.P.C.; Danchin, A.; Viari, A. Translation in *Bacillus subtilis*: Roles and trends of initiation and termination, insights from a genome analysis. *Nucleic Acids Res.* **1999**, *27*, 3567–3576. [CrossRef] [PubMed]

111. Lovett, P.S.; Ambulos, N.P.; Mulbry, W.; Noguchi, N.; Rogers, E.J. UGA can be decoded as tryptophan at low efficiency in *Bacillus subtilis*. *J. Bacteriol.* **1991**, *173*, 1810–1812. [CrossRef]

112. Xu, J.; Hendrix, R.W.; Duda, R.L. Conserved translational frameshift in dsDNA bacteriophage tail assembly genes. *Mol. Cell* **2004**, *16*, 11–21. [CrossRef]

113. Xu, J.; Hendrix, R.W.; Duda, R.L. A balanced ratio of proteins from gene G and frameshift-extended gene GT is required for phage Lambda tail assembly. *J. Mol. Biol.* **2013**, *425*, 3476–3487. [CrossRef]

114. Xu, J.; Hendrix, R.W.; Duda, R.L. Chaperone–protein interactions that mediate assembly of the bacteriophage lambda tail to the correct length. *J. Mol. Biol.* **2014**, *426*, 1004–1018. [CrossRef] [PubMed]

115. Juhala, R.J.; Ford, M.E.; Duda, R.L.; Youlton, A.; Hatfull, G.F.; Hendrix, R.W. Genomic sequences of bacteriophages HK97 and HK022: Pervasive genetic mosaicism in the lambdoid bacteriophages. *J. Mol. Biol.* **2000**, *299*, 27–51. [CrossRef] [PubMed]

116. Hendrix, R.W.; Lawrence, J.G.; Hatfull, G.F.; Casjens, S. The origins and ongoing evolution of viruses. *Trends Microbiol.* **2000**, *8*, 504–508. [CrossRef]

117. Cumby, N.; Davidson, A.R.; Maxwell, K.L. The moron comes of age. *Bacteriophage* **2012**, *2*, e23146. [CrossRef]

118. Chai, S.; Kruft, V.; Alonso, J.C. Analysis of the *Bacillus subtilis* bacteriophages SPP1 and SF6 gene 1 product: a protein involved in the initiation of headful packaging. *Virology* **1994**, *202*, 930–939. [CrossRef]

119. Umene, K.; Shiraishi, A. Complete nucleotide sequence of *Bacillus subtilis* (natto) bacteriophage PM1, a phage associated with disruption of food production. *Virus Genes* **2013**, *46*, 524–534. [CrossRef] [PubMed]

120. Marks, T.J.; Hamilton, P.T. Characterization of a thermophilic bacteriophage of *Geobacillus kaustophilus*. *Arch. Virol.* **2014**, *159*, 2771–2775. [CrossRef] [PubMed]

121. Botstein, D. A theory of modular evolution for bacteriophages. *Ann. N. Y. Acad. Sci.* **1980**, *354*, 484–491. [CrossRef]

122. Hendrix, R.W.; Smith, M.C.; Burns, R.N.; Ford, M.E.; Hatfull, G.F. Evolutionary relationships among diverse bacteriophages and prophages: All the world's a phage. *Proc. Natl. Acad. Sci. USA* **1999**, *96*, 2192–2197. [CrossRef]

123. Alonso, J.C.; Tavares, P.; Lurz, R.; Trautner, T.A. Bacteriophage SPP1. In *The Bacteriophages*; Calendar, R., Ed.; Oxford University Press: New York, NY, USA, 2006; Volume 54, p. 746, ISBN 9780195148503.

124. Tavares, P.; Zinn-Justin, S.; Orlova, E.V. Genome gating in tailed bacteriophage capsids. In *Viral Molecular Machines*; Rossmann, M.G., Rao, V.B., Eds.; Advances in Experimental Medicine and Biology; Springer: Boston, MA, USA, 2012; Volume 726, pp. 585–600, ISBN 978-1-4614-0979-3.

125. Oliveira, L.; Tavares, P.; Alonso, J.C. Headful DNA packaging: bacteriophage SPP1 as a model system. *Virus Res.* **2013**, *173*, 247–259. [CrossRef]

126. Tavares, P. The bacteriophage head-to-tail interface. In *Subcellular Biochemistry*; Harris, J.R., Bhella, D., Eds.; Springer: Singapore, 2018; Volume 88, pp. 305–328, ISBN 978-981-10-8456-0.

127. Hendrix, R.W.; Duda, R.L. Bacteriophage lambda PaPa: Not the mother of all lambda phages. *Science* **1992**, *258*, 1145–1148. [CrossRef] [PubMed]

128. Shen, B.W.; Landthaler, M.; Shub, D.A.; Stoddard, B.L. DNA binding and bleavage by the HNH Homing Endonuclease I-HmuI. *J. Mol. Biol.* **2004**, *342*, 43–56. [CrossRef] [PubMed]

129. Edgell, D.R.; Gibb, E.A.; Belfort, M. Mobile DNA elements in T4 and related phages. *Virol. J.* **2010**, *7*, 290. [CrossRef]
130. Landthaler, M.; Lau, N.C.; Shub, D.A. Group I intron homing in *Bacillus* phages SPO1 and SP82: A gene conversion event initiated by a nicking homing endonuclease. *J. Bacteriol.* **2004**, *186*, 4307–4314. [CrossRef] [PubMed]

131. Quiles-Puchalt, N.; Carpena, N.; Alonso, J.C.; Novick, R.P.; Marina, A.; Penadés, J.R. Staphylococcal pathogenicity island DNA packaging system involving cos-site packaging and phage-encoded HNH endonucleases. *Proc. Natl. Acad. Sci. USA* **2014**, *111*, 6016–6021. [CrossRef] [PubMed]

132. Kala, S.; Cumby, N.; Sadowski, P.D.; Hyder, B.Z.; Kanelis, V.; Davidson, A.R.; Maxwell, K.L. HNH proteins are a widespread component of phage DNA packaging machines. *Proc. Natl. Acad. Sci. USA* **2014**, *111*, 6022–6027. [CrossRef]

133. Neamah, M.M.; Mir-Sanchis, I.; López-Sanz, M.; Acosta, S.; Baquedano, I.; Haag, A.F.; Marina, A.; Ayora, S.; Penadés, J.R. Sak and Sak4 recombinases are required for bacteriophage replication in *Staphylococcus aureus*. *Nucleic Acids Res.* **2017**, *45*, 6507–6519. [CrossRef]

134. Ogunleye, A.; Irorere, V.U.; Williams, C.; Hill, D.; Bhat, A.; Radecka, I. Poly-γ-glutamic acid: production, properties and applications. *Microbiology* **2015**, *161*, 1–17. [CrossRef]

135. Mamberti, S.; Prati, P.; Cremaschi, P.; Seppi, C.; Morelli, C.F.; Galizzi, A.; Fabbi, M.; Calvio, C. γ-PGA hydrolases of phage origin in *Bacillus subtilis* and other microbial genomes. *PLoS ONE* **2015**, *10*, e1–17. [CrossRef]

136. Kimura, K.; Itoh, Y. Characterization of poly-γ-glutamate hydrolase encoded by a bacteriophage genome: possible role in phage infection of *Bacillus subtilis* encapsulated with poly-γ-glutamate. *Appl. Environ. Microbiol.* **2003**, *69*, 2491–2497. [CrossRef]

137. Ghosh, K.; Senevirathne, A.; Kang, H.; Hyun, W.; Kim, J.; Kim, K.-P.; Ghosh, K.; Senevirathne, A.; Kang, H.S.; Hyun, W.; et al. Complete nucleotide sequence analysis of a novel *Bacillus subtilis*-infecting bacteriophage BSP10 and its effect on poly-gamma-glutamic acid degradation. *Viruses* **2018**, *10*, 240. [CrossRef] [PubMed]

138. Molshanski-Mor, S.; Yosef, I.; Kiro, R.; Edgar, R.; Manor, M.; Gershovits, M.; Laserson, M.; Pupko, T.; Qimron, U. Revealing bacterial targets of growth inhibitors encoded by bacteriophage T7. *Proc. Natl. Acad. Sci. USA* **2014**, *111*, 18715–18720. [CrossRef] [PubMed]

139. Bondy-Denomy, J.; Qian, J.; Westra, E.R.; Buckling, A.; Guttman, D.S.; Davidson, A.R.; Maxwell, K.L. Prophages mediate defense against phage infection through diverse mechanisms. *ISME J.* **2016**, *10*, 2854–2866. [CrossRef] [PubMed]

140. Dedrick, R.M.; Jacobs-Sera, D.; Bustamante, C.A.G.; Garlena, R.A.; Mavrich, T.N.; Pope, W.H.; Reyes, J.C.C.; Russell, D.A.; Adair, T.; Alvey, R.; et al. Prophage-mediated defence against viral attack and viral counter-defence. *Nat. Microbiol.* **2017**, *2*, 16251. [CrossRef] [PubMed]

141. Seco, E.M.; Ayora, S. *Bacillus subtilis* DNA polymerases, PolC and DnaE, are required for both leading and lagging strand synthesis in SPP1 origin-dependent DNA replication. *Nucleic Acids Res.* **2017**, *45*, 8302–8313. [CrossRef]

142. Alonso, J.C.; Ayora, S.; Canosa, I.; Weise, F.; Rojo, F. Site-specific recombination in Gram-positive theta-replicating plasmids. *FEMS Microbiol. Lett.* **1996**, *142*, 1–10. [CrossRef]

143. Stahl, F.W.; Murray, N.E. The evolution of gene clusters and genetic circularity in microorganisms. *Genetics* **1966**, *53*, 569–576.

144. Young, R. Phage lysis: Do we have the hole story yet? *Curr. Opin. Microbiol.* **2013**, *16*, 790–797. [CrossRef]

145. Krogh, S.; Jørgensen, S.T.; Devine, K.M. Lysis genes of the *Bacillus subtilis* defective prophage PBSX. *J. Bacteriol.* **1998**, *180*, 2110–2117.

146. Bläsi, U.; Young, R. Two beginnings for a single purpose: The dual-start holins in the regulation of phage lysis. *Mol. Microbiol. Lett.* **1996**, *21*, 675–682. [CrossRef]

147. Abrescia, N.G.A.A.; Bamford, D.H.; Grimes, J.M.; Stuart, D.I. Structure unifies the viral universe. *Annu. Rev. Biochem.* **2012**, *81*, 795–822. [CrossRef] [PubMed]

148. Veesler, D.; Cambillau, C. A common evolutionary origin for tailed-bacteriophage functional modules and bacterial machineries. *Microbiol. Mol. Biol. Rev.* **2011**, *75*, 423–433. [CrossRef] [PubMed]

149. Jamet, A.; Touchon, M.; Ribeiro-Gonçalves, B.; Carriço, J.A.; Charbit, A.; Nassif, X.; Ramirez, M.; Rocha, E.P.C. A widespread family of polymorphic toxins encoded by temperate phages. *BMC Biol.* **2017**, *15*, 75. [CrossRef] [PubMed]

150. Jamet, A.; Charbit, A.; Nassif, X. Antibacterial toxins: Gram-positive bacteria strike back! *Trends Microbiol.* **2018**, *26*, 89–91. [CrossRef] [PubMed]
151. Deichelbohrer, I.; Messer, W.; Trautner, T.A. Genome of *Bacillus subtilis* bacteriophage SPP1: structure and nucleotide sequence of *pac*, the origin of DNA packaging. *J. Virol.* **1982**, *42*, 83–90. [PubMed]

152. Bravo, A.; Alonso, J.C.; Trautner, T.A. Functional analysis of the *Bacillus subtilis* bacteriophage SPP1 pac site. *Nucleic Acids Res.* **1990**, *18*, 2881. [CrossRef] [PubMed]

153. Dröge, A.; Tavares, P. In vitro packaging of DNA of the *Bacillus subtilis* bacteriophage SPP1. *J. Mol. Biol.* **2000**, *296*, 103–115. [CrossRef] [PubMed]

154. Oliveira, L.; Alonso, J.C.; Tavares, P. A defined in vitro system for DNA packaging by the bacteriophage SPP1: Insights into the headful packaging mechanism. *J. Mol. Biol.* **2005**, *353*, 529–539. [CrossRef]

155. Tavares, P.; Dröge, A.; Lurz, R.; Graeber, I.; Orlova, E.; Dube, P.; Van Heel, M. The SPP1 connection. *FEMS Microbiol. Rev.* **1995**, *17*, 47–56. [CrossRef]

156. Bönemann, G.; Pietrosiuk, A.; Mogk, A. Tubules and donuts: a type VI secretion story: MicroReview. *Mol. Microbiol. Rev.* **2010**, *76*, 815–821. [CrossRef]

157. Scholl, D. Phage Tail-Like Bacteriocins. *Annu. Rev. Virol.* **2017**, *4*, 453–467. [CrossRef] [PubMed]

158. São-José, C.; Baptista, C.; Santos, M.A. *Bacillus subtilis* operon encoding a membrane receptor for bacteriophage SPP1. *J. Bacteriol.* **2004**, *186*, 8337–8346. [CrossRef] [PubMed]

159. King, J. Assembly of the tail of bacteriophage T4. *J. Mol. Biol.* **1968**, *32*, 231–262. [CrossRef]

160. Katsura, I.; Kühl, P.W. Morphogenesis of the tail of bacteriophage lambda: III. Morphogenetic pathway. *J. Mol. Biol.* **1975**, *91*, 257–273. [CrossRef]

161. Ghosh, N.; McKillop, T.J.; Jowitt, T.A.; Howard, M.; Davies, H.; Holmes, D.F.; Roberts, I.S.; Bella, J. Collagen-like proteins in pathogenic *E. coli* strains. *PLoS ONE* **2012**, *7*, e37872. [CrossRef] [PubMed]

162. Brüssow, H.; Desiere, F. Comparative phage genomics and the evolution of Siphoviridae: Insights from dairy phages. *Mol. Microbiol.* **2001**, *39*, 213–222. [CrossRef] [PubMed]

163. Casjens, S.R. Comparative genomics and evolution of the tailed-bacteriophages. *Curr. Opin. Microbiol.* **2005**, *8*, 451–458. [CrossRef]

164. Stockdale, S.R.; Collins, B.; Spinelli, S.; Douillard, F.P.; Mahony, J.; Cambillau, C.; van Sinderen, D. Structure and assembly of TP901-1 virion unveiled by mutagenesis. *PLoS ONE* **2015**, *10*, e0131676. [CrossRef]

165. Murphy, J.; Bottacini, F.; Mahony, J.; Kelleher, P.; Neve, H.; Zomer, A.; Nauta, A.; van Sinderen, D. Comparative genomics and functional analysis of the 936 group of lactococcal Siphoviridae phages. *Sci. Rep.* **2016**, *6*, 21345. [CrossRef]

166. Lopes, A.; Tavares, P.; Petit, M.-A.; Guérois, R.; Zinn-Justin, S. Automated classification of tailed bacteriophages according to their neck organization. *BMC Genom.* **2014**, *15*, 1027. [CrossRef]

167. Zivanovic, Y.; Confalonieri, F.; Ponchon, L.; Lurz, R.; Chami, M.; Flayhan, A.; Renouard, M.; Huet, A.; Decottignies, P.; Davidson, A.R.; et al. Insights into bacteriophage T5 structure from analysis of its morphogenesis genes and protein components. *J. Virol.* **2014**, *88*, 1162–1174. [CrossRef] [PubMed]

168. Fraser, J.S.; Yu, Z.; Maxwell, K.L.; Davidson, A.R. Ig-Like Domains on bacteriophages: a tale of promiscuity and deceit. *J. Mol. Biol.* **2006**, *359*, 496–507. [CrossRef] [PubMed]

169. Spinelli, S.; Veesler, D.; Bebeacua, C.; Cambillau, C. Structures and host-adhesion mechanisms of lactococcal siphophages. *Front. Microbiol.* **2014**, *5*, 3. [CrossRef] [PubMed]

170. 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. [CrossRef] [PubMed]

<figure><img src="image_5.png" /></figure>

© 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).