Next Article in Journal
The Preclinical and Clinical Progress of Bacteriophages and Their Lytic Enzymes: The Parts are Easier than the Whole
Next Article in Special Issue
The Characteristics and Genome Analysis of vB_AviM_AVP, the First Phage Infecting Aerococcus viridans
Previous Article in Journal
New Viruses from the Ectoparasite Mite Varroa destructor Infesting Apis mellifera and Apis cerana
Previous Article in Special Issue
Isolation of A Novel Bacillus thuringiensis Phage Representing A New Phage Lineage and Characterization of Its Endolysin
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 11
Issue 2
10.3390/v11020095
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

Metagenomic Characterization of the Viral Community of the South Scotia Ridge

by
Qingwei Yang
1,
Chen Gao
1,
Yong Jiang
1,2,3,*,
Min Wang
1,2,3,*,
Xinhao Zhou
1,
Hongbing Shao
1,
Zheng Gong
1 and
Andrew McMinn
1,4
1
College of Marine Life Sciences, Ocean University of China, Qingdao 266003, China
2
Key Lab of Polar Oceanography and Global Ocean Change, Ocean University of China, Qingdao 266003, China
3
Institute of Evolution and Marine Biodiversity, Ocean University of China, Qingdao 266003, China
4
Institute for Marine and Antarctic Studies, University of Tasmania, Hobart, Tasmania 7001, Australia
*
Authors to whom correspondence should be addressed.
Viruses 2019, 11(2), 95; https://doi.org/10.3390/v11020095
Submission received: 24 October 2018 / Revised: 28 December 2018 / Accepted: 22 January 2019 / Published: 24 January 2019
(This article belongs to the Special Issue Diversity and Evolution of Phage Genomes)
Browse Figures
Versions Notes

Abstract

:
Viruses are the most abundant biological entities in aquatic ecosystems and harbor an enormous amount of genetic diversity. Whereas their influence on marine ecosystems is widely acknowledged, current information about their diversity remains limited. We conducted a viral metagenomic analysis of water samples collected during the austral summer of 2016 from the South Scotia Ridge (SSR), near the Antarctic Peninsula. The taxonomic composition and diversity of the viral communities were investigated, and a functional assessment of the sequences was performed. Phylotypic analysis showed that most viruses belonged to the order Caudovirales, especially the family Podoviridae (41.92–48.7%), which is similar to the situation in the Pacific Ocean. Functional analysis revealed a relatively high frequency of phage-associated and metabolism genes. Phylogenetic analyses of phage TerL and Capsid_NCLDV (nucleocytoplasmic large DNA viruses) marker genes indicated that many sequences associated with Caudovirales and NCLDV were novel and distinct from known phage genomes. High Phaeocystis globosa virus virophage (Pgvv) signatures were found and complete and partial Pgvv-like were obtained, which influence host–virus interactions. Our study expands existing knowledge of viral communities and their diversities from the Antarctic region and provides basic data for further exploring polar microbiomes.

1. Introduction

Viruses exist wherever life is found, including in extreme locations such as the deep ocean and polar areas. Arguably, viruses are by far the most numerous, genetically diverse, and pervasive biological entities on earth [1,2], including in aquatic ecosystems [3,4]. They are critical mortality agents of both eukaryotes and prokaryotes, affecting the abundance and diversity of microbial communities as well as global biogeochemical processes and energy fluxes, by causing lysis of a large proportion of both autotrophic and heterotrophic prokaryotes, shunting nutrients between particulate and dissolved phases [5,6,7,8,9,10], and modifying the efficiency of the carbon pump [11]. The genetic diversity of bacteria and protists was shaped by virus-mediated horizontal gene transfer, allowing viral genes to spread widely [4,12].
Both the ecology of Antarctic prokaryotes and protists [13,14,15,16,17,18] and the major role of viruses in prokaryotic and eukaryotic phytoplankton mortality [19,20,21,22] have been well studied. However, due to the geographical isolation and difficulty of accessing the Antarctic and Sub-Antarctic and the culturing viral hosts, an understanding of virus diversity and viral community structures in these regions is still lacking. There are comparatively few studies based on culture-independent methods, such as metagenomics and single-cell genomics, on the DNA and RNA of viral communities in Antarctica, although there are a few from freshwater habitats [23,24,25,26,27], the Southern Ocean close to the Western Antarctic Peninsula [28], and sediment soils [29,30]. These studies all identified a high viral biodiversity in these Antarctic ecosystems. However, despite the virome diversity information derived from these special habitats in the Antarctic, few studies have been completed in the open sea near the Antarctic Peninsula.
In this study, we conducted an analysis of viromes from three South Scotia Ridge (SSR) seawater samples including two from the surface and one from the bottom (water depth = 521 m) in an area influenced by Antarctic Circumpolar Current flow (ACC) [31]. The taxonomic composition of these viromes and the dominant viral species were identified and compared with viromes from other areas and habitats. A phylogenetic and/or genomic analysis of the representatives was undertaken.

2. Methods and Materials

2.1. Sample Collection and Sequencing

Seawater samples, including two surface waters and one from bottom water, at 5 m above the sediment–water interface (Table S1), were collected during the austral summer (December 2016) from two sites (D39 close to the edge of the Powell Basin and DA4 near the Clarence and Elephant Islands, Figure S1) on the southern flank of the SSR. Seawater temperature and salinity were recorded with a CTD profiler (SBE9/11 plus V5.2, Sea-Bird Inc., Newport, RI, USA). The temperatures of the seawater samples ranged from −0.04 to −0.57 °C and the salinity from 34.37‰ to 34.57‰ (Figure S1, Table S1). Water for biological and chemical analysis was collected with Niskin bottles attached to the CTD profiler and was prefiltered with 20-μm mesh to remove large particles.
The virome samples were processed immediately, according to the process used by Sun et al. [32]. Briefly, the samples were sequentially filtered through 3-μm and 0.22-μm pore size filters to remove any microorganisms, and then a two-step tangential flow filtration (TFF) with a 50-kDa cartridge (Pellicon® XL Cassette, Biomax® 50 kDa; polyethersulfone, Millipore Corporation, Billerica, MA, USA) was used to concentrate the viruses to a final volume of ca. 50 mL and stored at −80 °C. The samples were further concentrated by polyethylene glycol (PEG-8000) precipitation (10% w/v) and incubated at 4 °C overnight. The concentrate was then centrifuged at 8000× g for 80 min at 4 °C and suspended in 200 µL of SM buffer (0.05 M TRIS, 0.1 M NaCl, 0.008 M MgSO4, 0.01% (w/v) gelatin pH 7.5). Finally, DNA was extracted using the phenol/chloroform/isoamyl method and precipitated with ethanol without random amplification. High-throughput sequencing was performed by Novogene (Beijing, China) using Illumina Hiseq X ten (Paired End sequencing, 2 × 150 bp).

2.2. Virome Composition Analysis

The paired-end reads were quality trimmed by adopting the following conditions: (1) they contained more than 10% N, (2) had an adapter, and (3) were of low quality (70% read length, Q ≤ 30, and 80% read length, Q ≤ 20). All clean sequence data with quality-controlled were submitted to the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) under the following accession numbers: PRJNA505984 [33].
In order to avoid chimeras, SSR virome sequences were analyzed without assembly and queried by Diamond [34] against the NCBI non-redundant (nr) protein database [35] and the RefSeq complete viral genomes protein references (viral RefSeq) database [36], setting a maximum E-value of 10−3. Taxonomic identification was assigned based on best similarities and the relative taxonomy was normalized against complete viral genome length and sequencing depth which is the number of reads annotated to the virus divided by the viral complete length and the total number of reads sequenced in this sample.

2.3. Virome Comparison Analysis

Twenty previously published viromes taxonomic compositions, with the same maximum E-value based on read number, were selected from MetaVir to compare with this study [37]. These were obtained from a variety of habitats, including six temperate freshwater lakes (Lake Bourget and Lake Pavin [38], Antarctic lakes [23], Lough Neagh [39], Tilapia Channel [40]), nine seawater sites from the eastern tropical South Pacific oxygen minimum zones (ETSP-OMZ) [41], the Indian Ocean [42], the high salinity Jiulong River Estuary [43], Dunk Island, Fitzroy Island, LA26S (near Vancouver Island in British Columbia, Canada), M1CS (Monterey Bay, CA, USA) of the Pacific Ocean Virome (POV) [44], the Arctic Ocean, and Sargasso Sea (SAR) [45], three deep-sea surface sediment samples (Arctic Ocean, Black Sea and Mediterranean Sea [46]), soil, and hypolithon [47] (Table 3). The relative taxonomic composition of each community was normalized as described above. The similarity search algorithm BLAST was performed on the three SSR viromes against the 20 viromes obtained in MetaVir. The taxonomic composition distance matrix, based on relative abundances, was used in the non-metric multidimensional scaling(nMDS) analysis to plot viromes (the metaMDS function with a Bray–Curtis dissimilarity index using the VEGAN package in R software [48]) and a PERMANOVA (Permutational multivariate analysis of variance) test (p-test) was also performed.
In order to completely compare viromes rather than only their small known fraction, a qualitative comparison of viromes based on sequence similarity (tBLASTX comparison) was computed as described in MetaVir [37]. Briefly, we completed tBLASTX searches of the sub-sample set of sequences (that is, 50,000 sequences of 100 bp were randomly extracted from each entire virome library) from every single library versus all other subsamples. A similarity score between virome A and virome B was then computed as the sum of the top High Scoring Pairs (HSPs) scores of virome A reads against virome B reads (E-value < 10−5). Finally, the resulting score matrix (i.e., similarity scores for all virome pairs) was used in the nMDS analysis to plot viromes as described above.

