Previous Article | Next Article 
Applied and Environmental Microbiology, April 2008, p. 2544-2549, Vol. 74, No. 8
0099-2240/08/$08.00+0 doi:10.1128/AEM.02477-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
Tracing of Norovirus Outbreak Strains in Mussels Collected near Sewage Effluents
Nancy P. Nenonen,1*
Charles Hannoun,1
Peter Horal,1
Bodil Hernroth,2 and
Tomas Bergström1
Department of Infectious Diseases/Clinical Virology, University of Gothenburg, Guldhedsgatan 10b, 413 46 Gothenburg, Sweden,1
The Royal Swedish Academy of Sciences, Sven Loven Centre for Marine Sciences, Kristineberg 566, SE-450 34 Fiskebäckskil, Sweden2
Received 2 November 2007/
Accepted 20 February 2008

ABSTRACT
Noroviruses from mussels collected near sewage effluents were
compared with local patient outbreak strains. Sequence analyses
of RNA polymerase-capsid-poly(A)-3' (3.1-kilobase) regions confirmed
the 99.9% similarity between genotype I.1 strains from mussels
and patient strains from recreational-bathing outbreaks, indicating
the potential usefulness of sentinel norovirus mussel studies
in tracing human norovirus contamination of coastal waters.

INTRODUCTION
Noroviruses (NoV) are the most commonly identified cause of
acute nonbacterial gastroenteritis in all age groups in the
developed world (
21). Symptoms are short lived, but virus levels
in patient stools are high (
12), and NoV is excreted directly
into the sewage system (
20). Although real-time studies (
4)
show improved water quality following sewage plant treatment,
NoV breakthrough is common (
24,
30,
32), particularly after
heavy rainfall (
23). The high prevalence of NoV genogroup I
(GI) in wastewaters and treated effluents (
4,
16) raises questions
of environmental stability and resistance to wastewater treatment.
Less is known about the circulation of NoV in other waters (
22,
28,
30), although the importance of viral contamination from
sewage effluents in shellfish growing areas is well recognized
(
2,
18).
The efficient filtering and concentrating capacities of filter-feeding bivalves (1, 11) support the use of these bivalves as sentinels for tracing NoV in polluted waters. Using bivalves to trace recent contamination with human sewage may provide information on NoV dispersal in effluents and coastal waters. To assess this potential, we examined wild mussels gathered close to the effluent plume from the Gothenburg sewage treatment plant, correlating the molecular findings with those of human strains from local outbreaks of acute gastroenteritis.

Mussel samples.
Blue mussels (
Mytilus edulis) were collected from Fotö
Island in the northern archipelago, 13 km from the Rya wastewater
treatment plant (RyaWWTP), the Gothenburg sewage station. Situated
at the mouth of the Göta Álv River, the plant serves
a population of 650,000. Fotö (population, 617) lies directly
in the path of prevailing northerly currents carrying effluent
plume from the RyaWWTP (
7) and was therefore chosen for the
NoV studies with wild mussels (Fig.
1). In November 2004, a
batch of 40 mussels was collected and processed in the laboratory
assigned for bivalve analysis, following the principles of Kwok
and Higuchi (
15). Each mussel was examined individually by emulsifying
the whole digestive gland in 4 ml Hanks medium. Mussel gland
suspensions aliquoted in 500-µl volumes were stored at
–80°C. Suspensions (500 µl) were homogenized
in 1,050 µl RNA III tissue lysis buffer (Roche Diagnostic
GmbH, Germany) with a MagNA lyser instrument (Roche), using
ceramic beads and oscillation (two 60-s cycles at 6,500 rpm).
RNA was extracted from the tissue lysate (three 350-µl
volumes), using a MagNA Pure LC instrument according to the
RNA III protocol (elution volume, 50 µl), giving three
extracts.

Clinical samples.
Stool samples were obtained from symptomatic patients during
a surveillance of community and institutional outbreaks of acute
gastroenteritis in West Götaland. Samples emulsified as
10% suspensions in 4 ml Hanks medium were centrifuged at 700
x g for 7 min. Supernatants (130 µl) were added to MagNA
Pure RNA III lysis buffer (220 µl), and RNA (100 µl)
was extracted with the MagNA Pure LC instrument.

