Skip to main content

Advertisement

Functional and comparative genome analysis of novel virulent actinophages belonging to Streptomyces flavovirens

Article metrics

Abstract

Background

Next Generation Sequencing (NGS) technologies provide exciting possibilities for whole genome sequencing of a plethora of organisms including bacterial strains and phages, with many possible applications in research and diagnostics. No Streptomyces flavovirens phages have been sequenced to date; there is therefore a lack in available information about S. flavovirens phage genomics. We report biological and physiochemical features and use NGS to provide the complete annotated genomes for two new strains (Sf1 and Sf3) of the virulent phage Streptomyces flavovirens, isolated from Egyptian soil samples.

Results

The S. flavovirens phages (Sf1 and Sf3) examined in this study show higher adsorption rates (82 and 85%, respectively) than other actinophages, indicating a strong specificity to their host, and latent periods (15 and 30 min.), followed by rise periods of 45 and 30 min. As expected for actinophages, their burst sizes were 1.95 and 2.49 virions per mL. Both phages were stable and, as reported in previous experiments, showed a significant increase in their activity after sodium chloride (NaCl) and magnesium chloride (MgCl2.6H2O) treatments, whereas after zinc chloride (ZnCl2) application both phages showed a significant decrease in infection.

The sequenced phage genomes are parts of a singleton cluster with sizes of 43,150 bp and 60,934 bp, respectively. Bioinformatics analyses and functional characterizations enabled the assignment of possible functions to 19 and 28 putative identified ORFs, which included phage structural proteins, lysis components and metabolic proteins.

Thirty phams were identified in both phages, 10 (33.3%) of them with known function, which can be used in cluster prediction. Comparative genomic analysis revealed significant homology between the two phages, showing the highest hits among Sf1, Sf3 and the closest Streptomyces phage (VWB phages) in a specific 13Kb region. However, the phylogenetic analysis using the Major Capsid Protein (MCP) sequences highlighted that the isolated phages belong to the BG Streptomyces phage group but are clearly separated, representing a novel sub-cluster.

Conclusion

The results of this study provide the first physiological and genomic information for S. flavovirens phages and will be useful for pharmaceutical industries based on S. flavovirens and future phage evolution studies.

Background

Bacteriophages (phages), natural viral predators of bacteria, are engaged in a constant evolutionary arms race with their hosts [1], playing major roles in the ecological balance of microbial life and in microbial diversity.

Most double-stranded DNA (dsDNA) phages share the same gene pool [2]; however, sequence comparisons reveal a widespread horizontal exchange of sequences among genomes, mediated by both non-homologous and homologous recombination. High frequency exchange among phages occupying similar ecological niches results in a high rate of mosaic diversity in local populations [3]. Studies confirm that phage genomes are mosaics and represent a large common genetic pool due to horizontal exchange [4, 5].

The screening of microbial natural products continues to constitute an important route to the discovery of chemicals for developing new therapeutic agents and evaluating the therapeutic potential of bacterial taxa [68]. In this respect, actinomycetes are a group of microorganisms mostly used in biotechnology for handling bioactive compounds. [9, 10]. Moreover, bacteriophages can be used to detect antiviral compound production by actinomycetes. Finally, actinophages are isolated and investigated because they can influence antibiotic production in bacterial strains, causing problems in the pharmaceutical industry. The vast majority of actinophages were isolated from sediments, but direct isolation from soil generally yields extremely low titers [11, 12]. However, although it is difficult to grow bacteriophages from soil without enrichment, a wide range of counts has been reported [13, 14].

Recently, there has been expanding interest in bacteriophages that infect Streptomyces species, since the phages can support the development of cloning vectors [15]. Such vectors could open the way for genetic manipulation as an important tool for Streptomyces improvement. Moreover, the mechanisms of the system for phage infection and multiplication could be useful in the fermentation industry and lead to the development of phage cloning vectors [16]. To date, no studies on phages isolated from S. flavovirens, an important source for several pharmaceutical drugs, such as actinomycin complex, mureidomycin and pravastatin [17, 18], have been carried out.

The development of high-throughput NGS (Next Generation Sequencing) technologies [19, 20] and the possibility to sequence entire genomes or transcriptomes more efficiently and economically than with first generation sequencing strategies permitted the collection of large amounts of information and the analysis of sequences from hundreds of thousands of species. Therefore, the dawn of next generation sequencing technologies has opened up exciting possibilities for whole genome sequencing in a wide range of organisms and the bacterial viruses have not been excluded from this revolution, despite the fact that their genomes are orders of magnitude smaller in size compared with bacteria and other organisms.