2.4. Metagenomic Assembly and Function Analysis

SSR virome assemblies were performed via a random subsampling approach as previously described [49]. They were designed to obtain as the longest possible contigs by reducing the microdiversity within the samples [50,51]. Briefly, the assembling strategy was based on random selection of a subset of the reads: 1% (75×), 5% (50×), 10% (50×), 25% (25×), 75% (25×), and 100% (1×) from each sample and then assembling these subsets individually with IDBA_UD (v 1.1.2) [52] using the default parameters. We combined contigs derived from all the assemblies of the same samples and removed those <500 bp. To this end, contigs were clustered at 90% global average nucleotide identity with cd-hit-est (v 4.7, options: -c 0.9 -n 8) [53]. The relative abundance of each non-redundancy (nr) contigs was determined based on the mapping of the quality-filtered reads to the contigs, computed with bowtie2 (v 2.3.3.1) [54] and SAMtools [55], using the default parameters (the total length of reads mapping to the contig divided by the contig length). Then, the nr contigs were uploaded to the IMG system [56,57] and analyzed with the standard operating procedure of the DOE-JGI Metagenome Annotation Pipeline (MAP v.4) [58]. Finally, the IMG genomes 3300028548, 3300028550, and 3300028925 were obtained. The functional content was further characterized using MG-RAST [59] (with MG-RAST accession number 4808192.3, 4808195.3 and 4808193.3 respectively), an online metagenome annotation service [60], which was used to compare data to the SEED Subsystems database using a maximum E-value of 10−5, a minimum identity of 60%, and a minimum alignment length of 15.

2.5. Phylogenetic Analysis

Two double-stranded DNA (dsDNA) markers were present: the phage terminase large-subunit domain (TerL), which was present in phages of the order Caudovirales (Terminase_6, PF03237), and the major capsid protein (MCP) gene, which was present in large eukaryotic DNA viruses (Capsid_NCLDV, PF04451). Both of these were used to construct the phylogenetic trees and the TerL sequences were dereplicated at the 97% nucleotide level using cd-hit [53]. These markers from the SSR virome genes were screened by the DOE-JGI Metagenome Annotation Pipeline and compared to the viral RefSeq database using BLASTP (E-value < 10−5) to recruit relevant reference sequences. All sequences were aligned at the amino acid level using MUSCLE [61] (using default parameters), manually inspected and trimmed as necessary, and both maximum likelihood (ML) trees (MCP and TerL) with 1000 bootstraps were constructed with the program FastTree (v2.1.10) [62] using a JTT + CAT model and an estimation of the gamma parameter. Finally, the data were visualized and displayed using iTOL (Interactive Tree of Life) [63].

2.6. Genomic Comparison

The Phaeocystis globosa virus virophage (Pgvv)-like genomes were annotated with RAST and predicted open reading frames (ORFs) were searched against the NCBI reference viral protein (taxid:10239) with online BLASTP [64]. The partial functional annotations of the Pgvv reference sequence was obtained from Yutin et al. [65]. Visualization of the genomes’ map comparisons was completed EasyFig [66].

3. Results

3.1. Overview of SSR Viromes

After extraction, a total of 129,710,606 paired-end 150 bp sequences, with 109,923,264 (84.75%) reads passing the quality screening (Table 1), were obtained. The best BLAST Hit (E-value < 10−3) affiliations of unassembled high-quality reads from the three data sets are consistent with viral metagenomes published so far, as more than three-quarters of the reads (75.7–88.24%) did not show any significant sequence similarity to current NCBI nr data (Figure 1a). According to the NCBI nr and viral RefSeq annotation, the reads classified as viruses were 3.31–10.87% and 2.68–6.61% respectively (Figure 1b). A comparison of the annotation results of NCBI nr and viral RefSeq, showed that the virus sequences annotated with virus in the NCBI nr database were more abundant than those in viral RefSeq (Figure 1c), indicating that a certain proportion of sequences belong to an unidentified virus that viral RefSeq excluded, such as the uncultured Mediterranean phage uvMED.

3.2. Taxonomic Diversity Analysis

The BLAST data results (against viral RefSeq) of the virome composition were visualized using the Krona tool (Figures S2–S5) [67], which showed that, as expected, the majority of viral reads (93.69–95.16%) with significant hits belonged to double-stranded DNA (dsDNA) viruses with no RNA stage. These were largely comprised of members of the Caudovirales comprising the families Podoviridae, Siphoviridae, and Myoviridae, with similarities to single-stranded DNA (ssDNA) viruses also observed (Table 2 and Table 3). Podoviridae sequences (41.92–48.7%) were the most abundant in all three viromes, followed by Myoviridae (22.92–29.46%) and Siphoviridae (11.92–14.08%). Viruses from Phycodnaviridae (infecting algae) and Mimiviridae (infecting amoebas and algae) were more abundant in surface waters than the bottom water (D39s: 3.57% and 0.16%, DA4s: 2.22table% and 0.22%, DA4b: 1.32% and 0.10%, respectively). There was a significant proportion of virophages that prey on phycodnaviruses in the surface water samples: approximately 2% in D39s (Table 2). The top 10 most abundant viral species (Figure S6) included Puniceispirillum phage HMO-2011 [68] (Podoviridae, circular genome), which infects a bacterium of the SAR116 clade and was the most abundant in the SSR virome (18.50–25.75%), accounting for ~25% in DA4s, and Pelagibacter phage HTVC008M [69] (Myoviridae, linear), a T4-like myovirus infecting a SAR11 bacteriophage, which was the second most abundant (8.6–11.11%).

3.3. Comparison with Other Viromes

To compare viromes from the present study with previously published data sets, 20 viromes from different habitats were selected from MetaVir (see Materials and Methods for details). The result showed that the three SSR viromes were most closely related to ocean surface samples, except for the samples from ETSP-OMZ and SAR (Figure 2 and Figure S5) (p < 0.001). At the ocean surface, virome composition at the family level was dominated by the Caudovirales (Myoviridae, Siphoviridae, and Podoviridae), which collectively contributed 43.74–92.03% of the genomes. Viromes from within special habitats, including deep-ocean surface sediments, ETSP-OMZs, Antarctic freshwater, soil, and hypolithon, are dominated by Ciroviridae and Microviridae, members of ssDNA viruses, which contributed 25.45–88.45% of the genomes (Figure 3, Table 3). Less than 5% of these viromes’ sequences showed any similarity (E-value < 10−3) to the SSR viromes (Viromes were highlighted in bold shown in Table 4).

3.4. Contigs and Function Analysis

As the contigs assembled by the random subsampling approach could still contain redundant sequences derived from the same (or closely related) populations contigs derived from the same population were merged into clusters with 90% global average nucleotide identity using cd-hit-est. This resulted in 145,023 (D39s), 135,910 (DA4s), and 234,648 (DA4b) non-redundant genome fragments (>500 bp) (Table 1). Of these, 43.32% (D39s), 35.07% (DA4s), and 55.49% (DA4b) quality-filtered reads were assigned to nr contigs.
The putative functions of the annotated ORFs from the nr contigs dataset were predicted using MG-RAST, which assigns sequences to metabolic categories based on their Best BLAST Hit against the SEED database (E-value < 10−5). Using the subsystems approach, nearly 25% (17.54–26.46%) of the annotated proteins fell into ‘Phage, Prophage, Transposable elements, or Plasmids’ (Figure 4). Phage structural, integration/excision, and DNA metabolism-related proteins were most commonly identified and 10–11.96% of them were classified into “Clustering-based subsystems”, with phage endolysin commonly found in this category. The other SEED functional annotation categories showed that the metabolism of amino acids, carbohydrates, cofactors, vitamins, proteins, RNA, DNA, and nucleosides/nucleotides were the dominant annotations. In these categories, many proteins, such as DNA polymerases and helicase, could be phage-related (or of possible cellular origin). These hits were also found in the Pfam and COGs databases (see IMG system), with ‘Replication, recombination and repair’ being the most common protein categories identified.

3.5. Phylogenetic Tree Analysis

3.5.1. Terminase Phylogeny

