Genetic Variability of the Noncoding Control Region of Cutaneous Merkel Cell Polyomavirus: Identification of Geographically Related Genotypes

Merkel Cell Polyomavirus in the Context of Oral Squamous Cell Carcinoma and Oral Potentially Malignant Disorders

Despite recent advances in prevention, detection and treatment, oral squamous cell carcinoma (OSCC) remains a global health concern, strongly associated with environmental and lifestyle risk factors and infection with oncogenic viruses. Merkel Cell Polyomavirus (MCPyV), well known to be the causative agent of Merkel Cell Carcinoma (MCC) has been found in OSCC, suggesting its potential role as a co-factor in the development of oral cavity cancers. To improve our understanding about MCPyV in oral cavities, the detection and analysis of MCPyV DNA, transcripts and miRNA were performed on OSCCs and oral potentially malignant disorders (OPMDs). In addition, the cellular miR-375, known to be deregulated in tumors, was examined. MCPyV DNA was found in 3 out of 11 OSCC and 4 out of 12 OPMD samples, with a viral mean value of 1.49 × 102 copies/mL. Viral integration was not observed and LTAg and VP1 transcripts were detected. Viral miRNAs were not detected whereas the cellular miR-375 was found over expressed in all MCPyV positive oral specimens. Our results reported evidence of MCPyV replication in both OSCC and OPMD suggesting the oral cavity as a site of replicative MCPyV infection, therefore underscoring an active role of this virus in the occurrence of oral lesions.

Detection of Merkel cell polyomavirus in multiple primary oral squamous cell carcinomas

Oral microbiome studies have mainly focussed on bacteria, with the relationship between viruses and oral cancers remaining poorly understood. Oral cancers can develop even in the absence of any history of daily smoking or drinking. Oral cancer patients frequently have multiple primary cancers in the oral cavity and other organs, such as the upper gastrointestinal tract. Merkel cell polyomavirus (MCPyV) is a novel oncovirus identified from a subtype of skin cancer in 2008. In this study, we investigated the potential involvement of MCPyV in the pathogenesis of oral squamous cell carcinoma (OSCC). Participants comprised 115 Japanese patients with OSCC (single primary: 109 tumours in 109 patients; multiple primaries: 16 tumours in 6 patients) treated in our department between 2014 and 2017. DNA was extracted from formalin-fixed paraffin-embedded specimens of primary lesions. MCPyV DNA copy counts were analysed by quantitative real-time polymerase chain reaction. Twenty-four of the 115 patients (20.9%) were positive for MCPyV DNA. No association was found between presence or absence of MCPyV DNA and clinical characteristics other than number of primary lesions. The MCPyV DNA-positive rate was significantly higher for multiple primary OSCCs (62.5%, 10/16 tumours) than for single primary OSCCs (16.5%, 18/109 tumours; P < 0.001). Furthermore, MCPyV DNA load was significantly higher for patients with multiple primaries (P < 0.05). MCPyV was observed more frequently and DNA load was significantly higher with multiple primary OSCCs than with single primary OSCC. MCPyV may play some role as an oncovirus for multiple primary OSCCs.

Cutavirus on the skin in an Asian cohort: identification of a novel geographically related genotype

BACKGROUND

Cutavirus (CuV) is the newest human parvovirus and is currently receiving increasing attention because of its possible association with cutaneous T-cell lymphoma. Despite the pathogenetic potential of CuV, it has been detected in normal skin; however, little is known about the prevalence, infection levels, and genetic variations of this virus in the skin of the general population.

METHODS

We investigated the CuV DNA prevalence and viral loads concerning age, sampling location, and gender using 678 skin swabs collected from the normal-appearing skins of 339 Japanese participants aged 2-99 years. Phylogenetic analyses were also conducted based on the near-full-length CuV sequences identified in this study.

RESULTS

Both the CuV DNA prevalence and viral loads were significantly higher in the skin of elderly persons aged ≥60 years compared with those of persons aged < 60 years. CuV DNA tended to persist in the skin of elderly individuals. No significant difference in viral loads was observed between the skin of the upper arm and the skin of the forehead in CuV DNA-positive specimens. Significantly higher viral loads were evident in men vs. women, although no gender-associated differences in viral prevalence were noted. Phylogenetic analyses demonstrated the existence of Japanese-specific viruses that were genetically distinct from viruses prevalent in other areas, especially Europe.

