Quantitative Viral Community DNA Analysis Reveals the Dominance of Single-Stranded DNA Viruses in Offshore Upper Bathyal Sediment from Tohoku, Japan

Previous studies on marine environmental virology have primarily focused on double-stranded DNA (dsDNA) viruses; however, it has recently been suggested that single-stranded DNA (ssDNA) viruses are more abundant in marine ecosystems. In this study, we performed a quantitative viral community DNA analysis to estimate the relative abundance and composition of both ssDNA and dsDNA viruses in offshore upper bathyal sediment from Tohoku, Japan (water depth = 500 m). The estimated dsDNA viral abundance ranged from 3 × 10⁶ to 5 × 10⁶ genome copies per cm³ sediment, showing values similar to the range of fluorescence-based direct virus counts. In contrast, the estimated ssDNA viral abundance ranged from 1 × 10⁸ to 3 × 10⁹ genome copies per cm³ sediment, thus providing an estimation that the ssDNA viral populations represent 96.3–99.8% of the benthic total DNA viral assemblages. In the ssDNA viral metagenome, most of the identified viral sequences were associated with ssDNA viral families such as Circoviridae and Microviridae. The principle components analysis of the ssDNA viral sequence components from the sedimentary ssDNA viral metagenomic libraries found that the different depth viral communities at the study site all exhibited similar profiles compared with deep-sea sediment ones at other reference sites. Our results suggested that deep-sea benthic ssDNA viruses have been significantly underestimated by conventional direct virus counts and that their contributions to deep-sea benthic microbial mortality and geochemical cycles should be further addressed by such a new quantitative approach.

Introduction

Viruses are abundant components in aquatic ecosystems and play important roles in microbial mortality, gene transfer, and biogeochemical cycles (Fuhrman, 1999; Wommack and Colwell, 2000; Suttle, 2007). For many years, most viral ecology studies have focused on double-stranded DNA (dsDNA) viruses that infect bacteria (bacteriophages; also known as phages) (Mann et al., 2003; Weinbauer, 2004). Recent viral metagenome (virome) studies have greatly expanded knowledge of marine dsDNA viruses, particularly their hidden diversity and global distribution (Brum and Sullivan, 2015). In addition, virome studies have also suggested a previously unseen diversity and co-occurrence of single-stranded DNA (ssDNA) viruses in the marine virosphere (Angly et al., 2006; Tucker et al., 2011; Labonté and Suttle, 2013; Yoshida et al., 2013). However, these studies have always performed whole-genome amplification by using the multiple displacement amplification (MDA) technique (Edwards and Rohwer, 2005). It has been noted that this technique has a substantial risk of preferential amplification of circular ssDNA to dsDNA, and thus, the relative abundance of ssDNA virome sequences and associated viruses may have been overestimated (Lizardi et al., 1998; Dean et al., 2001). To solve this issue, a less biased approach for DNA viral community analysis is required to quantify the relative abundance of viruses with different nucleic acid types in various environments. Recently, a new virome library preparation protocol using a commercial kit for both ssDNA and dsDNA viral templates has been developed to estimate the respective viral abundance on the basis of the metagenomic read abundance of relevant virus-associated sequences (Roux et al., 2016). Otherwise, a mass apportioning approach to viral nucleic acids can be used to estimate the abundance ratio of RNA and DNA viruses (Steward et al., 2013; Miranda et al., 2016); however, this approach has not been tested for ssDNA and dsDNA viral quantification.

Interactions between viruses and their hosts (e.g., bacteria and archaea) play important ecological and biogeochemical roles in deep-sea surface sediments (down to 50 cm below the seafloor [cmbsf]) (Danovaro et al., 2008a,b, 2016). Yoshida et al. (2013) have reported that the surface sediments from bathyal to hadal depths in the northwest Pacific harbor highly diverse ssDNA viromes, including sequences related to Microviridae phages and Circoviridae viruses. However, the abundance and diversity of ssDNA viruses in marine benthic environments have not been quantitatively evaluated. Here, we present a mass apportioning-based virus community DNA analysis to investigate the relative abundance and composition of ssDNA and dsDNA viral assemblages in deep-sea surface sediments (down to 15 cmbsf) off Tohoku, Japan (water depth = 500 m), which is located in the eutrophic northwest Pacific region, by using the quantification of ssDNA and dsDNA viral genomes and subsequent metagenomic analysis for each of the viral genome assemblages. Because we first separated the total sediment viral DNA into ssDNA and dsDNA components and then prepared and analyzed the respective sequencing libraries, this analysis would make it possible to mitigate the preferential amplification of ssDNA in the total DNA metagenomic sample preparation during the MDA reaction and obtain more reliable data in the DNA viral composition analysis, as well as independently, provide unbiased estimates in terms of both the total viral abundance.

Materials and Methods

Sediment Samples

Short sediment cores (24 cm in length) were obtained using a multiple corer offshore from the Tohoku region in Japan (station N-500: 39°15′ N, 142°13′ E; water depth of 500 m) during the Japan Agency for Marine-Earth Science and Technology (JAMSTEC) MR12-E02 Leg3 cruise of the R/V Mirai in March 2012 (Nomaki et al., 2016a). The sediment cores were sliced at 1–5 cm intervals as previously described (Nomaki et al., 2016a), and the subsamples were stored at -80°C. Among them, three sections (0–2, 5–8, and 10–15 cmbsf) were used for the following analyses.

Direct Viral and Prokaryotic Cell Counts