The ML phylogenetic analysis of the phage large terminase subunits identified in this work is shown in Figure 5. The topology of the phylogenetic tree clearly shows that the majority of the SSR viromes’ TerL amino acid sequences were widely distributed among the Myo-, Sipho-, and Podoviridae. Among them, several branches (black dotted line) were moderately related to known members of the Guernseyvirinae family, T5virus, Luz24likevirus, and T4virus genera, and a few sequences were relatively closely related to those of known cultured representatives, including Pelagibacter phage HTVC010P, Rhodothermus phage RM378, and Cyanophage P-RSM6. Most sequences, however, were phylogenetically distant to known complete phage genomes (black solid lines). Notably, six groups (bright blue solid line) did not cluster with any known species and formed novel phylogenetic clusters. This separation is supported by high bootstrap values, which highlight the important but uncharacterized diversity of the Caudovirales in SSR.

3.5.2. Capsid_NCLDV Phylogeny

An ML phylogenetic tree, based on the MCP, which includes a group of putative MCP of Pgvv-like infected Phaeocystis globosa virus (Pgv), is shown in Figure 6. The MCP tree shows that several sequences from the SSR viromes are closely related to known NCLDV—mainly those belonging to Phycodnaviruses, and these can be classified into Prasinovirus, Pgv, and Pgvv. The three clades differed from Phycodnaviruses and Mimivirus MCPs and formed three distinct groups with the well-supported clades. One of the clades, marked as Group3, was only found in the surface ocean of DA4 station. The Pgv group, which included five new Pgvv-like MCPs, was distantly related to the Pgv group and had a higher relative abundance in the surface samples than the bottom samples.

3.6. Novel Pgvv Group

From the MCP phylogenetic tree, one distinct group of virophages was defined for which there is one known related virophage genome. However, this group still different from the known Pgvv. An alignment of Pgvv-like group genomes is shown in Figure 7. The Pgvv-like group genomes appear to have a relatively high GC content (37.36–38.17%), which was expected as the GC content of the Pgvv-like 04 genome (GC, 35.85%) was similar to Pgvv (GC, 35.8%). All virophages share four homologous proteins or domains: (1) packaging ATPase (ATPase), (2) lipase, (3) major capsid protein (MCP), (4) minor capsid protein (mCP). In addition, Pgvv-like 02 also contains the OLV11-like tyrosine recombinase (Yrec) gene, which is distantly related to the OLV11-like family [65]. Three genes with functional annotation (shown in yellow), which were absent in the Pgvv genome, were carried by the Pgvv-like sequences, including putative primase-helicase and DNA methyltransferase genes in the Pgvv-like 02 and recombination endonuclease VII gene in Pgvv-like 04. These characteristics further indicate that these viruses may belong to a new Pgvv-like group.

4. Discussion

Marine viral communities are still largely undescribed and many basic features, such as their global ocean distribution and their actual genetic and species richness, remain unknown [4,11,70]. With the advent of metagenomic methods, an association with high-depth sequencing, and meta-analyses of bioinformatics, an increasing number of studies have been conducted [71,72]. So far, only a few of these have focused on viral communities from the Antarctic region and most of these are from unusual habitats, such as freshwater lakes [23], hydrothermal vents [73], and soils [30].The important role of marine dsDNA viral communities viruses (that is, those capable of both lysogeny and lytic replication) in the Western Antarctic Peninsula has recently been demonstrated [28]. Also, the major differences in viral community composition between the subtropical Indian and the Southern Oceans have been identified [74].
The number of reads identified as either bacteria or eukaryote was similar to that reported in viral metagenomes of other environments [70,75]. In addition, the relatively low number of rRNA and tRNA genes (<1%) matching sequences (Table S2) indicates a certain degree of bacterial and eukaryotic contamination of the metagenomes; this has previously been reported to occur with TFF-based concentration methods [76]. One possible reason is that bacterial genes can be packaged into generalized transduced phage particles [77,78]. The bacterial-like sequences might have originated from excised prophages, mistakenly annotated as bacterial, and/or from genes of bacterial origins that were transferred to their phages [75]. Another problem is that cd-hit-est does not de-replicate full circular genomes that are linearized at different assembly sites, leading to circular viruses of the same source potentially merging into different clusters, which does not affect the reads-based taxonomic identification.
BLAST searches showed that more than 75% of the sequences before assembly did not have homologs in current sequence databases. This is consistent with the results of previously published viral metagenomic projects [44,74,79,80]. The SSR viromes were mostly dominated by Caudovirales, including Myoviruses, Siphoviruses, and Podoviruses, which are the dominant viral types recovered during metagenomic analyses of most marine environments [4,79]. In the three SSR viromes investigated here, the largest number of reads (>40%) were related to podoviruses and ~13% of reads were of siphoviruses (viruses that infect photosynthetic bacteria such as Prochlorococcus and Synechococcus) (in bold, Table S3). Consistent with a previous investigation [68,81], Puniceispirillum phage HMO-2011, which infects Candidatus Puniceispirillum marinum strain IMCC1322 of the SAR116 clade, and the Pelagibacter phage group (HTVC008M, HTVC010P, HTVC011P, and HTVC019P), infects SAR11 populations were widespread and most abundant in the SSR. Both SAR11 and SAR 116 clades play important roles in oceanic dimethylsulfide (DMS) production and biogeochemical sulfur cycles, especially via bacteria-mediated dimethylsulfoniopropionate (DMSP) degradation [82,83]. Marine viruses are likely to indirectly influence the global sulfur cycle by mediating the death of both hosts. Interestingly, Pgv is the tenth most abundant viral species in the SSR region (2.53% in D93s), infecting the temperate algal species Phaeocystis globosa [84]. In the Antarctic, however, the most abundant Phaeocystis species is P. antarctica [85], but a P. antarctica-specific virus has not yet been isolated or identified, which may suggest a high genome similarity between P. antarctica virus and Pgv. Compared with the surface viromes, a relatively smaller number of Phycodnaviridae and their virophage were still identified in the bottom virome where eukaryotic algae cannot photosynthesize. The origin and activity of these viruses needs further study.
Despite being in a cold marine environment with an average temperature below 0 °C, the SSR viral community had a similar structure to those found in the Pacific Ocean. However, there were still significant differences in nucleic acid levels. Is likely that the genotype of many viruses changed, allowing them to infect psychrophiles and thus evolve into new viral groups. The previously studied viromes from deep-ocean surface sediments, ETSP-OMZs, Antarctic freshwater, soil and hypolithon, in which ssDNA viruses played dominated roles, were clearly different from those of the SSR. However, all of those viromes, except those from the deep-ocean surface sediment, were amplified using multiple displacement amplification (MDA) with phi29 polymerase. In these genomes of the ssDNA, viruses were selectively amplified [86,87], leading to an overestimation of the role of ssDNA viruses. Although existence bias from MDA in these studies and the prevalence of Caudovirales sequences has been observed in most marine viromes, previously published research on global morphological analysis of marine viruses, conducted by the Tara Oceans Expedition, showed that non-tailed viruses (largely ssDNA and RNA) numerically dominate the upper oceans [88], and small, non-tailed viruses were undoubtedly underestimated in the SSR region.
The deep sequencing method, combined with a random subsampling assembly approach, enable obtaining a nearly complete viral genome and undertaking phylogenetic analyses on marker genes. Analysis of the major viral groups found in the SSR viromes showed broad diversity with many previously unknown virotypes. The terminase gene, which is responsible for DNA recognition and initiation of DNA packaging, is an essential component of all head-tail phages (Caudovirales), as it encodes the molecular movements that translocate DNA into empty capsids [89]. There is a large diversity of terminases that can be used to resolve different Caudoviruses groups [90]. The NCLDV comprises a monophyletic group of viruses infecting both animals and a diverse range of unicellular eukaryotes, including the Phycodna-, Mimi-, Asco-, Asfar-, Irido-, and Poxviridae families. The MCP of NCLDV (capsid_NCLDV), a redox protein that encodes complex DNA replication and transcription systems and involved in the formation of disulfide bond in virion membrane proteins, is relatively conserved among NCLDVs evolution [91,92,93]. Using phylogenetic trees based on these two viral marker genes (TerL and MCP), a high diversity among Caudovirales and NCLDV was identified. In a single habitat, the high diversity of viruses not only expressed in genotypes, but also in morphological and biological properties such as plaque morphology, originated from a broad range of hosts and propagation temperatures [73]. A high proportion of TerL sequences was distributed both far from the reference and far from each other, highlighting both the richness of Caudovirales in the SSR communities and the absence of closely-related reference sequences. In addition, some SSR virome sequences appear to have formed a new clade (Group 6) related to the T4 viruses, one of the best described Caudovirales families.
The topology of the MCP tree and genomic comparisons strongly suggest that the five putative virophage genomes are more closely related to the Pgvv than to other NCDLV families, including the Pgvv host. The Pgvv-like group also has a high relative abundance. The Lotka–Volterra simulation demonstrated that virophages promote secondary production through the microbial loop by reducing overall mortality of the algal cell after a bloom and increasing the frequency of blooms during the summer [20]. According to the above model, it can be inferred that the Pgvv-like group plays a previously unrecognized role in regulating virus-host interactions in the SSR area during summer.