CONCLUSIONS

This large study suggests that high levels of CuV DNA are prevalent on the skin of elderly adults. Our findings also indicated the prevalence of geographically related CuV genotypes. A follow-up study of this cohort should provide helpful information on whether CuV may become pathogenic.

Detection Analysis and Study of Genomic Region Variability of JCPyV, BKPyV, MCPyV, HPyV6, HPyV7 and QPyV in the Urine and Plasma of HIV-1-Infected Patients

Since it was clearly established that HIV/AIDS predisposes to the infection, persistence or reactivation of latent viruses, the prevalence of human polyomaviruses (HPyVs) among HIV-1-infected patients and a possible correlation between HPyVs and HIV sero-status were investigated. PCR was performed to detect and quantify JCPyV, BKPyV, MCPyV, HPyV6, HPyV7 and QPyV DNA in the urine and plasma samples of 103 HIV-1-infected patients. Subsequently, NCCR, VP1 and MCPyV LT sequences were examined. In addition, for MCPyV, the expression of transcripts for the LT gene was investigated. JCPyV, BKPyV and MCPyV's presence was reported, whereas HPyV6, HPyV7 and QPyV were not detected in any sample. Co-infection patterns of JCPyV, BKPyV and MCPyV were found. Archetype-like NCCRs were observed with some point mutations in plasma samples positive for JCPyV and BKPyV. The VP1 region was found to be highly conserved among these subjects. LT did not show mutations causing stop codons, and LT transcripts were expressed in MCPyV positive samples. A significant correlation between HPyVs' detection and a low level of CD4+ was reported. In conclusion, HPyV6, HPyV7 and QPyV seem to not have a clinical relevance in HIV-1 patients, whereas further studies are warranted to define the clinical importance of JCPyV, BKPyV and MCPyV DNA detection in these subjects.

Seroprevalence of Four Polyomaviruses Linked to Dermatological Diseases: New Findings and a Comprehensive Analysis

Our aim was to study the seroprevalence of human polyomaviruses (HPyV) linked to skin diseases. A total of 552 serum samples were analysed by the enzyme-linked immunosorbent assay to detect IgG antibodies against Merkel cell polyomavirus (MCPyV), HPyV6, HPyV7 and Trichodysplasia spinulosa-associated polyomavirus (TSPyV) using recombinant major capsid proteins of these viruses. The individuals (age 0.8−85 years, median 33) were sorted into seven age groups: <6, 6−10, 10−14, 14−21, 21−40, 40−60 and >60 years. The adulthood seroprevalence was 69.3%, 87.7%, 83.8% and 85% for MCPyV, HPyV6, HPyV7 and TSPyV, respectively. For all four polyomaviruses, there was increasing seropositivity with age until reaching the adulthood level. There was a significant increase in seroreactivity for those age groups in which the rate of already-infected individuals also showed significant differences. The adulthood seropositvity was relatively stable with ageing, except for TSPyV, for which elevated seropositivity was observed for the elderly (>60 years) age group. Since seroepidemiological data have been published with wide ranges for all the viruses studied, we performed a comprehensive analysis comparing the published age-specific seropositivities to our data. Although the cohorts, methods and even the antigens were variable among the studies, there were similar results for all studied polyomaviruses. For MCPyV, geographically distinct genotypes might exist, which might also result in the differences in the seroprevalence data. Additional studies with comparable study groups and methods are required to clarify whether there are geographical differences.

High prevalence of Merkel cell polyomavirus is associated with dysregulation in transcript levels of TLR9 and type I IFNs in a large cohort of CF patients from the Italian (Lazio) reference center for cystic fibrosis