NoV RNA detection and comparative analysis of RdRp.
Screening of RNA extracts was based on nested combined reverse
transcription-PCR (RT-PCR) (
25) of the NoV RNA-dependent RNA
polymerase (RdRp) gene, using the primers (
8,
17,
33,
35) and
amplification programs outlined in Table
1. Serial dilutions
(10
0, 10
–1, and 10
–2) of mussel RNA (10 µl)
were amplified in NVp36/d110 RT-PCR master mixes (40 µl).
Three RNA extracts were assayed. Water controls were included
between samples (
15). Amplified RT-PCR mixtures (5 µl
and 8 µl) were nested into PCR master mixes (45 µl)
containing JV12Y/JV13I primers. The specificities of gel-electrophoresed
328-bp products were confirmed by purification (
25), cycle sequencing
of both strands with a BigDye Terminator kit (version 1.1; Applied
Biosystems, Foster City, CA), and sequence analysis on an Avant-3100
gene analyzer (Applied Biosystems, Foster City, CA). Three aliquots
of each mussel suspension were homogenized and screened for
NoV RdRp RNA.
View this table:
[in this window]
[in a new window]
|
TABLE 1. Primers and RT-PCR conditions used in screening mussel digestive gland and patient stool samples for NoV RNA
|
In diagnostic testing of patient samples, RNA (10 µl)
was added to three RT-PCR master mixes (40 µl). Each contained
one forward primer (NI, NVp36, or NVp69) and the reverse primer
NVp110 (Table
1). When amplicons were detected on gel electrophoresis
of NI/NVp110 or NVp69/NVp110 RT-PCR mixtures, the NVp36/NVp110
RT-PCR mixture (5 µl) was selected for nesting into JV12Y/JV13I
PCR master mixes (45 µl). Amplification gave standard-length
328-bp products for sequence confirmation and phylogenic analysis.
Gel products, when detected in first-round NVp36/NVp110 RT-PCR
mixtures, were sequenced with the respective primers.
NoV was detected in 23 of 40 mussels (57%), as confirmed by RdRp sequence analyses. Eighteen strains from mussels (n = 18) aligned with GI.1 reference sequences from GenBank, with one strain from another mussel aligning as GI.2. NoV GII strains were detected in four other mussels (two GII.2 and two GII.4 strains).
Mussel NoV RdRp sequences from bivalves collected in November were compared with patient outbreak sequences from the same year to trace the possible source of mussel contamination. NoV GI.1 strains from November mussels showed 100% sequence similarity (285 nucleotides [nt]) with GI.1 strains (n = 9) from Lake Delsjö bathers infected in August, earlier that year, during outbreaks associated with recreational bathing in inland lakes (27). More than 300 individuals were affected, and two distinct GI strains were identified in clinical samples examined in August, NoV GI.1 in Delsjö patients and GI.3 in Lake Aspen bathers (Fig. 1). None of the mussel strains aligned with GI.3 strains (n = 17) from Aspen bathers.

Long-fragment amplification.
Comparative sequence analysis of RdRp-capsid-poly(A)-3' (3.1-kb)
regions of NoV GI.1 strains from November mussels and Delsjö
patients was used to trace the origins of mussel NoV strains.
RNA from three mussel and five patient samples testing positive
for RdRp was examined by seminested RT-PCR with primers as shown
in Table
2 (
6,
9,
19,
35). RT-PCR master mixes (
25) were modified
to include 0.4 mM betaine (Sigma-Aldrich), 5 units avian myeloblastosis
virus reverse transcriptase (Promega, Madison, WI), 1 unit AmpliTAQ
DNA polymerase (Applied Biosystems), and NoV primers. Viral
RNA was amplified from mussel extracts by RT-PCR with primers
NVp36/N7610R, followed by two seminested PCRs using primer pairs
NVp36/GV7 and GIFFN/N7610R. RNA extracts from patient samples
were amplified by RT-PCR with primers NVp36/T
25VN-3' and seminested
PCR using NVp36/N7610R (Table
2). Purified gel products were
cycle sequenced with primers giving contiguous sequences in
both directions, as shown in Table
3 (
6,
9,
12,
14,
17,
19,
33-
35).
View this table:
[in this window]
[in a new window]
|
TABLE 2. Primers and RT-PCR conditions used to amplify long-fragment NoV RdRp-capsid-poly(A)-3' regions from mussel and patient samples
|

