Next Article in Journal
EWI-2 Inhibits Cell–Cell Fusion at the HIV-1 Virological Presynapse
Next Article in Special Issue
Promises and Pitfalls of In Vivo Evolution to Improve Phage Therapy
Previous Article in Journal
Characterization of Klebsiella pneumoniae ST11 Isolates and Their Interactions with Lytic Phages
Previous Article in Special Issue
Look Who’s Talking: T-Even Phage Lysis Inhibition, the Granddaddy of Virus-Virus Intercellular Communication Research
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Benthic Phage Infecting Shewanella with Strong Replication Ability

1
Key Laboratory of Biofuels, Shandong Provincial Key Laboratory of Energy Genetics, Qingdao Institute of Bioenergy and Bioprocess Technology, Chinese Academy of Sciences, Qingdao 266101, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
CAS Key Lab of Separation Sciences for Analytical Chemistry, National Chromatographic Research and Analysis Center, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
*
Author to whom correspondence should be addressed.
Submission received: 4 November 2019 / Accepted: 17 November 2019 / Published: 19 November 2019
(This article belongs to the Special Issue Phage Ecology)

Abstract

:
The coastal sediments were considered to contain diverse phages playing important roles in driving biogeochemical cycles based on genetic analysis. However, till now, benthic phages in coastal sediments were very rarely isolated, which largely limits our understanding of their biological characteristics. Here, we describe a novel lytic phage (named Shewanella phage S0112) isolated from the coastal sediments of the Yellow Sea infecting a sediment bacterium of the genus Shewanella. The phage has a very high replication capability, with the burst size of ca. 1170 phage particles per infected cell, which is 5–10 times higher than that of most phages isolated before. Meanwhile, the latent period of this phage is relatively longer, which might ensure adequate time for phage replication. The phage has a double-stranded DNA genome comprising 62,286 bp with 102 ORFs, ca. 60% of which are functionally unknown. The expression products of 16 ORF genes, mainly structural proteins, were identified by LC-MS/MS analysis. Besides the general DNA metabolism and structure assembly genes in the phage genome, there is a cluster of auxiliary metabolic genes that may be involved in 7-cyano-7-deazaguanine (preQ0) biosynthesis. Meanwhile, a pyrophosphohydrolase (MazG) gene being considered as a regulator of programmed cell death or involving in host stringer responses is inserted in this gene cluster. Comparative genomic and phylogenetic analysis both revealed a great novelty of phage S0112. This study represents the first report of a benthic phage infecting Shewanella, which also sheds light on the phage–host interactions in coastal sediments.

1. Introduction

In their natural environment, phages are important agents influencing the physiology and survival of prokaryotes and have significant impacts on the microbially driven biogeochemical cycles [1,2,3]. In recent decades, extensive studies on marine phages in seawater have been carried out via field investigation, laboratory isolation, gene-based analysis, and metaproteomic methods. Now we have a certain understanding of the biological characteristics, community composition, and ecological functions of the planktonic phages. The total amount of bacteriophage in the ocean is estimated to be ca. 1030, about 15-fold times than prokaryotes [2]. Phage infections significantly control the bacterial mortality and, through viral shunt, affect nutrient cycling and food web dynamics [2,4]. Recently, the most abundant plankton in the oceans, SAR11 bacteria, formerly considered immune to phage predation, was found also susceptible to dominant phages (i.e., pelagiphages) in seawater [5]. Meanwhile, phage infection can promote the coevolution of bacteria and phages [4]. Many phages have been found to contain host genes via lateral gene transfer and participate in the host metabolism, such as the genes (psbA and psbD) encoding for photosynthetically important proteins, which were recently found in a number of cultured cyanophage genomes [6]. By contrast, so far there is still a great lack of knowledge about benthic phages, particularly the benthic phages in coastal sediments. For example, although previous gene-based analysis has revealed that the coastal sediments contain extremely diverse phages and may play important roles in driving carbon cycles [7,8], benthic phages have been rarely isolated from coastal sediments in the laboratory, which largely limits our understanding of their biological characteristics. At present, isolation and characterization of benthic phages infecting ecologically important bacteria is especially vital for us to explore the mystery of phage activities in coastal sediments.
Shewanella is an ecologically important bacterial genus, whose members are widely distributed in freshwater, seawater, sediment, and deep ocean [9]. Members of the genus Shewanella have multiple metabolic capabilities and have potentials in bioremediation and microbial energy generation [10]. Meanwhile, they were considered to be model organisms to study the microbially driven carbon-cycling process [11]. Till now, only nine Shewanella phage genomes have been sequenced, including four siphoviruses (i.e., 1/44, 3/49, Spp001, and SppYZU05) and five mycoviruses (i.e., 1/4, 1/40, 1/41, SFcil, and SppYZU01) [12,13,14,15]. The nine Shewanella phages were isolated from various environments, but none was from marine sediments. Phage SppYZU01 infecting S. baltica and phages SppYZU05, Spp001 infecting S. putrefaciens were isolated from waste effluents of a seafood market [12,13]. Phage SFCi1 infecting S. fidelis was isolated from the seawater of Mission Bay [14]. The remaining five Shewanella phages were isolated from Baltic Sea ice [15]. As we know, members of Shewanella are widespread in marine sediments [9] and play important roles in extracellular electron transfer [16], reduction of humic acids and iron [10], and degradation of persistent organic compounds, such as petroleum, chlorinated solvents, tiazine, and tetrachloroethene [17]. However, phages infecting Shewanella spp. from marine sediments had never been reported. Even sediment phages are rarely reported. So far, to the best of our knowledge, only two phages isolated from coastal sediments have been comprehensively identified (including the genome annotation), whose hosts are Vibrio parahaemolyticus and Acinetobacter baumannii, respectively [18,19].
In this study, Shewanella indica KJW27, which was isolated from coastal sediments of the Arabian Sea [20], was used as the host bacterium, to isolate its phages from the coastal sediments of Yellow Sea, China. Via biological identification and genome analysis, we found that certain unique characteristics of this benthic phage might be related to phage life strategy and environmental adaptability in coastal sediments. The isolation and characterization of this phage open a window for future explorations of the interactions between phage and the ecologically important Shewanella spp. in marine sediments.

2. Materials and Methods

2.1. Isolation and Purification of Benthic Phages

The bacterial strain used in this study was S. indica KJW27, which was isolated in 2011 from the coastal sediments of Arabian Sea, India [20]. The strain was grown on RO medium (yeast extract 1 g·L−1, peptone 1 g·L−1, sodium acetate 1 g·L−1, artificial seawater 1 L, pH 7.8–8.0), at 30 °C, in a shaking incubator. Coastal sediment samples (5 g) were collected from Aoshan Bay of the Yellow Sea, China (120°79′ N; 36°44′ E), in September 2018, and immediately transferred into a 100 mL of RO culture medium, for 7 days, in order to promote phage reproduction and increase the probability of phage isolation. Then, the filtrate containing phage particles were collected by filtration through 0.22 µm pore-size filters (Millipore, Bedford, MA, USA). One milliliter of the phage filtrate was added to S. indica KJW27 culture (100 mL) and incubated for 24 h, at 30 °C, in a shaker. The culture was filtered again and serially diluted to determine phage activity via the double-layer agar method [21]. Single phage plaques were picked up from the plate by using a sterile pipette and then purified five times, using the double-layer agar method. The purified phages were cultivated for expansion in liquid RO medium. The mixture was centrifuged at 10,000× g at 4 °C, for 15 min, to collect the phage-containing supernatant. Supernatants were collected after filtration through 0.22 μm pore-size filters, in order to remove the remaining cells and debris. Phage lysates were then treated with 2 ng·L−1 of DNase I and 2 ng·L−1 RNase A, at room temperature, for 1 h. After digestion, NaCl was added to the treated lysates (final concentration: 1 M) and incubated at 4 °C for 1 h. The treated lysates were centrifuged at 10,000× g for 10 min at 4 °C. The supernatant was filtered through a 0.22 µm pore-size filter to remove the debris. The filtrate was concentrated by polyethylene glycol 8000 precipitation (final concentration: 100 g·L−1) overnight, at 4 °C. The mixture was centrifuged at 10,000× g, for 60 min, at 4 °C, and the precipitates were resuspended in 6 mL of SM buffer (100 mM of NaCl, 8 mM of MgSO4, 50 mM of Tris-HCl, at pH 7.5). Afterward, the phages were purified by CsCl gradient ultra-centrifugation (gradient-density: 1.5 g·mL−1, 200,000× g, 8 h, 4 °C; CP-100WX, Hitachi Limited, Tokyo, Japan). The purified phage particles were collected and dialyzed three times in SM buffer, using 30 kDa super-filters (UFC5030, Millipore), and they were then used for morphologic observation and phage genome extraction.

