ORIGINAL RESEARCH article
Volume 7 - 2016 | https://doi.org/10.3389/fmicb.2016.01807
Comparative Functional Genomic Analysis of Two Vibrio Phages Reveals Complex Metabolic Interactions with the Host Cell
- 1Laboratory of Molecular Biology, Department of Biotechnology, School of Food, Biotechnology and Development, Agricultural University of Athens, Athens, Greece
- 2Institute of Marine Biology, Biotechnology and Aquaculture, Hellenic Centre for Marine Research, Heraklion, Crete, Greece
- 3Marine Biological Section, University of Copenhagen, Helsingør, Denmark
Sequencing and annotation was performed for two large double stranded DNA bacteriophages, φGrn1 and φSt2 of the Myoviridae family, considered to be of great interest for phage therapy against Vibrios in aquaculture live feeds. In addition, phage–host metabolic interactions and exploitation was studied by transcript profiling of selected viral and host genes. Comparative genomic analysis with other large Vibrio phages was also performed to establish the presence and location of homing endonucleases highlighting distinct features for both phages. Phylogenetic analysis revealed that they belong to the “schizoT4like” clade. Although many reports of newly sequenced viruses have provided a large set of information, basic research related to the shift of the bacterial metabolism during infection remains stagnant. The function of many viral protein products in the process of infection is still unknown. Genome annotation identified the presence of several viral open reading frames (ORFs) participating in metabolism, including a Sir2/cobB (sirtuin) protein and a number of genes involved in auxiliary NAD+ and nucleotide biosynthesis, necessary for phage DNA replication. Key genes were subsequently selected for detail study of their expression levels during infection. This work suggests a complex metabolic interaction and exploitation of the host metabolic pathways and biochemical processes, including a possible post-translational protein modification, by the virus during infection.
The ever growing demand for fishery products and seafood has led to intensification of aquaculture. The overuse and abuse of antibiotics has resulted in the selection of resistant bacteria that are considered one of the biggest risk for humanity. Bacteriophage therapy has been suggested as a potential alternative method for both treatment and prophylaxis of bacterial infections, including aquaculture, showing very promising results (Stone, 2002; Sulakvelidze, 2011; Jassim and Limoges, 2014). The advancement of sequencing technology has boosted the genomic characterization of isolated phages providing fascinating insights to their biology and interaction with their host. Furthermore, genomic analysis may provide crucial information for their safety as potential therapeutics as it is of imperative importance to exclude temperate phages or phages carrying toxin or antibiotic resistance genes (Muniesa et al., 2004; Modi et al., 2013).
Although the T4 bacteriophage’s genome has been fully sequenced, annotated, and characterized for many years now (Miller et al., 2003b), the need for T4-like phages genomic characterization from a varied host range remains high, with the number of newly isolated bacteriophages and their genomic characterization increasing (Klumpp et al., 2012).
By using protein sequence similarity in the past, large Vibrio phages have been distinct from the T4-like phages. Known as “schizoT4like” they have been also proposed to be classified as KVP40-like phages (Lavigne et al., 2009) mainly because of the size of the head and the high genome and protein similarities with the KVP40 phage (Miller et al., 2003a). According to GenBank this group now contains the Vibrio phages KVP40, φpp2 and nt-1 all of which seem to be both morphologically and genetically similar, with homing endonucleases (HEs), NAD+ biosynthesis and nucleotide metabolism related genes standing out.
HEs are enzymes found in all forms of microbial life, like phages or even mitochondria and chloroplasts of eukaryotes (Stoddard, 2014). HEs are known as transposable elements that prefer to duplicate into specific genomic regions, which can be either intron, intron-less, or intein sites by a phenomenon called homing (Dujon, 1989). Seg-like and mob-like phage HEs are known as free standing genes and they have been reported as responsible for horizontal gene transfer among host and phage neighboring genes (Zeng et al., 2009).
Additionally the presence of NAD+ biosynthesis and nucleotide metabolism related genes in bacteriophages have been mentioned in various works (Mesyanzhinov et al., 2002; Miller et al., 2003a; Petrov et al., 2010; Javed et al., 2014). By being obligate parasites, bacteriophages are known to carry genes toward their own benefit (Gazzaniga et al., 2009) for host metabolic manipulation. Sirtuins are able to deacetylase acetyl-lysine in proteins (Burgos et al., 2013) and their presence in large bacteriophage genomes is still enigmatic. It is natural for phages to carry molecular tools for metabolic reprogramming of their host, since they need them in order to replicate their genome, before capsid packaging.
Biological characteristics such as morphology, burst size, latency period, adsorption time, and in vitro and in vivo lytic efficacy have been previously reported for two isolated and partially characterized Vibrio phages. These Vibrio alginolyticus lytic phages have also been proposed as promising agents for disinfecting live feeds in aquaculture (Kalatzis et al., 2016), but information over genomic features was lacking. In this study, the phages were sequenced, annotated, and compared, revealing the largest known double stranded DNA Vibrio phage until now with various notable features. Combined study of presence and position of HEs provided insights, such as evolutionary relationships. We also attempted to describe metabolic interactions between the viral and the host biochemical processes by studying the relationship of presumable NAD+ increased biosynthesis, Sir2 viral gene and the necessity of increased ATP accumulation for quality phage DNA replication.
Expanding the knowledge of genomic features and better understanding of the complex phage–host biochemical interactions can provide valuable insights for the efficient application of phage therapy.
Materials and Methods
Phages and Bacterial Host
Both phages φGrn1 and φSt2 belong to Myoviridae family and have been isolated from coastal seawater in Crete, Greece (Kalatzis et al., 2016). The bacterial host was a clinical V. alginolyticus strain isolated from sick gilthead seabream (Sparus aurata) and has been fully sequenced (Castillo et al., 2015).
Amplification, Precipitation, and DNA Extraction of Bacteriophages
Two liquid bacterial cultures of V. alginolyticus strain V1 in the exponential phase of growth were infected separately by bacteriophages φSt2 and φGrn1. The infection was performed with a multiplicity of infection (MOI) of 10 and both tubes were incubated overnight at 25°C with reciprocal shaking. The following day, the cultures were centrifuged and their supernatants were filtered (0.22 μm), tittered and stored at 4°C. Having an optimal titer of 1010 PFU ml-1, phages were concentrated using a standard poly-ethylene glycol/NaCl precipitation (Supplemental Data Sheet S1). DNA extraction was conducted using a Qiagen protocol of the QIAamp DNA Blood Mini Kit (QIAGEN, Hilden, Germany) with the addition of ethanol 100% before the first column wash. A yield of at least 10 μg of DNA was retrieved. Finally, polymerase chain reaction (PCR) (VerityTM Thermal Cycler, Thermo Fisher Scientific, Waltham, MA, USA) with 16S universal primers (forward: 5′-AGAGTTTGATCCTGGCTCAG-3′, reverse: 5′-GACGGGCGGTGTGTACAAG-3′) was conducted before and after DNase RQ1 (Promega, Madison, WI, USA) treatment in order to verify the absence of the host’s or other contaminant DNA (Turner et al., 1999; Glockner et al., 2000). DNA quality was evaluated with NanoDrop (Thermo Fisher Scientific, Waltham, MA, USA) measurements and agarose gel before library construction.
DNA Sequencings and Annotations
Five micrograms of the extracted DNA were used for the construction of a pair-end library with an insert size of 800 bp following Illumina sequencing using an Illumina Hi Seq 2000 (Illumina, San Diego, CA, USA) sequencer. Sequencing was conducted at the Beijing Genomic Institute (Shenzhen, Guangdong, China) according to the manufacturer’s protocol. Possible contaminated reads, primers, N-terminus, and 3′-, 5′-low quality reads were trimmed off with an error rate threshold of 0.05. De novo assembly was conducted using Velvet software (Zerbino and Birney, 2008) under the Geneious platform (R8 version; Biomatters Ltd, Auckland, New Zealand). Finally, assembling resulted in single contigs in both occasions. Annotations were made using ab initio gene predictor Glimmer 3 (Delcher et al., 1999) and Rapid Annotation Subsystem Technology (R.A.S.T.; Aziz et al., 2008; Overbeek et al., 2014) where tRNAs were also identified. Hypothetical proteins were identified by using the B2Go (BioBam, Valencia, Spain) platform against non-reductant protein database and UniProt database with an E-value threshold of ≤10-6 which allowed to identify and manually annotate more coding DNA sequences (CDSs) for φGrn1 and φSt2. Verification of tRNAs took place with tRNAscan-se software1 (version 1.21; Lowe and Eddy, 1997). Synteny was studied by using MAUVE software (Darling et al., 2004) and a list of high homologous genes (blastp threshold: 90) was generated with the online software CoreGene (Zafar et al., 2002). Kyoto Encyclopedia of Genes and Genomes2 (K.E.G.G.; Ogata et al., 1999; Kanehisa et al., 2016) was used for verification of protein products involved in the metabolic processes described in the study.
Whole Genome Alignment and Phylogeny
Whole genome alignment was carried out using the LastZ algorithm (Harris and Pierpoint, 2012). Similarity between the phages was recorded as the highest identity of distances between the different alignments. Whole genome neighbor-joining consensus tree with free end gaps and Tamura–Nei method (bootstrap: 10, consensus method threshold: 87%) was generated, after the alignment, with Geneious (Biomatters Ltd, New Zealand) software. A maximum-likelihood phylogenetic tree for HEs was generated with the MEGA 6 software (Tamura et al., 2013), Jones–Taylor–Thornton substitution model and nearest-neighbor-interchange model of tree interference (100 bootstrap). DELTA-BLAST algorithm was used to compare, identify, and characterize HEs amino acid sequencings and domains. Visualization of alignments was performed with the BioEdit software (version 7.2.5) (Hall, 1999). Relative expression patterns were generated with SigmaPlot (Systat Software Inc., San Jose, CA, USA). Heat maps were created with the FiRe 2.2 Microsoft Excel add-on (University of Fribourg, Fribourg, Switzerland).
Protein Structure and Modeling
The “beta-lactamase domain” open reading frame (ORF) was studied with InterProScan (Jones et al., 2014) and Prosite (de Castro et al., 2006) bioinformatics tools to explore putative domains and active sites. The viral Sir2/cobB protein was studied and visualized using the SwissPdb software viewer3 (version 4.1; Guex and Peitsch, 1997). Packing, solvent exposure, and stereochemical structure were evaluated with Verify 3D4 (Molecular Biology Institute, UCLA, Los Angeles, CA, USA; Bowie et al., 1991; Luthy et al., 1992) and Prosa II (Sippl, 1993; Wiederstein and Sippl, 2007). Investigation of zinc ligand in a finger-like binding site was conducted with ZincExplorer (Chen et al., 2013). Salt bridges were evaluated with ESBRI online software (Costantini et al., 2008).
Transcriptional Study of Bacterial and Viral Genes
Gene expression (Supplemental Data Sheet 1) was studied in wild type (uninfected control) and phage-treated bacteria (MOI:100). Bacteriophages were incorporated into bacterial cultures during exponential phase. Three biological replicates were used for both treatments. After 1 min of vigorous shaking at 25°C, 5 ml of each phage-treated culture were harvested. Harvest was repeated at 5, 10, 20, and 30 min post-infection (p.i.). For the control treatments only one harvest per replicate was performed at 30 min. Cells were immediately centrifuged at 4°C and washed with 150 mM NaCl prior to RNA extraction. The duration of the experiment was set based on the latency time which is 30 min for both phages (Kalatzis et al., 2016).
RNA extraction was performed using a standard TRIzolTM reagent (Thermo Fisher Scientific, Waltham, MA, USA) according to manufacturer protocol. This method resulted in at least 13 μg of RNA per sample. 6 μg of RNA per sample were treated with DNase RQ1 (Promega, Madison, WI, USA) according to manufacturer’s protocol. Samples were tested with PCR to verify purity from bacterial and viral DNA. RNA was then extracted by using a phenol:chloroform protocol. A 70% yield of RNA was retrieved after DNase treatment. Approximately 1 μg of RNA was used per cDNA synthesis by using Superscript II (Thermo Fisher Scientific, Waltham, MA, USA) enzyme according to manufacturer’s protocol. Both bacterial and viral primers for cDNA amplification were designed using Geneious software (Supplementary Table S1) and were tested against both genomic DNAs to confirm that a single amplicon of 70 bp would result from quantitative real-time PCR (qPCR). qPCR was performed on a StepOnePlusTM Real-Time PCR System (Applied Biosystems, Foster City, CA, USA) using SYBR Select Master Mix (Applied Biosystems, Austin, TX, USA), gene-specific primers at a final concentration of 0.2 μM each, and 1 μl of the cDNA as template. Primer specificity and formation of primer dimers were monitored by dissociation curve analysis. The expression levels of V. alginolyticus gyrase A (gyrA) and the HSP70 protein (dnaK) were used as housekeeping (HK) genes to normalize cDNA templates. In order to evaluate the experiment’s and the HK gene’s reliability, two viral gene expression motifs were studied at first after normalization. The first one is a glutaredoxin gene (grx) which is considered as a late early transcribed phage gene and the second one the major capsid protein (MCP), which is considered as a late transcribed one (Luke et al., 2002; Supplementary Figure S1).
The GenBank accession numbers for these two new genomes correspond to KT919972 for φGrn1 and KT919973 for φSt2.
Results and Discussion
Genome Sequencing and Annotation
The genomic sequences of the two bacteriophages were determined after the assembly of short raw reads. All nucleotides were aligned in one single contig of 248,605 bp and a GC content of 38.8% for φGrn1 and 250,485 bp and a GC content of 42.6% for φSt2, which correspond to larger currently known Vibrio bacteriophage genomes according to GenBank. GC content of bacteriophages has been strongly linked to the host’s GC content in Staphylococcus aureus phages (Kwan et al., 2005), although disagreements have been mentioned in other species (Wittmann et al., 2014). The presented phages’ GC content is placed between other Vibrio phages like φpp2 (42.55%) and KVP40 (42.6%) and Enterobacteria T4 phage (35.3%). Host’s GC content (V. alginolyticus V1 strain) is 44.5%. In order to identify the origin of replication and the terminus point of the phage genomes, cumulative GC skews were generated as described before (Grigoriev, 1998, 1999; Uchiyama et al., 2008; Jin et al., 2014). φGrn1 appears to have a putative origin of replication at 51,089 and a putative terminus location at 24,676, whereas φSt2 at nucleotide 1 and 242,751, respectively (Figure 1). Both phage GC skews agree with the transcriptional direction of most of the CDSs. Purine excess is strongly correlated with the leading transcriptional strand (Freeman, 1998).
FIGURE 1. Cumulative GC skews of phages φGrn1 and φSt2. The global minimum and maximum are displayed in the cumulative graph. Putative origin of replication and putative terminus location are highlighted.
In total, 410 and 412 genes were annotated for φGrn1 and φSt2 phages, respectively (Supplementary Tables S2 and S3; Supplementary Figures S2 and S3). By using the online software CoreGene, we were able to detect 77 proteins with high identity among the T4 and large Vibrio phages contained in the GenBank; when we included only Vibrio “schizoT4like” phages in Coregene, identical proteins increased to 271 out of 381 of KVP40’s genome. This is strong indication that the newly characterized phages are “schizoT4like” and can been characterized as KVP40-like, and that VH7D should be included as well in that clade, increasing the number of known and characterized “schizoT4like” Vibrio phages from 3 to 6. When evaluating candidates for phage therapy, it is important to study thoroughly their genome for potential presence of known genes involved in bacterial resistance to antibiotics (Balcazar, 2014), especially when they are physically associated with transposable elements, like HEs. Comparative genome analysis revealed the presence of a small ORF in both studied genomes initially annotated as “beta-lactamase domain protein” (ALP47273 for φGrn1 and ALP47653 for φSt2; Supplementary Tables S2 and S3). Interestingly a similar ORF is present in the KVP40 genome (NP_899337.1) and all Vibrio large phages cited in this work. Protein analysis and amino acid sequence comparison with characterized bacterial beta-lactamases, revealed that these polypeptides exhibited a low degree of similarity, while at the same time residue domains (data not shown) important for catalysis are absent. More specifically metallo-beta-lactamase domain, essential for catalysis (Moali et al., 2003), was not detected using InterProScan, while also no active site was detected using Prosite analysis. Finally, similarity with a well-known and characterized beta-lactamase of the Gram-negative bacteria Stenotrophomonas maltophilia (EC: 184.108.40.206) was lower than 1.4%. Thus, the results of these analyses do not support the automatic in silico annotation as this ORF does not appear to code for a functional beta lactamase protein. Further work is needed for the functional characterization of these ORFs, as the presence of active beta-lactamases on phages could represent a drawback for their application in phage therapy.
Multiple tRNA Genes Are Present in Both Genomes
We verified the structure and the presence of tRNAs for φGrn1 and φSt2, respectively, arranged in clusters in a region of approximately 10,000 bp (37,088 to 46,099) for φGrn1 and approximately 8,000 bp (96,267 to 106,279) for φSt2. Both phages contained two pseudo-forms for GCA and TGC anticodons. Ten and seven hypothetical proteins are scattered inside the tRNA clusters of φGrn1 and φSt2, respectively. tRNAs have only been found in double stranded DNA phages. The perception supported by T4 sequence that the average number of phage’s tRNAs was approximately 10 was challenged after the sequencing of many large phages, showing in addition that the presence of high number of tRNAs is a characteristic of virulent phages (Wilson, 1973); φGrn1 a member of large phages presented here contains 28 (without pseudoforms), one of the highest number of tRNAs identified until today. Deletion of eight T4’s tRNAs has led to a decline in burst size and in phage protein synthesis, highlighting their significance and the evolutionary pressure that favors their conservation through natural selection (Freeman, 1998). Although phages harbor their own tRNAs, they are strongly host-dependent for the efficient translation of their proteins (Kunisawa, 2002). In our case, V. alginolyticus strain V1 possesses a number of at least 67 tRNAs. Codon usage of φGrn1 is able to utilize its tRNAs and incorporate at least 41.8% of the amino acids in its proteins, while this percentage is 39.6% for φSt2. Codon usage from phage’s tRNAs is also associated with low and late expressed viral genes for which host’s tRNAs are not frequent (Kunisawa, 1992).
Whole Genome Phylogenetic and Synteny Study
Both bacteriophages have the highest identity with two large Vibrio phages: ValKK3 (Lal et al., 2016; unclassified, possible “schizoT4like” Vibrio phage as well) and VH7D (Luo et al., 2015). Following whole genome alignment, identity was reported as the highest percentage between more similar alignments. The two phages had a 99.3% nucleotide identity. Highest genomic similarity was identified with VH7D (KC131129) reaching 98.74 and 99.24% for φGrn1 and φSt2, respectively, with the ValKK3 (K671755) at 98.53%. Similarity with large phages KVP40 (AY283928), φpp2 (JN849462), nt-1 (HQ317393), and reference phage T4 (AF138101) was also examined. φGrn1 had a 94.13% similarity along KVP40 and a 94.5% with φpp2, while φSt2 had 93.31 and 93.57%, respectively (Figure 2). Bacteriophage T4 having only a 168,903 bp genomic size, presents 55.99% similarity with φGrn1 and 61.94% with φSt2.
FIGURE 2. Whole genome consensus tree. Neighbor-joining consensus tree of all known Vibrio “schizoT4like” phages. Numbers next to branches represent consensus support. For improved visualization only closer whole genome alignments to the reference are shown. KVP40 bacteriophage was used as reference genome.
Synteny of the genome organization between these phages was also examined as it has been proposed in Acinetobacter T4-like phages (Jin et al., 2014). The alignment of φGrn1, φSt2, KVP40, φpp2, nt-1, and VH7D resulted (Figure 3) into two small synteny local collinear blocks (LCBs) with 3,327 bp (purple) and 13,585 bp (blue) and three large LCBs with 57,120 bp (red), 66,636 bp (green), and 106,660 bp (yellow), indicating DNA regions which are homologous among the genomes. Graphs inside the blocks show very high similarity between the genomes, however, there are some non-identical genomic regions that are represented with white color inside the blocks. Although there seems to be a genomic rearrangement, the block sequence remains the same across the genomes of all phages, which can be speculated as a result of the possible circularly permuted linear double stranded DNA genomes of T4-like phages and their conserved genome organization (Petrov et al., 2006, 2010).
FIGURE 3. Multiple genome alignment of bacteriophage genomes. Genomes of all known Vibrio “schizoT4like” bacteriophages were compared using Mauve software. Local collinear blocks (LCBs) are highlighted with different colors. Same colored blocks indicate high synteny between genomes without genomic rearrangements. Graphs inside the blocks represent the level of synteny. White regions represent unique genomic regions. For improved visualization lines connecting the high syntenic regions have been obliterated.
Identification of Viral Lysozymes and Endolysins
Genome analysis also led to two ORFs for each phage belonging to the large lysozyme superfamily. The first one is a tail lysozyme at position 57,157 for φGrn1 (ALP46994) and 116,942 for φSt2 (ALP47375). Both tail lysozymes contain the N-terminal of the OB domain, which was firstly thoroughly described for the T4 bacteriophage (gp5) (Nakagawa et al., 1985; Kanamaru et al., 1999; Arisaka et al., 2003). Additionally, one transglycosylase for each phage was detected at position 8,121 for φGrn1 (ALP47078) and 69,098 for φSt2 (ALP47458). Both transglycosylases harbor the SLT domain characteristic of murein hydrolases (Van Asselt et al., 1999). Application of recombinant endolysins can be a major biotechnological weapon against infectious bacteria. Generally, difficulties toward successful endolysin utilization have been reported in Gram-negative bacteria, mainly due to their limited access to the interior peptidoglycan (Hermoso et al., 2007; Lai et al., 2011). However, application of endolysins in a slightly acidic environment has tackled this obstacle (Oliveira et al., 2014) and reports of successful extracellular applications of recombinant endolysins against Gram-negative bacteria have been recently reported (Lim et al., 2014; Oliveira et al., 2016).
Unique HEs Are Present in Both Viral Genomes
Both phages presented here bear HEs. With the increasing number of T4-like phages being sequenced and annotated, it becomes clear that the 15 HEs of the T4 phage is a feature unique to that phage (Edgell et al., 2010). Represented in our work, φGrn1 has only one HE (like KVP40), while φSt2 has three (like φpp2). We were able to identify conserved domains, compare and study syntenic relationships with other phage-associated HEs. The HE of φGrn1 (ALP47050) is located at 90,429 and encodes for a 238 aa protein. It belongs to the seg-like HE family as it encodes the GIY_YIG domain of the homonymous superfamily. DELTA-BLAST showed a 43% similarity with Salmonella phage S16 HE (YP_007501215) with a 51% query coverage (E-value = 9 × e-23). Concerning Vibrio phages, ValKK3 presents the highest similarity with the seg-like HE (AJT61075) with a query coverage of 52% (E-value = 2 × e-22). BlastX algorithm could not match any nucleotide sequences with this particular HE. Synteny study showed that it is established at the same genomic area as the T4 phage segD HE (NP_049788.2), but with opposite orientation (Figure 4A). These data support the hypothesis that this is a unique HE, although the genomic region after the MCP is a frequent site of reported HEs in Enterobacteria phages (like T4). For φSt2, three HEs were identified and annotated. The first HE is a mob-like HE containing the HNHc domain (ALP47430) and is located at 240,309. It has a 98% similarity with the unique HE (Lin and Lin, 2012) of the φpp2 bacteriophage (AFN37352) and a query coverage of 100% (E-value = 6 × e-150). It is also established at the same genomic area like φpp2, revealing a possible evolutionary relationship between them (Figure 4B; Supplementary Figure S4). The other two HEs are characterized as seg-like endonucleases, because they both contain the GIY_YIG domain. Specifically, the second HE (ALP47453) at 216,143 appears to have only a 46% similarity with the HE of ValKK3 bacteriophage (AJT61075) with a query coverage of 97% (E-value = 3 × e-43) and subsequently should be enlisted as unique as well. Bacteriophage φpp2 also has a HE established just after the regA gene, while it is located upstream in φSt2 (Figure 4C). Finally, the third HE (ALP47432) at 6,477 contains three main domains. Apart from the GIY_YIG domain, it also has two tandem repeats of the NUMOD3 domain (Sitbon and Pietrokovski, 2003) and presents the highest similarity (90% similarity with a 100% query coverage, E-value = 4 × e-80) with the segD HE of the KVP40 Vibrio phage (NP_899393), although it was previously described as a unique one. It is established at the same area as KVP40’s after the rIIA and rIIB early lysis protectors (Figure 4D; Supplementary Figure S4) suggesting an evolutionary relationship between the bacteriophages. In order to avoid deletion, transposition of HEs can take place within a phage genome without compromising the viability of the phage and can be responsible for genomic shuﬄing (Sandegren et al., 2005) because they can also cause mobility of the surrounding components (Belfort, 1991). Lin and Lin (2012) described the occurrence of unique HEs and their indels in bacteriophage φpp2 as evidence of being a distinct new species different from KVP40 within the T4-like phage family. Although they act as selfish DNA features (Jin et al., 2014), reports and characterizations of viral HEs can highlight putative evolutionary processes. For instance, φSt2 although geographically distant, may be an evolutionary link between φpp2 and KVP40. This can be also supported by the fact that φSt2 can also infect V. parahaemolyticus (strain V2; Kalatzis et al., 2016), which is also the host-species of φpp2 and KVP40. Apart from the evolutionary impact that HEs have in phagic genomes and the population diversity that can confer (Lin and Lin, 2012), they are also responsible for regulating surrounding genes, highlighting the importance of locating and characterizing them in comparative genomic studies (Edgell et al., 2010; Stoddard, 2014).
FIGURE 4. Neighbor genes of HEs from various T4-like bacteriophages compared to φGrn1 and φSt2. Same colored arrows represent homologous genes. Light yellow arrows represent putative HE of compared bacteriophages. The genomic areas were aligned according to the gene represented from the red arrow. (A) Neighbor gene products of unique φGrn1 segD HE. 1 for prohead assembly protein, 2 for major capsid protein, and 3 for inhibitor of prohead lysis. (B) Neighbor genes of φSt2 HE, mob-like. 2 for GTP cyclohydrolase I. 4 for hypothetical protein. (C) Neighbor genes of φSt2 unique HE, seg-like. 2 and 3 DNA polymerase clamp loader subunit. 4 for regA. (D) Neighbor genes of φSt2 segD HE. 1 for rIIA protector, 2 for rIIB protector, and 3 for hypothetical protein.
Viral Infection Results in Complex Metabolic Interactions with the Host
The existence of NAD+ salvage enzymes in the genome of bacteriophage KVP40 has been previously described and thought to be of interest (Miller et al., 2003a). All of the Vibrio phages cited in the present publication bear a number of enzymes associated with the conversion of nicotinamide to NAD+, so it seems that large phages have the molecular tools to increase the cellular NAD+ content. Additionally, several bacteriophages carry genes with products involved in pyrimidine and purine biosynthesis. It is believed that large genome size imposes a strong need for DNA replication. Thus, the phage can increase the bacterial capacity for nucleotide biosynthesis in order to enhance the pathways and therefore its gene dosage (Karam and Drake, 1994; De Smet et al., 2016). This is also the case for the bacteriophages presented in this work. In the interest of pointing out the significance of the described genes and their effect during the infection, relative expression levels of both bacterial and bacteriophagic genes were studied (Figure 5; Supplementary Data Sheet S2; Supplementary Table S1), where the viral relative transcript levels outnumbered the host’s, as expected (Chevallereau et al., 2016). Apart from the nicotinamide mononucleotide adenyltransferase (NMNAT; ALP47012 for φGrn1 and ALP47393 for φSt2), an enzyme present in both phages and the host, φGrn1 and φSt2 also possess a nicotinamide phosphoribosyltransferase (NAMPT; ALP46980 for φGrn1 and ALP47363 for φSt2), which is absent from V. alginolyticus and is able to utilize nicotinamide as a substrate and convert it to nicotinamide d-ribonucleotide, making a shortcut in the NAD+ biosynthesis pathway. No fluxes in the expression levels of the bacterial genes of pyrazinamidase (pnac) and NMNAT were noticed, showing that the pathway remains unaffected. On the other hand, φSt2’s NAMPT gene is transcribed almost immediately after the infection and NMNAT shows that the bacteriophage is probably utilizing intracellular nicotinamide instantly for NAD+ biosynthesis by using a quick two-step pathway. This suggests that bacteriophages try to enhance NAD+ production. Although many NAD+-dependent enzymes important for DNA replication have been reported (DNA ligase) in T4-like bacteriophages (Hertveldt et al., 2005), the characterized enzymes in the bacteriophages presented in this work are ATP-dependent or NADPH-dependent, except for a sirtuin, a deacetylase protein (Sir2/CobB protein, Sir2) (ALP47040 for φGrn1 and ALP47418 for φSt2). This sirtuin is a NAD+-dependent conserved enzyme among large Vibrio phages and it was initially reported in bacteriophages after the sequencing of the KVP40 phage, where it was described as having a NAD+ hydrolysis role at the time. The homologous eukaryotic Sirt2 and Sirt3 proteins have been connected to increased lifespan and cell growth (Frye, 2000; Chang and Min, 2002). Both eukaryotic and prokaryotic sirtuins are known for post-translational modifications using NAD+ as a co-substrate. Specifically, deacetylation of acetyl-lysine by the sirtuins produces nicotinamide as a byproduct, which can also be recycled for NAD+ biosynthesis as mentioned above (Burgos et al., 2013). Deacetylation of acetyl-lysine has been strongly correlated with the activation of acetyl-coA synthetase (ACS) in prokaryotes, which is characterized as a Sir2-dependent enzyme (Starai, 2002). Lysine acetylation is a major post-translational modification in both prokaryotic and eukaryotic proteins and is a frequent regulatory phenomenon (Ouidir et al., 2015). V. alginolyticus strain V1 already possesses a Sir2/cobB deacetylase protein (Sir2). The presence of an additional highly divergent Sir2/cob protein in bacteriophages may suggest bacterial protein activation (such as the activation of ACS) by post-translational modifications by the virus. The strong relationship between NAMPT and sirtuins has been well described in prokaryotes and eukaryotes, along with the recycling of nicotinamide and the importance of the de novo synthesis of NAD+ (Imai et al., 2000; Lin et al., 2010; Burgos et al., 2013). The relative expression levels of the Sir2 gene are high 20 min p.i., implying a possible priority of NAMPT and NMNAT regulation, and therefore NAD+ production before Sir2 transcription. Sir2/cobB protein is able to activate ACS protein by deacetylation of acetyl-lysine and by consuming 1 ATP molecule, and therefore possibly advance toward the synthesis of acetyl-coenzyme A (AcoA; Gulick et al., 2003). Interestingly, a statistically significant increase in transcription levels is observed in the corresponding bacterial Sir2 during late infection and in the ACS gene, which may be part of the phage’s metabolic manipulation and the biotic stress the cell is experiencing. AcoA can then be incorporated in the bacterial citric cycle and be one of the sources of increased intracellular ATP content by consuming the abundant NAD+. The presence of multiple ATP-dependent enzymes could indicate a high ATP demand during the lytic cycle. Direct evidence of increased accumulation of ATP during phage infection has been recently provided by Chevallereau et al. (2016; Supplementary Table S3). Additionally in our case, at least 97 genomes of 250,485 bp have to be synthetized during phage infection, which could dictate even higher nucleotide metabolic demand and subsequently energy in the form of ATP. Specifically, in association to purine metabolism, both phages carry the two subunits of a ribonucleoside diphosphate reductase (nrdAB; ALP46965, ALP46998 for φGrn1 and ALP47345, ALP47379 for φSt2) and a ribonucleoside triphosphate reductase (nrdD; ALP46970 for φGrn1 and ALP47350 for φSt2), both involved in the final steps of dATP and dGTP biosynthesis. In addition to nrdAB and nrdD, more enzymes are involved in pyrimidine metabolism; a dCMP deaminase (ALP47131 for φGrn1 and ALP47515 for φSt2), a thymidine kinase (ALP47080 for φGrn1 and ALP47460 for φSt2), and a thymidylate synthase (thyA; ALP47026 for φGrn1 and ALP47405 for φSt2) are also present, enhancing dTTP biosynthesis from dCTP, and dUTP, a well-established fact during phage infection. Host and viral RNA decay could potentially be a source of free nucleoside diphosphates during T4-like phage infection (Carpousis et al., 1989, 1994; Ueno and Yonesaki, 2004; Uzan, 2009) and along with the presence of viral nrdAB and nrdD ATP-dependent ribonucleotidases, the phage might be able to enhance the much needed deoxyribonucleotide biosynthesis for DNA replication (Chevallereau et al., 2016). These enzymes reach their transcript levels plateau in 20 min p.i. for φSt2, while T4 bacteriophage—having a 20-min latency time—reaches these levels at 10 min p.i. (Luke et al., 2002). Statistically significant differences are noted in the two nrdD host genes, a phenomenon also observed recently during viral infection in a Pseudomonas aeruginosa strain (Chevallereau et al., 2016), an obligatory anaerobic enzyme. Bacterial ribonucleotide reductases are known for allosteric and transcriptional regulation depending on the balance of NTPs present in the cell (Torrents, 2014). Increased mutation rates during DNA replication can take place if uneven presence of NTPs is spotted (Wheeler et al., 2005). This upregulation of the two host nrdD reductases hint toward an imbalance of ATP content in the cell, which can be justified as described in this section. Interestingly, both phages contain a dUTP pyrophosphatase (DUT; AL47106 for φGrn1 and AP47489 for φSt2), which has been reported in all Vibrio “schizoT4like” viruses and most bacterial species, but has not been reported or annotated in any V. alginolyticus bacterial strains (including V1). This host possibly lacks a dedicated enzyme for diphosphatase activity in order to one-step hydrolyze dUTP to dUMP, or that a possible chimeric protein has this enzymatic activity (Moroz et al., 2005). Nonetheless, V. alginolyticus can initially convert dUTP to dUDP with a nucleoside-diphosphate kinase (NDK) and then to dUMP with a dTMP kinase (TMK) (Kielley, 1970; Chakrabarty, 1998). DUTP pyrophosphatase has been extensively studied in yeast and proven to be efficient in preventing the incorporation of uracil into DNA during the replication stage (Gadsden et al., 1993). This suggests that whereas dUTP pyrophosphatase is absent in V. alginolyticus, lytic Vibrio bacteriophages carry it to possibly satisfy the need for quick uracil hydrolysis, which can interfere during the rolling circle replication if misused by DNA polymerase as a building block. This hypothesis is also supported by the presence of dUTPase in many retroviruses and the enzyme’s role in circumventing the deleterious effects of high uracil presence during the reverse transcription of the viral RNA (Hizi and Herzig, 2015). It is noteworthy that dUTP pyrophosphatase is found in most sequenced bacterial genomes, with Escherichia coli knock-out mutants resulting in accretion of putative short Okazaki fragments and subsequent errors in DNA replication (Tye and Lehman, 1977; Shlomai and Kornberg, 1978), while the T4 phage carries a bifunctional homologous dCTPase-dUTPase gene (gp56), which also takes part in forming 5-hydroxy-methyl-cytosine (Gary et al., 1998). Viral DNA replication peaks at 20 min p.i. (last one-third of the latency period), at the same time point that we were able to identify maximum expression levels for the DUT gene, showing a threefold increase in comparison to 10 min p.i. This indicates the possible need to shift the nucleotide biosynthesis balance toward DNA replication, rather than RNA production (transcription), in order to prevent DNA polymerase from using uracil as a substrate for DNA synthesis. At 10 min p.i. we also noticed fluxes in the transcription levels of the NDK (statistically significant) and TMK bacterial genes, possibly mirroring the high uracil content and the need to hydrolyze it. The importance of hydrolyzing dUTP and converting it to dTMP, and later dTTP, is also reflected by the presence of the thyA gene in the phage and its induced transcription levels, along with DUT. Fluxes of the two bacterial thyAs are also observed, but not statistically significant.
FIGURE 5. Schematic representation of biochemical processes during the infection of V. alginolyticus V1 from bacteriophage φSt2. Heat maps represent gradient changes in relative transcript levels of bacterial genes (black bolded) for control (C), 1 min post-infection (p.i.) (1′), 5 min p.i. (5′), 10 min p.i. (10′), 20 min p.i. (20′), and 30 min p.i. (30′) treatments. White asterisks represent statistically significant differences compared to control treatments (p < 0.05). Bars represent relative transcript levels of bacteriophage genes (purple bolded; ±SE) for 1 min p.i., 5 min p.i., 10 min p.i., 20 min p.i., 30 min p.i. treatments. Arrows represent bacterial processes and dotted arrows represent possible phage processes. Dark red arrows represent extracellular compounds that are inserted in the cell. Green cycle represents possible nicotinamide recycling.
In silico Functionality Study of Sir2/cobB Protein
In an attempt to provide insights in the relationship between NAD+ production and nucleotide biosynthesis, we noted the possible neuralgic role Sir2/cobB protein may have and tried to partially characterize it. Post-translational protein modifications are of high research interest, especially if they are taking place during a host–parasite interaction. Phages bear a large number of non-functional ORFs and studies aiming for functional verification can prove valuable. Both phages possess the same Sir2/cobB protein. A molecular model was constructed based on the Sir2/cobB crystal structure of the homologous protein from E. coli (PDB ID: 1SP5P; Zhao et al., 2004). Modeling predicted 10 protein sheets and 9 helices, along with a large Rossmann fold domain, a small Zinc binding domain, and the loops connecting the two (Supplementary Figure S5). Verify3D and PROSA II profiles (Z score -5.32) of packing, solvent exposure, and stereochemical structure proved that the final model was of high overall quality. Despite the low protein identity between E. coli and bacteriophage φSt2 (23.5 and 24.1% with V. alginolyticus), the overall predicted secondary structure is similar (Figure 6A). Unlike the bacterial protein, superposition of the phage’s protein shows absence of ligand zinc in the small finger domain, a phenomenon also observed in the Sir2/cobB protein of KVP40 bacteriophage. Furthermore, verification of the absence of a zinc ligand site was carried out by ZincExplorer software. This raises the question of functionality of the enzyme due to a possibly unstable secondary structure. Thus we examined the possibility of a salt bridge forming instead in the finger-like domain (Kumar and Nussinov, 2002). Although zinc binding sites can be extremely important for a stable secondary structure (Webster et al., 1991), salt bridges are able to substitute them efficiently in prokaryotic proteins (Baglivo et al., 2009). Prediction of salt bridges resulted in the detection of two possible sites between the aspartic acid and arginine residues in the 112 and 135 positions, respectively (Figure 6B) showing that the finger domain could remain stable. An additional salt bridge might also take place between lysine 116 and glutamic acid 129, with its presence needing further examination. Active site of the viral Sir2/cobB protein is detected at the histidine 112 position, a conserved region of bacterial and viral sirtuins. Additionally, the binding site of acetyl-lysine is also conserved, with domains FNE and INP creating a similar to E. coli tunnel for acetyl-lysine to bind close to the active site (Figure 7). Specifically, the strictly conserved phenylalanine 190 and proline 221 (of E. coli Sir2/cobB protein) are present, in contrary the bacterial conserved tyrosine 220 is replaced by the polarly neutral amino acid Asparagine. This is a somewhat usual replacement also present in at least Archaea Archaeoglobus fulgidus’ Sir2Af1 (Ringel et al., 2014) and the human SIRT2 protein (Feldman et al., 2015), with the polarly neutral glutamine replacing tyrosine and remaining functional. Finally, valine 219 of E. coli is replaced with isoleucine, both being non-polar aliphatic amino acids. This information suggests that viral Sir2/cobB proteins of Vibrio “schizoT4like” phages can act similarly to the bacterial and are able to deacetylase acetyl-lysines of enzymes, like ACS, and subsequently activate them. The Sir2/cobB protein is also conserved at the genomes of E. coli bacteriophage T5 (Wang et al., 2005), at Salmonella phage SPC35 (Kim and Ryu, 2011), at Cronobacter phage vB_CsaM_GAP32 (Abbasifar et al., 2014), at Pectobacter phage My1 (Lee et al., 2012) and at Klebsiella phage JD001 (Cui et al., 2012).
FIGURE 6. Schematic representation of φSt2 and φGrn1 Sir2/cobB protein. (A) Similar (1) and dissimilar (2) amino acid residues of Sir2/cobB protein of φSt2 (magenta) are shown as superposition of structural models, with Sir2/cobB protein of V. alginolyticus (turquoise). Protein of E. coli is also shown (white). Ligand zinc of the finger domain is also represented (yellow sphere). (B) Ribbon diagram represents the finger domain of Sir2/cobB protein. Sheets (yellow), helices (red), and coils (green) are highlighted. Carbon (white), nitrogen (black), and oxygen (red) of the putative salt bridge residues are shown (represented as sticks). Distances between atoms are shown with yellow numbers (Armstrong).
FIGURE 7. Alignment of Sir2/cobB proteins of E. coli, V. alginolyticus, φSt2, and KVP40. Red asterisk highlights the histidine active site. Boxes represent conserved regions of active sites (red), acetyl-lysine binding sites (orange), salt bridges of finger-like viral protein domains (blue), zinc ligand sites of finger-like bacterial protein domains (magenta).
The knowledge of the genomic organization of bacteriophages provides a valuable insight into their interactions with the host cell and facilitates their efficient application for phage therapy. Although phage therapy is a relatively old technique, it still lacks in basic research in order to understand processes of auxiliary metabolism during infection. This work underlines the significance of clarifying biochemical processes and interactions during bacteriophage infection and also, host metabolic hijacking, as a result of features unmasked after DNA sequencing. By determining the genome of these two phages and describing genomic features we know six (possibly seven including ValKK3) fully characterized large “schizoT4like” Vibrio bacteriophages with a wide spectrum of bacterial hosts of Vibrio species, which threat fisheries and aquaculture, increasing our “armory” against vibriosis without the use of antibiotics. Both importance of detailed genomic study and characterization of endolysins as potential antibacterial agents are highlighted in the manuscript. Additionally, the characterization and localization of HEs resulted to the description of viral evolutionary relationships. In an attempt to expand our knowledge in phage infection and lysis efficacy, we monitored both bacterial and viral gene regulation related to NAD+ biosynthesis and nucleotide metabolism, during the latency period of the φSt2 phage. The results hint toward possible post-translational modifications by the viral Sir2 gene in order to activate inert bacterial ACS protein, a binary model never previously described in detail. Although there is a high variability in the Sir2/cobB protein, its partial characterization indicates that it can act similarly to the bacterial ones and also contribute to the increased cell needs in ATP for an efficient phage DNA replication. Overall our data raise the possibility that the ability of large phages to maximize biochemical host exploitation could render phage therapy more efficient. Future experiments including the biochemical characterization of the Sir2/cobB protein, its deletion and the monitoring of intracellular ATP during large-size phage infection can strengthen our assumptions. This information could maybe apply in the future toward creating more efficient molecularly engineered virions in the battle against bacterial drug-resistant infections.
EF, PK, and DS conceived the study and designed the research. DS and PGK performed molecular work. DS performed bioinformatics work. DS and PGK analyzed the data. EF and DS wrote the manuscript and discussed the results and all the authors commented on the manuscript.
This work was kindly funded and supported by the Greek National Strategic Reference Framework 2007–2013 of General Secretariat for Research and Technology (co-funded by European Social Fund and Greek National Funds), FISHPHAGE project 131, European Union, FP7 Marie Curie, IRSES 2010, AQUAPHAGE project 269175, and PROAQUA project 12-132390 (Danish Committee for Strategic Research in Health, Food and Welfare). The funders had no role in study, design, analysis, decision to publish, or preparation of the manuscript.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2016.01807/full#supplementary-material
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
- ^ http://lowelab.ucsc.edu/tRNAscan-SE/
- ^ http://www.genome.jp/kegg/
- ^ http://www.expasy.org/spdbv/
- ^ http://services.mbi.ucla.edu/Verify_3D/
Abbasifar, R., Griffiths, M. W., Sabour, P. M., Ackermann, H. W., Vandersteegen, K., Lavigne, R., et al. (2014). Supersize me: Cronobacter sakazakii phage GAP32. Virology 46, 138–146. doi: 10.1016/j.virol.2014.05.003
Arisaka, F., Kanamaru, S., Leiman, P., and Rossmann, M. G. (2003). The tail lysozyme complex of bacteriophage T4. Int. J. Biochem. Cell Biol. 35, 16–21. doi: 10.1016/S1357-2725(02)00098-5
Aziz, R. K., Bartels, D., Best, A. A., DeJongh, M., Disz, T., Edwards, R. A., et al. (2008). The RAST server: rapid annotations using subsystems technology. BMC Genomics 9:75. doi: 10.1186/1471-2164-9-75
Baglivo, I., Russo, L., Esposito, S., Malgieri, G., Renda, M., Salluzzo, A., et al. (2009). The structural role of the zinc ion can be dispensable in prokaryotic zinc-finger domains. Proc. Natl. Acad. Sci. U.S.A. 106, 6933–6938. doi: 10.1073/pnas.0810003106
Balcazar, J. L. (2014). Bacteriophages as vehicles for antibiotic resistance genes in the environment. PLoS Pathog. 10:e1004219. doi: 10.1371/journal.ppat.1004219
Belfort, M. (1991). Self-splicing lntrons in prokaryotes. Cell 64, 9–11. doi: 10.1016/0092-8674(91)90201-9
Bowie, J. U., Luthy, R., and Eisenberg, D. (1991). A method to identify protein sequences that fold into a known three- dimensional structure. Science 253, 164–170. doi: 10.1126/science.1853201
Burgos, E. S., Vetticatt, M. J., and Schramm, V. L. (2013). Recycling nicotinamide. The transition-state structure of human nicotinamide phosphoribosyltransferase. J. Am. Chem. Soc. 135, 3485–3493. doi: 10.1021/ja310180c
Carpousis, A. J., Mudd, E. A., and Krisch, H. M. (1989). Transcription and messenger RNA processing upstream of bacteriophage T4 gene 32. Mol. Gen. Genet. 219, 39–48. doi: 10.1007/BF00261155
Carpousis, A. J., Van Houwe, G., Ehretsmann, C., and Krisch, H. M. (1994). Copurification of E. coli RNAase E and PNPase: evidence for a specific association between two enzymes important in RNA processing and degradation. Cell 76, 889–900. doi: 10.1016/0092-8674(94)90363-8
Castillo, D., D’Alvise, P., Kalatzis, P. G., Kokkari, C., Middelboe, M., Gram, L., et al. (2015). Draft genome sequences of Vibrio alginolyticus strains V1 and V2, opportunistic marine pathogens. Genome Announc. 3:e00729-15. doi: 10.1128/genomeA.00729-15
Chakrabarty, A. M. (1998). Nucleoside diphosphate kinase: role in bacterial growth, virulence, cell signalling and polysaccharide synthesis. Mol. Microbiol. 28, 875–882. doi: 10.1046/j.1365-2958.1998.00846.x
Chang, K. T., and Min, K. T. (2002). Regulation of lifespan by histone deacetylase. Ageing Res. Rev. 1, 313–326. doi: 10.1016/S1568-1637(02)00003-X
Chen, Z., Wang, Y., Zhai, Y.-F., Song, J., and Zhang, Z. (2013). ZincExplorer: an accurate hybrid method to improve the prediction of zinc-binding sites from protein sequences. Mol. Biosyst. 9, 2213–2222. doi: 10.1039/c3mb70100j
Chevallereau, A., Blasdel, B. G., De Smet, J., Monot, M., Zimmermann, M., Kogadeeva, M., et al. (2016). Next-generation “-omics” approaches reveal a massive alteration of host RNA metabolism during bacteriophage infection of Pseudomonas aeruginosa. PLoS Genet. 12:e1006134. doi: 10.1371/journal.pgen.1006134
Costantini, S., Colonna, G., and Facchiano, A. M. (2008). Bioinformation ESBRI: A web server for evaluating salt bridges in proteins. Bioinformation 3, 137–138. doi: 10.6026/97320630003137
Cui, Z., Shen, W., Wang, Z., Zhang, H., Me, R., Wang, Y., et al. (2012). Complete genome sequence of Klebsiella pneumoniae phage JD001. J. Virol. 86:13843.
Darling, A. C. E., Mau, B., Blattner, F. R., and Perna, N. T. (2004). Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 14, 1394–1403. doi: 10.1101/gr.2289704
de Castro, E., Sigrist, C. J. A., Gattiker, A., Bulliard, V., Langendijk-Genevaux, P. S., Gasteiger, E., et al. (2006). ScanProsite: detection of PROSITE signature matches and ProRule-associated functional and structural residues in proteins. Nucleic Acids Res. 4, 362–365. doi: 10.1093/nar/gkl124
De Smet, J., Zimmermann, M., Kogadeeva, M., Ceyssens, P.-J., Vermaelen, W., Blasdel, B., et al. (2016). High coverage metabolomics analysis reveals phage-specific alterations to Pseudomonas aeruginosa physiology during infection. ISME J. 10, 1823–1835. doi: 10.1038/ismej.2016.3
Delcher, A. L., Harmon, D., Kasif, S., White, O., and Salzberg, S. L. (1999). Improved microbial gene identification with GLIMMER. Nucleic Acids Res. 27, 4636–4641. doi: 10.1093/nar/27.23.4636
Dujon, B. (1989). Group I introns as mobile genetic elements: Facts and mechanistic speculations – a review. Gene 82, 91–114. doi: 10.1016/0378-1119(89)90034-6
Edgell, D. R., Gibb, E. A., and Belfort, M. (2010). Mobile DNA elements in T4 and related phages. Virol. J. 7, 290–304. doi: 10.1186/1743-422X-7-290
Feldman, J. L., Dittenhafer-Reed, K. E., Kudo, N., Thelen, J. N., Ito, A., Yoshida, M., et al. (2015). Kinetic and structural basis for Acyl-group selectivity and NAD+ dependence in sirtuin-catalyzed deacylation. Biochemistry 54, 3037–3050. doi: 10.1021/acs.biochem.5b00150
Freeman, J. M. (1998). Patterns of genome organization in bacteria. Science 279:1827.
Frye, R. A. (2000). Phylogenetic classification of prokaryotic and eukaryotic Sir2-like proteins. Biochem. Biophys. Res. Commun. 273, 793–798. doi: 10.1006/bbrc.2000.3000
Gadsden, M. H., McIntosh, E. M., Game, J. C., Wilson, P. J., and Haynes, R. H. (1993). dUTP pyrophosphatase is an essential enzyme in Saccharomyces cerevisiae. EMBO J. 12, 4425–4431.
Gary, T. P., Colowick, N. E., and Mosig, G. (1998). A species barrier between bacteriophages T2 and T4: Exclusion, join- copy and join-cut-copy recombination and mutagenesis in the dCTPase genes. Genetics 148, 1461–1473.
Gazzaniga, F., Stebbins, R., Chang, S. Z., McPeek, M. A., and Brenner, C. (2009). Microbial NAD metabolism: lessons from comparative genomics. Microbiol. Mol. Biol. Rev. 73, 529–541. doi: 10.1128/MMBR.00042-08
Glockner, F. O., Zaichikov, E., Belkova, N., Denissova, L., Pernthaler, J., Pernthaler, A., et al. (2000). Comparative 16S rRNA analysis of lake bacterioplankton reveals globally distributed phylogenetic clusters including an abundant group of actinobacteria. Appl. Environ. Microbiol. 66, 5053–5065. doi: 10.1128/AEM.66.11.5053-5065.2000
Grigoriev, A. (1998). Analyzing genomes with cumulative skew diagrams. Nucleic Acids Res. 26, 2286–2290. doi: 10.1093/nar/26.10.2286
Grigoriev, A. (1999). Strand-specific compositional asymmetries in double-stranded DNA viruses. Virus Res. 60, 1–19. doi: 10.1016/S0168-1702(98)00139-7
Guex, N., and Peitsch, M. C. (1997). SWISS-MODEL and the Swiss-PdbViewer: an environment for comparative protein modeling. Electrophoresis 18, 2714–2723. doi: 10.1002/elps.1150181505
Gulick, A. M., Starai, V. J., Horswill, A. R., Homick, K. M., and Escalante-Semerena, J. C. (2003). The 1.75 Å crystal structure of acetyl-CoA synthetase bound to adenosine-5’-propylphosphate and coenzyme A. Biochemistry 42, 2866–2873. doi: 10.1021/bi0271603
Hall, T. (1999). BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl. Acids Symp. Ser. 41, 91–98.
Harris, F., and Pierpoint, L. (2012). Photodynamic therapy based on 5-aminolevulinic acid and its use as an antimicrobial agent. Med. Res. Rev. 29, 1292–1327. doi: 10.1002/med.20251
Hermoso, J. A., García, J. L., and García, P. (2007). Taking aim on bacterial pathogens: from phage therapy to enzybiotics. Curr. Opin. Microbiol. 10, 461–472. doi: 10.1016/j.mib.2007.08.002
Hertveldt, K., Lavigne, R., Pleteneva, E., Sernova, N., Kurochkina, L., Korchevskii, R., et al. (2005). Genome comparison of Pseudomonas aeruginosa large phages. J. Mol. Biol. 354, 536–545. doi: 10.1016/j.jmb.2005.08.075
Hizi, A., and Herzig, E. (2015). dUTPase: the frequently overlooked enzyme encoded by many retroviruses. Retrovirology 12, 70–84. doi: 10.1186/s12977-015-0198-9
Imai, S., Armstrong, C. M., Kaeberlein, M., and Guarente, L. (2000). Transcriptional silencing and longevity protein Sir2 is an NAD-dependent histone deacetylase. Nature 403, 795–800. doi: 10.1038/35001622
Jassim, S. A. A., and Limoges, R. G. (2014). Natural solution to antibiotic resistance: bacteriophages “The Living Drugs.” World J. Microbiol. Biotechnol. 30, 2153–2170. doi: 10.1007/s11274-014-1655-7
Javed, M. A., Ackermann, H. W., Azeredo, J., Carvalho, C. M., Connerton, I., Evoy, S., et al. (2014). A suggested classification for two groups of Campylobacter myoviruses. Arch. Virol. 159, 181–190. doi: 10.1007/s00705-013-1788-2
Jin, J., Li, Z.-J., Wang, S.-W., Wang, S.-M., Chen, S.-J., Huang, D.-H., et al. (2014). Genome organisation of the Acinetobacter lytic phage ZZ1 and comparison with other T4-like Acinetobacter phages. BMC Genomics 15:793. doi: 10.1186/1471-2164-15-793
Jones, P., Binns, D., Chang, H. Y., Fraser, M., Li, W., McAnulla, C., et al. (2014). InterProScan 5: genome-scale protein function classification. Bioinformatics 30, 1236–1240. doi: 10.1093/bioinformatics/btu031
Kalatzis, P. G., Bastías, R., Kokkari, C., and Katharios, P. (2016). Isolation and characterization of two lytic bacteriophages, φSt2 and φGrn1; phage therapy application for biological control of Vibrio alginolyticus in aquaculture live feeds. PLoS ONE 11:e0151101. doi: 10.1371/journal.pone.0151101
Kanamaru, S., Gassner, N. C., Ye, N., Takeda, S., and Arisaka, F. (1999). The C-terminal fragment of the precursor tail lysozyme of bacteriophage T4 stays as a structural component of the baseplate after cleavage. J. Bacteriol. 181, 2739–2744.
Kanehisa, M., Sato, Y., Kawashima, M., Furumichi, M., and Tanabe, M. (2016). KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 44, D457–D462. doi: 10.1093/nar/gkv1070
Karam, J. D., and Drake, J. W. (1994). Molecular Biology of Bacteriophage. Washington, DC: American Society for Microbiology.
Kielley, R. K. (1970). Purification and properties of thymidine monophosphate kinase from mouse hepatoma. J. Biol. Chem. 245, 4204–4212.
Kim, M., and Ryu, S. (2011). Characterization of a T5-like coliphage, SPC35, and differential development of resistance to SPC35 in Salmonella enterica serovar typhimurium and Escherichia coli. Appl. Environ. Microbiol. 77, 2042–2050. doi: 10.1128/AEM.02504-10
Klumpp, J., Fouts, D. E., and Sozhamannan, S. (2012). Next generation sequencing technologies and the changing landscape of phage genomics. Bacteriophage 2, 190–199. doi: 10.4161/bact.22111
Kumar, S., and Nussinov, R. (2002). Close range electrostatic interactions in proteins close-range electrostatic interactions in proteins. Chembiochem 3, 604–617. doi: 10.1002/1439-7633(20020703)3:7<604::AID-CBIC604>3.0.CO;2-X
Kunisawa, T. (1992). Synonymous codon preferences in bacteriophage T4: a distinctive use of transfer RNAs from T4 and from its host Escherichia coli. J. Theor. Biol. 159, 287–298. doi: 10.1016/S0022-5193(05)80725-8
Kunisawa, T. (2002). Functional role of bacteriophage transfer RNAs: codon usage analysis of genomic sequences stored in the GENBANK/EMBL/DDBJ databases TT – Functional role of bacteriophage transfer RNAs: codon usage analysis of genomic sequences stored in the GENBANK/EMBL/. Data Sci. J. 1, 216–228. doi: 10.2481/dsj.1.216
Kwan, T., Liu, J., DuBow, M., Gros, P., and Pelletier, J. (2005). The complete genomes and proteomes of 27 Staphylococcus aureus bacteriophages. Proc. Natl. Acad. Sci. U.S.A. 102, 5174–5179. doi: 10.1073/pnas.0501140102
Lai, M. J., Lin, N. T., Hu, A., Soo, P. C., Chen, L. K., Chen, L. H., et al. (2011). Antibacterial activity of Acinetobacter baumannii phage ΦaB2 endolysin (LysAB2) against both gram-positive and gram-negative bacteria. Appl. Microbiol. Biotechnol. 90, 529–539. doi: 10.1007/s00253-011-3104-y
Lal, T. M., Sano, M., Hatai, K., and Ransangan, J. (2016). Complete genome sequence of a giant Vibrio phage ValKK3 infecting Vibrio alginolyticus. Genomics Data 8, 37–38. doi: 10.1016/j.gdata.2016.03.002
Lavigne, R., Darius, P., Summer, E. J., Seto, D., Mahadevan, P., Nilsson, A. S., et al. (2009). Classification of Myoviridae bacteriophages using protein sequence similarity. BMC Microbiol. 9:224. doi: 10.1186/1471-2180-9-224
Lee, J.-H., Shin, H., Ji, S., Malhotra, S., Kumar, M., Ryu, S., et al. (2012). Complete genome sequence of phytopathogenic Pectobacterium carotovorum subsp. carotovorum bacteriophage PP1. J. Virol. 86, 8899–8900. doi: 10.1128/JVI.01283-12
Lim, J. A., Shin, H., Heu, S., and Ryu, S. (2014). Exogenous lytic activity of SPN9CC endolysin against gram-negative Bacteria. J. Microbiol. Biotechnol. 24, 803–811. doi: 10.4014/jmb.1403.03035
Lin, H., Kwan, A. L., and Dutcher, S. K. (2010). Synthesizing and salvaging NAD+: lessons learned from Chlamydomonas reinhardtii. PLoS Genet. 6:e1001105. doi: 10.1371/journal.pgen.1001105
Lin, Y. R., and Lin, C. S. (2012). Genome-wide characterization of vibrio phage ϕpp2 with unique arrangements of the mob-like genes. BMC Genomics 13:224. doi: 10.1186/1471-2164-13-224
Lowe, T. M., and Eddy, S. R. (1997). tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 25, 955–964. doi: 10.1093/nar/25.5.0955
Luke, K., Radek, A., Liu, X., Campbell, J., Uzan, M., Haselkorn, R., et al. (2002). Microarray analysis of gene expression during bacteriophage T4 infection. Virology 299, 182–191. doi: 10.1006/viro.2002.1409
Luo, Z. H., Yu, Y. P., Jost, G., Xu, W., and Huang, X. L. (2015). Complete genome sequence of a giant Vibrio bacteriophage VH7D. Mar. Genomics 24, 293–295. doi: 10.1016/j.margen.2015.10.005
Luthy, R., Bowie, J. U., and Eisenberg, D. (1992). Assessment of protein models with three-dimensional profiles. Nature 356, 83–85. doi: 10.1038/356083a0
Mesyanzhinov, V. V., Robben, J., Grymonprez, B., Kostyuchenko, V. A., Bourkaltseva, M. V., Sykilinda, N. N., et al. (2002). The genome of bacteriophage φKZ of Pseudomonas aeruginosa. J. Mol. Biol. 317, 1–19. doi: 10.1006/jmbi.2001.5396
Miller, E. S., Heidelberg, J. F., Eisen, J. A., Nelson, W. C., Durkin, A. S., Ciecko, A., et al. (2003a). Complete genome sequence of the broad-host-range vibriophage KVP40: comparative genomics of a T4-related bacteriophage. J. Bacteriol. 185, 5220–5233. doi: 10.1128/JB.185.17.5220-5233.2003
Miller, E. S., Kutter, E., Mosig, G., Kunisawa, T., Rüger, W., Arisaka, F., et al. (2003b). Bacteriophage T4 genome bacteriophage T4 genome. Microbiol. Mol. Biol. Rev. 67, 86–156. doi: 10.1128/MMBR.67.1.86-156.2003
Moali, C., Anne, C., Lamotte-Brasseur, J., Groslambert, S., Devreese, B., Van Beeumen, J., et al. (2003). Analysis of the importance of the metallo-β-lactamase active site loop in substrate binding and catalysis. Chem. Biol. 10, 319–329. doi: 10.1016/S1074-5521(03)00070-X
Modi, S. R., Lee, H. H., Spina, C. S., and Collins, J. J. (2013). Antibiotic treatment expands the resistance reservoir and ecological network of the phage metagenome. Nature 499, 219–222. doi: 10.1038/nature12212
Moroz, O. V., Murzin, A. G., Makarova, K. S., Koonin, E. V., Wilson, K. S., and Galperin, M. Y. (2005). Dimeric dUTPases, HisE, and MazG belong to a new superfamily of all-α NTP pyrophosphohydrolases with potential “house-cleaning” functions. J. Mol. Biol. 347, 243–255. doi: 10.1016/j.jmb.2005.01.030
Muniesa, M., García, A., Miró, E., Mirelis, B., Prats, G., Jofre, J., et al. (2004). Bacteriophages and diffusion of beta-lactamase genes. Emerg. Infect. Dis. 10, 1134–1137. doi: 10.3201/eid1006.030472
Nakagawa, H., Arisaka, F., and Ishii, S. (1985). Isolation and characterization of the bacteriophage T4 tail-associated lysozyme. J. Virol. 54, 460–466.
Ogata, H., Goto, S., Sato, K., Fujibuchi, W., Bono, H., and Kanehisa, M. (1999). KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 27, 29–34. doi: 10.1093/nar/27.1.29
Oliveira, H., Boas, D. V., Mesnage, S., Kluskens, L. D., Lavigne, R., Sillankorva, S., et al. (2016). Structural and enzymatic characterization of ABgp46, a novel phage endolysin with broad anti-gram-negative bacterial activity. Front. Microbiol. 7:208. doi: 10.3389/fmicb.2016.00208
Oliveira, H., Thiagarajan, V., Walmagh, M., Sillankorva, S., Lavigne, R., Neves-Petersen, M. T., et al. (2014). A thermostable Salmonella phage endolysin, Lys68, with broad bactericidal properties against gram-negative pathogens in presence of weak acids. PLoS ONE 9:e108376. doi: 10.1371/journal.pone.0108376
Ouidir, T., Kentache, T., and Hardouin, J. (2015). Protein lysine acetylation in bacteria: current state of the art. Proteomics 16, 301–309. doi: 10.1002/pmic.201500258
Overbeek, R., Olson, R., Pusch, G. D., Olsen, G. J., Davis, J. J., Disz, T., et al. (2014). The SEED and the rapid annotation of microbial genomes using subsystems technology (RAST). Nucleic Acids Res. 42, 206–214. doi: 10.1093/nar/gkt1226
Petrov, V. M., Nolan, J. M., Bertrand, C., Levy, D., Desplats, C., Krisch, H. M., et al. (2006). Plasticity of the gene functions for DNA replication in the T4-like phages. J. Mol. Biol. 361, 46–68. doi: 10.1016/j.jmb.2006.05.071
Petrov, V. M., Ratnayaka, S., Nolan, J. M., Miller, E. S., and Karam, J. D. (2010). Genomes of the T4-related bacteriophages as windows on microbial genome evolution. Virol. J. 7, 292–310. doi: 10.1186/1743-422X-7-292
Ringel, A. E., Roman, C., and Wolberger, C. (2014). Alternate deacylating specificities of the archaeal sirtuins Sir2Af1 and Sir2Af2. Protein Sci. 23, 1686–1697. doi: 10.1002/pro.2546
Sandegren, L., Nord, D., and Sjöberg, B. M. (2005). SegH and Hef: Two novel homing endonucleases whose genes replace the mobC and mobE genes in several T4-related phages. Nucleic Acids Res. 33, 6203–6213. doi: 10.1093/nar/gki932
Shlomai, J., and Kornberg, A. (1978). Deoxyuridine triphosphatase of Escherichia coli. J. Biol. Chem. 253, 3305–3312.
Sippl, M. J. (1993). Recognition of errors in three-dimensional structures of proteins. Proteins 17, 355–362. doi: 10.1002/prot.340170404
Sitbon, E., and Pietrokovski, S. (2003). New types of conserved sequence domains in DNA-binding regions of homing endonucleases. Trends Biochem. Sci. 28, 473–477. doi: 10.1016/S0968-0004(03)00170-1
Starai, V. J. (2002). Sir2-dependent activation of acetyl-CoA synthetase by deacetylation of active lysine. Science 298, 2390–2392. doi: 10.1126/science.1077650
Stoddard, B. L. (2014). Homing endonucleases from mobile group I introns: discovery to genome engineering. Mob. DNA 5:7. doi: 10.1186/1759-8753-5-7
Stone, R. (2002). Stalin’s forgotten cure. Science 298, 63–66. doi: 10.1126/science.298.5594.728
Sulakvelidze, A. (2011). The challenges of bacteriophage therapy. Ind. Pharm. 45, 14–18.
Tamura, K., Stecher, G., Peterson, D., Filipski, A., and Kumar, S. (2013). MEGA6: molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 30, 2725–2729. doi: 10.1093/molbev/mst197
Torrents, E. (2014). Ribonucleotide reductases: essential enzymes for bacterial life. Front. Cell Infect. Microbiol. 4:52. doi: 10.3389/fcimb.2014.00052
Turner, S., Pryer, K. M., Miao, V. P. W., and Palmer, J. D. (1999). Investigating deep phylogenetic relationships among Cyanobacteria and Plastids by small subunit rRNA sequence analysis. J. Eukaryot. Microbiol. 46, 327–338. doi: 10.1111/j.1550-7408.1999.tb04612.x
Tye, B. K., and Lehman, I. R. (1977). Excision repair of uracil incorporated in DNA as a result of a defect in dUTPase. J. Mol. Biol. 117, 293–306. doi: 10.1016/0022-2836(77)90128-0
Uchiyama, J., Rashel, M., Takemura, I., Wakiguchi, H., and Matsuzaki, S. (2008). In silico and in vivo evaluation of bacteriophage φEF24C, a candidate for treatment of Enterococcus faecalis infections. Appl. Environ. Microbiol. 74, 4149–4163. doi: 10.1128/AEM.02371-07
Ueno, H., and Yonesaki, T. (2004). Phage-induced change in the stability of mRNAs. Virology 329, 134–141. doi: 10.1016/j.virol.2004.08.001
Uzan, M. (2009). RNA processing and decay in bacteriophage T4. Prog. Mol. Biol. Transl. Sci. 85, 43–89. doi: 10.1016/S0079-6603(08)00802-7
Van Asselt, E. J., Dijkstra, A. J., Kalk, K. H., Takacs, B., Keck, W., and Dijkstra, B. W. (1999). Crystal structure of Escherichia coli lytic transglycosylase Slt35 reveals a lysozyme-like catalytic domain with an EF-hand. Structure 7, 1167–1180. doi: 10.1016/S0969-2126(00)80051-9
Wang, J., Jiang, Y., Vincent, M., Sun, Y., Yu, H., Wang, J., et al. (2005). Complete genome sequence of bacteriophage T5. Virology 332, 45–65. doi: 10.1016/j.virol.2004.10.049
Webster, L. C., Zhang, K., Chance, B., Ayene, I., Culp, J. S., Huang, W. J., et al. (1991). Conversion of the E1A Cys4 zinc finger to a nonfunctional His2,Cys2 zinc finger by a single point mutation. Proc. Natl. Acad. Sci. U.S.A. 88, 9989–9993. doi: 10.1073/pnas.88.22.9989
Wheeler, L. J., Rajagopal, I., and Mathews, C. K. (2005). Stimulation of mutagenesis by proportional deoxyribonucleoside triphosphate accumulation in Escherichia coli. DNA Repair 4, 1450–1456. doi: 10.1016/j.dnarep.2005.09.003
Wiederstein, M., and Sippl, M. J. (2007). ProSA-web: interactive web service for the recognition of errors in three-dimensional structures of proteins. Nucleic Acids Res. 35, 407–410. doi: 10.1093/nar/gkm290
Wilson, J. H. (1973). Function of the bacteriophage T4 transfer RNA’s. J. Mol. Biol. 74, 753–757. doi: 10.1016/0022-2836(73)90065-X
Wittmann, J., Dreiseikelmann, B., Rohde, M., Meier-Kolthoff, J. P., Bunk, B., and Rohde, C. (2014). First genome sequences of Achromobacter phages reveal new members of the N4 family. Virol. J. 11:14. doi: 10.1186/1743-422X-11-14
Zafar, N., Mazumder, R., and Seto, D. (2002). CoreGenes: a computational tool for identifying and cataloging “core” genes in a set of small genomes. BMC Bioinformatics 3:12. doi: 10.1186/1471-2105-3-12
Zeng, Q., Bonocora, R. P., and Shub, D. A. (2009). A free-standing homing endonuclease targets an intron insertion site in the psbA gene of cyanophages. Curr. Biol. 19, 218–222. doi: 10.1016/j.cub.2008.11.069
Zerbino, D. R., and Birney, E. (2008). Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 18, 821–829. doi: 10.1101/gr.074492.107
Zhao, K., Chai, X., and Marmorstein, R. (2004). Structure and substrate binding properties of cobB, a Sir2 homolog protein deacetylase from Escherichia coli. J. Mol. Biol. 337, 731–741. doi: 10.1016/j.jmb.2004.01.060
Keywords: bacteriophages (phages), phage therapy, Vibrio, sirtuins, comparative genomics, phage–host interaction, nucleotide metabolism, NAD+-dependent deacetylation
Citation: Skliros D, Kalatzis PG, Katharios P and Flemetakis E (2016) Comparative Functional Genomic Analysis of Two Vibrio Phages Reveals Complex Metabolic Interactions with the Host Cell. Front. Microbiol. 7:1807. doi: 10.3389/fmicb.2016.01807
Received: 25 August 2016; Accepted: 27 October 2016;
Published: 14 November 2016.
Edited by:Akio Adachi, Tokushima University, Japan
Reviewed by:Jumpei Uchiyama, Azabu University, Japan
Weili Liang, National Institute for Communicable Disease Control and Prevention, Chinese Center for Disease Control and Prevention, China
Andrei A. Zimin, Institute of Biochemistry and Physiology of Microorganisms (Russian Academy of Sciences), Russia
Copyright © 2016 Skliros, Kalatzis, Katharios and Flemetakis. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Emmanouil Flemetakis, firstname.lastname@example.org