5. Conclusions

Analysis of the SSR viromes has shown that novel, oceanic-related viromes are present. A high proportion of sequence reads was classified as unknown, with only 3.31–10.87% having known virus counterparts. Among these, members of the order Caudovirales were most abundant. This pattern is consistent with previously described viromes from the Pacific Ocean as well as from a range of different biomes. The diversity of the Caudovirales and NCLDV in the SSR viromes is high, suggesting that viral diversity is high in gelid environments. However, the abundance and diversity of ssDNA and RNA viruses require further research. The strong signatures of Pgvv found in the SSR may indicate that virophages play an important role in regulating virus-host interactions.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4915/11/2/95/s1, Figure S1: Map of area around the South Scotia Ridge (SSR) indicating where samples for metagenomic analysis were collected; Figure S2: Krona chart representing taxonomic composition of the sequence reads in the D39s from SSR surface seawater. Relative abundance of the sequence reads classified by the taxonomic grouping based on BLASTX similarity search (E-value <10−3); Figure S3: Krona chart representing taxonomic composition of the sequence reads in the DA4s from SSR surface seawater; Figure S4: Krona chart representing taxonomic composition of the sequence reads in the D39s from SSR bottom seawater; Figure S5: Comparison of viromes between SSR area and other environmental viromes depending on known taxonomic composition (stress value 6.76%); Figure S6: Top 10 relative abundance of viral species in the three SSR virome; Table S1: Details of the sampling sites where the viral metagenomes have been collected in the South Scotia Ridge; Table S2: General statistics of the gene annotation of the three SSR viromes; Table S3: Most represented viral genotypes among the viral hits according to complete normalized viral genome length and the total number of reads sequenced.

Author Contributions

M.W. and Y.J. conceived and designed the experiments; H.S. sampled the seawater; Q.Y., X.Z. and Z.G. performed the experiments; Q.Y., C.G. and X.Z. analyzed the data; Q.Y., Y.J. and A.M. wrote and modified the paper.

Funding

