Next Article in Journal
Attenuated Influenza Virions Expressing the SARS-CoV-2 Receptor-Binding Domain Induce Neutralizing Antibodies in Mice
Next Article in Special Issue
VIRIDIC—A Novel Tool to Calculate the Intergenomic Similarities of Prokaryote-Infecting Viruses
Previous Article in Journal
Screening for the Proteins That Can Interact with Grouper Nervous Necrosis Virus Capsid Protein
Previous Article in Special Issue
The Phage-Encoded N-Acetyltransferase Rac Mediates Inactivation of Pseudomonas aeruginosa Transcription by Cleavage of the RNA Polymerase Alpha Subunit
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 12
Issue 9
10.3390/v12090986
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploring the Remarkable Diversity of Culturable Escherichia coli Phages in the Danish Wastewater Environment

by
Nikoline S. Olsen
1,2,
Laura Forero-Junco
2,
Witold Kot
2,* and
Lars H. Hansen
2,*
1
Department of Environmental Science, Aarhus University, Frederiksborgvej 399, 4000 Roskilde, Denmark
2
Department of Plant and Environmental Sciences, University of Copenhagen, Thorvaldsensvej 40, 1871 Frederiksberg C, Denmark
*
Authors to whom correspondence should be addressed.
Viruses 2020, 12(9), 986; https://0-doi-org.brum.beds.ac.uk/10.3390/v12090986
Submission received: 17 July 2020 / Revised: 24 August 2020 / Accepted: 31 August 2020 / Published: 4 September 2020
(This article belongs to the Special Issue Viruses of Microbes 2020: The Latest Conquests on Viruses of Microbes)
Browse Figures
Versions Notes

Abstract

:
Phages drive bacterial diversity, profoundly influencing microbial communities, from microbiomes to the drivers of global biogeochemical cycling. Aiming to broaden our understanding of Escherichia coli (MG1655, K-12) phages, we screened 188 Danish wastewater samples and isolated 136 phages. Ninety-two of these have genomic sequences with less than 95% similarity to known phages, while most map to existing genera several represent novel lineages. The isolated phages are highly diverse, estimated to represent roughly one-third of the true diversity of culturable virulent dsDNA Escherichia phages in Danish wastewater, yet almost half (40%) are not represented in metagenomic databases, emphasising the importance of isolating phages to uncover diversity. Seven viral families, Myoviridae, Siphoviridae, Podoviridae, Drexlerviridae, Chaseviridae, Autographviridae, and Microviridae, are represented in the dataset. Their genomes vary drastically in length from 5.3 kb to 170.8 kb, with a guanine and cytosine (GC) content ranging from 35.3% to 60.0%. Hence, even for a model host bacterium, substantial diversity remains to be uncovered. These results expand and underline the range of coliphage diversity and demonstrate how far we are from fully disclosing phage diversity and ecology.

Graphical Abstract

1. Introduction

Phages are important ecological contributors, renewing organic matter supplies in nutrient cycles and driving bacterial diversity by enabling co-existence of competing bacteria by “Killing the winner” and by serving as genomic reservoirs and transport units [1,2]. Phage genomes are known to contain auxiliary metabolism genes (AMGs), toxins, and virulence factors [3,4,5,6,7]. Through lysogeny and transduction, they can transfer metabolic traits including antibiotic resistance to their hosts and can confer immunity against homologous phages [1].
Still, despite their ecological role, potential as antimicrobials and the fact that they carry a multitude of unknown genes with great potential for biotechnological applications, phages are vastly understudied. Around 10,000 phage genomes have now been published [8]. Though the number increases rapidly, we may have merely scratched the surface of the expected diversity. It is estimated that at least one billion bacterial species exist [9]. Hence, only phages targeting a tiny fraction of potential hosts have been reported. Efforts to disclose the range and diversity of phages targeting a single host have revealed a stunning display of diversity. The most scrutinized phage host is the Mycobacterium smegmatis, for which the Science Education Alliance Phage Hunters program has isolated more than 4700 phages and fully sequenced 680. These represent 30 distinct phage clusters [10,11]. This endeavour has provided a unique insight into viral and host diversity, evolution, and genetics [12,13,14,15]. No other phage host has been equally targeted, but numerous Escherichia coli phages (coliphages) have been isolated. The International Committee on Taxonomy of Viruses (ICTV) currently recognises 276 phage species originally isolated from Escherichia species [16], while 733 Escherichia phage assemblies are listed in the European Nucleotide Archive (ENA) [17].
As phages are expected to have evolutionary potential to migrate across microbial populations, host species may not be an ideal indicator of relatedness, but it serves as an excellent starting point to explore phage diversity. The hierarchical classification of phages is complicated by the high degree of horizontal gene transfer. Consequently, several classification systems have been proposed [18,19,20]. We may not yet have reached a point where it is reasonable to establish the criteria for a universal system [20]. Nonetheless, a system enabling a mutual understanding and exchange of knowledge is needed. Accordingly, we have classified our phages as per the ICTV guidelines [21].
Here, we aim to expand our understanding of coliphage diversity, by screening for coliphages targeting a single strain of E. coli. Earlier studies on coliphage diversity have explored other aspects. Grose & Casjens (2014) did a comprehensive in silico study on 337 genomes of tailed phages infecting Enterobacteriaceae to characterise the known diversity [8]. Jurczak-Kurek et al. (2016) isolated 60 coliphages from sewage samples on a single E. coli strain and thoroughly assessed the physiological diversity but did not sequence any of the coliphages [22]. Korf et al. (2019) used 29 individual E. coli strains from various sources to isolate, characterise, and sequence 50 diverse tailed coliphages including representatives of novel phage lineages, verifying that there is still something to discover [23]. Mathieu et. al., (2020) explored the presence of virulent and temperate coliphages in 648 faecal samples from 1-year-old children, revealing interesting compositional trends likely to impact gut microbiota dynamics [24].
In this study we applied the High Throughput Screening (HiTS) screening method on 188 wastewater samples using the K-12 MG1655 strain as host. This approach favours virulent, culturable dsDNA phages with a large burst size and a short latency period [25]. We hypothesised that using a non-pathogenic lab-strain like MG1655 would provide a broad diversity of phages as opposed to wildtype strains which are likely better equipped for avoiding a wide array of successful phage infections. Hence, we expected the screening to yield coliphages that were easy to work with but also distinct enough to expand the number of known coliphages.

2. Materials and Methods

The screening for coliphages was performed with the HiTS method as described in [25], though instead of direct plaque sequencing [26], lysates of wells giving rise to plaques were sequenced. In short, an overnight enrichment (37 °C) was performed in microplates with E. coli, media, and wastewater (0.5 mL/well); the next day, the enrichments were filtrated (0.45 µm), re-inoculated (∼1 μL), and re-incubated overnight (37 °C); on the third day, a second filtration (0.45 µm) and a spot-test (soft-agar overlay) were performed to indicate positive wells.

2.1. Sample Bacteria and Media

Inlet wastewater samples (188) were collected (40–50 mL) at two- to four time-points during July and August 2017 from 48 Danish wastewater treatment plants (WWTPs) geographically distributed in Denmark. The samples were centrifuged (9000× g, 4 °C, 10 min) and the supernatant filtered (0.45 μm) before storage in aliquots (−20 °C) until screening. The host bacterium was E. coli (MG1655, K-12), and the medium, Lysogeny Broth (LB), was amended with CaCl2 and MgCl2 (final concentration 10 mM).

2.2. Sequencing and Genome Characterisation

DNA extractions, clean-up (ZR-96 Clean and Concentrator kit, Zymo Research, Irvine, CA, USA), and sequencing libraries (Nextera® XT DNA kit, Illumina, San Diego, CA, USA) were performed according to manufacturer’s protocols with minor modifications as described in Kot et al. (2014) [26]. The libraries were sequenced as paired-end reads on an Illumina NextSeq platform with the Mid Output Kit v2 (300 cycles). The obtained reads were trimmed and assembled in CLC Genomics Workbench version 10.1.1. (CLC BIO, Aarhus, Denmark). Overlapping reads were merged with the following settings: mismatch cost: 2, minimum score: 15, gap cost: 3 and maximum unaligned end mismatches: 0, and then assembled de novo. Additional assemblies were constructed using SPAdes version 3.12.0 [27]. Gene prediction and annotation were performed using a customized RASTtk version 2.0 [28] workflow with GeneMark [29], with manual curation and verification using BLASTP [30], HHpred [31], and Pfam version 32.0 [32], or de novo annotation was performed using VIGA version 0.11.0 [33] based on DIAMOND searches (RefSeq Viral protein database) and HMMer searches (pVOG HMM database). All genomes were assessed for antibiotic resistance genes, bacterial virulence genes, restriction-modification genes, and auxiliary metabolism genes (AMGs) using ResFinder 3.1 [34,35], VirulenceFinder 2.0 [36], Restriction-ModificationFinder 1.1 (REBASE) [37], and VIBRANT version 1.0.1 [38], respectively.

2.3. Bioinformatics

Nucleotide (NT) and amino acid (AA) similarities were calculated using tools recommended by the ICTV [21], i.e., BLAST [30] for identification of the closest relative (BLASTn when possible, discontinuous megaBLAST (word size 16) for larger genomes), and Gegenees version 2.2.1 [39] for assessing phylogenetic NT (BLASTn) and AA (tBLASTx) distances of multiple genomes, with fragment size 200 bp and step size 100 bp. Intergenomic nucleotide sequence similarity and aligned genome fractions between all isolated phage species were plotted with VIRIDIC [40]. NT similarity was determined as percentage query cover multiplied by percentage identity. Novel phages were categorised according to ICTV taxonomy. The criterion of 95% DNA sequence similarity for demarcation of species was applied to identify novel species representatives and to determine uniqueness within the dataset. Evolutionary analyses for phylogenomic trees were conducted in MEGA7 version 2.1 (default settings) [41]. These were based on the large terminase subunit (Caudovirales), a gene commonly applied for phylogenetic analysis [42,43] and on the DNA replication gene (gpA) (Microviridae). The NT sequences were aligned by MUSCLE [44] and the evolutionary history inferred by the Maximum Likelihood method based on the Tamura-Nei model [45]. The trees with the highest log-likelihood and are shown. Pairwise whole genome comparisons were performed with Easyfig 2.2.2 [46] (BLASTn), curated by adding color-codes and identifiers in Inkscape version 0.92.2. The R package iNEXT [47,48] in R studio version 1.1.456 [49] was used for rarefaction, species diversity (q = 0, datatype: incidence_raw), extrapolation thereof (estimadeD), and estimation of sample coverage. The visualisation of genome sizes and GC contents was prepared in Excel version 16.31. Blast+ 2.9.0 [50] was used to perform a NT search of the coliphages (queries) against a database with the IMG/VR v2.0 database sequences [51] and the human gut virome database (GVD) v 1.7 [52]. Reads from metagenomes and metaviromes were mapped using bbmap 38.22 [53]. Genome breadth and depth coverage was calculated using genomecov from BEDtools 2.28.0 [54] and BamM 1.7.3 [55], respectively.

