INTRODUCTION
Enterically transmitted infections are common and can be caused by a large number of different pathogens, including viruses. Many of these viruses cause subclinical infections, which is why their prevalence is difficult to estimate. Increased urbanization and high frequency of global travel favor the emergence and reemergence of enteric pathogens and outbreaks of gastrointestinal infections.
Most of the fecal-oral-transmitted viruses are highly resistant in the water environment, where they may persist at high levels despite the decontamination processes commonly used for drinking water and sewage treatment (
1,
2). Most of them are nonenveloped viruses, such as norovirus, enterovirus, Aichi virus, parechovirus, hepatitis A and E viruses (HAV and HEV, respectively), astrovirus, rotavirus, and adenovirus (Ad). These viruses cause outbreaks or sporadic cases with a wide range of symptoms from mild to severe gastroenteritis to meningitis, respiratory disease, conjunctivitis, paralysis, or hepatitis.
Human norovirus (HuNoV), a positive-sense RNA virus, is one of the leading gastroenteritis agents and has become an increasing public concern worldwide. It belongs to the
Caliciviridae family and is classified into 5 genogroups (GI to GV), of which GI, GII, and GIV infect humans. While GII causes the majority of reported outbreaks, GI is frequently found in environmental water settings (
3,
71).
Aichi virus, parechovirus, and hepatitis A virus all belong to the
Picornaviridae family and are spread globally. They possess a single-stranded positive-sense RNA genome of 7 to 8 kb. Aichi virus causes gastroenteritis and infects both children and adults. Human parechoviruses infect mainly infants, causing mild or severe gastroenteritis, meningitis, encephalitis, and paralysis (
4–6). Hepatitis A and E viruses may cause sporadic hepatitis as well as large outbreaks. Hepatitis A virus (HAV) is classified into 6 genotypes (I to VI), with somewhat different geographical distributions, with types IA and IB being common in Europe (
7).
Hepatitis E virus (HEV) belongs to the
Hepeviridae family and is a positive-sense RNA virus. HEV infecting humans is classified into four genotypes (
1,
2,
4). Types 1 and 2 infect only humans, and the infection has a high mortality rate among pregnant females and young children (
8). These two types may cause large waterborne outbreaks (
8). Types 3 and 4 have zoonotic spread, and have been identified in samples from many different mammals apart from humans, such as pigs, wild boar, deer, moose, and mongooses (
8,
9).
Astroviruses belong to the
Astroviridae family and are also positive-sense RNA viruses. They cause predominately gastroenteritis in children under the age of 2 years. Elderly persons and immunocompromised individuals are involved in up to 20% of the clinical cases and 0.5 to 15% of the outbreaks (
10).
Rotavirus belongs to the
Reoviridae family. Its genome consists of 11 segments of double-stranded RNA (
11). Globally nearly every child has been infected by rotavirus by the age of 5 years (
12), and it is the most common cause of severe diarrhea among infants. However, it also causes diarrhea among elderly persons. Although rotavirus infection is prevalent worldwide, most deaths from this infection occur in developing countries, mostly due to dehydration (
13).
There are seven species of human adenoviruses (Ads), classified A through G in the
Adenoviridae family. Each species is further classified into different serotypes (
14–16). The genome is a nonsegmented DNA 26 to 45 kbp in length. Gastroenteritis is caused mainly by Ad-F types 40 and 41 and Ad-G type 52. Most of the Ads belonging to the other species are also spread fecally and orally. They may cause respiratory diseases, conjunctivitis, and, more rarely hepatitis, pancreatitis, or encephalitis.
It is possible to detect these viruses in sewage, since each individual sheds around 100 g fecal matter per day (
17,
18), and 10
5 to 10
9 enteric virus particles per gram of stool are released daily from an infected individual (
19).
During the study period, there was an ongoing hepatitis A outbreak caused by two different IB strains in Scandinavia (
20). By 17 March 2013, there were 52 hepatitis A cases infected with IB strains in Scandinavia, and at least 28 of these persons were infected by the outbreak strains (
20). This outbreak is thought to have been caused by contaminated frozen strawberries.
Several techniques have been developed to detect viruses in sewage (
21–24). Many of these have been developed for poliovirus surveillance as part of the WHO polio eradication program (
25). The techniques have also been used to detect noroviruses and other viruses in water (
26–28). In this study, human viruses in sewage were concentrated and detected by PCR to investigate if this technique could be used as an early warning system for incipient outbreaks of enterically transmitted viruses.
MATERIALS AND METHODS
Sewage treatment plant and sampling.
Sewage samples were collected at the wastewater treatment plant Ryaverket (owned and run by Gryaab AB), Gothenburg, Sweden. The plant receives and treats wastewater from a population of 690,000 in the Gothenburg region and industrial wastewater and storm water. While the inflow of wastewater from households is fairly constant during the year, the total flow increases during periods of heavy rain due the inflow of storm water. Such daily variations in the flow during the sampling period are shown in
Fig. 1.
The sampling in the present study encompassed every second week, between 14 January and 15 April 2013 (weeks 3, 5, 7, 9, 11, 13, and 15). Wastewater was sampled over 24 h by a fixed-site sampler (model Maxx SP II; Mess-u. Probenahmetechnik GmbH, Rangendingen, Germany) routinely used at Ryaverket for monitoring and control of the sewage. The volume collected was proportional to the flow of the influent wastewater. The sampler was adjusted to withdraw 30 ml wastewater per cycle from the inlet pipe of the sewage plant. Combined weekly samples were prepared by pooling daily collected samples. All daily samples were stored in a refrigerator at +4°C before weekly pooling, and the pooled samples were kept frozen at −20C before analysis.
Virus concentration.
The concentration of viruses from the weekly sample of wastewater was performed by virus adsorption to milk powder, mainly as previously described (
29). Briefly, viruses in 1,000 ml sewage were adsorbed to 10 g acidified skim milk powder (pH 3.5 [Difco]) by stirring at room temperature for 8 h. The skim milk powder was precipitated by centrifugation at 4,500 ×
g. The pellet was dissolved by addition of 3 ml phosphate buffer (pH 7.5) and further treated with 6 ml 0.25 M glycine buffer (pH 9.5) for 45 min. An additional 10 ml phosphate buffer was added, and centrifugation was performed at 180,000 ×
g for 2 h. The pellets obtained were dissolved in 200 μl phosphate buffer (pH 7.5) and stored at −20°C.
To estimate the efficiency of the technique for virus concentration from sewage, ∼0.5 cm3 of three norovirus genogroup II (GII)-positive fecal samples was dissolved in 2.5 ml physiologic NaCl and added to 1,000 ml of the first sewage sample that was collected during week 3 prior to the addition of skimmed milk powder and concentration of viruses. The added fecal extract had a threshold cycle (CT) value of 14.7 in the PCR assay for norovirus GII, which corresponded to 3.5 ×109 genomes/ml cDNA, based on CT values obtained from dilutions of a known concentration of a pUC57plasmid with all targeted regions cloned into the EcoRV site (pUC57cl; GenScript HK, Ltd., Hong Kong). The fecal extract was negative for the other investigated fecal-oral-transmitted viruses. The values obtained after concentration were compared with those obtained from concentration of virus in sewage without addition of fecal suspensions.
Extraction of nucleic acid.
The DNeasy blood and tissue kit (Qiagen) was used for extraction of nucleic acids from 100 μl of concentrated virus solution from sewage or patient samples according to the manufacturer's instructions. Extraction was also performed directly from 100 μl raw sewage without prior concentration of virus by adsorption to milk powder. The nucleic acids were eluted with 100 μl RNase-free H2O (Sigma) and stored at −70°C.
Synthesis of cDNA.
For cDNA synthesis, 5 μl of extracted nucleic acids was added to 6.6 μl of RNase-free H
2O (Sigma), 1 μl (4 μg/μl) random primer (Roche), 0.4 μl (25 mM) deoxynucleoside triphosphate (dNTP) (Roche), 4 μl 5× First Strand buffer (Invitrogen, Carlsbad, CA), 1 μl dithiothreitol (DTT) (Invitrogen), 1 μl RNase OUT (Invitrogen), and 1 μl Superscript III (Invitrogen), which gave a total reaction mixture volume of 20 μl. For detection of hepatitis A virus (HAV), the antisense primer used in PCR, HA3381AS (
Table 1), was used in the cDNA reaction instead of random primers. The cDNA synthesis was performed at 50°C for 60 min followed by inactivation of the reverse transcriptase at 70°C for 15 min. The cDNA was stored at −20°C until use.
Detection of viruses by real-time qPCR.
Primers for hepatitis A virus (HAV), hepatitis E virus (HEV), rotavirus, adenovirus, Aichi virus, astrovirus, and parechovirus were used at a concentration of 600 nM, and the primers for norovirus detection were used at a concentration of 900 nM. A 200 nM concentration of the probes was used in the quantitative PCR (qPCR) for seven of the viruses. For HEV, 300 nM probe was used in the qPCR. The sequences of the primers and probes used are given in
Table 1. A total volume of 50 μl reaction mixture contained 5 μl of cDNA or diluted plasmid pUC57cl, 25 μl universal master mix containing AmpliTaq Gold DNA polymerase, AmpErase UNG, dNTPs, dUTP, passive reference, and optimized buffer components (Applied Biosystems), as well as the primers and probes at the concentrations given above. The qPCRs were optimized with regard to annealing and extension temperatures and to primer and probe concentrations by using dilution series of plasmid pUC57cl containing the targeted regions. For HAV, HEV, and astrovirus, the qPCR had also been optimized for detection of these viruses in patient samples, and this qPCR has since long been used routinely for diagnosis of these viruses. Attempts were made to perform multiple qPCRs for detection of several viruses simultaneously in the same assay, but since the sensitivity of the assay then decreased 10- to 100-fold, single qPCRs were performed. The qPCR for HAV, HEV, rotavirus, and adenovirus was performed by an initial incubation at 50°C for 2 min and 95°C for 10 min followed by amplification consisting of 45 cycles of 95°C for 15 s and at 55°C for 1 min. The qPCRs for Aichi virus, astrovirus, and parechovirus were performed as for the hepatitis A and E viruses, but with extension at 60°C for 1 min. The qPCRs were performed in a 7300 real-time PCR system (Applied Biosystems). All samples were tested in triplicate with each qPCR.
Calculation of the number of viral genomes and the approximate number of infected individuals shedding 1011 virus particles per day required for detection of viral genomes in sewage.
The detection limit for the different viruses by qPCR was 6 to 10 viral genomes/50-μl reaction mixture, based on dilution series of a known amount of the pUC57cl plasmid with all targeted regions cloned into the EcoRV site.
The CT values obtained for sewage samples in the qPCR were used to calculate the number of viral genomes (Ci) by performing linear regression of the CT values obtained from serial dilutions of pUC57cl in relation to the dilutions of the plasmid.
An infected individual is assumed to excrete between 10
7 and 10
11 norovirus, HAV, enterovirus, and adenovirus particles per day (
19). The number of virus particles excreted daily is, to our knowledge, not known for the other investigated viruses. An infected individual was therefore assumed to excrete similar amounts of all investigated viruses. Calculations of the number of infected persons shedding viruses in this study were based on the maximum amount of virus that is shed per day by a newly infected person, 10
11 virus particles/day, thereby estimating the lowest number of individuals shedding virus into the wastewater.
The number of virus particles that was expected to be present in the sewage on the daily basis from one infected person excreting (Cexp) was calculated for each week according to the following equation: Cexp = 1011/{[Σ(respective daily flow)]/7}.
The number of potentially infected individuals (Ninfected) based on the presence of the respective virus in the sewage was estimated as Ninfected = Ci/Cexp.
Patient samples.
Between the end of December 2012 and mid-May 2013, there were nine patients notified for hepatitis A infection in the county of Västra Götaland. Five of these patients lived in the Gothenburg area, and wastewater from their households passed through Ryaverket. Serum samples from these patients were sent for hepatitis A diagnosis to the Clinical Microbiology Laboratory (CML) at Sahlgrenska University Hospital. During the study period, there was also an outbreak of norovirus in Gothenburg. Throughout the sewage sampling period, a total of 970 fecal samples from patients with gastroenteritis were sent to the CML for diagnosis of norovirus, sapovirus, rotavirus, astrovirus, and adenovirus. No patient was identified with HEV IgM or had detectable HEV RNA during the study period.
PCR amplification of HAV RNA.
For the first PCR, 5 μl of cDNA with HAV-specific primers was added to 34.4 μl RNase-free water (Sigma), 5 μl 10× TaqMan buffer (Applied Biosystems, Foster City, CA), 3 μl 25 mM MgCl
2 (Applied Biosystems), 0.4 μl (25 mM) dNTP (Roche Diagnostics, Bromma, Sweden), 200 nM sense primer, 200 nM antisense primer, and 1 U (5 U/μl)
Taq polymerase in a total 50-μl reaction volume for each PCR. The sequences of the primers covering the N-terminal region of VP1 are given in
Table 1. The PCR was performed by using an initial step at 94°C for 3 min and then 40 cycles at 94°C for 20 s, 54°C for 30 s, and 72°C for 1 min. From the product of this initial PCR, 5 μl was used in a nested PCR, which was performed as described above but with other antisense primers (
Table 1). The amplified fragments were visualized on a 1.5% agarose gel with DNA GelRed (Biotium, Hayward, CA).
Sequencing.
Amplified PCR fragments were purified using the QIAquick PCR purification kit (Qiagen) according to the manufacturer's instructions. The eluted DNA was sequenced with the BigDye terminator cycle sequencing ready reaction kit (Applied Biosystems) using the same primers as in the second PCR. After the cycle sequence reaction, the amplified DNA fragments were precipitated by adding 25 μl 95% ethanol and 1 μl 3 M NaAc to each tube and incubated at room temperature for 10 min prior to centrifugation for 45 min. The supernatant was removed, and the pellet was washed with 70% ethanol. Each pellet was dissolved in 15 μl formamide before loading for separation and detection on an Avanti3130XL genetic analyzer (Applied Biosystems).
Sequence analysis.
The sequences obtained were checked manually against the chromatograms with the program SeqMan in the DNASTAR program package (DNASTAR, Inc., Madison, WI). The corrected sequences were aligned with the corresponding region of 500 sequences obtained from GenBank by use of the ClustalX2 program (
30). After manual correction of the alignment, phylogenetic analysis was performed with the PHYLIP program package, version 3.65 (
31). The F84 algorithm with a gamma correction and a transition/transversion ratio of 3.76 was used for genomic distance determination by the DNADIST program. Phylogenetic trees were constructed using the unweighted pair-group method using arithmetic averages (UPGMA) in the NEIGHBOR program. The trees were visualized by the TreeView program (
32).
Nucleotide sequence accession numbers.
The HAV sequences obtained in this study have been deposited in GenBank under accession no. KM486802 to KM486807.
DISCUSSION
Analysis of seven weekly pooled untreated sewage samples from Ryaverket (a large wastewater treatment plant in Gothenburg) performed in the present study clearly demonstrated the presence of seven different enteric viruses: norovirus, astrovirus, rotavirus, adenovirus, Aichi virus, hepatitis A virus, and hepatitis E virus. The amount of virus genomes per liter sewage of each virus was quantified, and even if there was a high loss of viral particles, with only about 10% of the original amount of viruses detected, all viruses except parechovirus could be identified in the samples. The low level of recovery of virus was determined in one experiment only and needs to be confirmed. The magnitude of viral loss in this study is, however, in accordance with those of other studies on virus detection in wastewater (
21,
22,
33). In addition, the concentration seemed to have removed inhibitors of qPCR present in raw sewage, as has been shown also by others (
21,
22,
33,
34). Since an infected person excretes 10
7 to 10
13 virus particles per day (
19), analysis of incoming sewage may be a useful tool to reveal the presence of and quantify excreted fecal human pathogens and thereby give an estimation of the number of infected persons. Most studies on enteric viruses in sewage only detect the viruses and do not relate the virus sequences identified in sewage to those from patients from the same sampling time and region (
27,
35–42). Several of the techniques for virus detection in sewage have also been developed for detection of poliovirus in line with the WHO polio eradication program (
36,
43–49).
There was a peak in number of detected viral genomes during week 7 for most of the investigated viruses in this study. This was the week for winter school holidays in this region of Sweden. The assays to detect the different virus types were not performed on the same day or even during the same week, which means that the identified elevated amount of virus particles in the wastewater probably did not reflect a systematic technical error but may indicate that there was an influx of persons from other regions who shed virus during these holidays.
Analysis of wastewater for enterically transmitted viruses has a number of advantages as it makes it possible to monitor a large population by analyzing samples collected in one place, since wastewater from hundreds of thousands of households gathers at the same wastewater treatment plant. Wastewater contains enteric viruses excreted from persons who are ill as well as from cases of subclinical infection, which reflects the real magnitude of circulating virus in the community. In addition, virus can be detected before an outbreak occurs, as shown in this study for norovirus, since the virus may be excreted in feces before onset of symptoms (
50), which is 1 to 2 days for norovirus GII and 4 to 5 days for astrovirus, (
51), while for hepatitis A and E viruses, the excretion period is longer and occurs for up to 7 weeks (
8,
52).
Detection of norovirus in sewage has been extensively studied, since it is known to cause large waterborne and seasonal outbreaks, sometimes with severe consequences for vulnerable persons (
3,
26,
53–57). In this study, the amount of norovirus GII in wastewater peaked 2 to 3 weeks before the infection spread in hospital wards and nursing homes, indicating that the outbreak was ongoing at least 2 weeks before most persons severely affected by the norovirus infection came to medical attention.
Genetic comparison of viruses in sewage with virus from clinical sporadic cases and from outbreaks may provide a model for understanding the epidemiology of enteric viral pathogens in the population. This was highlighted by the Scandinavian hepatitis A outbreak, which began in Denmark during autumn 2012 and spread during winter 2012 and spring 2013 throughout Scandinavia, with at least 28 persons infected in 2012 (
20). One of the two strains causing this outbreak was detected in sewage in weeks 7 and 9, together with samples from three patients infected with these strains. Two of the patients were diagnosed with acute hepatitis A in weeks 10 and 19 and were both infected with the strain found in the wastewater from week 7. This sewage sample had a large amount of HAV RNA, indicating that there were probably more infected persons than those identified. The third patient, who was infected with the second outbreak strain, had onset of hepatitis in week 3 and was infected by the strain that was found in the wastewater in week 9. None of the reported hepatitis A patients during the study period was infected with the strain that was identified in the wastewater in week 5 and which resembled a strain originating from Iraq. This suggests that there may have been several patients with subclinical cases of hepatitis A who did not seek health care.
There was a constant low concentration of HEV in all but one of the sewage samples, indicating that this virus circulated in the population during the study period. Since HEV is common among Swedish pigs and wild boars (
58,
59), the strains in sewage need to be sequenced to identify if they derive from infected humans, animals, or contaminated food products. The seroprevalence of anti-HEV is comparably high in Sweden (
60,
61), but there are only 10 to 15 notified cases yearly, indicating that most cases are subclinical. However, hepatitis E may also be underdiagnosed, since the awareness of endemic hepatitis E is still low among treating physicians in Sweden.
The surprisingly large amount of Aichi virus found in this study indicates that Sweden could have as high a prevalence of this virus as has been described for other European countries (
62). However, no routine investigation for Aichi virus in patients with diarrhea has been performed in Sweden, and to date only three cases of infection with Aichi virus have been described (
63). This virus is known to cause diarrhea in adults as well as in children and has been suggested to play an etiological role in gastroenteritis, especially in outbreaks associated with contaminated seafood (
64–67). The finding in this study indicates that patients in Sweden with diarrhea should also be investigated for Aichi virus. Investigation of patient samples and sequencing of the strains identified from both patients and sewage will give a better picture of which strains of Aichi virus are circulating in the area and whether only certain strains will cause disease.
This study illustrates the value of typing strains identified in sewage for the communicable disease officer, since regular monitoring of sewage for these viruses may give an early warning of a possible upcoming outbreak. Typing can also be used as a tool for estimating the burden of infection in a specific community. When these viruses are detected in sewage, it is known that there are several persons infected in the community and that these viruses are circulating in the population. Our results indicated that there were more infected persons excreting viruses than those diagnosed with these infections during the study period. However, the number of infected persons is approximate and provides only an estimate of the extent and weekly variation in the number of people who excrete the viruses in the community. According to an experiment performed in Helsinki, Finland, poliovirus could be detected in sewage if 1 only in 10,000 inhabitants excreted the virus (
68). The approach used in this study with molecular-based qPCR may increase the level of detection of the investigated viruses. However, it may also detect genomes from noninfectious virus particles (
69,
70). Even if only viral fragments are identified, it will still give a good indication as to which viruses circulate in the population and which viruses infect persons that do not seek medical attention. Further testing of the material needs to be done to establish the sensitivity of the method and to find out if the detected viruses are infectious viruses that can escape treatment of the sewage and spread further. Future studies should also include viruses that are excreted in urine and investigate the presence of additional enterically transmitted viruses in order to get a powerful tool that can be used for both epidemiological purposes and early warning of outbreaks.