Merkel cell polyomavirus (MCPyV) has been detected in respiratory specimens including those from Cystic Fibrosis (CF) patients, raising questions about its immunological and clinical relevance in the respiratory tract. MCPyV might promote an inappropriate antiviral response contributing to a chronic inflammatory response and resulting in detrimental effects in CF. Respiratory samples (n = 1138) were randomly collected from respiratory tract of CF patients (n = 539) during July 2018-October 2019. MCPyV-DNA detection was performed by Real Time-PCR and positive samples were characterized by sequencing of the NCCR genomic region. The transcript levels of Toll-like receptor 9 (TLR9) and type I interferon (IFN-I) genes (IFNα, IFNβ and IFNε) were examined by RT/Real Time-PCR assays. MCPyV-DNA was detected in 268 out of 1138 respiratory specimens (23.5%) without any difference in the prevalence of MCPyV-DNA according to age, gender or bacteriological status of CF individuals. Thirteen out of 137 CF patients remained positive for MCPyV-DNA over the time (a median follow-up period of 8.8 months). Detection of MCPyV-DNA in respiratory specimens was not associated with the occurrence of exacerbation events. Both MCPyV positive adolescents (11-24 years) and adults (>25 years) had lower mRNA levels of TLR9, IFNβ, IFNε and IFNα than the negative patients of the same age group, while MCPyV positive children produced increased levels of TLR9 and IFN-I genes (p < 0.05 for TLR9, IFNβ, IFNε) with respect to the negative ones. There were significant differences in TLR9 levels (p < 0.01), but not in those of IFNs, between MCPyV-DNA positive and negative patients with S. aureus, P. aeruginosa or both. Overall, these results indicate that MCPyV-DNA is frequently detected in the respiratory samples of CF patients and might influence the expression levels of IFN-related genes in an age dependent manner. The concomitant detection of MCPyV together with S. aureus and/or P. aeruginosa correlated with alterations in TLR9 levels suggesting that virus-bacteria coinfections might contribute to affect antiviral immunity in CF patients.

Evaluation of Merkel Cell Polyomavirus DNA in Tissue Samples from Italian Patients with Diagnosis of MCC

Because the incidence of Merkel cell carcinoma (MCC) has increased significantly during the last 10 years and it is recognized that Merkel cell polyomavirus (MCPyV) and ultraviolet (UV) radiation represent two different etiological inputs sharing clinical, histopathological, and prognostic similar features, although with different prognosis, this study investigated the detection of MCPyV in skin and lymph nodes with histological diagnosis of MCC. Formalin-fixed paraffin-embedded tissue (FFPE) were retrieved from archived specimens and MCPyV non-coding control region (NCCR) and viral capsid protein 1 (VP1) sequences were amplified and sequenced. Results provide an interesting observation concerning the discrepancy between the MCPyV DNA status in primary and metastatic sites: in fact, in all cases in which primary and metastatic lesions were investigated, MCPyV DNA was detected only in the primary lesions. Our data further support the “hit-and-run” theory, also proposed by other authors, and may lead to speculation that in some MCCs the virus is only necessary for the process of tumor initiation and that further mutations may render the tumor independent from the virus. Few point mutations were detected in the NCCR and only silent mutations were observed in the VP1 sequence compared to the MCPyV MCC350 isolate. To unequivocally establish a role of MCPyV in malignancies, additional well-controlled investigations are required, and larger cohorts should be examined.

Genetic Diversity of the Noncoding Control Region of the Novel Human Polyomaviruses

The genomes of polyomaviruses are characterized by their tripartite organization with an early region, a late region and a noncoding control region (NCCR). The early region encodes proteins involved in replication and transcription of the viral genome, while expression of the late region generates the capsid proteins. Transcription regulatory sequences for expression of the early and late genes, as well as the origin of replication are encompassed in the NCCR. Cell tropism of polyomaviruses not only depends on the appropriate receptors on the host cell, but cell-specific expression of the viral genes is also governed by the NCCR. Thus far, 15 polyomaviruses have been isolated from humans, though it remains to be established whether all of them are genuine human polyomaviruses (HPyVs). The sequences of the NCCR of these HPyVs show high genetic variability and have been best studied in the human polyomaviruses BK and JC. Rearranged NCCRs in BKPyV and JCPyV, the first HPyVs to be discovered approximately 30 years ago, have been associated with the pathogenic properties of these viruses in nephropathy and progressive multifocal leukoencephalopathy, respectively. Since 2007, thirteen novel PyVs have been isolated from humans: KIPyV, WUPyV, MCPyV, HPyV6, HPyV7, TSPyV, HPyV9, HPyV10, STLPyV, HPyV12, NJPyV, LIPyV and QPyV. This review describes all NCCR variants of the new HPyVs that have been reported in the literature and discusses the possible consequences of NCCR diversity in terms of promoter strength, putative transcription factor binding sites and possible association with diseases.

