Abstract
Fragmented data are available on the human polyomaviruses (HPyVs) prevalence in the gastrointestinal tract. Rearrangements in the non-coding control region (NCCR) of JCPyV and BKPyV have been extensively studied and correlated to clinical outcome; instead, little information is available for KIPyV, WUPyV and MCPyV NCCRs. To get insights into the role of HPyVs in the gastrointestinal tract, we investigated JCPyV, BKPyV, KIPyV, WUPyV and MCPyV distribution among hematological patients in concomitance with gastrointestinal symptoms. In addition, NCCRs and VP1 sequences were examined to characterize the strains circulating among the enrolled patients. DNA was extracted from 62 stool samples and qPCR was carried out to detect and quantify JCPyV, BKPyV, KIPyV, WUPyV and MCPyV genomes. Positive samples were subsequently amplified and sequenced for NCCR and VP1 regions. A phylogenetic tree was constructed aligning the obtained VP1 sequences to a set of reference sequences. qPCR revealed low viral loads for all HPyVs searched. Mono and co-infections were detected. A significant correlation was found between gastrointestinal complications and KIPyV infection. Archetype-like NCCRs were found for JCPyV and BKPyV, and a high degree of NCCRs stability was observed for KIPyV, WUPyV and MCPyV. Analysis of the VP1 sequences revealed a 99% identity with the VP1 reference sequences. The study adds important information on HPyVs prevalence and persistence in the gastrointestinal tract. Gastrointestinal signs were correlated with the presence of KIPyV, although definitive conclusions cannot be drawn. HPyVs NCCRs showed a high degree of sequence stability, suggesting that sequence rearrangements are rare in this anatomical site.
Similar content being viewed by others
Introduction
Human polyomaviruses (HPyVs) are small non-enveloped DNA viruses with a circular double-stranded DNA genome divided into three functional regions: early, late and, interposed between these regions, the non-coding control region (NCCR) [1]. HPyVs have a worldwide distribution and each member has a restrict host range. JC virus (JCPyV), isolated in 1971 from the brain of a patient with Hodgkin’s disease, was the first HPyV described [2]. In the same year, BK virus (BKPyV) was isolated from the urine sample of a renal transplant patient [3]. Since then, other HPyVs have been identified including KI polyomavirus (KIPyV) [4] and WU polyomavirus (WUPyV) [5], isolated from respiratory samples, and Merkel cell polyomavirus (MCPyV) identified in Merkel cell carcinoma (MCC) [6]. Following primary infection, which occurs usually in childhood through the respiratory or oral route, HPyVs persist in the body in a latent state [7]. Reactivation is frequent in immunocompromised patients with possible severe consequences [7, 8]. JCPyV can cause progressive multifocal leukoencephalopathy (PML) [9], whereas BKPyV is associated with hemorrhagic cystitis and interstitial nephropathy in transplant patients [10]. MCPyV is associated with MCC, a rare and aggressive skin cancer [6]. Conversely, disease association has not been documented for KIPyV and WUPyV [11]. Both viruses were detected in the respiratory secretions of children with acute respiratory symptoms [4, 5, 11] and rarely in those with gastrointestinal disorders [4, 12].
Prevalence studies detected HPyVs in stool specimens of hospitalized children [13,14,15,16] and JCPyV and BKPyV also in feces of healthy adults [17, 18]. Furthermore, reactivation and gastrointestinal excretion of KIPyV and WUPyV have been reported in immunosuppressed transplant patients [19]. Environmental virology studies detected HPyVs in sewage samples and polluted water supporting the possibility of an oral–fecal route of spread [7, 20,21,22]. Despite these studies, the data on the prevalence of HPyVs in the gastrointestinal tract are still scarce.
Concerning the NCCR’s genetic variability, while rearrangements in the NCCR of JCPyV and BKPyV have been extensively studied and linked to clinical outcomes [23,24,25], little information is available on KIPyV, WUPyV and MCPyV NCCRs [26, 27].
To get insight into the role played by HPyVs in the gastrointestinal tract, we investigated the distribution of JCPyV, BKPyV, KIPyV, WUPyV and MCPyV in a group of patients with hematological disorders. The presence of the viruses was evaluated in concomitance with the gastrointestinal symptoms, and the genetic characterization of the strains circulating among patients was performed analyzing the VP1 and NCCR sequences.
Materials and methods
Patients, sample collection and processing
The study includes 62 transplanted patients with hematological disorders afferent to the University Hospital Tor Vergata (Rome, Italy). Patients were divided into two sub-groups with respect to the disease: a sub-group of 31 thalassemic patients (15 females and 16 males; mean age: 9.25 years; median age: 9 years) and a sub-group of 31 patients affected by acute myeloid leukemia (AML) (16 females and 15 males, mean age: 41.38 years; median age: 38 years). A total of 62 stool samples was collected from April to December 2018. Signed informed consent was obtained in accordance with the Declaratio7n of Helsinki and the Ethic Committee guidelines (RTN-CL-MOD020).
Total DNA was extracted from 0.2 mg of stool using the Stool DNA Isolation Kit following the manufacturer’s instructions (NORGEN BIOTEK, Canada). The extraction and quality of nucleic acid were checked amplifying the human β-globin gene as previously described [28].
Screening of HPyVs DNA by quantitative real-time PCR
Quantitative amplification assay (qPCR) was carried out to detect and quantify JCPyV, BKPyV, KIPyV, WUPyV, MCPyV DNA, using the 7300 Real-Time PCR System (Applied Biosystems, CA, USA) and published protocols [6, 29,30,31]. All samples were tested in triplicate and the number of viral copies was calculated from an external standard curve obtained using serial dilutions at known titer (range: 102–105 copies) of JCPyV pCY/cl1 plasmid (ATCC®VRMC-1™), WUPyV pcDWUER (#37093) and KIPyV pcDKIER (#37094) plasmids, and MCPyV pcDNA6 sTco plasmid (#40201) (Addgene, Cambridge, MA). Quantitative determination of BKPyV DNA was performed using Thermo Scientific™ AcroMetrix™ BKPyV Panel containing intact, encapsidated viral particles (VP1). Standard precautions were taken to prevent contamination during amplification procedures. The lower detection limit of the assay was 10 DNA copies of the target gene per amplification reaction, corresponding to 10 genome equivalents per reaction (10 gEq/reaction).
Amplification, sequencing and analysis of HPyVs NCCRs
All HPyVs DNA-positive samples were subsequently amplified for NCCR following published protocols [21, 32,33,34], on the 9700 GeneAmp® PCR System (AB Applied Biosystems, CA, USA). The amplified products were purified using the MinElute PCR Purification Kit (QIAGEN, Milan, Italy) and sequenced in a dedicated facility (Bio-Fab research s.r.l., Rome, Italy). The obtained sequences were compared to reference sequences deposited in GenBank (AB081613, AB263926, EF127906, EF444549, EU375803). Sequence alignments were performed with ClustalW2 at the European Molecular Biology Laboratory–European Bioinformatics Institute (EMBL-EBI) website using default parameters [35].
Amplification, sequencing and analysis of HPyVs VP1 regions
All HPyVs DNA-positive samples were also subjected to PCR amplification of the VP1 region, following published protocols [4, 5, 34, 36, 37]. The amplicon size was as follows: BKPyV, 327 bp; JCPyV, 215 bp; KIPyV, 207 bp; WUPyV, 369 bp; MCPyV, 589 bp.
The amplified products were purified using the MinElute PCR Purification Kit (QIAGEN, Milan, Italy) and sequenced in a dedicated facility (Bio-Fab research s.r.l., Rome, Italy). BKPyV subtypes/sub-groups and JCPyV genotypes/subtypes were classified based on the single nucleotide polymorphisms found within the amplified VP1 region. Sequence alignment for all HPyVs VP1 isolates were performed with ClustalW2 at the EMBL-EBI website using default parameters [35].
Amplification and sequencing of reference strains
The quality of our sequences was also checked amplifying the VP1 and NCCR regions of the reference strains AB081613, AB263926, EF127906, EF444549, EU375803 cloned into a luciferase reporter plasmid. PCR conditions were the same as for the DNA extracted from stool samples. Amplicons processing and sequences analyses were as reported above.
HPyVs phylogenetic analysis
A phylogenetic tree was generated using Molecular Evolutionary Genetics Analysis (Mega) version 6.0 software program [38] after aligning the VP1 sequences isolated from the patients to those of the reference strains: JCPyV AB081613, BKPyV AB263926, KIPyV EF127906, WUPyV EF444549, MCPyV EU375803. A bootstrap test with 1000 replicates was performed to evaluate the confidence of the branching pattern of the tree.
Statistical analysis
HPyVs detection was summarized by counts and proportions. If continuous variables were normally distributed, they were expressed as mean ± SD, if not, they were expressed by median and range. The χ2 test was performed to evaluate differences in the viral detection among patient groups, while Mann–Whitney U test was applied for non-normally distributed continuous variables to analyze differences between patients. A p value below 0.05 was considered statistically significant.
Results
HPyVs detection
KIPyV, WUPyV, BKPyV, JCPyV and MCPyV DNAs were detected in both thalassemic and leukemic patients, although with different prevalences. Overall, KIPyV DNA was found in 20/62 samples (32.3%), WUPyV DNA in 20/62 (32.3%), BKPyV DNA in 20/62 (32.3%), JCPyV DNA in 29/62 (47%) and MCPyV DNA in 25/62 (40%). qPCR results showed low amounts of viral DNA with an average value of 5 × 102 gEq/mg (95% CI 4.3–5.8) for KIPyV, 1.5 × 102 gEq/mg (95% CI 1–2) for WUPyV, 6.3 × 102 gEq/mg (95% CI 6.2–6.5) for BKPyV, 7 × 102 gEq/mg (95% CI 6.9–7.2) for JCPyV, and 4.5 × 102 gEq/mg (95% CI 3–8) for MCPyV.
HPyVs mono- and co-infections
Considering HPyVs mono-infection, KIPyV DNA was detected in 12/31 thalassemic patients (6 females/6 males) and in 8/31 (4 females/4 males) AML patients. WUPyV DNA was identified in 11/31 thalassemic patients (6 females/5 males) and in 9/31 AML patients (6 females/3 males). BKPyV DNA was detected in 9/31 thalassemic patients (4 females/5 males) and 11/31 AML patients (4 females/7 males). JCPyV DNA was present in 13/31 thalassemic (8 females/5 males) and in 16/31 AML patients (9 females/7 males). Finally, MCPyV DNA was found in 14/31 thalassemic patients (8 females/6 males) and in 11/31 AML patients (6 females/5 males) (Fig. 1). Examining HPyVs co-infection patterns, the most frequent combination was JCPyV/MCPyV, found in nine samples (three from thalassemic patients and six from patients with AML). All co-infections are summarized in Fig. 1.
No significant association was found between the presence of HPyV DNA and age or gender (Table 1). Instead, a significant association was found between the presence of gastrointestinal symptoms and KIPyV infection, both in thalassemic (p = 0.014) and AML patients (p = 0.015) (Table 1).
Sequencing analysis of HPyVs NCCRs
Sequence analysis of 29 JCPyV NCCRs revealed an archetype CY-like NCCR, with the conservative regions A, B, C, D, E and F [39] in all samples. However, single nucleotide variations were observed in two specimens: G95C and G217A in an AML patient, and G108A and G217A in a thalassemic patient. The G95C and G217A involve the cellular transcription-binding sites TAR and NF1, respectively, while the G108A falls outside transcription-binding sites (Fig. 2).
NCCRs amplified from 20 BKPyV-positive samples were compared to the archetype BKPyV NCCR sequence, arbitrarily divided into O, P, Q, R, and S blocks [40, 41]. All BKPyV analyzed regions were identical to the archetype NCCR. In one thalassemic patient and in one AML patient, point mutations at nucleotide positions 4 (R-block) and 18 (S-block) and a thymidine deletion at nucleotide position 50 (S-block) were identified. All mutations were located outside the cellular transcriptional factor-binding sites (Fig. 3). Regarding KIPyV NCCRs sequences, no nucleotide variations were detected in our KIPyV-positive samples. Sequences were identical to Stockholm 60 reference sequence (GenBank: EF127906) [4]. Finally, NCCRs recovered from 20 WUPyV- and 25 MCPyV-positive samples were identical to those deposited in GenBank under the accession numbers EF444549 and EU375803 [5, 6], except for some nucleotide variations already described [34, 42] (data not shown). Overall, sequence similarity with WUPyV and MCPyV reference strains ranged from 99.5 to 100%.
To confirm that HPyVs NCCR mutations did not originate from DNA polymerase-induced mistakes, we also amplified the NCCRs of the reference strains AB081613, AB263926, EF127906, EF444549 and EU375803, and sequenced them. Analysis of the sequences revealed that they were identical to those of the reference strains, indicating that the PCR did not introduce mutations.
Analysis of HPyVs VP1 region
Sequences analysis of KIPyV, WUPyV and MCPyV VP1 amplicons showed 1–3 nucleotides difference with respect to the reference strains EF127906, EF444549, EU375803, although these variations did not produce any amino acid change in the derived protein sequence. Among the 29 JCPyV-positive samples, there was a prevalence of the European genotypes 1A and 1B. Instead, all BKPyV-positive samples contained the BKPyV subtype I/subgroup b-2.
HPyVs phylogenetic analysis
The phylogenetic analysis carried out on the VP1 sequences obtained from the thalassemic and AML patients revealed that the isolates were 99% identical to the reference sequences (JCPyV AB081613, BKPyV AB263926, KIPyV EF127906, WUPyV EF444549, MCPyV EU375803). All isolates of each investigated virus clustered together and with the corresponding reference strains (Fig. 4).
Discussion
Infections in transplanted patients are related to the kinetics of immune reconstitution and are frequently due to reactivation of latent viruses [19]. BKPyV and JCPyV reactivation in hematological transplanted patients are a common observation because of the immunosuppressive status induced by the conditioning regimen and/or immunosuppressive drugs administered after transplant. For this reason, in the post-transplant phase, BKPyV and JCPyV are regularly monitored by qPCR in urine and blood. The presence of both viruses in these compartments, along with a progressive increase in viral load, may be associated with hemorrhagic cystitis [43], or an increased risk of developing PML [7], respectively. While the association with these diseases has been established, few and fragmented data are currently available about the persistence and excretion of BKPyV and JCPyV as well as KIPyV, WUPyV and MCPyv in the gastrointestinal tract. Our data provides evidence that HPyVs are shed in stool specimens of hematological patients, although at low level, supporting the hypothesis that gastrointestinal tract may be a site of HPyVs persistence and reactivation in humans.
There was no association between polyomaviruses shedding in fecal matter and gender as well as HPyVs shedding in fecal matter and age as reported also by Vanchieri et al. [44].
Among the HPyVs investigated, JCPyV and MCPyV were the polyomaviruses most frequently detected. The higher detection rate of the two polyomaviruses compared to the others could be explained by the wide circulation of JCPyV and MCPyV among the populations [45] or by a higher excretion rate of the two HPyVs.
Several studies have examined the prevalence of HPyVs in stool samples with conflicting results. Bergallo et al. [46] reported detection rates for KIPyV (31%) and WUPyV (25%) similar to ours, while a lower prevalence was reported by other authors [47,48,49,50]. Detection rate ranged from 2.7 to 11% for KIPyV [48, 49] and from 0.9 to 14% for WUPyV [47,48,49,50]. Similarly, the prevalence reported in this study for MCPyV, BKPyV and JCPyV differs from that reported in other studies. In our leukemic patients, MCPyV was detected at a higher rate compared to the study of Kantola et al. [49] who investigated the prevalence of MCPyV in stool samples from leukemic children. The difference in the age group might explain these results. Also, the detection rate of BKPyV and JCPyV is much higher than that reported among healthy adults [44], but similar to that found in hematopoietic stem cell transplant patients, strengthening the hypothesis that immunosuppression increases fecal shedding and the gastrointestinal tract may be a site of latency of HPyVs [51].
After transplantation, gastrointestinal symptoms may be caused by infectious agents, drug effects, metabolic conditions or mechanical complications of surgery. In this study, the presence of gastrointestinal signs was statistically correlated with KIPyV infection. Since gastrointestinal complications are often associated with worse outcomes, determining the cause may suggest whether to initiate an antiviral therapy and the duration of treatment. Reactivation of HPyVs might influence treatment efficacy and have a synergistic effect on the reactivation of other viruses. Thus, studying in depth the contemporary presence of HPyVs might be significant and improve the understanding of the pathogenic potential of these viral species in immunocompromised hosts.
NCCR structure organization is important because of the association of the rearranged variants with specific human diseases, such as PML for JCPyV and nephropathy for BKPyV [52,53,54,55]. In light of that, we investigated the sequence variability of the HPyVs NCCRs in the gastrointestinal tract of hematological patients, in effort of advancing our understanding of HPyVs biology. The NCCRs isolated from JCPyV-positive samples showed an archetype structure with nucleotide changes corresponding to known sequence polymorphisms [56,57,58]. A G217 to A nucleotide transition was identified in the NCCR sequences of patients with the European genotype, corroborating that this mutation represents a common feature of the European strains [58]. Similarly to JCPyV, alignment of BKPyV NCCRs showed an archetype structure with point mutations that did not involve any cellular transcriptional factor-binding sites. Alignment of KIPyV, WUPyV and MCPyV NCCRs showed a high degree of sequence stability, suggesting that rearrangements in this anatomical site are probably rare and providing also an explanation for the low viral load measured. In fact, it is well known that NCCR rearrangements can increase viral replication and influence gene expressions and virulence properties [59].
HPyVs VP1 phylogenetic analysis showed that all isolates grouped into a main cluster together with corresponding reference strain. Since the circulation and the genetic evolution of HPyVs were influenced by virus infectivity and/or virus antigenic variability, monitoring amino acid changes could be useful to improve the understanding of the epidemiological and clinical features of HPyVs.
In conclusion, our results add new information on the prevalence, replication and persistence of HPyVs in the gastrointestinal tract, although the duration of HPyVs shedding in fecal matter after primary infection has not been settled yet. Larger studies that will include serial sample collections, to control for intermittent excretion, will be necessary to clarify this point. Moreover, cell culture systems will be useful to study the effects of NCCR variants on viral replication and virulence in this anatomical site [60].
References
Decaprio JA, Garcea RL (2013) A cornucopia of human polyomaviruses. Nat Rev Microbiol 11:264–276
Padgett BL, Zu Rhein GM, Walker DL, Echroade R, Dessel B (1971) Cultivation of papova-like virus from human brain with progressive multifocal leukoencephalopathy. Lancet 1:1257–1260
Gardner SD, Field AM, Coleman DV, Hulme B (1971) New human papovavirus (B.K.) isolated from urine after renal transplantation. Lancet 1:1253–1257
Allander T, Andreasson K, Gupta S, Bjerkner A, Bogdanovic G, Persson MA et al (2007) Identification of a third human polyomavirus. J Virol 81:4130–4136
Gaynor AM, Nissen MD, Whiley DM, Mackay IM, Lambert SB, Wu G et al (2007) Identification of a novel polyomavirus from patients with acute respiratory tract infections. PLoS Pathog 3:e64
Feng H, Shuda M, Chang Y, Moore PS (2008) Clonal integration of a polyomavirus in human Merkel cell carcinoma. Science 319:1096–1100
Pietropaolo V, Prezioso C, Bagnato F, Antonelli G (2018) John Cunningham virus: an overview on biology and disease of the etiological agent of the progressive multifocal leukoencephalopathy. New Microbiol 41:179–186
Costa C, Bergallo M, Astegiano S, Terlizzi ME, Sidoti F, Segoloni GP et al (2008) Monitoring of BK virus replication in the first year following renal transplantation. Nephrol Dial Transpl 23:3333–3336
Koralnik IJ (2006) Progressive multifocal leukoencephalopathy revisited: has the disease outgrown its name? Ann Neurol 60:162–173
Pinto M, Dobson S (2014) BK and JC virus: a review. J Infect 68:S2–S8
Babakir-Mina M, Ciccozzi M, Bonifacio D, Bergallo M, Costa C, Cavallo R (2009) Identification of the novel KI and WU polyomaviruses in human tonsils. J Clin Virol 46:75–79
Ren L, Gonzalez R, Xu X, Li J, Zhang J, Vernet G (2009) WU polyomavirus in fecal specimens of children with acute gastroenteritis, China. Emerg Infect Dis 15:134–135
Vanchiere JA, Carillo B, Morrow AL, Jiang X, Ruiz-Palacios GM, Butel JS (2016) Fecal polyomavirus excretion in infancy. J Pediatric Infect Dis Soc 5:210–213
Venter M, Visser A, Lassauniere R (2009) Human polyomaviruses, WU and KI in HIV exposed children with acute lower respiratory tract infections in hospitals in South Africa. J Clin Virol 44:230–234
Vanchiere JA, Nicome RH, Greer JM, Demmler GJ, Butel JS (2005) Frequent detection of polyomaviruses in stool samples from hospitalized children. J Infect Dis 192:658–664
Li K, Guo J, Zhao R, Xue Y, Chen L, Yang J (2013) Prevalence of 10 human polyomaviruses in fecal samples from children with acute gastroenteritis: a case-control study. J Clin Microbiol 51:3107–3109
Bialasiewicz S, Whiley DM, Lambert SB, Nissen MD, Sloots TP (2009) Detection of BK, JC, WU, or KI polyomaviruses in faecal, urine, blood, cerebrospinal fluid and respiratory samples. J Clin Virol 45:249–254
Vanchiere JA, Abudayyeh S, Copeland CM, Lu LB, Graham DY, Butel JS (2009) Polyomavirus shedding in the stool of healthy adults. J Clin Microbiol 47:2388–2391
Babakir-Mina M, Ciccozzi M, Alteri C, Polchi P, Picardi A, Greco F (2009) Excretion of the novel polyomaviruses KI and WU in the stool of patients with hematological disorders. J Med Virol 81:1668–1673
Bofill-Mas S, Albinana-Gimenez N, Clemente-Casares P, Hundesa A, Rodriguez-Manzano J, Allard A (2006) Quantification and stability of human adenoviruses and polyomavirus JCPyV in wastewater matrices. Appl Environ Microbiol 72:7894–7896
Bofill-Mas S, Pina S, Girones R (2000) Documenting the epidemiologic patterns of polyomaviruses in human populations by studying their presence in urban sewage. Appl Environ Microbiol 66:238–245
Di Bonito P, Della Libera S, Petricca S, Iaconelli M, Accardi L, Muscillo M et al (2015) Frequent and abundant Merkel cell polyomavirus detection in urban wastewaters in Italy. Food Environ Virol 7:1–6
Gosert R, Rinaldo CH, Funk GA, Egli A, Ramos E, Drachenberg CB (2008) Polyomavirus BK with rearranged noncoding control region emerge in vivo in renal transplant patients and increase viral replication and cytopathology. J Exp Med 205:841–852
Olsen GH, Hirsch HH, Rinaldo CH (2009) Functional analysis of polyomavirus BK non-coding control region quasispecies from kidney transplant recipients. J Med Virol 81:1959–1967
Gosert R, Kardas P, Major EO, Hirsch HH (2010) Rearranged JC virus noncoding control regions found in progressive multifocal leukoencephalopathy patient samples increase virus early gene expression and replication rate. J Virol 84:10448–10456
Song XB, Van Ghelue M, Ludvigsen M, Nordbo SA, Ehlers B, Moens U (2016) Characterization of the non-coding control region of polyomavirus KI isolated from nasopharyngeal samples from patients with respiratory symptoms or infection and from blood from healthy blood donors in Norway. J Gen Virol 97:1647–1657
Ajuh ET, Wu Z, Kraus E, Weissbach FH, Bethge T, Gosert R (2018) Novel Human Polyomavirus non-coding control regions differ in bi-directional gene expression according to host cell, large T-antigen expression, and clinically occurring rearrangements. J Virol 92:e02231-17
Saiki RK, Bugawan TL, Horn GT, Mullis KB, Erlich HA (1986) Analysis of enzymatically amplified beta-globin and HLA-DQ alpha DNA with allele-specific oligonucleotide probes. Nature 324:163–166
Bialasiewicz S, Whiley DM, Lambert SB, Gould A, Nissen MD, Sloots TP (2007) Development and evaluation of real-time PCR assays for the detection of the newly identified KI and WU polyomaviruses. J Clin Virol 40:9–14
Prezioso C, Scribano D, Bellizzi A, Anzivino E, Rodio DM, Trancassini M (2017) Efficient propagation of archetype JC polyomavirus in COS-7 cells: evaluation of rearrangements within the NCCR structural organization after transfection. Arch Virol 162:3745–3752
Tremolada S, Delbue S, Larocca S, Carloni C, Elia F, Khalili K (2010) Polymorphisms of the BK virus subtypes and their influence on viral in vitro growth efficiency. Virus Res 149:190–196
Flaegstad T, Sundsfjord A, Arthur RR, Pedersen M, Traavik T, Subramani S (1991) Amplification and sequencing of the control regions of BK and JC virus from human urine by polymerase chain reaction. Virology 180:553–560
Markowitz RB, Thompson HC, Mueller JF, Cohen JA, Dynan WS (1993) Incidence of BK virus and JC virus viruria in human immunodeficiency virus-infected and uninfected subjects. J Infect Dis 167:13–20
Hashida Y, Higuchi T, Matsui K, Shibata Y, Nakajima K, Sano S et al (2018) Genetic variability of the noncoding control region of cutaneous Merkel Cell polyomavirus: identification of geographically related genotypes. J Infect Dis 217:1601–1611
ClustalW2–Multiple Sequence Alignment: http://www.ebi.ac.uk/Tools/msa/clustalw2/
Agostini HT, Yanagihara R, Davis V, Ryschkewitsch CF, Stoner GL (1997) Asian genotypes of JC virus in Native Americans and in a Pacific Island population: markers of viral evolution and human migration. Proc Natl Acad Sci USA 94:14542–14546
Jin L, Gibson PE, Knowles WA, Clewley JP (1993) BK virus antigenic variants: sequence analysis within the capsid VP1 epitope. J Med Virol 39:50–56
Kumar S, Nei M, Dudley J, Tamura K (2008) MEGA: a biologist-centric software for evolutionary analysis of DNA and protein sequences. Brief Bioinform 9:299–306
Yogo Y, Kitamura T, Sugimoto C, Ueki T, Aso Y, Hara K (1990) Isolation of a possible archetypal JC virus DNA sequence from non-immunocompromised individuals. J Virol 64:3139–3143
Moens U, Van Ghelue M (2005) Polymorphism in the genome of non-passaged human polyomavirus BK: implications for cell tropism and the pathological role of the virus. Virology 331:209–231
Yogo Y, Zhong S, Xu Y, Zhu M, Chao Y, Sugimoto C (2008) Conserved archetypal configuration of the transcriptional control region during the course of BK polyomavirus evolution. J Gen Virol 89:1849–1856
Csoma E, Lengyel G, Bányai K, Takács P, Ánosi N, Marton S et al (2018) Study of Karolinska Institutet and Washington University polyomaviruses in tonsil, adenoid, throat swab and middle ear fluid samples. Future Microbiol 13:1719–1730
Gilis L, Morisset S, Billaud G, Ducastelle-Leprêtre S, Labussière-Wallet H, Nicolini FE et al (2014) High burden of BK virus-associated hemorrhagic cystitis in patients undergoing allogeneic hematopoietic stem cell transplantation. Bone Marrow Transpl 49:664–670
Vanchiere JA, Abudayyeh S, Copeland CM, Lu LB, Graham DY, Butel JS (2009) Polyomavirus shedding in the stool of healthy adults. J Clin Microbiol 47:2388–2391
Kamminga S, van der Meijden E, Feltkamp MCW, Zaaijer HL (2018) Seroprevalence of fourteen human polyomaviruses determined in blood donors. PLoS ONE 13:e0206273
Bergallo M, Terlizzi ME, Astegiano S, Ciotti M, Babakir-Mina M, Perno CF, Cavallo R, Costa C (2009) Real time PCR TaqMan assays for detection of polyomaviruses KIV and WUV in clinical samples. J Virol Methods 162:69–74
Neske F, Blessing K, Pröttel A, Ullrich F, Kreth HW, Weissbrich B (2009) Detection of WU polyomavirus DNA by real-time PCR in nasopharyngeal aspirates, serum, and stool samples. J Clin Virol 44:115–118
Motamedi N, Mairhofer H, Nitschko H, Jäger G, Koszinowski UH (2012) The polyomaviruses WUPyV and KIPyV: a retrospective quantitative analysis in patients undergoing hematopoietic stem cell transplantation. Virol J 9:209
Kantola K, Sadeghi M, Lahtinen A, Koskenvuo M, Aaltonen LM, Möttönen M et al (2009) Merkel cell polyomavirus DNA in tumor-free tonsillar tissues and upper respiratory tract samples: implications for respiratory transmission and latency. J Clin Virol 45:292–295
Wen TK, Tang MB, Yu CP (2016) Detection of WU polyomavirus and norovirus genogroup II in stools of children with acute gastroenteritis in Taiwan. Southeast Asian J Trop Med Public Health 47:23–30
Wong AS, Cheng VC, Yuen KY, Kwong YL, Leung AY (2009) High frequency of polyoma BK virus shedding in the gastrointestinal tract after hematopoietic stem cell transplantation: a prospective and quantitative analysis. Bone Marrow Transpl 43:43–47
White MK, Safak M, Khalili K (2009) Regulation of gene expression in primate polyomaviruses. J Virol 83:10846–10856
Helle F, Brochot E, Handala L, Martin E, Castelain S, Francois C et al (2017) Biology of the BKPyV: an update. Viruses 9:E327
Reid CE, Li H, Sur G, Carmillo P, Bushnell S, Tizard R et al (2011) Sequencing and analysis of JC virus DNA from natalizumab-treated PML patients. J Infect Dis 204:237–244
Seppälä H, Virtanen E, Saarela M, Laine P, Paulín L, Mannonen L et al (2017) Single-molecule sequencing revealing the presence of distinct JC polyomavirus populations in patients with progressive multifocal leukoencephalopathy. J Infect Dis 215:889–895
Agostini HT, Yanagihara R, Davis V, Ryschkewitsch CF, Stoner GL (1997) Asian genotypes of JC virus in Native Americans and in a Pacific Island population: markers of viral evolution and human migration. Proc Natl Acad Sci USA 94:14542–14546
Agostini HT, Ryschkewitsch CF, Brubaker GR, Shao J, Stoner GL (1997) Five complete genomes of JC virus type 3 from Africans and African Americans. Arch Virol 142:637–655
Agostini HT, Deckhut A, Jobes DV, Girones R, Schlunck G, Prost MG et al (2001) Genotypes of JC virus in East, Central and Southwest Europe. J Gen Virol 82:1221–1331
Ferenczy MW, Marshall LJ, Nelson CD, Atwood WJ, Nath A, Khalili K, Major EO (2012) Molecular biology, epidemiology, and pathogenesis of progressive multifocal leukoencephalopathy, the JC virus-induced demyelinating disease of the human brain. Clin Microbiol Rev 25:471–506
Anzivino E, Rodio DM, Mischitelli M, Bellizzi A, Sciarra A, Salciccia S (2015) High frequency of JCV DNA detection in prostate cancer tissues. Cancer Genomics Proteomics 12:189–200
Funding
Anna Teresa Palamara: Funding acquisition (MIUR-PRIN Grant 2015W729WH); Valeria Pietropaolo: Funding acquisition (MIUR Research Grant RP11715C210E7374).
Author information
Authors and Affiliations
Contributions
CP: Conceptualization, methodology, data curation, writing the original draft, review and editing; FO, DA and DMR: Methodology; AM, AN: Statistical analysis; LC, JG, CF, WA: Patients’ enrollment and clinical data collection; VP, MC, ATP: Conceptualization, data curation, writing the original draft, review and editing. VP: Supervision.
Corresponding author
Ethics declarations
Conflict of interest
We certify that we have no affiliations with or involvement in any organization or entity with any financial interest, or non-financial interest regarding the subject matter or materials discussed in this manuscript.
Ethical approval
Ethical approval for the study was granted by the Ethical Committee of the University Hospital Tor Vergata (Rome, Italy). In case of minors, signed informed consent was obtained from the parents/LAT.
Additional information
Edited by Matthias J. Reddehase.
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Prezioso, C., Ciotti, M., Obregon, F. et al. Polyomaviruses shedding in stool of patients with hematological disorders: detection analysis and study of the non-coding control region’s genetic variability. Med Microbiol Immunol 208, 845–854 (2019). https://doi.org/10.1007/s00430-019-00630-9
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00430-019-00630-9