Gene studies and phylogenetic analyses.
Contiguous sequences of NoV GI.1 strains (
n = 3) from November
mussels showed 99.9% similarity, across 3,082 nt, with GI.1
strains (
n = 5) from Delsjö patients infected the previous
August. Mussel NoV strains differed from human strains at one
position, where a synonymous change (C7444T) was noted (GenBank
accession no. M87661). Identical sequences (3,082 nt) were obtained
from three of the five Delsjö patients; sequences from
the two other patients showed one synonymous change (nt position
7276). NoV sequences from mussels and patients were most similar
to GI.1 reference sequences and showed unique GI.1 deletions
in capsid amino acids 198 to 199, 289 to 299, 347 to 351, 375
to 377, and 442 to 445 compared with sequences from other GI
subgroup strains. These gene study analyses indicate bioaccumulation
of NoV strains of human origin in the mussels and provide evidence
of the tracing of clinical strains in the environment.
Phylogenetic trees of 285- and 3,082-nt sequences were constructed based on neighbor-joining methods by using MacVector 7.2 software. The results confirmed the high similarity between NoV GI.1 strains from November mussels and those from Delsjö patients of August, as shown by clustering of strains in Fig. 2A and B. The GI.1 Delsjö strain showed 96% similarity with patient strain SRSV-KY-89/89/J, from a Japanese oyster outbreak in 1989 (GenBank accession no. L23828). Bootstrapping (1,000 replicates) gave reliable values of >75 on all nodes of constructed trees.

Significance of molecular findings.
We present molecular evidence of contaminant human NoV GI.1
RNA in wild mussels collected close to the effluent plume from
a sewage treatment plant following major outbreaks of NoV GI.1
in the region. Sequence similarity between NoV GI.1 strains
from mussels collected in November and from bathers sampled
in August was high (>99.9%; 3,082 nt). With a 3-month time
scale between the GI.1 Delsjö (August) outbreak and detection
of NoV GI.1 in mussels sampled later that year (November), the
sequence identity indicates the human origin of the mussel NoV
strains. During that time period, heavy rainfall in Gothenburg
was recorded (
5), conditions associated with high concentrations
of NoV in effluents from sewage plants (
20,
32). Although it
is unlikely, we cannot exclude the possibility that illicit
waste dumping from recreational or commercial shipping was a
cause of the gross mussel contamination (
2).
These mussel findings suggest the environmental stability of NoV GI strains in water and bivalves (22, 26, 31). Recent European real-time studies show erratic peaking of GI strains in waters before and after sewage plant treatment, although NoV GI is rarely represented in medical reporting (4, 16, 24). In the present study, NoV GI.1 and GI.3 strains caused extensive outbreaks affecting children and adults, confirming the importance of GI infections in community outbreaks. The subsequent detection of NoV GI.1 (3.1 kb) in Fotö mussels collected more than 3 months after the Delsjö outbreak provides evidence of the stability of GI strains in bivalves and coastal waters.
This stability of the noncultivable NoV, combined with the unique filter-feeding habits of bivalves, may provide epidemiological information in "compacted time." Here, we focused on the potential of bivalve studies in tracing strains from past or recent-past NoV outbreaks circulating in a region, recording events that cannot be traced by direct examination of sewage effluents. Similar sentinel studies of hepatitis A virus in clams from Mozambique showed correlation with strains from patients of South and Northeast African origin (25). Bivalves, therefore, may prove a useful tool in tracing other viral agents.
In assessing molecular findings, the absence of Aspen NoV GI.3 strains from Fotö mussels at first seemed incongruous. This variant, causing outbreaks in neighboring Lerum, had circulated with GI.1 (Delsjö) during August (Fig. 1). However, wastewater from Lerum was processed not at the RyaWWTP but at a smaller plant with an outfall pipe close to the Aspen bathing area. As Lake Aspen lies approximately 34 km from Fotö (Fig. 1), disparity in regional population sizes (Gothenburg, 650,000; Lerum, 36,000) and factors such as distance, dilution, and rapid river flow (10) could explain the absence of Aspen strains in Fotö mussels. The negative NoV GI.3 (Aspen) findings further emphasize the significance of the GI.1 mussel findings from Fotö.
This study focused on tracing clinical strains in the environment, where environmental samples, mussels, were taken 3 months after the clinical samples were obtained from patients affected in extensive NoV outbreaks. To date, nucleotide information on NoV strains detected in shellfish has been limited to short sequences (
400 nt), mainly from bivalves implicated in outbreaks (3, 13, 29). In contrast, here we have succeeded in sequencing 3.1-kb NoV regions from wild mussels collected close to the effluent plume from a sewage station by using long-fragment, classical nested RT-PCR methods in well-controlled systems (15). Resultant gene study analyses of mussel and outbreak strains provide evidence of bioaccumulation of NoV strains of human origin in the mussels.
We conclude that sentinel mussel studies can provide evidence of NoV circulation in coastal waters and contribute to our understanding of NoV epidemiology. Moreover, long-fragment analyses may provide an exploratory tool for studying virus-bivalve interactions.