The Actinophage Sequence Databases (http://phagesdb.org/) currently include 5861 genomes from putative actinophages, 120 of which infect Streptomyces species and sixty-five of which are sequenced, but no genomes of phages isolated form S. flavovirens are currently available. The NCBI genome database contains around 600 Caudovirales genomes to date but the number of complete bacteriophage genomes published is growing slowly [21].

Until now, no phages belonging to S. flavovirens have been sequenced and relatively little is known about S. flavovirens phage genomics. In the present work, we report the first whole genome sequencing study and annotation of two S. flavovirens virulent phages. The results will provide an important genomic resource for future investigations in the bacteriophages related to S. flavovirens and for phage evolution studies.

Methods

Source of lytic actinophages

Two isolates of Streptomyces flavovirens phages, named Sf1 and Sf3, were obtained from the virology lab, Agric. Microbiology Department, Faculty of Agriculture, Ain Shams University, Cairo, Egypt. Phages were isolated from soil and the morphological properties were analyzed by standard methodology and reported in Marei and Elbaz (2013) [22].

Purification of lytic actinophages

The high titer phage suspension of each isolated phage was prepared using a liquid culture enrichment technique. The high titer phage suspension of each phage was ultra-centrifuged at 30000 rpm for 90 min. at 4 °C in a Beckman L7-35 ultracentrifuge. The pellet was gently resuspended in 0.5 ml of 0.2 M phosphate buffer pH 7.2 [23].

Adsorption rate and one-step growth experiments

The adsorption experiments were carried out with two isolated phage suspensions added to spores of their indicator host (S. flavovirens). Suspensions of each phage were incubated at 30 °C with gentle shaking. Samples were withdrawn at regular intervals after inoculation.

The mycelial fragments of the indicator strain were removed by centrifugation and the concentration of phage remaining in the supernatant was counted. The adsorption rates of the two phages were determined by measuring residual plaque-forming ability in membrane-filtered samples of an attachment mixture [24] and the adsorption rate constant k (mL/min) was calculated [25]. The one-step growth experiment was performed as described by Dowding (1973) [24].

Physiochemical stability

To evaluate the phages’ stability three different chemicals (NaCl, MgCl2.6H2o and ZnCl2), were used. Five concentrations (0.1, 0.2, 0.3, 0.4 and 0.5 mM) for each salt were employed [26]. To test the effect of different treatments phage solutions for both tested strains with final concentrations of 107 PFU/ml were utilized. The mixture was incubated for 10 min at room temperature (RT). The number of plaques was determined using the double layer method (plaque assay test) [27]. A control test was prepared by mixing bacterial suspension with phage without the tested chemicals.

DNA isolation, library preparation and whole genome sequencing

Genomic DNA was isolated from the propagated phages according to the procedure described by Kieser et al. [28]. DNA quality was assessed using a Nanodrop Bioanalyzer ND1000 (ThermoScientific). Sequencing libraries were prepared by shearing 1 μg of DNA in blunt-ended fragments by linking the Ion adapters using an Ion XpressTM Plus Fragment Library Kit (Life Technologies, Carlsbad, USA) according to the manufacturer’s specifications. The sized and ligated fragments were amplified by emulsion-PCR using the Ion OneTouch 200 Template kit (Life Technologies, Carlsbad, USA). Quality and insert size distribution were assessed using an Agilent Bioanalyzer DNA 1000 chip. Libraries were sequenced on an Ion Torrent PGM semiconductor sequencer (Life Technologies, Carlsbad, USA) using the 200 bp protocol and an Ion Torrent 314 chip following the manufacturer instructions (Life Technologies, Carlsbad, USA).

Assembly and bioinformatics analyses

Raw reads resulting from Sf1 and Sf3 sequencing were trimmed using Trimmomatic with single end mode (no quality encoding was specified to allow the program to determine it automatically [29]) and assembled separately using the gsAssembler (Roche Applied Science, Indianapolis, IN); the Graphical User Interface (GUI) version was used with the default parameters. The collected contigs were visualized and validated using Hawkeye [30]. Resulting contigs for each phage showed approximately 60-fold sequence read coverage. The expected sequence accuracy was 95% with a statistical error of less than 1 in 10,000 bp. Sequence homologies were determined by using BLASTn against the actinophage database to assign the phages to a cluster [31].

Open reading frame (ORF) analysis and gene prediction

Open reading frames (ORFs) were identified and the genome sequences of each phage were annotated as described previously in Dobbins et al., 2004 by using DNA Master (J. G. Lawrence) (http://cobamide2.bio.pitt.edu) software and visual inspection [32]. For a genome-wide viewpoint an association with the annotation refinement, functional analysis and other explorations was developed using Phamerator. Protein sequence relationships and conserved domains within genes were also studied. Gene products were grouped into “Phamilies” generally referred to as “Phams”, or groups of proteins with a high degree of similarity to one another. The pairwise alignment scores and significant rate were determined using BLASTp and ClustalW [33].

Genomic comparisons between the sequenced and the close related phages

Sequence comparisons were performed by using the BLAST algorithm available at NCBI [34] and Mauve software [35]. A comparison map among Sf1 and Sf3 Streptomyces phages and closely related phages (VWB and SV1) with available genomes in the National Center for Biotechnology Information (NCBI) nucleotide database (https://www.ncbi.nlm.nih.gov/) was generated by Circoletto (http://tools.bat.infspire.org/circoletto/) [34, 36]. For pictogram construction, bit-score values were used to describe the quality of the alignment at a given point. The bit-score is a normalized version of the score value returned by the BLAST searches, expressed in bits [37].

The phylogenetic tree of Major Capsid Protein (MCP) genes from two new isolated phages (Sf1 and Sf3) and 20 related Streptomyces phages available in the NCBI database was constructed with Geneious software version (R8) (http://www.geneious.com) [38] based on the Neighbor-Joining (NJ) algorithm.

Results and discussion

Adsorption rate constant and growth characteristics of isolated phages

Adsorption of Sf1 and Sf3 was determined using S. flavovirens cells grown in phage medium to the early exponential phase of growth (15-h cultures). About 82 and 85% of all infective Sf1 and Sf3 particles, respectively, were adsorbed within 20 min of contact. The adsorption reached a maximum after 30 min. for both phages. The adsorption constant K was 3.66 pL/min for Sf1 and 3.80 pL/min for Sf3, determined by the Adams’s formula [27]. The phages adsorption rates were higher than other actinophages [39], which was probably due to the strong specificity of the Sf1 and Sf3 phages to their host.

The production of Sf1 and Sf3 phages were determined in a one-step growth experiment at 30 °C. Results revealed that the latent periods of Sf1 and Sf3 were approximately 15 and 30 mins, respectively. After 30 and 45 mins the maximum rise period was shown and the burst sizes were 1.95 and 2.49 PFU/mL for Sf1 and Sf3, respectively (Fig. 1). The present results are in agreement with the data obtained from a study on 24 actinophages [40], underlining that under controlled cultural conditions the infection of isolated Streptomycetes cells by phages was varied.

Fig. 1
figure1

One-step growth experiment for Sf1 and Sf3 phages development on S. flavovirens at 30 °C

Physiochemical stability of isolated actinophages

Sodium and magnesium chloride treatments yielded a significant increase in both phages’ activity for all concentrations used compared with the control, while zinc chloride application with concentrations > 0.3 mM caused a significant decrease of activity for Sf1 and Sf3 (Fig. 2). Similar results were reported in previous studies [4143]. Absence of calcium and magnesium ions prevents adsorption and the lysis cycle, while their presence stimulates a significant increase in phage activity, probably due to the increase of adsorption and penetration rates. On the contrary, zinc and aluminum chloride showed significant loss of infectivity in both phages. This is in accordance with the experiments performed by Robert and Charles, which suggested that aluminum caused viral inactivation related to the dissociation of viral capsid proteins [44].

Fig. 2
figure2

Effect of three different chemicals on the Sf1 and Sf3 infectivity

Genome organization of phages Sf1 and Sf3

Genome sequencing generated 69,719 and 107,273 reads for each phage with around 60-fold coverage and 43,150 bp, and 60,934 bp assembled sequences for Sf1 and Sf3, respectively. The pair-wise alignment [45] revealed that the genomes of Sf1 and Sf3 shared an overall high level of similarity, with conserved regions of high identity (100% identity) interspersed between regions with high variability (ranging from 23.9% to 87.5%) (Fig. 3a). A similar mosaic genome structure has been observed in most other phage genomes, indicating extensive horizontal genetic exchange among phages [4649]. No close relatives (Singleton) from modeling of both genome construction were revealed (Fig. 1).

Fig. 3
figure3

Genomic organization of Sf1 and Sf3 phages. Phages were mapped using Phamerator; the purple lines between phages underline the regions with high similarity, while the ruler corresponds to genome base pairs. The predicted genes are shown as boxes either above or below the genome (ruler), depending on whether are rightwards or leftwards transcribed, respectively. Gene numbers are shown within each box; pink boxes refereed to the genes with high similarity between two phages while the blue boxes refereed to the genes that show low similarity. a The phages maps showed by cluster conservation. b The phages maps showed by phams; genes are colored according to their function categories “phams”

Annotation of Sf1 and Sf3 genomes revealed 52 and 91 putative ORFs, respectively. According to their homology, 19 out of 52 ORFs (36.5%) from Sf1 and 28 out of 91 (30.8%) from the Sf3 genome have been assigned functions compared to known conserved domains [50] (Tables 1 and 2). Isolated genes were mainly involved in DNA replication and repair, nucleotide metabolism, lysis, phage structural proteins and other enzymes. The results obtained are in agreement with other bacteriophage studies [5153]. Phage Sf1 showed 52 ORFs (Table 1), named gp1 - gp52, while 91 ORF were identified from Phage Sf3, from gp1 to gp91 (Table 2). The majority of members of identified families are bacteriophage proteins, while others (75%) have unknown function [54, 55].

Table 1 Overview of Sf1 phage ORFs, summary of homology searches and annotations
Table 2 Overview of Sf3 phage, ORFs, summary of homology searches and annotations

Phage structure and assembly genes

Several genes code for terminase subunit proteins, such as gp1 and 2 which code for terminase_4 (pfam05119) and terminase_1 (pfam03354) super-families, respectively. The gp3 and gp23 genes encode for the phage portal protein (pfam05133), an important protein involved in DNA transport during its packaging and ejection. Another relevant gene is gp6 which, together with gp27,codes for the major capsid protein (PHA00665) [56] and the major capsid protein E domain (pfam03864) [57], respectively, involved in the stabilization of the condensed form of DNA in phage heads. Some genes involved in tail development, gp14 (pfam10145) and gp17 (pfam13550), were also identified.

In Sf3we found a gene (gp3) encoding phage portal protein (pfam05133), crucial for DNA migration and building the junction between head and tail proteins [58], and others, such as gp7 and gb16, that encode for the major capsid protein E domain (pfam03864) [57] or for lyase (gp21), like pectate lyase_3 superfamily protein (pfam12708). A phage putative head morphogenesis protein (TIGR01641) of 110 amino acids found exclusively in phage-related proteins, was encoded by gp84. Putaive head morphogenesis proteins such as gp85, which encodesthe transcriptional activator RhaR (PRK13502), and gp89, involved in the phage terminase_3 (COG1783) synthesis, are activated during the beginning of double-stranded viral DNA packaging [59].

DNA replication and metabolic genes

The gp44 gene encodes YabA (COG4467), a protein that interacts with the DnaA initiator and the DnaN sliding clamp and drives the control of DNA replication initiation [60, 61]. gp46 and gp52 encode for helix-turn-helix XRE-family like proteins (cd00093) [62] and histidine kinase-like ATPases (cd00075) [63], respectively, two important binding proteins with roles in the replication, repair, storage and modification of DNA. gp4 encodes a protein belonging to the MATE family (cd13126), which functions as a translocase for lipopolysaccharides [64], while gp5 codes for the golgin subfamily protein A5, a protein responsible for maintaining Golgi structure in intra-Golgi retrograde transport [65].

ORFs with the same biological roles were also identified in Sf3 phage. Indeed gp35 encodes for a HhH-GPD superfamily base excision DNA repair protein (pfam00730). This group includes endonuclease III, 8-oxoguanine DNA glycosylases and DNA-3-methyladenine glycosylase II [66]. Other members include different types of DNA and RNA exonucleases such as RNase T, oligoribonuclease, and RNA exonuclease (REX) [67]; Holliday junction resolvases (HJRs) (cd00529), endonucleases structurally similar to RNase H and Hsp70, which specifically resolve Holliday junction DNA intermediates during homologous recombination was encoded by gp52 [68]. Gp76 encodes for HNH nucleases (cd00085), an endonuclease signature which is found in viral, prokaryotic and eukaryotic proteins [69].

Cell lysis genes

Crucial genes implicated in lysis activities, such as the cell wall degradation process in bacteria during host infection, were identified in the Sf1 genome. Indeed, gp36 encodes for the lytic transglycosylase (LT) (cd00254) that catalyzes the cleavage of the beta-1,4-glycosidic bond between N-acetylmuramic acid and N-acetyl-D-glucoseamine, similar to “goose-type” lysozymes. gp42 encodespeptidoglycan recognition proteins (PGRPs) (cd06583), namely receptors that bind and hydrolyze peptidoglycans of bacterial cell walls, and contains two conserved histidines and a cysteine, typical residues of zinc binding sites [70].

While gp21 is included in the pectate lyase superfamily (pfam12708), proteins with a beta helical structure like pectate lyase and most closely related to glycosyl hydrolase family and gp22 encodes to Peptidoglycan recognition proteins (PGRPs) (cd06583) [70], were identified in Sf3 genome.

Both phage genomes show up to bring a modular organization, with genes of related function clustered together (Fig. 3a and b). DNA sequences of the first 13 kb in Sf3 are highly similar to the last DNA sequences in Sf1 and encode for DNA packaging structural proteins (Fig. 3b).

On the basis of the amino acid sequence similarity between the gene products, the conserved pfam05133 motif and the gene locations, orf3 is predicted to encode a portal protein in both phages. No small terminase-encoding gene could be identified in either genome. The largest gene in Sf1 genome is located in orf36 (3.5 kb) encoding the lytic transglycosylase (LT), while the largest one in Sf3 genome with the same length is orf16, encoding the major capsid protein E domain. [48, 71, 72]. A possible lyase gene is positioned distinctively in both phage genomes (orf41 for Sf1 and orf21 for Sf3). Those genes located downstream in both phage genomes encode proteins involved in DNA synthesis, metabolism and repair (Fig. 3b).

Evolutionary relationship of Sf1 and Sf3

Sf1 and Sf3 phages show 30 phams, where 29 out of 30 phams contain two members (Table 3), while three members belong to pham number 12. Ten phams (33.3%) were assigned with known functionality; the others are unknown. Therefore, some of these phams are informative and can be used in evolutionary studies. Indeed, as reported for mycobacteriophages [73], single, ubiquitous, semi-conserved genes can be utilized for cluster prediction, useful when the whole genome sequence is unavailable. The 30 identified phams, which include important genes (see below), underline a close phylogenetic relationship between the two isolated phages and provide important information that can be used in future evolutionary relationship studies by comparing the genes identified in the Streptomyces flavovirens phages and homologous genes in other bacteriophages.

Table 3 Related Conserved Domains (CD) to the detected Phamilies using Phamerator

orf27 (Sf1) and orf7 (Sf3) as members of pham n.7 were assigned as phage major capsid protein (MCP) E domains; this important class of genes was also used as a single gene prediction system for the mycobacteriophage clusters analysis [73]. orf23 (Sf1) and orf3 (Sf3), members of pham n. 3, were classified as phage portal proteins. These proteins were used in some previous investigations as a marker of diversity indicating, in some cases, the connections between habitat properties, microbial community structure and phage community composition [74]. orf29 (Sf1) and orf9 (Sf3) are the members of pham n.9, were assigned to phage protein gp19, an important tail component. Most of the proteins forming the phage tail components as well as other needle-like assemblies (e.g. secretion systems and bacteriocins) have a common origin from a single protein module [74]. This evidence emphasizes the importance of phage protein diversification and specialization in the evolution of different and complex bacterial systems and in bacterial adaptation, developing new functions and providing a distinct selective advantage [74].

As expected, the virulent phages developed phams involved in lysogenic pathways. Indeed, orf41 (Sf1) and orf21 (Sf3), grouped in pham n.20, showed high homology to the pectate lyase superfamily protein that can modify the properties of polysaccharides. Since the pectinolytic protein family is commonly represented in prokaryotic and eukaryotic microorganisms and, in plants, is involved in remodelling cell walls, it is clear that the divergence from the ancestral protein over time has allowed different micro-organisms to target a range of pectin-like substrates while the overall structure has been maintained [75]. orf42 (Sf1) and orf22 (Sf3) are members of pham n.21 and classified as peptidoglycan recognition proteins (PGRPs), an innate class of immunity molecules present in insects, mollusks, echinoderms, and vertebrates that by interacting with peptidoglycan in the cell wall, rather than permeabilizing bacterial membranes, kills bacteria. These proteins were reported, at least in one carboxy-terminal domain, as homologous in bacteriophage and bacteria [76]. orf46 (Sf1) and orf26 (Sf3) are grouped in pham n.24 and were identified as helix-turn-helix (HTH) XRE-family-like proteins, one of the early studied regulatory DNA-binding proteins involved in metabolic regulation in bacteria. This class of genes encodes components to process environmental metabolites (e.g. lactose) and to produce interacting constituents in the development of a lytic or lysogenic pathway in phages. A common ancestor for all DNA-binding domains was suggested and, through its duplication and divergence, the diversity of transcription regulators that drive bacterial and phage genes was generated. The HTH fold investigations confirmed the significance of this module in DNA–protein interactions across a wide phylogenetic spectrum including a wide variety of phages [77].

orf26 (Sf1) and orf6 (Sf3), members of pham n. 6, were classified as bacteriophage lambda head decoration protein D. Since the protein allows for the display of many copies of a foreign protein, which is advantageous for displaying weak ligands for affinity selection, a useful platform for phage polypeptide display was recently developed [78]. Interestingly, orf32 in Sf1 and orf12 in Sf3 were not assigned functions previously, although they belong to the pham n. 12 together with orf 1 (Sf1) which is classified as terminase_4.

A standard Nucleotide NCBI BLAST (blastn) search was developed using both Sf1 and Sf3 phage whole genome sequences as a query against a non-redundant nucleotide sequences database. Starting from a whole phage dataset (https://www.ncbi.nlm.nih.gov/) the available phage genomes with the best identity percentages (VWB and SV1) were chosen and a pictogram was developed (Fig. 4). Seventy-eight percent identity for both S. flavovirens phages compared to the complete genome of bacteriophage VWB, isolated from S. venezuelae strain ETH 14630 (AY320035.2), was exhibited (with 29% and 36% of coverage for Sf1 and Sf3, respectively), while 75% of identity for both studied phages with S. venezuelae phage SV1 (JX182371.1) was reported, but with low query coverage (11% for Sf1 and 14% for Sf3), probably due to the phylogenetic distance between the compared phages.

Fig. 4
figure4

Sequence similarities among Sf1, Sf3, VWB and SV1 phages. The picture shows the results of the BLAST local alignments using Sf1 and Sf3 as a query against the VWB and SV1 phages sequences. The different colours (blue, green, orange and red) represent the overall quality of the aligned segments along the phage sequences, evaluated on the basis of the bit-score values from the worst to the best score (blue to red). The bit-score is a normalized version of the score value obtained by BLAST searches, expressed in bits. The height of the coloured bars in the histogram shows how many times each colour hits a specific fragment of the other phage sequences. A twist in a ribbon indicates that the local alignment is inverted (query and database sequence on opposite strands)

The alignment of both Sf1 and Sf3 genomes against the sequences of VWB phage, carried out by Mauve software, revealed that most hits occurred around a 13Kb region (Fig. 4). The approximate location of this region were (18000–31000) within the Sf1 genome, (1–13000) in the Sf3 genome and (23000–36000) in the VWB genome. On the contrary, the alignment of both S. flavovirens phage genomes versus the sequences of SV1 showed only a short region (~1Kb) with moderate bit score ranging from 9691–10707 and 10300–11208 in the genomes of Sf1 and Sf3, respectively, consistent with the low sequence coverage obtained.

The MCPs diversity between Sf1, Sf3 and 20 related Streptomyces phages, due to a combination of illegitimate and homologous recombination [79] and mutational drift, was also evaluated. The current investigation highlighted the hybrid generation between phage genera [80] or phage families [81]. Twenty-two Streptomyces phages were grouped in five main branches (Fig. 5). The Lannister MCP shared a close evolutionary relationship with the Izzy, Aaronocolus, and Caliburn sequences, demonstrating that phages may undergo genetic exchange by horizontal gene transfer from a large shared pool [4] and that horizontal gene transfer between phages is a component of their evolution. Numerous gene exchanges within each major clade and core phage functions do not appear to have co-evolved with specific hosts [82].

Fig. 5
figure5

Phylogenetic analysis of studied phages and other members (20) of the Streptomyces phages group based on MCPs. Bootstrap values indicate the number of times a node was supported in 1000 resampling replications

Our phylogenetic analysis is useful for further studies, since both Sf1 and Sf3 were recovered in a clade that included phages that infect Streptomyces species but most of these phages (Maih, YDN12, Xkcd426 and TP1604) were members of the BG phage cluster; this clustering does not represent a phylogenetic or taxonomic grouping but rather provides a framework for reflecting their overall genome relationships and for identifying genes that have been recently exchanged and their genomic context [83, 84]. Moreover, Sf1 and Sf3 grouped in a separate branch, indicating that isolated phages belong to the BG phage cluster but represent a different sub-cluster.

Conclusion

Recently, large advances have occurred in phage genomics; nevertheless,the full extent of phage diversity and evolutionary pathways are yet unknown. With the advent of NGS technologies a much greater volume of transcriptome and genome sequences is available and we can therefore expect an increased flow of new data in upcoming years. Current assessment suggests that more than 1031 phages exist on earth, representing more than ten million phage “species”. Of these, less than 6000 have been observed using electron microscopy and fewer than 1000 genomes have been sequenced. The available sequences show that the majority of phages analyzed are tailed phages belonging to the family Siphoviridae, but less is known about the degree of their genetic diversity. The genomic characterization of phages is necessary to evaluate their important ecological impact. In spite of their ubiquity, phages have not yet been characterized for many bacterial genera. In the present study, biological, physiochemical and genome sequences of two new virulent Streptomyces phages are presented, representing the first genomic report of S. flavovirens phages which may represent a new sub-cluster of the BG Streptomyces phage cluster.

Abbreviations

dsDNA:

Double-stranded DNA

HJRs:

Holliday junction resolvases

HTH:

Helix-turn-helix

LT:

Lytic transglycosylase

MCP:

Major capsid protein

NCBI:

National Center for Biotechnology Information

NGS:

New Generation Sequencing

ORFs:

Open reading frames

PGRPs:

Peptidoglycan recognition proteins

Phages:

Bacteriophages

REX:

RNA exonuclease

TDP:

Thymidine diphosphate

TMP:

Thymidine monophosphate

TMPK:

Thymidine monophosphate kinase

TTP:

Thymidine triphosphate

References

  1. 1.

    Klumpp J, Fouts DE, Sozhamannan S. Next generation sequencing technologies and the changing landscape of phage genomics. Bacteriophage. 2012;2(3):190–9.

  2. 2.

    Sandmeier H. Acquisition and rearrangement of sequence motifs in the evolution of bacteriophage tail fibres. Mol Microbiol. 1994;12(3):343–50.

  3. 3.

    Juhala RJ, Ford ME, Duda RL, et al. Genomic sequences of bacteriophages HK97 and HK022: pervasive genetic mosaicism in the lambdoid bacteriophages. J Mol Biol. 2000;299:27–51.

  4. 4.

    Hendrix RW, Smith MCM, Burns RN, et al. Evolutionary relationships among diverse bacteriophages and prophages: all the world’s a phage. Proc Natl Acad Sci U S A. 1999;96:2192–7.

  5. 5.

    Hendrix RW. Bacteriophage genomics. Current Opinion in Microbiol. 2003;6:506–11.

  6. 6.

    Bradley SG, Ritzi D. Electron microscopical studies of actinophage multiplication. J Gen Virol. 1967;1(3):285–90.

  7. 7.

    Bull AT, Ward AC, Goodfellow M. Search and discovery strategies for biotechnology: the paradigm shift. Microbiol Mol Biol Rev. 2000;64:573–606.

  8. 8.

    Ipek KD. Actinophages as indicators of actinomycete taxa in marine environments. Antonie Van Leeuwenhoek. 2005;87:19–28.

  9. 9.

    Korn-Wendish F, Schneider J. Phage typing-a useful tool in actinomycete systematics. Gene. 1992;115:243–7.

  10. 10.

    Williams ST, Locci R, Beswick A, et al. Detection and identification of novel actinomycetes. Res Microbiol. 1993;144(8):653–6.

  11. 11.

    Ackerman HW, Berthiaume L, Jones LA. New actinophage species. InterVirology. 1985;23:121–30.

  12. 12.

    Othman BA, Askora AA, Awny NM, et al. Virulent Bacteriophages for Streptomyces griseoflsvus Isolated from Soil. Pak J Biotechnol. 2008;5:109–18.

  13. 13.

    Ashelford KE, Day MJ, Bailey MJ, et al. In situ population dynamics of bacterial viruses in a terrestrial environment. Appl Environ Microbiol. 1999;65:169–74.

  14. 14.

    Ashelford KE, Day MJ, Fry JC. Elevated Abundance of Bacteriophage Infecting Bacteria in Soil. Appl Environ Microbiol. 2003;69(1):285–9.

  15. 15.

    Chater KF, Hopwood DA, Kieser T, et al. Gene cloning in Streptomyces. Curr Top Microbiol Immunol. 1982;96:69–95.

  16. 16.

    Harris JE, Chater KF, Bruton CJ, et al. The restriction mapping of c gene deletions in Streptomyces bacteriophage 4C31 and their use in cloning vector development. Gene. 1983;22(2–3):167–74.

  17. 17.

    Oliver K. Grundwissen Pharmazeutische Biotechnologie. Vieweg + Teubner Verlag: Wiesbaden; 2002. ISBN 3-322-80124-1.

  18. 18.

    Satinder KB, Gurpreet SD, Carlos RS. Biotransformation of Waste Biomass into High Value Biochemicals. New York, NY: Springer New York; 2014.

  19. 19.

    Eisenstein M. Oxford Nanopore announcement sets sequencing sector abuzz. Nat Biotechnol. 2012;30:295–6.

  20. 20.

    Srivatsan A, Han Y, Peng J, et al. High-precision, whole-genome sequencing of laboratory strains facilitates genetic studies. PLoS Genet 2008, doi:10.1371/journal.pgen.1000139.

  21. 21.

    Hatfull GF. Bacteriophage genomics. Curr Opin Microbiol. 2008;11:447–53.

  22. 22.

    Marei EM, Elbaz RM. Isolation and Molecular Characterization of Three Virulent Actinophages Specific for Streptomyces flavovirens. J Virol Res. 2013;2:12–7.

  23. 23.

    Lanning S, Williams ST. Methods for the direct isolation and enumeration of actinophages in soil. J Gen Microbiol. 1982;128:2063–71.

  24. 24.

    Dowding JE. Characterization of a bacteriophage virulent for Streptomyces coelicolor A3(2). J Gen Microbiol. 1973;76:163–76.

  25. 25.

    ADAMS MH. Methods of study of bacterial viruses. Methods in Med Research. 1950;2:1–73.

  26. 26.

    Paunikar WN, Sanmukh SG, Ghosh TK. Effect of metal ions and chemical solvents on the adsorption of Salmonella phage on Salmonella choleraesuis subspecies indica. Int J of Pharma and Biosci. 2012;3(1):181–90.

  27. 27.

    Adams MH. Bacteriophages. New York: Interscience; 1959.

  28. 28.

    Kieser T, Bibb MJ, Buttner MJ, et al. Practical Streptomyces Genetics. Norwich: The John Innes Foundation; 2000.

  29. 29.

    Lohsem M, Bolger AM, Nagel A, et al. RobiNA: a user-friendly, integrated software solution for RNA-Seq-based transcriptomics. Nucleic Acids Res 2012, doi:10.1093/nar/gks540.

  30. 30.

    Schatz MC, Phillippy AM, Shneiderman B, et al. Hawkeye: an interactive visual analytics tool for genome assemblies. Genome Biology 2007, doi:10.1186/gb-2007-8-3-r34.

  31. 31.

    Hatfull GF, Jacobs-Sera D, Lawrence JG, et al. Comparative genomic analysis of sixty mycobacteriophage genomes: Genome clustering, gene acquisition and gene size. J Mol Biol. 2010;397(1):119–43.

  32. 32.

    Dobbins AT, George Jr M, Basham DA, et al. Complete Genomic Sequence of the Virulent Salmonella Bacteriophage SP6. J Bacteriol. 2004;186(7):1933–44.

  33. 33.

    Cresawn SG, Bogel M, Day N, et al. Phamerator: a bioinformatic tool for comparative bacteriophage genomics. BMC Bioinformatics. 2011;12:395.

  34. 34.

    Altschul SF, Gish W, Miller W, et al. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.

  35. 35.

    Darling AE, Mau B, Perna NT. ProgressiveMauve: multiple genome alignment with gene gain, loss and rearrangement. PLoS One 2010, doi:10.1371/journal.pone.0011147.

  36. 36.

    Krzywinski M, Schein J, Birol I, et al. Circos: an information aesthetic for comparative genomics. Genome Res. 2009;19(9):1639–45.

  37. 37.

    Darzentas N. Circoletto: visualizing sequence similarity with Circos. Bioinformatics. 2010;26(20):2620–1.

  38. 38.

    Kearse M, et al. Geneious basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28:1647–9.

  39. 39.

    Anné J, Van Mellaert L, Decock B, et al. Further biological and molecular characterization of actinophage VWB. J Gen Microbiol. 1990;136:1365–72.

  40. 40.

    Pringsulaka O, Chvanich S, Doi K, et al. Characteristics and Distribution of Actinophages Isolated from Thai Soil Samples. Sci Asia. 2004;30:223–30.

  41. 41.

    AlKhazindar M, Sayed ETA, Khalil MS, et al. Isolation and Characterization of Two Phages Infecting Streptomyces scabies. Res J Pharm, Biol Chem Sci. 2016;7(2):1351–63.

  42. 42.

    Binetti A, Quiberoni A, Reinheimer J. Phage adsorption to Streptococcus thermophilus. Influence of environmental factors and characterization of cell-receptors. Food Res Int. 2002;35:73–83.

  43. 43.

    Quiberoni A, Guglielmotti D, Binetti A, et al. Characterization of three Lactobacillus delbrueckii sub sp. bulgaricus phages and the physiochemical analysis of phage adsorption. J Appl Microbiol. 2004;96:340–51.

  44. 44.

    Robert B, Charles PJ. Characterization of the effect of aluminum metal on poliovirus. J of Indus Microbiol. 1988;3:33–8.

  45. 45.

    Brudno M, Do CB, Cooper GM, et al. NISC Comparative sequencing program LAGAN and multi-LAGAN: efficient tools for large-scale multiple alignment of genomic DNA. Genome Res. 2003;13(4):721–31.

  46. 46.

    Casjen SR, Thuman-Commike PA. Evolution of mosaically related bacteriophage genomes seen through the lens of phage P22 virion assembly. Virology. 2011;411:393–415.

  47. 47.

    Hendrix RW, Lawrence JG, Hatfull GF, et al. The origins and ongoing evolution of viruses. Trends Microbiol. 2000;8(11):504–8.

  48. 48.

    Petrovski S, Tillett D, Seviour RJ. Genome sequences and characterization of the Gordonia phages GTE5 and GRU1 and their use as potential biocontrol agents. Appl Environ Microbiol. 2012;78:42–7.

  49. 49.

    Rohwer F, Edwards R. The phage proteomic tree: a genome-based taxonomy of phage. J Bacteriol. 2002;184:4529–5.

  50. 50.

    Marchler-Bauer A, Derbyshire MK, Gonzales NR, et al. CDD: NCBI's conserved domain database. Nucleic Acids Res. 2015;43(Database issue):222–6.

  51. 51.

    Van Dessel W, Van Mellaert L, Liesegang H, et al. Complete genomic nucleotide sequence and analysis of the temperate bacteriophage VWB. J Anne Virology. 2005;331:325–37.

  52. 52.

    Grande L, Michelacci V, Tozzoli R, et al. Whole genome sequence comparison of vtx2-converting phages from Enteroaggregative Haemorrhagic Escherichia coli strains. BMC Genomics. 2014. doi:10.1186/1471-2164-15-574.

  53. 53.

    Leon-Velarde CG, Kropinski AM, Chen S, et al. Complete genome sequence of bacteriophage vB_YenP_AP5 which infects Yersinia enterocolitica of serotype O:3. Virol J. 2014. doi:10.1186/1743-422X-11-188.

  54. 54.

    Schouler C, Ehrlich SD, Chopin MC. Sequence and organization of the lactococcal prolate-headed bIL67 phage genome. Microbiology. 1994;140:3061–9.

  55. 55.

    Ford ME, Sarkis GJ, Belanger AE, et al. Genome structure of mycobacteriophage D29: implications for phage evolution. J Mol Biol. 1998;279(1):143–64.

  56. 56.

    Kropinski AM, Kovalyova IV, Billington SJ, et al. The genome of epsilon15, a serotype-converting, Group E1 Salmonella enterica-specific bacteriophage. Virology. 2007;369(2):234–44.

  57. 57.

    Katsura I. Structure and inherent properties of the bacteriophage lambda head shell. VI. DNA-packaging-defective mutants in the major capsid protein. J Mol Biol. 1989;205(2):397–405.

  58. 58.

    Lurz R, Orlova EV, Günther D, et al. Structural organisation of the head-to-tail interface of a bacterial virus. J Mol Biol. 2001;310(5):1027–37.

  59. 59.

    Gual A, Camacho AG, Alonso JC. Functional analysis of the terminase large subunit, G2P, of Bacillus subtilis bacteriophage SPP1. J Biol Chem. 2000;275(45):35311–9.

  60. 60.

    Noirot-Gros MF, Velten M, Yoshimura M, et al. Functional dissection of YabA, a negative regulator of DNA replication initiation in Bacillus subtilis. Proc Natl Acad Sci U S A. 2006;103(7):2368–73.

  61. 61.

    Hayashi M, Ogura Y, Harry EJ, et al. Bacillus subtilis YabA is involved in determining the timing and synchrony of replication initiation. FEMS Microbiol Lett. 2005;247(1):73–9.

  62. 62.

    Luscombe NM, Austin SE, Berman HM, et al. An overview of the structures of protein-DNA complexes. Genome Biol 2000, doi:10.1186/gb-2000-1-1-reviews001.

  63. 63.

    Dutta R, Inouye M. GHKL, An emergent ATPase/kinase superfamily. Trends Biochem Sci. 2000;25(1):24–8.

  64. 64.

    Hvorup RN, Winnen B, Chang AB, et al. The multidrug/oligosaccharidyl-lipid/polysaccharide (MOP) exporter superfamily. Eur J Biochem. 2003;270(5):799–813.

  65. 65.

    Satoh A, Wang Y, Malsam J, et al. Golgin-84 is a rab1 binding partner involved in Golgi structure. Traffic. 2003;4(3):153–61.

  66. 66.

    Bruner SD, Norman DP, Verdine GL. Structural basis for recognition and repair of the endogenous mutagen 8-oxoguanine in DNA. Nature. 2000;403(6772):859–66.

  67. 67.

    Ibrahim H, Wilusz J, Wilusz CJ. RNA recognition by 3’-to-5’ exonucleases: the substrate perspective. Biochim Biophys Acta. 2008;1779(4):256–65.

  68. 68.

    Aravind L, Makarova KS, Koonin EV. SURVEY AND SUMMARY: holliday junction resolvases and related nucleases: identification of new families, phyletic distribution and evolutionary trajectories. Nucleic Acids Res. 2000;28(18):3417–32.

  69. 69.

    Ku WY, Liu YW, Hsu YC, et al. The zinc ion in the HNH motif of the endonuclease domain of colicin E7 is not required for DNA binding but is essential for DNA hydrolysis. Nucleic Acids Res. 2002;30(7):1670–8.

  70. 70.

    Cheng X, Zhang X, Pflugrath JW, et al. The structure of bacteriophage T7 lysozyme, a zinc amidase and an inhibitor of T7 RNA polymerase. Proc Natl Acad Sci U S A. 1994;91(9):4034–8.

  71. 71.

    Petrovski S, Seviour RJ, Tillett D. Genome characterization of the polyvalent lytic bacteriophage GTE2 with the potential for biocontrol of Gordonia, Rhodococcus and Nocardia stabilized foams in activated sludge plants. Appl Environ Microbiol. 2011;77(12):3923–9.

  72. 72.

    Petrovski S, Dyson ZA, Seviour RJ, et al. Small but sufficient: the Rhodococcus phage RRH1 has the smallest known Siphoviridae genome at 14.2 kb. J Virol. 2012;86:358–63.

  73. 73.

    Smith KC, Castro-Nallar E, Fisher JNB, et al. Phage cluster relationships identified through single gene analysis. BMC Genomics. 2013;14:410.

  74. 74.

    Sullivan MB, Coleman ML, Quinlivan V, et al. Portal protein diversity and phage ecology. Environ Microbiol. 2008;10(10):2810–23.

  75. 75.

    Hooton SP, Timms AR, Rowsell J, et al. Salmonella Typhimurium-specific bacteriophage FSH19 and the origins of species specificity in the Vi01-like phage family. Virol J. 2011;8:498.

  76. 76.

    Dziarski R, Gupta D. Protein family review The peptidoglycan recognition proteins (PGRPs). Genome Biol. 2006;7:232.

  77. 77.

    Aravind L, Anantharaman V, Balaji S, et al. The many faces of the helix-turn-helix domain: Transcription regulation and beyond. FEMS Microbiol Rev. 2005;29:231–62.

  78. 78.

    Wendt JL, Feiss M. A fragile lattice: replacing bacteriophage E’s head stability gene D with the shp gene of phage 21 generates the Mg2 + −dependent virus, E shp. Virology. 2004;326:41–6.

  79. 79.

    Clark AJ, Inwood W, Cloutier T, Dhillon TS. Nucleotide sequence of coliphage HK620 and the evolution of lambdoid phages. J Mol Biol. 2001;311:657–79.

  80. 80.

    Canchaya C, Fournous G, Chibani-Chennoufi S, Dillmann ML, Brussow H. Phage as agents of lateral gene transfer. Curr Opin Microbiol. 2003;6:417–24.

  81. 81.

    Recktenwald J, Schmidt H. The nucleotide sequence of Shiga toxin (Stx) 2e-encoding phage phiP27 is not related to other Stx phage genomes, but the modular genetic structure is conserved. Infect Immun. 2002;70:1896–908.

  82. 82.

    Dekel-Bird NP, Avrani S, Sabehi G, Pekarsky I, Marston MF, Kirzner S, Lindell D. Diversity and evolutionary relationships of T7-like podoviruses infecting marine cyanobacteria. Environ Microbiol. 2013;15:1476–91.

  83. 83.

    Lawrence JG, Hatfull GF, Hendrix RW. Imbroglios of viral taxonomy: genetic exchange and failings of phenetic approaches. J Bacteriol. 2002;184:4891–905.

  84. 84.

    Lima-Mendez G, Van Helden J, Toussaint A, Leplae R. Reticulate representation of evolutionary and functional relationships between phage genomes. Mol Biol Evol. 2008;25:762–77.

Download references

Acknowledgement

The authors would like to thank the Center for Research in Agricultural Genomics (CRAG) service laboratory, Barcelona, Spain for providing the sequencing instruments and reagents used in the study, Ezio Fontana, IBBR-CNR, Palermo, Italy for his advice and discussion about the whole genome data analysis and Heather Esson, Biology Center ASCR, Institute of Parasitology, Czech Republic for assistance with language editing.

Funding

This work does not obtained any fund.

Availability of data and materials

The complete genome sequences of Sf1 and Sf3 phages were deposited in the National Center for Biotechnology Information (NCBI) GenBank under accession numbers (KT221033 and KT221034), respectively.

Authors’ contributions

All authors conceived and designed the experiments; AS carried out the experiments and performed the bioinformatics and statistical analysis; AS and FM compiled the results and drafted the manuscript. All the authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

This article does not contain any studies with human participants or animals performed by any of the authors. Informed consent was obtained from all individual participants included in the study.

Author information

Correspondence to A. Sharaf.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Sharaf, A., Mercati, F., Elmaghraby, I. et al. Functional and comparative genome analysis of novel virulent actinophages belonging to Streptomyces flavovirens . BMC Microbiol 17, 51 (2017) doi:10.1186/s12866-017-0940-7

Download citation

Keywords

  • Bacteriophage
  • Biological stability
  • Whole genome sequence
  • NGS
  • Comparative genomics