We are grateful for the funding from the National Natural Science Foundation of China (Nos. 41676178, 41076088, and 31500339), National Key Research and Development Program of China (2017YFA0603200), Scientific and Technological Innovation Project Financially Supported by Qingdao National Laboratory for Marine Science and Technology (No. 2016ASKJ14), and Fundamental Research Funds from the Central University of Ocean University of China (Grant Nos. 201812002, 201762017, 201562018).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rohwer, F. Global phage diversity. Cell 2003, 113, 141. [Google Scholar] [CrossRef]
  2. Lipkin, W.I. Microbe hunting. Microbiol. Mol. Biol. Rev. 2010, 74, 363–377. [Google Scholar] [CrossRef] [PubMed]
  3. Fuhrman, J.A. Marine viruses and their biogeochemical and ecological effects. Nature 1999, 399, 541–548. [Google Scholar] [CrossRef] [PubMed]
  4. Suttle, C.A. Viruses in the sea. Nature 2005, 437, 356–361. [Google Scholar] [CrossRef] [PubMed]
  5. Proctor, L.M.; Fuhrman, J.A. Viral mortality of marine-bacteria and cyanobacteria. Nature 1990, 343, 60–62. [Google Scholar] [CrossRef]
  6. Suttle, C.A. The significance of viruses to mortality in aquatic microbial communities. Microb. Ecol. 1994, 28, 237–243. [Google Scholar] [CrossRef] [PubMed]
  7. Brussaard, C.P.D.; Wilhelm, S.W.; Thingstad, F.; Weinbauer, M.G.; Bratbak, G.; Heldal, M.; Kimmance, S.A.; Middelboe, M.; Nagasaki, K.; Paul, J.H.; et al. Global-scale processes with a nanoscale drive: The role of marine viruses. ISME J. 2008, 2, 575–578. [Google Scholar] [CrossRef]
  8. Rohwer, F.; Thurber, R.V. Viruses manipulate the marine environment. Nature 2009, 459, 207–212. [Google Scholar] [CrossRef]
  9. 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]
  10. Jover, L.F.; Effler, T.C.; Buchan, A.; Wilhelm, S.W.; Weitz, J.S. The elemental composition of virus particles: Implications for marine biogeochemical cycles. Nat. Rev. Microbiol. 2014, 12, 519–528. [Google Scholar] [CrossRef]
  11. Suttle, C.A. Marine viruses—Major players in the global ecosystem. Nat. Rev. Microbiol. 2007, 5, 801–812. [Google Scholar] [CrossRef] [PubMed]
  12. Mann, N.H.; Cook, A.; Millard, A.; Bailey, S.; Clokie, M. Marine ecosystems: Bacterial photosynthesis genes in a virus. Nature 2003, 424, 741. [Google Scholar] [CrossRef] [PubMed]
  13. Kuparinen, J.; Bjornsen, P.K. Spatial-distribution of bacterioplankton production across the weddell-scotia confluence during early austral summer 1988–1989. Pol. Biol. 1992, 12, 197–204. [Google Scholar] [CrossRef]
  14. Archer, S.D.; Leakey, R.J.G.; Burkill, P.H.; Sleigh, M.A. Microbial dynamics in coastal waters of east antarctica: Herbivory by heterotrophic dinoflagellates. Mar. Ecol. Prog. Ser. 1996, 139, 239–255. [Google Scholar] [CrossRef]
  15. Vaque, D.; Calderon-Paz, J.I.; Guixa-Boixereu, N.; Pedros-Alio, C. Spatial distribution of microbial biomass and activity (bacterivory and bacterial production) in the northern weddell sea during the austral summer (January 1994). Aquat. Microb. Ecol. 2002, 29, 107–121. [Google Scholar] [CrossRef]
  16. Christaki, U.; Obernosterer, I.; Van Wambeke, F.; Veldhuis, M.; Garcia, N.; Catala, P. Microbial food web structure in a naturally iron-fertilized area in the southern ocean (kerguelen plateau). Deep-Sea Res. Part II 2008, 55, 706–719. [Google Scholar] [CrossRef]
  17. Montes-Hugo, M.; Doney, S.C.; Ducklow, H.W.; Fraser, W.; Martinson, D.; Stammerjohn, S.E.; Schofield, O. Recent changes in phytoplankton communities associated with rapid regional climate change along the western antarctic peninsula. Science 2009, 323, 1470–1473. [Google Scholar] [CrossRef]
  18. Moschonas, G.; Gowen, R.J.; Paterson, R.F.; Mitchell, E.; Stewart, B.M.; McNeill, S.; Glibert, P.M.; Davidson, K. Nitrogen dynamics and phytoplankton community structure: The role of organic nutrients. Biogeochemistry 2017, 134, 125–145. [Google Scholar] [CrossRef]
  19. Guixa-Boixereu, N.; Vaque, D.; Gasol, J.M.; Sanchez-Camara, J.; Pedros-Alio, C. Viral distribution and activity in antarctic waters. Deep-Sea Res. Part II 2002, 49, 827–845. [Google Scholar] [CrossRef]
  20. Yau, S.; Lauro, F.M.; DeMaere, M.Z.; Brown, M.V.; Thomas, T.; Raftery, M.J.; Andrews-Pfannkoch, C.; Lewis, M.; Hoffman, J.M.; Gibson, J.A.; et al. Virophage control of antarctic algal host-virus dynamics. Proc. Natl. Acad. Sci. USA 2011, 108, 6163–6168. [Google Scholar] [CrossRef] [Green Version]
  21. Malits, A.; Christaki, U.; Obernosterer, I.; Weinbauer, M.G. Enhanced viral production and virus-mediated mortality of bacterioplankton in a natural iron-fertilized bloom event above the kerguelen plateau. Biogeosciences 2014, 11, 6841–6853. [Google Scholar] [CrossRef]
  22. Vaque, D.; Boras, J.A.; Torrent-Llagostera, F.; Agusti, S.; Arrieta, J.M.; Lara, E.; Castillo, Y.M.; Duarte, C.M.; Sala, M.M. Viruses and protists induced-mortality of prokaryotes around the antarctic peninsula during the austral summer. Front. Microbiol. 2017, 8, 241. [Google Scholar] [CrossRef] [PubMed]
  23. Lopez-Bueno, A.; Tamames, J.; Velazquez, D.; Moya, A.; Quesada, A.; Alcami, A. High diversity of the viral community from an antarctic lake. Science 2009, 326, 858–861. [Google Scholar] [CrossRef] [PubMed]
  24. Zawar-Reza, P.; Arguello-Astorga, G.R.; Kraberger, S.; Julian, L.; Stainton, D.; Broady, P.A.; Varsani, A. Diverse small circular single-stranded DNA viruses identified in a freshwater pond on the mcmurdo ice shelf (antarctica). Infect. Genet. Evol. 2014, 26, 132–138. [Google Scholar] [CrossRef] [PubMed]
  25. De Carcer, D.A.; Lopez-Bueno, A.; Pearce, D.A.; Alcami, A. Biodiversity and distribution of polar freshwater DNA viruses. Sci Adv. 2015, 1, e1400127. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. De Carcer, D.A.; Lopez-Bueno, A.; Alonso-Lobo, J.M.; Quesada, A.; Alcami, A. Metagenomic analysis of lacustrine viral diversity along a latitudinal transect of the antarctic peninsula. FEMS Microbiol. Ecol. 2016, 92. [Google Scholar]
  27. Cavicchioli, R.; Erdmann, S. The discovery of antarctic rna viruses: A new game changer. Mol. Ecol. 2015, 24, 4809–4811. [Google Scholar] [CrossRef]
  28. Brum, J.R.; Hurwitz, B.L.; Schofield, O.; Ducklow, H.W.; Sullivan, M.B. Seasonal time bombs: Dominant temperate viruses affect southern ocean microbial dynamics. ISME J. 2016, 10, 437–449. [Google Scholar] [CrossRef]
  29. Hopkins, M.; Kailasan, S.; Cohen, A.; Roux, S.; Tucker, K.P.; Shevenell, A.; Agbandje-McKenna, M.; Breitbart, M. Diversity of environmental single-stranded DNA phages revealed by pcr amplification of the partial major capsid protein. ISME J. 2014, 8, 2093–2103. [Google Scholar] [CrossRef]
  30. Adriaenssens, E.M.; Kramer, R.; Van Goethem, M.W.; Makhalanyane, T.P.; Hogg, I.; Cowan, D.A. Environmental drivers of viral community composition in antarctic soils identified by viromics. Microbiome 2017, 5, 83. [Google Scholar] [CrossRef]
  31. Howe, J.A.; Livermore, R.A.; Maldonado, A. Mudwave activity and current-controlled sedimentation in powell basin, northern weddell sea, antarctica. Mar. Geol. 1998, 149, 229–241. [Google Scholar] [CrossRef]
  32. Sun, G.; Xiao, J.; Wang, H.; Gong, C.; Pan, Y.; Yan, S.; Wang, Y. Efficient purification and concentration of viruses from a large body of high turbidity seawater. MethodsX 2014, 1, 197–206. [Google Scholar] [CrossRef] [PubMed]
  33. SRA. Available online: https://www.ncbi.nlm.nih.gov/sra/PRJNA505984 (accessed on 20 November 2018).
  34. Buchfink, B.; Xie, C.; Huson, D.H. Fast and sensitive protein alignment using diamond. Nat. Methods 2015, 12, 59–60. [Google Scholar] [CrossRef] [PubMed]
  35. Download data from NCBI FTP site. Available online: ftp://ftp.ncbi.nlm.nih.gov/blast/db/FASTA/ (accessed on 18 July 2018).
  36. Download data from NCBI FTP site. Available online: ftp://ftp.ncbi.nlm.nih.gov/refseq/release/viral/ (accessed on 18 July 2018).
  37. Roux, S.; Faubladier, M.; Mahul, A.; Paulhe, N.; Bernard, A.; Debroas, D.; Enault, F. Metavir: A web server dedicated to virome analysis. Bioinformatics 2011, 27, 3074–3075. [Google Scholar] [CrossRef] [PubMed]
  38. Roux, S.; Enault, F.; Robin, A.; Ravet, V.; Personnic, S.; Theil, S.; Colombet, J.; Sime-Ngando, T.; Debroas, D. Assessing the diversity and specificity of two freshwater viral communities through metagenomics. PLoS ONE 2012, 7, e33641. [Google Scholar] [CrossRef] [PubMed]
  39. Skvortsov, T.; de Leeuwe, C.; Quinn, J.P.; McGrath, J.W.; Allen, C.C.R.; McElarney, Y.; Watson, C.; Arkhipova, K.; Lavigne, R.; Kulakov, L.A. Metagenomic characterisation of the viral community of lough neagh, the largest freshwater lake in Ireland. PLoS ONE 2016, 11, e0150361. [Google Scholar] [CrossRef] [PubMed]
  40. Rodriguez-Brito, B.; Li, L.L.; Wegley, L.; Furlan, M.; Angly, F.; Breitbart, M.; Buchanan, J.; Desnues, C.; Dinsdale, E.; Edwards, R.; et al. Viral and microbial community dynamics in four aquatic environments. ISME J. 2010, 4, 739–751. [Google Scholar] [CrossRef] [Green Version]
  41. Cassman, N.; Prieto-Davo, A.; Walsh, K.; Silva, G.G.Z.; Angly, F.; Akhter, S.; Barott, K.; Busch, J.; McDole, T.; Haggerty, J.M.; et al. Oxygen minimum zones harbour novel viral communities with low diversity. Environ. Microbiol. 2012, 14, 3043–3065. [Google Scholar] [CrossRef]
  42. Williamson, S.J.; Allen, L.Z.; Lorenzi, H.A.; Fadrosh, D.W.; Brami, D.; Thiagarajan, M.; McCrow, J.P.; Tovchigrechko, A.; Yooseph, S.; Venter, J.C. Metagenomic exploration of viruses throughout the Indian ocean. PLoS ONE 2012, 7, e42047. [Google Scholar] [CrossRef]
  43. Cai, L.L.; Zhang, R.; He, Y.; Feng, X.Y.; Jiao, N.Z. Metagenomic analysis of virioplankton of the subtropical jiulong river estuary, China. Viruses 2016, 8, 35. [Google Scholar] [CrossRef]
  44. Hurwitz, B.L.; Sullivan, M.B. The pacific ocean virome (POV): A marine viral metagenomic dataset and associated protein clusters for quantitative viral ecology. PLoS ONE 2013, 8, e57355. [Google Scholar] [CrossRef] [PubMed]
  45. Angly, F.E.; Felts, B.; Breitbart, M.; Salamon, P.; Edwards, R.A.; Carlson, C.; Chan, A.M.; Haynes, M.; Kelley, S.; Liu, H.; et al. The marine viromes of four oceanic regions. PLoS Biol. 2006, 4, 2121–2131. [Google Scholar] [CrossRef] [PubMed]
  46. Danovaro, R.; Dell’Anno, A.; Corinaldesi, C.; Magagnini, M.; Noble, R.; Tamburini, C.; Weinbauer, M. Major viral impact on the functioning of benthic deep-sea ecosystems. Nature 2008, 454, 1084–1087. [Google Scholar] [CrossRef] [PubMed]
  47. Zablocki, O.; van Zyl, L.; Adriaenssens, E.M.; Rubagotti, E.; Tuffin, M.; Cary, S.C.; Cowan, D. High-level diversity of tailed phages, eukaryote-associated viruses, and virophage-like elements in the metaviromes of antarctic soils. Appl. Environ. Microbiol. 2014, 80, 6888–6897. [Google Scholar] [CrossRef]
  48. Dixon, P. Vegan, a package of r functions for community ecology. J. Veg. Sci. 2003, 14, 927–930. [Google Scholar] [CrossRef]
  49. Coutinho, F.H.; Silveira, C.B.; Gregoracci, G.B.; Thompson, C.C.; Edwards, R.A.; Brussaard, C.P.D.; Dutilh, B.E.; Thompson, F.L. Marine viruses discovered via metagenomics shed light on viral strategies throughout the oceans. Nat. Commun. 2017, 8, 15955. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Nagarajan, N.; Pop, M. Sequence assembly demystified. Nat. Rev. Genet. 2013, 14, 157–167. [Google Scholar] [CrossRef] [PubMed]
  51. De Crcer, D.A.; Angly, F.E.; Alcami, A. Evaluation of viral genome assembly and diversity estimation in deep metagenomes. Bmc Genomics 2014, 15, 989. [Google Scholar]
  52. Peng, Y.; Leung, H.C.M.; Yiu, S.M.; Chin, F.Y.L. IDBA-UD: A de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth. Bioinformatics 2012, 28, 1420–1428. [Google Scholar] [CrossRef]
  53. Li, W.; Godzik, A. Cd-hit: A fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 2006, 22, 1658–1659. [Google Scholar] [CrossRef]
  54. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [PubMed]
  55. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R.; Proc, G.P.D. The sequence alignment/map format and samtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [PubMed]
  56. Chen, I.M.A.; Markowitz, V.M.; Chu, K.; Palaniappan, K.; Szeto, E.; Pillay, M.; Ratner, A.; Huang, J.H.; Andersen, E.; Huntemann, M.; et al. Img/m: Integrated genome and metagenome comparative data analysis system. Nucleic Acids Res. 2017, 45, D507–D516. [Google Scholar] [CrossRef]
  57. DOE Joint Genome Insitute-IMG. Available online: https://img.jgi.doe.gov/ (accessed on 23 August 2018).
  58. Huntemann, M.; Ivanova, N.N.; Mavromatis, K.; Tripp, H.J.; Paez-Espino, D.; Tennessen, K.; Palaniappan, K.; Szeto, E.; Pillay, M.; Chen, I.M.A.; et al. The standard operating procedure of the DOE-JGI Metagenome Annotation Pipeline (MAP v.4). Stand. Genomic Sci. 2016, 11, 17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Meyer, F.; Paarmann, D.; D’Souza, M.; Olson, R.; Glass, E.M.; Kubal, M.; Paczian, T.; Rodriguez, A.; Stevens, R.; Wilke, A.; et al. The metagenomics rast server—A public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformatics 2008, 9, 386. [Google Scholar] [CrossRef] [PubMed]
  60. MG-RAST metagenomics analysis server. Available online: http://metagenomics.anl.gov/ (accessed on 9 September 2018).
  61. Edgar, R.C. Muscle: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [PubMed]
  62. Price, M.N.; Dehal, P.S.; Arkin, A.P. Fasttree 2-approximately maximum-likelihood trees for large alignments. PLoS ONE 2010, 5, e9490. [Google Scholar] [CrossRef] [PubMed]
  63. Letunic, I.; Bork, P. Interactive tree of life (ITOL) v3: An online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016, 44, W242–W245. [Google Scholar] [CrossRef]
  64. Altschul, S.F.; Madden, T.L.; Schaffer, A.A.; Zhang, J.H.; Zhang, Z.; Miller, W.; Lipman, D.J. Gapped blast and psi-blast: A new generation of protein database search programs. Nucleic Acids Res. 1997, 25, 3389–3402. [Google Scholar] [CrossRef]
  65. Yutin, N.; Shevchenko, S.; Kapitonov, V.; Krupovic, M.; Koonin, E.V. A novel group of diverse polinton-like viruses discovered by metagenome analysis. BMC Biol. 2015, 13, 95. [Google Scholar] [CrossRef]
  66. Sullivan, M.J.; Petty, N.K.; Beatson, S.A. Easyfig: A genome comparison visualizer. Bioinformatics 2011, 27, 1009–1010. [Google Scholar] [CrossRef]
  67. Ondov, B.D.; Bergman, N.H.; Phillippy, A.M. Interactive metagenomic visualization in a web browser. BMC Bioinformatics 2011, 12, 385. [Google Scholar] [CrossRef] [PubMed]
  68. Kang, I.; Oh, H.M.; Kang, D.; Cho, J.C. Genome of a sar116 bacteriophage shows the prevalence of this phage type in the oceans. Proc. Natl. Acad. Sci. USA 2013, 110, 12343–12348. [Google Scholar] [CrossRef]
  69. Eggleston, E.M.; Hewson, I. Abundance of two pelagibacter ubique bacteriophage genotypes along a latitudinal transect in the north and south atlantic oceans. Front. Microbiol. 2016, 7, 1534. [Google Scholar] [CrossRef] [PubMed]
  70. Edwards, R.A.; Rohwer, F. Viral metagenomics. Nat. Rev. Microbiol. 2005, 3, 504–510. [Google Scholar] [CrossRef] [PubMed]
  71. Mokili, J.L.; Rohwer, F.; Dutilh, B.E. Metagenomics and future perspectives in virus discovery. Curr. Opin. Virol. 2012, 2, 63–77. [Google Scholar] [CrossRef]
  72. Brum, J.R.; Sullivan, M.B. Rising to the challenge: Accelerated pace of discovery transforms marine virology. Nat. Rev. Microbiol. 2015, 13, 147–159. [Google Scholar] [CrossRef] [PubMed]
  73. Millard, A.D.; Pearce, D.; Zwirglmaier, K. Biogeography of bacteriophages at four hydrothermal vent sites in the antarctic based on g23 sequence diversity. FEMS Microbiol. Lett. 2016, 363, fnw043. [Google Scholar] [CrossRef] [PubMed]
  74. Flaviani, F.; Schroeder, D.C.; Lebret, K.; Balestreri, C.; Highfield, A.C.; Schroeder, J.L.; Thorpe, S.E.; Moore, K.; Pasckiewicz, K.; Pfaff, M.C.; et al. Distinct oceanic microbiomes from viruses to protists located near the antarctic circumpolar current. Front. Microbiol. 2018, 9, 1474. [Google Scholar] [CrossRef]
  75. Fancello, L.; Trape, S.; Robert, C.; Boyer, M.; Popgeorgiev, N.; Raoult, D.; Desnues, C. Viruses in the desert: A metagenomic survey of viral communities in four perennial ponds of the mauritanian sahara. ISME J. 2013, 7, 359–369. [Google Scholar] [CrossRef]
  76. Hurwitz, B.L.; Deng, L.; Poulos, B.T.; Sullivan, M.B. Evaluation of methods to concentrate and purify ocean virus communities through comparative, replicated metagenomics. Environ. Microbiol. 2013, 15, 1428–1440. [Google Scholar] [CrossRef] [PubMed]
  77. Beumer, A.; Robinson, J.B. A broad-host-range, generalized transducing phage (sn-t) acquires 16s rRNA genes from different genera of bacteria. Appl. Environ. Microb. 2005, 71, 8301–8304. [Google Scholar] [CrossRef] [PubMed]
  78. Ghosh, D.; Roy, K.; Williamson, K.E.; White, D.C.; Wommack, K.E.; Sublette, K.L.; Radosevich, M. Prevalence of lysogeny among soil bacteria and presence of 16s rRNA and tRZN genes in viral-community DNA. Appl. Environ. Microb. 2008, 74, 495–502. [Google Scholar] [CrossRef] [PubMed]
  79. Breitbart, M.; Salamon, P.; Andresen, B.; Mahaffy, J.M.; Segall, A.M.; Mead, D.; Azam, F.; Rohwer, F. Genomic analysis of uncultured marine viral communities. Proc. Natl. Acad. Sci. USA 2002, 99, 14250–14255. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Roux, S.; Enault, F.; Ravet, V.; Colombet, J.; Bettarel, Y.; Auguet, J.C.; Bouvier, T.; Lucas-Staat, S.; Vellet, A.; Prangishvili, D.; et al. Analysis of metagenomic data reveals common features of halophilic viral communities across continents. Environ. Microbiol. 2016, 18, 889–903. [Google Scholar] [CrossRef] [PubMed]
  81. Zhao, Y.L.; Temperton, B.; Thrash, J.C.; Schwalbach, M.S.; Vergin, K.L.; Landry, Z.C.; Ellisman, M.; Deerinck, T.; Sullivan, M.B.; Giovannoni, S.J. Abundant SAR11 viruses in the ocean. Nature 2013, 494, 357–360. [Google Scholar] [CrossRef] [PubMed]
  82. Howard, E.C.; Sun, S.L.; Biers, E.J.; Moran, M.A. Abundant and diverse bacteria involved in dmsp degradation in marine surface waters. Environ. Microbiol. 2008, 10, 2397–2410. [Google Scholar] [CrossRef]
  83. Choi, D.H.; Park, K.T.; An, S.M.; Lee, K.; Cho, J.C.; Lee, J.H.; Kim, D.; Jeon, D.; Noh, J.H. Pyrosequencing revealed SAR116 clade as dominant dddp-containing bacteria in oligotrophic nw pacific ocean. PLoS ONE 2015, 10, e0116271. [Google Scholar] [CrossRef]
  84. Baudoux, A.C.; Brussaard, C.P.D. Characterization of different viruses infecting the marine harmful algal bloom species phaeocystis globosa. Virology 2005, 341, 80–90. [Google Scholar] [CrossRef]
  85. DiTullio, G.R.; Grebmeier, J.M.; Arrigo, K.R.; Lizotte, M.P.; Robinson, D.H.; Leventer, A.; Barry, J.B.; VanWoert, M.L.; Dunbar, R.B. Rapid and early export of phaeocystis antarctica blooms in the ross sea, antarctica. Nature 2000, 404, 595–598. [Google Scholar] [CrossRef]
  86. Lasken, R.S.; Stockwell, T.B. Mechanism of chimera formation during the multiple displacement amplification reaction. BMC Biotechnol. 2007, 7, 19. [Google Scholar] [CrossRef] [PubMed]
  87. Kim, K.H.; Bae, J.W. Amplification methods bias metagenomic libraries of uncultured single-stranded and double-stranded DNA viruses. Appl. Environ. Microb. 2011, 77, 7663–7668. [Google Scholar] [CrossRef] [PubMed]
  88. Brum, J.R.; Schenck, R.O.; Sullivan, M.B. Global morphological analysis of marine viruses shows minimal regional variation and dominance of non-tailed viruses. ISME J. 2013, 7, 1738–1751. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  89. Maniloff, J.; Ackermann, H.W. Taxonomy of bacterial viruses: Establishment of tailed virus genera and the order caudovirales. Arch. Virol. 1998, 143, 2051–2063. [Google Scholar] [CrossRef] [PubMed]
  90. Mizuno, C.M.; Rodriguez-Valera, F.; Garcia-Heredia, I.; Martin-Cuadrado, A.B.; Ghai, R. Reconstruction of novel cyanobacterial siphovirus genomes from mediterranean metagenomic fosmids. Appl. Environ. Microb. 2013, 79, 688–695. [Google Scholar] [CrossRef] [PubMed]
  91. Tidona, C.A.; Schnitzler, P.; Kehm, R.; Darai, G. Is the major capsid protein of iridoviruses a suitable target for the study of viral evolution? Virus Genes 1998, 16, 59–66. [Google Scholar] [CrossRef] [PubMed]
  92. Iyer, L.M.; Aravind, L.; Koonin, E.V. Common origin of four diverse families of large eukaryotic DNA viruses. J. Virol. 2001, 75, 11720–11734. [Google Scholar] [CrossRef]
  93. Larsen, J.B.; Larsen, A.; Bratbak, G.; Sandaa, R.A. Phylogenetic analysis of members of the phycodnaviridae virus family, using amplified fragments of the major capsid protein gene. Appl. Environ. Microb. 2008, 74, 3048–3057. [Google Scholar] [CrossRef]