Nucleotide sequence accession numbers.
The sequences described in this study have GenBank accession
numbers EU085483 to EU085490 and EU085497 to EU085529.

ACKNOWLEDGMENTS
This work was supported by the Swedish International Development
Cooperation Agency, Sida (grants SWE-2003-108 and 709 603),
and the LUA-ALF Foundation (grant 7346), Sahlgren's University
Hospital.
We thank Birgitta Bidefors, Mona Groot-Jensen, and members of the Department of Virus Detection, Virology, Sahlgren's University Hospital, Gothenburg, for skilled support and Jan Mattson from the RyaWWTP for advice on the hydrodynamics of the effluent plume.

FOOTNOTES
* Corresponding author. Mailing address: Department of Infectious Diseases/Clinical Virology, University of Gothenburg, Guldhedsgatan 10b, 413 46 Gothenburg, Sweden. Phone: 46 736 601 740. Fax: 46 31827032. E-mail:
nancy.nenonen{at}microbio.gu.se 
Published ahead of print on 29 February 2008. 

REFERENCES
1 - Bayne, B., J. Widdows, and R. J. Thompson. 1976. Physiological integrations, p. 261-292. In B. Bayne (ed.), Marine mussels: their ecology and physiology. Cambridge University Press, Cambridge, United Kingdom.
2 - Berg, D. E., M. A. Kohn, T. A. Farley, and L. M. McFarland. 2000. Multi-state outbreaks of acute gastroenteritis traced to fecal-contaminated oysters harvested in Louisiana. J. Infect. Dis. 181(Suppl. 2):S381-S386.[CrossRef][Medline]
3 - Boxman, I. L., J. J. Tilburg, N. A. Te Loeke, H. Vennema, K. Jonker, E. de Boer, and M. Koopmans. 2006. Detection of noroviruses in shellfish in the Netherlands. Int. J. Food Microbiol. 108:391-396.[Medline]
4 - da Silva, A. K., J. C. Le Saux, S. Parnaudeau, M. Pommepuy, M. Elimelech, and F. S. Le Guyader. 2007. Evaluation of removal of noroviruses during wastewater treatment, using real-time reverse transcription-PCR: different behaviors of genogroups I and II. Appl. Environ. Microbiol. 73:7891-7897.[Abstract/Free Full Text]
5 - Ewert, L., and L. Elander. 2007. Weather and environment, p. 24. In L. Ewart and L. Elander (ed.), Statistisk Årsbok Göteborg 2008. Göteborg Stadskansliet, Gothenburg, Sweden. (In Swedish.)
6 - Gallimore, C. I., J. S. Cheesbrough, K. Lamden, C. Bingham, and J. J. Gray. 2005. Multiple norovirus genotypes characterised from an oyster-associated outbreak of gastroenteritis. Int. J. Food Microbiol. 103:323-330.[CrossRef][Medline]
7 - Garde, K., I. S. Hansen, F. Mohlenberg, E. K. Rasmussen, and C. Murray. 2005. Assessment of effects of phosphate discharge from Rya wastewater treatment plant 3D modelling. Report for Bohuskustens Vattenvårdsförbund 53105-2005. Dansk Institut for vand og miljö (DHI), Hørsholm, Denmark. (In Swedish.)
8 - Green, J., C. I. Gallimore, J. P. Norcott, D. Lewis, and D. W. Brown. 1995. Broadly reactive reverse transcriptase polymerase chain reaction for the diagnosis of SRSV-associated gastroenteritis. J. Med. Virol. 47:392-398.[Medline]
9 - Green, J., J. Vinje, C. I. Gallimore, M. Koopmans, A. Hale, D. W. Brown, J. C. Clegg, and J. Chamberlain. 2000. Capsid protein diversity among Norwalk-like viruses. Virus Genes 20:227-236.[CrossRef][Medline]
10 - Isaeus, M., P. Stålnacke, P. Stenström, J. Magnusson, and K.-M. Norderhaug. 2005. Report on the effects of reduced phosphate release from Rya wastewater treatment plant. Report no. 4947-2005. Norsk institutt for vannforskning (NIVA), Oslo, Norway. (In Swedish.)
11 - Jørgensen, C. B. 1990. Bivalve filter feeding: hydrodynamics, bioenergetics, physiology and ecology. Olsen and Olsen, Fredensborg, Denmark.
12 - Kageyama, T., S. Kojima, M. Shinohara, K. Uchida, S. Fukushi, F. B. Hoshino, N. Takeda, and K. Katayama. 2003. Broadly reactive and highly sensitive assay for Norwalk-like viruses based on real-time quantitative reverse transcription-PCR. J. Clin. Microbiol. 41:1548-1557.[Abstract/Free Full Text]
13 - Kingsley, D. H., G. K. Meade, and G. P. Richards. 2002. Detection of both hepatitis A virus and Norwalk-like virus in imported clams associated with food-borne illness. Appl. Environ. Microbiol. 68:3914-3918.[Abstract/Free Full Text]
14 - Kojima, S., T. Kageyama, S. Fukushi, F. B. Hoshino, M. Shinohara, K. Uchida, K. Natori, N. Takeda, and K. Katayama. 2002. Genogroup-specific PCR primers for detection of Norwalk-like viruses. J. Virol. Methods 100:107-114.[CrossRef][Medline]
15 - Kwok, S., and R. Higuchi. 1989. Avoiding false positives with PCR. Nature 339:237-238.[CrossRef][Medline]
16 - La Rosa, G., S. Fontana, A. Di Grazia, M. Iaconelli, M. Pourshaban, and M. Muscillo. 2007. Molecular identification and genetic analysis of norovirus genogroups I and II in water environments: comparative analysis of different reverse transcription-PCR assays. Appl. Environ. Microbiol. 73:4152-4161.[Abstract/Free Full Text]
17 - Le Guyader, F., M. K. Estes, M. E. Hardy, F. H. Neill, J. Green, D. W. Brown, and R. L. Atmar. 1996. Evaluation of a degenerate primer for the PCR detection of human caliciviruses. Arch. Virol. 141:2225-2235.[CrossRef][Medline]
18 - Le Guyader, F. S., F. Bon, D. DeMedici, S. Parnaudeau, A. Bertone, S. Crudeli, A. Doyle, M. Zidane, E. Suffredini, E. Kohli, F. Maddalo, M. Monini, A. Gallay, M. Pommepuy, P. Pothier, and F. M. Ruggeri. 2006. Detection of multiple noroviruses associated with an international gastroenteritis outbreak linked to oyster consumption. J. Clin. Microbiol. 44:3878-3882.[Abstract/Free Full Text]
19 - Liu, B. L., I. N. Clarke, E. O. Caul, and P. R. Lambden. 1995. Human enteric caliciviruses have a unique genome structure and are distinct from the Norwalk-like viruses. Arch. Virol. 140:1345-1356.[CrossRef][Medline]
20 - Lodder, W. J., J. Vinje, R. van De Heide, A. M. de Roda Husman, E. J. Leenen, and M. P. Koopmans. 1999. Molecular detection of Norwalk-like caliciviruses in sewage. Appl. Environ. Microbiol. 65:5624-5627.[Abstract/Free Full Text]
21 - Lopman, B. A., D. W. Brown, and M. Koopmans. 2002. Human caliciviruses in Europe. J. Clin. Virol. 24:137-160.[CrossRef][Medline]
22 - Maunula, L., I. T. Miettinen, and C. H. von Bonsdorff. 2005. Norovirus outbreaks from drinking water. Emerg. Infect. Dis. 11:1716-1721.[Medline]
23 - Miossec, L., F. Le Guyader, L. Haugarreau, and M. Pommepuy. 2000. Magnitude of rainfall on viral contamination of the marine environment during gastroenteritis epidemics in human coastal population. Rev. Epidemiol. Sante Publique 48(Suppl. 2):2S62-2S71. (In French.)[Medline]
24 - Myrmel, M., E. M. Berg, B. Grinde, and E. Rimstad. 2006. Enteric viruses in inlet and outlet samples from sewage treatment plants. J. Water Health 4:197-209.[Medline]
25 - Nenonen, N. P., B. Hernroth, A. A. Chauque, C. Hannoun, and T. Bergstrom. 2006. Detection of hepatitis A virus genotype IB variants in clams from Maputo Bay, Mozambique. J. Med. Virol. 78:896-905.[CrossRef][Medline]
26 - Rzezutka, A., and N. Cook. 2004. Survival of human enteric viruses in the environment and food. FEMS Microbiol. Rev. 28:441-453.[Medline]
27 - Sartorius, B., Y. Andersson, I. Velicko, B. De Jong, M. Lofdahl, K. O. Hedlund, G. Allestam, C. Wangsell, O. Bergstedt, P. Horal, P. Ulleryd, and A. Soderstrom. 2007. Outbreak of norovirus in Vastra Gotaland associated with recreational activities at two lakes during August 2004. Scand. J. Infect. Dis. 39:323-331.[CrossRef][Medline]
28 - Schaub, S. A., and R. K. Oshiro. 2000. Public health concerns about caliciviruses as waterborne contaminants. J. Infect. Dis. 181(Suppl. 2):S374-S380.[CrossRef][Medline]
29 - Shieh, Y., S. S. Monroe, R. L. Fankhauser, G. W. Langlois, W. Burkhardt III, and R. S. Baric. 2000. Detection of norwalk-like virus in shellfish implicated in illness. J. Infect. Dis. 181(Suppl. 2):S360-S366.[CrossRef][Medline]
30 - Ueki, Y., D. Sano, T. Watanabe, K. Akiyama, and T. Omura. 2005. Norovirus pathway in water environment estimated by genetic analysis of strains from patients of gastroenteritis, sewage, treated wastewater, river water and oysters. Water Res. 39:4271-4280.[Medline]
31 - Vainio, K., and M. Myrmel. 2006. Molecular epidemiology of norovirus outbreaks in Norway during 2000 to 2005 and comparison of four norovirus real-time reverse transcriptase PCR assays. J. Clin. Microbiol. 44:3695-3702.[Abstract/Free Full Text]
32 - van den Berg, H., W. Lodder, W. van der Poel, H. Vennema, and A. M. de Roda Husman. 2005. Genetic diversity of noroviruses in raw and treated sewage water. Res. Microbiol. 156:532-540.[Medline]
33 - Vennema, H., E. de Bruin, and M. Koopmans. 2002. Rational optimization of generic primers used for Norwalk-like virus detection by reverse transcriptase polymerase chain reaction. J. Clin. Virol. 25:233-235.[CrossRef][Medline]
34 - Vinje, J., R. A. Hamidjaja, and M. D. Sobsey. 2004. Development and application of a capsid VP1 (region D) based reverse transcription PCR assay for genotyping of genogroup I and II noroviruses. J. Virol. Methods 116:109-117.[CrossRef][Medline]
35 - Wang, J., X. Jiang, H. P. Madore, J. Gray, U. Desselberger, T. Ando, Y. Seto, I. Oishi, J. F. Lew, K. Y. Green, et al. 1994. Sequence diversity of small, round-structured viruses in the Norwalk virus group. J. Virol. 68:5982-5990.[Abstract/Free Full Text]
Applied and Environmental Microbiology, April 2008, p. 2544-2549, Vol. 74, No. 8
0099-2240/08/$08.00+0 doi:10.1128/AEM.02477-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
This article has been cited by other articles:
-
Rosario, K., Symonds, E. M., Sinigalliano, C., Stewart, J., Breitbart, M.
(2009). Pepper Mild Mottle Virus as an Indicator of Fecal Pollution. Appl. Environ. Microbiol.
75: 7261-7267
[Abstract]
[Full Text]
-
Symonds, E. M., Griffin, D. W., Breitbart, M.
(2009). Eukaryotic Viruses in Wastewater Samples from the United States. Appl. Environ. Microbiol.
75: 1402-1409
[Abstract]
[Full Text]