Next Article in Journal
Importance of Viral Late Domains in Budding and Release of Enveloped RNA Viruses
Previous Article in Journal
Post-COVID-19 Arthritis and Sacroiliitis: Natural History with Longitudinal Magnetic Resonance Imaging Study in Two Cases and Review of the Literature
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 13
Issue 8
10.3390/v13081557
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Isolation and Characterisation of the Bundooravirus Genus and Phylogenetic Investigation of the Salasmaviridae Bacteriophages

by
Cassandra R. Stanton
1,
Daniel T. F. Rice
1,
Michael Beer
2,
Steven Batinovic
1,† and
Steve Petrovski
1,*,†
1
Department of Physiology, Anatomy & Microbiology, La Trobe University, Melbourne, VIC 3086, Australia
2
Department of Defence Science and Technology, Port Melbourne, VIC 3207, Australia
*
Author to whom correspondence should be addressed.
†
These authors contributed equally.
Viruses 2021, 13(8), 1557; https://0-doi-org.brum.beds.ac.uk/10.3390/v13081557
Submission received: 29 June 2021 / Revised: 28 July 2021 / Accepted: 2 August 2021 / Published: 6 August 2021
(This article belongs to the Section Bacterial Viruses)
Browse Figures
Versions Notes

Abstract

:
Bacillus is a highly diverse genus containing over 200 species that can be problematic in both industrial and medical settings. This is mainly attributed to Bacillus sp. being intrinsically resistant to an array of antimicrobial compounds, hence alternative treatment options are needed. In this study, two bacteriophages, PumA1 and PumA2 were isolated and characterized. Genome nucleotide analysis identified the two phages as novel at the DNA sequence level but contained proteins similar to phi29 and other related phages. Whole genome phylogenetic investigation of 34 phi29-like phages resulted in the formation of seven clusters that aligned with recent ICTV classifications. PumA1 and PumA2 share high genetic mosaicism and form a genus with another phage named WhyPhy, more recently isolated from the United States of America. The three phages within this cluster are the only candidates to infect B. pumilus. Sequence analysis of B. pumilus phage resistant mutants revealed that PumA1 and PumA2 require polymerized and peptidoglycan bound wall teichoic acid (WTA) for their infection. Bacteriophage classification is continuously evolving with the increasing phages’ sequences in public databases. Understanding phage evolution by utilizing a combination of phylogenetic approaches provides invaluable information as phages become legitimate alternatives in both human health and industrial processes.

1. Introduction

Organisms that belong to the Bacillus genus are Gram-positive, aerobic, endospore forming rods [1]. They are diverse and important environmental microbes; however, members of this genus have been implicated with human disease [2,3]. These organisms most commonly include Bacillus cereus and Bacillus anthracis, which are capable of causing severe foodborne illnesses and anthrax, respectively [4]. Bacillus pumilus has also been known to intermittently cause food borne illness and contamination in assumed sterile areas [5]. This is largely due to the ability of B. pumilus endospores to survive extreme environments, including hydrogen peroxide treatment, which is a common method of equipment sterilization [6]. B. pumilus genomes also contain an arsenal of genes able to survive oxidative stress and antimicrobial compounds [7]. This alarming robustness, along with the looming threat of antibiotic resistance, suggests that alternative treatments to chemical sterilization and antibiotics are required.
Bacteriophages (or phages) are viruses that propagate by infecting and lysing bacterial cells. Phages are predicted to be the most diverse and abundant biological entities on the planet [8,9]. Due to their ability to alter bacterial genomes through horizontal gene transfer and impact the population dynamics within microbial communities, they play a major role in microbial ecology and evolution [10,11]. Phages have also gained great interest as potential therapeutic agents to be used as an alternative to antibiotics [12,13]. With the emergence of next generation sequencing, an abundance of phage sequences are continuously deposited into public databases, providing comprehensive information on phage genetic diversity and taxonomy [14,15,16,17]. The abundance of sequence data has led to the creation of the Mycobacteriophage databases that has subsequently become the Actinobacteriophage database (phagesDB.org) [18]. This website has expanded to also include a Bacillus phage database with over 1400 Bacillus phage sequences uploaded. This useful resource has provided us with information on the genetic diversity and evolution of Bacillus phages. It is becoming increasingly important that the more Bacillus phage genome sequences we have, the better we understand their evolution and interactions with the environment and their hosts [19].
Grose et al. [20] conducted a thorough analysis of 93 Bacillus phage sequences and generated clusters based on genomic similarities. The groupings were based on that described previously for the Mycobacterium phages [21] and the phi29-like phages [22]. Cluster B was noted for phages relating to phi29, now classified in the Salasvirus genus. At the time, a small subset of phage sequences was available and was classified into three sub-clusters, including phi29 and PZA (grouped in sub-cluster B1), B103 and Nf (sub-cluster B2), and GA-1 (sub-cluster B3) [20]. In 2018, Schilling, Hoppert, and Hertel [9] discussed 21 phi29-like phage sequences available in GenBank, providing preliminary insights into this group’s genetic relatedness [9]. However, due to advances in viral classification techniques and the addition of other phi29-like phages in Genbank, this work is outdated. Recently, Li, et al. [23] isolated one of the newest members of this group, DLc1, and concluded that the now present 30 phi29-like phages should be included within the Salasvirus genus. However, this study failed to use multiple reticulate phylogeny methods for this classification and relied only on basic genome comparison methods. As of earlier this year, the International Committee on Virus Taxonomy (ICTV) reviewed the phi29-like phages and classified them in the new Salasmaviridae family [24]. However, within this report there were no detailed conclusions about these new taxonomic rankings, therefore, validation of this work is needed.
In this study we isolated two novel B. pumilus phages, PumA1 and PumA2, from Australian soil samples. These phages were characterized based on morphology, host range, and genome sequence. Annotation and sequence analysis revealed that both phages are novel at the DNA sequence but share conserved protein families, genome structure, and phenotypic characteristics similar to phi29 and other Salasmaviridae phages. This study also provides an in-depth genomic analysis of the newly classified Salasmaviridae phages and insights into their evolution and diversity.

2. Materials and Methods

2.1. Bacterial Strains and Media

In this study, Bacillus pumilus LTU1 strain was used, isolated by Dr Darryl Reanney, from a soil sample collected from Victoria, Australia. Bacterial cultures were grown on LB (1% tryptone (Oxoid, Adelaide, Australia), 0.5% yeast extract (Oxoid), and 1% sodium chloride (Sigma)) broth or agar (LB plus 1.4% agar (Oxoid)) at 28 °C. All chemicals were obtained from Sigma, Sydney, Australia) unless otherwise stated.

2.2. Isolation and Purification of Phages and DNA Extraction

Bacteriophages were isolated from soil samples taken from multiple points within 30 km of Darwin, Australia. The soil samples (1 g) were suspended in 2 mL of sterile water. The mixture was vortexed for 60 s and then centrifuged for 5 min before being filtered through a cellulose acetate membrane filter (0.2-µm pore size) to remove cells and other solid matter. Following incubations, the remaining bacterial cells were removed by centrifugation and filtration through a 0.2-µm cellulose acetate membrane filter. Lawn plates of B. pumilus were prepared and 20 µL aliquots of enriched supernatants were applied onto the lawn plates and allowed to dry. Plates were incubated overnight and visually inspected for the presence of the plaques the following day. Single plaques were observed and isolated from two different soil samples in the Northern Territory in Australia. Plaques were purified through six rounds of dilution and re-isolation to ensure their purity.

2.3. Phage DNA Isolation, Genome Sequencing, and Annotation