To enumerate the viral and prokaryotic cell abundances, previously described direct cell and viral count methods were applied with minor modifications (Middelboe et al., 2011; Yanagawa et al., 2014). Briefly, approximately 2 cm³ of the frozen sediment was promptly suspended in 10 ml of modified SM buffer (50 mM Tris-HCl, pH 7.5; 10 mM MgSO₄; 100 mM NaCl) containing 2% (v/v) formaldehyde in a 50-ml centrifuge tube. The slurry was shaken with a ShakeMaster Auto (BioMedical Science, Tokyo, Japan) for 1 min at maximum speed and then sonicated for 1 min with an ultrasonic homogenizer (UH-50; SMT company, Tokyo, Japan) to detach viral particles from sediment matrices (Middelboe et al., 2011). After centrifugation for 5 min at 700 g, microbial cells in the supernatant were filtered with a 0.2-μm pore size black polycarbonate filter (Advantec, Tokyo, Japan). For viral counts, viral particles in the filtrate were further filtered with a 0.02-μm pore size Anodisc filter (GE Healthcare, Piscataway, NJ, United States). Direct cell and viral counts were obtained with 5× SYBR Gold (Invitrogen, Eugene, OR, United States; Chen et al., 2001) for 15 min at dark, by using an Olympus BX53 epifluorescence microscope at a magnification of 1,500× (Olympus, Tokyo, Japan). Each filter was mounted on a glass slide with a phosphate-buffered saline-glycerol mixture (1:1) containing 0.1% (w/v) p-phenylenediamine (Sigma–Aldrich, St. Louis, MO, United States) to minimize fading (Noble and Fuhrman, 1998). Typically, 10–20 microscopic fields were randomly selected, and >200 viral particles and >200 prokaryotic cells were counted for each filter.

Small Subunit (SSU) rRNA Gene Tag Sequencing

Environmental DNA was extracted from each sediment sample (approximately 5 g) by using a PowerMax Soil DNA Isolation kit (Mo Bio Lab, Carlsbad, CA, United States). SSU rRNA genes were amplified by LA Taq polymerase (Takara Bio, Kusatsu, Japan) using a universal primer set of U530F and U907R (Nunoura et al., 2012). The Illumina adaptor sequence (5′-ACA CTC TTT CCC TAC ACG ACG CTC TTC CGA TCT-3′; Illumina, San Diego, CA, United States) and Illumina Multiplexing PCR Primer 2.0 (5′-GTG ACT GGA GTT CAG ACG TGT GCT CTT CCG ATC T-3′; Illumina) were added at the 5′ ends of both primers (U530F and U907R). The PCR mixture contained LA Taq polymerase (final concentration of 0.1 U μl^-1), 1 × GC Buffer I (Takara Bio), dNTPs (final concentration of 0.25 mM), each primer (final concentration of 0.2 mM), and template DNA. The PCR conditions were as follows: 25 to 35 cycles of 20 s at 96°C, 45 s at 52°C, and 1 min at 72°C. Library preparation and sequencing with Illumina MiSeq are described in Hirai et al. (2017). For the PCR cycle number of each sample, we used the minimum cycle number that provided enough amplified products for pyrosequencing based on the preliminary PCR amplification experiments using the same templates.

The paired-end reads of 16S rRNA gene tags were merged together using the program PEAR v0.9.10 (Zhang et al., 2014). Low-quality reads that included 3% or more nucleotides with low-quality values (quality score <30) were removed using the FASTX-Toolkit v0.0.14¹. PCR primers were removed from the processed sequences using Cutadapt v1.10 (Martin, 2011). The resulting SSU rRNA gene tag libraries were analyzed with QIIME v1.9.1 pipelines (Caporaso et al., 2010). Briefly, the sequences flagged as chimeras were removed and were then clustered into OTUs by 97% sequence identity in each library. The taxonomic position of each OTU was automatically assigned on the basis of BLAST analysis, using the SILVA128 database (Quast et al., 2013) as a reference data set of SSU rRNA genes. The sequences closely related to the potential laboratory contaminants described above were omitted from further analyses.

Viral DNA Extraction and Quantification of ssDNA and dsDNA

In this study, sediment viral fraction prepared for the DNA extraction and the subsequent virome analysis was also used as sample for electron microscopic observation. Thus, to avoid damage to viral structure by sonication, we adopted a conventional and versatile virus extraction method (Casas and Rohwer, 2007) for analyzing the viruses in sediment sample. Approximately 50 cm³ frozen sediment from each core section was immediately suspended in 125 ml of modified SM buffer and dispensed into 50-ml centrifuge tubes. The slurry was shaken for 1 min with a ShakeMaster Auto instrument at the maximum speed and centrifuged at 1,500 g for 2 min. The supernatant was then centrifuged at 8,000 g for 10 min. After collection of the supernatant, 125 ml of fresh modified SM buffer was added to the remaining sediments, and the extraction procedure was repeated. The resulting supernatants were pooled (a total of approximately 250 ml) and filtered with a 0.22-μm pore size filter (Millipore). Viral particles in the filtrates were collected by ultracentrifugation at 235,000 g (Type 45 Ti rotor, Optima L-90K ultracentrifuge; Beckman Coulter, Fullerton, CA, United States) for 3 h at 4°C. The pellet was suspended in TM buffer (50 mM Tris-HCl, pH 7.5; 10 mM MgCl₂) and was subjected to cesium chloride density ultracentrifugation at 207,000 g (SW41Ti rotor, Optima L-90K ultracentrifuge; Beckman Coulter) for 48 h at 15°C (initial CsCl density of 0.45 g ml^-1 in TM buffer). Then, a fraction (∼1.2 g ml^-1) potentially containing small microbial cells that passed through 0.2-μm pore size filter (Kuhn et al., 2014) was removed. Finally, the virus particles were collected in a density range from 1.2 to 1.5 g ml^-1 (Thurber et al., 2009; Holmfeldt et al., 2012). The viral fraction was dialyzed in fresh TM buffer for more than 3 h at 4°C, and was ultracentrifuged and collected again under the conditions described above. The pellet was suspended in SM buffer up to 200 μl and treated with DNase I (final concentration of 0.05 U μl^-1; Roche, Mannheim, Germany) for 1 h at 37°C to remove naked extracellular DNA. The addition of EDTA solution (pH 8.0; final concentration of 0.5 M) terminated the reaction. Viral particles were lysed with SDS and proteinase K, and the viral DNA was obtained by the phenol-chloroform extraction and ethanol precipitation (Yoshida et al., 2008). RNA contaminants in the sample were digested with RNase A (final concentration of 0.5 μg μl^-1; Invitrogen, Carlsbad, CA, United States). For the purification of ssDNA, a QIAprep Spin M13 kit (Qiagen, Hilden, Germany) was used. For the purification of dsDNA, ssDNA in the sample was digested with S1 nuclease (final concentration of 0.5 U μl^-1; Takara Bio) at 37°C for 30 min, and the reaction mixture was processed with a QIAamp DNA mini kit (Qiagen). The purified ssDNA and dsDNA were quantified using a Qubit 2.0 Fluorometer (Life Technologies, Eugene, OR, United States) with a Qubit ssDNA assay kit (Life Technologies) and a Qubit dsDNA HS Assay kit (Life Technologies), respectively.

