Human polyomaviruses JC and BK (JCPyV and BKPyV) are two members of the Polyomaviridae family that persistently infect humans and cause disease in immunocompromised individuals. These viruses have been potentially implicated in certain cancers 12. Both respiratory and oral routes have been postulated for their transmission 345. A high frequency of excretion of JCPyV and BKPyV has been reported, and both viruses have been detected in urban sewage from various geographical areas 67. This observation indicates that they could be transmitted by water or food.

In 2007 and 2008, three new polyomaviruses, KI WU and Merkel cell polyomavirus (KIPyV, WUPyV and MCPyV), were reported in humans 8910. KIPyV and WUPyV have been detected mainly in respiratory tract specimens from children and also immunocompromised individuals. In 4 continents these viruses showed equivalent prevalence and highly conserved nucleotide sequences. KIPyV and WUPyV have also been co-detected with other viruses in patients with respiratory and, in some cases, gastrointestinal disorders. Both viruses have been detected in feces 1112 and their role in the etiology of respiratory infections has recently been questioned 13.

MCPyV, which has also been described in respiratory secretions 141516, is strongly associated with Merkel cell carcinomas (MCC) 17. This association strongly supports an etiological role for MCPyV in the development of MCC 18. Recent serological data show that KIPyV, WUPyV and MCPyV are prevalent in the healthy population 19.

Antibodies against lymphotropic polyomavirus (LPyV), a virus of simian origin, have been found in human blood samples 1920. Moreover, LPyV has been reported in human peripheral blood from patients with leukoencephalopathies as well as in immunocompromised and healthy subjects 2122.

Here we assessed KIPyV, WUPyV, MCPyV and LPyV in urban wastewater to determine whether these viruses are prevalent in the environment, as reported for JCPyV and BKPyV 7. For this purpose, we performed nested-PCR (nPCR) assays and compared our results with the nucleotide sequences available in data banks. Wastewater samples collected over the last 6 years from a treatment plant processing domestic and industrial wastewater from a population of 175,000 inhabitants were tested for the presence of KIPyV, WUPyV and MCPyV (8 sewage samples) and also for LPyV (13 sewage and 9 biosolid and sludge samples). In addition, 7 samples collected in 2009 from river water used to source a drinking water treatment plant were also analyzed for the presence of KIPyV, WUPyV and MCPyV. The presence of JCPyV and human adenoviruses (HAdVs) was evaluated by quantitative PCR (qPCR) as a control of the procedures applied and as an index of the level of fecal pollution of human origin present in the samples 6.

Urine samples collected from 4 healthy pregnant women were also tested for WUPyV, KIPyV and MCPyV.

Viral particles were concentrated using methods developed in a previous study using JCPyV as a model. Metods were based on: ultracentrifugation and elution of samples with glycine buffer pH 9.5 for sewage 7 and sludge or biosolids 6, glass wool columns filtration and glycine buffer elution for river water 23 and on ultracentrifugation for urine 3. Negative controls were established for each batch of samples. Nucleic acids were extracted with the QIAamp Viral RNA kit (QIAGEN, Inc.). Oligonucleotide primers (Table 1) were designed based on existing polyomaviral sequences and their specificity against other known polyomaviruses (JCPyV, BKPyV, SV40, LPyV) was checked by nPCR. Samples were analyzed by nPCR in final 50-μL reaction volumes. Briefly, 10 μL of the extracted nucleic acids (corresponding to 2 mL of sewage, 2.5 mL of sludge, 1 g of biosolids, 13.5 mL of river water, and 2 mL of urine) and a 10-fold dilution (to prevent enzymatic inhibition) of each nucleic acid extraction were analyzed in a 40-μL reaction mixture containing 1xPCR Buffer, MgCl₂ at 1.5 mM, 0.025 mM of each dNTP, 0.5 μM of primers and 2 units of TaqGold DNA polymerase (Applied Biosystems). After a first-round PCR, 1 μL of the product was added to 49 μL of the nPCR mixture containing the same components as the first-round PCR mixture. The conditions for the first-round and nPCR reaction conditions were as follows: 95°C for 10 min, 30 cycles of 94°C for 60 sec, 60 sec at the corresponding annealing temperature (Table 1) and extension at 72°C for 60 sec. Amplification was completed with a 7-min extension step at 72°C. Amplicons of the expected size were purified (QIAquick PCR purification kit, QIAGEN, Inc) and sequenced (BigDye sequencing kit and ABI Prism 377 genetic analyzer; Applied Biosystems).

