Volume 20, Number 6—June 2014
Human Polyomavirus 9 Infection in Kidney Transplant Patients
Several human polyomaviruses of unknown prevalence and pathogenicity have been identified, including human polyomavirus 9 (HPyV9). To determine rates of HPyV9 infection among immunosuppressed patients, we screened serum samples from 101 kidney transplant patients in the Netherlands for HPyV9 DNA and seroreactivity. A total of 21 patients had positive results for HPyV9 DNA; positivity rates peaked at 3 months after transplantation, but the highest viral loads were measured just after transplantation. During 18 months of follow-up, HPyV9 seroprevalence increased from 33% to 46% among transplant patients; seroprevalence remained stable at ≈30% in a control group of healthy blood donors in whom no HPyV9 DNA was detected. Further analysis revealed an association between detection of HPyV9 and detection of BK polyomavirus but not of cytomegalovirus. Our data indicate that HPyV9 infection is frequent in kidney transplant patients, but the nature of infection—endogenous or donor-derived—and pathogenic potential of this virus remain unknown.
The Polyomaviridae constitute a family of small DNA viruses that infect a variety of hosts. BK polyomavirus (BKPyV) and JC polyomavirus (JCPyV), discovered in 1971 (1,2), are well-known examples of human polyomaviruses (HPyVs) that cause severe disease in immunocompromised patients. Serologic data have revealed that most polyomaviruses are ubiquitous (3–6). In case of JCPyV and BKPyV, primary infection occurs early in life, without apparent symptoms, and persists throughout life as latent infection in the kidneys, accompanied by occasional virus shedding in urine (7). When immunity is decreased, these viruses can reactivate with detectable viremia and manifestation of disease, which poses a threat to, among others, patients who receive solid-organ transplants. For kidney transplant patients, BKPyV infection is considered the most common viral complication and causes nephropathy and graft loss in 1%–10% of cases if left untreated (8). It is not known what determines the severity of BKPyV infection and whether co-infection is involved in the pathogenesis.
Since 2007, at least 10 novel HPyVs have been discovered (9–20); of these, Merkel cell polyomavirus (MCPyV) and trichodysplasia spinulosa–associated polyomavirus (TSPyV) have been shown to be associated with disease (11,19,21,22). Human polyomavirus 9 (HPyV9), so far without a disease association, was identified in 2011 from a serum sample from a kidney transplant patient (17). Overall seroprevalence of HPyV9 has been determined to be 25% to 50% (23–26).
Because HPyV9 was originally isolated from a kidney transplant recipient (17), we aimed to systematically study the presence of HPyV9 infection in kidney transplant patients and investigate a possible association with the known nephropathogenic BKPyV. We analyzed a cohort of 101 transplant patients who received either a kidney transplant or a simultaneous kidney–pancreas transplant for the appearance of markers for HPyV9 infection during the 18 months after transplantation. We assessed the presence of HPyV9 DNA and IgG seroresponses in serum samples. The HPyV9 findings in the transplant cohort were compared with those obtained for an age- and sex-matched cohort of healthy blood donors. Co-infection with BKPyV was investigated by comparing observed HPyV9 and BKPyV viremia levels in the transplant cohort. For comparative purposes, we also tested for cytomegalovirus (CMV), which, like polyomaviruses, frequently reactivates during immunosuppressive drug use after transplantation.
The cohort study consisted of 101 patients who received kidney (n = 83) or kidney–pancreas (n = 18) transplants during 2002–2004 at Leiden University Medical Center (LUMC), Leiden, the Netherlands (Table 1). This study population is part of a larger prospective European multicenter study designed to investigate the role of human papillomavirus infection in the development of skin cancer in solid-organ transplant patients (27). The study adhered to the Declaration of Helsinki Principles, and the medical ethical committee of the LUMC approved of the study design (Medical Ethical Committee no. P02.111). Participants gave written informed consent.
All patients received induction with interleukin-2 receptor blocker daclizumab (100 mg/d) on the day of transplantation and 10 days after transplantation or basiliximab (40 mg at days 0 and 4), followed by triple therapy with prednisone, tacrolimus, or cyclosporine and mycophenolate mofetil. For kidney transplant patients, the dose of the calcineurin inhibitor (tacrolimus or cyclosporine) was tapered at 6 weeks after transplantation, whereas for kidney–pancreas transplant patients, the calcineurin inhibitor was reduced at 3 months after transplantation.
The time points of serum sample collection and the number of samples collected per time point are summarized in Table 2 and shown in relation to the date of transplantation. The baseline samples were obtained in the days immediately after transplantation (time point 0, T0). To collect the subsequent samples, patients were asked to visit the LUMC outpatient clinic for follow-up sample collection at the preferred time points of 3, 6, 9, 12, and 18 months after transplantation (T3–T18). A total of 58 patients provided a sample at all 6 time points: 31, 6, 4, and 2 patients provided 5, 4, 3, or 1 samples, respectively.
Stored pretransplantation serum samples, if available, were retrieved and tested for HPyV9 DNA (n = 65; Table 2) and antibodies (n = 45, 40 of which were also included in pretransplantation DNA testing). The average dates of obtaining samples for DNA and antibody testing were 8 and 2 days before transplantation, respectively.
To obtain a healthy control population, we analyzed anonymized samples from 87 random unpaid blood donors (Table 2). For each donor, 2 follow-up serum samples were studied, collected 1 year apart (T0 and T12). The donors were matched for age and sex with the transplant patient study population (Table 1).
Viral DNA Detection and Quantification
Total nucleic acids were extracted from 200 µL serum by using the MagNA Pure LC Total Nucleic Acid Isolation Kit–High Performance and MagNA Pure LC Instrument (Roche Diagnostics, Indianapolis, IN, USA). To monitor the quality of DNA extraction and potential PCR inhibition, we added low concentrations of phocine herpesvirus (28) to the lysis buffer. DNA was eluted in a final volume of 100 μL elution buffer, of which 10 μL was used as input for real-time quantitative PCR (qPCR).
Primers and Taqman probes were designed by using Beacon Designer software (Premier Biosoft, Palo Alto, CA, USA). For HPyV9, we used the following primers and probe, located in the viral protein (VP) 1 gene amplifying a product of 109 nt: sense primer 5′- CCTGTAAGCTCTCTCCTTA-3′, antisense primer 5′- CCTGATAAATTCTGACTTCTTC-3′, and probe FAM-5′- CTTGTTCTCTGGTCTTATGCCTCA-3′-BHQ-1. For BKPyV, we used the following primers and probe, located in the VP1 gene amplifying a product of 90 nt: sense primer 5′-GAAAAGGAGAGAGTGTCCAGGG-3′, antisense primer 5′-GAACTTCTACTCCTCCTTTTATTAGT-3′, and probe FAM-5′-CCAAAAAGCCAAAGGAACCC-3′-BHQ1.
The BKPyV qPCR and phocine herpesvirus PCR were duplexed for DNA quality and potential PCR inhibition monitoring. Furthermore, the BKPyV qPCR was validated to detect BKPyV genotypes I–IV. qPCR reactions were performed in a total volume of 50 μL, containing 25 μL HotStart Taq mastermix (QIAGEN, Hilden, Germany), 0.5 µmol/L of each primer, 0.35 µmol/L BKPyV probe or 0.4 µmol/L HPyV9 probe, and 3.5 mmol/L MgCl2. Reactions were performed by using a CFX96 real-time detection system (Bio-Rad, Hercules, CA, USA) with the following cycle conditions: 15 min at 95°C followed by 45 cycles of amplification (30 s at 95°C; 30 s at 60°C for HPyV9 qPCR and 55°C for BKPyV qPCR; 30 s at 72°C). For quantification, a standard of pGEX 5×3 HPyV9 VP1 plasmid (Genscript, Piscataway, NJ, USA) and of a quantified BKPyV-positive urine sample were used. Analytical sensitivity of the HPyV9 and BKPyV qPCRs was ≈10 copies/mL. CMV load was measured as described (29) with minor adjustments: 0.5 µmol/L of each primer, 0.2 µmol/L probe, and HotStart Taq mastermix (QIAGEN) was used with the following cycle conditions: 15 min at 95°C, followed by 45 cycles of amplification (5 s at 95°C, 15 s at 63°C, 15 s at 72°C). Analytical sensitivity of the CMV qPCR was ≈100 copies/mL. On each plate, 3 negative controls were included; these controls tested negative in all PCR assays. PCR results with a cycle threshold ≥40 were considered negative.
HPyV9 DNA Sequencing
HPyV9-positive PCR samples were confirmed by sequencing. PCR products were cloned by using the TOPO TA Cloning Kit (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions and subsequently sequenced. Sequence reactions were performed by using the BigDye Terminator Kit (Applied Biosystems, Foster City, CA, USA) and analyzed on an ABI Prism 3130 Genetic Analyzer (Applied Biosystems).
HPyV9 Serologic Testing
To detect IgG seroresponses against the major capsid protein VP1 of HPyV9, we performed an antibody-binding assay using Luminex xMAP technology (30), as described (26). Briefly, the assay is based on cross-linking of glutathione to casein, which is subsequently coupled to fluorescent polystyrene beads (Bio-Rad). Glutathione S-transferase HPyV9 VP1 fusion protein was affinity purified on the beads. Serum samples were tested in a 1:100 dilution, and VP1-bound antibodies were detected with biotinylated goat anti-human IgG (H+L; Jackson Immuno Research, West Grove, PA, USA), followed by streptavidine-R-phycoerythrin (Invitrogen). Finally, the beads and the phycoerythrin signal were analyzed in a Bio-Plex 100 Analyzer (Bio-Rad), which gave results in median fluorescent intensity (MFI). For background correction, MFI values of glutathione S-transferase alone were subtracted to obtain HPyV9 VP1–specific signals. Quality control was performed on each plate with a serum pool consisting of 4 serum samples that had been analyzed in a 1:4 serial dilution, starting with a dilution of 1:100 up to 1:409,600. Little interplate variance was observed.
Cutoff Value Determination
The cutoff value of the antibody-binding assay was defined on the basis of a group of healthy children 0.5–2 years of age and determined as described by van der Meijden et al. (26). The transplant patients and the healthy blood donors were analyzed in 2 independent antibody-binding assays. HPyV9 cutoff values of 252 MFI and 311 MFI were determined for the transplant patient group and the blood donor group, respectively.
Differences between groups in terms of HPyV9 DNA or seroprevalence were assessed by using the Pearson χ2 or Fisher exact test, as appropriate for population size. Independent Student t tests or analyses of variance were used for comparisons of mean values between groups. Occurrence of HPyV9 viremia at transplantation was calculated by using the Kaplan-Meier method, with the time from transplantation to the next detected HPyV9 DNA as the outcome variable for HPyV9-seronegative and -positive patients at baseline (T0). HPyV9 seroconversion at transplantation was calculated by using the Kaplan-Meier method with the time from transplantation to the next seropositive sample as the outcome variable for HPyV9 viremic and nonviremic patients during follow-up. For all tests, 2-tailed p values ≤0.05 were considered significant. The statistical analyses in this study were performed by using SPSS 20 (IBM, Armonk, NY, USA) and Prism 3 statistical software (GraphPad Software Inc., San Diego, CA, USA).
HPyV9 Viremia in Kidney Transplant Patients
To determine HPyV9 viremia in the transplant patients, we assessed the presence of HPyV9 DNA in the complete sample set. During the 18 months after transplantation, HPyV9 DNA was detected at some point in 21 (20.8%) of the 101 patients (Table 2). No significant difference in the detection of HPyV9 DNA was observed between kidney and kidney–pancreas transplant patients (20.5% vs. 22.2%, respectively). For 3 (3.0%) patients, results were positive for consecutive serum samples; persistent HPyV9 DNA detection throughout the follow-up period was observed for 1 patient. Shortly after transplantation (T0, on average 11 days after transplantation), HPyV9 DNA was detected in 3.0% of patients. Detection of HPyV9 DNA peaked 3 months after transplantation (7.1% positivity) and gradually decreased to 2.5% 18 months after transplantation (Figure 1; Table 2).
Lack of HPyV9 Viremia before Transplantation and in Controls
To explore the possibility that the observed viremia was not related to the transplant and immunosuppression but to the underlying cause of the kidney disease (e.g., diabetes), we retrieved and analyzed pretransplantation serum samples for 65 (64%) of the 101 transplant patients. In addition, the group of 87 healthy blood donors was analyzed for the presence of HPyV9 DNA. No DNA was detected in either of these sample sets. To confirm the HPyV9-specificity of our PCR findings in the transplant patients, we cloned and sequenced 13 of the 27 HPyV9-positive PCR products, 109 nt in length. The results revealed a complete match with the described HPyV9 DNA sequence in GenBank (accession no. NC_015150) for 12 of 13 samples. In 1 sample, a single nonsynonymous nucleotide mismatch was observed (A→G at position 2403), resulting in an I321V amino acid mutation in the VP1 capsid protein.
Peak of HPyV9 Viral Load Immediately after Transplantation
The mean HPyV9 DNA load after transplantation was 137 copies/mL (range 25–530 copies/mL). On average, the highest viral loads were observed immediately after transplantation (Figure 1). The kidney–pancreas transplant patients tended to show higher HPyV9 DNA loads than did kidney transplant patients (mean values of 250 and 135 copies/mL, respectively; Table 2), but this difference was not statistically significant (p = 0.123 by Student t test).
HPyV9 Seroreactivity Increase in Transplant Patients but Not in Controls
HPyV9 seroresponses were analyzed for the complete sample set. At baseline, just after transplantation (T0), 33% of transplant patients were HPyV9 seropositive. This percentage corresponds to the percentages that we measured in healthy blood donors (29%) and in 45 pretransplantation serum samples (31%). However, at 1 year after transplantation (T12), the seropositivity rate for transplant patients rose to 46%. This percentage differed significantly from the rate measured for blood donors, which remained stable at ≈30% during 1 year of follow-up (p = 0.029 by χ2 test) (Figure 2, panel A). In total, 15 (15%) of 101 transplant patients seroconverted during follow-up (Technical Appendix Figure 1); these patients represent 23% (15/66) of the patients who were seronegative at baseline.
The intensity of measured HPyV9 serologic responses also increased after transplantation, whereas HPyV9 seroreactivity in blood donors was lower at baseline and remained low within a comparable follow-up period of 1 year (Figure 2, panel B). Mean HPyV9 seroreactivity in the 1-year follow-up samples was significantly higher for the transplant group than for blood donors (p = 0.008 by Student t test). Further analysis of the transplant population revealed that kidney–pancreas transplant patients in particular were responsible for the observed increase in HPyV9 seroreactivity after transplantation (Figure 2, panels C, D). Kidney–pancreas transplant patients had lower mean seroreactivity at baseline than did kidney transplant patients (Figure 2, panel D); the relative increase of seroreactivity in kidney–pancreas transplant patients was confirmed by analyzing the complete dataset with a mixed model analysis (p = 0.003; data not shown).
No Correlation between HPyV9 Viremia and Seroreactivity
Because HPyV9 DNA detection and seroresponses increased after transplantation, we investigated the correlation between these parameters. Comparable proportions of patients who were HPyV9-seropositive and -seronegative at baseline became HPyV9 DNA-positive during follow-up (6/35 [17%] and 15/66 [23%], respectively; p = 0.510 by χ2 test) (Figure 3, panel A), and the measured mean viral loads were comparable for the 2 groups (169 and 152 copies/mL, respectively; p = 0.798 by Student t test). Furthermore, we analyzed whether the presence of HPyV9 DNA influenced HPyV9 seroreactivity during follow up and found no association (Figure 3, panel B); we also found no association when we compared HPyV9 DNA positivity between high and low seroresponders (above and below median MFI) and seroconverters (data not shown). Stratified analyses for kidney–pancreas and kidney transplant patients did not alter the lack of association. HPyV9 DNA and seroreactivity profiles for patients who seroconverted and/or became viremic during follow-up are shown in the online Technical Appendix.
Correlation between HPyV9 and BKPyV Viremia
Of 541 samples tested, 225 (42%) were BKPyV DNA positive; these samples came from 86 (85%) of the 101 transplant patients. HPyV9 DNA was detected significantly more frequently in BKPyV DNA–positive samples than in BKPyV DNA–negative samples (9.8% vs. 1.6%, respectively; p<0.001 by χ2 test) (Figure 4, panel A). During follow-up, HPyV9 DNA was more often detected in BKPyV DNA–positive patients than in BKPyV-negative patients (23.3% vs. 6.7%%, respectively; p = 0.185 by Fisher exact test) (Figure 4, panel B). Furthermore, we divided BKPyV viremic patients into 2 groups, those with high (>103 copies/mL) and low (<103 copies/mL) BKPyV DNA loads, and found HPyV9 DNA–positive patients were overrepresented among patients with high BKPyV loads (p = 0.001 by Fisher exact test; Figure 4, panel C). For 11 (55%) of 20 co-infected patients, BKPyV viremia coincided with HPyV9 viremia; for 8 (40%), BKPyV viremia preceded HPyV9 viremia.
We additionally assessed the presence of viremia caused by CMV, a herpes virus that is not phylogenetically related to HPyV9 and not particularly related to urinary tract infections but that frequently reactivates during immunosuppressive drug use after transplantation. CMV DNA was detected at some point after transplantation in 27 (27%) of the 101 patients (43 [8%] of 541 tested samples). The proportion of HPyV9 DNA–positive samples was similar among CMV DNA positive and negative samples (7% and 5%, respectively; Figure 4, panel A), and no associations were found when comparing HPyV9 DNA positivity and seropositivity among patients who were negative or positive for CMV DNA (Figure 4, panel B; data not shown).
We systematically assessed the presence of HPyV9 DNA and IgG responses in posttransplantation serum samples from kidney transplant patients. These markers of viremia and seroreactivity were shown to increase after transplantation, indicative of active HPyV9 infection, whereas the levels remained stable in matched healthy blood donors.
HPyV9 viremia was detected in 21% of transplant patients at some point within 18 months after transplantation. Most patients were viremic at a single time point, predominantly 3 months after transplantation. The highest mean viral loads were observed immediately after transplantation and decreased gradually over time, but overall HPyV9 loads were low (25–530 copies/mL). Repeat analysis of the complete sample set of 541 serum samples reconfirmed HPyV9 viremia in the same patients (data not shown). At the same time, reanalysis of our cohort showed that the time of a positive finding sometimes differed within viremic patients, compatible with the idea that the viral loads are generally low in persons with HPyV9 viremia and sometimes fall below the PCR detection limit.
Since the identification of HPyV9 in 2011 (17), one study has reported detection of the virus in blood from 2% of immunosuppressed patients (31), whereas other studies did not find HPyV9 (32,33). These studies did not report the time of sampling in relation to transplantation and immunosuppression. Our data suggest that active HPyV9 infection is particularly found in the first year after transplantation. After 18 months, only 2.5% of our transplant patients were HPyV9 DNA–positive, with a mean viral load of 51 copies/mL. The use of different primer sets (and probes) in different studies, with different specificity and sensitivity for the detection of HPyV9 DNA, hampers an accurate comparison among studies.
We observed a peak in HPyV9 DNA detection and load in the first 3 months after transplantation, which coincides with the highest dose of immunosuppressive medication administered to these patients. HPyV9 DNA was not detected in serum samples from patients before transplantation or in serum samples from healthy blood donors. Taken together, these observations indicate a close relationship between active HPyV9 infection and transplantation and/or immunosuppression. The higher mean viral load detected in patients who received a combined kidney–pancreas transplant might be the result of the more intensified immunosuppressive regime applied to these patients. Alternatively, the underlying cause of kidney failure might predispose patients for more frequent HPyV9 infection: 94% of kidney–pancreas transplant patients had diabetes, compared with only 6% of kidney transplant patients.
During follow-up, HPyV9 seroprevalence significantly increased among transplant patients, from 33% to 46%, but remained stable at ≈30% in a control group of healthy blood donors among whom no HPyV9 DNA was detected. A previous cross-sectional study observed a comparable difference in HPyV9 seroprevalence between kidney transplant patients (65%) and healthy persons (45%) (25).
The detection of HPyV9 DNA and the increase in HPyV9 seroreactivity observed after transplantation could reflect primary infection and reactivation. Polyomavirus infections after transplantation and immunosuppression could result from endogenous reactivation, but proof of this concept is lacking. Infection/reactivation originating from the transplanted organ, as suggested for BKPyV (34), should be considered in the case of HPyV9, especially because HPyV9 viremia and baseline HPyV9 seroreactivity were not correlated in this study. HPyV9 viremia was frequently detected in baseline samples from seronegative transplant patients, which suggests donor-derived infection rather than endogenous reactivation.
Additional analysis of our findings showed that HPyV9 viremia was more prevalent in BKPyV DNA–positive samples and in BKPyV-viremic patients than in their BKPyV DNA–negative equivalents; this association reached statistical significance in BKPyV DNA–positive samples. Stratified analysis revealed that HPyV9 DNA positivity was correlated with high BKPyV load. Taken together, these observations suggest that these related viruses benefit from a joint risk factor present in immunosuppressed kidney transplant patients. The observation that CMV and HPyV9 viremia were not associated, however, suggests that the joint risk factor for the polyomaviruses is not simply explained by immunosuppression, which is a well-known risk factor for CMV.
Although we provide strong evidence for emergence of HPyV9 infection in kidney transplant patients and for association between HPyV9 and BKPyV infection, this study has its limitations. The cohort of kidney transplantation patients we tested was small (n = 101) and was formed >10 years ago. Confirmation of our observations in a more recent and larger cohort will strengthen our findings. Because we were not able to investigate whether donor HPyV9 serostatus correlated with HPyV9 viremia in the recipient, future research might explore the possibility of the donor organ as the source of HPyV9 infection. Furthermore, studies that include urine samples that were not available for our analyses could investigate urinary excretion of HPyV9 in the infected patients and might confirm the epidemiologic correlations we found in serum samples. Finally, the sensitivity of detecting HPyV9 viremic episodes was limited by the 3-month sampling interval. Future studies using a shorter sampling interval would increase the number of measurements.
In conclusion, we identified HPyV9 as an emerging infection in immunosuppressed kidney transplant patients. The observed prevalence of HPyV9 DNA in serum samples (21%) considerably exceeded detection rates of HPyV9 found by others in skin samples of immunocompetent (0.9%) and immuncompromised (2.0%) persons (35), which suggests that HPyV9 causes systemic rather than skin infection. Whether HPyV9 is pathogenic in immunocompromised patients, alone or in concert with the well-known pathogen BKPyV, deserves further study. In this context, it would be worthwhile to investigate the course of BKPyV viremia and development of BK virus–associated nephropathy in HPyV9-positive and -negative patients.
Ms van der Meijden is a PhD student at the Leiden University Medical Center. Her research interests include epidemiologic and molecular characterization of human polyomaviruses.
We thank Eric Claas for designing primers and probes of the HPyV9 and BKPyV PCR and Ann Vossen for fruitful discussions.
This work was supported in part by Dutch Kidney Foundation grant WO 2013.06/13A1D302.
- Gardner SD, Field AM, Coleman DV, Hulme B. New human papovavirus (B.K.) isolated from urine after renal transplantation. Lancet. 1971;297:1253–7.
- Padgett BL, Walker DL, ZuRhein GM, Eckroade RJ, Dessel BH. Cultivation of papova-like virus from human brain with progressive multifocal leucoencephalopathy. Lancet. 1971;297:1257–60 .
- Carter JJ, Paulson KG, Wipf GC, Miranda D, Madeleine MM, Johnson LG, Association of Merkel cell polyomavirus–specific antibodies with Merkel cell carcinoma. J Natl Cancer Inst. 2009;101:1510–22.
- Egli A, Infanti L, Dumoulin A, Buser A, Samaridis J, Stebler C, Prevalence of polyomavirus BK and JC infection and replication in 400 healthy blood donors. J Infect Dis. 2009;199:837–46.
- Kean JM, Rao S, Wang M, Garcea RL. Seroepidemiology of human polyomaviruses. PLoS Pathog. 2009;5:e1000363.
- van der Meijden E, Kazem S, Burgers MM, Janssens R, Bouwes Bavinck JN, de Melker H, Seroprevalence of trichodysplasia spinulosa–associated polyomavirus. Emerg Infect Dis. 2011;17:1355–63 .
- Gjoerup O, Chang Y. Update on human polyomaviruses and cancer. Adv Cancer Res. 2010;106:1–51.
- Hirsch HH, Randhawa P. BK polyomavirus in solid organ transplantation. Am J Transplant. 2013;13(Suppl 4):179–88.
- Allander T, Andreasson K, Gupta S, Bjerkner A, Bogdanovic G, Persson MA, Identification of a third human polyomavirus. J Virol. 2007;81:4130–6 .
- Buck CB, Phan GQ, Raiji MT, Murphy PM, McDermott DH, McBride AA. Complete genome sequence of a tenth human polyomavirus. J Virol. 2012;86:10887.
- Feng H, Shuda M, Chang Y, Moore PS. Clonal integration of a polyomavirus in human Merkel cell carcinoma. Science. 2008;319:1096–100.
- Gaynor AM, Nissen MD, Whiley DM, Mackay IM, Lambert SB, Wu G, Identification of a novel polyomavirus from patients with acute respiratory tract infections. PLoS Pathog. 2007;3:e64.
- Korup S, Rietscher J, Calvignac-Spencer S, Trusch F, Hofmann J, Moens U, Identification of a novel human polyomavirus in organs of the gastrointestinal tract. PLoS ONE. 2013;8:e58021.
- Lim ES, Reyes A, Antonio M, Saha D, Ikumapayi UN, Adeyemi M, Discovery of STL polyomavirus, a polyomavirus of ancestral recombinant origin that encodes a unique T antigen by alternative splicing. Virology. 2013;436:295–303.
- Sauvage V, Foulongne V, Cheval J, Ar GM, Pariente K, Dereure O, Human polyomavirus related to African green monkey lymphotropic polyomavirus. Emerg Infect Dis. 2011;17:1364–70 .
- Schowalter RM, Pastrana DV, Pumphrey KA, Moyer AL, Buck CB. Merkel cell polyomavirus and two previously unknown polyomaviruses are chronically shed from human skin. Cell Host Microbe. 2010;7:509–15.
- Scuda N, Hofmann J, Calvignac-Spencer S, Ruprecht K, Liman P, Kuhn J, A novel human polyomavirus closely related to the African green monkey–derived lymphotropic polyomavirus (LPV). J Virol. 2011;85:4586–90.
- Siebrasse EA, Reyes A, Lim ES, Zhao G, Mkakosya RS, Manary MJ, Identification of MW polyomavirus, a novel polyomavirus in human stool. J Virol. 2012;86:10321–6 .
- van der Meijden E, Janssens RW, Lauber C, Bouwes Bavinck JN, Gorbalenya AE, Feltkamp MC. Discovery of a new human polyomavirus associated with trichodysplasia spinulosa in an immunocompromized patient. PLoS Pathog. 2010;6:e1001024.
- Yu G, Greninger AL, Isa P, Phan TG, Martinez MA, de la Luz SM, Discovery of a novel polyomavirus in acute diarrheal samples from children. PLoS ONE. 2012;7:e49449.
- Kazem S, van der Meijden E, Kooijman S, Rosenberg AS, Hughey LC, Browning JC, Trichodysplasia spinulosa is characterized by active polyomavirus infection. J Clin Virol. 2012;53:225–30.
- Shuda M, Arora R, Kwun HJ, Feng H, Sarid R, Fernandez-Figueras MT, Human Merkel cell polyomavirus infection I. MCV T antigen expression in Merkel cell carcinoma, lymphoid tissues and lymphoid tumors. Int J Cancer. 2009;125:1243–9.
- Nicol JT, Touze A, Robinot R, Arnold F, Mazzoni E, Tognon M, Seroprevalence and cross-reactivity of human polyomavirus 9. Emerg Infect Dis. 2012;18:1329–32 .
- Nicol JT, Robinot R, Carpentier A, Carandina G, Mazzoni E, Tognon M, Age-specific seroprevalence of Merkel cell polyomavirus, human polyomaviruses 6, 7, and 9 and trichodysplasia Spinulosa–associated polyomavirus. Clin Vaccine Immunol. 2013;20:363–8.
- Trusch F, Klein M, Finsterbusch T, Kuhn J, Hofmann J, Ehlers B. Seroprevalence of human polyomavirus 9 and cross-reactivity to African green monkey–derived lymphotropic polyomavirus. J Gen Virol. 2012;93:698–705 .
- van der Meijden E, Bialasiewicz S, Rockett RJ, Tozer SJ, Sloots TP, Feltkamp MC. Different serologic behavior of MCPyV, TSPyV, HPyV6, HPyV7 and HPyV9 polyomaviruses found on the skin. PLoS ONE. 2013;8:e81078.
- Antonsson A, Waterboer T, Bouwes Bavinck JN, Abeni D, de Koning M, Euvrard S, Longitudinal study of seroprevalence and serostability of 34 human papillomavirus types in European organ transplant recipients. Virology. 2013;436:91–9.
- Niesters HG. Quantitation of viral load using real-time amplification techniques. Methods. 2001;25:419–29 .
- Kalpoe JS, Kroes AC, de Jong MD, Schinkel J, de Brouwer CS, Beersma MF, Validation of clinical application of cytomegalovirus plasma DNA load measurement and definition of treatment criteria by analysis of correlation to antigen detection. J Clin Microbiol. 2004;42:1498–504.
- Waterboer T, Sehr P, Michael KM, Franceschi S, Nieland JD, Joos TO, Multiplex human papillomavirus serology based on in situ–purified glutathione S–transferase fusion proteins. Clin Chem. 2005;51:1845–53.
- Csoma E, Sapy T, Meszaros B, Gergely L. Novel human polyomaviruses in pregnancy: higher prevalence of BKPyV, but no WUPyV, KIPyV and HPyV9. J Clin Virol. 2012;55:262–5.
- Imajoh M, Hashida Y, Taniguchi A, Kamioka M, Daibata M. Novel human polyomaviruses, Merkel cell polyomavirus and human polyomavirus 9, in Japanese chronic lymphocytic leukemia cases. J Hematol Oncol. 2012;5:25 .
- Rockett RJ, Sloots TP, Bowes S, O’Neill N, Ye S, Robson J, Detection of novel polyomaviruses, TSPyV, HPyV6, HPyV7, HPyV9 and MWPyV in feces, urine, blood, respiratory swabs and cerebrospinal fluid. PLoS ONE. 2013;8:e62764.
- Hirsch HH, Steiger J, Polyomavirus BK. Lancet Infect Dis. 2003;3:611–23.