3. Results and Discussion

3.1. Wastewater Coliphages Are Remarkably Diverse

The sequenced coliphages were analysed strictly in silico, focusing on their relatedness to known phages, taxonomy, and distinctive characteristics. The genome assemblies had a coverage of ×20-12122 with an average of ×390.5 (Table S1). The genome screening algorithms identified no homologs of known virulence or antibiotic resistance genes. Though not a definitive exclusion, this is interpreted as a reduced risk of presence, a preferable trait for phage therapy. The majority of genes identified when screening for AMGs code for phage DNA modification pathways (Table S2).
The isolation method (HiTS) favours easily culturable plaque-forming virulent phages [25]. Still, even though we screened wastewater samples, which is a commonly used source for isolation of coliphages, we identified 136 coliphages of which 92 differed by ≥5% from published phage genomes and some with nucleotide (NT) similarities as low as 29% (Table 1). Based on Blastn analyses and the 95% nucleotide similarity demarcation, 104 of the coliphages are unique phage species (Table 1, Figure 1). Based on DNA homology and phylogeny, the 104 unique coliphages group into 14 distinct clusters and 7 single phages (Figure 1 and Figure S1). Coliphages were identified in samples from 44 of the 48 WWTPs (Table S1). There was no substantial difference in phage diversity distribution between samples of urban or rural origin (Figure S2). Samples without coliphages likely reflect the crude nature of the screening method and in some cases sequence or assembly issues and not actual absence. From the majority of positive samples (n = 58) a single phage was sequenced, though some lysates held more than one phage (28 lysates: 2 phages, 6 lysates: 3 phages, 1 lysate: 4 phages, Table S1).
The 95% nucleotide identity demarcation of species is an arbitrary delimitation. It does not consider the biological importance of the non-identical sequence parts and imposes a discrepancy between the demarcation of species depending on genome size. However, it provides a means to quantify and compare relatedness enabling estimations of, e.g., culturable virulent coliphage species richness in the Danish wastewater environment (Figure 2. An extrapolation of species richness (q = 0) predicts a total of 311 distinct species (requiring a sample size of ∼900 phages) (Table S3). The relatively small sample size in this study (n = 136) may subject the estimation to a large prediction bias. The sampling-method also introduces bias by selecting for abundance, latency, and burst size, thereby potentially underestimating diversity. Sequencing and assembly methods as well as the choice of a host further reduce the number of detected phage genomes. Nonetheless, the results indicate the minimal diversity of culturable virulent dsDNA coliphages (MG1655, K-12) in Danish wastewater, estimated to be as a minimum in the range of 183 to 350 unique phage species (Figure 2, Table S3).
The diversity of tailed dsDNA coliphages is well documented [8,22,23] and it is to be expected that a screening of nearly 200 wastewater samples would yield hitherto unknown phages. However, considering the use of only a single host strain and a crude isolation method ensuring that only a single or the few most successful phage(s) from each sample were sequenced, then the degree of novelty and diversity revealed is remarkable and verifies our hypothesis, as well as the efficiency of the HiTS method for exploring diverse phages of a single host [25].

3.2. Taxonomy of the 104 Novel Coliphages

Based on the confirmed morphology of closely related phages, six different families of the order Caudovirales are represented; Myoviridae (57.4%), Drexlerviridae former Siphoviridae (18.4%), Autographviridae former Podoviridae (7.4%), Siphoviridae (6.6%), Podoviridae (3.7%) and Chaseviridae (0.7%), as well as Microviridae (5.9%), order Petitvirales (Figure 1). A similar distribution of coliphages from surface water, manure, sewage, and animal faeces was found by Korf et al. (2019); 70% Myoviridae, 22% Siphoviridae, and 8% Podoviridae [23]. Grose & Casjens (2014) also identified more clusters belonging to the Myoviridae, than the Siphoviridae and the fewest for the Podoviridae when analysing genomes of Caudovirales infecting Enterobacteriaceae [8]. Jurczak-Kurek et al. (2016) found more Siphoviridae than Myoviridae, but also found the Podoviridae to be the least abundant coliphages in sewage [22]. However, these distributions likely reflect abundance distributions of culturable phages and not necessarily natural abundances.

3.2.1. Fifty-Five Myoviridae Species

The 55 Myoviridae phage species represent the greatest span in genome sizes, from the Suspvirus mistaenkt (86.7 kb) to the Dhakavirus dhaeg (170.8 kb) (Figure 2B), and all, except the Krischvirus, code for tRNAs (Table 1). The Myoviridae group into eight distinct clusters and one single phage (mistaenkt), representing three subfamilies; Tevenvirinae, Vequintavirinae and Ounavirinae, in addition to the Phapecoctaviruses and a cluster of six unclassified Myoviridae (Figure 1, Table 1). The Tevenvirinae represent four genera, two krischviruses, five tequatrovirus, two dhakaviruses and two mosigviruses notable for their ability to perform arabinosylation of hmC [56].
The isolated Vequintavirinae are all vequintaviruses closely related (91.1–93.8%, BLAST) to classified species and were identified in samples from 12 of the 48 WWTPs. All of the Ounavirinae but the Suspvirus mistaenkt are felixounaviruses (89.7–93.9%, BLAST). The Felixounavirus is a relatively large genus with 16 recognized species isolated from Escherichia and Salmonella. In this study, felixounaviruses were identified 33 times in samples from no less than 23 WWTPs, indicating that they are ubiquitous in the Danish wastewater environment and that they are easily cultivated, though the method prevents assessment of relative preponderance. Felixounaviruses often have broad within-genus host ranges, and isolates have been shown to be able to rapidly expand their host range when challenged, co-coinciding with mutations in the long tail gene [57,58]. Five of the Myoviridae are members of the newly announced genus Phapecoctavirus with substantial similarity (86–90%, BLAST) to the type species Escherichia phage phAPEC8 (JX561091) [23,59]. The five phages in the last of the Myoviridae clusters are an even more homogeneous group than the phapecoctaviruses (Figure 3 and Figure S1). All five are closely related (92–95%, BLAST) to the same five unclassified Enterobacteriaceae phages vB_Ecom_PHB05 (MF805809), vB_vPM_PD06 (MH816848), ECGD1 (KU522583), phi92 (NC_023693), and vB_vPM_PD114 (MH675927) [60,61]; this group of nine phages is distinct (<44% NT similarity, BLAST) from all other described phages and thus represents a yet to be classified genus, presumably with the first sequenced phage phi92 as type species. Phi92 was isolated in 1982 and has been thoroughly characterised; it has a broad across-genus (Salmonella, Escherichia) host range enabled by multiple divergent tail fibres and can infect both non-capsulated and encapsulated hosts as it has a unique endosialidase tailspike encoded by gene 143 [60,61]. Interestingly, this gene appears to be unique for phi92, though other versions of a putative tailspike are present at the same position in the genomes of alia, PHB05, ECGD1, and the two PD06 and PD114 phages (Figure 3). Both the phapecoctaviruses and the unclassified Myoviridae genomes code for a complete dTDP-rhamnose biosynthesis pathway. The presence of a dTDP-rhamnose biosynthesis pathway in the DNA metabolism region of phage genomes is peculiar; one possible explanation is that these phages utilize rhamnose for glycosylation of hydroxy-methylated NTs in the same manner as the T4-generated glucosyl-hmC [56].

3.2.2. A New Addition to the Small Family Chaseviridae

The distinctive phage flopper only shares NT similarity (38.5–87%, BLAST) with ten other phages; it belongs to the newly established Carltongylesvirus (80.8–87% NT similarity, BLAST) of the new family Chaseviridae. This family currently has only nine species and the Carltongylesvirus only two species, Escherichia phage phiEcoM_GJ1 (EF460875) and Escherichia phage ST32 (MF044458). Both type species GJ1 and ST32 have broad within-genus host ranges [62,63]. NT similarity between flopper and GJ1 is partially low in the gene for the putative tail tape measure and also low between all three phage genomes in a tail fiber gene (Figure 4). The carltongylesviruses are unique in having characteristic Myoviridae morphology, i.e., icosahedral head, neck, and a contractile tail with tail fibres and also code for RNA polymerases, a feature otherwise characteristic to the T7-like phages of the Autographiviridae family [62,64,65].

3.2.3. Six Microviridae of Two Genera

The single phage Lilleven and the five gequatroviruses belong to the subfamily Bullavirinae, family Microviridae, order Petitvirales, characterised by ssDNA non-enveloped icosahedral phages (Table 1). Lilleven is a novel species of the genus Alphatrevirus, closely related to (93.9% NT similarity, BLAST, 89–90% AA similarity, Gegenees) the Alphatrevirus Enterobacteria phage St1 (NC_012868) (Figure S3). The five gequatroviruses only differ from one another by single NT polymorphisms and in non-coding regions (Figure S3). They cluster and share genomic organisation and extensive NT similarity (92.6–94.5%, BLAST) with the unclassified Microviridae Escherichia phage SECphi17 (LT960607), but only have 59.1–67.9% NT similarity (BLAST) with recognised Gequatrovirus species, with which they have almost no sequence similarity in the region coding for the major spike protein (gpG), a distinctive marker of the subfamily Bullavirinae involved in host attachment (Figure S3, Table 1) [66]. However, considering the pronounced gene synteny between their relatively small genomes and a conserved AA similarity (62–64%, Gegenees), they are considered gequatroviruses.
The sequencing of the Microviridae is peculiar, as library preparation with the Nextera® XT DNA kit applies transposons targeting dsDNA. However, during Microviridae infection, the host polymerase converts the viral ssDNA into an intermediate state of covalently closed dsDNA, which is then replicated in a rolling circle by viral replication proteins transcribed by the host RNA polymerase [67]. This intermediate state may have enabled the library preparation. The presence of host DNA (2.8–39.1% of reads) in the sequence results of these samples indicates an insufficient initial DNase I treatment (Table S4), which can be attributed to chemical inhibition or inactivation of the enzyme by adhesion to the sides of wells. Hence, it is reasonable to assume that the extracted microvirus DNA was captured as free dsDNA inside host cells during ongoing infections.

3.2.4. Twenty Drexlerviridae Phages Including a New Linage Representative

The 20 species of the new family Drexlerviridae represent a considerable expansion of the new subfamily Tempevirinae [68]. Eight of the Drexlerviridae belong to the new genus Warwickvirus (five species) with Escherichia virus swan01 as type species (LT841308), as they have ≥84.9% NT similarity (BLASTS) to recognised species thereof. The other eight are of the genus Hanrivervirus (NT: 86–90%, BLAST and AA: 77–85% Gegenees, Figure S4), currently consisting of only the type species Shigella virus pSf-1 (NC_021331) isolated from the Han River in Korea [69]. The warwickviruses and hanriverviruses isolated in this study all have comparable genome sizes, GC contents, and gene organisation with the respective type species (Figure 5, Table 1). During their differentiation, many deletions and insertions of small hypothetical genes have occurred; most notable is a unique version of a putative tail-spike protein in seven of the new Hanrivervirus species and all of the new Warwickvirus species, indicating a variety of divergent host ranges (Figure 5). All the hanriverviruses code for (putative) dam, and Psf-1 is resistant against at least six restriction endonucleases [69], suggesting that these phages employ DNA methylation as a defence strategy.
The last Drexlerviridae is Jahat. The warwickviruses and hanriverviruses form a monophyletic clade together with Jahat (Figure 1 and Figure S1). Even though Jahat has its own branch, this phage has gene synteny, slightly higher but comparable GC content, and shares an equal degree of NT similarity ≤68.7% with phages of both the Hanrivervirus and Warwickvirus (Figure 5). Hence, Jahat cannot with confidence be assigned to either genus but falls in between, barely different enough to represent its own genus—an indicator of the genetic continuum of phages challenging taxonomic delimitations.

3.2.5. Eight Siphoviridae Species and a Novel Genus Representative

The eight Siphoviridae species vary greatly in GC content, ranging from 44.6% (Skure) to 54.6% (welsh), but are quite similar in genome sizes, 49.7–54.6 kb (Figure 2). Five of these phages are of the genus Dhillonvirus as they have substantial NT similarity (77–80%, BLAST) and pronounced gene synteny with the type species Escherichia virus HK578. As with the hanriverviruses and warwickviruses, their genomes only differ in minor hypothetical genes and have limited NT similarity in a gene of highly variable length coding for a tail fiber (gp26 in HK578) (Figure 6), a phenomenon also observed in the dhillonviruses isolated by Korf et al. (2019), which correspondingly had divergent host ranges [23]. Each of the three remaining Siphoviridae represents a different genus. Based on NT similarity and the presence of the canonical 7-deazaguanine operon, Skure is of the 13-species genus Seuratvirus, while buks is assigned to the two-species genus Jerseyvirus, subfamily Guernseyvirinae.
Interestingly, the Siphoviridae Halfdan has only minuscule similarity with described phages (12–29%, BLAST). These entail two Pseudomonas phages vB_PaeS_SCUT-S3 (MK165657) and Ab26 (HG962376) [70], both Septimatreviruses, two Acinetobacter phages of the Lokivirus IMEAB3 (KF811200) and type species Acinetobacter virus Loki [71], and to a lesser degree the unclassified Achromobacter phage phiAxp-1 (KP313532) [72]. They have a common gene organization, yet their intra-Gegenees scores are low (≤1% BLASTn, <43% BLASTx, Figure S5), and NT similarity is negligent in roughly one-third of Halfdan’s 57 CDSs (Figure 7). The TerL-based phylogeny and AA similarity also indicate a distant relation, although grouping Halfdan closer with the lokiviruses (40–43%, Gegenees BLATSx) than the septimatreviruses (33–34%, Gegenees BLATSx) (Figure 7, Figure S5). Clearly Halfdan is distinct from all other described phages and hence the first phage sequenced of a new Siphoviridae genus.

3.2.6. Nine Autographiviridae Species

The nine Autographiviridae all have the hallmarks of this new family, i.e., unidirectionally encoded genes and RNA polymerases [65,68,73]. They belong to the genus Bonnellvirus, as they have conserved gene organisation, a similar GC content, and also share considerable NT similarity (69–93%, BLAST) with the type species Enterobacteria phage J8-65 (NC_025445) (Figure S6). The genomes of the nine new bonnellviruses and J8-65 are highly similar with differences primarily in small hypothetical genes, though Lidtsur codes for a unique version of tailspike colanidase (Figure S6). Lidtsur was deposited to the NCBI GenBank before the others and is currently the only one which is an ICTV-approved species representative.

3.2.7. Four Podoviridae Species Including Two Novel Genus Representatives

The four Podoviridae all have high (>59%) GC contents and represent no less than three distinct genera (Figure 2, Figure 8). Skarpretter is the type and only species of the genus Skarprettervirus [74]. Skarpretter is distinct from all described phages sharing only 38% NT similarity (BLAST) with the Giessenvirus Escherichia phage C130_2 (MH363708) isolated from cheese [75] (Figure 8 and Figure S7). Sortsne is the type species of the genus Sortsnevirus [74], currently consisting of only Sortsne and Klebsiella phage vB_KpnS_IME279 (MF614100); however, based on high NT similarity and conserved gene organization with IME279 (89.8%, BLAST), we suggest that sortkaff also belongs in Sortsnevirus (Figure 8 and Figure S7). The last Podoviridae sortsyn is of the new 2-species genus Murrayvius [76], as it shares a high degree of NT similarity and has conserved gene organization with the type species Enterobacteria phage IME_EC2 (KF591601) isolated from hospital sewage [77] (Figure 8).

3.3. The Wastewater Coliphages Are Largely Absent in Metaviromes

In order to investigate the prevalence of the 104 coliphage species in different environments we mapped the reads of 510 metagenomes from studies of primarily Danish wastewater, pig, and human gut samples (Table S5) [78,79]. The threshold for significant hits was set as mapped reads covering ≥70% of a coliphage genome, and the distribution of the mapped reads was assessed to verify that this threshold ensured identification of closely related phages (Figure S8). No hits were found for any of the coliphages. This is likely a consequence of sequencing depth and sample preparation, as prior to sequencing, these metagenome samples were concentrated by centrifugation as a pellet or by CsCl gradient and the supernatant was either discarded or stored for future studies, and as a result, a large proportion of potential phage reads was omitted. Subsequently, we also searched for the coliphages in hundreds of metavirome datasets (Table S5) from Irish and Chinese faecal, human, animal, and water samples using the same read mapping method (Table S6). There were no hits to the human faecal contigs from Ireland [80], while 22 of the 104 coliphage genomes (21%) representing ten genera were covered by >70% by reads from 10 (mammals and birds) of the 38 (26%) Chinese Wang study libraries (Figure 9) [80]. For most phage genera, only reads from a single sample matched, though reads from five metaviromes (pet dog, pig, yak, and flamingo faeces) matched (>70% read coverage) with the Alphatrevirus Lilleven, and reads from seven metaviromes (dog, red panda, giant panda, non-human primate, masked civet, pig, and chicken faeces) matched (>70% read coverage) with the Carltongylesvirus flopper (Table S6). Finally, the genome sequences of the 104 coliphages (queries) were blasted against a database of 735,106 uncultured viral genomes (UVIGs) from the Integrated Microbial Genomes/Virus (IMG/VR) database, derived from a wide range of sample types including marine, freshwater, terrestrial, and hosts [51], as well as 13,203 UVIGs from human gut samples retrieved from the human gut virome (GVD) database (Table S7) [52]. The coliphage genomes were also blasted against the 8392 isolated virus genomes (iVGs) of the IMG/VR database and based on the observed alignment coverage distribution (Figure S9), significant matches were defined as those covering >80% of coliphage genomes. With this threshold, there were significant matches for 23 of the 104 (22%) coliphage genomes to four of the 735,106 (0.0005%) IMG/VR UVIG sequences (Figure 9, Table S7).
Only phages from 14 of the 24 taxonomic groups of coliphages from this study had matches in the virome databases assessed. For only 14 of the 62 coliphages with matches, a closely related phage could be identified in more than one virome. Although the coliphages are omnipresent and culturable in Danish wastewater, they are for a large part not represented in metagenomic data, and therefore these coliphage genomes provide valuable information. A lack of representation in metagenomic data could be caused by low natural abundance as this would result in insufficient sequencing depth for genome assembly within metagenomes/-viromes. The Siphoviridae Halfdan, the Myoviridae mistaenkt, and the phages of the Drexlerviridae genus Warwickvirus, the Autographviridae genus Bonnellvirus, all the Podoviridae genera, Sortsnevirus, Murrayvirus and Skarprettervirus, as well as the Microviridae genus Gequatrovirus did not match any UVIG sequences, nor did the reads from any virome cover ≥70% of their genomes. Surprisingly, the gequatroviruses were not detected in any of the gut-viromes, even though that apart from temperate and crAss-like phages, Microviridae dominates in human, mammal, and bird gut-microbiomes and that 860 Microviridae genomes were assembled from the assorted Wang et al. (2019) metaviromes [80,81]. However, due to the relatively small size of the Microviridae genomes, substitution of a single gene is enough to warrant a ∼5–33% difference in NT similarity, putting them below the set threshold for identification. The fact that both Halfdan and the bonnellviruses and all the Podoviridae represent novel genera with very few close relatives suggests that these lineages are under-sampled and not sufficiently abundant in the environments explored by metagenomic sequencing and assessed in this study. These findings underline the importance of isolating and sequencing individual phages in order to uncover diversity. It is plausible that phages selected for by plating techniques are not those that are naturally abundant; however, this cannot be concluded based on these results. Future studies should compare the diversity obtained by isolating to metagenomic sequencing of metaviromes of identical samples in order to establish the degree of discrepancy between these two methods.

4. Conclusions

By screening 188 wastewater samples, we identified 104 coliphages species (MG1655–K12), enabling us to predict the species richness of culturable virulent dsDNA coliphages in Danish wastewater, which is predicted to be at least 183–350 and expected to fluctuate drastically over time. The true species richness is likely even higher as the isolation, DNA extraction, library construction, and genome assembly method as well as the choice of a host all are liable to reduce the number of phages detected. Ninety-two of the newly isolated coliphages represent novel species of seven families; Myoviridae, Siphoviridae, Podoviridae, Drexlerviridae, Chaseviridae, and Microviridae. Though most of them distribute into 18 established genera, the diversity of these many phages isolated from a single strain is notable. They vary greatly in genome size and have a broad GC content range.
Apart from the analyses applied, the main difference between this and the comprehensive Korf et al., study from 2019 [23] is the isolation approach. Korf et al., isolated 50 phages from various sample types over several years from a wide collection of clinical E. coli isolates, whereas the wastewater sample collection and phage isolation in this study were performed in a matter of weeks on a single strain of E. coli. Still, the distribution of phage types including many of the same genera and the discovery of a handful of phages with limited similarity to known phages are in many aspects comparable, suggesting that the method of isolation (plaque purification) is perhaps the key limiting factor for uncovering the diversity of coliphages. However, fewer than 60% of the 104 coliphages are represented in the assessed metaviromes, emphasising the importance of cultivating phages to uncover the true diversity.
These findings add to our understanding of phage ecology and diversity, and through classification of these many phages we come yet another step closer to a more refined taxonomic understanding of phages. Furthermore, the numerous and diverse phages isolated in this study, all lytic to the same single strain, serve as an excellent opportunity to learn important phage-host interactions in future studies. These include, but are not limited to, lysogen-induced phage immunity, host-range, and anti-RE systems.
Finally, the first genus representative for at least three novel genera was sequenced in this study. Skarprettervirus and Sortsnevirus of the Podoviridae have been accepted by the ICTV. We propose that Halfdan is the type species of a new Siphoviridae genus, that the four novel Myoviridae species muut, alia, outra, and inny together with five unclassified Myoviridae also represent a new genus, and as the Drexlerviridae Jahat cannot with confidence be assigned to any described genera, Jahat may also represent its own lineage. In conclusion, this study shows that uncharted territory remains for even well-studied phage hosts and that cultivation approaches uncover vital genomes that seem absent from metagenomic studies.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1999-4915/12/9/986/s1, Figure S1: VIRIDIC Plot of intergenomic sequence similarity, aligned genome fraction and genome length ratio for the 104 coliphage species, Figure S2: Distribution of the 136 coliphage species per WWTP location (urban vs rural), Figure S3: Phylogenetic and genomic analyses of the six Microviridae coliphages, Figure S4: Phylogenomic nucleotide distances (Gegenees, BLASTn) for the Hanrivervirus coliphages, Figure S5: Phylogenomic nucleotide distances (Gegenees, BLASTn) for phage Halfdan, Figure S6: Phylogenetic and genomic analyses of the Bonnellvirus coliphages, Figure S7: Phylogenomic nucleotide distances (Gegenees, BLASTn) for the Podoviridae coliphages, Figure S8: Mapped reads from the Wang et al. (2019) study to the genomes of flopper and Lilleven, Figure S9: Alignment coverage distribution of IMG/VR iVGs to the 104 coliphage genomes and significant hits to these, Table S1: List of phages per wastewater sample and corresponding sequence result metrics, Table S2: AMG screening results, Table S3: Species richness estimations, Table S4: Read mapping results for samples with Microviridae, Table S5: List of sequence data assessed for matches with the 104 coliphage genomes, Table S6: Mapped reads from the Wang et al. (2019) study and the IMG/VR microbiomes, Table S7: IMG/VR and GVD blast results.

Author Contributions

Conceptualization, N.S.O., W.K. and L.H.H.; Data curation, N.S.O. and L.F.-J.; Formal analysis, N.S.O. and L.F.-J.; Funding acquisition, W.K. and L.H.H.; Investigation, N.S.O.; Methodology, N.S.O., W.K. and L.H.H.; Project administration, N.S.O., W.K. and L.H.H.; Resources, N.S.O., W.K. and L.H.H.; Supervision, W.K. and L.H.H.; Validation, N.S.O. and W.K.; Visualization, N.S.O.; Writing–original draft, N.S.O.; Writing–review & editing, N.S.O., L.F.-J., W.K. and L.H.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Human Frontier Science Program Grant RGP0024/2018, Villum Experiment Grant 17595, and Aarhus University Research Foundation AUFF Grant E-2015-FLS-7-28.

Acknowledgments

The authors thank the many members of the Danish Water and Wastewater Association (DANVA) who kindly supplied us with time-series of wastewater samples from their WWTPs. A special thank you is extended to the Billund and Grindsted WWTP of Billund Vand & Energi, Lynetten, Avedøre and Damhusåen WWTP of BIOFOS, Kolding WWTP of BlueKolding, Esbjerg Øst, Esbjerg Vest, Ribe, Varde og Skovlund WWTPs of DIN Forsyning, Drøsbro, Hadsten, Hammel, Hinnerup and Voldum WWTPs of Favrskov Forsyning, Haerning WWTP of Herning Vand, Hillerød WWTP of Hillerød Forsyning, Lemvig WWTP of Lemvig Vand og Spildevand, Bevtoft, Gram, Haderslev, Halk, Jegerup, Nustrup, Over Jerstal, Skovby, Skrydstrup, Sommersted, Vojens and Årøsund WWTPs of Provas, Marselisborg, Egå and Viby WWTPs of Aarhus vand, Hedensted, Juelsminde and Tørring WWTPs of Hedensted Spildevand, Ejby Mølle, Bogense, Hofmansgave, Hårslev, Nordvest, Nordøst, Otterup and Søndersø WWTPs of VandCenterSyd, Øst and Vest WWTPs of Aalborg Forsyning and finally Helsingør WWTP of Forsyning Helsingør.

Conflicts of Interest

The authors declare no competing interests. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Weinbauer, M.G.; Rassoulzadegan, F. Are viruses driving microbial diversification and diversity? Environ. Microbiol. 2003, 6, 1–11. [Google Scholar] [CrossRef] [PubMed]
  2. Weitz, J.S.; Wilhelm, S.W. Ocean viruses and their effects on microbial communities and biogeochemical cycles. F1000 Biol. Rep. 2012, 4, 17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Crummett, L.T.; Puxty, R.J.; Weihe, C.; Marston, M.F.; Martiny, J.B. The genomic content and context of auxiliary metabolic genes in marine cyanomyoviruses. Virology 2016, 499, 219–229. [Google Scholar] [CrossRef] [PubMed]
  4. Boyd, E.F. Bacteriophage-Encoded Bacterial Virulence Factors and Phage–Pathogenicity Island Interactions. Adv. Virus Res. 2012, 82, 91–118. [Google Scholar] [CrossRef] [PubMed]
  5. Fortier, L.-C.; Sekulovic, O. Importance of prophages to evolution and virulence of bacterial pathogens. Virulence 2013, 4, 354–365. [Google Scholar] [CrossRef]
  6. Volkova, V.V.; Lu, Z.; Besser, T.; Gröhn, Y.T. Modeling the infection dynamics of bacteriophages in enteric Escherichia coli: Estimating the contribution of transduction to antimicrobial gene spread. Appl. Environ. Microbiol. 2014, 80, 4350–4362. [Google Scholar] [CrossRef] [Green Version]
  7. Bearson, B.L.; Allen, H.K.; Brunelle, B.W.; Lee, I.S.; Casjens, S.R.; Stanton, T.B. The agricultural antibiotic carbadox induces phage-mediated gene transfer in Salmonella. Front. Microbiol. 2014, 5, 52. [Google Scholar] [CrossRef] [Green Version]
  8. Grose, J.H.; Casjens, S.R. Understanding the enormous diversity of bacteriophages: The tailed phages that infect the bacterial family Enterobacteriaceae. Virology 2014, 468, 421–443. [Google Scholar] [CrossRef] [Green Version]
  9. Dykhuizen, D. Species Numbers in Bacteria. Proc. Calif. Acad. Sci. 2005, 56, 62. [Google Scholar]
  10. Hatfull, G.F.; Pedulla, M.L.; Jacobs-Sera, D.; Cichon, P.M.; Foley, A.; Ford, M.E.; Gonda, R.M.; Houtz, J.M.; Hryckowian, A.J.; Kelchner, V.A.; et al. Exploring the Mycobacteriophage Metaproteome: Phage Genomics as an Educational Platform. PLoS Genet. 2006, 2, e92. [Google Scholar] [CrossRef] [Green Version]
  11. Hatfull, G.F. Mycobacteriophages: Windows into Tuberculosis. PLoS Pathog. 2014, 10, e1003953. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Hatfull, G.F.; Jacobs-Sera, D.; Lawrence, J.G.; Pope, W.H.; Russell, D.A.; Ko, C.C.; Weber, R.J.; Patel, M.C.; Germane, K.L.; Edgar, R.H.; et al. Comparative Genomic Analysis of 60 Mycobacteriophage Genomes: Genome Clustering, Gene Acquisition, and Gene Size. J. Mol. Biol. 2010, 397, 119–143. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Dedrick, R.M.; Jacobs-Sera, D.; Bustamante, C.A.; Garlena, R.A.; Mavrich, T.N.; Pope, W.H.; Reyes, J.C.; Russell, D.A.; Adair, T.; Alvey, R.; et al. Prophage-mediated defence against viral attack and viral counter-defence. Nat. Microbiol. 2017, 2, 16251. [Google Scholar] [CrossRef] [PubMed]
  14. Jacobs-Sera, D.; Marinelli, L.J.; Bowman, C.; Broussard, G.W.; Bustamante, C.G.; Boyle, M.M.; Petrova, Z.O.; Dedrick, R.M.; Pope, W.H.; Advancing, S.E.; et al. On the nature of mycobacteriophage diversity and host preference. Virology 2012, 434, 187–201. [Google Scholar] [CrossRef]
  15. Hatfull, G.F. The Secret Lives of Mycobacteriophages. Adv. Virus Res. 2012, 82, 179–288. [Google Scholar] [CrossRef]
  16. Lefkowitz, E.J.; Dempsey, D.M.; Hendrickson, R.C.; Orton, R.J.; Siddell, S.G.; Smith, D.B. Virus taxonomy: The database of the International Committee on Taxonomy of Viruses (ICTV). Nucleic. Acids Res. 2018, 46, D708–D717. [Google Scholar] [CrossRef] [Green Version]
  17. European Nucleotide Archive ENA Browser. Available online: https://www.ebi.ac.uk/ena/browser/text-search?query=escherichiaphage (accessed on 29 May 2020).
  18. Lawrence, J.G.; Hatfull, G.F.; Hendrix, R.W. Imbroglios of viral taxonomy: Genetic exchange and failings of phenetic approaches. J. Bacteriol. 2002, 184, 4891–4905. [Google Scholar] [CrossRef] [Green Version]
  19. Aiewsakun, P.; Adriaenssens, E.M.; Lavigne, R.; Kropinski, A.M.; Simmonds, P. Evaluation of the genomic diversity of viruses infecting bacteria, archaea and eukaryotes using a common bioinformatic platform: Steps towards a unified taxonomy. J. Gen. Virol. 2018, 99, 1331–1343. [Google Scholar] [CrossRef]
  20. Nelson, D. Phage taxonomy: We agree to disagree. J. Bacteriol. 2004, 186, 7029–7031. [Google Scholar] [CrossRef] [Green Version]
  21. Adriaenssens, E.; Brister, J.R. How to Name and Classify Your Phage: An Informal Guide. Viruses 2017, 9, 70. [Google Scholar] [CrossRef] [Green Version]
  22. Jurczak-Kurek, A.; Gąsior, T.; Nejman-Faleńczyk, B.; Bloch, S.; Dydecka, A.; Topka, G.; Necel, A.; Jakubowska-Deredas, M.; Narajczyk, M.; Richert, M.; et al. Biodiversity of bacteriophages: Morphological and biological properties of a large group of phages isolated from urban sewage. Sci. Rep. 2016, 6, 34338. [Google Scholar] [CrossRef] [PubMed]
  23. Korf, I.H.; Meier-Kolthoff, J.P.; Adriaenssens, E.M.; Kropinski, A.M.; Nimtz, M.; Rohde, M.; Van Raaij, M.J.; Wittmann, J. Still Something to Discover: Novel Insights into Escherichia coli Phage Diversity and Taxonomy. Viruses 2019, 11, 454. [Google Scholar] [CrossRef] [Green Version]
  24. Mathieu, A.; Dion, M.; Deng, L.; Tremblay, D.; Moncaut, E.; Shah, S.A.; Stokholm, J.; Krogfelt, K.A.; Schjørring, S.; Bisgaard, H.; et al. Virulent coliphages in 1-year-old children fecal samples are fewer, but more infectious than temperate coliphages. Nat. Commun. 2020, 11, 1–12. [Google Scholar] [CrossRef] [PubMed]
  25. Olsen, N.S.; Hendriksen, N.B.; Hansen, L.H.; Kot, W. A New High-throughput Screening (HiTS) Method for Phages—Enabling Crude Isolation and Fast Identification of Diverse Phages with Therapeutic Potential. bioRxiv 2020. [Google Scholar] [CrossRef] [Green Version]
  26. Kot, W.; Vogensen, F.K.; Sørensen, S.J.; Hansen, L.H. DPS—A rapid method for genome sequencing of DNA-containing bacteriophages directly from a single plaque. J. Virol. Methods 2014, 196, 152–156. [Google Scholar] [CrossRef]
  27. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef] [Green Version]
  28. Brettin, T.; Davis, J.J.; Disz, T.; Edwards, R.A.; Gerdes, S.; Olsen, G.J.; Olson, R.; Overbeek, R.; Parrello, B.; Pusch, G.D.; et al. RASTtk: A modular and extensible implementation of the RAST algorithm for building custom annotation pipelines and annotating batches of genomes. Sci. Rep. 2015, 5, 8365. [Google Scholar] [CrossRef] [Green Version]
  29. Besemer, J.; Borodovsky, M. GeneMark: Web software for gene finding in prokaryotes, eukaryotes and viruses. Nucleic. Acids Res. 2005, 33, W451–W454. [Google Scholar] [CrossRef] [Green Version]
  30. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  31. Soding, J.; Biegert, A.; Lupas, A.N. The HHpred interactive server for protein homology detection and structure prediction. Nucleic. Acids Res. 2005, 33, W244–W248. [Google Scholar] [CrossRef] [Green Version]
  32. Finn, R.D.; Coggill, P.; Eberhardt, R.Y.; Eddy, S.R.; Mistry, J.; Mitchell, A.L.; Potter, S.C.; Punta, M.; Qureshi, M.; Sangrador-Vegas, A.; et al. The Pfam protein families database: Towards a more sustainable future. Nucleic. Acids Res. 2016, 44, D279–D285. [Google Scholar] [CrossRef] [PubMed]
  33. González-Tortuero, E.; Sutton, T.D.; Velayudhan, V.; Shkoporov, A.N.; Draper, L.A.; Stockdale, S.R.; Ross, R.P.; Hill, C. VIGA: A sensitive, precise and automatic de novo VIral Genome Annotator. bioRxiv 2018. [Google Scholar] [CrossRef] [Green Version]
  34. Zankari, E.; Hasman, H.; Cosentino, S.; Vestergaard, M.; Rasmussen, S.; Lund, O.; Aarestrup, F.M.; Larsen, M.V. Identification of acquired antimicrobial resistance genes. J. Antimicrob. Chemother. 2012, 67, 2640–2644. [Google Scholar] [CrossRef] [PubMed]
  35. Kleinheinz, K.A.; Joensen, K.G.; Larsen, M.V. Applying the ResFinder and VirulenceFinder web-services for easy identification of acquired antibiotic resistance and E. coli virulence genes in bacteriophage and prophage nucleotide sequences. Bacteriophage 2014, 4, e27943. [Google Scholar] [CrossRef] [Green Version]
  36. Joensen, K.G.; Scheutz, F.; Lund, O.; Hasman, H.; Kaas, R.S.; Nielsen, E.M.; Aarestrup, F.M. Real-Time Whole-Genome Sequencing for Routine Typing, Surveillance, and Outbreak Detection of Verotoxigenic Escherichia coli. J. Clin. Microbiol. 2014, 52, 1501–1510. [Google Scholar] [CrossRef] [Green Version]
  37. Roberts, R.J.; Vincze, T.; Posfai, J.; Macelis, D. REBASE-a database for DNA restriction and modification: Enzymes, genes and genomes. Nucleic. Acids Res. 2015, 43, D298–D299. [Google Scholar] [CrossRef]
  38. Kieft, K.; Zhou, Z.; Anantharaman, K. VIBRANT: Automated recovery, annotation and curation of microbial viruses, and evaluation of virome function from genomic sequences. bioRxiv 2019. [Google Scholar] [CrossRef] [Green Version]
  39. Ågren, J.; Sundström, A.; Håfström, T.; Segerman, B. Gegenees: Fragmented Alignment of Multiple Genomes for Determining Phylogenomic Distances and Genetic Signatures Unique for Specified Target Groups. PLoS ONE 2012, 7, e39107. [Google Scholar] [CrossRef]
  40. Moraru, C.; Varsani, A.; Kropinski, A.M. VIRIDIC-a Novel Tool to Calculate the Intergenomic 1 Similarities of Prokaryote-Infecting Viruses 2; Cold Spring Harbor Laboratory: Cold Spring Harbor, NY, USA, 2020. [Google Scholar]
  41. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [Green Version]
  42. Lopes, A.; Tavares, P.; Petit, M.A.; Guérois, R.; Zinn-Justin, S. Automated classification of tailed bacteriophages according to their neck organization. BMC Genom. 2014, 15, 1027. [Google Scholar] [CrossRef] [Green Version]
  43. Mercanti, D.J.; Rousseau, G.M.; Capra, M.L.; Quiberoni, A.; Tremblay, D.M.; Labrie, S.J.; Moineau, S. Genomic Diversity of Phages Infecting Probiotic Strains of Lactobacillus paracasei. Appl. Environ. Microbiol. 2016, 82, 95–105. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic. Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Tamura, K.; Nei, M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol. Biol. Evol. 1993, 10, 512–526. [Google Scholar] [CrossRef] [PubMed]
  46. Sullivan, M.J.; Petty, N.K.; Beatson, S.A. Easyfig: A genome comparison visualizer. Bioinformatics 2011, 27, 1009. [Google Scholar] [CrossRef] [PubMed]
  47. Hsieh, T.C.; Ma, K.H.; Chao, A. iNEXT: An R package for rarefaction and extrapolation of species diversity (Hill numbers). Methods Ecol. Evol. 2016, 7, 1451–1456. [Google Scholar] [CrossRef]
  48. Chao, A.; Gotelli, N.J.; Hsieh, T.C.; Sander, E.L.; Ma, K.H.; Colwell, R.K.; Ellison, A.M. Rarefaction and extrapolation with Hill numbers: A framework for sampling and estimation in species diversity studies. Ecol. Monogr. 2014, 84, 45–67. [Google Scholar] [CrossRef] [Green Version]
  49. RStudio Team. RStudio: Integrated Development for R; RStudio Inc.: Boston, MA, USA, 2016. [Google Scholar]
  50. Camacho, C.; Coulouris, G.; Avagyan, V.; Ma, N.; Papadopoulos, J.; Bealer, K.; Madden, T.L. BLAST+: Architecture and applications. BMC Bioinform. 2009, 10, 421. [Google Scholar] [CrossRef] [Green Version]
  51. Paez-Espino, D.; Chen, I.M.; Palaniappan, K.; Ratner, A.; Chu, K.; Szeto, E.; Pillay, M.; Huang, J.; Markowitz, V.M.; Nielsen, T.; et al. IMG/VR: A database of cultured and uncultured DNA viruses and retroviruses. Nucleic. Acids Res. 2017, 45, D457–D465. [Google Scholar] [CrossRef]
  52. Gregory, A.C.; Zablocki, O.; Howell, A.; Bolduc, B.; Sullivan, M.B. The human gut virome database. bioRxiv 2019. [Google Scholar] [CrossRef]
  53. Bushnell, B. BBMap Download|SourceForge.net. 2013. Available online: https://sourceforge.net/projects/bbmap/ (accessed on 16 June 2020).
  54. Quinlan, A.R.; Hall, I.M. BEDTools: A flexible suite of utilities for comparing genomic features. Bioinformatics 2010, 26, 841–842. [Google Scholar] [CrossRef] [Green Version]
  55. GitHub—Ecogenomics/BamM: Metagenomics-Focused BAM File Manipulation. 2020. Available online: https://github.com/ecogenomics/BamM (accessed on 24 August 2020).
  56. Thomas, J.A.; Orwenyo, J.; Wang, L.X.; Black, L.W. The odd “RB” phage—Identification of arabinosylation as a new epigenetic modification of DNA in T4-like phage RB69. Viruses 2018, 10, 313. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Kropinski, A.M.; Sulakvelidze, A.; Konczy, P.; Poppe, C. Salmonella phages and prophages—Genomics and practical aspects. Methods Mol. Biol. 2007, 394, 133–175. [Google Scholar]
  58. Rivera, D.; Hudson, L.K.; Denes, T.G.; Hamilton-West, C.; Pezoa, D.; Moreno-Switt, A.I. Two phages of the genera felixunavirus subjected to 12 h challenge on salmonella infantis showed distinct genotypic and phenotypic changes. Viruses 2019, 11, 586. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Tsonos, J.; Adriaenssens, E.M.; Klumpp, J.; Hernalsteens, J.P.; Lavigne, R.; De Greve, H. Complete genome sequence of the novel Escherichia coli phage phAPEC8. J. Virol. 2012, 86, 13117–13118. [Google Scholar] [CrossRef] [Green Version]
  60. Schwarzer, D.; Buettner, F.F.; Browning, C.; Nazarov, S.; Rabsch, W.; Bethe, A.; Oberbeck, A.; Bowman, V.D.; Stummeyer, K.; Mühlenhoff, M.; et al. A Multivalent Adsorption Apparatus Explains the Broad Host Range of Phage phi92: A Comprehensive Genomic and Structural Analysis. J. Virol. 2012, 86, 10384–10398. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Schwarzer, D.; Browning, C.; Stummeyer, K.; Oberbeck, A.; Mühlenhoff, M.; Gerardy-Schahn, R.; Leiman, P.G. Structure and biochemical characterization of bacteriophage phi92 endosialidase. Virology 2015, 477, 133–143. [Google Scholar] [CrossRef] [Green Version]
  62. Jamalludeen, N.; Kropinski, A.M.; Johnson, R.P.; Lingohr, E.; Harel, J.; Gyles, C.L. Complete genomic sequence of bacteriophage phiEcoM-GJ1, a novel phage that has myovirus morphology and a podovirus-like RNA polymerase. Appl. Environ. Microbiol. 2008, 74, 516–525. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Liu, H.; Geagea, H.; Rousseau, G.M.; Labrie, S.J.; Tremblay, D.M.; Liu, X.; Moineau, S. Characterization of the Escherichia coli Virulent Myophage ST32. Viruses 2018, 10, 616. [Google Scholar] [CrossRef] [Green Version]
  64. Kropinski, A.M.; Anany, H.; Adriaenssens, E.M. Create One New Family (Chaseviridae) Including Six New Genera in the Order Caudovirales 2019.047B. Available online: https://ictv.global/ICTV/proposals/2019.047B.zip (accessed on 17 May 2020).
  65. Lavigne, R.; Seto, D.; Mahadevan, P.; Ackermann, H.W.; Kropinski, A.M. Unifying classical and molecular taxonomic classification: Analysis of the Podoviridae using BLASTP-based tools. Res. Microbiol. 2008, 159, 406–414. [Google Scholar] [CrossRef]
  66. Labrie, S.J.; Dupuis, M.È.; Tremblay, D.M.; Plante, P.L.; Corbeil, J.; Moineau, S. A new Microviridae phage isolated from a failed biotechnological process driven by Escherichia coli. Appl. Environ. Microbiol. 2014, 80, 6992–7000. [Google Scholar] [CrossRef] [Green Version]
  67. Doore, S.M.; Fane, B.A. The microviridae: Diversity, assembly, and experimental evolution. Virology 2016, 491, 45–55. [Google Scholar] [CrossRef] [PubMed]
  68. Adriaenssens, E.M.; Sullivan, M.B.; Knezevic, P.; van Zyl, L.J.; Sarkar, B.L.; Dutilh, B.E.; Alfenas-Zerbini, P.; Łobocka, M.; Tong, Y.; Brister, J.R.; et al. Taxonomy of prokaryotic viruses: 2018-2019 update from the ICTV Bacterial and Archaeal Viruses Subcommittee. Arch. Virol. 2020, 165, 1253–1260. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Jun, J.W.; Kim, J.H.; Shin, S.P.; Han, J.E.; Chai, J.Y.; Park, S.C. Characterization and complete genome sequence of the Shigella bacteriophage pSf-1. Res. Microbiol. 2013, 164, 979–986. [Google Scholar] [CrossRef] [PubMed]
  70. Essoh, C.; Latino, L.; Midoux, C.; Blouin, Y.; Loukou, G.; Nguetta, S.P.; Lathro, S.; Cablanmian, A.; Kouassi, A.K.; Vergnaud, G.; et al. Investigation of a Large Collection of Pseudomonas aeruginosa Bacteriophages Collected from a Single Environmental Source in Abidjan, Côte d’Ivoire. PLoS ONE 2015, 10, e0130548. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Turner, D.; Wand, M.E.; Briers, Y.; Lavigne, R.; Sutton, J.M.; Reynolds, D.M. Characterisation and genome sequence of the lytic Acinetobacter baumannii bacteriophage vB_AbaS_Loki. PLoS ONE 2017, 12, e0172303. [Google Scholar] [CrossRef] [PubMed]
  72. Ma, Y.; Li, E.; Qi, Z.; Li, H.; Wei, X.; Lin, W.; Zhao, R.; Jiang, A.; Yang, H.; Yin, Z.; et al. Isolation and molecular characterisation of Achromobacter phage phiAxp-3, an N4-like bacteriophage. Sci. Rep. 2016, 6, 24776. [Google Scholar] [CrossRef] [Green Version]
  73. Turner, D.; Kropinski, A.M.; Alfernas-Zerbini, P.; Buttimer, C.; Lavigne, R.; Bister, J.R.; Tolstoy, I.; Morozova, V.V.; Babkin, I.V.; Kozlova, Y.N.; et al. Create one new family (Autographiviridae) including nine subfamilies and one hundred and thirty-two genera in the order Caudovirales. Microb. Biotechnol. 2019, 13, 1428–1445. [Google Scholar]
  74. Kropinski, A.M.; Adriaenssens, E.M. Create Four New Genera Including Six New Species in the Family Podoviridae 2019.076B. Available online: https://ictv.global/ICTV/proposals/2019.076B.zip (accessed on 17 May 2020).
  75. Sváb, D.; Falgenhauer, L.; Rohde, M.; Chakraborty, T.; Tóth, I. Complete genome sequence of C130_2, a novel myovirus infecting pathogenic Escherichia coli and Shigella strains. Arch. Virol. 2019, 164, 321–324. [Google Scholar] [CrossRef] [Green Version]
  76. Kropinski, A.M.; Adriaenssens, E.M. Create One New Genus (Murrayvirus) Including of Two New Species in the Family Siphoviridae 2019.069B. Available online: https://ictv.global/ICTV/proposals/2019.069B.zip (accessed on 17 May 2020).
  77. Hua, Y.; An, X.; Pei, G.; Li, S.; Wang, W.; Xu, X.; Fan, H.; Huang, Y.; Zhang, Z.; Mi, Z.; et al. Characterization of the morphology and genome of an Escherichia coli podovirus. Arch. Virol. 2014, 159, 3249–3256. [Google Scholar] [CrossRef]
  78. Hendriksen, R.S.; Munk, P.; Njage, P.; Van Bunnik, B.; McNally, L.; Lukjancenko, O.; Röder, T.; Nieuwenhuijse, D.; Pedersen, S.K.; Kjeldgaard, J.; et al. Global monitoring of antimicrobial resistance based on metagenomics analyses of urban sewage. Nat. Commun. 2019, 10, 1–12. [Google Scholar] [CrossRef]
  79. Xiao, L.; Estellé, J.; Kiilerich, P.; Ramayo-Caldas, Y.; Xia, Z.; Feng, Q.; Liang, S.; Pedersen, A.Ø.; Kjeldsen, N.J.; Liu, C.; et al. A reference gene catalogue of the pig gut microbiome. Nat. Microbiol. 2016, 1, 1–6. [Google Scholar] [CrossRef] [PubMed]
  80. Wang, H.; Ling, Y.; Shan, T.; Yang, S.; Xu, H.; Deng, X.; Delwart, E.; Zhang, W. Gut virome of mammals and birds reveals high genetic diversity of the family Microviridae. Virus Evol. 2019, 5, 1–8. [Google Scholar] [CrossRef] [PubMed]
  81. Shkoporov, A.N.; Clooney, A.G.; Sutton, T.D.; Ryan, F.J.; Daly, K.M.; Nolan, J.A.; McDonnell, S.A.; Khokhlova, E.V.; Draper, L.A.; Forde, A.; et al. The Human Gut Virome Is Highly Diverse, Stable, and Individual Specific. Cell Host. Microbe. 2019, 26, 527–541.e5. [Google Scholar] [CrossRef] [PubMed]