Virome Composition Analysis

The viral ssDNA and dsDNA assemblages were amplified using a REPLI-g Mini Kit (Qiagen, Valencia, CA, United States) for 13 h and were purified by phenol/chloroform/isoamyl alcohol extraction. The amplified DNA products were digested with S1 nuclease at 37°C for 30 min and purified with phenol/chloroform/isoamyl alcohol extraction. Then, the amplified DNA was sheared using a Covaris S220 instrument (Covaris, Woburn, MA, United States) under the manufacturer’s recommended conditions for generating a 350-bp peak. Then, sequencing libraries were constructed using an Ion Plus Fragment Library Kit and an Ion Xpress Barcode Adaptors (Life Technologies, Carlsbad, CA, United States). Emulsion PCR was performed using an Ion PGM Hi-Q OT2 Kit (Life Technologies), and sequencing was performed using an Ion Torrent PGM with an Ion PGM Hi-Q Sequencing Kit and an Ion 318 Chip Kit V2 (Life Technologies). The raw sequence reads were processed using CLC Genomics Workbench ver. 8.0.1 (CLC Bio, Aarhus, Denmark) to trim low-quality reads (with the parameters of a minimum Phred score of 20 and a minimum read length of 50). The processed reads of each virome were subjected to BLASTx analysis (Altschul et al., 1997) against the NCBI GenBank non-redundant (nr) protein database. To identify the viral populations by the virome analysis, the RefseqVirus database including the complete virus genome sequences was also used in the BLASTx analysis. MEGAN (MEtaGenome Analyzer; version 4.70.4) was used to assign taxonomic groups of viruses and cellular organisms (bacteria, archaea, and eukaryota) to query sequences based on BLAST search (E-values>10^-5) (Huson et al., 2011). The rRNA genes in viromes were also evaluated by BLAST analysis against SILVA123 database (Quast et al., 2013). Principal component analysis (PCA) for the comparison of the viromes in this study and other references was performed using the prcomp function in R.

Transmission Electron Microscopy

In this study, sediment viral fraction prepared for the DNA extraction and the subsequent virome analysis was also used as sample for electron microscopic observation. Thus, to avoid damage to viral structure by sonication, we adopted a conventional and versatile virus extraction method (Casas and Rohwer, 2007) for analyzing the viruses in sediment sample. A portion of the viral fraction taken from the 10–15 cmbsf sediment sample was fixed using 1% (v/v) glutaraldehyde. An aliquot of the phage suspension was absorbed onto Formvar/carbon-coated copper grids, followed by staining with 2% (w/v) uranyl acetate and analysis using a TECNAI 20 transmission electron microscope (FEI, Tokyo, Japan) at 120 kV. Images were obtained using a CCD Eagle digital camera (FEI).

Nucleotide Sequence Accession Numbers

The metagenomic data obtained by the multiple virome analyses and SSU rRNA gene tag sequencing of the deep-sea sediment samples have been submitted to the DDBJ Sequence Read Archive (DRA)² under accession number DRA005987.

Results and Discussion

Characteristics of the Sediment Samples

Three sections (0–2, 5–8 and 10–15 cmbsf) of an upper bathyal sediment core (15 cm in total length) from the eutrophic northwest Pacific Ocean were used in this study to obtain a vertical profile of ssDNA and dsDNA viruses in deep-sea sediment (Nomaki et al., 2016a). The geochemical properties (Supplementary Figure S1) and biomasses in the sediment cores have been reported previously (Nomaki et al., 2016a,b). The direct virus and cell counts were reexamined for the virome samples analyzed in this study (Table 1), and they ranged from 5.1 × 10⁵ to 5.3 × 10⁶ particles per cm³ sediment and from 1.6 × 10⁶ to 9.1 × 10⁶ cells per cm³ sediment, respectively. The virus-to-cell ratio (VCR) decreased with increasing sediment depth, ranging from 0.32 to 0.58 (Table 1). The compositions of the small subunit (SSU) rRNA gene tag libraries are summarized in Figure 1. The predominant taxa were shared among the three sections, and their relative abundance changed along the sediment depth. The operational taxonomic units (OTUs) affiliated with Gammaproteobacteria, Deltaproteobacteria, and Bacteroidetes dominated the tag libraries in the 0–2 and 5–8 cmbsf sediments, whereas the Deltaproteobacteria and Chloroflexi OTUs were predominant at deeper depths (10–15 cmbsf). The archaeal OTUs belonging to Thaumarchaeota, Woesearchaeota (DHVEG-6), and Euryarchaeota were observed as minor components throughout all depths (3.8-8.1%).