2.2. Phage Morphological Observation Using Transmission Electron Microscopy (TEM)

The CsCl-purified phages were prepared for imaging on a 200-mesh copper grid, using a negative stain with 2% aqueous uranyl acetate, as described previously [22]. Samples were viewed at 80 kV voltage, using an H-7650 transmission electron microscope (Hitachi Limited, Tokyo, Japan). Images were taken using GATAN INC CCD image transmission system (Gatan Inc., Pleasanton, CA, USA).

2.3. Detection of Phage Host Range

The host range of phage S0112 was tested by the spot test method with the ability to form plaques on bacterial lawn culture. Besides S. indica KJW27, the tested bacteria included seven strains from marine sediments, i.e., Alginatibacterium sediminis ALS 81, Woeseia oceani SDUM189001, Sediminicola luteus SDUM 701001, Kordiimonas sediminis N39, Shewanella chilikensis JC5, Shewanella basaltis CJW-54, and Shewanella japonica KCTC22435; ten strains from seawater or saline water, i.e., Marinobacter vinifirmus D7035, Halomonas denitrificans D7027, Vibrio alginolyticus CIP 82.01, Vibrio neocaledonicus NC 470, Shewanella algae JCM 21037, Vibrio azureus NBRC 104587, Vibrio harveyi NBRC 15634, Vibrio parahaemolyticus NBRC 12711, Vibrio campbellii CAIM 519, and Roseobacter denitrificans OCH114. Briefly, 2 µL of phage suspension (~106 phage particles) was spotted onto the surface of double-layer agar plate inoculated with the tested strains at 30 °C, for two weeks, and plaque formation was examined each day. Three replicates were tested for each bacterial strain.

2.4. Lipid Detection in the Viral Capsid

To determine whether there were lipids in the viral capsid, 1 mL of phage suspension was mixed with 0, 20, and 200 µL of pure chloroform, vigorously shaken for 1 min, and then kept at room temperature for 30 min. The samples were centrifuged at a low speed to ensure that the phages were kept in the supernatant. The supernatant containing phages were dropped onto an S. indica KJW27 plate, in order to observe the emergence of plaques [23].

2.5. The Influence of the External Factors on Phage Particles Stability

To investigate the effect of external factors on the stability of the phage particles, the following external factors were tested: pH (2, 4, 10, and 12), temperature (−20, 15, 20, 28, 30, 37, 40, 50, 60, 95, and 100 °C), detergents (0.1% CTAB, 0.09% SDS, and 0.1% Sarkosyl), and organic solvents (63% ethanol, 90% acetone, 50% DMSO, and chloroform). The effect of osmotic shock on phage particles production was also analyzed, as previously described [24].

2.6. Lysis Profile Assay

Host culture was inoculated with phage lysate at a multiplicity of infection (MOI) of 0.01. Their density was monitored by OD600 measurement at the 30 min intervals up to 8 h. The bacterial density (OD600), survival of host bacteria after phage infection (CFU·mL−1) and lysate titer (PFU·mL−1) were analyzed during the experiment. To investigate the number of surviving cells after phage infection, 100 μL samples were harvested at times indicated above, and serial dilutions (10-fold each) were made in RO medium. Then, 50 μL of each dilution was spread onto RO agar plates. The CFU of each sample was calculated by counting the colonies. To estimate the lysate titer (PFU·mL−1), 100 μL samples were removed every 30 min, and their serial 10-fold dilutions were prepared in SM buffer. Then, the PFU of each dilution was calculated by counting the plaques on the bacterial lawn.

2.7. One-Step Growth Curve

One-step growth curve experiments were performed, as previously described [25]. Briefly, phages were added to 1 mL aliquots of the bacterial culture at mid-exponential growth phase at MOI = 0.01. After 20 min, cells were centrifuged at 6000× g for 10 min at 4 °C, to remove the non-adsorbed phages in the supernatant, and the pellets were resuspended in 1 mL of RO medium. This process was repeated twice. Then, 50 µL of the resuspended culture (the bacteria and adsorbed phages) was transferred to 50 mL of RO medium and incubated over 8 h at 30 °C [26]. Two sets of duplicate samples were taken every 15 min, for eight hours, and chloroform (1% final concentration) was added to the second set, in order to release the intracellular phage. The two samples were serially diluted (10-fold each) and immediately plated for phage titration, using the double-layer agar plate method [27]. Three replicates of this experiment were carried out in this study.
To investigate whether the phage replication was independent of the host, 50 mL of the host culture was added, with rifampicin, at various concentration, including 0, 5, 10, and 20 μg·mL−1 (final concentration) [28]. Twenty min after addition of the rifampicin, the phage lysate was added to the host culture at MOI of 1.0 and incubated for 20 h at 30 °C. The phage titer was assayed as described above.

2.8. Phage Genome Extraction

The CsCl-purified phages were incubated with proteinase K (100 mg·mL−1), SDS (100 μg·mL−1), and EDTA (0.5 mol·L−1, pH 8.0) and incubated at 55 °C, in water, for 3 h. The digested sample was extracted by using an equal volume of phenol/chloroform/isoamyl alcohol (25:24:1). The samples were completely mixed before centrifugation at 12,000× g and 4 °C for 5 min. The aqueous phase was purified by adding chloroform/isoamyl alcohol (24:1) before centrifugation at 12,000× g and 4 °C for 10 min. This step was repeated twice. The aqueous phase was mixed with isopropyl alcohol and incubated overnight at −20 °C, before centrifugation at 12,000× g for 15 min. The precipitate was washed twice with cold ethanol (70%) and then resuspended in 100 μL of TE buffer (10 mM of Tris-HCl, 1 mM of EDTA, pH 8.0). The purified phage genomic DNA was stored at −20 °C before analysis.

2.9. Genomic Analysis