Viral Genomic Characterization and Replication Pattern of Human Polyomaviruses in Kidney Transplant Recipients

Human Polyomavirus (HPyV) infections are common, ranging from 60% to 100%. In kidney transplant (KTx) recipients, HPyVs have been associated with allograft nephropathy, progressive multifocal leukoencephalopathy, and skin cancer. Whether such complications are caused by viral reactivation or primary infection transmitted by the donor remains debated. This study aimed to investigate the replication pattern and genomic characterization of BK Polyomavirus (BKPyV), JC Polyomavirus (JCPyV), and Merkel Cell Polyomavirus (MCPyV) infections in KTx. Urine samples from 57 KTx donor/recipient pairs were collected immediately before organ retrieval/transplant and periodically up to post-operative day 540. Specimens were tested for the presence of BKPyV, JCPyV, and MCPyV genome by virus-specific Real-Time PCR and molecularly characterized. HPyVs genome was detected in 49.1% of donors and 77.2% of recipients. Sequences analysis revealed the archetypal strain for JCPyV, TU and Dunlop strains for BKPyV, and IIa-2 strain for MCPyV. VP1 genotyping showed a high frequency for JCPyV genotype 1 and BKPyV genotype I. Our experience demonstrates that after KTx, HPyVs genome remains stable over time with no emergence of quasi-species. HPyVs strains isolated in donor/recipient pairs are mostly identical, suggesting that viruses detected in the recipient may be transmitted by the allograft.

Human Polyomavirus 6 Detected in Cases of Eosinophilic Pustular Folliculitis

BACKGROUND

Human polyomaviruses (HPyVs) have been associated with several cutaneous inflammatory conditions. More investigation is needed to identify further presentations of cutaneous pathology associated with HPyVs. Our aim was to investigate the possible association of skin-tropic HPyVs with folliculitis, particularly eosinophilic pustular folliculitis (EPF).

METHODS

This study included 55 Japanese patients, comprising 13 patients with EPF and 42 patients with suppurative folliculitis. HPyV DNAs were detected by quantitative polymerase chain reaction. Expression of viral antigen and geographically related viral genotypes were also assessed.

RESULTS

Human polyomavirus 6 (HPyV6) DNA was found in 9 of 13 (69%) patients with EPF, a rate significantly higher than that found in suppurative folliculitis (1/42; 2%). Of the 7 HPyV6 DNA-positive EPF specimens analyzed, 4 were positive for HPyV6 small tumor antigen. All the HPyV6 strains detected in this study were of the Asian/Japanese genotype.

CONCLUSIONS

The predominant detection of HPyV6 DNA and the expression of viral antigen suggest a possible association between HPyV6 infection and EPF in a subset of patients. Worldwide studies are warranted to determine whether Asian/Japanese genotype HPyV6 is associated preferentially with the incidence and pathogenesis of this eosinophil-related skin disease that has an ethnic predilection for the East Asian population.

Regulation of Polyomavirus Transcription by Viral and Cellular Factors

Polyomavirus infection is widespread in the human population. This family of viruses normally maintains latent infection within the host cell but can cause a range of human pathologies, especially in immunocompromised individuals. Among several known pathogenic human polyomaviruses, JC polyomavirus (JCPyV) has the potential to cause the demyelinating disease progressive multifocal leukoencephalopathy (PML); BK polyomavirus (BKPyV) can cause nephropathy in kidney transplant recipients, and Merkel cell polyomavirus (MCPyV) is associated with a highly aggressive form of skin cancer, Merkel cell carcinoma (MCC). While the mechanisms by which these viruses give rise to the relevant diseases are not well understood, it is clear that the control of gene expression in each polyomavirus plays an important role in determining the infectious tropism of the virus as well as their potential to promote disease progression. In this review, we discuss the mechanisms governing the transcriptional regulation of these pathogenic human polyomaviruses in addition to the best-studied simian vacuolating virus 40 (SV40). We highlight the roles of viral cis-acting DNA elements, encoded proteins and miRNAs that control the viral gene expression. We will also underline the cellular transcription factors and epigenetic modifications that regulate the gene expression of these viruses.

