1. Introduction
Worldwide, non-typhoidal Salmonella (NTS) enterica subspecies enterica serovars Typhimurium and Enteritidis (S. Typhimurium and S. Enteritidis) cause self-limiting enterocolitis in humans. However, in recent years novel lineages of NTS serovars Typhimurium and Enteritidis have been associated with increasing levels of invasive non-typhoidal Salmonella (iNTS) disease, mostly in Sub-Saharan Africa [1]. These bloodstream infections occur in immunocompromised individuals, resulting in about 50,000 deaths each year in Africa [1]. These lineages are characterised by genomic degradation, a distinct prophage repertoire, and multidrug resistance [2,3]. S. Typhimurium sequence type (ST) 313 has been widely studied [4,5] due to its association with iNTS disease, and D23580 is the representative strain for Lineage 2. Although S. Enteritidis is the other main serovar responsible for iNTS in Africa, this serovar has been primarily studied in relation to the implications in food industry. Several distinct and geographically linked lineages have been identified that cause iNTS (the West African, and Central and Eastern African lineages) or the global enterocolitis epidemic [3].
Phage-resistance mechanisms and other frequently carried genes involved in bacterial metabolism and pathogenesis increase bacterial fitness and support adaptation to different environments [6,7]. Prophages found in iNTS have been described in S. Typhimurium [8,9], and partially described in S. Enteritidis [3], for example Salmonella prophages Gifsy-1, Gifsy-2, BTP1, and Fels2-BT. Prophages play an important role in Salmonella virulence, increasing fitness and preventing infection by other phages [6,10,11]. Indeed, a novel phage exclusion system was identified encoded in S. Typhimurium D23580 prophage BTP1 [10], supporting the host to survive phage predation in the environment and drive niche specialisation. There are diverse anti-phage systems widely spread in bacteria [12,13], stressing the need to find phages that are capable of escaping these anti-phage mechanisms that protect African iNTS.
Moreover, the increase in antimicrobial resistant strains is driving the search for alternatives to widely used antibiotics. Here, we aimed to isolate and characterise phages capable of lysing the novel lineages of S. Typhimurium and S. Enteritidis associated with bloodstream infection. We sought phages on two continents, in both the UK and Malawi, with the objective of isolating a pool of phages with a different host range. 2. Materials and Methods 2.1. Cloning Procedures
All the bacterial strains, plasmids, and bacteriophages used in this study are described in Table 1, and oligonucleotides (primers) are described in Table S1. DNA manipulations were carried out according to standard protocols [14] and to enzyme/kit manufacturer’s recommendations. Restriction enzymes and T4 DNA ligase were purchased from Thermo Fisher Scientific or New England Biolabs (NEB) and DNA purification and plasmid isolation kits were obtained from Bioline (Cat#BIO-52060 and Cat#BIO-52057, respectively).
For the construction of pEMG-based suicide plasmids [26], the DNA inserts of interest were digested with EcoRI and BamHI and ligated into pEMG, using the corresponding sites.
For the construction of pNAW1, two DNA fragments (~750 bp each) flanking the Gifsy-2 attachment site (attBGifsy-2 4/74) of S. Typhimurium strain 4/74 were PCR amplified with primer pairs NW_1/NW_2 and NW_3/NW_4. The resulting fragments were fused together by overlap extension PCR [30], as primers NW_2 and NW_3 share a complementary stretch of 25 nucleotides. Finally, the resulting fragment (1529 bp) was digested and ligated into pEMG. Similarly, a suicide plasmid was constructed to knock out the galE gene (pNAW34): the flanking regions of galE were amplified with primer pairs NW_82/NW_83 and NW_84/NW_85 and the resulting fragments (~450 bp each) were fused by overlap extension PCR, digested and ligated into pEMG.
S. Typhimurium strain 4/74 carries the SopEΦ prophage that is flanked by a P4-like remnant prophage (hereafter called P4-like) [31]. Those prophage elements are absent in strain 14028s that was used to amplify the corresponding empty attB site with primers NW_65 and NW_65. After digestion, the resulting fragment (1531 bp) was ligated into pEMG. All the pEMG-derived plasmids were checked by Sanger sequencing using primers M13_-40_long and M13_rev_long and the Lightrun sequencing service (Eurofins Genomics).
2.2. Construction of Mutant Hosts
The anti-phage brex operon [32] (genes STMMW_44361-STMMW_44431, GenBank accession number FN424405.1) was deleted in the prophage-free strain S. Typhimurium D23580 ΔΦ (JH3949) by λ red recombination [27]. Electro-competent D23580 ΔΦ carrying the temperature inducible λ red plasmid pSIM5-tet [28] were prepared, as previously described [33]. Competent bacteria were electroporated with a PCR fragment obtained by PCR with primers Brex_del_Fw and Brex_del_Rv and pKD4 as template. After selection on kanamycin (Km) medium, recombinant D23580 ΔΦ ∆brex::Km mutants were obtained. The mutation was transferred into D23580 ΔΦ by transduction using the P22 HT 105/1 int-201 transducing phage [8,16]. Finally, the Km R cassette was removed with the flippase expressing plasmid pCP20-Gm [29], resulting in strain D23580 ΔΦ ∆brex::frt (JH4314). Finally, the mutation was confirmed by Sanger sequencing using the primers Brex_up and Brex_down.
For the scarless excision of prophages and for the deletion of galE, the technique based on the pEMG suicide was used [26]. Suicide plasmids carrying the appropriate empty attB or the ∆galE fragment were mobilised from Escherichia coli S17-1 λpir into S. Typhimurium by conjugation and the resulting Km resistant meroploids were resolved with the pSW-2 plasmid, as described previously [8,27,34]. After successive prophage deletions, the prophage-free strain S. Typhimurium 4/74 ΔΦ (JH4180) was obtained: this strain lacked the Gifsy-1, Gifsy-2, ST64B, SopEΦ and P4-like prophage elements and was sequenced using the Illumina whole-genome sequencing service at MicrobesNG (Birmingham, UK).
2.3. Phage Isolation
Water samples were enriched with either S. Typhimurium or S. Enteritidis hosts (Table 2) for initial bacteriophage isolation. Further propagation and storage of isolated phages was carried out by using a susceptible Salmonella host (S. Typhimurium D23580 ∆Φ ∆brex) constructed from the iNTS strain D23580 [4] with all prophages removed [8], lacking in addition the BREX bacteriophage exclusion defence system [32].
Several locations were selected in the UK and Malawi for field sampling. The sampling locations in Blantyre, Malawi were around Mudi river [35], Naperi river, drains near food markets and Queen Elizabeth Central Hospital. The sampling locations in the UK were pond and river water near parks (Sefton Park, Liverpool, and Botanic Gardens, Southport) and the outflow from sewage treatment plants (River Alt, Liverpool, and River Wear, Durham).
Briefly, we collected 50 mL water samples that were first centrifuged at 7000 rpm (Eppendorf Centrifuge 5430R) at room temperature for 5 min. The enrichments were prepared by filtering 10 mL supernatant through 0.22 µm pore-size membranes (Star Lab Cat#E4780-1226) into 1:1 v/v LB 2X (2% tryptone, 1% yeast extract, 1% NaCl; pH 7.0) inoculated with the Salmonella host (1:20 v/v) grown overnight, and incubating with agitation (200 rpm) for 3 days at 37 °C. Then, 5 mL aliquots from each enrichment were centrifuged at 7000 rpm (Eppendorf Centrifuge 5430R) at room temperature for 15 min, and the supernatant filtered through 0.22 µm pore-size membranes. A 10-fold serial dilution was prepared using sterile LB for double-layer plaque assay [10,36] on either 90 mm or 12 × 12 cm Petri dishes, incubated overnight at 37 °C. Individual plaques with different morphologies were selected and purified repeating the serial dilution and plaque assay process twice.
Phage titration and storage were carried out using the susceptible Salmonella Typhimurium strain D23580 ∆Φ ∆brex. For storage, plaques were selected from plates and the phages resuspended in sterile LB with 1% chloroform, followed by vigorous shaking and centrifugation at 14,000 rpm (Eppendorf Centrifuge 5424) for 5 min at room temperature, and stored at 4 °C in the dark. 2.4. Strain Selection and Phage Host Range
The S. Typhimurium D23580 strain was used for phage isolation because it is a representative of ST313 lineage 2, the Salmonella variant responsible for iNTS epidemics in Malawi. Other relevant Salmonella strains were also included to expand the likelihood of phage isolation. The S. Enteritidis strains selected were representatives of four of the major clades, namely the Global Epidemic, Central/Eastern, and West African clades [3].
To determine the host range of the isolated bacteriophages, a spot test [37] was performed using phage stocks that contained between 1010 and 1011 pfu/mL. Briefly, bacterial lawns of different Salmonella hosts (Table 1) were prepared by diluting overnight bacterial cultures 40 folds with top-agar (0.5% agar in LB, 50 °C), and pouring over 1.5% LB-agar plates (4 mL top agar for 90 mm diameter or 8 mL for 12 × 12 cm Petri dishes). The plates were incubated for 24 h at 37 °C, and lysis formation recorded (Figure S1). This assay was also performed using characterised published (reference) phages (Table 1) and using between two and six biological replicates.
2.5. RAPD PCR
Randomly amplified polymorphic DNA polymerase chain reaction (RAPD PCR) was performed using the phage lysate as template. The phage lysate was prepared with 85 µL phage stock (1010–1011 pfu/mL), treated with 10 µL 10XDNase buffer (100 mM Tris-HCl, 25 mM MgCl2, pH 7.5) and 5 µL 10 mg/mL DNase I (SIGMA DN25), incubated at 37 °C shaking (300 rpm) for 60 min, followed by a second incubation at 75 °C for 10 min. Four oligonucleotides were designed [38] for PCR using MyTaq Red Mix 1X (Bioline, Cat#BIO-25043), according to the manufacturer’s recommendations. The oligonucleotides used were OPL5, RAPD5, P1 and P2 (Table S1), as previously described [38]. The PCR products were analysed by 0.8% agarose-gel electrophoresis with 1 kb ladder (Bioline, Cat#H1-819101A), run at 100 V for 60 min, and visualised with Midori Green DNA/RNA staining (Nippon Genetics, Cat#MG06).
2.6. Phage DNA Extraction and Sequencing Phage DNA was extracted from 410 µL of phage stock (5 × 1010–1011 pfu/mL) with 50 µL 10X DNase buffer (100 mM Tris-HCl, 25 mM MgCl2, pH 7.5), 20 µL 10 mg/mL DNase I (SIGMA DN25), and 20 µL of 1 mg/mL RNase A (SIGMA), incubated at 37 °C shaking (300 rpm) for 60 min, followed by a second incubation at 75 °C for 10 min. The lysate was treated with 10 µL 20 mg/mL Proteinase K (BIOLINE, Cat#BIO-37037) and purified using the Norgen Phage DNA Isolation Kit (cat. 46850) following manufacturer’s instructions. The purified DNA was eluted in 75 µL molecular grade water. DNA integrity was assessed by 1% agarose-gel electrophoresis, and DNA concentration was determined fluorometrically with dsDNA HS Assay kit (Qubit, ThermoFisher Scientific).
Phage DNA sequencing was carried out at Nu-Omics DNA sequencing research facility (Northumbria University, UK) using MiSeq (Illumina). Generated reads were quality checked using FastQC v0.11.5 (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/, accessed on 9 May 2019), and quality trimmed using SEQTK v1.3-r106 (https://github.com/lh3/seqtk, accessed on 9 May 2019). Phage genomes were assembled using Unicycler v0.3.0b [39] with default parameters, and the quality of the assemblies were assessed with QUAST v4.6.3 [40].
Genomes that had more than one contiguous sequence were visualised in Artemis v10.2 [41], and circularised contigs were extracted for viral species analysis using BLAST [42]. Phages BPS7 and ER24 had multiple linear contigs, therefore BLAST was used to identify the viral contig for extraction and analysis. Phage assemblies were annotated using Prokka v1.12 [43] with--proteins/blastdb/Viral_Genomes/all_viral.faa, visualised using SnapGene v5.2 (Insightful Science, snapgene.com) and manually curated (Tables S2–S5). PhageAI [44] is a tool used to determine lifecycle of phages with high accuracy. Phage genome assemblies were analysed using PhageAI v0.1.0 LifeCycleClassifier using default parameters.
All genomes were uploaded to GenBank with phage names assigned based on Adriaenssens and Brister (2017) [45]; however, common names were used for simplicity (Table 2).
2.7. Comparative Genomics and Phylogenetic Analyses
Gene sequence alignment of terminase large subunit (terL) was performed by ClustalO v1.2.4 [46] using default parameters, and a maximum likelihood tree was constructed using RAxML v8.2.8 [47], visualised in the interactive Tree of Life (iTOL) v4.2 [48] rooted at midpoint. Genomic similarities were calculated by pairwise nucleotide comparison with the Virus Intergenomic Distance Calculator (VIRIDIC) [49] with default parameters.
Genomes of phages belonging to the same cluster were analysed. A gene alignment was performed using Roary v3.11.0 [50] with default parameters for each cluster. A maximum likelihood phylogenetic tree was built using RAxML-NG v0.4.1 beta [47] with 100 bootstrap replicates to assess support, and visualised in iTOL v4.2 unrooted. Genome comparisons were performed using nucleotide alignment with MEGA v6.06 [51] and BLASTn [42], and visualised with Brig v0.95 [52].
2.8. Transmission Electron Microscopy (TEM)
Morphology of selected phages was determined by TEM [8], using 10 µL of phage stock (1010–1011 pfu/mL) pipetted onto carbon-coated copper 200-mesh grid (TAAB F077/050), and incubated at room temperature for 10 min. The grids were washed twice in distilled water for 2 min. To negative-stain the phages, 10 µL of 2% uranyl-acetate was added and incubated for 1 min. Grids were allowed to dry and TEM was performed on EI 120 kV Tecnai G2 Spirit BioTWIN transmission electron microscope.
3. Results and Discussion 3.1. Phage Isolation
To isolate a diverse pool of phages, we used water samples collected in the UK and Malawi, close to areas where Salmonella or other related bacteria would be expected to coexist, such as sewage plants or outlets. We enriched for phages that plated on S. Typhimurium or S. Enteritidis hosts (Table 2), and we selected plaques with different morphologies. The RAPD PCR assay was used to differentiate the phages, and 32 distinct phages were selected for whole genome sequencing (Figure S2). The 32 bacteriophages included 20 phages isolated from Malawi and 12 phages from the UK.
The majority of the phages (25 phages, 78%) were isolated using the phage-susceptible S. Typhimurium D23580 ΔΦ Δbrex mutant as a host. This mutant was constructed by removing all known ‘anti-phage’ encoding genetic features, which include all five prophages [8] and BREX bacteriophage exclusion defence system [32], from a well-characterised S. Typhimurium ST313 associated bloodstream infection in Malawi [4].
The remaining seven phages were isolated with several other host strains, including the S. Typhimurium ST19 representative strain 4/74, S. Typhimurium D23580 ΔΦ [8] (with the BREX system intact), and S. Enteritidis strains.
3.2. Host Range Analysis
The host range of all 32 phages was tested using bacterial strains representing S. Typhimurium ST19 and ST313, and S. Enteritidis Global Epidemic and African clades (Figure 1). No significant difference in host range was observed between phages isolated from the UK and those from Malawi (Figure 1). A small proportion of phages (9 phages, 28%) that originated from UK or Malawi were able to infect S. Enteritidis from the African clades, of which only two (BPS5 and BPS6) were isolated in Liverpool with S. Enteritidis as the host. BPS5 was the only phage capable of infecting all hosts tested, including mutants. The phages BPS3, BPS6 and BPS7 were able to infect all wild-type Salmonella strains tested. These results indicate that phages isolated in the UK had a broader host range than phages isolated from Malawi.
Of the 32 phages isolated, seven were unable to infect S. Enteritidis. The Global Epidemic S. Enteritidis strains P125109 and A1636 were susceptible to infection by the same phages, with the exception of phage ER8 that only infected A1636 (Figure 1). This result reflects the fact that these S. Enteritidis strains were very similar at the genomic level, differing by ~40 SNPs. S. Enteritidis from the African clades were resistant to infection of most (>70%) of the 32 phages. S. Enteritidis strain 4030/15, a representative from the West African clade, was susceptible to infection of the same phages as D7795 and CP255, representatives of the Central/Eastern African clade, with the addition of phages BPS1, ER2, ER4, ER5 and ER10.
S. Enteritidis strains, especially from the African clades, are also resistant to being infected by other commonly used bacteriophages, such as the Salmonella phages 9NA, BTP1 and P22 (Figure 1). The P22 transducing phage is commonly used for genetic manipulation of Salmonella. However, P22 infects S. Enteritidis from the Global Epidemic clade at a lower level than S. Typhimurium, and did not cause clearing of infected S. Enteritidis cultures due to an unknown mechanism. Interestingly, Salmonella phage Det7 was capable of infecting both S. Typhimurium and S. Enteritidis hosts, sharing a similar host range profile to phages BPS3, BPS5, BPS6, and BPS7.
More phages infected S. Typhimurium ST19 (Figure 1) strains 14028s (32 phages) and LT2 (29 phages), than S. Typhimurium 4/74 (eight phages), possibly due to prophage-encoded “anti-phage” defence mechanisms. Host range was tested on a S. Typhimurium 4/74 mutant that lacked prophages (4/74 ΔΦ), and the number of infecting phages increased from eight to 13 (Figure 1).
All isolated phages, except for ER21, were able to infect S. Typhimurium D23580 (Figure 1). However, in comparison to the original host D23580 strains that were deleted for either prophages (D23580 ΔΦ) or prophages and BREX (D23580 ΔΦΔbrex), some plaque sizes were smaller on D23580 wild type (Figure 1). This suggested phages ER22, ER23, ER25, BPS3, BPS6 and BPS7 are susceptible to a resistance mechanism encoded by one of the D23580 prophages. BREX did not appear to impact these phages (Figure 1).
Phage BTP1 is one of the key prophages carried by S. Typhimurium ST313 that has very high spontaneous induction (~109 pfu/mL) [8]. BTP1 is active and highly conserved in both the lineages 1 and 2 of ST313 that are responsible for about two-thirds of the iNTS cases in Sub-Saharan Africa [53]. Here, we used a S. Typhimurium D23580 mutant (D23580 ΔBTP1) to show that BTP1 conferred resistance against phage ER21, possibly due to the removal of the novel BstA defence system [10] or the GtrAC LPS-modification enzymes [11].
Cell-surface properties frequently mediate phage-resistance by blocking attachment of tailed-bacteriophages, either by modification of the lipopolysaccharide (LPS), outer membrane proteins or flagella [54]. The LPS oligo-polysaccharide residues (O-antigen) are common phage receptors [55]. To determine whether the O-antigen was required for phage infection, a S. Typhimurium D23580 mutant lacking all prophages, the BREX system, and galE (D23580 ΔΦ Δbrex ΔgalE) was tested. The gene galE encodes the UDP-glucose 4-epimerase component of the galactose biosynthetic pathway, and is required for the synthesis of the LPS O-antigen of Salmonella. Phages ER25, BPS3, BPS6, BPS7, ER6 and ER18 were unable to infect the mutant D23580 ΔΦ Δbrex ΔgalE (Figure 1), suggesting that LPS is the receptor for attachment of these bacteriophages.
3.3. Comparative Genomics and Phylogeny of Isolated Phages
We used whole genome sequencing to investigate the taxonomy and relatedness of all the Malawian and UK phages. Pairwise genome comparison of nucleotide similarity between the 32 isolated phages revealed three main clusters (Figure S3). The relatedness of the bacteriophages was also assessed by alignment of the gene encoding the terminase large sub-unit (terL) that was present in all phages (Figure 2).
Cluster 1 contained phages isolated in the UK that were able to infect all wild-type Salmonella strains tested. Cluster 2 included a mixture of bacteriophages isolated in the UK and Malawi that were only capable of infecting S. Typhimurium (except phage BPS1). Cluster 3 was divided into two sub-clusters: 3.a contained phages isolated from the UK, and 3.b from Malawi. Generally, the Cluster 3 phages infected S. Typhimurium 14028s, LT2 and D23580, and S. Enteritidis from the Global Epidemic clade, with a few of the phages infecting the representative West African clade strain 4030/15. Phage ER25 was the only one that did not belong to Cluster 1, 2 or 3. Phages belonging to Cluster 1 had a relatively large genome of ~158 kb, compared to the ~60 kb genomes of Cluster 2 phages, and ~43 kb genomes of Cluster 3 phages.
Comparative genomic analysis of Cluster 1 phages (Figure 3A) showed that BPS3, BPS5 and BPS6 were closely related to S117 (>98% nucleotide identity; GenBank MH370370.1), a phage originally isolated on S. Typhimurium LT2 [56]. Figure 3B shows that most phage-encoded genes are conserved throughout Cluster 1. BPS3 lacked one of the tail-spike genes present in BPS6, BPS5, BPS7, and reference phage S117, a gene that is upstream of the predicted virulence protein VriC, also found in other Salmonella phages [57]; however, it is not clear whether VriC is a functional effector protein. Other genes were missing from Cluster 1 phages such as BPS7, that encode hypothetical proteins preventing any inference of their function. The comparative genomic analysis identified 26 SNP differences between BPS5 and BPS6, all of which are in intergenic regions upstream hypothetical proteins and tRNA genes. None of these nucleotide differences explained why BPS5 is capable of infecting all hosts tested, including the mutant S. Typhimurium D23580 strain with a short LPS (D23580 ΔΦ Δbrex ΔgalE), but BPS6 could not. Interestingly, BPS3 is the only phage in Cluster 1 with a different tail-spike protein (tail-spike 1).
Cluster 2 (Figure 4A) included phages isolated from Malawi and the UK that were closely related (>93% identity) to S. Typhimurium phages iEPS5 (GenBank KC677662.1) [58] and øχ (Chi, GenBank NC_025442.1) [59], which uses the flagellum as the receptor, and phage 37 [60] isolated from sewage in India. Based on the comparative genomic analysis (Figure 4B), seven structural genes including a tape measure protein, head vertex, and a portal protein were identified. Additionally, genes encoding putative endolysins proteins A and B and terminase large and small subunits (terSL) were identified.
Cluster 3 was the largest and most diverse cluster, with different host range and genome composition. Genome comparison (Figure 5A) indicated a clear division between sub-clusters 3.a and 3.b, with sub-cluster 3.b only containing phages isolated from Malawi. Phages ER21, ER22 and ER23 (sub-cluster 3.a) (Figure 5B) were closely related (>96% identity) to phage S134 (GenBank MH370381) [56], originally isolated using S. Enteritidis PT1 strain E2331, which might indicate their ability to infect S. Enteritidis from the Global Epidemic clade. Sub-cluster 3.b included very closely related phages, with differences only in two hypothetical proteins (Figure 5C). Several structural genes were identified in phages from sub-cluster 3.a, however no structural genes were found in sub-cluster 3.b isolates, highlighting the novelty of these phages. Genes related to replication, like DNA polymerases and topoisomerases, and recombination were identified in phages from sub-cluster 3.b. Additionally, the lysozyme encoded by gene rrrD was identified; however, most genes were identified as hypothetical proteins.
Finally, phage ER25 did not belong to any of the identified clusters. Comparative genomic analysis demonstrated that phage ER25 is very similar (>99% identity) to P22 (GenBank AF217253.1) [61], differing by only 49 SNPs (Figure S4). This finding, complemented by analysis with PhageAI, suggests that ER25 was the only temperate phage identified, whereas all other isolated phages were likely to have a virulent life cycle.
Four bacteriophages were selected as representatives of each Cluster for morphological analysis by transmission electron microscopy. All phages were from the order Caudovirales, with two main morphotypes representing different phage families. Phages from Clusters 2 and 3 had a long (>100 nm) non-contractile tail and belonged to the Siphoviridae family, whereas the phages from Cluster 1 had a contractile tail, belonging to the Ackermannviridae family, confirmed by genome comparison with phage Det7 that also groups in Cluster 1. 3.4. Perspective
This study represents the first isolation of virulent phages that target Salmonella variants associated with bloodstream infection in Sub-Saharan Africa. The phages target multidrug resistant African lineages of iNTS that are particularly difficult to control by many antimicrobials and are not susceptible to well-characterised phages. Our data also demonstrate that prophage repertoire is critical for phage resistance in Salmonella, suggesting that novel anti-phage systems are likely to be discovered in the genomes of clinically relevant Salmonella. We have identified phages with different host ranges, including generalist and specialist phages against S. Typhimurium and S. Enteritidis. This phage collection has the potential for wider applications, such as decontamination of food, swine slurry [62], and hospital or food preparation surfaces [63], which is particularly relevant in low and middle-income countries [64], and could have utility against multidrug resistant iNTS infections in clinical settings.
Enlarge this image.Figure 1. Host range of isolated phages. Host range analysis was determined by spot assay on S. Enteritidis or S. Typhimurium host strains, with phages isolated in the UK (blue) or Malawi (red). Reference phages are displayed in black. Host range of biological replicates was recorded as defined lysis (dark green), turbid edge or smaller plaques (light green), and no lysis (grey). Examples of the different lysis observed are shown in Figure S1. The original Salmonella host used for phage isolation is indicated in purple. NT: not tested.
Enlarge this image.Figure 2. Phylogenetic and phenotypic characterisation of related isolated phages. Gene alignment of terL gene visualised in iTOL, with numbers indicating bootstraps on branches leading to clusters. The three clusters of related phages identified are presented in lilac (Cluster 1), green (Cluster 2), light blue (Cluster 3.a) and dark blue (Cluster 3.b). The host range for S. Typhimurium and S. Enteritidis is presented as a filled circle (lysis) or a clear circle (no lysis) from Figure 1, genome size (kb), and morphology observed by transmission electron microscopy (TEM) with scale bars of 100 nm. Inferred taxonomy is indicated next to TEM images. Phage names in blue were isolated from the UK and red from Malawi. Names in black indicate reference phages, and stars indicate representative phages used for TEM.
Enlarge this image.Figure 3. Comparative genomics of phages from Cluster 1. (A) Maximum-likelihood tree of phages in Cluster 1, unrooted and visualised in iTOL with numbers on branches indicating bootstraps. Scale bar represents number of single nucleotide polymorphisms. Phage names in blue were isolated from the UK and black font indicates reference phages. (B) Whole-genome comparison of phages in Cluster 1. The local similarity of each genome at the nucleotide level was calculated based on BLASTn high scoring pairs and plotted against a circular map of BPS6 represented as the inner circle. Blue circles represent reference genomes. Genome annotation is shown as arrows in the inner circle, colour-coded based on gene categories. Full annotation can be found in Table S2.
Enlarge this image.Figure 4. Comparative genomics of phages from Cluster 2. (A) Maximum-likelihood tree of phages in Cluster 2, unrooted and visualised in iTOL with numbers on branches indicating bootstraps. Scale bar represents number of single nucleotide polymorphisms. Phage names in blue were isolated from the UK and red from Malawi. Names in black font indicate reference phages. (B) Whole-genome comparison of phages in Cluster 2. The local similarity of each genome at the nucleotide level was calculated based on BLASTn high scoring pairs and plotted against a circular map of ER6 represented as the inner circle. Blue circles represent reference genomes. Genome annotation is shown as arrows in the inner circle, colour-coded based on gene categories. Full annotation can be found in Table S3.
Enlarge this image.Figure 5. Comparative genomics of phages from Cluster 3. (A) Maximum-likelihood tree of phages in Cluster 3, unrooted and visualised in iTOL with numbers on branches indicating bootstraps. Scale bar represents number of SNPs. Phage names in blue were isolated from the UK and red from Malawi. Names in black font indicate reference phages. (B) Whole-genome comparison of phages in Cluster 3.a. The local similarity of each phage genome at the nucleotide level was calculated based on BLASTn high scoring pairs and plotted against a circular map of ER22, represented as the inner circle. Blue circles represent reference genomes. Genome annotation is shown as arrows in the inner circle, colour-coded based on gene categories. Full annotation can be found in Table S4. (C) Whole-genome comparison of phages in Cluster 3.b. The local similarity of each genome was calculated based on BLASTn high scoring pairs and plotted against a circular map of ER20 represented as the inner circle. Blue circles represent reference genomes. Genome annotation is shown as arrows in the inner circle, colour-coded based on gene categories. Full annotation can be found in Table S5.
| Strain | Description a | Reference |
|---|---|---|
| Bacteriophages | ||
| P22 | Podoviridae, temperate, wild type from strain LT2 | [8,15] |
| P22 HT 105/1 int-201 | Podoviridae, virulent transducing phage | [16] |
| BTP1 | Podoviridae, temperate, wild type from strain D23580 | [8] |
| Det7 | Ackermannviridae, virulent, wild type | [17] |
| 9NA | Siphoviridae, virulent, wild type | [18] |
| S. Enteritidis | ||
| P125109 | Wild type | [19] |
| A1636 | Wild type | [3] |
| 4030_15 | Wild type | [20] |
| D7795 | Wild type | [3] |
| CP255 | Wild type | [21] |
| S. Typhimurium | ||
| 14028s | Wild type | [22] |
| LT2 | Wild type | [23] |
| 4/74 derivatives | ||
| 4/74 | Wild type | [24] |
| JH4180 | 4/74 ΔΦ (ΔGifsy-1 ΔGifsy-2 ΔST64B ΔSopEΦ ΔP4-like) | This study |
| D23580 derivatives | ||
| D23580 | Wild type | [4] |
| JH3877 | D23580 ΔBTP1 | [8] |
| JH3949 | D23580 ΔΦ (ΔBTP1 ΔBTP5 ΔGifsy-1 ΔGifsy-2 ΔST64B) | [8] |
| JH4314 | D23580ΔΦ Δbrex::frt (ΔSTMMW_44361- STMMW_44431) | This study |
| JH4655 | D23580ΔΦ Δbrex::frt ΔgalE | This study |
| Escherichia coli | ||
| S17-1 λpir | pro thi hsdR recA chromosome::RP4-2 Tc::Mu Km::Tn7/λpir; Tp R, Sm R | [25] |
| Plasmids | ||
| pEMG | Suicide plasmid; Km R | [26] |
| pNAW1 | pEMG::attBGifsy-2 4/74; Km R | This study |
| pNAW15 | pEMG::attBGifsy-1; Km R | [8] |
| pNAW19 | pEMG::attBSopEΦ&P4-like; Km R | This study |
| pNAW34 | pEMG::∆galE; Km R | This study |
| pNAW42 | pEMG::attBST64-B; Km R | [8] |
| pSW-2 | I-SceI expressing vector, m-toluate inducible; Gm R | [26] |
| pKD4 | frt-aph-frt (Km R cassette flanked with frt) template plasmid; Km R | [27] |
| pSIM5-tet | λ Red recombination plasmid, temperature-inducible; Tc R | [28] |
| pCP20-Gm | Flippase recombinase expression plasmid; Gm R | [29] |
a Relevant antibiotic resistances are indicated by R: Gm, gentamicin; Km, kanamycin; Sm, streptomycin; Tc, tetracycline; Tp, trimethoprim. Natural resistances of wild type strains are not indicated. The flippase recognition target sequence from pKD4 is indicated by “frt”.
| Common Name | Phage Name | Place of Isolation | Original Host | Cluster a | GenBank Accession Number |
|---|---|---|---|---|---|
| ER1 | vB_SenS_ER1 | Blantyre, Malawi | S. Typhimurium D23580ΔΦbrex | 3.b | MW355461 |
| ER2 | vB_SenS_ER2 | Blantyre, Malawi | S. Typhimurium D23580ΔΦbrex | 3.b | MW355466 |
| ER3 | vB_SenS_ER3 | Blantyre, Malawi | S. Typhimurium D23580ΔΦbrex | 2 | MW355467 |
| ER4 | vB_SenS_ER4 | Blantyre, Malawi | S. Typhimurium D23580ΔΦbrex | 3.b | MW355468 |
| ER5 | vB_SenS_ER5 | Blantyre, Malawi | S. Typhimurium D23580ΔΦbrex | 3.b | MW355469 |
| ER6 | vB_SenS_ER6 | Blantyre, Malawi | S. Typhimurium D23580ΔΦbrex | 2 | MW355470 |
| ER7 | vB_SenS_ER7 | Blantyre, Malawi | S. Typhimurium D23580ΔΦbrex | 3.b | MW355471 |
| ER8 | vB_SenS_ER8 | Blantyre, Malawi | S. Typhimurium D23580ΔΦbrex | 3.b | MW355472 |
| ER9 | vB_SenS_ER9 | Blantyre, Malawi | S. Typhimurium D23580ΔΦbrex | 3.b | MW355473 |
| ER10 | vB_SenS_ER10 | Blantyre, Malawi | S. Typhimurium D23580ΔΦbrex | 3.b | MW355451 |
| ER11 | vB_SenS_ER11 | Blantyre, Malawi | S. Typhimurium D23580ΔΦbrex | 3.b | MW355452 |
| ER12 | vB_SenS_ER12 | Blantyre, Malawi | S. Typhimurium D23580ΔΦbrex | 3.b | MW355453 |
| ER13 | vB_SenS_ER13 | Blantyre, Malawi | S. Typhimurium D23580ΔΦbrex | 3.b | MW355454 |
| ER14 | vB_SenS_ER14 | Blantyre, Malawi | S. Typhimurium D23580ΔΦbrex | 3.b | MW355455 |
| ER15 | vB_SenS_ER15 | Blantyre, Malawi | S. Typhimurium D23580ΔΦbrex | 3.b | MW355456 |
| ER16 | vB_SenS_ER16 | Blantyre, Malawi | S. Typhimurium D23580ΔΦbrex | 3.b | MW355457 |
| ER17 | vB_SenS_ER17 | Blantyre, Malawi | S. Typhimurium D23580ΔΦbrex | 3.b | MW355458 |
| ER18 | vB_SenS_ER18 | Blantyre, Malawi | S. Typhimurium D23580ΔΦbrex | 3.b | MW355459 |
| ER19 | vB_SenS_ER19 | Blantyre, Malawi | S. Typhimurium D23580ΔΦbrex | 2 | MW355460 |
| ER20 | vB_SenS_ER20 | Blantyre, Malawi | S. Typhimurium D23580ΔΦbrex | 3.b | MW355462 |
| ER21 | vB_SenS_ER21 | Liverpool, UK | S. Typhimurium D23580ΔΦbrex | 3.a | MW355463 |
| ER22 | vB_SenS_ER22 | Liverpool, UK | S. Typhimurium D23580ΔΦbrex | 3.a | MW355464 |
| ER23 | vB_SenS_ER23 | Southport, UK | S. Typhimurium D23580ΔΦbrex | 3.a | MW355465 |
| ER24 | vB_SenS_ER24 | Durham, UK | S. Typhimurium D23580ΔΦbrex | 2 | MW355479 |
| ER25 | vB_SenP_ER25 | Durham, UK | S. Typhimurium D23580ΔΦbrex | ND | MW355480 |
| BPS1 | vB_SenS_BPS1 | Liverpool, UK | S. Typhimurium D23580ΔΦ | 2 | MW355474 |
| BPS2 | vB_SenS_BPS2 | Liverpool, UK | S. Typhimurium D23580ΔΦ | 2 | MW355449 |
| BPS3 | vB_SenAc_BPS3 | Liverpool, UK | S. Typhimurium 4/74 | 1 | MW355475 |
| BPS4 | vB_SenS_BPS4 | Liverpool, UK | S. Typhimurium D23580ΔΦ | 2 | MW355476 |
| BPS5 | vB_SenAc_BPS5 | Liverpool, UK | S. Enteritidis D7795 | 1 | MW355477 |
| BPS6 | vB_SenAc_BPS6 | Liverpool, UK | S. Enteritidis D7795 | 1 | MW355450 |
| BPS7 | vB_SenAc_BPS7 | Liverpool, UK | S. Enteritidis P125109 | 1 | MW355478 |
a Clusters determined based on whole genome alignment. ND: not determined.
Supplementary Materials
The following are available online at https://www.mdpi.com/1999-4915/13/3/478/s1, Figure S1: Plaque morphology used for defining host range, Figure S2: Random amplified polymorphic DNA (RAPD), Figure S3: Genomic pairwise similarities calculated using VIRIDIC, Figure S4: Comparative genomic analysis of phages ER25 and P22, Table S1: Oligonucleotide used in this study, Table S2: BPS6 genome annotation, Table S3: ER6 genome annotation, Table S4: ER22 genome annotation, Table S5: ER20 genome annotation.
Author Contributions
Conceptualization, B.M.P.-S.; methodology, E.V.R., C.V.P., N.W. and B.M.P.-S.; formal analysis, E.V.R., C.V.P. and B.M.P.-S.; investigation, E.V.R., C.V.P., N.W., Y.C., A.B.-B. and B.M.P.-S.; resources T.R.B., N.A.F. and J.C.D.H.; data curation, E.V.R. and B.M.P.-S.; writing-original draft preparation, B.M.P.-S.; writing-review and editing, E.V.R., N.W., C.V.P., Y.C., A.B.-B., A.B., J.R., D.M.P., T.R.B., N.A.F., J.C.D.H., and B.M.P.-S.; visualization, B.M.P.-S.; supervision, A.B., J.R., D.M.P., T.R.B., N.A.F., J.C.D.H., and B.M.P.-S.; funding acquisition, J.C.D.H. All authors have read and agreed to the published version of the manuscript.
Funding
This research was supported by the Wellcome Trust Senior Investigator Award, grant number 106914/Z/15/Z, to Jay C. D. Hinton, and a BBSRC NLD Doctoral Training Partnership studentship number BB/M011186/1, to David M. Picton.
Institutional Review Board Statement
Not applicable.
Informed Consent Statement
Not applicable.
Data Availability Statement
The data presented in this study are publicly available in GenBank, with accession numbers listed in Table 2.
Acknowledgments
We are grateful to present and former members of the Hinton laboratory for helpful discussions. We are indebted to Tom Heyes for technical support and to Josep Casadesús for kindly sharing their bacteriophages.
Conflicts of Interest
Ella V. Rodwell is affiliated to the National Institute for Health Research Health Protection Research Unit (NIHR HPRU) in Gastrointestinal Infections at University of Liverpool in partnership with Public Health England (PHE), in collaboration with University of Warwick. Ella V. Rodwell is based at Public Health England. The views expressed are those of the author(s) and not necessarily those of the NHS, the NIHR, the Department of Health and Social Care or Public Health England. The other authors declare no conflict of interest.
1. Stanaway, J.D.; Parisi, A.; Sarkar, K.; Blacker, B.F.; Reiner, R.C.; Hay, S.I.; Nixon, M.R.; Dolecek, C.; James, S.L.; Mokdad, A.H.; et al. The global burden of non-typhoidal Salmonella invasive disease: A systematic analysis for the Global Burden of Disease Study 2017. Lancet Infect. Dis. 2019, 19, 1312-1324.
2. Okoro, C.K.; Kingsley, R.A.; Quail, M.A.; Kankwatira, A.M.; Feasey, N.A.; Parkhill, J.; Dougan, G.; Gordon, M.A. High-resolution single nucleotide polymorphism analysis distinguishes recrudescence and reinfection in recurrent invasive nontyphoidal salmonella typhimurium disease. Clin. Infect. Dis. 2012, 54, 955-963.
3. Feasey, N.A.; Hadfield, J.; Keddy, K.H.; Dallman, T.J.; Jacobs, J.; Deng, X.; Wigley, P.; Barquist Barquist, L.; Langridge, G.C.; Feltwell, T.; et al. Distinct Salmonella Enteritidis lineages associated with enterocolitis in high-income settings and invasive disease in low-income settings. Nat. Genet. 2016, 48, 1211-1217.
4. Kingsley, R.A.; Msefula, C.L.; Thomson, N.R.; Kariuki, S.; Holt, K.E.; Gordon, M.A.; Harris, D.; Clarke, L.; Whitehead, S.; Sangal, V.; et al. Epidemic multiple drug resistant Salmonella Typhimurium causing invasive disease in sub-Saharan Africa have a distinct genotype. Genome Res. 2009, 19, 2279-2287.
5. Okoro, C.K.; Kingsley, R.A.; Connor, T.R.; Harris, S.R.; Parry, C.M.; Al-Mashhadani, M.N.; Kariuki, S.; Msefula, C.L.; Gordon, M.A.; de Pinna, E.; et al. Intracontinental spread of human invasive Salmonella Typhimurium pathovariants in sub-Saharan Africa. Nat. Genet. 2012, 44, 1215-1221.
6. Brüssow, H.; Canchaya, C.; Hardt, W.D. Phages and the evolution of bacterial pathogens: From genomic rearrangements to lysogenic conversion. Microbiol. Mol. Biol. Rev. 2004, 68, 560-602.
7. Fortier, L.-C.; Sekulovic, O. Importance of prophages to evolution and virulence of bacterial pathogens. Virulence 2013, 4, 354-365.
8. Owen, S.V.; Wenner, N.; Canals, R.; Makumi, A.; Hammarlöf, D.L.; Gordon, M.A.; Aertsen, A.; Feasey, N.A.; Hinton, J.C.D. Characterization of the prophage repertoire of African Salmonella Typhimurium ST313 reveals high levels of spontaneous induction of novel phage BTP1. Front. Microbiol. 2017, 8, 235.
9. Ashton, P.M.; Owen, S.V.; Kaindama, L.; Rowe, W.P.M.; Lane, C.R.; Larkin, L.; Nair, S.; Jenkins, C.; de Pinna, E.M.; Feasey, N.A.; et al. Public health surveillance in the UK revolutionises our understanding of the invasive Salmonella Typhimurium epidemic in Africa. Genome Med. 2017, 9, 1-13.
10. Owen, S.V.; Wenner, N.; Dulberger, C.L.; Rodwell, E.V.; Bowers-Barnard, A.; Quinones-Olvera, N.; Rigden, D.J.; Rubin, E.J.; Garner, E.C.; Baym, M.; et al. Prophage-encoded phage defence proteins with cognate self-immunity. bioRxiv 2020.
11. Kintz, E.; Davies, M.R.; Hammarlöf, D.L.; Canals, R.; Hinton, J.C.D.; van der Woude, M.W. A BTP1 prophage gene present in invasive non-typhoidal Salmonella determines composition and length of the O-antigen of the lipopolysaccharide. Mol. Microbiol. 2015, 96, 263-275.
12. Bernheim, A.; Sorek, R. The pan-immune system of bacteria: Antiviral defence as a community resource. Nat. Rev. Microbiol. 2020, 18, 113-119.
13. Rousset, F.; Dowding, J.; Bernheim, A.; Rocha, E.P.C.; Bikard, D. Prophage-encoded hotspots of bacterial immune systems. bioRxiv 2021.
14. Sambrook, J.; Russell, D.W. Molecular Cloning: A Laboratory Manual, 3rd ed.; Cold Spring Harbor Laboratory Press: Cold Spring Harbor, NY, USA, 2001.
15. Pedulla, M.L.; Ford, M.E.; Karthikeyan, T.; Houtz, J.M.; Hendrix, R.W.; Hatfull, G.F.; Poteete, A.R.; Gilcrease, E.B.; Winn-Stapley, D.A.; Casjens, S.R. Corrected sequence of the bacteriophage P22 genome. J. Bacteriol. 2003.
16. Schmieger, H. Phage P22-mutants with increased or decreased transduction abilities. Mol. Gen. Genet. MGG 1972, 119, 75-88.
17. Casjens, S.R.; Jacobs-Sera, D.; Hatfull, G.F.; Hendrix, R.W. Genome sequence of Salmonella enterica phage Det7. Genome Announc. 2015.
18. Casjens, S.R.; Leavitt, J.C.; Hatfull, G.F.; Hendrix, R.W. Genome sequence of Salmonella phage 9NA. Genome Announc. 2014.
19. Thomson, N.R.; Clayton, D.J.; Windhorst, D.; Vernikos, G.; Davidson, S.; Churcher, C.; Quail, M.A.; Stevens, M.; Jones, M.A.; Watson, M.; et al. Comparative genome analysis of Salmonella Enteritidis PT4 and Salmonella Gallinarum 287/91 provides insights into evolutionary and host adaptation pathways. Genome Res. 2008, 18, 1624-1637.
20. Darboe, S.; Bradbury, R.; Phelan, J.; Kanteh, A.; Muhammad, A.-K.; Worwui, A.; Yang, S.; Nwakanma, D.; Perez-Sepulveda, B.; Kariuki, S.; et al. Genomic diversity and antimicrobial resistance among non-typhoidal Salmonella associated with human disease in The Gambia. bioRxiv 2020.
21. Perez-Sepulveda, B.M.; Predeus, A.V.; Fong, W.Y.; Parry, C.M.; Cheesbroughb, J.; Wigley, P.; Feasey, N.A.; Hinton, J.C.D. Complete genome sequences of African Salmonella Enteritidis clinical isolates associated with bloodstream infection. Microbiol. Resour. Announc. 2021, in press.
22. Fields, P.I.; Swanson, R.V.; Haidaris, C.G.; Heffron, F. Mutants of Salmonella typhimurium that cannot survive within the macrophage are avirulent. Proc. Natl. Acad. Sci. USA 1986, 83, 5189-5193.
23. Zinder, N.D.; Lederberg, J. Genetic exchange in Salmonella. J. Bacteriol. 1952, 64, 679-699.
24. Rankin, J.D.; Taylor, R.J. The estimation of doses of Salmonella typhimurium suitable for the experimental production of disease in calves. Vet. Rec. 1966, 78, 706-707.
25. Simon, R.; Priefer, U.; Pühler, A. A broad host range mobilization system for in vivo genetic engineering: Transposon mutagenesis in gram negative bacteria. Bio/Technology 1983.
26. Martínez-García, E.; de Lorenzo, V. Engineering multiple genomic deletions in Gram-negative bacteria: Analysis of the multi-resistant antibiotic profile of Pseudomonas putida KT2440. Environ. Microbiol. 2011.
27. Datsenko, K.A.; Wanner, B.L. One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products. Proc. Natl. Acad. Sci. USA 2000.
28. Koskiniemi, S.; Pränting, M.; Gullberg, E.; Näsvall, J.; Andersson, D.I. Activation of cryptic aminoglycoside resistance in Salmonella enterica. Mol. Microbiol. 2011.
29. Doublet, B.; Douard, G.; Targant, H.; Meunier, D.; Madec, J.Y.; Cloeckaert, A. Antibiotic marker modifications of λ Red and FLP helper plasmids, pKD46 and pCP20, for inactivation of chromosomal genes using PCR products in multidrug-resistant strains. J. Microbiol. Methods 2008.
30. Heckman, K.L.; Pease, L.R. Gene splicing and mutagenesis by PCR-driven overlap extension. Nat. Protoc. 2007.
31. Pelludat, C.; Mirold, S.; Hardt, W.D. The SopEφ phage integrates into the ssrA gene of Salmonella enterica serovar Typhimurium A36 and is closely related to the Fels-2 prophage. J. Bacteriol. 2003.
32. Goldfarb, T.; Sberro, H.; Weinstock, E.; Cohen, O.; Doron, S.; Charpak-Amikam, Y.; Afik, S.; Ofir, G.; Sorek, R. BREX is a novel phage resistance system widespread in microbial genomes. EMBO J. 2015, 34, 169-183.
33. Hammarlöf, D.L.; Kröger, C.; Owen, S.V.; Canals, R.; Lacharme-Lora, L.; Wenner, N.; Schager, A.E.; Wells, T.J.; Henderson, I.R.; Wigley, P.; et al. Role of a single noncoding nucleotide in the evolution of an epidemic African clade of Salmonella. Proc. Natl. Acad. Sci. USA 2018, 115, E2614-E2623.
34. Canals, R.; Hammarlöf, D.L.; Kröger, C.; Owen, S.V.; Fong, W.Y.; Lacharme-Lora, L.; Zhu, X.; Wenner, N.; Carden, S.E.; Honeycutt, J.; et al. Adding function to the genome of African Salmonella Typhimurium ST313 strain D23580. PLoS Biol. 2019, 17, e3000059.
35. Kumwenda, S.; Kambala, C. Determination of Biological, Physical and Chemical Pollutants in Mudi River, Blantyre, Malawi RBF4MNH Evaluation View project. J. Basic Appl. Sci. Res. 2012, 2, 6833-6839.
36. Kropinski, A.M.; Mazzocco, A.; Waddell, T.E.; Lingohr, E.; Johnson, R. Enumeration of Bacteriophages by Double Agar Overlay Plaque Assay. In Bacteriophages: Methods and Protocols, Volume 1: Isolation, Characterization, and Interactions; Kropinski, A.M., Clokie, M.R.J., Eds.; Humana Press: Totowa, NJ, USA, 2009; pp. 69-76. ISBN 978-1-58829-682-5.
37. Mazzocco, A.; Waddell, T.E.; Lingohr, E.; Johnson, R.P. Enumeration of bacteriophages using the small drop plaque assay system. In Bacteriophages: Methods in Molecular Biology, Volume 1: Isolation, Characterization, and Interactions; Kropinski, A.M., Clokie, M.R.J., Eds.; Humana Press: Totowa, NJ, USA, 2009; pp. 81-85. ISBN 978-1-58829-682-5.
38. Gutiérrez, D.; Martín-Platero, A.M.; Rodríguez, A.; Martínez-Bueno, M.; García, P.; Martínez, B. Typing of bacteriophages by randomly amplified polymorphic DNA (RAPD)-PCR to assess genetic diversity. FEMS Microbiol. Lett. 2011, 322, 90-97.
39. Wick, R.R.; Judd, L.M.; Gorrie, C.L.; Holt, K.E. Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput. Biol. 2017, 13, e1005595.
40. Gurevich, A.; Saveliev, V.; Vyahhi, N.; Tesler, G. QUAST: Quality assessment tool for genome assemblies. Bioinformatics 2013, 29, 1072-1075.
41. Carver, T.; Harris, S.R.; Berriman, M.; Parkhill, J.; McQuillan, J.A. Artemis: An integrated platform for visualization and analysis of high-throughput sequence-based experimental data. Bioinformatics 2012, 28, 464-469.
42. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990.
43. Seemann, T. Prokka: Rapid prokaryotic genome annotation. Bioinformatics 2014, 30, 2068-2069.
44. Tynecki, P.; Guziński, A.; Kazimierczak, J.; Jadczuk, M.; Dastych, J.; Onisko, A. PhageAI-Bacteriophage life cycle recognition with machine learning and natural language processing. bioRxiv 2020.
45. Adriaenssens, E.; Brister, J.R. How to Name and Classify Your Phage: An Informal Guide. Viruses 2017, 9, 70.
46. Sievers, F.; Wilm, A.; Dineen, D.; Gibson, T.J.; Karplus, K.; Li, W.; Lopez, 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.
47. Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 2014, 30, 1312-1313.
48. Letunic, I.; Bork, P. Interactive Tree of Life (iTOL) v4: Recent updates and new developments. Nucleic Acids Res. 2019, 47, W256-W259.
49. Moraru, C.; Varsani, A.; Kropinski, A.M. VIRIDIC-A Novel Tool to Calculate the Intergenomic Similarities of Prokaryote-Infecting Viruses. Viruses 2020, 12, 1268.
50. Page, A.J.; Cummins, C.A.; Hunt, M.; Wong, V.K.; Reuter, S.; Holden, M.T.G.; Fookes, M.; Falush, D.; Keane, J.A.; Parkhill, J. Roary: Rapid large-scale prokaryote pan genome analysis. Bioinformatics 2015, 31, 3691-3693.
51. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0. Mol. Biol. Evol. 2013, 30, 2725-2729.
52. Alikhan, N.F.; Petty, N.K.; Ben Zakour, N.L.; Beatson, S.A. BLAST Ring Image Generator (BRIG): Simple prokaryote genome comparisons. BMC Genom. 2011, 12, 402.
53. Reddy, E.A.; Shaw, A.V.; Crump, J.A. Community-acquired bloodstream infections in Africa: A systematic review and meta-analysis. Lancet Infect. Dis. 2010, 10, 417-432.
54. Labrie, S.J.; Samson, J.E.; Moineau, S. Bacteriophage resistance mechanisms. Nat. Rev. Microbiol. 2010, 8, 317-327.
55. Broeker, N.K.; Barbirz, S. MicroCommentary Not a barrier but a key: How bacteriophages exploit host's O-antigen as an essential receptor to initiate infection. Mol. Microbiol. 2017, 105, 353-357.
56. Gencay, Y.E.; Gambino, M.; Prüssing, T.F.; Brøndsted, L. The genera of bacteriophages and their receptors are the major determinants of host range. Environ. Microbiol. 2019, 21, 2095-2111.
57. Fong, K.; Tremblay, D.M.; Delaquis, P.; Goodridge, L.; Levesque, R.C.; Moineau, S.; Suttle, C.A.; Wang, S. Diversity and Host Specificity Revealed by Biological Characterization and Whole Genome Sequencing of Bacteriophages Infecting Salmonella enterica. Viruses 2019, 11, 854.
58. Choi, Y.; Shin, H.; Lee, J.H.; Ryua, S. Identification and characterization of a novel flagellum-dependent Salmonella-infecting bacteriophage, iEPS5. Appl. Environ. Microbiol. 2013, 79, 4829-4837.
59. Meynell, E.W. A phage, øχ, which attacks motile bacteria. J. Gen. Microbiol. 1961, 25, 253-290.
60. Karpe, Y.A.; Kanade, G.D.; Pingale, K.D.; Arankalle, V.A.; Banerjee, K. Genomic characterization of Salmonella bacteriophages isolated from India. Virus Genes 2016, 52, 117-126.
61. Byl, C.V.; Kropinski, A.M. Sequence of the genome of Salmonella bacteriophage P22. J. Bacteriol. 2000, 182, 6472-6481.
62. Grygorcewicz, B.; Grudziński, M.; Wasak, A.; Augustyniak, A.; Pietruszka, A.; Nawrotek, P. Bacteriophage-mediated reduction of Salmonella Enteritidis in swine slurry. Appl. Soil Ecol. 2017, 119, 179-182.
63. Woolston, J.; Parks, A.R.; Abuladze, T.; Anderson, B.; Li, M.; Carter, C.; Hanna, L.F.; Heyse, S.; Charbonneau, D.; Sulakvelidze, A. Bacteriophages lytic for Salmonella rapidly reduce Salmonella contamination on glass and stainless steel surfaces. Bacteriophage 2013, 3, e25697.
64. Nagel, T.E.; Chan, B.K.; De Vos, D.; El-Shibiny, A.; Kang'ethe, E.K.; Makumi, A.; Pirnay, J.P. The developing world urgently needs phages to combat pathogenic bacteria. Front. Microbiol. 2016, 7, 882.
© 2021. This work is licensed under http://creativecommons.org/licenses/by/3.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.