Sequencing libraries were constructed by using TruSeq DNA Library Prep Kits (Illumina, San Diego, CA, USA). Whole genome was sequenced by using the Illumina Hiseq X Ten (2 × 150 bp), and the sequencing data reached coverage of 50,000× and were assembled by using SPAdes genome assembler (v3.12.0), with default k-mer lengths. Open reading frames (ORFs) were predicted with the GeneMarkS online server (http://exon.gatech.edu/Genemark/genemarks.cgi) and were then verified manually, using the NCBI ORF Finder online server (http://www.ncbi/nlm.nih.gov/orffinder/). The putative function of translated products was analyzed and annotated, using BLAST searches against the NCBI non-redundant (nr) protein database with E-value ≤ 10−3. A termini analysis and PhageTerm were used to identify the phage’s termini and genome packaging [29,30], while tRNAscan-SE was used to scan tRNA sequences [31]. Progressive Mauve algorithm was used to compare the phage S0112 with the other phages infecting Shewanella spp. isolated from different environments and phage VpKK5 isolated from coastal sediments [32]. Gene maps were created by using the Java Operon combined with the genome annotations. To test any spacers of CRISPR array within the phage genome, the phage genome was searched against viral spacer database of the Integrated Microbial Genome/Virus (IMG/VR) database (https://img.jgi.doe.gov/cgi-bin/vr/main.cgi) [33], as well as in spacers CRISPRs loci of its host S. indica KJW27, the isolates used in this study [34].

2.10. Proteomic Analysis

Proteins of phage S0112 were prepared according to the filter-aided sample preparation procedure [35]. All samples were analyzed on a Dionex UltiMate 3000 nanoLC System (Thermo Fisher Scientific, Waltham, MA, USA) coupled online to an LTQ Orbitrap Velos mass spectrometer (Thermo Fisher Scientific, Waltham, MA, USA). The data were acquired in a data-dependent mode. The obtained spectrometry information was searched against the database created from translated open reading frames that have been found in the genome of phage S0112 by using Mascot 2.3.02 software (Matrix Science Inc., Boston, MA, USA). The important parameter settings for protein identification were as follows: peptide mass tolerance = 20 ppm, MS/MS tolerance = 0.05 Da, enzyme = trypsin, missed cleavage = 2, and fixed modification: Carbamidomethyl (C), and variable modification (M).

2.11. Construction of Phylogenetic Tree

The amino acid sequences of DNA polymerase subunit B of phage S0112 was aligned with those of the closely related bacteriophages, deposited in the NCBI database, using MUSCLE implanted in the MEGA 7.0 [36]. The maximum likelihood phylogenetic tree with 1000 bootstrap replications for DNA polymerase subunit B amino acid sequences was constructed by the model of Poisson. To find the closest relatives of phage S0112, the ViPTree was used to construct the viral proteomic tree [37].

2.12. Genome Recruitment

In order to evaluate the similarity between the phage S0112 genome and the virioplankton communities, the sequences of the phage ORFs and the databases of the Pacific Ocean Virome (POV), Global Ocean Survey (GOS) were compared by using reciprocal best-hit BLAST [5,38,39]. Samples in POV and GOS were collected from aquatic environments during various seasons, at multiple depths, and from different habitats. In addition, to completely compare whole genome rather than only the ORF, we search for the phage S0112 genome in the IMG/VR database, which is the largest available database related to viral genomics, including the genome of cultivated viruses, prophages derived from cultivated microbial isolates, and viral fragments generated by the metagenomes [33].

2.13. Nucleotide Sequence Accession Number

The genome sequence of the phage S0112 is available under accession number MK675901. The raw reads were submitted to the SRA and can be retrieved via SRA PRJNA573990.

3. Results and Discussion

3.1. Phage Basic Characteristics: Infection Mode, Morphology, Host Ranges, and Existence or Absence of Lipids in Phage Capsid

The phage S0112 is a lytic phage infecting S. indica KJW27, which can form small clear round plaques on plates with an average diameter of 0.5–1 mm after two days of incubation with bacterial hosts (Figure 1a). Morphological analysis via TEM showed that phage S0112 had a polyhedral capsid (62 ± 3 nm in diameter) and a long, flexible, non-contractile tail (97 ± 1 nm in length and 10 ± 1 nm in width) (Figure 1b).
Host range analysis revealed that phage S0112 could not infect any other bacterium except S. indica KJW27 among the eighteen tested strains, which include four different Shewanella spp. from marine sediments and seawater that have similar evolutionary status with S. indica KJW27 (Table 1). After being treated with chloroform, the phage can still survive on a bacterial lawn, indicating the absence of lipids in phage capsid. As introduced above, only two phages (i.e., VpKK5 and AB1) isolated from coastal sediments have been identified comprehensively in the past [18,19]. Based on the few reports on coastal sediment phages, the three known benthic phages including phage S0112 in this study were found to affect only very limited species within a specific genus. However, we believe that the benthic phages in coastal sediments must have much more complex and diverse biological characteristics than what we know now and remain to be explored. Isolating more phages from coastal sediments is of great significance to discover their largely unknown characteristics.

3.2. Lysis Profile Assay and Sensitivity of Virions to Physical and Chemical Factors

The entire course of the lysis of host bacteria spends almost five hours (Figure 2). Meanwhile, we found that phage S0112 could infect S. indica KJW27 and form clear plaques on plates within a pH range of 4–12, temperature from −20 to 60 °C, at the presence of organic solvents (63% ethanol, 90% acetone, 50% DMSO, and chloroform), detergents (0.1% Sarkosyl and 0.1% CTAB), and osmotic shock, whereas when the pH value was lower than 2, or the temperature was higher than 95 °C, or there was 0.09% SDS, the phage infection activity disappeared (Table 2).

3.3. High Replication Capability Reflected by the High Burst Size

The one-step growth curve of the phage S0112 revealed that its eclipse and latent period of phage infection occurred during 2 h 45 min and 3 h post-infection, respectively (Figure 3), which is much longer than most phages isolated before (Figure 4). We randomly counted sixty-eight phages isolated from aquatic environments infecting strains such as Roseobacter spp., Vibrio spp., Pseudomonas spp., Aeromonas spp., Acinetobacter spp., Escherichia spp., etc., and found that 91% of the viruses have the latent period of less than 2 h (Table S1). Comparatively, the latent periods of the phages which were isolated from seawater or sewage and have similar hosts with that of phage S0112 (i.e., Shewanella spp. strains) are only 12.7–60 min [12,13,14]. We speculate that the physiological state or doubling time of bacterial hosts may vary in different environments, which possibly have influences on the phage latent period. The hosts listed in Table S1., such as Vibrio spp., Escherichia spp., and Acinetobacter spp., usually have a rapid doubling time of <30 min, whereas the doubling time of the phage S0112′s host Shewanella indica KJW27 is about 70 min. At present, it is not clear about the potential relationship between the host doubling time and the phage latent period, which is worth exploring in the future.
Concurrent with the long phage latency, the burst size of phage S0112 was as high as 1170 phage particles per infected cell, which is 5–10 times higher than that of most phages. It has been suggested that the increase of latent period might result in the formation of more progeny phages, i.e., higher burst size [40]. Based on this, it is not abnormal for phage S0112 to have such a high burst size, because the longer latent period can provide ample time for the replication and packaging of progeny phages. Moreover, the host metabolic activity and its survival environment were considered to be closely related to phage latency and burst size [41]. Fast-growing cells under favorable environments can produce more progeny phages than in tough environments [41]. In this study, phage S0112 and its host were all isolated from nutrient-replete coastal sediments. High content of organic matter in coastal surface sediments can boost bacterial growth. When vigorous bacteria are infected by phages, their stronger metabolic activities might provide more energy for viral replication. This is likely why phage release or phage burst size is often higher in laboratory culture than in natural environment. Based on the very limited reports of laboratory-cultured benthic phages, the other two coastal sediment phages, i.e., VpKK5 and AB1, infecting Vibrio parahaemolyticus and Acinetobacter baumannii respectively, also showed relatively higher burst sizes, which is 180 and 409 phage particles per infected cell [18,19].
In addition, viral production in marine sediments was reported to be higher than in the water column [42]. Although the nutrient-rich coastal sediments are beneficial for bacterial growth, the complexity of sediment components, such as the presence of abundant extracellular enzymes, may be detrimental to phages and promote phage decay and degradation [42]. Therefore, we propose that the strong replication ability of benthic phages might be a survival strategy adapting to the coastal sediment environment to ensure that enough progeny phages can survive the tough conditions. Meanwhile, the fluidity of phages in sediments is relatively weaker than that in water columns, which may lead to a lower contact rate between phages and their hosts. The benthic phages with high burst size can produce more progeny phages, to ensure the adequate occurrence of phage–host contacts and maintain their own survival via infections.
To test if the phage S0112 replication is independent of host, we analyzed the sensitivity of phage infection to rifampicin, which can inhibit host RNA polymerase (RNAP) by binding to the β subunits [28]. The addition of 5 μg·mL−1 of rifampicin to host culture at prior (20 min) infection completely abolished the production of phage progeny, suggesting that phage S0112 cannot complete the infection cycle independently without host RNAP (Table S2).

3.4. Genomic Features of the Phage S0112

The phage S0112 genome consisted of 62,286 bp of linear ds DNA, with the average G+C content of 44.7%, which is lower than that of its host (51.2%) [20]. PhageTerm analysis showed that the DNA termini and packaging mechanisms of phage S0112 is similar with that of T5-like phages. No tRNA genes were found in the phage genome. Furthermore, the genome of phage S0112 was searched against viral spacer database of IMG/VR and spacers within CRISPRs of its host. The result showed no match between phage S0112 and viral spacers sequences within CRISPRs. The phage genome was predicted to encode 102 putative open reading frames (ORFs), of which only 42 ORFs were assigned putative functions based on their amino acid sequence homology to known proteins. About 59% of the ORFs (60 ORFs) are functionally unknown, among which 21 ORFs were annotated as hypothetical proteins and 39 ORFs had no any homology with protein databases (Table S3). Among the 42 ORFs with known functions, eight ORFs shared similarities with that of Vibrio phage VpKK5, which is also a siphovirus isolated from the marine sediments [18], with percent identities of 42–59%; six ORFs share similarities with the Pantoea phage vB_PagS_Vid5 and Pseudomonas phage MP1412; and eight ORFs shared low identity with Alphaproteobacteria phage Jl001, Roseophage RDJL Ф 1 and Burkholderia phage BcepGomr. All the genes with homology to S0112 are from siphoviruses.
So far, nine Shewanella phage genomes were reported [12,13,14,15]. These phages vary significantly in phage size, GC content, protein amount, and isolation sources (Table 3). Multiple genome alignments of the nine Shewanella phages and phage S0112 showed that the genome of phage S0112 has a very small amount of similarity with others [32] (Figure 5). Among phage S0112 genome, only three ORFs of phage S0112 show certain similarities with other Shewanella phages, including S0112_005 (DNA/RNA helicase), S0112_018 (primase/helicase), and S0112_053 (hypothetical protein).
The encoded proteins of phage S0112 can be divided into four functional modules: DNA packaging module, structure/morphogenesis module, DNA replication/recombination/metabolism module, and queuosine biosynthesis module (Figure 6). The DNA packaging module mainly includes the gene product S0112_063 (terminase small subunit) and S0112_064 (terminase large subunit, sharing 49.3% amino acid identity with the corresponding proteins of Fodinicurvate sediminis). It is reported that terminase generally consists of one larger submit and one small submit. The role of the terminase small submit is to specifically recognize the packaging initiation site (cos or pac) and helps assemble a holoenzyme complex with a larger submit which can cleave genomic DNA at special site, to ensure the whole genome into preassembled empty capsids, using its endonuclease and ATPase activities [43].
The structure/morphogenesis module includes head–tail connector protein (S0112_066), major capsid protein (S0112_071 and _076), tail completion protein (S0112_074), tail terminator protein (S0112_075), tail chaperonin (S0112_078), tail length tape measure protein (S0112_079), virion structural protein (S0112_081), tail assembly structural protein (S0112_082, _084, _085, and _086), minor tail protein (S0112_083 and _087), peptidase C39 family protein (S0112_034), alpha/beta hydrolase (S0112_060), and putative endolysin (S0112_062). Among this module, the tape measure protein usually exists in siphoviruses and is involved in the assembly of phage tails and the determination of the tail length [44]. Endolysin is a lytic enzyme used to lyse the cell wall of the host. The endolysin consists of two domains: One is the N-terminal domain, which is responsible for degrading bacterial peptidoglycan, and the other is the C-terminal domain, called holin, which is mainly bound to the cell membrane [45].
The DNA replication/recombination/metabolism module in phage genome includes ten ORFs, whose gene products include DNA polymerase I (S0112_092), DNA polymerase subunit B (S0112_093), DNA/RNA helicase (S0112_005), primase/helicase (S0112_018), DNA-binding domain protein (S0112_017), exonuclease (S0112_006), endonuclease (S0112_036), ATPase (S0112_009), P-loop containing nucleoside triphosphate hydrolase (S0112_011), and DprA-like protein (S0112_023). Among them, DNA polymerase is vital for phage DNA replication, and its coding gene is relatively conservative in many tailed phages [46]. Primase and helicase found in the genome of phage S0112 are of key importance for DNA replication. The primosome formed by primase and helicase can synthesize a short RNA primer, which can be elongated by DNA polymerase [47]. Exonucleases are enzymes that participate in the process of breaking phosphodiester bonds at either 3′ or 5′ end [48]. Other enzymes including endonuclease, DNA-binging protein, and DprA-like protein play important roles in DNA repair and transcription regulation. It is worth noting that ATPase encoded by ORF 9 is similar to the ATPase of the benthic Vibrio phage VpKK5, which was also isolated from coastal sediments. ATPase is important for DNA replication, DNA recombination, and DNA packaging [18].
Interestingly, we found ORF 99 in the phage genome of S0112 is an auxiliary metabolic gene encoding the MazG nucleotide pyrophosphohydrolase domain [49]. MazG was implicated as a regulator of programmed cell death in Escherichia coli by interfering with the MazEF toxin–antitoxin system through reducing the intracellular guanosine 3′,5′-bispyrophosphate (ppGpp) levels. The small nucleotide ppGpp, as a global regulator of gene expression, whose accumulation can cause the reorientation of bacterial transcription, and only genes important for starvation survival can be expressed [50,51]. Given the role of MazG in E. coli, it was proposed that phages containing mazG gene may trick the host to mimic a nutrient replete cell state by reducing the accumulation of ppGpp in the host cells, so that phages can utilize the host RNA transcription system to amplify the progeny phages [52]. However, a recent study found that the purified MazG from cyanophage has no binding or hydrolysis activity to ppGpp [49]. Instead, cyanophage encoded MazG protein can bind and hydrolyze dGTP and dCTP deoxyribonucleotides, and it was speculated that it may indirectly regulate the host stringent response by hydrolyzing other nucleotide substrates [49]. Here, the specific function of MazG gene in the phage S0112 remains to be studied.
Another feature of phage S0112 is the presence of a gene cluster consisting of five functional genes encoding QueC (7-cyano-7-deazaguanine synthase), QueD (6-carboxy-5,6,7,8-tetrahydropterin synthase), FolE (GTP cyclohydrolase I), queuosine tRNA-ribosyltransferase, and organic radical activating enzyme. This kind of gene cluster has also been observed in other phages, such as Rhizobium phage RHEph04 [53], R5C [54], Vibrio phage VpKK5 [18], Escherichia coli phage 9g [55], and the genus Seuratvirus [56]. This gene cluster involves in the biosynthesis of a precursor, i.e., 7-cyano-7-deazaguanine (preQ0), of queuosine or archaeosine [57]. Queuosine is a well-known nucleoside derivative that modifies tRNAs by replacement of guanine at certain position [58]. Queuosine modification in bacteria was considered a way to improve reading frame maintenance [59], but this opinion is still controversial, as a recent study did not observe this function in Escherichia coli [60]. Meanwhile, a phage infecting Escherichia coli that contains homologues of queuosine biosynthesis genes, including organic radical activating enzyme, queuosine tRNA-ribosyltransferase, QueC, QueD, FolE, and (GAT)-QueC (glutamine amidotransferase class-II) was found to be resistant to a wide range of restriction endonucleases, by modifying its DNA (converting deoxyguanosine to 2′-dexoy-archaeosine) [55,56,57]. As reported, phages belonging to the genus Seuratvirus contain QueC, QueD, and QueE, as well as FolE, modify their DNA with queuosine, as they lack homologues of (GAT)-QueC, which is required for the insertion of archaeosine [56]. Our analysis showed that phage S0112 contains only QueC homologues without the homologues of (GAT)-QueC, and we propose that phage S0112 may modify their DNA in a similar manner to the genus Seuratvirus. Further work is needed to verify the biological role of the queuosine biosynthesis genes in phage S0112.

3.5. Mass-Spectrometric Identification of Phage Proteins

LC–MS/MS analysis identified 16 out of the 102 predicted ORF expression products, among which, eleven are phage structural proteins, including major capsid protein, minor tail protein, tail length tape measure protein, tail assembly structure protein, tail terminator protein, head-tail connector protein, and virion structural protein (Table 4). Another five proteins (encoded by ORF1, 65, 68, 69, and 80), having no significant homologs in the database, were classified as uncharacterized proteins. However, other predicted ORF gene products, including the functional proteins, were not detected via LC–MS/MS analysis. The reason might be that the expression level of these proteins is relatively low and did not reach the detection limit of LC–MS/MS in our experiment.

3.6. Phylogenetic Analysis

We used the gene-encoding DNA polymerase subunit B, which has been widely used as a marker gene to study phage diversity to determine the phylogeny of phage S0112. Here, the phylogenetic analysis based on phage DNA polymerase subunit B gene and the viral proteomic tree constructed by the ViPTree both indicated that the closest relative of phage S0112 was Vibrio phage VpKK5 (Figure 7 and Figure 8), while, in the whole genome level, phage S0112 and VpKK5 share very low similarity with each other (Figure S1). Phylogenetic analysis combined with morphological assessment indicated that phage S0112 could be a representative of a novel genus belonging to the family Siphoviridae.

3.7. Environmental Distribution

By using the phage genome of S0112 as queries to search against the IMG/VR database, which is the largest available database of viral genomics, with parameters described elsewhere (blastn -db ~/IMG/VR -query queryname -outfmt “6 std qlen slen” -evalue 1e-5 -max_target_seqs 1000) [61], we did not find any homologous contig and genome, further indicating the novelty of phage S0112. Reciprocal best-hit BLAST showed that the phage proteins have very low similarity to those in both POV and GOS databases, revealing that this kind of phage may rarely exists in other environments beside sediment (Figure S2).

4. Conclusions

The benthic phage S0112 infecting Shewanella indica KJW27 is a novel lytic phage that was isolated from coastal sediments and can be considered a representative of a novel genus under the family Siphoviridae. Fifty-nine percent of its ORFs are functionally unknown. Its strong replication ability, as reflected by the high burst size (i.e., 1170), and its relatively longer latency might be adaptive survival strategies for coastal sediment environments. Several auxiliary metabolic genes were observed in the phage genome, i.e., the pyrophosphohydrolase (MazG) gene and a gene cluster that may be involved in 7-cyano-7-deazaguanine (preQ0) biosynthesis. To the best of our knowledge, this study represents the first report of a benthic phage infecting Shewanella and sheds light on the phage-host interactions in coastal sediments. As, so far, very few phages have been isolated from coastal sediments, which largely limits our understanding of their biological characteristics, the isolation and characterization of more benthic phages infecting different hosts in coastal sediments are vital for unveiling the phage–host interactions in coastal sediments and their ecological significance, as well as the largely unknown metagenomic dark matters of sediment phages.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1999-4915/11/11/1081/s1. Table S1: Summary of latent period and burst size of sixty-eight phages belonging to the order Caudovirales. Table S2: Effects of rifampicin on the phage progeny production. Table S3: The Shewanella sp. phage S0112 genome annotations (MK675901). Figure S1: Mauve alignment of phage S0112 genome against phage VpKK5 genome. Figure S2: Prevalence of phage S0112-like ORFs in Pacific Ocean Virome (POV) and Global Ocean Survey (GOS) metagenomic datasets, using the reciprocal best-BLAST strategy.

Author Contributions

Y.Z. designed the experiments and revised the manuscript; J.Z., L.W., J.L., L.Z., and C.L. helped to analyze the data and revise the manuscript; Z.W. performed the experiments, analyzed the data, and wrote the manuscripts.

Funding

This research was funded by the open task of Qingdao National Laboratory for Marine Science and Technology (QNLM2016ORP0311), Natural Science Foundation of China (no. 41876174), the senior user’s project of RV KEXUE (KEXUE2019GZ03), DICP&QIBEBT (Grant: DICP&QIBEBT UN201803), QIBEBT (Grant: QIBEBT ZZBS201805), Dalian National Laboratory for Clean Energy (DNL), CAS. This study is a contribution to the international IMBER project.

Acknowledgments

We thank Zhenzhen Hu, Shailesh Nair of Qingdao Institute of Bioenergy and Bioprocess Technology for their useful suggestions and help, and the support from the Center for Ocean Mega-Science, Chinese Academy of Sciences.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Fuhrman, J.A. Marine viruses and their biogeochemical and ecological effects. Nature 1999, 399, 541. [Google Scholar] [CrossRef]
  2. Suttle, C.A. Marine viruses-major players in the global ecosystem. Nat. Rev. Microbiol. 2007, 5, 801–812. [Google Scholar] [CrossRef]
  3. Suttle, C.A. Viruses in the sea. Nature 2005, 437, 356–361. [Google Scholar] [CrossRef]
  4. Zhang, Y.; Huang, C.; Yang, J.; Jiao, N. Interactions between marine microorganisms and their phages. Chin. Sci. Bull. 2011, 56, 1770–1777. [Google Scholar] [CrossRef]
  5. Zhao, Y.; Temperton, B.; Thrash, J.C.; Schwalbach, M.S.; Vergin, K.L.; Landry, Z.C.; Ellisman, M.; Deerinck, T.; Sullivan, M.B.; Giovannoni, S.J. Abundant SAR11 viruses in the ocean. Nature 2013, 494, 357–360. [Google Scholar] [CrossRef]
  6. Zeidner, G.; Bielawski, J.P.; Shmoish, M.; Scanlan, D.J.; Sabehi, G.; Beja, O. Potential photosynthesis gene recombination between Prochlorococcus and Synechococcus via viral intermediates. Environ. Microbiol. 2005, 7, 1505–1513. [Google Scholar] [CrossRef]
  7. Glud, R.N.; Mathias, M. Virus and bacteria dynamics of a coastal sediment: Implication for benthic carbon cycling. Limnol. Oceanogr. 2004, 49, 2073–2081. [Google Scholar] [CrossRef]
  8. Breitbart, M.; Felts, B.; Kelley, S.; Mahaffy, J.M.; Nulton, J.; Salamon, P.; Rohwer, F. Diversity and population structure of a near-shore marine-sediment viral community. Proc. Biol. Sci. 2004, 271, 565–574. [Google Scholar] [CrossRef] [PubMed]
  9. Hau, H.H.; Gralnick, J.A. Ecology and biotechnology of the genus Shewanella. Annu. Rev. Microbiol. 2007, 61, 237–258. [Google Scholar] [CrossRef] [PubMed]
  10. Fredrickson, J.K.; Romine, M.F.; Beliaev, A.S.; Auchtung, J.M.; Driscoll, M.E.; Gardner, T.S.; Nealson, K.H.; Osterman, A.L.; Pinchuk, G.; Reed, J.L.; et al. Towards environmental systems biology of Shewanella. Nat. Rev. Microbiol. 2008, 6, 592–603. [Google Scholar] [CrossRef] [PubMed]
  11. Li, B.; Cheng, Y.; Fan, Y.; Liu, D.; Fang, C.; Wu, C.; Li, W.; Yang, Z.; Yu, H. Estimates of abundance and diversity of Shewanella genus in natural and engineered aqueous environments with newly designed primers. Sci. Total Environ. 2018, 637–638, 926–933. [Google Scholar] [CrossRef] [PubMed]
  12. Yang, Z.; Tao, X.; Zhang, H.; Rao, S.; Gao, L.; Pan, Z.; Jiao, X. Isolation and characterization of virulent phages infecting Shewanella baltica and Shewanella putrefaciens, and their application for biopreservation of chilled channel catfish (Ictalurus punctatus). Int. J. Food Microbiol. 2019, 292, 107–117. [Google Scholar] [CrossRef] [PubMed]
  13. Han, F.; Li, M.; Lin, H.; Wang, J.; Cao, L.; Khan, M.N. The novel Shewanella putrefaciens-infecting bacteriophage Spp001: Genome sequence and lytic enzymes. J. Ind. Microbiol. Biotechnol. 2014, 41, 1017–1026. [Google Scholar] [CrossRef] [PubMed]
  14. Leigh, B.; Karrer, C.; Cannon, J.P.; Breitbart, M.; Dishaw, L.J. Isolation and characterization of a Shewanella phage-host system from the gut of the tunicate, Ciona intestinalis. Viruses 2017, 9, 60. [Google Scholar] [CrossRef]
  15. Sencčilo, A.; Luhtanen, A.M.; Saarijӓrvi, M.; Bamford, D.H.; Roine, E. Cold-active bacteriophages from the Baltic Sea ice have diverse genomes and virus-host interactions. Environ. Microbiol. 2015, 17, 3628–3641. [Google Scholar] [CrossRef]
  16. Zhao, J.; Wang, L.; Tang, L.; Ren, R.; You, W.; Farooq, R.; Wang, Z.; Zhang, Y. Changes in bacterial community structure and humic acid composition in response to enhanced extracellular electron transfer process in coastal sediment. Arch. Microbiol. 2019, 201, 897–906. [Google Scholar] [CrossRef]
  17. Zhao, J.S.; Manno, D.; Beaulieu, C.; Paquet, L.; Hawari, J. Shewanella sediminis sp. nov., a novel Na+-requiring and hexahydro-1,3,5-trinitro-1,3,5-triazine-degrading bacterium from marine sediment. Int. J. Syst. Evol. Microbiol. 2005, 55, 1511–1520. [Google Scholar] [CrossRef]
  18. Lal, T.M.; Sano, M.; Ransangan, J. Genome characterization of a novel vibriophage VpKK5 (Siphoviridae) specific to fish pathogenic strain of Vibrio parahaemolyticus. J. Basic Microbiol. 2016, 56, 872–888. [Google Scholar] [CrossRef]
  19. Yang, H.; Liang, L.; Lin, S.; Jia, S. Isolation and characterization of a virulent bacteriophage AB1 of Acinetobacter baumannii. BMC Microbiol. 2010, 10, 131. [Google Scholar] [CrossRef]
  20. Verma, P.; Pandey, P.K.; Gupta, A.K.; Kim, H.J.; Baik, K.S.; Seong, C.N.; Patole, M.S.; Shouche, Y.S. Shewanella indica sp. nov., isolated from sediment of the Arabian Sea. Int. J. Syst. Evol. Microbiol. 2011, 61, 2058–2064. [Google Scholar] [CrossRef]
  21. Zhang, Y.; Jiao, N. Roseophage RDJL Ф1, infecting the aerobic anoxygenic phototrophic bacterium Roseobacter denitrificans OCh114. Appl. Environ. Microbiol. 2009, 75, 1745–1749. [Google Scholar] [CrossRef] [PubMed]
  22. Weinbauer, M.G.; Rowe, J.M.; Wilhelm, S. Determining rates of virus production in aquatic systems by the virus reduction approach. In Manual of Aquatic Viral Ecology; Suttle, C.A., Wilhelm, S.W., Weinbauer, M.G., Eds.; American Society of Limnology and Oceanography: Waco, TX, USA, 2010; pp. 1–8. [Google Scholar]
  23. Alonso, M.D.C.; Rodríguez, J.; Borrego, J.J. Characterization of marine bacteriophages isolated from the Alboran Sea (Western Mediterranean). J. Plankton Res. 2002, 24, 1079–1087. [Google Scholar] [CrossRef]
  24. Jurczak-Kurek, A.; Gasior, T.; Nejman-Falenczyk, B.; Bloch, S.; Dydecka, A.; Topka, G.; Necel, A.; Jakubowska-Deredas, M.; Narajczyk, M.; Richert, M.; et al. Biodiversity of bacteriophages: Morphological and biological properties of a large group of phages isolated from urban sewage. Sci. Rep. 2016, 6, 34338. [Google Scholar] [CrossRef] [PubMed]
  25. Middelboe, M.; Chan, A.M.; Bertelsen, S.K. Isolation and life cycle characterization of lytic viruses infecting heterotrophic bacteria and cyanobacteria. In Manual of Aquatic Viral Ecology; Wilhelm, S.W., Weinbauer, M.G., Suttle, C.A., Eds.; American Society of Limnology and Oceanography: Waco, TX, USA, 2010; pp. 118–133. [Google Scholar]
  26. Anderson, T.F. Bacteriophages. Annu. Rev. Microbiol. 1950, 4, 21–34. [Google Scholar] [CrossRef]
  27. Lu, L.; Cai, L.; Jiao, N.; Zhang, R. Isolation and characterization of the first phage infecting ecologically important marine bacteria Erythrobacter. Virol. J. 2017, 14, 104. [Google Scholar] [CrossRef]
  28. Matsui, T.; Yoshikawa, G.; Mihara, T.; Chatchawankanphanich, O.; Kawasaki, T.; Nakano, M.; Fujie, M.; Ogata, H.; Yamada, T. Replications of two closely related groups of Jumbo phages show different level of dependence on host-encoded RNA polymerase. Front. Microbiol. 2017, 8, 1010. [Google Scholar] [CrossRef]
  29. Zhang, X.; Wang, Y.; Li, S.; An, X.; Pei, G.; Huang, Y.; Fan, H.; Mi, Z.; Zhang, Z.; Wang, W.; et al. A novel termini analysis theory using HTS data alone for the identification of Enterococcus phage EF4-like genome termini. BMC Genom. 2015, 16, 414. [Google Scholar] [CrossRef]
  30. Garneau, J.; Depardieu, F.; Fortier, L.C.; Bikard, D.; Monot, M. PhageTerm: A fast and user-friendly software to determine bacteriophage termini and packaging mode using randomly fragmented NGS data. bioRxiv 2017, 108100. [Google Scholar] [CrossRef]
  31. Schattner, P.; Brooks, A.N.; Lowe, T.M. The tRNAscan-SE, snoscan and snoGPS web servers for the detection of tRNAs and snoRNAs. Nucl. Acids Res. 2005, 33, W686–W689. [Google Scholar] [CrossRef]
  32. Stajich, J.E.; Darling, A.E.; Mau, B.; Perna, N.T. progressiveMauve: Multiple genome alignment with gene gain, loss and rearrangement. PLoS ONE 2010, 5, e11147. [Google Scholar]
  33. Paez-Espino, D.; Roux, S.; Chen, I.A.; Palaniappan, K.; Ratner, A.; Chu, K.; Huntemann, M.; Reddy, T.B.K.; Pons, J.C.; Llabrés, M.; et al. IMG/VR v.2.0: An integrated data management and analysis system for cultivated and environmental viral genomes. Nucl. Acids Res. 2018, 47, D678–D686. [Google Scholar] [CrossRef] [PubMed]
  34. Grissa, I.; Vergnaud, G.; Pourcel, C. CRISPRcompar: A website to compare clustered regularly interspaced short palindromic repeats. Nucl. Acids Res. 2008, 36, W145–W148. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Wisniewski, J.R.; Zougman, A.; Nagaraj, N.; Mann, M. Universal sample preparation method for proteome analysis. Nat. Methods 2009, 6, 359–362. [Google Scholar] [CrossRef] [PubMed]
  36. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Nishimura, Y.; Yoshida, T.; Kuronishi, M.; Uehara, H.; Ogata, H.; Goto, S. ViPTree: The viral proteomic tree server. Bioinformatics 2017, 33, 2379–2380. [Google Scholar] [CrossRef] [PubMed]
  38. Hurwitz, B.L.; Sullivan, M.B. The Pacific Ocean Virome (POV): A marine viral metagenomic dataset and associated protein clusters for quantitative viral ecology. PLoS ONE 2013, 8, e57355. [Google Scholar] [CrossRef]
  39. Rusch, D.B.; Halpern, A.L.; Sutton, G.; Heidelberg, K.B.; Williamson, S.; Yooseph, S.; Wu, D.; Eisen, J.A.; Hoffman, J.M.; Remington, K.; et al. The sorcerer II global ocean sampling expedition: Northwest Atlantic through eastern tropical Pacific. PLoS Biol. 2007, 5, e77. [Google Scholar] [CrossRef]
  40. Abedon, S.T. Selection for bacteriophage latent period length by bacterial density: A theoretical examination. Microb. Ecol. 1989, 18, 79–88. [Google Scholar] [CrossRef]
  41. Børsheim, K.Y. Native marine bacteriophages. FEMS Microbiol. Ecol. 1993, 11, 141–159. [Google Scholar] [CrossRef]
  42. Mei, M.L.; Danovaro, R. Virus production and life strategies in aquatic sediments. Limnol. Oceanogr. 2004, 49, 459–470. [Google Scholar] [CrossRef]
  43. Hamdi, S.; Rousseau, G.M.; Labrie, S.J.; Tremblay, D.M.; Kourda, R.S.; Ben Slama, K.; Moineau, S. Characterization of two polyvalent phages infecting Enterobacteriaceae. Sci. Rep. 2017, 7, 40349. [Google Scholar] [CrossRef] [PubMed]
  44. Mahony, J.; Alqarni, M.; Stockdale, S.; Spinelli, S.; Feyereisen, M.; Cambillau, C.; Sinderen, D.V. Functional and structural dissection of the tape measure protein of lactococcal phage TP901-1. Sci. Rep. 2016, 6, 36667. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Smith, D.L.; Struck, D.K.; Scholtz, J.M.; Young, R. Purification and biochemical characterization of the lambda holin. J. Bacteriol. 1998, 180, 2531–2540. [Google Scholar] [PubMed]
  46. Filée, J.; Forterre, P.; Sen-Lin, T.; Laurent, J. Evolution of DNA polymerase families: Evidences for multiple gene exchange between cellular and viral proteins. J. Mol. Evol. 2002, 54, 763–773. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Matson, S.W.; Bean, D.W.; George, J.W. DNA helicases: Enzymes with essential roles in all aspects of DNA metabolism. Bioessays 1994, 16, 13–22. [Google Scholar] [CrossRef] [PubMed]
  48. Subramanian, K.; Rutvisuttinunt, W.; Scott, W.; Myers, R.S. The enzymatic basis of processivity in lambda exonuclease. Nucl. Acids Res. 2003, 31, 1585–1596. [Google Scholar] [CrossRef] [Green Version]
  49. Rihtman, B.; Bowman-Grahl, S.; Millard, A.; Corrigan, R.M.; Clokie, M.R.J.; Scanlan, D.J. Cyanophage MazG is a pyrophosphohydrolase but unable to hydrolyse magic spot nucleotides. Environ. Microbiol. Rep. 2019, 11, 448–455. [Google Scholar] [CrossRef] [Green Version]
  50. Magnusson, L.U.; Farewell, A.; Nystrom, T. ppGpp: A global regulator in Escherichia coli. Trends Microbiol. 2005, 13, 236–242. [Google Scholar] [CrossRef]
  51. Gross, M.; Marianovsky, I.; Glaser, G. MazG—A regulator of programmed cell death in Escherichia coli. Mol. Microbiol. 2006, 59, 590–601. [Google Scholar] [CrossRef]
  52. Clokie, M.R.; Mann, N.H. Marine cyanophages and light. Environ. Microbiol. 2006, 8, 2074–2082. [Google Scholar] [CrossRef]
  53. Santamaría, R.I.; Bustos, P.; Sepúlveda-Robles, O.; Lozano, L.; Rodríguez, C.; Fernández, J.L.; Juárez, S.; Kameyama, L.; Guarneros, G.; Dávila, G.; et al. Narrow-host-range bacteriophages that infect Rhizobium etli associate with distinct genomic types. Appl. Environ. Microbiol. 2014, 80, 446–454. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Yang, Y.; Cai, L.; Ma, R.; Xu, Y.; Tong, Y.; Huang, Y.; Jiao, N.; Zhang, R. A novel roseosiphophage isolated from the oligotrophic South China Sea. Viruses 2017, 9, 109. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Kulikov, E.E.; Golomidova, A.K.; Letarova, M.A.; Kostryukova, E.S.; Zelenin, A.S.; Prokhorov, N.S.; Letarov, A.V. Genomic sequencing and biological characteristics of a novel Escherichia coli bacteriophage 9g, a putative representative of a new Siphoviridae genus. Viruses 2014, 6, 5077–5092. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Sazinas, P.; Redgwell, T.; Rihtman, B.; Grigonyte, A.; Michniewski, S.; Scanlan, D.J.; Hobman, J.; Millard, A. Comparative genomics of bacteriophage of the genus Seuratvirus. Genome Biol. Evol. 2018, 10, 72–76. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Thiaville, J.J.; Kellner, S.M.; Yuan, Y.; Hutinet, G.; Thiaville, P.C.; Jumpathong, W.; Mohapatra, S.; Brochier-Armanet, C.; Letarov, A.V.; Hillebrand, R.; et al. Novel genomic island modifies DNA with 7-deazaguanine derivatives. Proc. Natl. Acad. Sci. USA 2016, 113, E1452–E1459. [Google Scholar] [CrossRef] [Green Version]
  58. Yuan, Y.; Hutinet, G.; Valera, J.G.; Hillebrand, R.; Gustafson, A.; Iwata-Reuyl, D.; Dedon, P.C.; de Crécy-Lagard, V. Identification of the minimal bacterial 2′-deoxy-7-amido-7-deazaguanine synthesis machinery. Mol. Microbiol. 2018, 110, 469–483. [Google Scholar] [CrossRef]
  59. Urbonavičius, J.; Qian, Q.; Durand, J.M.; Hagervall, T.G.; Björk, G.R. Improvement of reading frame maintenance is a common function for several tRNA modifications. EMBO. J. 2001, 20, 4863–4873. [Google Scholar] [CrossRef] [Green Version]
  60. Xu, J.; Chen, M.; He, L.; Zhang, S.; Ding, T.; Yao, H.; Lu, C.; Zhang, W. Isolation and characterization of a T4-like phage with a relatively wide host range within Escherichia coli. J. Basic Microbiol. 2016, 56, 405–421. [Google Scholar] [CrossRef]
  61. Toro, N.; Mestre, M.R.; Martinez-Abarca, F.; Gonzalez-Delgado, A. Recruitment of reverse transcriptase-cas1 fusion proteins by type VI-A CRISPR-Cas Systems. Front. Microbiol. 2019, 10, 2160. [Google Scholar] [CrossRef]
Figure 1. (a) Phage plaques formed in double-layer agar plates and (b) transmission electron microscopy image of the phage S0112.
Figure 1. (a) Phage plaques formed in double-layer agar plates and (b) transmission electron microscopy image of the phage S0112.
Viruses 11 01081 g001
Figure 2. Kinetics of lytic development of Shewanella phage S0112 in Shewanella indica. Results are shown as (a) bacterial culture density measured at OD600, (b) the number of host cells after phage S0112 infection per 1 mL (CFU·mL−1), and (c) the number of phages per 1 mL (PFU·mL−1). Error bars show standard deviations among triplicate samples.
Figure 2. Kinetics of lytic development of Shewanella phage S0112 in Shewanella indica. Results are shown as (a) bacterial culture density measured at OD600, (b) the number of host cells after phage S0112 infection per 1 mL (CFU·mL−1), and (c) the number of phages per 1 mL (PFU·mL−1). Error bars show standard deviations among triplicate samples.
Viruses 11 01081 g002
Figure 3. One-step growth curve of Shewanella phage S0112. Open circles represent non-chloroform-treated samples; closed circles represent chloroform-treated samples. Error bars show standard deviations among triplicate samples.
Figure 3. One-step growth curve of Shewanella phage S0112. Open circles represent non-chloroform-treated samples; closed circles represent chloroform-treated samples. Error bars show standard deviations among triplicate samples.
Viruses 11 01081 g003
Figure 4. The latent period of sixty-eight phages infecting Roseobacter spp., Vibrio spp., Pseudoalteromonas spp., Aeromonas spp., Acinetobacter spp., Escherichia spp., etc.
Figure 4. The latent period of sixty-eight phages infecting Roseobacter spp., Vibrio spp., Pseudoalteromonas spp., Aeromonas spp., Acinetobacter spp., Escherichia spp., etc.
Viruses 11 01081 g004
Figure 5. Multiple genome alignments of the whole genomes of all Shewanella phages, using Mauve algorithm. The Shewanella phages in the analysis are isolated from different environments, including Baltic Sea ice, seawater, sewage, waste effluents from fish market, and coastal sediments. Colored boxes indicate homologous DNA regions between genomes without genomic rearrangements. Genomic similarity is represented by the height of the bars, which corresponds to the average level of conservation in that region of the genome sequence. Completely white regions represent fragments that were not aligned or contained sequence elements specific to a particular genome.
Figure 5. Multiple genome alignments of the whole genomes of all Shewanella phages, using Mauve algorithm. The Shewanella phages in the analysis are isolated from different environments, including Baltic Sea ice, seawater, sewage, waste effluents from fish market, and coastal sediments. Colored boxes indicate homologous DNA regions between genomes without genomic rearrangements. Genomic similarity is represented by the height of the bars, which corresponds to the average level of conservation in that region of the genome sequence. Completely white regions represent fragments that were not aligned or contained sequence elements specific to a particular genome.
Viruses 11 01081 g005
Figure 6. Genomic map of phage S0112. The phage S0112 genome is presented schematically with ORFs indicated by leftward or rightward oriented arrows, based on the direction of transcription. Gene features (DNA packaging, structure/morphogenesis protein, DNA replication/recombination/metabolism, and queuosine biosynthesis) are color-coded according to the legend below the figure.
Figure 6. Genomic map of phage S0112. The phage S0112 genome is presented schematically with ORFs indicated by leftward or rightward oriented arrows, based on the direction of transcription. Gene features (DNA packaging, structure/morphogenesis protein, DNA replication/recombination/metabolism, and queuosine biosynthesis) are color-coded according to the legend below the figure.
Viruses 11 01081 g006
Figure 7. The maximum likelihood phylogenetic tree of phages based on their DNA polymerase subunit B amino acid sequences. The red star represents Shewanella phage S0112 isolated in this study. The reference sequences were collected from the NCBI database. The tree was constructed based on the MUSCLE alignment of MEGA 7.0. The bootstrap values were based on 1000 replicates.
Figure 7. The maximum likelihood phylogenetic tree of phages based on their DNA polymerase subunit B amino acid sequences. The red star represents Shewanella phage S0112 isolated in this study. The reference sequences were collected from the NCBI database. The tree was constructed based on the MUSCLE alignment of MEGA 7.0. The bootstrap values were based on 1000 replicates.
Viruses 11 01081 g007
Figure 8. The viral proteomic tree including phage S0112 and other closest relative phages constructed using the ViPTree. The red star represents the Shewanella phage S0112.
Figure 8. The viral proteomic tree including phage S0112 and other closest relative phages constructed using the ViPTree. The red star represents the Shewanella phage S0112.
Viruses 11 01081 g008
Table 1. Infection specificities of the phage S0112 against eighteen bacterial strains.
Table 1. Infection specificities of the phage S0112 against eighteen bacterial strains.
Genus and SpeciesStrainSourceStrains Lysed by
Phage S0112
Shewanella indicaKJW27 asediment+
Shewanella basaltisCJW-54sediment-
Shewanella chilikensisJC5sediment-
Shewanella japonicaKCTC22435sediment-
Shewanella algaeJCM 21037seawater-
Alginatibacterium sediminisALS 81sediment-
Woeseia oceaniSDUM189001sediment-
Sediminicola luteusSDUM701001sediment-
Kordiimonas sediminisN39sediment-
Roseobacter denitrificansOCH114seawater-
Marinobacter vinifirmusD7035wastewater-
Halomonas denitrificansD7027saline water-
Vibrio alginolyticusCIP 82.01seawater-
Vibrio neocaledonicusNC 470seawater-
Vibrio azureusNBRC 104587seawater-
Vibrio harveyiNBRC 15634seawater-
Vibrio parahaemolyticusNBRC 12711seawater-
Vibrio campbelliiCAIM 519seawater-
+, Lysed; -, not lysed. a Host strain used in this study.
Table 2. Stability of Shewanella phage S0112 treated with various physical and chemical agents.
Table 2. Stability of Shewanella phage S0112 treated with various physical and chemical agents.
Phage Name−20 °C15 °C20 °C28 °C30 °C37 °C40 °C50 °C60 °C95 °C100 °CpH 2pH 4pH 10pH 12Osmotic Shock0.1% CTAB0.09% SDS0.1% Sarkosyl63% Ethanol90% Acetone50% DMSOChloroform
phage S0112951.34.599.210086.772.91.40.040006.882.10.330.452.5012.71.80.430.991.9
Comparison of the effects of different physical and chemical agents on the survival of phages S0112. Percent of surviving phages under various conditions is shown.
Table 3. The genome information of ten Shewanella phages.
Table 3. The genome information of ten Shewanella phages.
Shewanella Phages aHost StrainsFamilySize (Kb)GC%ProteinNo. of tRNAsIsolation Source
S0112Shewanella indicaSiphoviridae62.28644.71020Coastal sediment
Spp001S. putrefaciensSiphoviridae54.78949.4670Sewage
3/49S. balticaSiphoviridae40.16142700Baltic Sea ice
1/44S. frigidimarinaSiphoviridae49.6439.8750Baltic Sea ice
SppYZU05S. putrefaciensSiphoviridae54.31950.63650Waste effluents
1/41S. balticaMyoviridae43.5142.7690Baltic Sea ice
SFCi1S. fidelisMyoviridae42.27959.1400Seawater
SppYZU01S. balticaMyoviridae43.56755.72490Waste effluents
1/4S. frigidimarinaMyoviridae133.82436.92353Baltic Sea ice
1/40S. balticaMyoviridae139.00436.92363Baltic Sea ice
a Phage S0112 was obtained in the present study. Data for the other nine phages are from the references [12,13,14,15].
Table 4. Mass spectrometry analysis of phage S0112 virion.
Table 4. Mass spectrometry analysis of phage S0112 virion.
Detected ProteinPredicted FunctionMolecular Mass (kDa)Number of PeptidesSequence Coverage (%)Protein Score
S0112_001Uncharacterized protein11.2333.7175
S0112_065Uncharacterized protein52.61746.73837
S0112_066Head-tail connector protein19.6749.71873
S0112_068Uncharacterized protein16.9530.1575
S0112_069Uncharacterized protein43.149.4220
S0112_071Major capsid protein37.31342.19723
S0112_074Tail completion protein13.7443.7418
S0112_075Tail terminator protein16.1644.3694
S0112_076Major capsid protein54.11758.715,158
S0112_079Tail length tape measure protein992431.45007
S0112_080Uncharacterized protein38.4527.11859
S0112_081Virion structural protein30.8730.8709
S0112_082Tail assembly structural protein62.2827.11781
S0112_083Minor tail protein30.8837.71024
S0112_086Tail assembly structural protein83.31223.91972
S0112_087Minor tail protein14.9445.5696

Share and Cite

MDPI and ACS Style

Wang, Z.; Zhao, J.; Wang, L.; Li, C.; Liu, J.; Zhang, L.; Zhang, Y. A Novel Benthic Phage Infecting Shewanella with Strong Replication Ability. Viruses 2019, 11, 1081. https://0-doi-org.brum.beds.ac.uk/10.3390/v11111081

AMA Style

Wang Z, Zhao J, Wang L, Li C, Liu J, Zhang L, Zhang Y. A Novel Benthic Phage Infecting Shewanella with Strong Replication Ability. Viruses. 2019; 11(11):1081. https://0-doi-org.brum.beds.ac.uk/10.3390/v11111081

Chicago/Turabian Style

Wang, Zengmeng, Jiulong Zhao, Long Wang, Chengcheng Li, Jianhui Liu, Lihua Zhang, and Yongyu Zhang. 2019. "A Novel Benthic Phage Infecting Shewanella with Strong Replication Ability" Viruses 11, no. 11: 1081. https://0-doi-org.brum.beds.ac.uk/10.3390/v11111081

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