Viruses 11 00095 g001 550
Figure 1. Taxonomic assignment of metagenomic reads (a) Percentage of the sequence reads classified by the taxonomic grouping based on BLASTX similarity search with NCBI nr database (E-value < 10−3). Sequences with no hits with E-value > 10−3 were regarded as unidentified reads (“no annotation” category in the table and red in the pie graphs). “LUCA” (last universal common ancestor) (green) denotes reads that could not unambiguously be assigned to a domain of life. (b) Taxonomic assignment of metagenomic reads based on BLASTX similarity search with viral RefSeq database (E-value < 10−3). (c) Comparison of the annotation results of NCBI nr and viral RefSeq.
Figure 1. Taxonomic assignment of metagenomic reads (a) Percentage of the sequence reads classified by the taxonomic grouping based on BLASTX similarity search with NCBI nr database (E-value < 10−3). Sequences with no hits with E-value > 10−3 were regarded as unidentified reads (“no annotation” category in the table and red in the pie graphs). “LUCA” (last universal common ancestor) (green) denotes reads that could not unambiguously be assigned to a domain of life. (b) Taxonomic assignment of metagenomic reads based on BLASTX similarity search with viral RefSeq database (E-value < 10−3). (c) Comparison of the annotation results of NCBI nr and viral RefSeq.
Viruses 11 00095 g001
Viruses 11 00095 g002 550
Figure 2. Comparison viromes between the SSR area and other environmental viromes depending on BLAST-based comparison. Twenty environmental viromes were available on MetaVir2, obtained from different habitats including freshwater, seawater, deep-sea surface sediments, soil, and hypolith. Bray–Curtis dissimilarity matrices of BLAST hits were calculated from virome data and used to represent the relative distances between individual viromes (stress value 10.7%).
Figure 2. Comparison viromes between the SSR area and other environmental viromes depending on BLAST-based comparison. Twenty environmental viromes were available on MetaVir2, obtained from different habitats including freshwater, seawater, deep-sea surface sediments, soil, and hypolith. Bray–Curtis dissimilarity matrices of BLAST hits were calculated from virome data and used to represent the relative distances between individual viromes (stress value 10.7%).
Viruses 11 00095 g002
Viruses 11 00095 g003 550
Figure 3. The relative abundance of viral sequences (normalized with genome length) mostly at the family level in each different habitat virome. Point size indicates the value of relative abundance percentage.
Figure 3. The relative abundance of viral sequences (normalized with genome length) mostly at the family level in each different habitat virome. Point size indicates the value of relative abundance percentage.
Viruses 11 00095 g003
Viruses 11 00095 g004 550
Figure 4. Composition of predicted functional genes of the SSR contigs. The Coding Sequences (CDSs) were compared with the SEED database using subsystems in MG-RAST. The metabolic categorization is based on the sequences Best BLAST Hits in the SEED database curated subsystems (E-value < 10−5).
Figure 4. Composition of predicted functional genes of the SSR contigs. The Coding Sequences (CDSs) were compared with the SEED database using subsystems in MG-RAST. The metabolic categorization is based on the sequences Best BLAST Hits in the SEED database curated subsystems (E-value < 10−5).
Viruses 11 00095 g004
Viruses 11 00095 g005 550
Figure 5. Terminase phylogeny. A maximum-likelihood phylogenetic tree of Caudovirales terminase large-subunit domains (PF03237) is shown (1000 iterations, JTT + G model). Only bootstrap values >50% are indicated at the nodes of the tree, and bootstrap scores greater than 90% are indicated with a black dot. Average branch length distance of leaves less than 0.4 were collapsed and are shown as triangles. Reference sequences are marked (see color legend at the top). Abbreviations are as follows. HMO-2011: Puniceispirillum phage HMO-2011; HTVC010P: Pelagibacter phage HTVC010P; RM378: Rhodothermus phage RM378; P-RSM6: Cyanophage P-RSM6. The black dotted line, solid line, and bright blue line indicate the sequences obtained from this study.
Figure 5. Terminase phylogeny. A maximum-likelihood phylogenetic tree of Caudovirales terminase large-subunit domains (PF03237) is shown (1000 iterations, JTT + G model). Only bootstrap values >50% are indicated at the nodes of the tree, and bootstrap scores greater than 90% are indicated with a black dot. Average branch length distance of leaves less than 0.4 were collapsed and are shown as triangles. Reference sequences are marked (see color legend at the top). Abbreviations are as follows. HMO-2011: Puniceispirillum phage HMO-2011; HTVC010P: Pelagibacter phage HTVC010P; RM378: Rhodothermus phage RM378; P-RSM6: Cyanophage P-RSM6. The black dotted line, solid line, and bright blue line indicate the sequences obtained from this study.
Viruses 11 00095 g005
Viruses 11 00095 g006 550
Figure 6. Capsid_NCLDV phylogeny. A maximum-likelihood phylogenetic tree drawn from the capsid_NCLDV (PF04451) and six virophage putative major capsid (MCP) protein multiple alignment is shown (1000 iterations, JTT + G model). Bootstrap scores greater than 90% are marked with black dots. Each MCP is associated with an abundance profile (right) that displays the relative abundance of the contig across the three SSR viromes (based on normalized coverage).
Figure 6. Capsid_NCLDV phylogeny. A maximum-likelihood phylogenetic tree drawn from the capsid_NCLDV (PF04451) and six virophage putative major capsid (MCP) protein multiple alignment is shown (1000 iterations, JTT + G model). Bootstrap scores greater than 90% are marked with black dots. Each MCP is associated with an abundance profile (right) that displays the relative abundance of the contig across the three SSR viromes (based on normalized coverage).
Viruses 11 00095 g006
Viruses 11 00095 g007 550
Figure 7. Overview of genomic synteny and similarities between Pgvv-like group. The complete Pgvv reference is covered by a grey shadow. A color scale for percent identity (nucleic) is shown at the bottom right. The name, percent GC content (GC%), and length for each genome are indicated. Genes are colored according to their functional affiliation. Tlr 6F: Toll-like receptor 6 family; GIY: GIY-YIG family nuclease; MCP: major capsid protein; mCP: minor capsid protein; and Yrec (OLV11): OLV11-like tyrosine recombinase.
Figure 7. Overview of genomic synteny and similarities between Pgvv-like group. The complete Pgvv reference is covered by a grey shadow. A color scale for percent identity (nucleic) is shown at the bottom right. The name, percent GC content (GC%), and length for each genome are indicated. Genes are colored according to their functional affiliation. Tlr 6F: Toll-like receptor 6 family; GIY: GIY-YIG family nuclease; MCP: major capsid protein; mCP: minor capsid protein; and Yrec (OLV11): OLV11-like tyrosine recombinase.
Viruses 11 00095 g007
Table
Table 1. The remaining number of quality-controlled reads and non-redundancy contigs. Q20: quality score of 20; Q30: quality score of 30.
Table 1. The remaining number of quality-controlled reads and non-redundancy contigs. Q20: quality score of 20; Q30: quality score of 30.
ProcessSSR Viromes
D39sDA4sDA4b
Quality controlRaw reads4,2547,32443,791,90843,371,374
Cut adaptor40,670,620 (95.59%)41,883,752 (95.64%)41,353,540 (95.35%)
Q20 >20%35,357,306 (83.10%)37,567,224 (85.79%)38,397,082 (88.53%)
Q30 >30%35,334,370 (83.05%)36,385,270 (83.09%)38,203,624 (88.08%)
AssembleAll assembled contigs (>500 bp)2,418,0813,699,5591,693,019
Non-redundancy contigs145,023135,910234,648
MappingMapped reads15,307,526 (43.32%)12,760,519 (35.07%)21,197,289 (55.49%)
Table
Table 2. Classification of reads from viromes hitting viral sequences. ssDNA: single-stranded DNA; dsDNA: double-stranded DNA.
Table 2. Classification of reads from viromes hitting viral sequences. ssDNA: single-stranded DNA; dsDNA: double-stranded DNA.
GroupOrderFamilyRelative Abundance (%)
D39sDA4sDA4b
dsDNACaudoviralesPodoviridae41.9248.7042.15
CaudoviralesMyoviridae28.3422.9229.46
CaudoviralesSiphoviridae11.9213.3814.08
Caudoviralesunclassified3.562.992.90
CaudoviralesPhycodnaviridae3.572.221.32
--Mimiviridae0.160.220.10
--Poxviridae0.140.221.13
--Iridoviridae0.370.440.30
unclassified 3.132.843.17
ssDNA--Inoviridae0.100.130.12
--Microviridae0.000.100.04
--Circoviridae0.010.110.01
--unclassified0.090.050.05
virophage-- 2.000.580.19
othersOrterviralesRetroviridae0.260.290.21
OrterviralesCaulimoviridae0.070.760.04
Unclassified phage/viruses-- 3.713.014.07
Table
Table 3. Taxonomic characterization of South Scotia Ridge (SSR) viromes and twenty previously published viromes as determined by the indicated BLAST comparison to viral Refseq (E-value < 10−3).
Table 3. Taxonomic characterization of South Scotia Ridge (SSR) viromes and twenty previously published viromes as determined by the indicated BLAST comparison to viral Refseq (E-value < 10−3).
ViromeMetaVir
ID.		Total Sequenced
Reads		Viral
Hits (%)		Relative Abundance (%)
dsDNA Viruses, No RNA StageRNAssDNASatelliteVirophage
AckermannBaculoMyoPodoSiphoHerpesIridoMimiPapillomaPhycodnPolyomPoxOthersCircoGeminiInoMicroNanoOthers
Antar Lake Spring1041,32220.51000.050.060.2200000.01000.030.1422.000.270.033.450.2223.396.820
Antar Lake Summer1138,47522.0700.010.320.300.720.030.010.3801.570.020.090.200.4061.960.490.111.310.9315.523.200.19
Lake Pavin6649,29025.04001.601.296.730.010.010.0100.060.010.020.8528.046.020.070.121.790.0617.070.930.25
Lake Bourget7593,08437.460.0102.302.9811.560.020.010.0100.140.020.061.860.281.770.030.0160.350.321.990.360.20
Tilapia Channel33264,8449.180.030.025.4011.4211.840.040.030.0400.090.0602.770.0219.8913.460.428.260.1612.313.800.03
Lough Neagh49252,295,05524.650.140.0810.9632.5936.710.110.080.020.060.380.050.3116.110.320.010.030.160.0400.3900.13
Jiulong River Estuary6305498,95731.400.050.0614.6732.3634.800.120.040.110.020.790.160.4913.580.2700.010.110.010.010.220.310.14
GS1171479480,37538.800.040.0316.5035.5927.050.060.040.060.010.530.060.4417.590.170.010.010.020.010.010.070.320.28
M1CS1440303,51932.930.050.0310.9228.0848.750.090.060.060.020.470.070.159.850.110.030.010.0100.040.110.290.11
Dunk Island13571,165,2563.610.010.0212.6140.7421.800.190.040.040.060.360.040.1522.430.0900.050.02000.070.530.20
Fitzroy Island135882,73932.660.020.0513.3348.0017.300.040.050.050.060.440.120.0119.360.2000.040.04000.0900.29
LA26S1396165,25628.450.040.1715.5144.0524.320.080.100.110.021.170.270.2711.760.1000.020.07000.210.760.19
D39s--17,667,18511.970.090.0628.3341.9011.920.030.330.1603.570.050.147.040.040.0100.07000.1002.00
DA4s--19,101,8125.090.060.2222.9148.6813.370.050.410.2202.210.060.226.330.100.1100.100.1000.0700.58
DA4b--18,192,63513.870.100.0529.4642.1414.080.020.260.1001.320.011.136.440.050.0100.090.0400.0700.19
Sargasso Sea12397,93911.800.0102.206.322.4800.020.0200.080.0402.880.0114.210043.690.0324.830.600.01
Arctic Vir15686,2091.300.040.049.919.6323.140.540.080.250.051.241.370.763.190.990.640.045.18001.8816.070.06
Arctic Ocean115879,64614.15000.120.120.39000.0100.030.010.020.090.3263.140.090.162.170.1913.934.250
Black Sea11557843616.47000.050.030.0800000.0200.010.020.1256.330.600.010.880.1714.245.790
Med Sea116165,34016.1300.012.290.822.490.040.050.020.010.160.041.580.350.2921.050059.0902.967.510
OMZst3viral10m897128,44115.560.010.012.904.002.120.010.020.0400.400.030.061.070.0860.7900.036.490.1314.891.880.16
Antar hypolith27261,057,53512.740.010.011.311.477.150.130.0100.050.020.030.020.360.130.0300.2188.4200.250.090
Antar open soil2727870,68716.990.080.059.254.9616.740.290.050.030.060.240.050.093.100.430.0100.0363.2800.400.070
Table
Table 4. The percentage of reads in other published viromes obtained from MetaVir with a significant similarity (BLASTN, E-value < 10−3) to the SSR viromes.
Table 4. The percentage of reads in other published viromes obtained from MetaVir with a significant similarity (BLASTN, E-value < 10−3) to the SSR viromes.
BiomeViromeMetaVir
Project ID.		Number of ReadsSouth Scotia Sea
D39sDA4sDA4b
Antarctic seawaterD39s-35,334,370100%50.35%49.60%
Antarctic seawaterDA4s-36,385,27048.90%100%28.37%
Antarctic seawaterDA4b-38,203,62445.88%27.02%100%
SeawaterOMZst3viral10m897128,4416.50%6.86%10.69%
SeawaterGS1171479480,3758.61%9.09%14.57%
Arctic seawaterArctic Vir15686,2091.45%1.33%2.10%
SeawaterSargasso Sea12397,9394.68%4.97%7.69%
POV seawaterDunk Island13571,165,2560.55%0.71%1.14%
POV seawaterFitzroy Island135882,7397.44%9.11%14.04%
POV seawaterLA26S1396165,25616.75%14.72%21.02%
POV seawaterM1CS1440303,51914.77%15.59%23.12%
Deep OceanArctic Ocean115879,6461.64%2.17%4.22%
Deep OceanBlack Sea115578,4360.57%0.62%0.69%
Deep OceanMed Sea116165,3400.86%0.89%2.04%
FreshwaterLake Bourget7593,0840.94%0.87%1.66%
FreshwaterLake Pavin6649,2900.25%0.23%0.44%
Antarctic freshwaterAntar Lake Spring1041,3220.07%0.06%0.16%
Antarctic freshwaterAntar Lake Summer1138,4750.41%0.43%1.87%
FreshwaterLough Neagh49252,295,0550.31%0.30%0.52%
FreshwaterJiulong River Estuary6305498,9575.74%6.06%9.61%
FreshwaterTilapia Channel33264,8440.14%0.15%0.31%
Antarctic soilAntar open soil2727870,6870.43%0.42%0.59%
Antarctic hypolithAntar hypolith27261,057,5350.26%0.42%0.36%