TABLE 1

TABLE 1. Direct cell and viral counts, the virus-to-cell ratio (VCR), the contents of viral DNAs and the estimated viral genome copy numbers of the ssDNA and dsDNA viral fractions, and the relative abundance of ssDNA viruses in the deep-sea surface sediments.

FIGURE 1

FIGURE 1. The depth profiles of prokaryotic community compositions based on the SSU rRNA gene sequences in deep-sea surface sediments from the off-Tohoku location (water depth = 500 m). The numbers of the sequences analyzed at each sediment depth are shown in parentheses. The “other bacteria” and “other archaea” categories represent the bacterial and archaeal taxa, respectively, and share less than 2% of the total sequences. The “others” category represents eukaryota and unclassified sequences.

Construction and Evaluation of ssDNA and dsDNA Virome Libraries

Both ssDNA and dsDNA were extracted from viral particles in the sediment samples. The content of the viral ssDNA in the sediment samples ranged from 1.50 to 2.87 ng per cm³ sediment and was approximately one order of magnitude higher than those of the viral dsDNA (0.17–0.29 ng per cm³ sediment) (Table 1).

To examine the viral ssDNA and dsDNA assemblages, a total of six metagenomic libraries were constructed (Table 2). Totals of 1,091,418; 725,380; and 213,547 sequence reads were obtained from the ssDNA virome libraries from the 0–2, 5–8, and 10–15 cmbsf sediments, respectively. For the virome libraries prepared from the dsDNA viral fraction, totals of 751,300; 647,462; and 271,392 sequence reads were taken from the 0–2, 5–8, and 10–15 cmbsf sediment samples, respectively.

TABLE 2

TABLE 2. Overview of quality-trimmed reads of the deep-sea sedimentary virome libraries analyzed in this study.

To verify the quality of the viromes, their sequence reads were subjected to BLASTx analysis against the GenBank nr protein database as a reference data set of known genes (Figure 2). Only 4–17% of the sequence reads in each library showed significant similarity (E-value < 10^-5 in BLASTx) to the sequences in the database and were subsequently classified into taxonomic groups by MEGAN analysis.

FIGURE 2

FIGURE 2. Summary of the MEGAN taxonomic assignments of the off-Tohoku sedimentary virome reads on the basis of BLASTx similarity search (E-value < 10^-5) against the nr database. The numbers of the sequenced reads in each virome library are shown in parentheses.

In the three ssDNA virome libraries, virus-associated reads were most abundant among the taxonomically assigned reads (9–13%), whereas more than 80% of each library exhibited no significant similarity to known sequences in the nr database (Figure 2). In contrast, in the dsDNA virome libraries, only 0.3–3% of sequence reads were classified as viral. The greater reads (3–5%) in these libraries showed significant homology to bacterial sequences, which was similar to the results of previous studies (e.g., Angly et al., 2006; Dinsdale et al., 2008). Most (75–94%) of the total annotated reads in each library were not assigned to any known KEGG metabolic functions (Kanehisa et al., 2012). Thus, the sequences assigned as non-viral are thought to be a result of the misclassification of viral sequences as host (bacteria, archaea, and eukaryota) genomic components, owing to both the lack of a comprehensive viral database and the horizontal gene transfer between viral and host genomes (Angly et al., 2006; Roux et al., 2012).

Among the reads classified as bacteria and archaea in each dsDNA virome library, the abundance of rRNA gene-associated reads was only 0–0.005% [for reference, length ratios of rRNA genes in genome sizes (1.5–5.5 Mb; Angly et al., 2009) inferred across aquatic bacterial and archaeal metagenomes: 0.03–0.1%], and it accounted for 0–0.0002% of the total virome reads (Table 2). According to the virome quality criteria proposed on the basis of the distribution of the rRNA gene-associated read ratios of the examined 67 public viromes (<0.02% of the total reads; Roux et al., 2013), our virome libraries produce relatively high quality with low contamination by cellular sequences.

Composition of ssDNA and dsDNA Viromes

To examine the composition of ssDNA and dsDNA viromes, all reads of the sedimentary virome libraries (Table 2) were reevaluated by BLASTx similarity search (E-value < 10^-5) against the RefseqVirus database including the complete virus genome sequences (Figure 3). Then, the reads that had significant similarity to the known viral sequences were defined as “viral reads” and were used for further analyses. The relative abundance of the viral reads in the ssDNA and dsDNA virome libraries ranged from 7.2 to 9.5% and from 0.4 to 4.1%, respectively.

FIGURE 3

FIGURE 3. Summary of the MEGAN taxonomic assignments of viral reads in the off-Tohoku sedimentary virome libraries. The viral reads were identified by BLASTx similarity search (E-value < 10^-5) against the viral sequence database constructed in this study. The numbers of the sequenced reads in each virome library are shown in parentheses. The “other ssDNA viruses” and “other dsDNA viruses” categories represent unclassified viruses and viral families constituting less than 2% of the total viral reads. The “others” comprises unclassified viruses and other components, single-stranded RNA viruses, retro-transcribing viruses, and virophages.

In all the ssDNA virome libraries, most of the viral reads were related to genes from ssDNA viral families, which include Circoviridae and Microviridae, and they consisted of 85–86% of the total viral reads in each library (Figure 3). The family Circoviridae viruses and the marine relatives are known to infect eukaryotic hosts (animals and marine diatoms, respectively), whereas Microviridae is a bacteriophage group (e.g., Rosario et al., 2009; Tomaru et al., 2011). The sequences associated with these ssDNA viral families have also been found in other deep-sea surface sediments (Yoshida et al., 2013; Tangherlini et al., 2016), as well as in oceanic waters (e.g., Angly et al., 2006). The predominant Inoviridae-like sequences which predominated in the viral sequences from deep subseafloor metatranscriptomes (Engelhardt et al., 2015) were detected as minor components of the surface sediment ssDNA viromes. The reads originating from dsDNA viruses were rarely detected in the three ssDNA libraries (1–2%).

Among the dsDNA libraries, the viral read assemblages in the 0–2 and 5–8 cmbsf sections were dominated (88–94%) by sequences related to dsDNA viral families, which includes Siphoviridae, Podoviridae, Myoviridae, and Phycodnaviridae (Figure 3). The first three families are typical “tailed bacteriophages” belonging to the order Caudovirales and commonly occur in various marine and freshwater environments (Suttle, 2005; Roux et al., 2012), whereas Phycodnaviridae are aquatic viruses that infect eukaryotic algae (Larsen et al., 2008). The Phycodnaviridae viruses were detected as relatively minor populations in the dsDNA libraries (2.3–2.7% of the total viral reads). Because those families include viruses larger than 0.2 μm, some of them are likely to be selectively removed by filtration. The reads associated with ssDNA viruses ranged from only 3 to 8% of the total viral reads in the dsDNA libraries. In contrast to these dsDNA libraries, the relative abundance of the sequences associated with ssDNA viruses in the 10–15 cmbsf dsDNA library was unexpectedly greater than that associated with dsDNA viruses, sharing 69% of the total viral reads. This result was probably because of potential contamination of ssDNA due to the incomplete S1 nuclease digestion of ssDNA during the sediment viral dsDNA sample preparation. The overall composition of the virome libraries suggested that the viral ssDNA and dsDNA were generally well prepared from the sediment viral communities, although contamination of DNA types other than the target DNA type was not completely excluded from each DNA fraction sample.

To examine the relationships among the different depth viromes in the present study and compare them with deep-sea sediment ones at other sites that we previously examined (Yoshida et al., 2013), the three data sets of the ssDNA and dsDNA viral sequence components in the relevant off-Tohoku sedimentary virome libraries (Figure 3) and the references were subjected to PCA (Figure 4). PCA for ssDNA viral components (Figure 4A) indicated that the different depth viral communities at the same site all exhibited similar profiles compared with those at other sites, while PCA for dsDNA viral components (Figure 4B) indicated that the different depth viral communities at the same site all revealed distinct relationships than with those at other sites.

FIGURE 4

FIGURE 4. Principal component analysis (PCA) analysis of (A) ssDNA and (B) dsDNA viral sequence components in the relevant deep-sea sedimentary virome libraries from the off-Tohoku and other reference sites (Ogasawara Trench, Mariana Trench, and off Shimokita Peninsula) based on the taxonomic composition of viral communities. The respective off-Tohoku viral communities with the three different depths are shown in blue and the references are shown in red.

Abundance of ssDNA and dsDNA Viruses

The abundance of ssDNA and dsDNA viruses in the samples was estimated as the content of the viral ssDNA and dsDNA extracts (Table 1) divided by each DNA content per virus particle, which was calculated on the basis of the presumed genome sizes of currently known ssDNA and dsDNA viral families. The genome sizes of the ssDNA viruses, ranging from 1.8 to 24.9 kb, were determined according to the latest virus taxonomy on the ViralZone website³ (Hulo et al., 2011), and the DNA content per virus particle ranges from 9.90 × 10^-19g (for 1.8 kb genome) to 1.37 × 10^-17g (for 24.9 kb genome). The abundance of ssDNA viruses in the sediments of our study site ranged from 1.1 × 10⁸ to 2.9 × 10⁹ genome copies per cm³ sediment (Table 1). For the calculation of dsDNA viral abundance, we assumed an average DNA content per virus particle of 5.5 × 10^-17 g (equivalent to 50 kb dsDNA) according to a previous estimate (Steward et al., 2013). The abundance of dsDNA viruses was estimated to be much lower than those of ssDNA viruses, ranging from 3.0 × 10⁶ to 5.3 × 10⁶ genome copies per cm³ sediment. Thus, we can provide estimation that the benthic ssDNA viral populations represent 96.3–99.8% of the benthic total DNA viral assemblages (Table 1). Further, morphological analysis of a sediment viral fraction (10-15 cmbsf sample) by transmission electron microscopy revealed the prevalence of non-tailed virus-like particles [diameter 31-62 nm (mean 47 nm), n = 20 particles] (e.g., Supplementary Figure S2) with sizes that were consistent with isolated marine ssDNA viral particles (diameter 30–73 nm; Brum et al., 2013), whereas there was no observation of tailed viruses, generally associated with dsDNA viruses. These data supported the results of the ssDNA and dsDNA viral quantification.

Ecological Implications of ssDNA Viral Prevalence

Direct virus counting with SYBR staining is generally applicable to dsDNA viruses, whereas ssDNA and RNA viruses are often insufficiently stained, owing to their small genomes (Tomaru and Nagasaki, 2007; Tucker et al., 2011; Holmfeldt et al., 2012). In this study, the quantification of viral dsDNA and the subsequent calculation of viral abundance suggested that the estimated dsDNA viral abundances were comparable to the direct virus counts (Table 1). In contrast, the estimated ssDNA viral abundances were at least 200-fold higher than those obtained by direct virus count. The results suggested that the abundance of deep-sea sedimentary ssDNA viruses would have been extraordinarily underestimated by fluorescence-based direct counts. The benthic viral infection is responsible for the abatement of up to ∼80% of heterotrophic carbon production by bacteria and archaea, thus causing the release of ∼0.37 to 0.63 gigatons of carbon per year on a global scale (Danovaro et al., 2008b). Our quantitative approach to both ssDNA and dsDNA viral communities provides new insights into the re-evaluation of viral abundance and ecological functions in deep-sea sediments and illuminates the potentially significant contribution of ssDNA viral populations to the bacterial and archaeal mortality and biogeochemical processes of deep-sea benthic microbial communities by the viral shunt. Future work needs to determine whether other diverse marine benthic virosphere, including abyssal and trench sediments and deep subseafloor sediments, reveal the relatively high abundance of ssDNA viruses compared to that of dsDNA viruses.

Author Contributions

MY conceived and designed the experiments, performed the experiments, analyzed the data, and wrote the manuscript. HN performed the onboard sample processing and critically edited the manuscript. TM, S-IU, and YY-T performed the laboratory experiments and critically edited the manuscript. MH performed the molecular analyses. SN and YT performed the bioinformatic analyses and critically edited the manuscript. TN and KT conceived and designed the experiments and wrote the manuscript.

Funding

This work was supported in part by a Grant-in-Aid for Exploratory Research (26660154) from the Japan Society for the Promotion of Science (JSPS), grants for Scientific Research on Innovative Areas (16H06429, 16K21723, 17H05830, and 16H06437) from the Ministry of Education, Culture, Science, Sports and Technology (MEXT) of Japan, and the research project Tohoku Ecosystem-Associated Marine Sciences from MEXT.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The reviewer KW and handling Editor declared their shared affiliation.

Acknowledgments

We would like to thank the captain, crew, and onboard scientists and technicians of the R/V Mirai (JAMSTEC) during the MR12-E02 cruise. We are grateful to Miwako Tsuda and Juliarni and for technical assistance with SSU rRNA gene tag sequencing. We are also grateful to Katsuyuki Uematsu for technical help with transmission electron microscopy.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.00075/full#supplementary-material

Footnotes

1. ^http://hannonlab.cshl.edu/fastx_toolkit/index.html
2. ^http://trace.ddbj.nig.ac.jp/dra/index_e.shtml
3. ^www.expasy.org/viralzone

References

Altschul, S. F., Madden, T. L., Schaffer, A. A., Zhang, J. H., Zhang, Z., Miller, W., et al. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389–3402. doi: 10.1093/nar/25.17.3389

PubMed Abstract | CrossRef Full Text | Google Scholar

Angly, F. E., Felts, B., Breitbart, M., Salamon, P., Edwards, R. A., Carlson, C., et al. (2006). The marine viromes of four oceanic regions. PLOS Biol. 4:e368. doi: 10.1371/journal.pbio.0040368

PubMed Abstract | CrossRef Full Text | Google Scholar

Angly, F. E., Willner, D., Prieto-Davó, A., Edwards, R. A., Schmieder, R., Vega-Thurber, R., et al. (2009). The GAAS metagenomic tool and its estimations of viral and microbial average genome size in four major biomes. PLOS Comput. Biol. 5:e1000593. doi: 10.1371/journal.pcbi.1000593

PubMed Abstract | CrossRef Full Text | Google Scholar

Brum, J. R., Schenck, R. O., and Sullivan, M. B. (2013). Global morphological analysis of marine viruses shows minimal regional variation and dominance of non-tailed viruses. ISME J. 7, 1738–1751. doi: 10.1038/ismej.2013.67

PubMed Abstract | CrossRef Full Text | Google Scholar

Brum, J. R., and Sullivan, M. B. (2015). Rising to the challenge: accelerated pace of discovery transforms marine virology. Nat. Rev. Microbiol. 13, 147–159. doi: 10.1038/nrmicro3404

PubMed Abstract | CrossRef Full Text | Google Scholar

Caporaso, J. S., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows integration and analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303

PubMed Abstract | CrossRef Full Text | Google Scholar

Casas, V., and Rohwer, F. (2007). Phage metagenomics. Methods Enzymol. 421, 259–268. doi: 10.1016/S0076-6879(06)21020-6

CrossRef Full Text | Google Scholar

Chen, F., Lu, J. R., Binder, B. J., Liu, Y. C., and Hodson, R. E. (2001). Application of digital image analysis and flow cytometry to enumerate marine viruses stained with SYBR Gold. Appl. Environ. Microbiol. 67, 539–545. doi: 10.1128/AEM.67.2.539-545.2001

PubMed Abstract | CrossRef Full Text | Google Scholar

Danovaro, R., Corinaldesi, C., Filippini, M., Fischer, U. R., Gessner, M. O., Jacquet, S., et al. (2008a). Viriobenthos in freshwater and marine sediments: a review. Freshw. Biol. 53, 1186–1213. doi: 10.1111/j.1365-2427.2008.01961.x

CrossRef Full Text | Google Scholar

Danovaro, R., Dell’Anno, A., Corinaldesi, C., Magagnini, M., Noble, R., Tamburini, C., et al. (2008b). Major viral impact on the functioning of benthic deep-sea ecosystems. Nature 454, 1084–1087. doi: 10.1038/nature07268

PubMed Abstract | CrossRef Full Text | Google Scholar

Danovaro, R., Dell’Anno, A., Corinaldesi, C., Rastelli, E., Cavicchioli, R., Krupovic, M., et al. (2016). Virus-mediated archaeal hecatomb in the deep seafloor. Sci. Adv. 2:e1600492. doi: 10.1126/sciadv.1600492

PubMed Abstract | CrossRef Full Text | Google Scholar

Dean, F., Nelson, J., Giesler, T., and Lasken, R. (2001). Rapid amplification of plasmid and phage DNA using phi29 DNA polymerase and multiply-primed rolling circle amplification. Genome Res. 11, 1095–1099. doi: 10.1101/gr.180501

PubMed Abstract | CrossRef Full Text | Google Scholar

Dinsdale, E. A., Edwards, R. A., Hall, D., Angly, F., Breitbart, M., Brulc, J. M., et al. (2008). Functional metagenomic profiling of nine biomes. Nature 452, 344–347. doi: 10.1038/nature06810

PubMed Abstract | CrossRef Full Text | Google Scholar

Edwards, R. A., and Rohwer, F. (2005). Viral metagenomics. Nat. Rev. Microbiol. 3, 504–510. doi: 10.1038/nrmicro1163

PubMed Abstract | CrossRef Full Text | Google Scholar

Engelhardt, T., Orsi, W. D., and Jørgensen, B. B. (2015). Viral activities and life-cycles in deep subseafloor sediment. Environ. Microbiol. Rep. 7, 868–873. doi: 10.1111/1758-2229.12316

PubMed Abstract | CrossRef Full Text | Google Scholar

Fuhrman, J. A. (1999). Marine viruses and their biogeochemical and ecological effects. Nature 399, 541–548. doi: 10.1038/21119

PubMed Abstract | CrossRef Full Text | Google Scholar

Hirai, M., Nishi, S., Tsuda, M., Sunamura, M., Takaki, Y., and Nunoura, T. (2017). Library construction from subnanogram DNA for pelagic sea water and deep-sea sediments. Microbes Environ. 32, 336–343. doi: 10.1264/jsme2.ME17132

PubMed Abstract | CrossRef Full Text | Google Scholar

Holmfeldt, K., Odić, D., Sullivan, M. B., Middelboe, M., and Riemann, L. (2012). Cultivated single-stranded DNA phages that infect marine Bacteroidetes prove difficult to detect with DNA-binding stains. Appl. Environ. Microbiol. 78, 892–894. doi: 10.1128/AEM.06580-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Hulo, C., de Castro, E., Masson, P., Bougueleret, L., Bairoch, A., Xenarios, I., et al. (2011). ViralZone: a knowledge resource to understand virus diversity. Nucleic Acids Res. 39, D576–D582. doi: 10.1093/nar/gkq901

PubMed Abstract | CrossRef Full Text | Google Scholar

Huson, D. H., Mitra, S., Ruscheweyh, H. J., Weber, N., and Schuster, S. C. (2011). Integrative analysis of environmental sequences using MEGAN4. Genome Res. 21, 1552–1560. doi: 10.1101/gr.120618.111

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanehisa, M., Goto, S., Sato, Y., Furumichi, M., and Tanabe, M. (2012). KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 40, D109–D114. doi: 10.1093/nar/gkr988

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuhn, E., Ichimura, A. S., Peng, V., Fritsen, C. H., Trubl, G., Doran, P. T., et al. (2014). Brine assemblages of ultrasmall microbial cells within the ice cover of Lake Vida, Antarctica. Appl. Environ. Microbiol. 80, 3687–3698. doi: 10.1128/AEM.00276-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Labonté, J. M., and Suttle, C. A. (2013). Previously unknown and highly divergent ssDNA viruses populate the oceans. ISME J. 7, 2169–2177. doi: 10.1038/ismej.2013.110

PubMed Abstract | CrossRef Full Text | Google Scholar

Larsen, J. B., Larsen, A., Bratbak, G., and Sandaa, R. A. (2008). Phylogenetic analysis of members of the Phycodnaviridae virus family, using amplified fragments of the major capsid protein gene. Appl. Environ. Microbiol. 74, 3048–3057. doi: 10.1128/AEM.02548-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Lizardi, P. M., Huang, X., Zhu, Z., Bray-Ward, P., Thomas, D. C., and Ward, D. C. (1998). Mutation detection and single-molecule counting using isothermal rolling-circle amplification. Nat. Genet. 19, 225–232. doi: 10.1038/898

PubMed Abstract | CrossRef Full Text | Google Scholar

Mann, N. H., Cook, A., Millard, A., Bailey, S., and Clokie, M. (2003). Marine ecosystems: bacterial photosynthesis genes in a virus. Nature 424, 741–742. doi: 10.1038/424741a

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin, M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 17, 10–12. doi: 10.14806/ej.17.1.200

CrossRef Full Text | Google Scholar

Middelboe, M., Glud, R. N., and Filippini, M. (2011). Viral abundance and activity in the deep sub-seafloor biosphere. Aquat. Microb. Ecol. 63, 1–8. doi: 10.3354/ame01485

CrossRef Full Text | Google Scholar

Miranda, J. A., Culley, A. I., Schvarcz, C. R., and Steward, G. F. (2016). RNA viruses as major contributors to Antarctic virioplankton. Environ. Microbiol. 18, 3714–3727. doi: 10.1111/1462-2920.13291

PubMed Abstract | CrossRef Full Text | Google Scholar

Noble, R. T., and Fuhrman, J. A. (1998). Use of SYBR Green I for rapid epifluorescence counts of marine viruses and bacteria. Aquat. Microb. Ecol. 14, 113–118. doi: 10.3354/ame014113

CrossRef Full Text | Google Scholar

Nomaki, H., Arai, K., Suga, H., Toyofuku, T., Wakita, M., Nunoura, T., et al. (2016a). Sedimentary organic matter contents and porewater chemistry at upper bathyal depths influenced by the 2011 off the Pacific coast of Tohoku Earthquake and tsunami. J. Oceanogr. 72, 99–111. doi: 10.1007/s10872-015-0315-3

CrossRef Full Text | Google Scholar

Nomaki, H., Mochizuki, T., Kitahashi, T., Nunoura, T., Arai, K., Toyofuku, T., et al. (2016b). Effects of mass sedimentation events after the 2011 off the Pacific coast of Tohoku Earthquake on benthic prokaryotes and meiofauna inhabiting the upper bathyal sediments. J. Oceanogr. 72, 113–128. doi: 10.1007/s10872-015-0293-5

CrossRef Full Text | Google Scholar

