Morphology of Phage vB_SenS_TUMS_E4
Phage E4 produces large, transparent, round plaques, 1.5 to 2 mm in size and slightly turbid halos on isolated Salmonella enteritidis host cultures (Fig.1a). It is suggested that halo formation results from endolysin secretion after the lysis of the host cells (29). According to Phage morphological analysis by transmitting electron microscopy, the phage exhibited a hexagonal capsid with a diameter of 82±2 nm and a non-contractile tail with a length of 150±5 nm (Fig.1b). Morphological assessment of phages indicated that they belonged to the Siphoviridae family.
Host range properties of the phage
Analysis of the host domain of phage E4 showed that this phage lyses all strains of Salmonella, including some drug-resistant strains (Table 1). At the same time, E4 did not show lytic ability on strains of other genera. Host specificity as a prerequisite for phage therapy is applied in the food industry, particularly for fermented foods or the gut microbiota in farm animals. These results indicate the high specificity of the vB_SenS_TUMS_E4.
Table 1
Host range of phage vB_SenS_TUMS_E4.
No
|
Species
|
strainID
|
Lytic Activity
|
1
|
Salmonella enteritidis ATCC 13076
|
ATCC 13076
|
+
|
2
|
Salmonella typhi ATCC 14023
|
ATCC 14023
|
+
|
3
|
Salmonella Paratyphi
|
Lab Collection
|
+
|
4
|
Salmonella entrica
|
Lab Collection
|
+
|
5
|
Salmonella typhi
|
Lab Collection
|
+
|
6
|
Salmonella Paratyphi C
|
Lab Collection
|
+
|
7
|
Salmonella clinical
|
Lab Collection
|
+
|
8
|
Salmonella clinical
|
Lab Collection
|
+
|
9
|
Salmonella clinical
|
Lab Collection
|
+
|
10
|
Salmonella clinical
|
Lab Collection
|
+
|
11
|
Escherichia coli PTCC 1330
|
PTCC 1330
|
-
|
12
|
Escherichia coli ATCC 25922
|
ATCC 25922
|
-
|
13
|
Staphylococcus aureus ATCC 6588
|
ATCC 6588
|
-
|
14
|
Staphylococcus aureus ATCC 23591
|
ATCC 23591
|
-
|
15
|
Staphylococcus aureus ATCC 33591
|
ATCC 33591
|
-
|
16
|
Staphylococcus aureus ATCC 25923
|
ATCC 25923
|
-
|
17
|
Pseudomonas syringae ATCC 10205
|
ATCC 10205
|
-
|
18
|
Pseudomonas aeruginosa ATCC 9027
|
ATCC 9027
|
-
|
19
|
Pseudomonas aeruginosa ATCC 27853
|
ATCC 27853
|
-
|
20
|
Methicillin-resistant Staphylococcus aureus (MRSA)
|
Lab Collection
|
-
|
21
|
Staphylococcus saprophyticus
|
Lab Collection
|
-
|
22
|
Proteus mirabilis
|
Lab Collection
|
-
|
23
|
Enterococcus faecalis
|
Lab Collection
|
-
|
24
|
Serratia marcescens
|
Lab Collection
|
-
|
25
|
Pseudomonas putida
|
Lab Collection
|
-
|
26
|
Enterobacter aerogenes
|
Lab Collection
|
-
|
27
|
Citrobacter freundi
|
Lab Collection
|
-
|
28
|
Shigella flexneri
|
Lab Collection
|
-
|
29
|
Citrobacter freundii
|
Lab Collection
|
-
|
30
|
Staphylococcus epidermidis
|
Lab Collection
|
-
|
31
|
Acinetobacter baumannii
|
Lab Collection
|
-
|
32
|
Klebsiella pneumoniae
|
Lab Collection
|
-
|
33
|
Lactobacillus salivarius
|
Lab Collection
|
-
|
34
|
Lactobacillus jensenii
|
Lab Collection
|
-
|
35
|
Lactobacillus plantarum
|
Lab Collection
|
-
|
36
|
Lactobacillus reuteri
|
Lab Collection
|
-
|
One-step growth of phage vB-SenS-TUMS-E4
A one-step growth test determined the latent period as well as phage size of vB_SenS_TUMS_E4. The latent period was about 15 min, and the burst size was approximately 287 plaque-forming units (PFU) per cell (Fig. 2a). The phage burst size is highly associated with its proliferation, and a right burst size is a favorable feature for an effective lytic bacteriophage. Thus, phages characterized by large burst size and short latency are suitable candidates for phage therapy (30). Hence, the found burst size and the short latency of VB_SenS_TUMS_E4 can be regarded as favorable features for its use in phage therapy.
Time-kill analysis of Salmonella enteritidis treated with phage
The phage lytic activity was assessed against S. enteritidis culture at its exponential developmental phase (OD600=0.4) with various MOIs of E4. After 2 h, using MOIs of 0.01, 0.1, 1, 10, and 100 led to a nearly 3.5 Log dcrease in the bacterium titer to approximately 102 CFU/ml (P<0.0001). Phage titers were similar at various MOIs, suggesting that the E4 phage titer was less affected by MOI. Nevertheless, lower MOIs (0.1 and 0.01) caused a moderate reduction in the bacterium titer in early hours (Fig. 2b). Thus, the application of higher phage levels led to a faster decrease in the bacterium count, possibly because of the elevated attachment rate at a higher phage titer. Also, incubation for 24 h significantly increased in the bacterium titer in both phage-treated and control samples. The recorded bacterium titers were about 2 Log less compared to the control (P<0.0001). Therefore, the observed decrease in the host bacterium cell count within the first three hours of being exposed to E4 phage and preserving such trend up to 8 h following incubation can be a marked feature for the possible use of this phage in phage therapy.
Temperature and pH Stability
Using phages as biocontrol agents to inhibit bacteria requires their stability under different stress conditions, like pH and temperature. pH stability testing of the E4 indicated that it is highly stable at pH 5 to 9 (P<0.0001). Nevertheless, under acidic conditions (pH=3), after 90 min, 4 logs of phage titer were observed (Fig. 2c).
Fig. 2d shows phage E4 is relatively stable at high temperatures without significant titer loss up to 50 ⁰C following incubation for 90 min (P<0.0001). After incubation of the phage at 60 ⁰C, PFU/mL decreased 2 logs in titer. Exposing phage E4 at temperatures more than 70 ⁰C caused rapid inactivation and titer reductions of 6 log PFU/ml.
Restriction profile
The DNA of phage E4 could be digested by both HindIII and EcoRI (Fig. 3), which indicated that the nucleic acid contained in phage E4 is double-stranded DNA.
vB_SenS_TUMS_E4 genome characterization
We obtained a complete assembled phage genome of vB_SenS_TUMS_E4 (43018 bp) and 49.7% C+G content from genome sequencing and annotation. The PhageTerm and PCR analysis results showed that the genome is circularly permuted (data not shown)(31). Based on RAST annotation analysis, a total of 60 protein encoding genes were predicted for E4, of which 38 genes are on the positive strand, and 22 genes are on the negative strand (Fig. 4). About 89% of genes start with ATG, and 11% with GTG (E4-2, 9, 49, and 54) and TTG (E4-5, 28, and 33). Three types of end codons were used to predict ORFs, including TAA (53%, 32 ORFs), TGA (33%, 20 ORFs), and TAG (13%, 8 ORFs). Functional analysis using Blastp revealed 43 genes encoding proteins with known functions and 17 genes with unknown functions (Additional File 1: Table S1). The vB_SenS_TUMS_E4 genome can be divided into four modules based on the predicted functions of its component genes: (i) phage assembly and structure, (ii) DNA packaging, (iii) replication and regulation proteins (iv) host lysis (Fig. 5). In addition, 21 transcriptional promoter sequences were recognized using Phage Promoter software (Additional File 1: Table S2). Thirty-seven Rho factor independent terminators were determined in the E4 genome by Arnold (Additional File 1: Table S3), which were evaluated based on the existence of a U-rich tail, their location, and strongly anticipated stem-loop secondary structure (ΔG≤−10 kcal/mol) determined using MFold (32). The completeness of the phage E4 was confirmed after mapping (99.94%) the reads with BWA-MEM (24). Our analysis showed that the E4 phage lacked known genes that encode antibiotic resistance, toxins, or virulence factors. Therefore, this phage is possibly a virulent phage against Salmonella enteritidis.
Phage structure and assembly
Based on Blastp bioinformatics analysis, the genome of Phage 4 contains 20 structural proteins. Fifteen genes related to tail and tail assembly (E4-3, 8, 14-21, 27-29, 33, and 34), five genes related to capsid and capsid assembly, and head morphogenesis were identified (E4-7, 11-13, and 46) (Additional File 1: Table S1). The structural module occupies the positive strand of genome, 20 genes including structural proteins and those involved in virion assembly.
Nine protein bands presenting the virion structural components were observed in the SDS polyacrylamide Coomassie-stained gel (Fig. 6). A dominant polypeptide bond of about 37 kDa indicates the presence of the original capsid protein resulting in a high number of copies of capsid proteins. The molecular mass identified corresponds to the protein anticipated molecular weight (Fig. 6). Major capsid protein is a highly conserved protein. A polypeptide bond of about 72 kDa indicates the presence of the tail spike protein (Fig. 6).
DNA packaging
Above the head module is the E4 DNA packaging module. The packaging machine usually consists of a portal ring and a terminase complex in tailed phages (29, 30). The most characteristic terminases are hetero-oligomers consisting of a small subunit (E4-1) involved in DNA detection and a large terminal subunit (E4-2) containing ATPase and endonuclease activities. E4-2 encodes terminase large subunit protein, which shows more than 99% identity (with query coverage 100%) with terminase large subunit protein encoded by VSe103, SE2 (Additional File 1: Table S1).
Replication and regulation proteins
E4 genome possessed four core genes associated with phage replication, like DNA polymerase I (E4-41), homing endonuclease (E4-24, 38, and 45), helicase (E4-48), and exonuclease (E4-36). Based on the analysis using Blastp, three regulatory core genes were identified. The gene coding for helix-turn-helix for regulating gene expression (E4-26,48, and 51), the gene coding for putative NinH-like protein transcriptional regulation, and phage site-specific recombination (E4-59), and the gene coding for putative DNA-binding protein (E4-6), were found.
Lysis
Lysis genes are usually found near each other and form a lysis cassette. In some siphoviridae members, including vB_SenS_SE1, wks13, FSL-SP031, K1H, and E4, the endolysin genes are upstream of the packaging genes (Fig. 5) (23). E4-56 encodes lysin protein, which shows more than 93% identity (with query coverage of 100%) with endolysin encoded by VSe103, SE2, and Ent1.
Phylogeny analysis
The E4 phage genome was compared with already sequenced phages in GenBank to study the evolutionary relationship of phage. Then, the E4 genome and the phages belonging to the Jerseyvirus were analyzed with JSpeciesWS, and the results showed that the ANI values for phages were above 80% (Additional File 1: Table S4)(33). Thus, phage E4 can be a member of the Jerseyvirus genus. Phylogenetic analysis based on the complete genome showed that phage E4 is closely related to bacteriophage 5sent1 (Fig. 7) belonging to the genus of Jerseyvirus. These findings are in line with those of Jersey phages who indicated low sequence similarity to those in the Siphoviridae family(7, 34).
Statistical Analysis
Data are presented as mean ± standard deviation. We applied GraphPad Prism 8.0.1 for statistical analysis. Two-way ANOVA with multiple comparisons was employed for indicating differences between groups at < 0.0001.