Share and Cite

MDPI and ACS Style

Yang, Q.; Gao, C.; Jiang, Y.; Wang, M.; Zhou, X.; Shao, H.; Gong, Z.; McMinn, A. Metagenomic Characterization of the Viral Community of the South Scotia Ridge. Viruses 2019, 11, 95. https://doi.org/10.3390/v11020095

AMA Style

Yang Q, Gao C, Jiang Y, Wang M, Zhou X, Shao H, Gong Z, McMinn A. Metagenomic Characterization of the Viral Community of the South Scotia Ridge. Viruses. 2019; 11(2):95. https://doi.org/10.3390/v11020095

Chicago/Turabian Style

Yang, Qingwei, Chen Gao, Yong Jiang, Min Wang, Xinhao Zhou, Hongbing Shao, Zheng Gong, and Andrew McMinn. 2019. "Metagenomic Characterization of the Viral Community of the South Scotia Ridge" Viruses 11, no. 2: 95. https://doi.org/10.3390/v11020095

APA Style

Yang, Q., Gao, C., Jiang, Y., Wang, M., Zhou, X., Shao, H., Gong, Z., & McMinn, A. (2019). Metagenomic Characterization of the Viral Community of the South Scotia Ridge. Viruses, 11(2), 95. https://doi.org/10.3390/v11020095

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