Merkel Cell Polyomavirus (MCPyV) in the Context of Immunosuppression: Genetic Analysis of Noncoding Control Region (NCCR) Variability among a HIV-1-Positive Population

Background: Since limited data are available about the prevalence of Merkel cell polyomavirus (MCPyV) and the genetic variability of its noncoding control region (NCCR) in the context of immunosuppression, this study aimed to investigate the distribution of MCPyV in anatomical sites other than the skin and the behavior of NCCR among an HIV-1-positive population. Methods: Urine, plasma, and rectal swabs specimens from a cohort of 66 HIV-1-positive patients were collected and subjected to quantitative real-time polymerase chain reaction (qPCR) for MCPyV DNA detection. MCPyV-positive samples were amplified by nested PCR targeting the NCCR, and NCCRs alignment was carried out to evaluate the occurrence of mutations and to identify putative binding sites for cellular factors. Results: MCPyV DNA was detected in 10/66 urine, in 7/66 plasma, and in 23/66 rectal samples, with a median value of 5 × 10² copies/mL, 1.5 × 10² copies/mL, and 2.3 × 10³ copies/mL, respectively. NCCR sequence analysis revealed a high degree of homology with the MCC350 reference strain in urine, whereas transitions, transversions, and single or double deletions were observed in plasma and rectal swabs. In these latter samples, representative GTT and GTTGA insertions were also observed. Search for putative binding sites of cellular transcription factors showed that in several strains, deletions, insertions, or single base substitutions altered the NCCR canonical configuration. Conclusions: Sequencing analysis revealed the presence of numerous mutations in the NCCR, including insertions and deletions. Whether these mutations may have an impact on the pathogenic features of the virus remains to be determined. qPCR measured on average a low viral load in the specimens analyzed, with the exception of those with the GTTGA insertion.

Promoter activity of Merkel cell Polyomavirus variants in human dermal fibroblasts and a Merkel cell carcinoma cell line

Background

Merkel cell polyomavirus (MCPyV) is a human polyomavirus that establishes a life-long harmless infection in most individuals, with dermal fibroblasts believed to be the natural host cell. However, this virus is the major cause of Merkel cell carcinoma (MCC), an aggressive skin cancer. Several MCPyV variants with polymorphism in their promoter region have been isolated, but it is not known whether these differences affect the biological properties of the virus.

Methods

Using transient transfection studies in human dermal fibroblasts and the MCC cell line MCC13, we compared the transcription activity of the early and late promoters of the most commonly described non-coding control region MCPyV variant and six other isolates containing specific mutation patterns.

Results

Both the early and late promoters were significantly stronger in human dermal fibroblasts compared with MCC13 cells, and a different promoter strength between the MCPyV variants was observed. The expression of full-length large T-antigen, a viral protein that regulates early and late promoter activity, inhibited early and late promoter activities in both cell lines. Nonetheless, a truncated large T-antigen, which is expressed in virus-positive MCCs, stimulated the activity of its cognate promoter.

Conclusion

The promoter activities of all MCPyV variants tested was stronger in human dermal fibroblasts, a cell line that supports viral replication, than in MCC13 cells, which are not permissive for MCPyV. Truncated large T-antigen, but not full-length large T-antigen stimulated viral promoter activity. Whether, the difference in promoter strength and regulation by large T-antigen may affect the replication and tumorigenic properties of the virus remains to be determined.

Intra-patient viral evolution in polyomavirus-related diseases

Human polyomaviruses show relatively little genetic polymorphism between isolates, indicating that these viruses are genetically stable between hosts. However, it has become increasingly clear that intra-host molecular evolution is a feature of some polyomavirus (PyV) infections in humans. Mutations inducing premature stop codons in the early region of the integrated Merkel cell PyV genome lead to the expression of a truncated form of the large tumour (LT) antigen that is critical for the transformation of Merkel cell carcinoma (MCC) cells. Non-coding control region (NCCR) rearrangements and point mutations in virion protein (VP) 1 have been described in both JCPyV and BKPyV infections. In the context of JCPyV infection, molecular evolution at both these loci allows the virus to replicate effectively in the central nervous system, thereby leading to the development of progressive multifocal leukoencephalopathy (PML). In BKPyV infection, NCCR rearrangements have been linked to higher rates of virus replication in the kidney, and are proposed to play a direct causal role in the development of PyV-associated nephropathy. In all three of these infections, therefore, intra-host viral evolution appears to be an essential component of the disease process. This article is part of the theme issue 'Silent cancer agents: multi-disciplinary modelling of human DNA oncoviruses'.

