Next Article in Journal
Begomovirus-Associated Satellite DNA Diversity Captured Through Vector-Enabled Metagenomic (VEM) Surveys Using Whiteflies (Aleyrodidae)
Next Article in Special Issue
Identification of Known and Novel Recurrent Viral Sequences in Data from Multiple Patients and Multiple Cancers
Previous Article in Journal
The Role of HBZ in HTLV-1-Induced Oncogenesis
Previous Article in Special Issue
HPV Population Profiling in Healthy Men by Next-Generation Deep Sequencing Coupled with HPV-QUEST
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 8
Issue 2
10.3390/v8020035
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 Analysis of Virioplankton of the Subtropical Jiulong River Estuary, China

by
Lanlan Cai
1,
Rui Zhang
1,*,
Ying He
2,
Xiaoyuan Feng
2 and
Nianzhi Jiao
1,*
1
State Key Laboratory of Marine Environmental Science, Institute of Marine Microbes and Ecospheres, Xiamen University (Xiang’an), Xiamen 361005, China
2
State Key Laboratory of Microbial Metabolism, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
*
Authors to whom correspondence should be addressed.
Viruses 2016, 8(2), 35; https://doi.org/10.3390/v8020035
Submission received: 29 November 2015 / Revised: 21 January 2016 / Accepted: 25 January 2016 / Published: 2 February 2016
(This article belongs to the Special Issue Next Generation Sequencing: New Developments and Discoveries in Virology)
Browse Figures
Versions Notes

Abstract

:
Viruses are the most abundant biological entities in the oceans, and encompass a significant reservoir of genetic diversity. However, little is known about their biodiversity in estuary environments, which represent a highly dynamic and potentially more diverse habitat. Here, we report a metagenomic analysis of the dsDNA viral community from the Jiulong River Estuary (JRE), China, and provide a comparative analysis with other closely related environments. The results showed that the majority of JRE virome did not show any significant similarity to the database. For the major viral group (Caudovirales) detected in the sample, Podoviridae (44.88%) were the most abundant family, followed by Siphoviridae (32.98%) and Myoviridae (17.32%). The two most abundant viruses identified in the virome were phages HTVC010P and HMO-2011, which infect bacteria belonging to marine SAR11 and SAR116 clades, respectively. Two contigs larger than 20 kb, which show similar overall genome architectures to Celeribacter phage P12053L and Thalosomonas phage BA3, respectively, were generated during assembly. Comparative analysis showed that the JRE virome was more similar to marine viromes than to freshwater viromes, and shared a relative coarse-grain genetic overlap (averaging 14.14% ± 1.68%) with other coastal viromes. Our study indicated that the diversity and community structure of the virioplankton found in JRE were mainly affected by marine waters, with less influence from freshwater discharge.

1. Introduction

Over the last three decades, viruses, particularly those infecting bacteria (bacteriophages), have received increasing recognition for their importance in marine environments. This recognition is based upon four main lines of evidence. First, viruses are the most abundant organisms in the oceans, comprising approximately 94% of the nucleic acid-containing particles, and encompassing a significant reservoir of genetic diversity [1,2]. Second, viruses keep microbial numbers at values less than the carrying capacity of the system through viral lysis. This is well known as the “kill the winner” hypothesis, and plays an important role in shaping microbial population structures and mediating genetic exchange between hosts [3,4]. Third, viruses are able to reprogram the metabolic output of their hosts through the expression of virally encoded auxiliary metabolic genes, which play critical roles in facilitating biochemical and metabolic processes [5,6]. For example, cyanophages can transcribe and express photosynthesis genes during lytic infection, helping to maintain photosynthesis during the phage lytic cycle and thereby promote viral propagation [7,8]. Finally, referred to as the “viral shunt”, the lytic processes of the host microbial cells infected by marine viruses will release dissolved organic matters into the water, which can be used by other heterotrophic bacteria. This means that virus-mediated mortality increases net respiration, dissolved organic matter recycling and nutrient transformation in aquatic systems [2,9,10,11].
Viral metagenome (virome) analysis overcomes the limitation of traditional methods and circumvents the drawback that there are no universal signature genes that can be used as phylogenetic markers to assess viral diversity. It has emerged as a promising method not only for studying viral diversity in a wide range of research samples, such as marine environments [12,13], fresh waters [14,15,16,17], plants [18], feces [19], blood [20], animal tissue [21], and other specialized habitats [22,23] but also for discovering novel viruses [13,24]. Combined findings from previous virome surveys suggest that viruses carry a tremendous genetic diversity. The majority (usually, 60%–99%) of sequences in viromes from any environment do not share homology with databases [25,26,27]. Recently, the Tara Oceans Expedition (2009–2013) carried out a global survey of viromes by collecting a wide variety of samples from the surface, deep chlorophyll maximum and mesopelagic zones of seven oceans and seas. The large-scale dataset confirms a high local viral diversity, which may result from localized pockets of viruses under small-scale environmental constraints [28]. However, so far, the majority of virome analyses have focused on marine or fresh waters and, to a lesser extent, transition zones, such as estuaries. Viral communities in estuary environments are likely to be more diverse than those of their marine and freshwater counterparts. The mixing zone between marine tidal and river outflow, together with nutritional surplus, is a highly dynamic environment that can trigger physiological, genetic, and ecological changes in microbes, and thus exert a potentially profound effect on food web interactions and estuarine biogeochemical cycling [29,30]. The Pacific Ocean Viromes dataset has revealed decreased functional richness with distance from shore in surface waters [31]. A previous study of viral diversity in a gradient of anthropogenic impacts demonstrated that the mixed sample produced the highest viral genotype richness [14]. Meanwhile, estuaries are of great ecological significance as they often serve as spawning sites for fish and shellfish species, and are the connection between watersheds and oceans. In spite of their significant importance, relatively little is known about the composition and structure of viral communities in estuaries compared with other aquatic environments.
The Jiulong River is the second longest river in Fujian Province, southeast China. The Jiulong River estuary (JRE) is a typical subtropical macro-tide estuary connected with Xiamen Sea harbor and on the southwest coast of the Taiwan Strait. A recently published study that investigated the distribution and diversity of bacterioplankton communities of Xiamen Sea harbor found that the sampling site near the JRE exhibited the highest bacterial diversity, and was significantly influenced by total nitrogen and total phosphorus [32]. Because of the nature of the host-specific interactions, we could expect the composition of viruses (most of which are phages) to follow that of the bacterial community. In this study, we present a double-stranded DNA (dsDNA) viral metagenome dataset from the JRE. Given the high level of plasticity of viral genomes, the current virome analysis improves our knowledge of aquatic viral ecology and assists in identifying viral populations specific to subtropical estuarine environment.