Viruses 12 00986 g001 550
Figure 1. Phylogenetic trees, bootstrap: 100, scalebar: substitutions per site. (A) Caudovirales based on large terminase subunit terL (Maximum log Likelihood: −1801.27). (B) Microviridae based on the DNA replication protein gene gpA (Maximum log Likelihood: −3922.55). (C) Taxonomic distribution of phages identified in 94 Danish wastewater samples, based on similarity to closest related and the ICTV Master Species list. 1 ≤95% similarity to other phages in the dataset. 2 ≤95% similarity to other phages in the dataset and the NCBI GenBank. 3 The Murrayvirus genus has by mistake been classified as Siphoviridae, but a proposed move to Podoviridae will be included in the 2021 ICTV ratification.
Figure 1. Phylogenetic trees, bootstrap: 100, scalebar: substitutions per site. (A) Caudovirales based on large terminase subunit terL (Maximum log Likelihood: −1801.27). (B) Microviridae based on the DNA replication protein gene gpA (Maximum log Likelihood: −3922.55). (C) Taxonomic distribution of phages identified in 94 Danish wastewater samples, based on similarity to closest related and the ICTV Master Species list. 1 ≤95% similarity to other phages in the dataset. 2 ≤95% similarity to other phages in the dataset and the NCBI GenBank. 3 The Murrayvirus genus has by mistake been classified as Siphoviridae, but a proposed move to Podoviridae will be included in the 2021 ICTV ratification.
Viruses 12 00986 g001
Viruses 12 00986 g002 550
Figure 2. (A) Number of phages and taxonomic family distribution per WWTP, as well as approximate geographical location thereof. (B) Bubble-diagram of the 104 unique coliphage species displaying genome size and GC content distribution; the area of the bubbles indicates the number of phages. (C) Sample completeness curve with confidence intervals (0.95). (D) Sample-size-based rarefaction and extrapolation curve with confidence intervals (0.95). Sampling units in figures C and D denote isolated phages (n = 136).
Figure 2. (A) Number of phages and taxonomic family distribution per WWTP, as well as approximate geographical location thereof. (B) Bubble-diagram of the 104 unique coliphage species displaying genome size and GC content distribution; the area of the bubbles indicates the number of phages. (C) Sample completeness curve with confidence intervals (0.95). (D) Sample-size-based rarefaction and extrapolation curve with confidence intervals (0.95). Sampling units in figures C and D denote isolated phages (n = 136).
Viruses 12 00986 g002
Viruses 12 00986 g003 550
Figure 3. Pairwise alignment of the unclassified Myoviridae phage species from this study (blue text) and their closest relatives (black text), the color bars between genomes indicate percent pairwise similarity (Easyfig, BlASTn). The genes marked with yellow code for a tailspike (gene 143 in phi92).
Figure 3. Pairwise alignment of the unclassified Myoviridae phage species from this study (blue text) and their closest relatives (black text), the color bars between genomes indicate percent pairwise similarity (Easyfig, BlASTn). The genes marked with yellow code for a tailspike (gene 143 in phi92).
Viruses 12 00986 g003
Viruses 12 00986 g004 550
Figure 4. Pairwise alignment of the new Carltongylesvirus phage flopper (blue) and the Carltongylesvirus phages GJ1 and ST32; the colour bars between genomes indicate percent pairwise similarity (Easyfig, BlASTn). Genomes have been modified to have similar starting points. Genes coding for the putative tail tape measure protein (red) and a tail fiber (orange) are colorised.
Figure 4. Pairwise alignment of the new Carltongylesvirus phage flopper (blue) and the Carltongylesvirus phages GJ1 and ST32; the colour bars between genomes indicate percent pairwise similarity (Easyfig, BlASTn). Genomes have been modified to have similar starting points. Genes coding for the putative tail tape measure protein (red) and a tail fiber (orange) are colorised.
Viruses 12 00986 g004
Viruses 12 00986 g005 550
Figure 5. Pairwise alignment of the new Hanrivervirus phages (green text) and type species pSf-1 (black text), the new Warwickvirus phages (blue text), and type species swan01 (black text) and Jahat (brown text); the color bars between genomes indicate percent pairwise similarity (Easyfig, BlASTn). Genomes have been modified to have similar starting points. Coloured genes (purple) code for a putative tail fiber.
Figure 5. Pairwise alignment of the new Hanrivervirus phages (green text) and type species pSf-1 (black text), the new Warwickvirus phages (blue text), and type species swan01 (black text) and Jahat (brown text); the color bars between genomes indicate percent pairwise similarity (Easyfig, BlASTn). Genomes have been modified to have similar starting points. Coloured genes (purple) code for a putative tail fiber.
Viruses 12 00986 g005
Viruses 12 00986 g006 550
Figure 6. Pairwise alignment of the new Dhillonvirus phage species (blue text) and the type species HK578 (black text); the colour bars between genomes indicate percent pairwise similarity (Easyfig, BlASTn). Genomes have been modified to have similar starting points. Highlighted genes (green) code for a putative tail fiber.
Figure 6. Pairwise alignment of the new Dhillonvirus phage species (blue text) and the type species HK578 (black text); the colour bars between genomes indicate percent pairwise similarity (Easyfig, BlASTn). Genomes have been modified to have similar starting points. Highlighted genes (green) code for a putative tail fiber.
Viruses 12 00986 g006
Viruses 12 00986 g007 550
Figure 7. Comparisons of the new phage lineage representative Halfdan (blue text) and closest relatives (black text). (A) Phylogenetic tree, Maximum log Likelihood: −7678.71, bootstrap 100, large terminase subunit TerL, scalebar: substitutions per site. (B) Pairwise alignment of phage genomes; color bars between genomes indicate percent pairwise similarity (Easyfig, BlASTn). Genomes have been modified to have similar starting points.
Figure 7. Comparisons of the new phage lineage representative Halfdan (blue text) and closest relatives (black text). (A) Phylogenetic tree, Maximum log Likelihood: −7678.71, bootstrap 100, large terminase subunit TerL, scalebar: substitutions per site. (B) Pairwise alignment of phage genomes; color bars between genomes indicate percent pairwise similarity (Easyfig, BlASTn). Genomes have been modified to have similar starting points.
Viruses 12 00986 g007
Viruses 12 00986 g008 550
Figure 8. Comparisons of the Podoviridae phage species from this study (blue text) with closest relatives (black text). (A) Phylogenetic tree, Maximum log Likelihood: −8023.43, bootstrap 100, large terminase subunit TerL, scalebar: substitutions per site. (B) Pairwise alignment of phage genomes; color bars between genomes indicate percent pairwise similarity (Easyfig, BlASTn). Genomes have been modified to have similar starting points.
Figure 8. Comparisons of the Podoviridae phage species from this study (blue text) with closest relatives (black text). (A) Phylogenetic tree, Maximum log Likelihood: −8023.43, bootstrap 100, large terminase subunit TerL, scalebar: substitutions per site. (B) Pairwise alignment of phage genomes; color bars between genomes indicate percent pairwise similarity (Easyfig, BlASTn). Genomes have been modified to have similar starting points.
Viruses 12 00986 g008
Viruses 12 00986 g009 550
Figure 9. Significant hits for finding close relatives of the phages isolated in this study in three databases. Hits are defined as hits ≥80% genome coverage when blasting the coliphage genomes against the IMG/VR and GVD databases, and as mapped reads covering ≥70% of individual coliphage genomes when mapping reads from the Wang study. The coliphages are grouped according to genera, and numbers in parentheses denote the number of coliphage species representing each genus; only genera with significant hits are shown. Color-codes of genera denote taxonomic family.
Figure 9. Significant hits for finding close relatives of the phages isolated in this study in three databases. Hits are defined as hits ≥80% genome coverage when blasting the coliphage genomes against the IMG/VR and GVD databases, and as mapped reads covering ≥70% of individual coliphage genomes when mapping reads from the Wang study. The coliphages are grouped according to genera, and numbers in parentheses denote the number of coliphage species representing each genus; only genera with significant hits are shown. Color-codes of genera denote taxonomic family.
Viruses 12 00986 g009
Table
Table 1. List of 104 unique Escherichia phage species identified in 94 Danish wastewater samples. Phages (n) denotes the number of phages isolated with more than 95% nucleotide sequence similarity. Similarity is sequence identity (%) times sequence coverage (%) to closest relative (Blastn). Taxonomy is based on similarity (BLASTn) to closest related.
Table 1. List of 104 unique Escherichia phage species identified in 94 Danish wastewater samples. Phages (n) denotes the number of phages isolated with more than 95% nucleotide sequence similarity. Similarity is sequence identity (%) times sequence coverage (%) to closest relative (Blastn). Taxonomy is based on similarity (BLASTn) to closest related.
PhageAccessionPhages (n)Genome (bp)ORFstRNAsGC (%)Family; GenusSimilarity (%)Closest RelativeAccession
tootikiMN850647188,2571282239Myoviridae; Felixounavirus90.2Escherichia phage vB_EcoM_Alf5NC_031082.1
mioMN850631183,4311211839.1Myoviridae; Felixounavirus89.7Salmonella virus VSe11MG251391.1
allfineMN850633186,9631252039Myoviridae; Felixounavirus91.2Escherichia phage vB_EcoM-AYO145ANC_028825.1
bumzenMN850635387,3601262039.1Myoviridae; Felixounavirus92.5Escherichia phage vB_EcoM_Alf5NC_031082.1
duneMN850636188,5111292039Myoviridae; Felixounavirus91.5Escherichia phage vB_EcoM_VpaE1NC_027337.1
warpigMN850637186,1061271739Myoviridae; Felixounavirus93Escherichia phage vB_EcoM_VpaE1NC_027337.1
radambzaMN850639186,7021271938.9Myoviridae; Felixounavirus91.6Escherichia phage vB_EcoM_VpaE1NC_027337.1
ekraMN850644187,2821282038.9Myoviridae; Felixounavirus92.9Escherichia phage vB_EcoM_Alf5NC_031082.1
humlepungMN850564385,3111191939.1Myoviridae; Felixounavirus92.1Escherichia phage vB_EcoM_VpaE1NC_027337.1
finnoMN850619187,5541292038.9Myoviridae; Felixounavirus89.7Escherichia phage vB_EcoM-AYO145ANC_028825.1
garusoMN850566285,7981302038.9Myoviridae; Felixounavirus90.9Escherichia phage vB_EcoM-AYO145ANC_028825.1
momoMN850580188,1681302039Myoviridae; Felixounavirus90.7Escherichia phage vB_EcoM-AYO145ANC_028825.1
heidMN850577687,5901262039Myoviridae; Felixounavirus91.2Escherichia phage vB_EcoM_Alf5NC_031082.1
skudenMN850585187,2631312038.9Myoviridae; Felixounavirus91.1Escherichia phage vB_EcoM_VpaE1NC_027337.1
pinkbiffMN850603188,8141292039Myoviridae; Felixounavirus93.9Escherichia phage vB_EcoM_Alf5NC_031082.1
fjerdesalMN850605387,7151282139Myoviridae; Felixounavirus90.6Escherichia phage vB_EcoM_AYO145ANC_028825.1
andreottiMN850610183,3911172039.2Myoviridae; Felixounavirus91.9Escherichia phage vB_EcoM_VpaE1NC_027337.1
nataliecMN850611289,1371342039Myoviridae; Felixounavirus90.3Escherichia phage vB_EcoM_AYO145ANC_028825.1
adrianhMN850614288,2261281938.9Myoviridae; Felixounavirus91.1Escherichia phage vB_EcoM_Alf5NC_031082.1
mistaenktMN850587186,6641282247.2Myoviridae; Suspvirus91.1Escherichia phage SUSP2NC_028935.2
nimiMN8506261137,039213543.7Myoviridae; Vequintavirus93.3Escherichia phage LL12MH491969.1
navnMN8506421141,707224443.6Myoviridae; Vequintavirus91.1Escherichia coli O157 typing phage 4KP869102.1
nomineMN8506491137,991220543.6Myoviridae; Vequintavirus91.5Escherichia phage LL12MH491969.1
naswaMN8505951138,583222543.6Myoviridae; Vequintavirus93.1Escherichia phage LL12MH491969.1
naamMN8506301137,129215543.7Myoviridae; Vequintavirus94.5Escherichia coli O157 typing phage 4KP869102.1
imeMN8505762137,114217543.6Myoviridae; Vequintavirus93.1Escherichia phage LL12MH491969.1
magacaMN8506121135,826217543.6Myoviridae; Vequintavirus96Escherichia phage slur12LN881735.1
nomMN8506461136,114213543.6Myoviridae; Vequintavirus92.6Escherichia phage LL12MH491969.1
isimMN8505971138,289219543.6Myoviridae; Vequintavirus93.8Escherichia phage LL12MH491969.1
nomoMN8505781137,702218543.7Myoviridae; Vequintavirus93.3Escherichia phage APCEc02NC_041869.1
inoaMN8505931138,710220543,6Myoviridae; Vequintavirus92Escherichia phage APCEc02NC_041869.1
pangalanMN8506273136,944215543.7Myoviridae; Vequintavirus94.8Escherichia phage vB_EcoM_FFH2NC_024134.1
tuntematonMN8506182150,4732791139.1Myoviridae; Phapecoctavirus89.6Escherichia phage phAPEC8NC_020079.1
anhysbysMN8506481149,3352711139.1Myoviridae; Phapecoctavirus91.5Escherichia phage phAPEC8NC_020079.1
ukendtMN8505651150,9472661139Myoviridae; Phapecoctavirus88.7Escherichia phage phAPEC8NC_020079.1
nepoznatoMN8505714151,5142651038.9Myoviridae; Phapecoctavirus85.6Escherichia phage phAPEC8NC_020079.1
nieznanyMN8505981144,9982541139.1Myoviridae; Phapecoctavirus88.9Escherichia phage phAPEC8NC_020079.1
muutMN8505731146,3072431337.4Myoviridae92Escherichia phage vB_EcoM_PHB05MF805809.1
aliaMN8506321147,0092461337.5Myoviridae93.1Enterobacteria phage ECGD1KU522583.1
outraMN8506451145,4822461337.4Myoviridae93.8Enterobacteria phage ECGD1KU522583.1
innyMN8506011147,4832471337.4Myoviridae92.4Enterobacteria phage ECGD1KU522583.1
arallMN8505841145,7152421337.4Myoviridae94.6Escherichia phage vB_vPM_PD06MH816848.1
kviMN8506151163,673266-40.5Myoviridae; Krischvirus94.2Escherichia phage ECD7NC_041936.1
kaaroeMN8505741163,719267-40.5Myoviridae; Krischvirus94.7Enterobacteria phage RB49NC_005066.1
dhabilMN8506211165,644266339.5Myoviridae; Dhakavirus87.5Enterobacteria phage JS10NC_012741.1
dhaegMN8506092170,817278339.4Myoviridae; Dhakavirus87.4Enterobacteria phage JS10NC_012741.1
mograMN8505791168,724263237.7Myoviridae; Mosigvirus91.1Escherichia phage vB_EcoM_PhAPEC2NC_024794.1
mobilluMN8506221163,063255237.7Myoviridae; Mosigvirus94.5Escherichia phage p000yMK047718.1
mohaMN8505901168,676267237.6Myoviridae; Mosigvirus94.8Escherichia phage APCEc01NC_029091.1
moskryMN8506511169,410269237.6Myoviridae; Mosigvirus93.2Escherichia virus vB_Eco_mar005P1LR027383.1
teqskovMN8954371165,017257635.4Myoviridae; Tequatrovirus91.7Yersinia phage phiD1NC_027353.1
teqdroesMN8954381166,8332691035.4Myoviridae; Tequatrovirus88.6Escherichia phage T2LC348380.1
teqhadMN8954341167,8922701035.3Myoviridae; Tequatrovirus90.1Escherichia phage T2LC348380.1
teqhalMN8954352168,0702661135.4Myoviridae; Tequatrovirus93.9Escherichia phage slur13LN881737.1
teqsoenMN8954361166,4682681035.5Myoviridae; Tequatrovirus91.7Yersinia phage phiD1NC_027353.1
flopperMN850594152,09278144.2Chaseviridae; Carltongylesvirus87Escherichia phage ST32NC_047830.1
damhausMN850602151,15489-44.1Drexlerviridae; Hanrivervirus85.8Shigella phage pSf-1 NC_021331.1
herniMN850640250,97189-44.1Drexlerviridae; Hanrivervirus87.6Shigella phage pSf-1 NC_021331.1
gramsMN850567149,53083-44.1Drexlerviridae; Hanrivervirus87.1Shigella phage pSf-1 NC_021331.1
aaroesMN850572151,66292-44.1Drexlerviridae; Hanrivervirus83Shigella phage pSf-1 NC_021331.1
aalborvMN850591146,66079-43.9Drexlerviridae; Hanrivervirus86.9Shigella phage pSf-1 NC_021331.1
haarsleMN850600248,61385-44Drexlerviridae; Hanrivervirus87.1Shigella phage pSf-1 NC_021331.1
egaaMN850607151,64387-44.1Drexlerviridae; Hanrivervirus89.7Shigella phage pSf-1 NC_021331.1
vojenMN850569150,70986-44.1Drexlerviridae; Hanrivervirus89.7Shigella phage pSf-1 NC_021331.1
tiwnaMN850643151,01485-44.6Drexlerviridae; Warwickvirus87.2Escherichia phage vB_Eco_Swan01NC_048202.1
tonijnMN850641251,62786-44.6Drexlerviridae; Warwickvirus88.4Escherichia phage vB_Eco_Swan01NC_048202.1
tonnikalaMN850613151,27786-44.8Drexlerviridae; Warwickvirus86.4Escherichia phage vB_Eco_Swan01NC_048202.1
atunaMN850620150,73288-44.6Drexlerviridae; Warwickvirus84.9Escherichia virus vB_Eco_mar001J1 NC_048204
tunusMN850638151,11187-44.8Drexlerviridae; Warwickvirus93.7Escherichia phage SECphi27NC_047938.1
orkinosMN850586249,79881-44.6Drexlerviridae; Warwickvirus91.3Escherichia phage SECphi27NC_047938.1
ityhunaMN850582150,76886-44.7Drexlerviridae; Warwickvirus93.3Escherichia phage SECphi27NC_047938.1
tonn MN850596251,01287-44.5Drexlerviridae; Warwickvirus94Escherichia phage vB_Eco_Swan01NC_048202.1
tinusoMN850634150,85686-44.8Drexlerviridae; Warwickvirus97.3Escherichia phage vB_Eco_Swan01NC_048202.1
tunzivisMN850604150,59684-44.6Drexlerviridae; Warwickvirus94.5Escherichia phage SECphi27NC_047938.1
tuinnMN850606150,50586-44.7Drexlerviridae; Warwickvirus94.8Escherichia phage vB_Eco_Swan01NC_048202.1
JahatMK552105151,10187-45.7Drexlerviridae68.5Escherichia phage vB_Eco_Swan01NC_048202.1
bobMN850628145,25263-54.5Siphoviridae; Dhillonvirus88.6Escherichia phage SECphi18LT960609.1
mckayMN850629144,44363-54.5Siphoviridae; Dhillonvirus83.8Escherichia phage slur05NC_028901.1
jatMN850650144,41763-54.5Siphoviridae; Dhillonvirus89.4Escherichia phage GluttonyNC_031113.1
rollingMN850575146,01764-54.2Siphoviridae; Dhillonvirus80.2Escherichia phage SlothKX534339.1
welshMN850589245,20762-54.6Siphoviridae; Dhillonvirus83.8Escherichia phage B2KX534339.1
buksMN850616140,30862-49.7Siphoviridae; Jerseyvirus91.3Salmonella phage vB_SenS-Ent1NC_019539.1
SkureMK672798159,47492-44.6Siphoviridae; Seuratvirus90.4Escherichia phage vB_Eco_SLUR25LT907986.1
HalfdanMH362766142,85857-53.7Siphoviridae28.8Pseudomonas phage vB_PaeS_SCUT-S3MK165657.1
LilleenMK629526153426-46.9Microviridae; Gequatrovirus93.8Escherichia phage SECphi17LT960607.1
LilleputMK629525154906-47Microviridae; Gequatrovirus93.4Escherichia phage SECphi17LT960607.1
LilletoMK629529354926-46.8Microviridae; Gequatrovirus92.7Escherichia phage SECphi17LT960607.1
LilleduMK791318154836-47.2Microviridae; Gequatrovirus92.6Escherichia phage SECphi17LT960607.1
lillemerMN850599154926-47.1Microviridae; Gequatrovirus94.5Escherichia phage SECphi17LT960607.1
LillevenMK629527160909-44.4Microviridae; Alphatrevirus93.9Enterobacteria phage St-1NC_012868.1
sortsynMN850623142,11661-59Podoviridae; Murrayvirus92.3Enterobacteria phage IME_EC2KF591601.1
sortregnMN850588138,20053-59.3Podoviridae; Murrayvirus97.3Salmonella phage LumpaelNC_048113.1
SkarpretterMK105855142,04263-55.8Podoviridae; Skarprettervirus37.9Escherichia phage C130_2MH363708.1
sortkaffMN850581142,53861-59.5Podoviridae; Sortsnevirus89.8Klebsiella phage vB_KpnS_IME279MF614100.1
SortsneMK651787141,91262-60Podoviridae; Sortsnevirus67.6Klebsiella phage vB_KpnS_IME279MF614100.1
aldrigsurMN850592142,37955-55.7Autographviridae; Bonnellvirus71.9Enterobacteria phage J8-65NC_025445.1
altidsurMN850568142,19753-55.7Autographviridae; Bonnellvirus71.8Enterobacteria phage J8-65NC_025445.1
forsurMN850617142,47656-55.4Autographviridae; Bonnellvirus72Enterobacteria phage J8-65NC_025445.1
glasurMN850583142,50756-55.4Autographviridae; Bonnellvirus72.3Enterobacteria phage J8-65NC_025445.1
LidtsurMK629528142,29156-54.6Autographviridae; Bonnellvirus69Enterobacteria phage J8-65NC_025445.1
megetsurMN850608142,13254-55.8Autographviridae; Bonnellvirus73.1Enterobacteria phage J8-65NC_025445.1
mellemsurMN850570140,77050-55.8Autographviridae; Bonnellvirus76.4Enterobacteria phage J8-65NC_025445.1
smaasurMN850625141,11050-55.4Autographviridae; Bonnellvirus93.3Enterobacteria phage J8-65NC_025445.1
usurMN850624241,90651-55.4Autographviridae; Bonnellvirus73.3Enterobacteria phage J8-65NC_025445.1

Share and Cite

MDPI and ACS Style

Olsen, N.S.; Forero-Junco, L.; Kot, W.; Hansen, L.H. Exploring the Remarkable Diversity of Culturable Escherichia coli Phages in the Danish Wastewater Environment. Viruses 2020, 12, 986. https://0-doi-org.brum.beds.ac.uk/10.3390/v12090986

AMA Style

Olsen NS, Forero-Junco L, Kot W, Hansen LH. Exploring the Remarkable Diversity of Culturable Escherichia coli Phages in the Danish Wastewater Environment. Viruses. 2020; 12(9):986. https://0-doi-org.brum.beds.ac.uk/10.3390/v12090986

Chicago/Turabian Style

Olsen, Nikoline S., Laura Forero-Junco, Witold Kot, and Lars H. Hansen. 2020. "Exploring the Remarkable Diversity of Culturable Escherichia coli Phages in the Danish Wastewater Environment" Viruses 12, no. 9: 986. https://0-doi-org.brum.beds.ac.uk/10.3390/v12090986

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop