microbial ecologists have devoted considerable effort to understanding the nature of the viruses in seawater , because viruses have key roles in the evolution , ecology and mortality of marine plankton ( rohwer and vega thurber , 2007 ; suttle , 2007 ) . for at least the past two decades , researchers have assumed that the pool of viruses in the ocean is dominated by bacteriophages with dna genomes ( steward et al . , 1992 ; breitbart et al . , 2002 ; weinbauer , 2004 ; comeau et al . , 2010 ; sullivan et al . , 2010 ) .
perhaps as a consequence , studies of the molecular diversity of marine viruses have most commonly ( exclusively , before 2003 ) focused on dna viruses ( edwards and rohwer , 2005 ; kristensen et al . ,
however , evidence that rna viruses are important contributors to marine plankton ecology has been steadily accumulating ( lang et al . , 2009 ) .
the isolation of a positive - sense , single - stranded rna ( ssrna ) virus that infects the raphidophyte heterosigma akashiwo ( harnav ) was the first recorded instance of an rna virus infecting a marine protist ( tai et al . , 2003 ) .
this was followed by reports of similar picorna - like ssrna viruses infecting diatoms ( nagasaki et al . , 2004 ; shirai et al . , 2008 ; tomaru et al . ,
molecular surveys using degenerate primers to target the rna - dependent rna polymerase gene of picorna - like viruses have shown that , in addition to the handful of isolates , a very diverse pool of uncultivated picornavirads exist in seawater ( culley et al . , 2003 ;
a double - stranded rna virus infecting the abundant marine prymnesiophyte micromonas pusilla represents a new genus in the family reoviridae ( brussaard et al . , 2004 ) . a positive - sense , ssrna virus that infects the dinoflagellate heterocapsa circularisquama ( tomaru et al . , 2004
) is only distantly related to existing viral families ( nagasaki et al . , 2005 ) and may represent a new family .
metagenomic surveys of the rna viral fraction of seawater are consistent with the view provided by the limited number of isolates available so far , and suggest that the positive - sense , single - stranded picornavirads dominate the marine rna virus pool , but other diverse rna viruses , including some containing double - stranded rna , are present as well ( culley et al . , 2006 ) . despite the emerging evidence that marine rna viruses are diverse and infect ecologically important members of the marine planktonic food web ,
there have been no reports that satisfactorily address the question of whether these viruses constitute a substantial fraction of the total virioplankton .
this question is of considerable ecological interest , because , unlike the pool of dna viruses in seawater , which is composed predominantly of bacteriophage - like sequences ( edwards and rohwer , 2005 ) , marine rna viruses are almost exclusively comprised of those that infect eukaryotes ( lang et al . , 2009 ) . directly quantifying the abundance of rna viruses in a mixed viral assemblage has proven difficult , because of technical limitations .
differences in counts using a dna - specific stain vs a nonspecific nucleic acid stain have been reported ( weinbauer and suttle , 1997 ; guixa - boixereu et al . , 1999 ;
bettarel et al . , 2000 ) , but the results are of uncertain significance , because the stains differ in their sensitivity ( weinbauer , 2004 ) .
in particular , the small genomes of single - stranded dna and rna viruses make the individual virions difficult to detect even with appropriate stains ( brussaard et al .
2012 ) . given the difficulties of obtaining reliable direct counts of rna viruses , we took a different approach . in the work reported here
, we measured the relative masses of rna and dna in natural assemblages of viruses purified from tropical coastal seawater and coupled this with estimates of the mass of nucleic acid per rna or dna virion to obtain the first estimates of the relative abundance of rna viruses in seawater .
the samples for this study were collected from a pier in the southern portion of kne'ohe bay on the windward side of o'ahu , hawaii ( 21 25 46.80 n , 157 47 31.51 w ) .
this tropical embayment is characterized by year - round warm temperatures ( 22  c to 28  c ) and salinities ranging from 32 to 35 , except during periods of heavy rain , when freshwater plumes from stream runoff can transiently suppress salinity to < 30 ( drupp et al . , 2011 ) .
concentrations of chlorophyll a in surface waters are typically low with average values recorded over a 2.5-year period ranging from 0.5 to 1.1  g  l in summer and from 1.3 to 3.4  g 
l in the winter depending on location ( drupp et al . , 2011 ) . following heavy rainfall ,
blooms have been reported with transient increases in chlorophyll a to ca 6  g 
l ( hoover et al . , 2006 ; de carlo et al . , 2007 )
seawater was collected in acid - washed , polycarbonate carboys on 1 august 2009 ( 35 liters ) and on 3 june 2010 ( 80 liters ) during non - storm conditions .
seawater was transported immediately to the laboratory ( < 1  h ) and filtered through 0.22  m pore - size filters ( sterivex , millipore , billerica , ma , usa ) . the sample from june
viruses in the filtrates were concentrated by iron flocculation ( john et al . , 2011 ) then concentrated further in a centrifugal 30  kda ultrafiltration device ( amicon 15 , millipore ) .
viruses in each concentrate were purified by a two - step process in cscl buoyant density gradients ( lawrence and steward , 2010 )
. fractions ( ca 0.5  ml each ) were collected from the final gradient for analysis with an auto densi - flow ( labconco , kansas city , mo , usa ) .
the density of each fraction was measured using a micropipet and an analytical balance ( lawrence and steward , 2010 ) .
cscl was exchanged with buffer ( 10  m tris , 1  m edta , ph 8) by centrifugal ultrafiltration ( steward , 2001 ) and total nucleic acids in a subsample of each fraction was extracted using spin columns ( qiaamp minelute viral spin kit , qiagen , valencia , ca , usa ) according to the manufacturer 's instructions .
samples for rna analysis were treated twice with dnase ( turbo dnase , life technologies , carlsbad , ca , usa ) to avoid nonspecific signal .
rna and dna contents of each fraction were measured separately in parallel subsamples by fluorometry ( quant - it dna and rna assays , life technologies ) in a cuvette fluorometer ( td-700 , turner designs , sunnyvale , ca , usa ) .
the rna content of the individual fractions from the june sample was low , so we also measured the rna content of the pooled putative viral fractions to obtain a more accurate estimate .
the nucleic acid masses ( 95% confidence interval ( ci ) ) were calculated by inverse prediction from the standard curves ( zar , 1996 ) .
the fluorometric assay was necessary to achieve the required sensitivity . to ensure that we did not overestimate viral rna with the fluorometric assay
, we compared the measurements of the mass of rna obtained from a purified rna virus by fluorometry and by spectrophometry ( model du 800 , beckman coulter , inc . , brea , ca , usa ) .
the virus was a positive - sense , single - stranded , picornavirad that infects a marine diatom and has a genome of 8800 nucleotides ( nt ) ( schvarcz et al .
we also tested for cross - reaction of dna in the rna fluorescence assay with and without the dnase digestion procedure described above ( figure 1 ) .
the fluorometer was calibrated with the quant - it kit standards ( escherichia coli rna ) .
the rna standards were then read as samples along with undigested and digested dna standards and the purified viral rna .
apparent rna concentration for each sample as determined by fluorometry was plotted as a function of the given rna or dna concentrations ( quant - it kit standards ) or , in the case of the viral rna , the concentration determined by absorption at 260  nm ( sambrook and russell , 2001 ) .
model i regression lines were calculated and a two - tailed t - test used to determine if the slopes were significantly different from zero ( zar , 1996 ) . to determine the composition of the rna in the selected rna peaks , we created and analyzed metagenomes from this material for the august sample and from one of the duplicate subsamples ( replicate 1 ) from june .
purified , dnase - treated rna from the indicated fractions ( figure 2 ) of each sample was pooled and amplified using random priming - mediated sequence - independent single - primer amplification as described previously ( djikeng et al . , 2008 ; culley et al . , 2010 ) .
tests of this method on rna viral genomes indicated that it results in coverage and redundancy similar to the ideal values that are predicted by the lander  waterman model ( djikeng et al . , 2008 ) , which suggests amplification biases are limited .
sequence libraries from the resulting amplified complementary dna samples were produced by pyrosequencing ( gs flx titanium , 454 life sciences , branford , ct , usa ) .
sequences that were > 2 s.d . from the mean length or had an average phred score of < 15 were discarded using the mg - rast qc pipeline ( http://metagenomics.anl.gov/ ) and sequence - independent single - primer amplification primers were trimmed from the remaining reads .
artificial ( or technical ) replicates ( gomez - alvarez et al . , 2009 )
were removed using the online tool mg - rast ( meyer et al . , 2008 ) .
processed sequences were assembled with clc genomics workbench version 5.0 ( clcbio , cambridge , ma , usa ) using global alignment with automatic word and bubble sizes , a minimum contig length of 200 , mismatch , insertion and deletion costs set to 3 , length fraction set to 0.5 , and the similarity threshold set to 0.8 .
the community composition of each sample was analyzed with megan ( huson et al . ,
2011 ) using the output from blastx comparisons ( altschul et al . , 1990 ) of the assembled metagenomes ( contigs plus singletons ) with the non - redundant ncbi sequence database .
the taxonomic assignment for a given contig was applied to all reads comprising that contig .
megan assignments were manually checked and some reassignments made in the case of annotation errors . as an independent check on possible contamination from cellular rna , sequences were compared with the silva database ( pruesse et al . , 2007 ) using blastn to search for ribosomal rna contamination , which is expected to be the major source of cellular rna contamination ( karpinets et al . , 2006 ) .
the mass of rna in the pooled viral fraction for each library was adjusted downward by the percentage of reads identified as being cell derived . to estimate the number of rna virus genomes
, we first summed the rna mass within a narrow buoyant density range ( between the solid vertical lines , figure 2 ) .
although rna viruses can have densities outside of that range , we used this conservative window because we analyzed the rna in only those fractions by metagenomic analysis . using the relative representation of different taxa in the metagenomic library and the average genome sizes for those taxa derived from data in the ninth report of the international committee on the taxonomy of viruses ( king et al . , 2012 ) , we calculated a weighted average rna mass per virion . for our prima facie estimates of the rna virus contribution , we assumed that the sequences having no significant hit to any sequences in genbank ( either directly or by association with other sequences in a contig ) had the same taxonomic distribution as the collection of sequences that did have significant blast hits .
the portion of the total rna determined to be viral was then divided by the weighted average mass of rna per virion to obtain the number of rna - containing viruses .
errors were calculated ( and propagated ) as the 95% ci of the inverse predictions of rna mass from the standard curves . to estimate abundance of dna viruses
, we summed the dna content in the fractions encompassing the main dna peak in the viral buoyant density range ( between the dashed vertical lines , figure 2 ) and converted the total dna mass to numbers of dna viruses assuming an average dna content per virion of 5.5  10  g ( equivalent to 50  kb double - stranded dna ) , an average that was found to be similar in a wide range of environments ( steward et al . , 2000 ) .
this value is conservative compared with some other estimates ( brum , 2005 ; angly et al . , 2009 ) and
we also calculated the contribution of rna viruses using more extreme assumptions to get a sense of how much higher and lower the percentage might be .
for the high estimate , we assumed that the dna per virion was twofold greater ( 100  kb , double - stranded dna ) as calculated for marine viral metagenomic data ( angly et al .
, 2009 ) . for the low estimate , we assumed an average dna content per virion twofold lower than expected ( 25  kb , double - stranded dna ) , and conservatively assumed that any rna sequences with no significant match in genbank were not viral ( which reduced the total mass of viral rna by about half ) .
the samples for this study were collected from a pier in the southern portion of kne'ohe bay on the windward side of o'ahu , hawaii ( 21 25 46.80 n , 157 47 31.51 w ) .
this tropical embayment is characterized by year - round warm temperatures ( 22  c to 28  c ) and salinities ranging from 32 to 35 , except during periods of heavy rain , when freshwater plumes from stream runoff can transiently suppress salinity to < 30 ( drupp et al . , 2011 ) .
concentrations of chlorophyll a in surface waters are typically low with average values recorded over a 2.5-year period ranging from 0.5 to 1.1  g  l in summer and from 1.3 to 3.4  g 
l in the winter depending on location ( drupp et al . , 2011 ) . following heavy rainfall ,
blooms have been reported with transient increases in chlorophyll a to ca 6  g 
l ( hoover et al . , 2006 ; de carlo et al . , 2007 )
seawater was collected in acid - washed , polycarbonate carboys on 1 august 2009 ( 35 liters ) and on 3 june 2010 ( 80 liters ) during non - storm conditions .
seawater was transported immediately to the laboratory ( < 1  h ) and filtered through 0.22  m pore - size filters ( sterivex , millipore , billerica , ma , usa ) . the sample from june
viruses in the filtrates were concentrated by iron flocculation ( john et al . , 2011 ) then concentrated further in a centrifugal 30  kda ultrafiltration device ( amicon 15 , millipore ) .
viruses in each concentrate were purified by a two - step process in cscl buoyant density gradients ( lawrence and steward , 2010 )
. fractions ( ca 0.5  ml each ) were collected from the final gradient for analysis with an auto densi - flow ( labconco , kansas city , mo , usa ) .
the density of each fraction was measured using a micropipet and an analytical balance ( lawrence and steward , 2010 ) .
cscl was exchanged with buffer ( 10  m tris , 1  m edta , ph 8) by centrifugal ultrafiltration ( steward , 2001 ) and total nucleic acids in a subsample of each fraction was extracted using spin columns ( qiaamp minelute viral spin kit , qiagen , valencia , ca , usa ) according to the manufacturer 's instructions .
samples for rna analysis were treated twice with dnase ( turbo dnase , life technologies , carlsbad , ca , usa ) to avoid nonspecific signal .
rna and dna contents of each fraction were measured separately in parallel subsamples by fluorometry ( quant - it dna and rna assays , life technologies ) in a cuvette fluorometer ( td-700 , turner designs , sunnyvale , ca , usa ) .
the rna content of the individual fractions from the june sample was low , so we also measured the rna content of the pooled putative viral fractions to obtain a more accurate estimate .
the nucleic acid masses ( 95% confidence interval ( ci ) ) were calculated by inverse prediction from the standard curves ( zar , 1996 ) .
the fluorometric assay was necessary to achieve the required sensitivity . to ensure that we did not overestimate viral rna with the fluorometric assay , we compared the measurements of the mass of rna obtained from a purified rna virus by fluorometry and by spectrophometry ( model du 800 , beckman coulter , inc . , brea , ca , usa ) .
the virus was a positive - sense , single - stranded , picornavirad that infects a marine diatom and has a genome of 8800 nucleotides ( nt ) ( schvarcz et al . , unpublished results ) .
we also tested for cross - reaction of dna in the rna fluorescence assay with and without the dnase digestion procedure described above ( figure 1 ) .
the fluorometer was calibrated with the quant - it kit standards ( escherichia coli rna ) .
the rna standards were then read as samples along with undigested and digested dna standards and the purified viral rna . apparent rna concentration for each sample as determined by fluorometry
was plotted as a function of the given rna or dna concentrations ( quant - it kit standards ) or , in the case of the viral rna , the concentration determined by absorption at 260  nm ( sambrook and russell , 2001 ) .
model i regression lines were calculated and a two - tailed t - test used to determine if the slopes were significantly different from zero ( zar , 1996 ) .
to determine the composition of the rna in the selected rna peaks , we created and analyzed metagenomes from this material for the august sample and from one of the duplicate subsamples ( replicate 1 ) from june .
purified , dnase - treated rna from the indicated fractions ( figure 2 ) of each sample was pooled and amplified using random priming - mediated sequence - independent single - primer amplification as described previously ( djikeng et al . , 2008 ; culley et al . , 2010 ) .
tests of this method on rna viral genomes indicated that it results in coverage and redundancy similar to the ideal values that are predicted by the lander 
waterman model ( djikeng et al . , 2008 ) , which suggests amplification biases are limited .
sequence libraries from the resulting amplified complementary dna samples were produced by pyrosequencing ( gs flx titanium , 454 life sciences , branford , ct , usa ) .
sequences that were > 2 s.d . from the mean length or had an average phred score of < 15 were discarded using the mg - rast qc pipeline ( http://metagenomics.anl.gov/ ) and sequence - independent single - primer amplification primers were trimmed from the remaining reads . artificial ( or technical ) replicates ( gomez - alvarez et al . , 2009 )
were removed using the online tool mg - rast ( meyer et al . , 2008 ) .
processed sequences were assembled with clc genomics workbench version 5.0 ( clcbio , cambridge , ma , usa ) using global alignment with automatic word and bubble sizes , a minimum contig length of 200 , mismatch , insertion and deletion costs set to 3 , length fraction set to 0.5 , and the similarity threshold set to 0.8 .
the community composition of each sample was analyzed with megan ( huson et al . ,
2011 ) using the output from blastx comparisons ( altschul et al . , 1990 ) of the assembled metagenomes ( contigs plus singletons ) with the non - redundant ncbi sequence database .
the taxonomic assignment for a given contig was applied to all reads comprising that contig .
megan assignments were manually checked and some reassignments made in the case of annotation errors . as an independent check on possible contamination from cellular rna , sequences were compared with the silva database ( pruesse et al . , 2007 ) using blastn to search for ribosomal rna contamination , which is expected to be the major source of cellular rna contamination ( karpinets et al . , 2006 ) .
the mass of rna in the pooled viral fraction for each library was adjusted downward by the percentage of reads identified as being cell derived .
to estimate the number of rna virus genomes , we first summed the rna mass within a narrow buoyant density range ( between the solid vertical lines , figure 2 ) .
although rna viruses can have densities outside of that range , we used this conservative window because we analyzed the rna in only those fractions by metagenomic analysis . using the relative representation of different taxa in the metagenomic library and the average genome sizes for those taxa derived from data in the ninth report of the international committee on the taxonomy of viruses ( king et al . , 2012 ) , we calculated a weighted average rna mass per virion . for our prima facie estimates of the rna virus contribution , we assumed that the sequences having no significant hit to any sequences in genbank ( either directly or by association with other sequences in a contig ) had the same taxonomic distribution as the collection of sequences that did have significant blast hits .
the portion of the total rna determined to be viral was then divided by the weighted average mass of rna per virion to obtain the number of rna - containing viruses .
errors were calculated ( and propagated ) as the 95% ci of the inverse predictions of rna mass from the standard curves . to estimate abundance of dna viruses
, we summed the dna content in the fractions encompassing the main dna peak in the viral buoyant density range ( between the dashed vertical lines , figure 2 ) and converted the total dna mass to numbers of dna viruses assuming an average dna content per virion of 5.5  10  g ( equivalent to 50  kb double - stranded dna ) , an average that was found to be similar in a wide range of environments ( steward et al . , 2000 ) .
this value is conservative compared with some other estimates ( brum , 2005 ; angly et al . , 2009 ) and
we also calculated the contribution of rna viruses using more extreme assumptions to get a sense of how much higher and lower the percentage might be .
for the high estimate , we assumed that the dna per virion was twofold greater ( 100  kb , double - stranded dna ) as calculated for marine viral metagenomic data ( angly et al . , 2009 ) .
for the low estimate , we assumed an average dna content per virion twofold lower than expected ( 25  kb , double - stranded dna ) , and conservatively assumed that any rna sequences with no significant match in genbank were not viral ( which reduced the total mass of viral rna by about half ) .
to achieve the sensitivity necessary to assay the nucleic acids in our cscl - purified viral fractions , we had to employ fluorometry .
we performed simple tests of the rna assay kit to ensure that it would provide reasonable estimates of viral rna concentration , and that our measurements would not be overestimated by cross - reaction from dna . with the fluorometer calibrated to the kit standards
, the standards themselves result in a linear curve with a significant slope of 1.000 as expected ( p<0.001 ) .
the assay of a purified ssrna virus falls very close to the calibration line ( low by 14% ) suggesting that the kit standards ( dilutions of e. coli rna ) are reasonably accurate for quantifying a ssrna viral genome ( figure 1 ) .
dna was found to have a limited cross - reaction with the rna stain with a significant ( p<0.001 ) slope of 0.084 , which suggests a cross - reaction of 8.4% on average .
however , this cross - reaction was effectively removed ( slope statistically indistinguishable from zero ; p>0.5 ) by the dnase treatment that we applied to all of our samples before rna assay ( figure 1 ) .
the majority of dna ( 8390% , depending on the sample ) fell within a range of buoyant densities from 1.33 to 1.53  g  ml with a peak around 1.45  g  ml ( figure 2 ) .
this buoyant density range and the peak location are similar to previous observations for dna - containing viruses from other marine environments ( steward et al . , 2000 ) and is well within the range for all known viruses ( 1.161.6 ; king et al . , 2012 ) .
rna concentrations displayed a local peak at the same density with 3269% of the total in the same viral range .
dna and rna were also found in higher density fractions nearer to , or at , the bottom of the gradient .
as these fractions are outside the known range of viral buoyant densities , data from them were not included in subsequent analyses .
for the metagenomic libraries prepared from rna within the viral density range ( figure 2 ) , the number of quality - controlled , de - replicated reads was 139  801 and 110  140 for the 2009 and 2010 samples , respectively ( table 1 ) .
the majority of reads ( 6978% ) in both libraries formed contigs with maximum lengths of 9378  bp for august and 9480  bp for june .
we found that , after assembly and classification of the sequences in the rna viral metagenomes , roughly half ( 5057% ) of the reads derived from the designated viral fractions were most similar to known eukaryote - infecting rna viruses , with the majority matching positive - sense , ssrna viruses in the order picornavirales ( table 2 ) .
the percentage of sequences matching known double - stranded rna viruses was very small ( 0.021.2% ) .
the calculated weighted average mass of rna per virion was 5.38  10  g in august and 5.25  10  g in june ( table 2 ) , which translate into weighted average genome sizes of 9528  nt and 9301  nt , respectively .
the average of these ( 9414  nt ) is very similar to the maximum contig lengths observed in each library and somewhat larger than the three complete rna virus genomes ( ranging from 4449 to 9212  nt ) assembled from a sample from coastal british columbia ( culley et al . , 2007 ) .
less than half of the reads ( 4143% ) had no significant matches in genbank .
we presume that most of these unidentified sequences are viral as well , because many otherwise unidentifiable sequences were found to assemble with virus - like sequences and the fraction of reads identified as cellular was small ( table 2 ) .
estimates of cellular contamination in each library based on blastn comparisons with the silva database ( 6.1% and 0.7% ) were nearly identical to those based on blastx comparisons with genbank nr database ( 6.7% and 0.8% of the sequences ) .
the reads identified as cell derived by the former method were a subset of those identified by the latter .
the identification of all of the putative ribosomal rna reads by blastx was somewhat surprising , but appears to be a result of mis - annotations of ribosomal rna genes as hypothetical proteins in genbank . as the blastx annotations used for megan were more comprehensive in the assignment of cell - derived reads , the larger percentages from that method were accepted as correct . after adjusting the rna mass estimates to discount the contribution from cellular rna , our estimates put the contribution of rna viruses to the total number of viruses at 3863% ( table 3 ) . by applying more conservative and more liberal assumptions
, we calculated extreme high and low estimates that ranged from 15% to 77% ( table 3 ) .
the absolute numbers for both types of viruses in our final samples are minimum estimates , because of losses during processing . assuming typical viral concentrations for these waters based on epifluorescence microscopy to be on the order of 0.5 to 1  10 per liter , this implies overall final yields of dna viruses on the order of 1733% ( june ) and 59% ( august ) , but these estimates have considerable uncertainty . in the absence of data to the contrary
, we assume these losses to be similar for rna and dna viruses . if dna viruses were preferentially lost , this would have led to overestimates of the contribution of rna viruses .
some of the largest viruses ( especially those > 0.2  m in diameter ) will be lost by the initial 0.2  m filtration and all of the largest viruses known so far contain double - stranded dna ( king et al . ,
2012 ) , suggesting that the purification procedure will have some bias against dna viruses .
however , the numerical contribution of viruses > 0.2  m to the total virioplankton appears to be low based on quantitative surveys using electron microscopy ( bratbak et al .
, 1992 ; cochlan et al . , 1993 ; wommack and colwell , 2000 ; auguet et al . , 2006 ) .
the viral concentration method we used has not been tested specifically on rna viruses , but is reported to be exceptionally efficient for dna viruses ( john et al . , 2011 ) .
we were also more conservative in our assignment of rna as viral compared with dna . from the above considerations , we feel it is unlikely that we have grossly overestimated the contribution of rna viruses , but we can not conclusively rule out the possibility that our procedure is significantly biased against dna viruses .
if rna viruses are as abundant as our data suggest , this would have important consequences for our understanding of marine viral ecology .
the data imply , for example , that eukaryotic viruses can be just as abundant as bacteriophages in coastal ocean waters , despite the fact that the eukaryotic plankton concentrations are orders of magnitude lower than those of bacteria .
it seems that the much larger burst sizes of eukaryotic rna viruses ( thousands to tens - of - thousands ; lang et al . , 2009 )
relative to those of bacteriophages ( tens to hundreds ; wommack and colwell , 2000 ) , compensates for the lower host abundances .
this is consistent with earlier theoretical work based on mass transport calculations ( murray and jackson , 1992 ) .
at the time of that report , only dna viruses with relatively small burst sizes were known to infect marine protists . as a consequence
they point out , however , that small viruses having a large burst size ( that is , those like the rna viruses that we now know exist ) , could have a large contribution to protistan mortality . as rna viruses and single - stranded dna viruses
are not reliably detected with the current routine methods for viral direct counts , our data also suggest that many rate estimates that depend on fluorescence - based viral direct counts , such as viral turnover times , virus  host contact rates and viral production rates , may be in need of revision .
development of new methods to directly count even the smallest viruses would be helpful in better constraining the rates of these important processes . in the meantime
, the approach we described here provides a means to estimate the relative contribution of rna viruses in natural aquatic habitats to determine whether our observations are more broadly representative of marine and freshwater habitats .
 viruses are abundant in the ocean and a major driving force in plankton ecology and evolution . it has been assumed that most of the viruses in seawater contain dna and infect bacteria , but rna - containing viruses in the ocean , which almost exclusively infect eukaryotes , have never been quantified . 
 we compared the total mass of rna and dna in the viral fraction harvested from seawater and using data on the mass of nucleic acid per rna- or dna - containing virion , estimated the abundances of each . 
 our data suggest that the abundance of rna viruses rivaled or exceeded that of dna viruses in samples of coastal seawater . 
 the dominant rna viruses in the samples were marine picorna - like viruses , which have small genomes and are at or below the detection limit of common fluorescence - based counting methods . 
 if our results are typical , this means that counts of viruses and the rate measurements that depend on them , such as viral production , are significantly underestimated by current practices . 
 as these rna viruses infect eukaryotes , our data imply that protists contribute more to marine viral dynamics than one might expect based on their relatively low abundance . 
 this conclusion is a departure from the prevailing view of viruses in the ocean , but is consistent with earlier theoretical predictions . 
Introduction
Materials and methods
Study site
Sample collection and processing
Fraction analysis
Metagenomic analysis
Calculation of DNA and RNA genome copies
Results and Discussion