Nucleotide sequences were analyzed using the basic BLAST program http://www.ncbi.nlm.nih.gov/BLAST/. Separate areas were used for the diverse steps of the procedures developed; non-template controls were included in each nPCR reaction. HAdV and JCPyV were tested as a control of the procedures applied as well as of the presence of enzymatic inhibitors in the samples.

We processed the samples as 3 separate batches at 3 separate periods of time. The samples showed typical levels of human fecal pollution, as shown by JCPyV and HAdV concentrations (Table 2). KIPyV and WUPyV were present in 1/8 and 2/8 sewage samples respectively while MCPyV was present in 7/8 sewage samples and was the unique newly described human polyomavirus found in the river water (Table 2). MCPyV was also detected in 1/4 urine samples. The VP1 and VP1/VP2/VP3 genes of the MCPyV genome were also amplified and sequenced in 3 sewage samples to confirm the presence of MCPyV genome (Table 2).

Although the detection technique used here was not quantitative, limiting-dilution nPCR experiments showed approximately 10-100 PCR units/mL of sewage for KIPyV, WUPyV and MCPyV. Samples showed positive results only after nPCR but not after the first-round PCR.

DNA cross contamination was ruled out since no viral strains or plasmids with the genomes of the viruses were available, only for LPyV was a plasmid available in the laboratory as positive control; however, all samples were found to be negative for this virus.

We found that the viruses showed a high degree of sequence stability. All but one sequenced MCPyV amplicon were identical and also identical to the reference sequence with GenBank accession number: EU375803, despite their distinct origins (sewage, river water or urine). This observation confirms the high level of conservation of the DNA of these viruses. Only one MCPyV VP1 amplicon showed a nucleotide that differed from the others and from strain EU375803 although it does not produce any change in the derived protein sequence.

The WUPyV amplicon sequenced was identical to reference strain EF444549 while the KIPyV amplicon sequenced showed one nucleotide of difference with reference strain EF127906.

The nucleotide sequences obtained were deposited in GenBank [GenBank: GQ376529 (WUPyV), GQ376528 (KIPyV), GQ376530 (MCPyV TAg region), GQ452776 (MCPyV VP1/VP2-VP3 region) and GQ390249/50 (MCPyV VP1 region)].

None of the 22 sewage, sludge and biosolid samples tested positive for LPyV although typical concentrations of JCPyV and HAdV indicated human fecal contamination (data not shown). The nPCR assay showed a sensitivity of 1-10 genomic copies/reaction when the complete LPyV genome 24 cloned in pBR322 and quantified spectrophotometrically was analyzed by limiting-dilution nPCR. Thus, LPyV was not detected in the tested samples by these methods.

The observation that MCPyV DNA was much more frequently detected than that of KIPyV or WUPyV might reflect that MCPyV is a more prevalent infection or that it is a highly excreted virus.

Our results on MCPyV in urine, urban sewage and river water strongly support the notion that this virus shows an excretion pattern that resembles that of JCPyV and BKPyV. Human excretion of new polyomaviruses, especially MCPyV, may lead to fecal (urine) contamination of water and food.

In this study we did not attempt the in vitro culture of the new polyomaviruses because no cell culture systems for these viruses are available at present. Furthermore, for other human polyomaviruses, such as JCPyV, the regulatory regions of strains excreted in urine present an archetypal structure and are inefficiently cultured.

To our knowledge, this is the first report of the presence of a virus strongly related to human cancer in sewage and river water samples. We propose that the methodology reported here is suitable to study the prevalence, excretion pattern and genetic variability of recently discovered human polyomaviruses in environmental matrices.

The authors declare that they have no competing interests.

SBM coordinated the study, concentrated urine samples and nucleic acid extractions of the urine samples, collaborated in PCR assays, typed the amplicons detected and drafted the manuscript. JRM concentrated the sewage and biosolid samples and performed the nucleic acid extractions; he also collaborated in the PCR analysis and in the sequencing of the resulting amplicons. BC concentrated river water samples and performed nucleic acid extraction of the same samples. AC collaborated in the production of standards for the quantification of HAdV and JCPyV and in the nucleotide sequence comparisons. RG participated in the development of the methodology, conception and coordination of the study and helped to draft the manuscript. All authors read and approved the final manuscript.

SBM is an assistant professor at the Department of Microbiology of the Faculty of Biology, University of Barcelona. Her main research interests are the epidemiology of human and animal polyomaviruses. She addresses their transmission through the environment and their potential as indicators of the presence of human or/and animal fecal contamination.