2. Materials and Methods

2.1. Sample Collection and Sequencing

About 150 L of surface water were collected from the JRE, China (118°16′54.56″ E, 24°30′17.00″ N). No specific permits or permissions were required for the sampling. Temperature, salinity, and pH were measured in situ with a YSI Professional Plus multiparameter meter (YSI Incorporated, Yellow Springs, OH, USA). Samples were processed immediately according to Thurber et al. [33]. Briefly, the samples were filtered through a 3-μm filter to remove large particles. Two-step tangential flow filtration (TFF) was then used to concentrate viruses. First, a 0.22-μm pore size filter (0.5 m2 Pellicon cartridge, Millipore Corp., Billerica, MA, USA) was used to remove bacteria and small eukaryotes. Viral particles were further concentrated using a TFF 30-kDa cartridge to a final volume of 2 L. The samples were precipitated using polyethylene glycol (PEG-8000) (10% w/v) and incubated at 4 °C overnight. Virus-like particles from the PEG pellet were purified using cesium chloride gradient ultracentrifugation (1.7, 1.5, and 1.35 g/mL layers). Virus-like particles were collected from the 1.35–1.5 g/mL fraction and confirmed by epifluorescence microscopy. After being treated with DNase I and RNase A at 4 °C for 1 h, the purified virus-like particles were treated with Proteinase K, EDTA (0.5 M) and 10% SDS at 55 °C for 3 h. DNA was extracted using phenol/chloroform/isoamylol method. To test for contaminating bacterial DNA, viral nucleic acids were screened by 16S rDNA PCR with 27F/1492R universal primers. After DNA quality and quantity checks by DNA agarose gel electrophoresis, DNA was sequenced using Roche 454 GS FLX system (Roche 454 Life Sciences, Branford, CT, USA).

2.2. Microbial Enumeration

For microbial counts, 2-mL aliquots were fixed in a final concentration of 0.5% glutaraldehyde at 4 °C for 20 min, then stored at −80 °C after snap freezing in liquid nitrogen. Viral and heterotrophic bacterial abundances were determined by flow cytometry (FCM) according to Marie et al. [34] and Bruaassrd [35]. Briefly, the fixed frozen samples were thawed at room temperature and stained with SYBR Green I (Invitrogen, Carlsbad, CA, USA). The stained particles were enumerated using a flow cytometer (Epics Altra II, Beckman Coulter, Miami, FL, USA) at event rates of 50–200 particles/s (bacteria) or 100–300 particles/s (virus). The autotrophic picoplankton abundance was determined according to the method of Jiao et al. [36]. Picoeukaryotes and Synechococcus were distinguished according to their positions in the cytometric plots. The analyses of microbial abundance by FCM for each sample were repeated twice. The data were analyzed with EXPOTM32 MultiCOMP software [37].

2.3. Bioinformatics Analysis

A series of filtering and trimming steps were undertaken to remove low quality reads using cd-hit 454 [38] and Mothur v.1.33.0 [39]. Clean reads were subsequently assembled using GS De novo Assembler (Newbler v2.5, Roche 454 Life Sciences) with default parameters. Both clean reads and assembled contigs from the JRE virome were uploaded to Metavir [40], an online tool for analyzing viral genomic data, and the Meta Genome Rapid Annotation using Subsystem Technology (MG-RAST) server [41], a general metagenomic webserver, with MG-RAST accession number 4,527,838.3. Metavir performs taxonomic composition analysis based on BLASTx searches against the National Center for Biotechnology Information (NCBI) viral refseq database. MG-RAST generates taxonomic assignments based on BLASTx searches against the M5NR database and functional categories based on BLASTx searches against the SEED-Subsystem database. Contigs larger than 20 kb obtained from the above assembly were examined further. Putative open reading frames (ORFs) within each contig were predicted using GeneMarkS [42], and then searched against the NCBI non-redundant protein database using BLASTp. BLAST homologies having an E-value less than 10−3 were considered putative hits. Comparisons to publicly available aquatic viromes were carried out in Metavir by performing tBLASTx comparisons after data normalization. To obtain an in-depth comparison of the virome obtained in this study with closely related datasets, we downloaded eight Pacific Ocean coastal surface virome datasets (PCSVs) from a previous study, including samples GFS, GDS, LJ4S, MACS, SFSS, SFDS, SFCS, and STCS [31]. All datasets were subsampled to a common depth (the minimum number of sequences in the GFS virome; 82,739), and then trimmed to 100-bp to homogenize the sequence length and effectively normalize the sampling effort. A cross-tblastx (E-value < 10−3) was undertaken to study the degree of putative coarse-grain genetic overlap between these environments [17].

3. Results

3.1. Environmental Characteristics

As an estuary environment, the JRE undergoes mixing of the Jiulong River outflow with tidal water from the Xiamen Sea harbor. The salinity of the JRE was 25.50, with a pH value of 7.88 and a sea surface temperature of approximately 23.3 °C. The abundance of virioplankton in the JRE, as estimated by FCM, was 3.38 × 107 particles/mL, which was more than an order of magnitude higher than that of bacterioplankton (2.21 × 106 particles/mL). The abundances of picoeukaryotes and Synechococcus were estimated as 2.09 × 104 cells/mL and 7.03 × 104 cells/mL, respectively (Figure 1).
Viruses 08 00035 g001 1024
Figure 1. Cytometric plot graph of bacterioplankton (A), virioplankton (B) and picophyplankton (C) populations of sample from the JRE.
Figure 1. Cytometric plot graph of bacterioplankton (A), virioplankton (B) and picophyplankton (C) populations of sample from the JRE.
Viruses 08 00035 g001