Purified phage particles were polyethylene glycol (PEG) precipitated followed by a proteinase K treatment to extract DNA as described previously [25]. Isolated phage DNA (100 ng) were prepared using the NEBNext® Ultra™ II DNA Library Prep Kit (NEB) (Australia) followed by whole genome sequencing on an Illumina MiSeq v3 600-cycle kit with 300 bp paired-end reads. Raw data were filtered using Trim Galore v0.6.4 with the default settings (Q scores of ≥20, with automatic adapter detection) [26]. Phage genomes were assembled with SPAdes v3.12.0 with default settings. The genome termini were corrected upon manual inspection of raw sequencing reads using CLC Genomics Workbench v9.5.4 (Qiagen, Melbourne, Victoria, Australia).
Putative open reading frames (ORFs) were predicted using Glimmer v3 and manually confirmed [27]. Sequence similarity searches were conducted using the predicted amino acid sequences against the GenBank database and the BLASTP algorithm was used with an E-value significance cut off of 10−4 [28,29]. Conserved domains and motifs were identified using the conserved domain database (CDD) (http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/Structure/cdd/cdd.shtml) (accessed on 7 March 2019) and the Pfam database (http://pfam.sanger.ac.uk) (accessed on 7 March 2019) [30]. The presence of genes encoding tRNAs was screened for using ARAGORN (http://130.235.244.92/ARAGORN/) (accessed on 7 March 2019) [31].

2.4. Electron Microscopy

Copper grids (ProSciTech, Townsville, Queensland, Australia) coated with carbon and formvar were subjected to a glow discharge treatment for 60 s. A total of 20 µL of high titer (>109 pfu/mL) phage filtrates were placed onto the grids, incubated for 10 min, and followed by removal of excess residue with filter paper. Grids were washed twice with 5 µL MilliQ, each wash being absorbed with filter paper. The grids were then negatively stained with 3 µL of 2% (w/v) uranyl acetate, followed by immediate removal with filter paper and one final MilliQ wash as outlined previously. The grids were then left to dry in a laminar flow for 30 min before imaging. The grids were examined under a JEOL JEM02010HC electron microscope.

2.5. Nucleotide Sequence

The nucleotide sequences for phages vB_BpuP_PumA1 and vB_BpuP_PumA2 have been deposited in GenBank under the accession numbers MN524844 and MN524845, respectively.

2.6. Identification of Phage Resistant B. pumilus Mutants

Lawn plates of B. pumilus and either PumA1 or PumA2 flooded in high titer were grown for 24 h. Colonies that emerged in the clearings were picked, streaked out for a total of three times, and re-spotted with phage to test that their resistance was stable. DNA extractions of the wild type B. pumilus and the strains showing phage resistance were prepared using the Wizard Genomic DNA Purification Kit (Promega, Sydney, Australia) as per manufacturer’s instructions. The DNA samples were then prepared for next generation sequencing as outlined above. For the wild type B. pumilus, the raw sequencing data was trimmed using Trim Galore v0.6.4 and the genome assembled using Unicycler v0.4.8. For SNP analysis, Snippy (Galaxy V.4.5.0) was used. The assembled wild type B. pumilus was the reference genome and each mutant strain was compared for differences.

2.7. Whole Genome Analysis and Clustering

The NCBI and Bacillus Phage (phagesdb.org) Databases were examined for related phi29-like phages. Thirty-four complete phage sequences were found on either database and used for this study (Table 1). A dot plot using the Genome Pair Rapid Dotter program (GEPARD) was chosen as a preliminary guide for clustering [32]. The genome sequences were collated into one FASTA file in the order of nucleotide identity found through BLASTN results. Clusters were validated using VIRIDIC v1.0 software [33].
Whole genome alignments of all phages and their clusters were undertaken using BLASTn comparisons and visualized with Easyfig [34]. For gene content analysis, a Clustal Omega alignment was performed and visualized in Splitstree4 [35,36]. Using the Neighbor-Joining method, a consensus network tree was developed. To create the reticulate gene sharing network, vConTACT v2.0 0.9.17 [37], a gene mapping program, was used. Gene2Genome was first used to assign each protein coding sequence of each phage and map it to its contig/genome ID. This output file was then combined with the collated FASTA file previously used and ran in vConTACT v.2.0. The network output file of vConTACT 2v.2.0 was then visualized in Cytoscape v3.8.1 [38]. Finally, a core protein phylogeny tree was created using both the DNA polymerase and DNA encapsidation ATPase in VICTOR [39]. All pairwise comparisons of the nucleotide sequences were conducted using the Genome-BLAST Distance Phylogeny (GBDP) method and branch support was inferred from 100 pseudo-bootstrap replicates each. The tree was then visualized using iTOL v6.2 [40]. The tree was rooted at midpoint and meta-data overlayed.

3. Results

3.1. Isolation of Bacteriophages and Their Morphological Features

After screening multiple soil samples collected in Darwin (Northern Territory, Australia), two samples from different locations produced plaques on lawn plates of B. pumilus (LTU1). They were isolated, purified, and named vB_BpuP_PumA1 (PumA1) and vB_BpuP_PumA2 (PumA2). For host range analysis, the phages were tested against other Bacillus species in our culture collection including B. anthracis, B. subtilis, B. mycoides, and B. cereus. Both phages exclusively lysed B. pumilus, suggesting a narrow host range (regarding available strains). Transmission electron microscopy imaging of negatively stained phages demonstrated that PumA1 and PumA2 displayed short tails and small, elongated capsids of 52 ± 5 nm (length) × 29 ± 6 nm (width) (Figure 1).

3.2. Sequencing and Genomic Features of PumA1 and PumA2

PumA1 and PumA2 were sequenced using Illumina sequencing technology. Genome assembly revealed both phages had linear genomes of 18,466 bp and 18,932 bp, respectively. Annotation of the PumA1 and PumA2 genomes revealed that they contained 26 and 28 putative open reading frames, respectively, and no tRNA genes (Figure 2). The genomes of both PumA1 and PumA2 share 82% similarity over 98% of the genomes. When compared to DNA sequences in the GenBank database, both phages are unique, with only 4% to 20% of the genomes sharing between 65% and 73% sequence identities with other phi29-like phage genomes. The phage genomes were flanked by 11 bp inverted repeat sequences (5′-AATGTAAAGGT-3′) consistent with phi29-like phages that all contain variations of inverted repeats [53]. The predicted amino acid sequence of each open reading frame was used in a BLASTp analysis to determine the closest homologue. Predicted functions can be assigned to sixteen gene products (Table S1). The genes were numbered consecutively with the exceptions of orf2.1 and orf2.2, which are present in PumA2 but not in PumA1. The genome organization and structure are similar to that of phage phi29 and its relatives, and the gene products share conserved similarities at the amino acid level [54].
The genomes of both PumA1 and PumA2 can be separated into three different modules based on the direction the genes are transcribed (Regions I, II, and III) (Figure 2). Regions I and III (also known as the early genes) contain genes that are transcribed in the same direction and are located at the 5′ and 3′ ends of the genomes. Region II (also referred to as the late genes) is located in the center of the genome and is transcribed in the opposite direction to the other genes [9].
Region I of the genomes encompasses eleven genes for PumA1 and thirteen genes for PumA2. The proteins encoded by the open reading frames orf1-orf6 in both phage genomes have no known predicted function and are noted as hypothetical proteins. However, orf7 to orf11 are described and associated with the phi29-like phages [9,44,45]. The genes encoding the DNA polymerase (orf7) are highly conserved between the phi29-like phages, and when compared using BLASTn, most of this gene (83% coverage) is 65% homologous to the DNA polymerase of phage phi29. The terminal binding protein in PumA1 and PumA2 is encoded by orf8, adjacent to the DNA polymerase [55]. The remaining genes within the first region encode a DNA transcriptional activator for the late genes (orf9), containing a characteristic conserved motif pfam05464 similar to that observed in phi29 [56,57], followed by a gene encoding single stranded binding protein (orf10) and double stranded binding protein (orf11). Region III contains four genes (orf23-orf26) transcribed in the same direction as the genes in region I. Three of the genes have unknown functions, and orf24 encodes a DNA replication organizer with a pfam06720 motif.
Region II contains genes orf12-orf22 that encode structural and morphogenesis genes. The head morphogenesis protein, orf12, contains a pfam11418 motif, similar to phi29 scaffolding protein [58]. This is followed by orf13, a putative major head protein containing a bacterial Ig-like domain (pfam02368) and orf14, a head fiber protein and motif (pfam11133). The major tail protein (orf15) contains a pfam16838 motif, conserved across groups of podoviruses [50,59]. The proteins that connect the phage head and tail are encoded by orf16 containing a pfam05352 [60], orf17 encoding a lower collar protein with the PHA00148 motif [61,62], and orf18 encoding a minor structural protein that is suspected to form the pre-neck appendage protein. Orf18’s closest homologues are found in other phi29-like phages WhyPhy and SRT01hs and a Staphylococcus phage ST134 with a conserved motif TIGR04523. This is followed by another morphogenesis protein, orf19, with the characteristic motif pfam01551. A putative holin is encoded by orf20 with the pfam05105 motif conserved in bacteriophage holin proteins. Orf21 encodes an endolysin with two motifs, pfam01520 and pfam01476 [63]. The last gene in the module orf22 encodes a protein that is predicted to encode a podovirus DNA encapsidation protein with the characteristic motif pfam05894 [64].

3.3. Whole Genome Comparison and Clustering of phi29-Like Phages

Since the genome organization and protein homologies of PumA1 and PumA2 were reminiscent of phi29, we next investigated the evolutionary relationship shared between PumA1, PumA2, and other phi29-like phages. The GenBank and Bacillus phage databases were searched for complete phage genomes that shared genomic similarities to phi29. This led to the identification of 34 phages, including the two isolated in this study (Table 1). Hatfull, et al. [65] previously described a classification system using a four-method approach to cluster 60 mycobacteriophage genomes. We employed their approach as a guide and included whole genome comparisons, network and clustering analysis, and candidate gene phylogenetics to organize the phi29-like phages into respective clusters (Figure S1). The clusters were expanded from the ones previously described by Grose, Jensen, Burnett, and Breakwell [20] and employ the same naming style of B and subclusters numbered (e.g., B1). The number of clusters were expanded from three [20] to seven and three singleton phages. New members were added to existing clusters and align with the ICTV classifications. While no new clusters or potential genera were formed from this analysis, this study provides a justification for the ICTV taxonomic rankings and framework for future clustering and classification of phages.

3.4. Dot Plot Analysis and Genomic Identities

Dot plot analysis of the 34 genomes revealed six clusters and two phages not pertaining to a cluster, referred to as singletons (Figure S2). To further define the clusters and provide a numerical value to their similarities, we used the Virus Intergenomic Distance Calculator (VIRIDIC) to calculate intergenomic similarities between each phage. VIRIDIC combines several similarity algorithms with genome alignment and length ratios to capture overall relatedness of prokaryotic viruses [33]. In correlation with the dot plot analysis, the heatmap presents the 34 phages in the same clustering pattern (Figure 3). Each cluster shares similarities of between 65.65–99.87% (Table S2). The clusters are now referred to as from B1 to B7, expanding from previous clusters [20].

3.5. Genome Map Alignments

Each cluster was aligned to show the dissimilarities between genomes (Figure 4). PumA1 and PumA2 contain an extra 796 bp and 1240 bp respectively compared to the ancestral phage phi29 in the 5′ early gene region. These hyperplastic regions are common throughout the phi29-like phages, where large insertions of up to 8256 bp are seen in the largest phi29-like phage, DLc1. Each cluster also shares the same pattern of insertions between them. The additional genes have no known functions but are presumed to be involved in the infection or replication processes since they are located within that region. Interestingly, the phages with genome sizes over 20 kb no longer contain the head fiber gene, a characteristic feature of phi29. The singletons show the least similarity to phi29 and other phages in this group, with small regions of similarity to their closest related phage.

3.6. Gene Sharing Networks

Reticulate networks have been recently shown to provide an accurate representation of phage relatedness versus traditional rooted phylogenetics, since phages undergo many recombination and horizontal gene transfer events [66]. To test if these methods aligned with the comparison methods previously mentioned, two reticulate methods were used. Firstly, an unrooted phylogenetic network was created in Splitstree4 using whole genome CLUSTAL Omega alignment (Figure 5). This network showed a consistent pattern of clustering in agreement with the other techniques.
vConTACT v.2.0 was then used in conjunction with the Splitstree network. vConTACT v.2.0 is a newly developed software for virus classification that extracts, aligns, and clusters all predicted input proteins [37]. The protein clusters are then used to calculate viral clusters (VC) by am “edge” weight or statistical confidence due to the amount of protein clusters that each phage shares. This is compared to a global network of phages in the GenBank database. Figure 6 depicts the global network produced from vConTACT v.2.0 with various phages color coded by their host genera. The phi29-like phages are not connected to the main network, showing little gene-sharing outside of this group of phages. When the phi29-like network is expanded, there is a clear differentiation of each phage cluster, denoted by the individual cluster colors. While these phages do not connect to the main network, they share proteins with two other phages that are not infective for Bacillus, Lactococcus phage Asccphi28 [67] and Weisella phage phiYS61 [68], which are genotypically similar to the phi29-like phages. This form of reticulate phylogenetics helped corroborate the established clusters previously outlined and is an accurate tool for investigating phage gene-sharing and evolution.

3.7. Candidate Gene Analysis

Finally, a phylogenetic tree was created combining two conserved and integral genes, DNA polymerase and DNA encapsidation ATPase (Figure 7). The tree correlates with the whole genome comparison approaches as the phages are seen to group with their clusters. The tree splits into two distinctive branches with singletons DLc1 and MG-B1 and clusters B6 and B7 diverging from the rest of the clusters. This correlates with the genome sizes of these phages, as they are the largest genomes of the phi29-like group. The branch distances are also short within the clusters, particularly B6 members, signifying small base pair substitutions. The other factors that are seen to contribute to the phages’ evolution are host species and country or region they were isolated. These clusters also agree with the ICTV taxonomic rankings as outlined by metadata on the tree.

3.8. PumA1 and PumA2 Host Receptor Site

PumA1 and PumA2 displayed a narrow host range, only able to infect the B. pumilus strain in our collection. This appears to be a characteristic of the B4 cluster or Bundooravirus genus. It has been shown that phi29 requires polymerized teichoic acid for its attachment to B. subtilis [69]. Given the observed host range differences between phi29 and the two phages isolated in this study, we pursued an investigation into the host receptor of PumA1 and PumA2. B. pumilus colonies that developed stable resistance to PumA1 and PumA2 were isolated, and whole genome sequenced to determine which genomic modifications were causing resistance to phages. Mutations were found in either the tagF or tagT genes, which are a part of the teichoic acid synthesis operon [70,71,72]. Four of the seven variant strains isolated had modifications to the tagF gene, and the other three had mutations in the tagT gene (Table S3). The majority of the mutants (A2M1, A2M11, A1M3, A2M14, and A1M5) contained frameshift mutations resulting in early termination of protein translation. Most of these frameshifts occurred near the N-terminus of the respective proteins likely resulting in non-functional TagT or TagF proteins. The remaining mutant had a single amino acid substitution in TagF (G688S). It is unclear how this mutation affects protein structure and function.

4. Discussion

Advances in whole genome sequencing techniques and the rise of antibiotic resistance has resulted in an abundance of publicly available sequencing data and reinvigoration of phage-based studies. In this study, we isolated two novel, narrow host range Bacillus pumilus phages, PumA1 and PumA2. Both phages were unable to form plaques on B. pumilus strains that had mutations in the tagT or tagF genes. Wall teichoic acid (WTA) is a major component of Gram-positive cell walls, with the tag operons encoding the necessary machine for its synthesis [70,71]. TagF is a poly(glycerol phosphate) polymerase that plays a key role in the formation of the glycerol phosphate chains in WTA. Single nucleotide polymorphisms have been shown to significantly alter the function of TagF, resulting in decreased polymerase activity [71]. Whereas tagT, which encodes a LytR-CspA-Psr (LCP) family protein, is implicated in the final stage of WTA synthesis, in which it catalyzes the attachment between the teichoic acid polymers and peptidoglycan. In B. subtilis mutants in which the tagTUV operon was knocked out, cells appear more rounded and lack the typical rod shape and teichoic acid present in growth cultures but are detached from peptidoglycan [70]. Given that both phages studied here and phi29 appear to share the same surface receptor, we hypothesize that the Bundooravirus phages have differentiated tail structures to specifically detect B. pumilus teichoic acid.
The Salasmaviridae (or phi29-like) group contain the smallest known genome sizes of phages that infect Bacillus sp. The highly conserved phages of this family are distributed throughout several continents, with PumA1 and PumA2 currently the only phages isolated from Australia. The multi-approach clustering indicated a correlation of Salasmaviridae phages by host range and genome size. This was aptly demonstrated with the Bundooravirus phages isolated in this study, which appeared to have a narrow host range against B. pumilus. We also noticed a correlation of cluster formation with the geographical location of phage isolation. This pattern is likely driven by the host bacterium’s biogeography and population; this may become clearer as more Salasmaviridae phages are isolated and sequenced [72,73].
Since the Salasmaviridae phages follow a strict lytic lifecycle with no evidence of lysogenic activity, there may be low gene-content flux and recombination event rates between themselves, their hosts, and other phages [74]. However, environmental and host pressures can naturally result in genomic mutations, and as a phage gains new adaptations such as expanded host range, this can result in enough genetic variance to be included as a new species and genus [19,66,75]. While the phage genomes remain well conserved in their arrangement and modules of essential genes, there appears to be three main regions where insertions are present, all in the “early” replication gene modules. These additional genes correlate to increased genome size and how the phages cluster, reminiscent of hyperplastic regions seen in other groups of phages [76]. While most of these early genes have unknown functions, it is hypothesized that early proteins are associated with phage-host interactions. This includes genes that encode for protection from host degradation and restriction, anti-CRISPR, and inhibition of host transcription [77]. Though they are not essential for phage function, and can be lost and gained readily, they have the potential to be advantageous for adaption to their host [78]. Mutations leading to increased genome and capsid sizes tend to be more favorable and have been conserved throughout the Northopvirinae phages [79]. In contrast, the structural or “late” gene modules remain almost identical in arrangement and constitution throughout the Salasmaviridae phages. One notable difference seen in the “late genes” region is that any phage with a genome larger than 20 kb had lost the “head fiber” gene. The head fibers are involved in sensing and interact with the bacteria cell wall but are not essential for phage viability [80,81]. Phages that have acquired new genes are seen to have an enlarged capsid and altered capsid architecture in response [79]. It is hypothesized that since the fibers have no essential functions and the addition of newer genes has either forced or allowed these phages to expand their capsids, the head fibers are no longer able to attach to this new capsid structure. Nevertheless, this needs to be investigated further.

5. Conclusions

PumA1 and PumA2 represent a novel genus in the newly formed Salasmaviridae family, Bundooravirus. Multiple clustering approaches, including reticulate networks, resulted in the clustering of all current phi29-like phages, including those recently classified by the ICTV. This clustering agrees with taxonomic rankings and allows for the addition of several phages into this family. Thornton and Baseball_field should be included in the Claudivirus genus, BSTP4 should be classified in the Salasvirus genus, and DLc1 should be classified into the Northopvirinae sub-family but in its own genus. WhyPhy should be classified into the Bundooravirus genus with PumA1 and PumA2. Current data suggests the Bundooravirus phages are unique in their host range and require B. pumilus specific teichoic acid residues for infection. The Salasmaviridae (or phi29-like) phages are globally distributed but remain well conserved in genome organization and protein domains. However, there are three patterns of clustering that contribute to their evolution and classification, including geographical location, host range, and genome size. This study demonstrated that a combination of whole genome comparisons and rooted and reticulate phylogenetic models can be used to our advantage to order and classify phages. Understanding phage evolution and their relationships with other phages and the environment will provide us with invaluable information into phage phylogenetics as their usage in medical and industrial processes continues.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/v13081557/s1, Figure S1: Workflow for whole genome comparisons; Figure S2: Whole genome dot plot of the Salasmaviridae/phi29-like phages; Table S1: PumA1 and PumA2 genome annotations; Table S2: VIRDIC intergenomic similarities between clusters; Table S3: B. pumilus mutant strains resistant to PumA1 and PumA2 infection.

Author Contributions

Conceptualization, C.R.S., S.B. and S.P.; methodology, C.R.S.; formal analysis, C.R.S., S.B. and S.P.; investigation, C.R.S., D.T.F.R., M.B., S.B. and S.P.; data curation, C.R.S., S.B. and S.P.; writing—original draft preparation, C.R.S., S.B. and S.P.; writing—review and editing, C.R.S., M.B., S.B. and S.P.; supervision, M.B., S.B. and S.P.; funding acquisition, S.P.; resources, S.P. All authors have read and agreed to the published version of the manuscript.

Funding

C.R.S. was supported by a La Trobe University postgraduate award and a Defence Science Institute (DSI) RHD student grant. This work was funded by a DSI CERA grant (VUP17-007), Defence Science and Technology (collaborative project agreement ID8047) and La Trobe University funds provided by the department.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The nucleotide sequences for phages vB_BpuP_PumA1 and vB_BpuP_PumA2 have been deposited in GenBank under the accession numbers MN524844 and MN524845, respectively.

Acknowledgments

The soil that was collected and used for the isolation of PumA1 and PumA2 was taken from Darwin, and we would like to acknowledge the Larrakia people as the traditional owners of the Darwin region and pay respects to Larrakia elders past and present. We would also like to acknowledge and pay respects to the Wurundjeri People of the Kulin nation (Melbourne) where the remainder of this work took place.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study, in the collection, analyses, or interpretation of data, in the writing of the manuscript, or in the decision to publish the results.

References

  1. Logan, N.A.; Halket, G. Developments in the taxonomy of aerobic, endospore-forming bacteria. In Endospore-forming Soil Bacteria; Springer: Heidelberg, Germany, 2011; pp. 1–29. [Google Scholar]
  2. Santoyo, G.; Orozco-Mosqueda, M.D.C.; Govindappa, M. Mechanisms of biocontrol and plant growth-promoting activity in soil bacterial species of Bacillus and Pseudomonas: A review. Biocontrol Sci. Technol. 2012, 22, 855–872. [Google Scholar] [CrossRef]
  3. Glasset, B.; Herbin, S.; Granier, S.A.; Cavalié, L.; Lafeuille, E.; Guérin, C.; Ruimy, R.; Casagrande-Magne, F.; Levast, M.; Chautemps, N.; et al. Bacillus cereus, a serious cause of nosocomial infections: Epidemiologic and genetic survey. PLoS ONE 2018, 13, e0194346. [Google Scholar] [CrossRef]
  4. Lorenz, L.; Lins, B.; Barrett, J.; Montgomery, A.; Trapani, S.; Schindler, A.; Christie, G.E.; Cresawn, S.G.; Temple, L. Genomic characterization of six novel Bacillus pumilus bacteriophages. Virology 2013, 444, 374–383. [Google Scholar] [CrossRef] [Green Version]
  5. From, C.; Hormazabal, V.; Granum, P.E. Food poisoning associated with pumilacidin-producing Bacillus pumilus in rice. Int. J. Food Microbiol. 2007, 115, 319–324. [Google Scholar] [CrossRef]
  6. Kempf, M.J.; Chen, F.; Kern, R.; Venkateswaran, K. Recurrent Isolation of Hydrogen Peroxide-Resistant Spores of Bacillus pumilus from a Spacecraft Assembly Facility. Astrobiology 2005, 5, 391–405. [Google Scholar] [CrossRef]
  7. Handtke, S.; Schroeter, R.; Jürgen, B.; Methling, K.; Schlüter, R.; Albrecht, D.; van Hijum, S.; Bongaerts, J.; Maurer, K.-H.; Lalk, M.; et al. Bacillus pumilus Reveals a Remarkably High Resistance to Hydrogen Peroxide Provoked Oxidative Stress. PLoS ONE 2014, 9, e85625. [Google Scholar] [CrossRef]
  8. Rohwer, F. Global Phage Diversity. Cell 2003, 113, 141. [Google Scholar] [CrossRef] [Green Version]
  9. Schilling, T.; Hoppert, M.; Hertel, R. Genomic Analysis of the Recent Viral Isolate vB_BthP-Goe4 Reveals Increased Diversity of φ29-Like Phages. Viruses 2018, 10, 624. [Google Scholar] [CrossRef] [Green Version]
  10. Salmond, G.P.C.; Fineran, P. A century of the phage: Past, present and future. Nat. Rev. Genet. 2015, 13, 777–786. [Google Scholar] [CrossRef]
  11. Levin, B.R.; Bull, J.J. Population and evolutionary dynamics of phage therapy. Nat. Rev. Genet. 2004, 2, 166–173. [Google Scholar] [CrossRef]
  12. Torres-Barceló, C.; Hochberg, M.E. Evolutionary Rationale for Phages as Complements of Antibiotics. Trends Microbiol. 2016, 24, 249–256. [Google Scholar] [CrossRef]
  13. Kutateladze, M.; Adamia, R. Bacteriophages as potential new therapeutics to replace or supplement antibiotics. Trends Biotechnol. 2010, 28, 591–595. [Google Scholar] [CrossRef]
  14. Hatfull, G.F. Bacteriophage genomics. Curr. Opin. Microbiol. 2008, 11, 447–453. [Google Scholar] [CrossRef] [Green Version]
  15. Hatfull, G.F.; Hendrix, R.W. Bacteriophages and their genomes. Curr. Opin. Virol. 2011, 1, 298–303. [Google Scholar] [CrossRef] [Green Version]
  16. Tolstoy, I.; Kropinski, A.M.; Brister, J.R. Bacteriophage Taxonomy: An Evolving Discipline. In Bacteriophage Therapy. Methods in Molecular Biology; Humana Press: New York, NY, USA, 2017; pp. 57–71. [Google Scholar]
  17. Aziz, R.K.; Ackermann, H.-W.; Petty, N.K.; Kropinski, A.M. Essential Steps in Characterizing Bacteriophages: Biology, Taxonomy, and Genome Analysis. In Bacteriophages. Methods in Molecular Biology; Humana Press: New York, NY, USA, 2017; pp. 197–215. [Google Scholar]
  18. Russell, D.A.; Hatfull, G.F. PhagesDB: The actinobacteriophage database. Bioinformatics 2016, 33, 784–786. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Mavrich, T.; Hatfull, G.F. Bacteriophage evolution differs by host, lifestyle and genome. Nat. Microbiol. 2017, 2, 1–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Grose, J.H.; Jensen, G.L.; Burnett, S.H.; Breakwell, N.P. Correction: Genomic comparison of 93 Bacillus phages reveals 12 clusters, 14 singletons and remarkable diversity. BMC Genom. 2014, 15, 1184. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Hatfull, G.F. Mycobacteriophages: Genes and Genomes. Annu. Rev. Microbiol. 2010, 64, 331–356. [Google Scholar] [CrossRef]
  22. Meijer, W.J.J.; Horcajadas, J.A.; Salas, M. φ29 Family of Phages. Microbiol. Mol. Biol. Rev. 2001, 65, 261–287. [Google Scholar] [CrossRef] [Green Version]
  23. Li, C.; Yuan, X.; Li, N.; Wang, J.; Yu, S.; Zeng, H.; Zhang, J.; Wu, Q.; Ding, Y. Isolation and Characterization of Bacillus cereus Phage vB_BceP-DLc1 Reveals the Largest Member of the Φ29-Like Phages. Microorganisms 2020, 8, 1750. [Google Scholar] [CrossRef]
  24. Walker, P.J.; Siddell, S.G.; Lefkowitz, E.J.; Mushegian, A.R.; Adriaenssens, E.M.; Dempsey, D.M.; Dutilh, B.E.; Harrach, B.; Harrison, R.L.; Hendrickson, R.C.; et al. Changes to virus taxonomy and the Statutes ratified by the International Committee on Taxonomy of Viruses (2020). Arch. Virol. 2020, 165, 2737–2748. [Google Scholar] [CrossRef]
  25. Petrovski, S.; Seviour, R.J.; Tillett, D. Genome Sequence and Characterization of the Tsukamurella Bacteriophage TPA2. Appl. Environ. Microbiol. 2011, 77, 1389–1398. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 2011, 17. [Google Scholar] [CrossRef]
  27. Kearse, M.; Moir, R.; Wilson, A.; Stones-Havas, S.; Cheung, M.; Sturrock, S.; Buxton, S.; Cooper, A.; Markowitz, S.; Duran, C.; et al. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 2012, 28, 1647–1649. [Google Scholar] [CrossRef]
  28. Altschul, S.F.; Madden, T.L.; Schäffer, A.A.; Zhang, J.; Zhang, Z.; Miller, W.; Lipman, D.J. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res. 1997, 25, 3389–3402. [Google Scholar] [CrossRef] [Green Version]
  29. Zhang, Z.; Schwartz, S.; Wagner, L.; Miller, W. A Greedy Algorithm for Aligning DNA Sequences. J. Comput. Biol. 2000, 7, 203–214. [Google Scholar] [CrossRef]
  30. Finn, R.D.; Coggill, P.; Eberhardt, R.; Eddy, S.R.; Mistry, J.; Mitchell, A.; Potter, S.; Punta, M.; Qureshi, M.; Sangrador, A.; et al. The Pfam protein families database: Towards a more sustainable future. Nucleic Acids Res. 2015, 44, D279–D285. [Google Scholar] [CrossRef]
  31. Laslett, D. ARAGORN, a program to detect tRNA genes and tmRNA genes in nucleotide sequences. Nucleic Acids Res. 2004, 32, 11–16. [Google Scholar] [CrossRef]
  32. Krumsiek, J.; Arnold, R.; Rattei, T. Gepard: A rapid and sensitive tool for creating dotplots on genome scale. Bioinformatics 2007, 23, 1026–1028. [Google Scholar] [CrossRef] [Green Version]
  33. Moraru, C.; Varsani, A.; Kropinski, A. VIRIDIC—A Novel Tool to Calculate the Intergenomic Similarities of Prokaryote-Infecting Viruses. Viruses 2020, 12, 1268. [Google Scholar] [CrossRef]
  34. Sullivan, M.J.; Petty, N.; Beatson, S.A. Easyfig: A genome comparison visualizer. Bioinformatics 2011, 27, 1009–1010. [Google Scholar] [CrossRef]
  35. Sievers, F.; Wilm, A.; Dineen, D.; Gibson, T.J.; Karplus, K.; Li, W.; López, R.; McWilliam, H.; Remmert, M.; Söding, J.; et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol. Syst. Biol. 2011, 7, 539. [Google Scholar] [CrossRef]
  36. Huson, D.H. SplitsTree: Analyzing and visualizing evolutionary data. Bioinformatics 1998, 14, 68–73. [Google Scholar] [CrossRef] [PubMed]
  37. Bolduc, B.; Bin Jang, H.; Doulcier, G.; You, Z.-Q.; Roux, S.; Sullivan, M.B. vConTACT: An iVirus tool to classify double-stranded DNA viruses that infect Archaea and Bacteria. PeerJ 2017, 5, e3243. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Lopes, C.T.; Franz, M.; Kazi, F.; Donaldson, S.L.; Morris, Q.; Bader, G.D. Cytoscape Web: An interactive web-based network browser. Bioinformatics 2010, 26, 2347–2348. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Meier-Kolthof, J.P.; Göker, M. VICTOR: Genome-based phylogeny and classification of prokaryotic viruses. Bioinformatics 2017, 33, 3396–3404. [Google Scholar] [CrossRef] [Green Version]
  40. Letunic, I.; Bork, P. Interactive Tree Of Life (iTOL) v5: An online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021, 49, W293–W296. [Google Scholar] [CrossRef]
  41. Redondo, R.A.F.; Kupczok, A.; Stift, G.; Bollback, J.P. Complete Genome Sequence of the Novel Phage MG-B1 Infecting Bacillus weihenstephanensis. Genome Announc. 2013, 1. [Google Scholar] [CrossRef] [Green Version]
  42. Ji, X.; Zhang, C.; Fang, Y.; Zhang, Q.; Lin, L.; Tang, B.; Wei, Y. Isolation and characterization of glacier VMY22, a novel lytic cold-active bacteriophage of Bacillus cereus. Virol. Sin. 2015, 30, 52–58. [Google Scholar] [CrossRef] [PubMed]
  43. Guo, X.; Zhang, T.; Jin, M.; Zeng, R. Characterization of Bacillus phage Gxv1, a novel lytic Salasvirus phage isolated from deep-sea seamount sediments. Mar. Life Sci. Technol. 2020, 3, 13–19. [Google Scholar] [CrossRef]
  44. Kong, L.; Ding, Y.; Wu, Q.; Wang, J.; Zhang, J.; Li, H.-Y.; Yu, S.; Yu, P.; Gao, T.; Zeng, H.; et al. Genome sequencing and characterization of three Bacillus cereus-specific phages, DK1, DK2, and DK3. Arch. Virol. 2019, 164, 1927–1929. [Google Scholar] [CrossRef]
  45. Willms, I.M.; Hertel, R. Phage vB_BsuP-Goe1: The smallest identified lytic phage of Bacillus subtilis. FEMS Microbiol. Lett. 2016, 363, fnw208. [Google Scholar] [CrossRef] [Green Version]
  46. Pecenkova, T.; Benes, V.; Paces, J.; Vlček, Č.; Pačes, V. Bacteriophage B103: Complete DNA sequence of its genome and relationship to other Bacillus phages. Gene 1997, 199, 157–163. [Google Scholar] [CrossRef]
  47. Bradley, D.E. The Isolation and morphology of Some New Bacteriophages Specific for Bacillus and Acetobacter species. J. Gen. Microbiol. 1965, 41, 233–241. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Castilla-Llorente, V.; Salas, M.; Meijer, W.J.J. Different responses to Spo0A-mediated suppression of the related Bacillus subtilis phages Nf and φ29. Environ. Microbiol. 2009, 11, 1137–1149. [Google Scholar] [CrossRef] [PubMed]
  49. Pačes, V.; Viček, Č.; Urbánek, P.; Hostomský, Z. Nucleotide sequence of the major early region of Bacillus subtilis phage PZA, a close relative of φ29. Gene 1985, 38, 45–56. [Google Scholar] [CrossRef]
  50. Anderson, D.L.; Hickman, D.D.; Reilly, B.E. Structure of Bacillus subtilis Bacteriophage φ29 and the Length of φ29 Deoxyribonucleic Acid. J. Bacteriol. 1966, 91, 2081–2089. [Google Scholar] [CrossRef] [Green Version]
  51. Erill, I.; Caruso, S.M. Complete Genome Sequences of Three phi29-like Bacillus cereus Group Podoviridae. Genome Announc. 2017, 5, e00701-17. [Google Scholar] [CrossRef] [Green Version]
  52. Duperier, J.; Bulpitt, M.; Bispo, F.; Greguske, E. Genome Annotations of Two Bacillus Phages, Tomato and BaseballField. Microbiol. Resour. Announc. 2021, 10, e01196-20. [Google Scholar] [CrossRef]
  53. Yoshikawa, H.; Ito, J. Terminal proteins and short inverted terminal repeats of the small Bacillus bacteriophage genomes. Proc. Natl. Acad. Sci. USA 1981, 78, 2596–2600. [Google Scholar] [CrossRef] [Green Version]
  54. Carstens, E.B. Ratification vote on taxonomic proposals to the International Committee on Taxonomy of Viruses (2009). Arch. Virol. 2009, 155, 133–146. [Google Scholar] [CrossRef] [PubMed]
  55. Paez, J.G. Genome coverage and sequence fidelity of φ29 polymerase-based multiple strand displacement whole genome amplification. Nucleic Acids Res. 2004, 32, e71. [Google Scholar] [CrossRef] [PubMed]
  56. Mencia, M.; Monsalve, M.; Rojo, F.; Salas, M. Transcription activation by phage phi29 protein p4 is mediated by interaction with the alpha subunit of Bacillus subtilis RNA polymerase. Proc. Natl. Acad. Sci. USA 1996, 93, 6616–6620. [Google Scholar] [CrossRef] [Green Version]
  57. Badia, D.; Camacho, A.; Pérez-Lago, L.; Escandón, C.; Salas, M.; Coll, M. The Structure of Phage ϕ29 Transcription Regulator p4-DNA Complex Reveals an N-Hook Motif for DNA Binding. Mol. Cell 2006, 22, 73–81. [Google Scholar] [CrossRef]
  58. Hawley, L.A.; Reilly, B.E.; Hagen, E.W.; Anderson, D.L. Viral Protein Synthesis in Bacteriophage φ29-Infected Bacillus subtilis. J. Virol. 1973, 12, 1149–1159. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Xu, J.; Wang, D.; Gui, M.; Xiang, Y. Structural assembly of the tailed bacteriophage φ29. Nat. Commun. 2019, 10, 2366. [Google Scholar] [CrossRef] [Green Version]
  60. Fu, C.-Y.; Prevelige, P.E. In vitro incorporation of the phage Phi29 connector complex. Virology 2009, 394, 149–153. [Google Scholar] [CrossRef] [Green Version]
  61. E Tosi, M.; E Reilly, B.; Anderson, D.L. Morphogenesis of bacteriophage phi29 of Bacillus subtilis: Cleavage and assembly of the neck appendage protein. J. Virol. 1975, 16, 1282–1295. [Google Scholar] [CrossRef] [Green Version]
  62. Tao, Y.; Olson, N.H.; Xu, W.; Anderson, D.L.; Rossmann, M.G.; Baker, T.S. Assembly of a Tailed Bacterial Virus and Its Genome Release Studied in Three Dimensions. Cell 1998, 95, 431–437. [Google Scholar] [CrossRef] [Green Version]
  63. Steiner, M.; Lubitz, W.; Bläsi, U. The missing link in phage lysis of gram-positive bacteria: Gene 14 of Bacillus subtilis phage phi29 encodes the functional homolog of lambda S protein. J. Bacteriol. 1993, 175, 1038–1042. [Google Scholar] [CrossRef] [Green Version]
  64. Chen, C.; Guo, P. Magnesium-induced conformational change of packaging RNA for procapsid recognition and binding during phage phi29 DNA encapsidation. J. Virol. 1997, 71, 495–500. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Hatfull, G.F.; Jacobs-Sera, D.; Lawrence, J.G.; Pope, W.H.; Russell, D.A.; Ko, C.-C.; Weber, R.J.; Patel, M.C.; Germane, K.; Edgar, R.H.; et al. Comparative Genomic Analysis of 60 Mycobacteriophage Genomes: Genome Clustering, Gene Acquisition, and Gene Size. J. Mol. Biol. 2010, 397, 119–143. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Lima-Mendez, G.; van Helden, J.; Toussaint, A.; Leplae, R. Reticulate Representation of Evolutionary and Functional Relationships between Phage Genomes. Mol. Biol. Evol. 2008, 25, 762–777. [Google Scholar] [CrossRef] [Green Version]
  67. Kotsonis, S.E.; Powell, I.B.; Pillidge, C.J.; Limsowtin, G.K.Y.; Hillier, A.J.; Davidson, B.E. Characterization and Genomic Analysis of Phage Asccφ28, a Phage of the Family Podoviridae Infecting Lactococcus lactis. Appl. Environ. Microbiol. 2008, 74, 3453–3460. [Google Scholar] [CrossRef] [Green Version]
  68. Kleppen, H.P.; Holo, H.; Jeon, S.-R.; Nes, I.F.; Yoon, S.-S. Novel Podoviridae Family Bacteriophage Infecting Weissella cibaria Isolated from Kimchi. Appl. Environ. Microbiol. 2012, 78, 7299–7308. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Young, F.E. Requirement of glucosylated teichoic acid for adsorption of phage in Bacillus subtilis 168. Proc. Natl. Acad. Sci. USA 1967, 58, 2377–2384. [Google Scholar] [CrossRef] [Green Version]
  70. Bhavsar, A.P.; Erdman, L.K.; Schertzer, J.W.; Brown, E.D. Teichoic Acid Is an Essential Polymer in Bacillus subtilis That Is Functionally Distinct from Teichuronic Acid. J. Bacteriol. 2004, 186, 7865–7873. [Google Scholar] [CrossRef] [Green Version]
  71. Kawai, Y.; Marles-Wright, J.; Cleverley, R.M.; Emmins, R.; Ishikawa, S.; Kuwano, M.; Heinz, N.; Bui, N.K.; Hoyland, C.N.; Ogasawara, N.; et al. A widespread family of bacterial cell wall assembly proteins. EMBO J. 2011, 30, 4931–4941. [Google Scholar] [CrossRef] [Green Version]
  72. Hanson, C.A.; Marston, M.F.; Martiny, J. Biogeographic Variation in Host Range Phenotypes and Taxonomic Composition of Marine Cyanophage Isolates. Front. Microbiol. 2016, 7, 983. [Google Scholar] [CrossRef] [Green Version]
  73. Focardi, A.; Ostrowski, M.; Goossen, K.; Brown, M.V.; Paulsen, I. Investigating the Diversity of Marine Bacteriophage in Contrasting Water Masses Associated with the East Australian Current (EAC) System. Viruses 2020, 12, 317. [Google Scholar] [CrossRef] [Green Version]
  74. Sousa, J.A.M.D.; Pfeifer, E.; Touchon, M.; Rocha, E.P.C. Causes and Consequences of Bacteriophage Diversification via Genetic Exchanges across Lifestyles and Bacterial Taxa. Mol. Biol. Evol. 2021, 38, 2497–2512. [Google Scholar] [CrossRef]
  75. Dion, M.B.; Oechslin, F.; Moineau, S. Phage diversity, genomics and phylogeny. Nat. Rev. Genet. 2020, 18, 125–138. [Google Scholar] [CrossRef] [PubMed]
  76. Salifu, S.P.; Rello, A.V.; Campbell, S.A.; Inglis, N.F.; Scortti, M.; Foley, S.; Vázquez-Boland, J.A. Genome and proteome analysis of phage E3 infecting the soil-borne actinomyceteRhodococcus equi. Environ. Microbiol. Rep. 2013, 5, 170–178. [Google Scholar] [CrossRef] [PubMed]
  77. Roucourt, B.; Lavigne, R. The role of interactions between phage and bacterial proteins within the infected cell: A diverse and puzzling interactome. Environ. Microbiol. 2009, 11, 2789–2805. [Google Scholar] [CrossRef] [PubMed]
  78. Kupczok, A.; Neve, H.; Huang, K.D.; Hoeppner, M.P.; Heller, K.J.; Franz, C.M.A.P.; Dagan, T. Rates of Mutation and Recombination in Siphoviridae Phage Genome Evolution over Three Decades. Mol. Biol. Evol. 2018, 35, 1147–1159. [Google Scholar] [CrossRef] [PubMed]
  79. Hua, J.; Huet, A.; Lopez, C.A.; Toropova, K.; Pope, W.H.; Duda, R.L.; Hendrix, R.W.; Conway, J.F. Capsids and Genomes of Jumbo-Sized Bacteriophages Reveal the Evolutionary Reach of the HK97 Fold. mBio 2017, 8, e01579-17. [Google Scholar] [CrossRef] [Green Version]
  80. Reilly, B.E.; Nelson, R.A.; Anderson, D.L. Morphogenesis of Bacteriophage φ29 of Bacillus subtilis: Mapping and Functional Analysis of the Head Fiber Gene. J. Virol. 1977, 24, 363–377. [Google Scholar] [CrossRef] [Green Version]
  81. Xiang, Y.; Rossmann, M.G. Structure of bacteriophage φ29 head fibers has a supercoiled triple repeating helix-turn-helix motif. Proc. Natl. Acad. Sci. USA 2011, 108, 4806–4810. [Google Scholar] [CrossRef] [Green Version]
Viruses 13 01557 g001 550
Figure 1. PumA1 and PumA2 morphologies. Transmission electron micrograph of (a) PumA1 and (b) PumA2. Scale bar = 100 nm.
Figure 1. PumA1 and PumA2 morphologies. Transmission electron micrograph of (a) PumA1 and (b) PumA2. Scale bar = 100 nm.
Viruses 13 01557 g001
Viruses 13 01557 g002 550
Figure 2. Genome maps of PumA1 and PumA2. Colored arrows represent open reading frames (ORFs). Purple arrows show the early region genes (regions I and II), and the green arrows depict the late region genes (region II). Key genes are named and full descriptions of the ORF functions are in Table S1.
Figure 2. Genome maps of PumA1 and PumA2. Colored arrows represent open reading frames (ORFs). Purple arrows show the early region genes (regions I and II), and the green arrows depict the late region genes (region II). Key genes are named and full descriptions of the ORF functions are in Table S1.
Viruses 13 01557 g002
Viruses 13 01557 g003 550
Figure 3. VIRIDIC heatmap of the 34 phi29-like phages. Intergenomic similarities of the phages are shown on the right side with the colored scale where more defined clusters are seen compared to the dot plot. The similarities scores are found in Table S2. The aligned genome fraction and genome length ratio values are shown on the left side with their corresponding scales. Overlayed is the color coding of the seven major clusters, where grey indicates singleton phages.
Figure 3. VIRIDIC heatmap of the 34 phi29-like phages. Intergenomic similarities of the phages are shown on the right side with the colored scale where more defined clusters are seen compared to the dot plot. The similarities scores are found in Table S2. The aligned genome fraction and genome length ratio values are shown on the left side with their corresponding scales. Overlayed is the color coding of the seven major clusters, where grey indicates singleton phages.
Viruses 13 01557 g003
Viruses 13 01557 g004 550
Figure 4. Genome map alignments of each cluster and singletons (shaded grey) with their closest relating phage. Color bar ranges from red to yellow, showing identities, and the blue to green shows homology in any inverted sections.
Figure 4. Genome map alignments of each cluster and singletons (shaded grey) with their closest relating phage. Color bar ranges from red to yellow, showing identities, and the blue to green shows homology in any inverted sections.
Viruses 13 01557 g004
Viruses 13 01557 g005 550
Figure 5. Neighbor-Joining network of the phi29-like phages. Clustal Omega alignments of the whole genomes were calculated and visualized in Splitstree. Clusters are displayed by corresponding colored overlays.
Figure 5. Neighbor-Joining network of the phi29-like phages. Clustal Omega alignments of the whole genomes were calculated and visualized in Splitstree. Clusters are displayed by corresponding colored overlays.
Viruses 13 01557 g005
Viruses 13 01557 g006 550
Figure 6. vConTACT2 reticulate genome network of phages in the Genbank and ICTV databases. A “node” or singular phage genome is connected by an “edge” or line which is calculated by how many proteins they share and is scored by significance. The Salasmaviridae/phi29-like phages cluster was extracted in the dotted box. The network is displayed using an edge-weight spring embedded layout that repels phages based on their lack of similarity in protein cluster scores.
Figure 6. vConTACT2 reticulate genome network of phages in the Genbank and ICTV databases. A “node” or singular phage genome is connected by an “edge” or line which is calculated by how many proteins they share and is scored by significance. The Salasmaviridae/phi29-like phages cluster was extracted in the dotted box. The network is displayed using an edge-weight spring embedded layout that repels phages based on their lack of similarity in protein cluster scores.
Viruses 13 01557 g006
Viruses 13 01557 g007 550
Figure 7. Phylogenetic tree using alignments of both DNA polymerase and encapsidation ATPase as candidate genes. The tree is rooted at midpoint and branches are supported by 100 bootstrap replicates. The colored shadings over each phage name indicate the clustering based off the other genomic comparisons; other metadata is defined in the legend. This clustering and color coding also correlate to the ICTV genus classifications. StevenHerd11 was eliminated from this tree due to its almost identical similarity to RadRaab.
Figure 7. Phylogenetic tree using alignments of both DNA polymerase and encapsidation ATPase as candidate genes. The tree is rooted at midpoint and branches are supported by 100 bootstrap replicates. The colored shadings over each phage name indicate the clustering based off the other genomic comparisons; other metadata is defined in the legend. This clustering and color coding also correlate to the ICTV genus classifications. StevenHerd11 was eliminated from this tree due to its almost identical similarity to RadRaab.
Viruses 13 01557 g007
Table
Table 1. Current Salasmaviridae (phi29-like) phages that have been sequenced and submitted to the Genbank and Bacillus phage database as of February 2021.
Table 1. Current Salasmaviridae (phi29-like) phages that have been sequenced and submitted to the Genbank and Bacillus phage database as of February 2021.
BacteriophageCountry IsolatedAccession NumberYear IsolatedGenome Size (kb)ReferenceHost
PumA1AustraliaMN524844201718,446This studyB. pumilus
PumA2AustraliaMN524845201718,932This studyB. pumilus
MG-B1Austria KC685370201127,190[41]B. weihenstephanensis
VMY22ChinaKT780304201518,609[42]B. cereus
Gxv1ChinaMT459794202021,781[43]Bacillus sp.
DK2ChinaMK284527201823,357[44]B. cereus
DK3ChinaMK284528201826,865[44]B. cereus
DK1ChinaMK284526201827,180[44]B. cereus
DLc1ChinaMW012634202028,950[23]B. cereus
Goe1GermanyKU831549201418,379[45]B. subtilis
Goe6GermanyMF407276201719,105UnpublishedB. velezensis
Goe4GermanyMH817022201825,722[9]B. thuringiensis kurstaki
B103PragueX99260198118,630[46]B. subtilis
SRT01hsRussiaMN857617202020,784UnpublishedB. altitudinis
GA-1ScotlandX96987196521,129[47]B. subtilis
BSTP4South KoreaMW354668202019,145UnpublishedB. subtilis
NfSpainEU622808200818,753[48]B. subtilis
PZAUSAPZACG197619,366[49]B. subtilis
Phi29USAEU771092196519,828[50]B. subtilis
KareziUSAMN082625201320,083UnpublishedB. thuringiensis kurstaki
BeachBumUSAKY921761201621,054[51]B. thuringiensis kutstaki
HarambeUSAKY821088201621,684[51]B. thuringiensis kutstaki
RadRaabUSAMF156580201623,946UnpublishedB. thuringiensis kurstaki
StevenHerd11USAMK084630201723,953UnpublishedB. thuringiensis kurstaki
StitchUSAKX349901201224,320UnpublishedBacillus sp.
JuanUSAMF156577201625,032UnpublishedB. thuringiensis kurstaki
AuroraUSAKX349899201025,908UnpublishedB. thuringiensis kurstaki
QCM11USAKX961631201626,054UnpublishedB. cereus group
KonjoTroubleUSAMF156578201626,061UnpublishedB. thuringiensis kurstaki
ClaudiUSAKX349900201426,504UnpublishedB. thuringiensis kurstaki
SerPounceUSAKY947509201627,206[51]B. thuringiensis kurstaki
WhyPhyUSAMW419775202018,642UnpublishedB. pumilus
ThorntonUSAMW348917201726,319UnpublishedB. thuringiensis kurstaki
Baseball_fieldUSAMT777452201526, 863[52]B. thuringiensis kurstaki
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Stanton, C.R.; Rice, D.T.F.; Beer, M.; Batinovic, S.; Petrovski, S. Isolation and Characterisation of the Bundooravirus Genus and Phylogenetic Investigation of the Salasmaviridae Bacteriophages. Viruses 2021, 13, 1557. https://0-doi-org.brum.beds.ac.uk/10.3390/v13081557

AMA Style

Stanton CR, Rice DTF, Beer M, Batinovic S, Petrovski S. Isolation and Characterisation of the Bundooravirus Genus and Phylogenetic Investigation of the Salasmaviridae Bacteriophages. Viruses. 2021; 13(8):1557. https://0-doi-org.brum.beds.ac.uk/10.3390/v13081557

Chicago/Turabian Style

Stanton, Cassandra R., Daniel T. F. Rice, Michael Beer, Steven Batinovic, and Steve Petrovski. 2021. "Isolation and Characterisation of the Bundooravirus Genus and Phylogenetic Investigation of the Salasmaviridae Bacteriophages" Viruses 13, no. 8: 1557. https://0-doi-org.brum.beds.ac.uk/10.3390/v13081557

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

Article Metrics

Back to TopTop