Previous Article | Next Article ![]()
Applied and Environmental Microbiology, March 2009, p. 1402-1409, Vol. 75, No. 5
0099-2240/09/$08.00+0 doi:10.1128/AEM.01899-08
Copyright © 2009, American Society for Microbiology. All Rights Reserved.

College of Marine Science, University of South Florida, St. Petersburg, Florida,1 Florida Integrated Science Center, U.S. Geological Survey, Tallahassee, Florida2
Received 15 August 2008/ Accepted 30 December 2008
|
|
|---|
|
|
|---|
Since fecal-associated bacteria are able to live in sediments in the absence of fecal pollution (18, 32, 55), their resuspension into the water column can result in false-positive results and mask correlations between their concentrations and the extent of recent fecal pollution. Another unfavorable characteristic of current bacterial indicators is their inability to predict or correlate with the presence of pathogenic viruses (25, 40, 41, 64, 80). Human-pathogenic viruses associated with feces are generally more robust than enteric bacteria and are not as easily eliminated by current methods of wastewater treatment (43, 80). For example, adenoviruses are more resilient to tertiary wastewater treatment and UV disinfection than are bacterial indicators of fecal pollution (74). Since bacterial indicators cannot accurately depict the risks to human health from fecal pollution, several studies have proposed the use of a viral indicator of wastewater contamination (35, 41, 61).
While it is impractical to monitor the presence of all viral pathogens related to wastewater pollution, the development of an accurate viral indicator of sewage contamination is needed for enhanced water quality monitoring. Enteric viruses (including viruses belonging to the families Adenoviridae, Caliciviridae, Picornaviridae, and Reoviridae) are transmitted via the fecal-oral route and are known to be abundant in raw sewage. These viruses have been used to identify fecal pollution in coastal environments throughout the world (27, 35, 39, 40, 48, 50, 56, 57, 63, 64, 67-69, 71, 80). To determine which viruses are effective indicators of fecal pollution, it is first necessary to establish a broad, baseline understanding of the many diverse groups of eukaryotic viruses in raw sewage. Several studies have identified adenoviruses, noroviruses, reoviruses, rotaviruses, and other enteroviruses (e.g., polioviruses, coxsackie viruses, and echoviruses) in raw sewage in Australia, Europe, and South Africa (30, 47, 58, 76-78). However, no broad baseline data on the presence of eukaryotic viruses in raw sewage in the United States currently exist.
This study determined the presence of 10 viral groups (adenoviruses, enteroviruses, hepatitis B viruses, herpesviruses, morbilliviruses, noroviruses, papillomaviruses, picobirnaviruses, reoviruses, and rotaviruses) in raw sewage samples collected throughout the United States. All viral groups that were detected in raw sewage were then examined further to determine if they were also present in final treated wastewater effluent. These 10 viral groups were chosen because of their potential to be transmitted via the fecal-oral route, suggesting that they might be found in raw sewage. Many of these viruses (excluding adenoviruses, enteroviruses, noroviruses, reoviruses, and rotaviruses) have not been studied in sewage despite their likely presence. Picobirnaviruses have been detected in individual fecal samples (12, 70, 79, 82); however, their presence has never been analyzed in collective waste, nor have they been proposed to be potential markers of fecal pollution. This study identified potential viral indicators of fecal pollution and will have important applications to water quality monitoring programs throughout the country.
|
|
|---|
PCR for eukaryotic viral groups.
PCR or reverse transcription-PCR for each of the targeted viral groups was executed using previously published primers and conditions (Table 1). All PCR mixtures had a total volume of 50 µl and contained 2 µl of target DNA, 1x REDTaq PCR buffer (10.0 mM Tris-HCl [pH 8.3], 50.0 mM KCl, 1.1 mM MgCl2, 0.01% gelatin; Sigma-Aldrich, St. Louis, MO), 0.25 mM each deoxynucleoside triphosphate, 1 µM of each primer, and 1 U REDTaq DNA polymerase (Sigma-Aldrich), unless otherwise noted.
|
View this table: [in a new window] |
TABLE 1. Primer nucleotide sequences used to PCR amplify 10 viral groups
|
Enteroviruses.
Nested PCR capable of detecting at least 25 different enteroviruses was used to amplify the 5'-untranslated region of the enterovirus genome (34). An additional 1.8 mM MgCl2 and an additional 1.4 mM MgCl2 were added to the first- and second-round PCR mixtures, respectively. The first-round PCR conditions were 40 cycles of 95°C for 30 s, 57.7°C for 30 s, and 72°C for 45 s, followed by 5 min at 72°C. One microliter of amplified PCR product was added to the second round of PCR, which was amplified by 40 cycles of 95°C for 30 s, 56.5°C for 30 s, and 72°C for 30 s, followed by 5 min at 72°C.
Hepatitis B virus.
Nested PCR was used to amplify the S gene of hepatitis B viruses (53). Both rounds of PCR had an additional 0.4 mM MgCl2 added to the mixture, and the PCR conditions were as follows: 5 min at 95°C, followed by 30 cycles of 95°C for 30 s, 55°C for 40 s, and 72°C for 40 s. Two microliters of product from the first reaction was used as a template in the final PCR.
Herpesviruses.
Nested PCR was used to amplify the DNA polymerase genes of eight human herpesvirus strains (11). The second round of PCR used 5 µl of the first-round product. The conditions for both rounds of PCR were 2 min at 94°C, followed by 55 cycles of 94°C for 20 s, 46°C for 30 s, and 72°C for 30 s, followed by incubation at 72°C for 10 min.
Morbilliviruses.
The phosphoprotein gene of morbilliviruses was amplified by PCR (10). The reaction mixture contained an additional 0.4 mM MgCl2, and reaction conditions were 35 cycles of 94°C for 1.5 min, 25°C for 2 min, and 72°C for 2 min, followed by incubation at 72°C for 7 min.
Noroviruses.
The RNA polymerase genes of noroviruses were amplified by PCR (51). The reaction mixture contained an additional 0.4 mM MgCl2 and 100 µg/ml bovine serum albumin. The PCR conditions were as follows: 3 min at 94°C, followed by 40 cycles of 94°C for 30 s, 49°C for 80 s, and 72°C for 1 min, and a final extension step of 72°C for 10 min.
Papillomaviruses.
The L1 gene of papillomaviruses was amplified by PCR that was capable of detecting 87% of human papillomavirus strains (36). The reaction mixture contained an additional 1.4 mM MgCl2 and was incubated for 10 min at 94°C, followed by 45 cycles of 94°C for 1 min, 50°C for 1 min, and 72°C for 1 min, followed by incubation at 72°C for 5 min.
Picobirnaviruses.
A portion of genomic segment 2 of genotype 1 picobirnaviruses was amplified by PCR (70). The picobirnavirus PCR mixture had an additional 0.9 mM MgCl2 and was incubated for 2 min at 94°C, followed by 40 cycles of 94°C for 1 min, 49°C for 2 min, and 72°C for 3 min and then incubation at 72°C for 7 min.
Reoviruses.
Nested PCR was used to amplify the L1 gene of reoviruses (54). The first PCR mixture contained an additional 1.4 mM MgCl2, and 1 µl of the first-round PCR product was added to the second PCR mixture, which contained an additional 0.9 mM MgCl2 and had a total volume of 25 µl. Both PCR mixtures were incubated for 1 min at 94°C, followed by 35 cycles of 94°C for 20 s, 50°C for 30 s, and 72°C for 30 s and then incubation at 72°C for 10 min.
Rotaviruses.
Nested PCR was used to amplify the VP7 gene of group A rotaviruses (38). The initial PCR mixture had an additional 0.4 mM MgCl2 and was incubated for 1 min at 94°C, followed by 25 cycles of 94°C for 30 s, 55°C for 30 s, and 72°C for 1 min, followed by incubation at 72°C for 3 min. The second reaction mixture, containing an additional 2.4 mM MgCl2 and 2 µg/ml bovine serum albumin, was incubated for 1 min at 94°C, followed by 40 cycles of 94°C for 30 s, 55°C for 30 s, and 72°C for 30 s and then incubation at 72°C for 3 min.
Detection and sequencing of PCR positives.
All PCR products were visualized using agarose gel (2%) electrophoresis with ethidium bromide. The presence of inhibitors in each sample was examined by spiking samples with a positive control and analyzing 1:10 serial dilutions to determine if the sample inhibited amplification. To verify the identity of positive assays, PCR products were sequenced. Positive PCR products with one distinct band were purified using the UltraClean PCR Clean-Up kit (Mo Bio Laboratories, Inc., Carlsbad, CA). If more than one band was present, the UltraClean Gelspin kit (Mo Bio Laboratories, Inc.) was utilized to gel purify the positive PCR product of the correct size. All purified positive PCR products were sequenced with their respective forward primers. Sequences were trimmed using Sequencher (Gene Codes Corporation, Ann Arbor, MI), and the identities of the positive PCR products were confirmed by comparing the sequences against the GenBank nonredundant database using BLASTN (3, 4).
Wastewater viral concentration efficiency.
Known quantities of adenoviruses were added to a raw sewage sample and a final effluent sample to determine the efficiency of the methods employed to concentrate viruses. Twenty milliliters of raw sewage and final effluent was autoclaved for 15 min at 121°C to eliminate adenovirus particles while minimizing the alteration of particulates in the sample. Ten milliliters of these wastewater samples was spiked with adenovirus 20 at a final concentration of 8 x 105 viruses per milliliter. The concentration of adenovirus 20 (strain AV-931; ATCC VR-1097) particles suspended in cell culture supernatant was determined through SYBR gold staining and epifluorescent microscopy (66). No aggregation of adenovirus 20 particles was observed. The remaining (unspiked) 10 ml of autoclaved raw sewage and final effluent served as a negative control throughout the viral concentration, nucleic acid isolation, and viral detection processes. Viruses were concentrated from all four samples as described above. PCR for adenoviruses was performed on 1:10 serial dilutions of DNA extracted from spiked and control samples using nested PCR as described above.
Determination of assay sensitivity.
A standard curve of target DNA was created for each assay to allow determinations of assay sensitivity. Positive PCR products were purified using the UltraClean PCR Clean-Up kit (Mo Bio Laboratories, Inc.), and the amount of nucleic acid in the purified product was quantified using a NanoDrop ND-1000 apparatus (NanoDrop Technologies, Wilmington, DE). In the case of nested PCRs, the positive control from the first round of amplification was used to create the standard curve. The number of targets was back calculated using the measured concentration of DNA and amplicon size using the following equation, where Y is the concentration of DNA measured by the NanoDrop instrument and Z is the length of the amplicon:
![]() |
Characterizing diversity of picobirnaviruses in raw sewage.
In order to gain a deeper understanding of the diversity of genotype I picobirnaviruses in raw sewage, positive PCR products from each state were cloned into pCR4-TOPO (Invitrogen), and transformants were screened for inserts by PCR with primers M13F and M13R. The UltraClean PCR Clean-Up kit (Mo Bio Laboratories, Inc.) was used to purify PCR products of the proper size, which were then sequenced with primer M13F. All sequences were trimmed using Sequencher (Gene Codes Corporation), and the identities of the positive PCR products were confirmed by comparing the sequences against the GenBank nonredundant database using BLASTN (3, 4).
All sequences with insignificant identities (E values of >0.001 and/or identities of less than 50 bp in length) to known picobirnaviruses were excluded from phylogenetic analyses. All of the remaining sequences were dereplicated to 99% sequence identity with gaps using the online program FastGroup II (81). The dereplication process and the representative sequences for each group were visually confirmed. FastGroup II was also used to calculate Chao1 (22, 23) and to perform rarefaction analyses (42, 46). The phylogenetic relationships of the picobirnaviruses detected in raw sewage to known picobirnaviruses were identified through aligning dereplicated raw sewage sequences with known sequences from GenBank using ClustalX v.1.8 (73). The resulting alignment was verified visually for accuracy and then imported into MEGA v. 4 for phylogenetic analyses (72).
The phylogenetic analysis of cloned picobirnavirus PCR products from raw sewage was executed in a manner similar to that used for previous picobirnavirus studies (12, 70). Prior to creating a phylogenetic tree, the average pairwise Jukes-Cantor distance was calculated to determine the appropriateness of creating a neighbor-joining tree. Since the average pairwise Jukes-Cantor distance was less than 1.0 (0.356), a neighbor-joining phylogenetic tree was calculated using a Jukes-Cantor model, and bootstrap analyses were performed with 2,000 replicates. Branches with insignificant bootstrap values (
50%) were condensed to create the final phylogenetic tree.
Nucleotide sequence accession numbers.
The nucleotide sequences of 207 genotype I picobirnavirus amplicons from raw sewage have been deposited in the GenBank database under accession numbers EU938707 to EU938913.
|
|
|---|
100 viruses. Negative (unspiked) controls were processed in parallel to the spiked samples, and the lack of PCR product in these controls ensured the elimination of preexisting adenovirus particles. Although this experiment was only semiquantitative, it demonstrated that the methods used here to concentrate wastewater samples were successful in concentrating viruses without major losses.
Viruses identified in wastewater samples.
The sensitivities of the assays used to identify the viruses under investigation ranged from 1 target to 100,000 targets and were identical to those reported previously (Table 1). The presence of 10 types of viruses was analyzed in raw sewage and final effluent samples from 12 wastewater treatment facilities located in 11 different states (Table 2). Adenoviruses and picobirnaviruses were found in 100% of the raw sewage samples. Enteroviruses were detected in 9 of 12 (75%) raw sewage samples, and noroviruses were detected in 7 of 12 (58%) raw sewage samples. Hepatitis B viruses, herpesviruses, morbilliviruses, papillomaviruses, reoviruses, and rotaviruses were not detected in any of the raw sewage samples; however, it is possible that they were present in concentrations below the detection limit of their assays (Table 1). Since adenoviruses, enteroviruses, noroviruses, and picobirnaviruses were detected in raw sewage, their presence was determined in final effluent samples collected from each wastewater treatment plant at the same time as raw sewage collection. Adenoviruses and picobirnaviruses were identified in 25% and 33% of the final effluent samples, respectively, while enteroviruses and noroviruses were identified in only one final effluent sample each.
|
View this table: [in a new window] |
TABLE 2. Viruses detected in raw sewage and final effluent samples collected from throughout the coastal United States
|
22 products per location) cloned picobirnavirus PCR products (200 bp) of the RNA-dependent RNA polymerase gene were sequenced by Agencourt (Beverly, MA), and 207 (72%) of these sequences had significant identities (E value of
0.001 and identities over 50 bp) to known picobirnaviruses in the GenBank database. The significant BLASTN sequence identities to known human and porcine picobirnaviruses ranged from 83% to 100% over regions ranging from 56 bp to 194 bp. Prior to alignment, a total of 70 groups were created from the 207 raw sewage picobirnavirus sequences using FastGroup II (81), which was set to group sequences with 99% identity with gaps. While the majority of the groups were composed of one or two sequences, seven groups contained three or more sequences. The dereplicated sequences along with the top BLASTN hits from GenBank were aligned with ClustalX v.1.8 (73), and phylogenetic analyses were carried out with MEGA v. 4 (72). Since the average pairwise Jukes-Cantor distance was 0.356, a neighbor-joining tree was created using the Jukes-Cantor model with a bootstrap replication of 2,000. The final phylogenetic tree with insignificant (bootstrap value of <50) branches condensed is displayed in Fig. 1.
![]() View larger version (29K): [in a new window] |
FIG. 1. Condensed neighbor-joining (Jukes-Cantor model) phylogenetic tree of a 200-bp segment of the picobirnavirus RNA-dependent RNA polymerase gene from 207 U.S. raw sewage sequences. Each branch represents a sequence or a group of sequences (99% identical with gaps) depending upon the block color. Every colored block is numbered from top to bottom, and the number corresponds to the adjacent table, which explains the locations represented by each block (Alabama [AL], California [CA], Connecticut [CT], Florida [FL], Florida Keys [FL K], Louisiana [LA], Maine [ME], Maryland [MD], North Carolina [NC], New Jersey [NJ], Oregon [OR], and Washington [WA]). Reference sequences of human and porcine picobirnaviruses from previously published studies are identified by their GenBank accession numbers, and porcine picobirnaviruses are indicated with a pig (8, 9, 13, 20, 59, 70, 79; J. Buesa, R. Montava, C. Tellez, and J. M. Ribes, unpublished data).
|
In addition to dereplicating sequences, FastGroup II (81) was also used to calculate the Chao1 value, a richness estimator (22, 23), and to plot a rarefaction curve, which displays the number of unique sequences versus the number of clones sequenced (42, 46). Chao1 predicted a minimum of 200 unique picobirnaviruses in the raw sewage, but only 70 were sampled in this study; thus, more clones need to be sequenced to adequately describe the diversity of picobirnaviruses. Since the rarefaction curve is not yet approaching an asymptote (data not shown), more clones need to be sequenced in order to ensure a complete analysis of picobirnavirus diversity in raw sewage throughout the United States.
|
|
|---|
Adenoviruses and picobirnaviruses are potentially useful viral indicators of fecal pollution since they were found in 100% of raw sewage samples collected throughout the United States. Since these viruses were also detected in 25% and 33% of final effluent samples, respectively, future studies will need to determine how the presence of these viruses correlates with important human pathogens. Enteroviruses were detected in 75% of the raw sewage samples, and noroviruses were detected in 58% of the raw sewage samples; thus, the use of these viruses as markers of fecal pollution could potentially underestimate the extent of fecal contamination. The results of this study support prior findings regarding the prevalence of adenoviruses and enteroviruses in raw sewage and final effluent (15, 26, 37, 47, 58, 74, 77, 78). Furthermore, the use of adenoviruses and enteroviruses to indicate fecal contamination has been demonstrated by numerous studies in a variety of environments (reviewed in references 35 and 49) and experimentally tested in coastal ecosystems (35, 39, 40, 56, 64, 69, 80). This study is the first to demonstrate the widespread prevalence of picobirnaviruses in raw sewage and to suggest that these viruses can potentially be used as indicators of sewage contamination.
Since many viral types under investigation were not detected in this study, it is important to acknowledge the limitations of this research. It is expected that the viral load in the human population differs for each virus and that the abundance of these viruses fluctuates daily and seasonally in raw sewage. The samples analyzed in this study represent only a single time point at each treatment facility; therefore, it is possible that the types of viruses found could vary if samples were collected at different times of the year. In addition, the presence of viruses in the final effluent will also vary in response to fluctuations in treatment efficiency. Furthermore, the results of this study may have been biased by the differential recovery efficiency of the concentration and nucleic acid isolation methods for different viral groups. For example, viruses adhering to particles may have been lost in the filtration step. It is also possible that inefficient reverse transcription or the degradation of RNA could have skewed the results. Further studies need to be completed to verify the absence of undetected or underdetected viruses.
Similarly, it is possible that undetected viruses were present at concentrations below the assay detection limits on the day of sampling. To detect smaller concentrations of these viruses, whole-genome amplification using GenomiPhi V2 (GE Healthcare, Piscataway, NJ) was performed according to the manufacturer's instructions prior to viral PCR for each raw sewage sample and each final effluent sample as a means of increasing the concentration of target nucleic acid. Whole-genome amplification increases the total amount of DNA up to 600-fold (reviewed in reference 14), thus increasing the sensitivity of each assay by more than 2 orders of magnitude. After executing whole-genome amplification, papillomaviruses were identified in 2 of 12 raw sewage samples, and enteroviruses were detected in an additional 2 raw sewage samples, raising the total rate of enterovirus detection to 11 of 12 samples. This suggests that papillomaviruses and enteroviruses were present in those samples but at a lower abundance.
Diversity of picobirnaviruses in raw sewage.
Picobirnaviruses are currently an unclassified group of viruses even though they have recently been proposed to belong to the family Picobirnaviridae (9). As their name suggests, these viruses have a bisegmented double-stranded RNA genome (21, 70). Picobirnavirus particles are fairly small (35 nm), nonenveloped, and spherical. They have been found in the feces of a wide range of mammals including humans (9, 12, 60, 70, 79). They have not been successfully cultured in the laboratory, and their pathogenesis is unknown (21). These viruses have been implicated as possible enteric pathogens and have occasionally been associated with gastroenteritis (12, 21). This is the first study to examine the presence and diversity of picobirnaviruses in raw sewage collected throughout the United States.
Several important conclusions can be drawn from the phylogenetic analysis of the picobirnavirus sequences from raw sewage. First, it appears that a complete understanding of the sequence diversity of picobirnaviruses remains unknown given the large number of single sequences. This is further supported by the Chao1 and rarefaction analyses. The high level of diversity of genotype I human picobirnaviruses, based on the RNA-dependent RNA polymerase gene, has been previously reported (8, 12). Similarly, extensive diversity has been observed in genotype I porcine picobirnaviruses, and their existence as quasispecies has been postulated (9). To fully understand the sequence diversity of genotype I human picobirnaviruses in raw sewage, it would be necessary to sequence more clones until an asymptote is reached in the rarefaction curve.
While the majority of picobirnavirus sequences from raw sewage grouped with known human sequences, some raw sewage sequences grouped most closely with porcine picobirnaviruses. Since the presence of porcine picobirnaviruses in human raw sewage is doubtful and because strong sequence identity between human and porcine picobirnaviruses was previously observed (9), the close grouping of a few raw sewage picobirnavirus sequences to known porcine picobirnavirus sequences is likely the result of the conserved nature of the RNA-dependent RNA polymerase gene.
Potential of adenoviruses, enteroviruses, noroviruses, and picobirnaviruses as indicators of water quality.
To protect public health, it is necessary to identify a practical method for assessing fecal pollution in recreational waters. This is often accomplished through the detection of indicators, which are used to approximate the presence of pathogens. The current bacterial indicators of fecal contamination are not good indicators of wastewater pollution or human health risk in the marine environment (18, 32, 40, 41, 55, 80). The results of this study suggest that adenoviruses, enteroviruses, noroviruses, and picobirnaviruses have potential as viral indicators of fecal pollution because they were detected in the majority of raw sewage samples collected throughout the United States. Since they were detected in 100% of raw sewage samples, adenoviruses and picobirnaviruses are the most promising viral indicators of fecal pollution. Numerous previous studies (34, 35) have used adenoviruses to monitor water quality; however, this is the first study to suggest picobirnaviruses as possible indicator viruses. To assess the potential of picobirnaviruses as indicators to monitor water quality, future studies will need to determine how the presence, abundance, and stability of these viruses correlate with pathogens of concern throughout the wastewater treatment process and in contaminated coastal environments. It will also be necessary to determine if these viruses occur naturally in the environment in the absence of fecal pollution. Our preliminary data indicated that picobirnaviruses were not naturally present at detectable levels in seawater (data not shown); however, further research is needed.
In conclusion, adenoviruses, enteroviruses, noroviruses, and picobirnaviruses were found in the majority of raw sewage samples collected from 12 wastewater treatment facilities throughout the coastal United States. However, adenoviruses and picobirnaviruses were the only viruses detected in 100% of the raw sewage samples. Adenoviruses and picobirnaviruses were also detected in 25% and 33% of final effluent samples, respectively. While adenoviruses, enteroviruses, and noroviruses are known human pathogens and were proposed as markers of fecal pollution in previous studies, the results of this research demonstrate the potential use of picobirnaviruses as an indicator of fecal pollution. Further research will be needed to determine if these candidate viruses have the necessary characteristics of a microbial water quality indicator.
We give special thanks to the personnel from the wastewater treatment facilities that supplied raw sewage and final effluent samples. We thank Mike Gray (USGS, St. Petersburg, FL) for his assistance in the phylogenetic analysis of picobirnaviruses and John Paul (University of South Florida, St. Petersburg, FL) for his insight, positive controls for norovirus and rotavirus assays, and reviewing the manuscript.
Published ahead of print on 5 January 2009. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»