3.2. Overview of the JRE Virome

Pyrosequencing of the JRE virome resulted in 527,507 raw reads, with an average length of 495 bp. A total of 498,957 reads with a 505-bp average length were generated after filtering and trimming. The average G + C content was 42.6% for clean data. Based on 90% sequence identity, the rarefaction curve generated by Metavir showed approximately 260,000 sequencing clusters. The number of clusters increased to 325,000 at 98% similarity (Figure 2). According to the current phage taxonomy, viral isolates are grouped into the same species at nucleotide identity levels of 90%–95% [43]. Given an average phage genome length of 50,000 bp and an average read length of 500 bp [12,44], the viral richness of JRE metavirome was estimated to be between 2600 and 3250 (# clusters/(genome length/read length)), which is almost twice higher than that of a recently published hypolithic virome [43].
Viruses 08 00035 g002 1024
Figure 2. Rarefaction curve of the JRE virome generated by Metavir. Clustering were set at 90% and 98% similarities.
Figure 2. Rarefaction curve of the JRE virome generated by Metavir. Clustering were set at 90% and 98% similarities.
Viruses 08 00035 g002

3.3. Taxonomic Composition

A total of 498,957 reads with a 505-bp average length were used for analysis of the JRE viral taxonomic composition. According to the MG-RAST annotation, 40.7% of reads had homology to known sequences. Reads were also assessed by Metavir for in-depth analysis of viral taxonomic composition using the GAAS tool, which normalizes the composition plot according to genome lengths [45], with an E-value threshold of 10−5. As with most viromes, the vast majority of reads from the JRE virome did not show any significant sequence similarity to the database. Only 21.42% of the reads produced significant hits. Most (97.60%) of the sequences with significant hits were recognized as belonging to dsDNA viruses (Figure 3), while similarities to single-stranded DNA phages and RNA viruses were also observed. The majority (84.81%) of sequences showing similarity to dsDNA viruses belonged to the order Caudovirales, with Podoviridae being the most abundant family (44.88%), followed by Siphoviridae (32.98%) and Myoviridae (17.32%). Furthermore, 4.82% of the sequences shared identity with “unclassified Caudovirales” sequences. Other dsDNA viruses identified at an abundance of 0.05%–1% included Phycodnaviridae, Ascoviridae and Mimiviridae. The rank-abundance curve for the JRE virome was characterized by a small number of abundant genomes and a large number of genomes with a relative abundance <1%. Among the 13 genomes with a relative abundance >1%, the five most abundant were identified as Pelagibacter phage HTVC010P (5.08%), Puniceispirillum phage HMO-2011 (3.91%) and Thalassomonas phage BA3 (2.07%), three Podoviridae viruses, and two members of unclassified dsDNA phages, Pseudomonas phage PA11 (2.64%) and Celeribacter phage P12053L (2.12%). Phages infecting enterobacteria, Staphylococcus species, and other pathogens infecting mammalian hosts were detected at >0.5% abundance. Two eukaryotic viruses, Tetraselmis viridis virus SI1 and Dunaliella viridis virus SI2 were also identified at relative abundances of 0.59% and 0.50%, respectively.
Viruses 08 00035 g003 1024
Figure 3. Taxonomic composition of the JRE virome, based on BLAST comparison with the NCBI Refseq complete viral genomes protein sequences database (release of 5 January 2015) by Metavir with an E-value threshold of 10−5.
Figure 3. Taxonomic composition of the JRE virome, based on BLAST comparison with the NCBI Refseq complete viral genomes protein sequences database (release of 5 January 2015) by Metavir with an E-value threshold of 10−5.
Viruses 08 00035 g003

3.4. Functional Analysis

Metabolic subsystems in the JRE virome were annotated using MG-RAST. The data were compared to subsystems using a maximum E-value of 1e−5 and a minimum identity of 60%. The database search resulted in 49,387 functional hits, the majority (31.26%) of which belonged to the subsystem “Phages, prophages, transposable elements, plasmids” with phage structural, integration/excision and DNA metabolism-related proteins. Other highly represented (>5%) functional categories included “Clustering-based subsystem”, “DNA metabolism”, “Nucleotides and nucleosides” and “Protein metabolism” (Figure 4).
Viruses 08 00035 g004 1024
Figure 4. Comparison of metabolic subsystems of the JRE virome with other aquatic viromes. The coding sequences were compared with SEED database using subsystems in MG-RAST. Metabolic subsystems from the JRE virome are colored by red. Green represents functional categories from temperate lakes (stations CMD, PM, and MO) [14]. Blue represents functional categories from the Pacific Ocean (stations LA26S, LF26S, LJ12S, LJ26S, M2MS, and M4OS) [31], and the Indian Ocean (stations GS122, GS108, GS117, and GS112) [46].
Figure 4. Comparison of metabolic subsystems of the JRE virome with other aquatic viromes. The coding sequences were compared with SEED database using subsystems in MG-RAST. Metabolic subsystems from the JRE virome are colored by red. Green represents functional categories from temperate lakes (stations CMD, PM, and MO) [14]. Blue represents functional categories from the Pacific Ocean (stations LA26S, LF26S, LJ12S, LJ26S, M2MS, and M4OS) [31], and the Indian Ocean (stations GS122, GS108, GS117, and GS112) [46].
Viruses 08 00035 g004
The functional categories of the JRE virome were then compared with those of other aquatic viromes, including temperate lakes (stations CMD, PM, and MO) [14], the Pacific Ocean (stations LA26S, LF26S, LJ12S, LJ26S, M2MS, and M4OS) [31], and the Indian Ocean (stations GS122, GS108, GS117, GS112) [46] on the MG-RAST server. Based on the comparison, we determined that subsystems related to “Nucleotides and nucleosides” and “Cofactors, vitamins, prosthetic groups, and pigments” were significantly over-represented in lake viromes compared with the JRE virome and other marine counterparts (Figure 4). However, some functions, including “Protein metabolism”, “Miscellaneous”, and “Carbohydrate” were more highly represented in the JRE virome and viromes from other marine environments.

3.5. Contig Analysis

Assembly using Newbler yielded 13,278 contigs, with an average length of 1126 bp accounting for a total of 15,729,985 bp for JRE virome. Because of the input restriction that only sequences longer than 300 bp are considered for projects of type contigs by Metavir, 10,246 contigs longer than 300 bp were selected for analysis. In assessing the composition of the assembled contigs, 39.15% of the contigs showed similarity to known sequences. A total of 31,476 genes were predicted, with a high proportion (78.4%) of predicted ORFs showing no significant homology to sequences present in the databases.
Of the 8367 contigs >500 bp, 28 contigs >20 kb were assembled, with the largest being 62 kb. Only two contigs produced significant hits (E-value < 10−3) to complete viral genomes in the publicly available database. The two contigs were 46,070 bp (contig 3) and 36,247 bp (contig 9) long, with the overall G + C content 44.70% and 40.35%, respectively. Sixty-two ORFs were predicted by GeneMarkS within contig 3, of which 15 showed similarity to Celeribacter phage P12053L (Figure 5), a marine member of the unclassified Caudovirales group [47], indicating that this contig almost certainly belonged to the order Caudovirales. The majority of genes in contig 3 for which a function could be assigned (11 of 17) encoded proteins related to DNA metabolism and replication, including endonuclease, DNA primase, DNA polymerase, single-strand DNA binding protein, thymidylate synthase and ribonucleotide reductase. Four ORFs were predicted to encode proteins involved in the structure and assembly of virions. The translated product of an ORF identified as the terminase large subunit gene showed 80% amino acid similarity to the corresponding protein of Celeribacter phage P12053L. Lysozyme genes were also predicted within contig 3. The genome maps showed that contig 9 present similar overall genome organization to that of phage BA3 (Figure 5), which lyses the coral pathogen Thalassomonas loyana. In total, 61.2% of ORFs identified (30 of 49) from contig 9 showed similarity to those of phage BA3, with the corresponding proteins showing 28%–100% amino acid identity.
Viruses 08 00035 g005 1024
Figure 5. Comparisons of contigs 3 and 9 with Celeribacter phage P12053L and Thalosomonas phage BA3, respectively. ORFs are depicted by leftward or rightward oriented arrows according to the direction of transcription. Colored arrows denote putative functions assigned by Blastp similarity. Homologous ORFs are connected by shadowing.
Figure 5. Comparisons of contigs 3 and 9 with Celeribacter phage P12053L and Thalosomonas phage BA3, respectively. ORFs are depicted by leftward or rightward oriented arrows according to the direction of transcription. Colored arrows denote putative functions assigned by Blastp similarity. Homologous ORFs are connected by shadowing.
Viruses 08 00035 g005

3.6. Comparison with Other Viromes

BLAST comparisons of the JRE virome with previously published aquatic viromes were carried out to examine the similarities between different aquatic viral communities. The results revealed that aquatic viromes could be classified into two representative groups, freshwater viromes and marine viromes (Figure 6A). Generally, the JRE virome clustered with the marine viromes, and was closely related to the Pacific Ocean surface samples from previous study [31]. To obtain a deeper comparison of related datasets, the JRE virome was further compared with eight PCSVs. Comparison of the dominant members of the JRE virome with the PCSVs showed that most viral genotypes were shared, but that there were large differences in the relative abundance of particular species between the samples. In addition, two Pseudomonas phages infecting P. aeruginosa, an opportunistic human pathogen, were exclusively found in the JRE sample. The results of the cross-tblastx analysis indicated that the JRE virome showed a coarse-grain genetic overlap of 14.14% ± 1.68% with the PCSVs, with the highest relative overlap observed with the virome of sample STCS, followed by GDS and GFS (Figure 6B). Samples SFSS, STCS, SFDS, and SFCS were derived from the same coastal surface water from Scripps Pier, San Diego, CA, USA [31]. Sites GDS and GFS, corresponding to the Great Barrier Reef, Australia, were under the influence of the Tully River, which appeared to explain the relatively high genetic overlap with the JRE virome.
Viruses 08 00035 g006 1024
Figure 6. (A) MDS plots of freshwater (green) and marine (cyan) publicly available viromes with the JRE virome (red) using Metavir; (B) Coarse-grain genetic overlap between the JRE virome and PCSVs. The heat-plot represents the results of cross-tblastx comparisons between samples. Values (% of total reads with hits) with E-value ≤ 10−3 are labeled within the boxes. Horizontal labels represent viromes acting as query, and vertical labels represent viromes acting as reference.
Figure 6. (A) MDS plots of freshwater (green) and marine (cyan) publicly available viromes with the JRE virome (red) using Metavir; (B) Coarse-grain genetic overlap between the JRE virome and PCSVs. The heat-plot represents the results of cross-tblastx comparisons between samples. Values (% of total reads with hits) with E-value ≤ 10−3 are labeled within the boxes. Horizontal labels represent viromes acting as query, and vertical labels represent viromes acting as reference.
Viruses 08 00035 g006

4. Discussion

4.1. The JRE Virome

Consistent with most other viromes examined [12,14,43,46,48], the investigation of dsDNA viral diversity in the JRE showed that the majority of sequence reads were classified as unknown, with only 21.42% having known virus matches. This is in line with identification rates for the value of the other aquatic dsDNA metaviromes [14,31,46,49], which range from 6.75% to 28.67%. These findings support the hypothesis that much of the global viral diversity remains uncharacterized. Analysis of taxonomic composition for JRE virome showed most viruses belonged to the order Caudovirales, with Podoviridae representing the major fraction, followed by Siphoviridae and Myoviridae. In addition, since the rarefaction curve of the JRE virome did not tend to approach a horizontal asymptote, the viral richness of the JRE estuary was potentially higher than that revealed here.
Based on clustering analysis in which viromes were compared to each other via Metavir through BLASTx, viral communities appeared to group according to habitat, with viromes from freshwater and marine samples forming distinct clusters. Similar results were also obtained from previous studies [14,16]. This reflected unique virome signatures contained in contrasting environments (seawater vs. freshwater). Therefore, viral communities in estuarine environments are expected to be more diverse due to the mixing of seawater and freshwater, the resuspension of sediments and the interference of human activities. Previous virome investigation of several environmental transects already suggested that surface viral richness decreased with coast-open ocean gradients [31], which mirrored the distribution pattern of their hosts (mainly bacteria).

4.2. Impacts of Marine Water on the Composition of the JRE Virome

Clustering analysis showed that the JRE virome clustered with other marine viromes and was closely related to the PCSVs, suggesting that it was mainly affected by marine waters. The JRE virome also shared a significant genetic overlap with viromes from costal and estuarine areas. Despite an overall similarity in the viral genotypes of the viromes examined, there were large differences in the relative abundance of particular families between the JRE virome and the PCSVs. This finding is consistent with the results of a previous study, which indicated that the change in abundance of the most common members is what sets different assemblages apart, rather than the exclusion of different viral genomes [13]. This idea forms the basis of a theoretical prediction made for marine virus communities, termed the “Seed-Bank model” [50].
Obvious marine signatures were observed following deeper inspection of phage genotypes. Viruses known to infect typical marine bacteria, such as Proteobacteria and Bacteroidetes, were widely present in the JRE sample. For example, the most abundant virus identified in the JRE virome was Pelagibacter phage (Pelagiphage) HTVC010P, which infect “Candidatus Pelagibacter ubique” of the SAR11 clade, the most abundant bacterium in surface seawater around the world [51]. In a more extensive metagenomic investigation including samples from coastal and open ocean areas, Pelagiphage HTVC010P was proposed as one of the most abundant virus subfamilies in the biosphere, in which 38.8% of successfully assigned reads were assigned to HTVC010P [52]. The second-most abundant virus genome was phage HMO-2011, a phage infecting “Candidatus Puniceispirillum marinum” of the SAR116 clade. Phage HMO-2011 was firstly isolated from the East Sea and was revealed to be widely distributed in the ocean, accounting for 10.3%–25.3% of reads in seven viromes from the coastal and open area of the Pacific and Indian Oceans [53]. Our study demonstrates that Pelagiphage HTVC010P and phage HMO-2011 may be abundant in estuarine environments as well.
In general, the prevalence of specific viral populations follows that of their hosts. The composition of the metavirome sequenced in the current study is consistent with the results of 16S rRNA gene analysis performed on samples from the same habitat [32]. The most dominant bacterial phylum was Proteobacteria, followed by Bacteroidetes, which is consistent with our assignment of the virome. At the operational taxonomic unit level, the most abundant bacterial sequences were affiliated with the SAR11 clade. Therefore, the prevalence of Pelagiphage HTVC010P and phage HMO-2011 sequences in the JRE virome corroborates the high abundance of their marine bacterial hosts.
Further evidence of the impact of marine water on the composition of the JRE virome was obtained from contig analysis. Two of the assembled contigs showed significant similarity to the complete genomes of Celeribacter phage P12053L and Thalosomonas phage BA3, respectively. Celeribacter phage P12053L, a lytic bacteriophage infecting bacteria of the the Roseobacter clade, was isolated from seawater collected off the coast of the Yellow Sea [47]. The Yellow Sea is close to the environment sampled in this study, confirming Celeribacter phage P12053L as a significant contributor to viral communities in coastal surface waters, and the high abundance of Roseobacteria in the ocean. Thalosomonas phage BA3, which lyses the coral pathogen Thalassomonas loyana, has previously been isolated from the Red Sea [54]. The presence of these two contigs in the JRE virome suggests the prevalence of Celeribacter phage P12053L and Thalosomonas phage BA3 in the JRE, and supports the importance of marine microbial populations in shaping the estuarine virome.

4.3. Impacts of River Outflow on the Composition of the JRE Virome

In general, relatively few impacts of river outflow were found in JRE virome, which agreed with the general environmental characteristics of our sampling station. The salinity and other environmental factors indicated that the sampling site was more similar to typical marine environments. Therefore, it is not surprising that only a few functional categories were consistent with lake viromes, showing higher relative abundance than in marine viromes. However, two P. aeruginosa phages were exclusively found in the JRE sample when compared with the PCSVs. P. aeruginosa is an opportunistic pathogen that is important in the etiology of many infectious diseases of humans, and is widely distributed in non-marine environments (e.g., soil, ground water, sewage, the mammalian gut). The apparent distribution of P. aeruginosa in the marine system has been restricted to river outfalls and shorelines [55,56,57]. Therefore, we deduced that Pseudomonas phages in the JRE virome probably originated from Jiulong River discharge. Phages infecting enterobacteria and other mammalian pathogens were also found in the JRE virome. Similar to their pathogen hosts, the occurrence of these non-marine phages in an estuary virome likely results from freshwater transfer from the river, which is a site of intense human activities.
The duality of the JRE virome can be explained by idea that the JRE is a sort of mixing chamber that receives Jiulong River outflow with the dash of marine water. Similar patterns for their hosts were frequently observed in estuarine area [32,58,59]. Previous research showed that the dynamics of viral abundance in the Yangtze River estuarine area were characterized by distinct seasonal and spatial variations caused by natural forces and anthropogenic impacts [58]. In addition, two eukaryotic viruses, SI1 and SI2, were identified in the current study. These viruses infect the eukaryotic microalgae Tetraselmis viridis Rouch and Dunaliella viridi. These two microalgae are commonly found in a wide range of habitats with varying degrees of salinity: from fresh water to salt lakes containing sodium chloride at concentrations up to saturation. The presence of these viruses in the JRE virome might reflect the adaptation of the eukaryotic microbial population to dynamic estuary environments with large temporal and spatial fluctuations of environmental parameters such as salinity.

4.4. Methodological Considerations and Limitations

Metagenomic analysis has become an invaluable tool to explore total viral taxonomic composition and diversity from environmental samples. However, the virome procedure from sample collection, sequencing preparation to bioinformatics analysis is experimental and informatics-challenging, which may induce biases in estimating viral diversity.
Due to the small individual size and biomass, virus enrichment steps are necessary to obtain sufficient experimental materials [60]. As with the majority of viromes published, a 0.22-μm TFF step was used in this study to remove cells. Large viruses such as Phycodnaviridae were likely removed during this process, resulting in the low representation in the virome. Recently, a new chemistry-based concentration method, FeCl3 precipitation, has emerged and showed nearly complete viral recovery, which may be widely applied for viral concentration in the future [61].
Although the prevalence of Caudovirales sequences was observed in most aquatic viromes, previous research on global viral morphology through the Tara Oceans Expedition has revealed that non-tailed viruses numerically dominated the upper oceans [62]. Non-tailed viruses are largely RNA and ssDNA viruses. It is possible that the relative lack of non-tailed viruses may be the results from the way of sample preparation here, where dsDNA viruses being selected in the CsCl ultracentrifugation and library preparation for 454-based sequencing [63]. Besides, it can also be inferred as the consequence of limited non-tailed viruses in the database [62].
Additionally, the paucity of viral reference databases affects the ability to identify viral sequences in the virome. Albeit the virome analyzed in this study was revealed to be composed primarily of dsDNA viruses, the BLAST results were limited to the small fraction of sequences with matches in reference database. The unknowns could represent a great treasure or a valuable blueprint for the discovery of novel viruses as the development of bioinformatics and the continuous supplement of viral sequences. Finally, further efforts to sequence more deeply will improve our population level analysis with better assembly.

5. Conclusions

The present study of the JRE virome has given us a glimpse into the complexity of a virus population from a typical subtropical estuary environment. Overall, the JRE shared a similar taxonomic composition, and presented some genetic overlap with viromes from coastal environments. However, it is distinct from other publicly available datasets, representing a degree of novelty. The endemic viruses in the JRE virome were shaped by local environmental conditions, including the river, coastal ocean, and anthropogenic activities. This study provides a baseline for subsequent systemic investigations.

Acknowledgments

We thank Yu Wang of the Institute of Marine Microbes and Ecospheres for his useful suggestions and help. This work was supported by the National Natural Science Foundation of China (41522603, 41376132), the 973 project (2013CB955700), and the State Oceanic Administration (SOA) Project (GASI-03-01-02-05). Rui Zhang was partially supported by the Fundamental Research Funds for the Central Universities (20720150156).

Author Contributions

Rui Zhang and Nianzhi Jiao conceived and designed the experiments; Lanlan Cai performed the experiments; Lanlan Cai, Ying He and Xiaoyuan Feng analyzed the data; Ying He contributed analysis tools; Lanlan Cai, Rui Zhang and Nianzhi Jiao wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Suttle, C.A. Marine viruses—Major players in the global ecosystem. Nat. Rev. Microbiol. 2007, 5, 801–812. [Google Scholar] [CrossRef] [PubMed]
  2. Jacquet, S.; Miki, T.; Noble, R.; Peduzzi, P.; Wilhelm, S. Viruses in aquatic ecosystems: Important advancements of the last 20 years and prospects for the future in the field of microbial oceanography and limnology. Adv. Oceanogr. Limnol. 2010, 1, 97–141. [Google Scholar] [CrossRef]
  3. Thingstad, T.F. Elements of a theory for the mechanisms controlling abundance, diversity, and biogeochemical role of lytic bacterial viruses in aquatic systems. Limnol. Oceanogr. 2000, 45, 1320–1328. [Google Scholar] [CrossRef]
  4. Rohwer, F.; Thurber, R.V. Viruses manipulate the marine environment. Nature 2009, 459, 207–212. [Google Scholar] [CrossRef] [PubMed]
  5. Breitbart, M. Marine viruses: Truth or dare. Ann. Rev. Mar. Sci. 2012, 4, 425–448. [Google Scholar] [CrossRef]
  6. Brum, J.R.; Sullivan, M.B. Rising to the challenge: Accelerated pace of discovery transforms marine virology. Nat. Rev. Microbiol. 2015. [Google Scholar] [CrossRef]
  7. Lindell, D.; Jaffe, J.D.; Johnson, Z.I.; Church, G.M.; Chisholm, S.W. Photosynthesis genes in marine viruses yield proteins during host infection. Nature 2005, 438, 86–89. [Google Scholar] [PubMed]
  8. Thompson, L.R.; Zeng, Q.; Kelly, L.; Huang, K.H.; Singer, A.U.; Stubbe, J.; Chisholm, S.W. Phage auxiliary metabolic genes and the redirection of cyanobacterial host carbon metabolism. Proc. Natl. Acad. Sci. USA 2011, 108, E757–E764. [Google Scholar] [CrossRef] [PubMed]
  9. Weinbauer, M.G. Ecology of prokaryotic viruses. FEMS Microbiol. Rev. 2004, 28, 127–181. [Google Scholar] [CrossRef] [PubMed]
  10. Zhang, R.; Wei, W.; Cai, L. The fate and biogeochemical cycling of viral elements. Nat. Rev. Microbiol. 2014, 12, 850–851. [Google Scholar] [CrossRef] [PubMed]
  11. Jiao, N.; Herndl, G.J.; Hansell, D.A.; Benner, R.; Kattner, G.; Wilhelm, S.W.; Kirchman, D.L.; Weinbauer, M.G.; Luo, T.; Chen, F.; et al. Microbial production of recalcitrant dissolved organic matter: Long-term carbon storage in the global ocean. Nat. Rev. Microbiol. 2010, 8, 593–599. [Google Scholar] [CrossRef]
  12. 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]
  13. 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, e368. [Google Scholar] [CrossRef] [PubMed]
  14. Green, J.C.; Rahman, F.; Saxton, M.A.; Williamson, K.E. Metagenomic assessment of viral diversity in lake matoaka, a temperate, eutrophic freshwater lake in southeastern Virginia, USA. Aquat. Microbial. Ecol. 2015, 75, 117–128. [Google Scholar] [CrossRef]
  15. Tseng, C.H.; Chiang, P.W.; Shiah, F.K.; Chen, Y.L.; Liou, J.R.; Hsu, T.C.; Maheswararajah, S.; Saeed, I.; Halgamuge, S.; Tang, S.L. Microbial and viral metagenomes of a subtropical freshwater reservoir subject to climatic disturbances. ISME J. 2013, 7, 2374–2386. [Google Scholar] [CrossRef] [PubMed]
  16. 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]
  17. Aguirre de Carcer, D.; Lopez-Bueno, A.; Pearce, D.A.; Alcami, A. Biodiversity and distribution of polar freshwater DNA viruses. Sci. Adv. 2015, 1. [Google Scholar] [CrossRef] [PubMed]
  18. Coetzee, B.; Freeborough, M.J.; Maree, H.J.; Celton, J.M.; Rees, D.J.; Burger, J.T. Deep sequencing analysis of viruses infecting grapevines: Virome of a vineyard. Virology 2010, 400, 157–163. [Google Scholar] [CrossRef] [PubMed]
  19. Breitbart, M.; Hewson, I.; Felts, B.; Mahaffy, J.M.; Nulton, J.; Salamon, P.; Rohwer, F. Metagenomic analyses of an uncultured viral community from human feces. J. Bacteriol. 2003, 185, 6220–6223. [Google Scholar] [CrossRef] [PubMed]
  20. Jones, M.S.; Kapoor, A.; Lukashov, V.V.; Simmonds, P.; Hecht, F.; Delwart, E. New DNA viruses identified in patients with acute viral infection syndrome. J. Virol. 2005, 79, 8230–8236. [Google Scholar] [CrossRef] [PubMed]
  21. Ng, T.F.; Suedmeyer, W.K.; Wheeler, E.; Gulland, F.; Breitbart, M. Novel anellovirus discovered from a mortality event of captive California sea lions. J. Gen. Virol 2009, 90, 1256–1261. [Google Scholar] [CrossRef] [PubMed]
  22. Schoenfeld, T.; Patterson, M.; Richardson, P.M.; Wommack, K.E.; Young, M.; Mead, D. Assembly of viral metagenomes from Yellowstone hot springs. Appl. Environ. Microbiol. 2008, 74, 4164–4174. [Google Scholar] [CrossRef]
  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. Edwards, R.A.; Rohwer, F. Viral metagenomics. Nat. Rev. Microbiol. 2005, 3, 504–510. [Google Scholar] [CrossRef] [PubMed]
  25. Kim, Y.; Aw, T.G.; Teal, T.K.; Rose, J.B. Metagenomic investigation of viral communities in ballast water. Environ. Sci. Technol. 2015, 49, 8396–8407. [Google Scholar] [CrossRef] [PubMed]
  26. Martinez Martinez, J.; Swan, B.K.; Wilson, W.H. Marine viruses, a genetic reservoir revealed by targeted viromics. ISME J. 2014, 8, 1079–1088. [Google Scholar] [CrossRef] [PubMed]
  27. Wommack, K.E.; Bhavsar, J.; Ravel, J. Metagenomics: Read length matters. Appl. Environ. Microbiol. 2008, 74, 1453–1463. [Google Scholar] [CrossRef]
  28. Brum, J.R.; Ignacio-Espinoza, J.C.; Roux, S.; Doulcier, G.; Acinas, S.G.; Alberti, A.; Chaffron, S.; Cruaud, C.; de Vargas, C.; Gasol, J.M.; et al. Patterns and ecological drivers of ocean viral communities. Science 2015, 348. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Helton, R.R.; Wommack, K.E. Seasonal dynamics and metagenomic characterization of estuarine viriobenthos assemblages by randomly amplified polymorphic DNA PCR. Appl. Environ. Microbiol. 2009, 75, 2259–2265. [Google Scholar] [CrossRef] [PubMed]
  30. Fortunato, C.S.; Herfort, L.; Zuber, P.; Baptista, A.M.; Crump, B.C. Spatial variability overwhelms seasonal patterns in bacterioplankton communities across a river to ocean gradient. ISME J. 2012, 6, 554–563. [Google Scholar] [CrossRef] [PubMed]
  31. 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]
  32. Shan, D.; Wei, G.; Li, M.; Wang, W.; Li, X.; Gao, Z.; Shao, Z. Distribution and diversity of bacterioplankton communities in subtropical seawater around Xiamen island, China. Microbiol. Res. 2015, 175, 16–23. [Google Scholar] [CrossRef] [PubMed]
  33. Thurber, R.V.; Haynes, M.; Breitbart, M.; Wegley, L.; Rohwer, F. Laboratory procedures to generate viral metagenomes. Nat. Protoc. 2009, 4, 470–483. [Google Scholar] [CrossRef] [PubMed]
  34. Marie, D.; Brussaard, C.P.D.; Thyrhaug, R.; Bratbak, G.; Vaulot, D. Enumeration of marine viruses in culture and natural samples by flow cytometry. Appl. Environ. Microbiol. 1999, 65, 45–52. [Google Scholar] [PubMed]
  35. Brussaard, C.P. Optimization of procedures for counting viruses by flow cytometry. Appl. Environ. Microbiol. 2004, 70, 1506–1513. [Google Scholar] [CrossRef] [PubMed]
  36. Jiao, N.Z.; Yang, Y.H.; Koshikawa, H.; Watanabe, M. Influence of hydrographie conditions on picoplankton distribution in the East China Sea. Aquat. Microbial. Ecol. 2002, 30, 37–48. [Google Scholar] [CrossRef]
  37. EXPO32; v1.2; Beckman Coulter: Miami, FL, USA, 2001.
  38. Niu, B.; Fu, L.; Sun, S.; Li, W. Artificial and natural duplicates in pyrosequencing reads of metagenomic data. BMC Bioinform. 2010, 11. [Google Scholar] [CrossRef] [PubMed]
  39. Mothur. Available online: http://www.mothur.org/ (accessed on 30 January 2016).
  40. 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]
  41. 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 Bioinform. 2008, 9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Besemer, J.; Lomsadze, A.; Borodovsky, M. Genemarks: A self-training method for prediction of gene starts in microbial genomes. Implications for finding sequence motifs in regulatory regions. Nucleic Acids Res. 2001, 29, 2607–2618. [Google Scholar] [CrossRef]
  43. Adriaenssens, E.M.; van Zyl, L.; de Maayer, P.; Rubagotti, E.; Rybicki, E.; Tuffin, M.; Cowan, D.A. Metagenomic analysis of the viral community in Namib Desert hypoliths. Environ. Microbiol. 2015, 17, 480–495. [Google Scholar] [CrossRef] [PubMed]
  44. Steward, G.F.; Montiel, J.L.; Azam, F. Genome size distributions indicate variability and similarities among marine viral assemblages from diverse environments. Limnol. Oceanogr. 2000, 45, 1697–1706. [Google Scholar] [CrossRef]
  45. Angly, F.E.; Willner, D.; Prieto-Davo, A.; Edwards, R.A.; Schmieder, R.; Vega-Thurber, R.; Antonopoulos, D.A.; Barott, K.; Cottrell, M.T.; Desnues, C.; et al. The GAAS metagenomic tool and its estimations of viral and microbial average genome size in four major biomes. PLoS Comput. Biol. 2009, 5, e1000593. [Google Scholar] [CrossRef] [PubMed]
  46. 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] [PubMed]
  47. Kang, I.; Jang, H.; Oh, H.M.; Cho, J.C. Complete genome sequence of Celeribacter bacteriophage P12053L. J. Virol. 2012, 86, 8339–8340. [Google Scholar] [CrossRef]
  48. Winter, C.; Garcia, J.A.; Weinbauer, M.G.; DuBow, M.S.; Herndl, G.J. Comparison of deep-water viromes from the Atlantic Ocean and the Mediterranean Sea. PLoS ONE 2014, 9, e100600. [Google Scholar]
  49. Chow, C.E.; Winget, D.M.; White, R.A., 3rd; Hallam, S.J.; Suttle, C.A. Combining genomic sequencing methods to explore viral diversity and reveal potential virus-host interactions. Front. Microbiol. 2015, 6. [Google Scholar] [CrossRef] [PubMed]
  50. Breitbart, M.; Rohwer, F. Here a virus, there a virus, everywhere the same virus? Trends Microbiol. 2005, 13, 278–284. [Google Scholar] [CrossRef] [PubMed]
  51. Morris, R.M.; Rappe, M.S.; Connon, S.A.; Vergin, K.L.; Siebold, W.A.; Carlson, C.A.; Giovannoni, S.J. SAR11 clade dominates ocean surface bacterioplankton communities. Nature 2002, 420, 806–810. [Google Scholar] [CrossRef] [PubMed]
  52. Zhao, Y.; Temperton, B.; Thrash, J.C.; Schwalbach, M.S.; Vergin, K.L.; Landry, Z.C.; Ellisman, M.; Deerinck, T.; Sullivan, M.B.; Giovannoni, S.J. Abundant SAR11 viruses in the ocean. Nature 2013, 494, 357–360. [Google Scholar] [CrossRef] [PubMed]
  53. 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] [PubMed]
  54. Efrony, R.; Loya, Y.; Bacharach, E.; Rosenberg, E. Phage therapy of coral disease. Coral Reefs 2007, 26, 7–13. [Google Scholar] [CrossRef]
  55. Velammal, A.; Aiyamperumal, B.; Venugopalan, V.K.; Ajmalkhan, S. Distribution of Pseudomonas aeruginosa in pondicherry coastal environs. Indian J. Mar. Sci. 1994, 23, 239–241. [Google Scholar]
  56. Mates, A. The significance of testing for Pseudomonas aeruginosa in recreational seawater beaches. Microbios 1992, 71, 89–93. [Google Scholar] [PubMed]
  57. Kimata, N.; Nishino, T.; Suzuki, S.; Kogure, K. Pseudomonas aeruginosa isolated from marine environments in Tokyo Bay. Microb. Ecol. 2004, 47, 41–47. [Google Scholar] [CrossRef]
  58. Jiao, N.Z.; Zhao, Y.L.; Luo, T.W.; Wang, X.L. Natural and anthropogenic forcing on the dynamics of virloplankton in the Yangtze river estuary. J. Mar. Biol. Assoc. UK 2006, 86, 543–550. [Google Scholar] [CrossRef]
  59. Liu, L.; Yang, J.; Zhang, Y. Genetic diversity patterns of microbial communities in a subtropical riverine ecosystem (Jiulong River, southeast China). Hydrobiologia 2011, 678, 113–125. [Google Scholar] [CrossRef]
  60. Cai, L.; Yang, Y.; Jiao, N.; Zhang, R. Evaluation of tangential flow filtration for the concentration and separation of bacteria and viruses in contrasting marine environments. PLoS ONE 2015, 10, e0136741. [Google Scholar] [CrossRef]
  61. John, S.G.; Mendez, C.B.; Deng, L.; Poulos, B.; Kauffman, A.K.; Kern, S.; Brum, J.; Polz, M.F.; Boyle, E.A.; Sullivan, M.B. A simple and efficient method for concentration of ocean viruses by chemical flocculation. Environ. Microbiol. Rep. 2011, 3, 195–202. [Google Scholar] [CrossRef]
  62. 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]
  63. Bellas, C.M.; Anesio, A.M.; Barker, G. Analysis of virus genomes from glacial environments reveals novel virus groups with unusual host interactions. Front. Microbiol. 2015, 6. [Google Scholar] [CrossRef] [PubMed]

Share and Cite

MDPI and ACS Style

Cai, L.; Zhang, R.; He, Y.; Feng, X.; Jiao, N. Metagenomic Analysis of Virioplankton of the Subtropical Jiulong River Estuary, China. Viruses 2016, 8, 35. https://doi.org/10.3390/v8020035

AMA Style

Cai L, Zhang R, He Y, Feng X, Jiao N. Metagenomic Analysis of Virioplankton of the Subtropical Jiulong River Estuary, China. Viruses. 2016; 8(2):35. https://doi.org/10.3390/v8020035

Chicago/Turabian Style

Cai, Lanlan, Rui Zhang, Ying He, Xiaoyuan Feng, and Nianzhi Jiao. 2016. "Metagenomic Analysis of Virioplankton of the Subtropical Jiulong River Estuary, China" Viruses 8, no. 2: 35. https://doi.org/10.3390/v8020035

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