Nunoura, T., Takaki, Y., Kazama, H., Hirai, M., Ashi, J., Imachi, H., et al. (2012). Microbial diversity in deep-sea methane seep sediments presented by SSU rRNA gene tag sequencing. Microbes Environ. 27, 382–390. doi: 10.1264/jsme2.ME12032

PubMed Abstract | CrossRef Full Text | Google Scholar

Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219

PubMed Abstract | CrossRef Full Text | Google Scholar

Rosario, K., Duffy, S., and Breitbart, M. (2009). Diverse circovirus-like genome architectures revealed by environmental metagenomics. J. Gen. Virol. 90, 2418–2424. doi: 10.1099/vir.0.012955-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Roux, S., Enault, F., Robin, A., Ravet, V., Personnic, S., Theil, S., et al. (2012). Assessing the diversity and specificity of two freshwater viral communities through metagenomics. PLOS ONE 7:e33641. doi: 10.1371/journal.pone.0033641

PubMed Abstract | CrossRef Full Text | Google Scholar

Roux, S., Krupovic, M., Debroas, D., Forterre, P., and Enault, F. (2013). Assessment of viral community functional potential from viral metagenomes may be hampered by contamination with cellular sequences. Open. Biol. 3:130160. doi: 10.1098/rsob.130160

PubMed Abstract | CrossRef Full Text | Google Scholar

Roux, S., Solonenko, N. E., Dang, V. T., Poulos, B. T., Schwenck, S. M., Goldsmith, D. B., et al. (2016). Towards quantitative viromics for both double-stranded and single-stranded DNA viruses. PeerJ 4:e2777. doi: 10.7717/peerj.2777

PubMed Abstract | CrossRef Full Text | Google Scholar

Steward, G. F., Culley, A. I., Mueller, J. A., Wood-Charlson, E. M., Belcaid, M., and Poisson, G. (2013). Are we missing half of the viruses in the ocean? ISME J. 7, 672–679. doi: 10.1038/ismej.2012.121

PubMed Abstract | CrossRef Full Text | Google Scholar

Suttle, C. A. (2005). Viruses in the sea. Nature 437, 356–361. doi: 10.1038/nature04160

PubMed Abstract | CrossRef Full Text | Google Scholar

Suttle, C. A. (2007). Marine viruses - major players in the global ecosystem. Nat. Rev. Microbiol. 5, 801–812. doi: 10.1038/nrmicro1750

PubMed Abstract | CrossRef Full Text | Google Scholar

Tangherlini, M., Dell’Anno, A., Allen, L. Z., Riccioni, G., and Corinaldesi, C. (2016). Assessing viral taxonomic composition in benthic marine ecosystems: reliability and efficiency of different bioinformatic tools for viral metagenomic analyses. Sci. Rep. 6:28428. doi: 10.1038/srep28428

PubMed Abstract | CrossRef Full Text | Google Scholar

Thurber, R. V., Haynes, M., Breitbart, M., Wegley, L., and Rohwer, F. (2009). Laboratory procedures to generate viral metagenomes. Nat. Protoc. 4, 470–483. doi: 10.1038/nprot.2009.10

PubMed Abstract | CrossRef Full Text | Google Scholar

Tomaru, Y., and Nagasaki, K. (2007). Flow cytometric detection and enumeration of DNA and RNA viruses infecting marine eukaryotic microalgae. J. Oceanogr. 63, 215–221. doi: 10.1007/s10872-007-0023-8

CrossRef Full Text | Google Scholar

Tomaru, Y., Takao, Y., Suzuki, H., Nagumo, T., Koike, K., and Nagasaki, K. (2011). Isolation and characterization of a single-stranded DNA virus infecting Chaetoceros lorenzianus Grunow. Appl. Environ. Microbiol. 77, 5285–5293. doi: 10.1128/AEM.00202-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Tucker, K. P., Parsons, R., Symonds, E. M., and Breitbart, M. (2011). Diversity and distribution of single-stranded DNA phages in the North Atlantic Ocean. ISME J. 5, 822–830. doi: 10.1038/ismej.2010.188

PubMed Abstract | CrossRef Full Text | Google Scholar

Weinbauer, M. G. (2004). Ecology of prokaryotic viruses. FEMS Microbiol. Rev. 28, 127–181. doi: 10.1016/j.femsre.2003.08.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Wommack, K. E., and Colwell, R. R. (2000). Virioplankton: viruses in aquatic ecosystems. Microbiol. Mol. Biol. Rev. 64, 69–114. doi: 10.1128/MMBR.64.1.69-114.2000

CrossRef Full Text | Google Scholar

Yanagawa, K., Morono, Y., Yoshida-Takashima, Y., Eitoku, M., Sunamura, M., Inagaki, F., et al. (2014). Variability of subseafloor viral abundance at the geographically and geologically distinct continental margins. FEMS Microbiol. Ecol. 88, 60–68. doi: 10.1111/1574-6941.12269

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoshida, M., Takaki, Y., Eitoku, M., Nunoura, T., and Takai, K. (2013). Metagenomic analysis of viral communities in (hado)pelagic sediments. PLOS ONE 8:e57271. doi: 10.1371/journal.pone.0057271

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoshida, M., Yoshida, T., Kashima, A., Takashima, Y., Hosoda, N., Nagasaki, K., et al. (2008). Ecological dynamics of the toxic bloom-forming cyanobacterium Microcystis aeruginosa and its cyanophages in freshwater. Appl. Environ. Microbiol. 74, 3269–3273. doi: 10.1128/AEM.02240-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Kober, K., Flouri, T., and Stamatakis, A. (2014). PEAR: a fast and accurate Illumina paired-end reAd mergeR. Bioinformatics 30, 614–620. doi: 10.1093/bioinformatics/btt593

PubMed Abstract | CrossRef Full Text | Google Scholar