A model for the formation of the duplicated enhancers found in polyomavirus regulatory regions

When purified from persistent infections, the genomes of most human polyomaviruses contain single enhancers. However, when isolated from productively infected cells from immunocompromised individuals, the genomes of several polyomaviruses contain duplicated enhancers that promote a number of polyoma-based diseases. The mechanism(s) that gives rise to the duplicated enhancers in the polyomaviruses is, however, not known. Herein we propose a model for the duplication of the enhancers that is based on recent advances in our understanding of; 1) the initiation of polyomavirus DNA replication, 2) the formation of long flaps via displacement synthesis and 3) the subsequent generation of duplicated enhancers via double stranded break repair. Finally, we discuss the possibility that the polyomavirus based replication dependent enhancer duplication model may be relevant to the enhancer-associated rearrangements detected in human genomes that are associated with various diseases, including cancers.

Prevalence and Viral Loads of Cutaneous Human Polyomaviruses in the Skin of Patients With Chronic Inflammatory Skin Diseases

BACKGROUND

Human skin microorganisms have been associated with various skin diseases. However, most studies have focused on bacterial communities, and little is known about normally resident skin viruses such as the Polyomaviridae and their association with cutaneous disorders.

METHODS

We investigated the infection levels of Merkel cell polyomavirus (MCPyV), human polyomavirus 6 (HPyV6), and human polyomavirus 7 (HPyV7), using triplet skin swabs collected from lesional and nonlesional skins of 86 Japanese patients with inflammatory skin diseases and mycosis fungoides, and from 149 healthy control individuals.

RESULTS

This age-matched case-control study provides the first analyses of the loads of polyomaviruses in association with various skin diseases. The viral loads were significantly higher for HPyV6/HPyV7 and lower for MCPyV in patients with psoriasis. The viral load variation was observed not only at lesion sites, but also at clinically unaffected skin sites in most of the patients. The viral strains tested were all of the Asian/Japanese genotype.

CONCLUSIONS

Our findings suggest a covariation in the infection levels of cutaneous polyomaviruses in certain inflammatory skin conditions. Worldwide prospective longitudinal studies are warranted to understand the influence of such alterations on the pathogenesis of inflammatory skin disorders.

Human Polyomaviruses in the Cerebrospinal Fluid of Neurological Patients

BACKGROUND

Central nervous system (CNS) infections by human polyomaviruses (HPyVs), with the exception of JC (JCPyV), have been poorly studied.

METHODS

In total, 234 cerebrospinal fluid (CSF) samples were collected from patients affected with neurological disorders. DNA was isolated and subjected to quantitative real-time PCR (Q-PCR) for the detection of six HPyVs: JCPyV, BKPyV, Merkel cell PyV (MCPyV), HPyV6, HPyV7, and HPyV9. Where possible, the molecular characterization of the viral strains was carried out by nested PCR and automated sequencing.

RESULTS

JCPyV was detected in 3/234 (1.3%), BKPyV in 15/234 (6.4%), MCPyV in 22/234 (9.4%), and HPyV6 in 1/234 (0.4%) CSF samples. JCPyV was detected at the highest (p < 0.05) mean load (3.7 × 10⁷ copies/mL), followed by BKPyV (1.9 × 10⁶ copies/mL), MCPyV (1.9 × 10⁵ copies/mL), and HPyV6 (3.3 × 10⁴ copies/mL). The noncoding control regions (NCCRs) of the sequenced viral strains were rearranged.

CONCLUSIONS

HPyVs other than JCPyV were found in the CSF of patients affected with different neurological diseases, probably as bystanders, rather than etiological agents of the disease. However, the fact that they can be latent in the CNS should be considered, especially in immunosuppressed patients.

Polyomaviruses shedding in stool of patients with hematological disorders: detection analysis and study of the non-coding control region's genetic variability

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.