Skip to main content

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



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.


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.


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.


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.


Cutavirus (CuV), a member of the Protoparvovirus genus, is a small nonenveloped virus with a single-stranded linear DNA genome of approximately 4.5 kb [1,2,3]. This virus is the newest human parvoviruses, discovered initially in fecal specimens from Brazilian children with diarrhea using metagenomics [1]. The CuV genome was subsequently identified in cutaneous T-cell lymphoma (CTCL) using in silico analysis in preexisting metagenomic libraries; moreover, it was detected in 4/17 (24%) specimens of CTCL from French patients by polymerase chain reaction (PCR) [1]. CuV DNA was also found in 1/10 (10%) biopsy specimens from Danish patients with cutaneous melanoma [4]. Subsequent studies reported the detection of CuV DNA in 6/117 (5%) and 4/25 (16%) biopsy specimens from patients with CTCL from Germany [5] and Finland [6], respectively. Thus, studies on CuV infection have mainly focused on malignant cutaneous tumors, with a possible association between CuV and CTCL.

To date, only a few studies have assessed the prevalence of CuV in the skin of healthy individuals. No viral DNA was found in the skin biopsy specimens from three French and 98 Finnish healthy adults [1, 6]. Using skin surface swabs, a study reported the detection rate of CuV DNA (3.8%, 9/237) in the skin of immunocompetent adults in Germany [7]. Thus, despite the skin tropism and the pathogenetic potential of CuV, the prevalence and viral loads of CuV in the skin among the general population have not been well delineated. Furthermore, epidemiological studies of CuV in the skin have all stemmed from individuals living in Europe. No data are available regarding the prevalence of CuV in the skin in other areas. Therefore, a worldwide survey including an Asian population will provide useful information on the geographical distribution of this novel human virus.

Previous studies have demonstrated that some skin-tropic viruses have distinct genotypes linked to the geographic origin of the infected individuals. Merkel cell polyomavirus (MCPyV) is etiologically linked to cutaneous Merkel cell carcinoma (MCC) [8]; moreover, it exists in healthy skin as a skin commensal [9, 10]. Several studies have shown that MCPyV has two major geographically related genotypes, namely European/North American and Asian/Japanese genotypes [11,12,13,14]. In turn, human polyomavirus 6 (HPyV6) and human polyomavirus 7 (HPyV7) are cutaneous viruses that are comprised of the human skin virome [15, 16]. We have also shown that the viral sequences of HPyV6 and HPyV7 vary between the prevalence in the skin of white and Asian populations [17]. Because it has been suggested that CuV may be part of the human skin virome [2], these findings highlight the need for studies to determine whether the CuV strains present in the skin also exhibit genotypic variation between white and Asian populations.

Based on this background, we were the first to investigate the age-specific prevalence and viral loads of CuV in 678 swab specimens collected from two different areas of the normal-appearing skin of individuals in a large Asian cohort. We also estimated the persistence and acquisition of CuV DNA in a subset of this cohort. Moreover, we conducted the first phylogenetic analyses of normal-skin-derived CuV strains using the near-full-length CuV genome, revealing the existence of Japanese-specific viruses genetically distinct from previously reported strains from Europe.

Materials and methods

Study participants and specimen collection

This study cohort included 339 Japanese individuals (persons of Japanese descent who resided in Japan) enrolled from 2014 to 2020. The participants ranged from 2 to 99 years (median age, 69 years); 213 were women, and 126 were men. In this study, gender refers to socially constructed roles, behaviors, and identities of women, men, and gender-diverse people, according to the Sex and Gender Equity in Research (SAGER) guidelines [18]. The authors interviewed the participants in this study for personal information, including current medical history of cancer and skin diseases and medication history. The participants and their families declared they were free of cancer, skin diseases, or having received immunosuppressive agents. The local skin areas from which the swabs were collected were carefully inspected and confirmed to appear healthy. However, the rest of the skin was not further examined. A total of 678 skin smear specimens were collected from about 50 cm2 of normal-appearing skin in two different regions in the same individuals, including the skin of the upper arm and the skin from the forehead to the top of the head. Sample collection was performed by rubbing the skin back and forth 5–10 times using sterile cotton swabs moistened with phosphate-buffered saline (PBS). All swabs were placed in a separate vial containing sterile PBS, and DNA was extracted within 2 h of sample collection. The subjects were divided according to age into six groups (Table 1).

Table 1 Age-specific prevalence of cutavirus DNA in swab specimens of normal skin of the arm and forehead

Quantification of the viral genome

DNA extraction was performed using the standard phenol: chloroform method. Concentrations of extracted DNAs were measured using BioPhotometer Plus (Eppendorf, Hamburg, Germany). Twenty-nanogram aliquots of the extracted DNA were analyzed for the detection and quantification of CuV using TaqMan-based real-time quantitative PCR (qPCR), as described elsewhere [19]. Primers and probes were prepared to amplify a region of the gene encoding the CuV viral capsid protein 2 (VP2). The forward and reverse primers are located at nucleotide positions 4245–4269 and 4316–4335 based on the GenBank CuV sequence [accession number NC_039050]), respectively (Additional file 1: Table S1) [6]. The reaction mixture was prepared in a total volume of 20 µl containing the TaqMan Gene Expression Master Mix (Thermo Fisher Scientific, Waltham, MA, USA), 900 nM of each primer, and 250 nM dual-labeled probe. Water was added to the PCR reaction mixture instead of extracted DNA as a negative PCR-amplification control in each PCR run. The negative control did not show any amplification in all of the specimens tested in this study. Precautions were also taken to prevent contamination in the PCR assays as described elsewhere [20], although negative DNA extraction controls were not included. The RNase P gene was amplified in separate PCR runs as a positive control to confirm the presence of PCR-amplifiable DNA. A PCR assay was performed using the same primers, and the PCR product (CuV VP2 sequence; nucleotide positions 4245–4335) was cloned into the pMD20-T vector (Takara Bio, Shiga, Japan). We prepared 10-fold serial dilutions using the cloned plasmid DNA to generate a standard curve, from which we calculated the viral copy number. The results are expressed as viral DNA copies/ng DNA. Because samples with cycle threshold values ≤ 41 were considered positive for targeted viral nucleic acids [13], we considered swab specimens with ≥ 3 × 10–2 copies/ng DNA to be positive.

Viral DNA sequencing analysis

The near-full-length CuV sequences (4455 bp, nucleotide positions 2–4456) were amplified by PCR using different combinations of 13 primer sets (Additional file 1: Table S1). The purified PCR products were sequenced directly. A total of 36 CuV sequences obtained in this study were deposited in the GenBank database under accession numbers LC744018–LC744023 and LC760811–LC760840.

Phylogenetic analysis and nucleotide identity analysis

The nucleotide sequences were aligned using Clustal W [21]. Phylogenetic trees were constructed using the maximum-likelihood method in MEGA X [22]. The bootstrap values were based on 1000 replicates for all trees. Nucleotide identity analysis was performed using the nucleotide BLAST program with default parameters [23].

Statistical analyses

Any correlations with CuV positivity rates were analyzed using Fisher’s exact test or Pearson’s chi-squared test if any values were < 5. The differences in viral loads were compared using the Mann–Whitney nonparametric U test. All statistical analyses were performed using R, version 1.61, with its graphical user interface, EZR [24]. Significance was set at P < 0.05.


CuV DNA prevalence and viral loads in normal skin

We tested 678 swab specimens obtained from the normal-appearing skin of the upper arm and forehead of 339 Japanese individuals, to detect and quantify the CuV genome. All specimens were tested twice independently in a blind manner. The results of the first and second tests were concordant for CuV DNA positivity or negativity in 657 of 678 (96.9%) specimens. The remaining samples with low CuV loads were all around the 3 × 10–2 copies/ng DNA borderlines and were determined to be positive or negative on the third test.

Overall, CuV DNA was detected in 53 of 339 (15.6%) arm skin specimens and 74 of 339 (21.8%) forehead skin specimens. The age-related prevalence was assessed. The detection rate according to the swab-sampling site in each age group is listed in Table 1. The positivity rates in groups aged 2–59 years were 7.9% (11/139) for the arm swabs and 7.2% (10/139) for the forehead swabs; the respective rate significantly increased to 21.0% (42/200; P = 0.002) and 32.0% (64/200; P < 0.001) in older groups (individuals aged 60–99 years). The rates of simultaneous detection of CuV DNA in both the arm and forehead skin specimens in the same individuals were also significantly higher in groups aged 60–99 years (13.0% [26/200]) than they were in younger groups (3.6% [5/139]; P = 0.003). No significant difference in site-specific detection rates between the arm and forehead was observed among the age groups (P = 0.354). The detection rate on either of the skins according to gender in each age group is listed in Table 2. Overall, no gender-associated differences in CuV DNA prevalence were noted (29.6% [63/213] for women and 26.2% [33/126] for men; P = 0.586).

Table 2 Gender-specific prevalence of cutavirus DNA in skin swab specimens according to age group

Next, we evaluated the CuV loads in skin swabs according to age group. Box plots of CuV DNA levels showed significantly higher viral loads in both the arm and forehead skins of persons aged ≥60 years compared with those of persons aged < 60 years (P < 0.001; Fig. 1). Thus, high loads of CuV DNA likely to exist in the skin of elderly individuals. For example, a CuV load of > 100 copies/ng DNA was detected in the arm and/or forehead skin in nine subjects; eight of whom (89%) were older than 80 years. Of note, a 91-year-old woman carried a very high CuV load of > 10,000 copies/ng DNA in her forehead skin without any apparent symptoms. We also conducted comparative analyses of the CuV loads between skin swabs obtained from the arms and foreheads and between those from women and men in CuV DNA-positive specimens (Fig. 2). No significant difference in site-specific viral loads was detected between the arm and forehead samples (P = 0.663). In turn, significantly higher viral loads were evident in men vs. women (P = 0.009).

Fig. 1
figure 1

Box plots showing age-specific cutavirus (CuV) DNA loads in skin swabs. (A) CuV DNA loads in skin swabs from the upper arm. (B) CuV DNA loads in skin swabs from the forehead. The horizontal bars extend the median by 1.5 times the interquartile ranges. The comparisons of the viral loads between age groups were performed using the Mann–Whitney nonparametric U test, and P values are presented above the box plots. The number of subjects is presented for each age group below the box plots

Fig. 2
figure 2

Box plots showing cutavirus (CuV) DNA loads in skin swabs positive for CuV DNA. (A) CuV DNA loads according to swab-sampling location. (B) CuV DNA loads according to gender. The horizontal bars extend the median by 1.5 times the interquartile ranges. The comparisons of the viral loads between age groups were performed using the Mann–Whitney nonparametric U test, and P values are presented above the box plots. The number of specimens is presented for each group below the box plots

Persistence and acquisition of CuV DNA in normal skin

We reassessed 46 subjects aged 32–98 years from whom both arm and forehead skin swabs were obtained approximately 6 months after collecting the first sample. The second specimens were collected from the same area as the first specimens, albeit with some minor variation in the exact site of collection. To investigate viral persistence, defined as the presence of CuV DNA at the two consecutive visits, we tested five and nine subjects whose first arm and forehead swabs contained viral DNA, respectively (Table 3). Of these individuals, three (60%) and four (44%), respectively, were found to have retained the CuV genome at a comparable level in the second skin swabs. Of note, all individuals with retained CuV DNA were older than 70. To investigate viral acquisition, defined as the new detection of CuV DNA that was not detected at the first visit, we tested 41 and 37 subjects whose first arm and forehead swabs were negative for CuV DNA, respectively (Table 3). The CuV viral genome was detected in the second skin swabs of five (12%) and 10 (27%) of these individuals, respectively. Similar to that observed for viral persistence, the acquisition of the CuV genome tended to occur more frequently in the skin of elderly individuals.

Table 3 Persistence and acquisition of cutavirus DNA among individuals in Japan

Phylogenetic analysis and geographic genetic diversity of CuV

Next, we analyzed the phylogenetic relationships between the CuV recovered from skin specimens in the Japanese population studied in this study and those reported previously in other populations. In this study, 36 near-full-length CuV sequences (4455 bp) were obtained successfully from the skin swabs of different individuals (nine sequences from the arm and 27 from the forehead), irrespective of the viral load of the specimens and the age/gender of the subjects (Table 4). As of December 1, 2022, four near-complete CuV sequences have been deposited in GenBank. These include sequences from one virus, termed BR-337, that was first identified in a fecal specimen from a patient in Brazil (GenBank accession no. KT868811) [1]; two viruses, termed FR-D and FR-F, that were recovered from skin biopsy specimens from patients with CTCL from France (KT868114 and KT868815, respectively) [1]; and one virus, termed CGC5-268, from a melanoma sample from a patient from Denmark (KX685945) [4]. We performed a phylogenetic study using these 40 CuV sequences (Fig. 3A). The phylogenetic tree generated using the maximum-likelihood method showed that all of the Japanese CuV sequences formed one major clade (designated clade 1) that included the Brazilian BR-337 sequence but was clearly distinct from a clade (clade 2) comprising an additional three sequences originating from European individuals. Slight sequence differences were seen in the phylogenetic tree of the 36 Japanese CuV sequences amplified from distinct individuals.

Table 4 Nucleotide identities between the near-full-length cutavirus sequences (4455 bp) identified in the current and previous studies
Fig. 3
figure 3

Phylogenetic trees generated using the maximum-likelihood method. (A) A phylogenetic tree was constructed based on 40 near-full-length cutavirus (CuV) sequences (4455 bp). These included 36 sequences that were successfully amplified from skin swab specimens from the Japanese individuals included in this study and four sequences retrieved from GenBank. The two major nucleotype clades (clades 1 and 2) are indicated. (B) A phylogenetic tree was constructed based on 55 CuV viral capsid protein 1 (VP1)/VP2 sequences with 538 bp (nucleotide positions 3778–4315 bp). These included the 36 sequences from the Japanese individuals and 19 sequences recovered from individuals of various geographic origins that were retrieved from GenBank as of December 1, 2022. The CuV reference sequence of the BR-337 strain (GenBank accession number no., NC_039050) is indicated with an asterisk. The sequences obtained from the Japanese individuals in the present study are presented in bold and colored in red. The CuV strains from Brazil, Botswana, and Europe are colored green, purple, and blue, respectively. The Brazilian and Botswanan strains were from fecal specimens, whereas the remaining strains were from skin specimens. The names of the sequences are indicated together with the GenBank accession numbers and country of origin in parentheses. The percentage bootstrap values calculated from 1000 replicates are indicated at the internal nodes. The scale bars represent the number of substitutions per site

We also conducted a phylogenetic study using 55 partial VP1/VP2 sequences (Fig. 3B). According to the data available on GenBank on December 1, 2022, a total of 19 partial sequences of CuV VP1/VP2 with 538 bp (nucleotide positions 3778–4315) have been reported from Brazil, Botswana, and Europe (including Finland, Denmark, and France), which were identified from skin or fecal specimens. As expected, all of the Japanese sequences belonged to clade 1. Importantly, the Japanese strains formed a paraphyletic group (tentatively designated Japan clade or clade 1a) distinct from other clades comprising the sequences originating from other regions (clades 1b and 2). As an exception, one strain from Brazil (BR-450, KT868812) was included in clade 1a. Thus, these phylogenetic analyses provided evidence of a specific Japanese geographic genotype of CuV. A phylogenetic tree constructed based on 44 complete CuV VP2 sequences (including the 36 Japanese sequences and eight sequences retrieved from GenBank), together with the VP2 sequences of gray fox amdovirus and bufavirus used as outgroups, also showed the presence of the Japan clade in the CuV group (Additional file 2: Fig. S1).

We compared the sequences of the Japanese strains with the CuV reference sequence of BR-337 in the genes encoding the nonstructural protein 1 (NS1) (nucleotide positions 1–1980), the middle open reading frame (ORF) (nucleotide positions 2022–2357), and the VP2 gene (nucleotide positions 2747–4456) which includes a part of the VP1 gene (Fig. 4). All of the Japanese strains had the following single-nucleotide substitutions: G→A at nucleotide position 84, G→A at position 621, G→A at position 696, T→C at position 1102, and T→C at position 1590 in NS1; C→A at position 2076 and T→A or C at position 2195 in middle ORF; and T→A at position 3406, T→C at position 3512, T→C at position 3586, T→C at position 4045, and T→A at position 4186 in VP2. Compared with the Japanese strains, the European FR-D, FR-F, and CGG5-268 strains, which belonged to clade 2, exhibited many different genetic alterations, including single-nucleotide substitutions with amino acid replacements and base deletions. These alterations appear to define the characteristic phylogenetic features of the Japanese strains.

Fig. 4
figure 4

Schematic diagrams showing the summaries of cutavirus (CuV) gene alterations. The sequences of the genes encoding the nonstructural protein 1 (NS1) and viral capsid protein 2 (VP2) and the sequence of the middle open reading frame (ORF) were compared with the CuV reference sequence of the BR-337 virus. The Japanese viruses included the 36 CuV sequences obtained in this study. The FR-D, FR-F, and CGG5-268 European CuV, which belong to clade 2, were also included. Filled ovals indicate alterations specific to the Japanese viruses. Open ovals indicate alterations common to the Japanese and European viruses. Filled rhombuses indicate alterations specific to the European viruses. Amino acid replacements are also shown. Nucleotide numbers refer to the sequences of the BR-337 strain (GenBank accession number no., NC_039050).

We also analyzed nucleotide identities of the near-full-length CuV sequences (4455 bp) between the 36 Japanese strains obtained in this study and four viruses, BR-337, FR-D, FR-F, and CGG5-268 (Table 4). We confirmed that our sequences of clade 1a were closer to the BR-337 sequence of clade 1b than the other sequences of clade 2. The Japanese CuV sequences were 97.40–98.00% and 93.77–94.71%, identical to the sequence of BR-337 and the sequences of FR-D, FR-F, and CGG5-268, respectively. The percent nucleotide identity between the Japanese CuV sequences was 97.53–99.98% (Additional file 3: Table S2).


CuV is receiving increasing attention because studies have suggested a causative role for this novel skin-tropic virus in CTCL [1, 2, 6, 25], although this warrants confirmation [5]. In this study, we examined the prevalence and viral load levels of CuV in the skin of 339 asymptomatic Japanese subjects who were free of skin diseases.

The present study had several strengths. One advantage was the availability of a large number of samples from participants. Because of the pathogenetic potential of CuV on the human skin, our survey was significant because it is thought that the presence of a virus at the site of disease development is essential for virus-associated disorders. Here, we provided the first age-specific analysis of the prevalence and quantification of CuV DNA in these skin samples. It should be noted that high burdens of CuV DNA were observed in persons aged ≥60 years. After 60 years, the viral prevalence and loads continued to increase, peaking in the groups aged ≥80 years. Very high viral loads of > 100 copies/ng DNA were predominantly observed in the skin specimens from individuals aged > 80 years. Thus, our assessment showed that the prevalence and viral loads of CuV increased with age. This may indicate a reduction of immunity with increasing age, which could lead to the activation of the virus in older individuals. In previous surveys, Väisänen et al. [6] did not detect CuV DNA among 159 skin biopsy specimens from 98 healthy immunocompetent adults (median age, 43 years; range, 18–67 years) from Finland, whereas they found it in four (2.9%) of 136 skin biopsy specimens from organ transplant recipients (median age, 62 years; range, 22–83 years), suggesting that CuV DNA is found more frequently in immunosuppressed patients. However, their finding related to healthy adults contrasts with our results. This might be attributed to the different approaches to sampling specimens because the viral detection rates obtained using surface skin swabs and dermis biopsies may vary [26]. Otherwise, because their survey of healthy adults was conducted on relatively younger individuals, CuV DNA could hardly be detected in their cohort. Wieland et al. [7] found CuV DNA more frequently in the skin swabs of human immunodeficiency virus-positive men (17.1%) than in those of healthy controls (3.8%) in Germany. Mohanraj et al. [27] analyzed fecal specimens in cohorts of patients with gastroenteritis in Finland. They showed that the CuV DNA prevalence among individuals older than 60 years (5.1%) was significantly higher than that among individuals younger than below 60 (0.2%). These findings support our observation that CuV DNA was more frequent on the skin of elderly individuals with waning immunity.

We also present the first data on the prevalence levels of the CuV genome in two different skin regions, including the skin of the forehead and the skin of the upper arm. There were no significant differences in the site-specific detection rates or viral loads, which was in contrast to the observation that the viral loads of MCPyV were higher on the forehead skin than on the arm skin [13]. In fact, MCPyV-associated MCC preferentially develops with a background of high viral load in sun-exposed skin areas such as forehead skin [10]. In contrast, mycosis fungoides, the most prevalent form of CTCL, often develop on the trunk, which is not commonly exposed to the sun [28]; however, other areas, such as the head and neck, are also affected [29, 30]. Risk factors for mycosis fungoides include older age and male sex [30, 31]. In this study, we evaluated the effect of gender differences on CuV DNA levels and found significantly higher viral loads in men. The CuV DNA prevalence patterns observed in this study should provide useful information to understand the CuV status on the skin.

We also evaluated the 6-month risk of CuV DNA persistence and acquisition. Our study showed that elderly subjects tended to retain CuV DNA in primarily infected regions. Similarly, the viral acquisition was more likely to occur in elderly individuals. These findings support the notion that the waning of the immune system with increasing age prevents the host from eradicating infecting viruses after the initial infection and/or preventing new infections. One limitation of the assay included the relatively small number of subjects for whom skin swabs were available for consecutive collection. Another limitation was that the short follow-up period of 6 months may not have allowed adequate assessment of viral persistence and may have underestimated viral acquisition. The classical human parvovirus B19 (B19V) can persist in the skin not only in symptomatic but also in healthy individuals [32,33,34,35]. Therefore, it is possible that CuV also persists after the initial infection in the skin of individuals, especially elderly persons.

Finally, this was the first report of the near-full-length CuV sequences identified in an Asian cohort. Our phylogenetic analyses demonstrated the existence of Japanese-specific CuV that were genetically distinct from the viruses recovered from fecal specimens of individuals from Brazil and Botswana and skin specimens of individuals from Europe. Based on the results of the phylogenetic analyses, we propose two major distinct genotypes, tentatively designated CuV genotypes 1 and 2 (Fig. 3). CuV genotype 1 can be divided into two subtypes; subtype 1a is predominantly composed of Japanese viruses that formed a paraphyletic clade, whereas subtype 1b includes a small number of viruses reported from other areas. CuV genotype 2 includes the majority of European viruses reported to date. Studies have shown that B19V is classified into three distinct genotypes (1, 2, and 3) and that the frequencies of the B19V genotypes vary in different continents [36,37,38]. Our findings suggest that different CuV genotypes are also prevalent worldwide. However, the global distribution patterns of CuV on normal skin have not been clarified; the information on CuV gene sequences available to date stems mainly from Europe, whereas that of North and South America, Africa, Oceania, and Asia is very limited or not available. Therefore, the findings should stimulate worldwide studies on geographically related variant CuV genotypes. Moreover, an important issue is whether a specific genotype is preferentially related to CuV-associated disorders. CTCL appears to be more frequent in Asian populations than in European populations [39]. In this context, whether the Japanese-specific CuV genotype represents a potential trigger of certain types of CTCL or merely indicates the presence of a Japan-related topotype warrant further investigation.


We have shown in this large study that both the CuV DNA prevalence rates and viral loads in the skin were increased among elderly individuals. Hence, a prospective follow-up study of this cohort should provide useful information on whether CuV might become pathogenic. We also detected Japanese-specific CuV strains. Our findings suggest that the viral genotype may vary according to ethnicity; therefore, further analyses of specimens from different countries are warranted to strengthen our original CuV genotype classification. Finally, our findings are expected to promote studies on whether specific geographically related CuV genotypes might predispose the development of CuV infection into CTCL.

Data Availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.


  1. Phan TG, Dreno B, da Costa AC, Li L, Orlandi P, Deng X, Kapusinszky B, Siqueira J, Knol AC, Halary F, Dantal J, Alexander KA, Pesavento PA, Delwart E. A new protoparvovirus in human fecal samples and cutaneous T cell lymphomas (mycosis fungoides). Virology. 2016;496:299–305.

    Article  CAS  PubMed  Google Scholar 

  2. Phan T, Nagaro K. Cutavirus: a newly discovered parvovirus on the rise. Infect Genet Evol. 2020;80:104175.

    Article  CAS  PubMed  Google Scholar 

  3. Väisänen E, Fu Y, Hedman K, Söderlund-Venermo M. Human protoparvoviruses. Viruses. 2017;9:354.

  4. Mollerup S, Fridholm H, Vinner L, Kjartansdóttir KR, Friis-Nielsen J, Asplund M, Herrera JA, Steiniche T, Mourier T, Brunak S, Willerslev E, Izarzugaza JM, Hansen AJ, Nielsen LP. Cutavirus in cutaneous malignant melanoma. Emerg Infect Dis. 2017;23:363–5.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Kreuter A, Nasserani N, Tigges C, Oellig F, Silling S, Akgül B, Wieland U. Cutavirus infection in primary cutaneous B- and T-cell lymphoma. JAMA Dermatol. 2018;154:965–7.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Väisänen E, Fu Y, Koskenmies S, Fyhrquist N, Wang Y, Keinonen A, Mäkisalo H, Väkevä L, Pitkänen S, Ranki A, Hedman K, Söderlund-Venermo M. Cutavirus DNA in malignant and nonmalignant skin of cutaneous T-cell lymphoma and organ transplant patients but not of healthy adults. Clin Infect Dis. 2019;68:1904–10.

    Article  PubMed  Google Scholar 

  7. Wieland U, Silling S, Hufbauer M, Mauch C, Zigrino P, Oellig F, Kreuter A, Akgül B. No evidence for role of cutavirus in malignant melanoma. Emerg Infect Dis. 2019;25:1600–2.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Feng H, Shuda M, Chang Y, Moore PS. Clonal integration of a polyomavirus in human Merkel cell carcinoma. Science. 2008;319:1096–100.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Foulongne V, Sauvage V, Hebert C, Dereure O, Cheval J, Gouilh MA, Pariente K, Segondy M, Burguière A, Manuguerra JC, Caro V, Eloit M. Human skin microbiota: high diversity of DNA viruses identified on the human skin by high throughput sequencing. PLoS ONE. 2012;7:e38499.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Hashida Y, Nakajima K, Nakajima H, Shiga T, Tanaka M, Murakami M, Matsuzaki S, Naganuma S, Kuroda N, Seki Y, Katano H, Sano S, Daibata M. High load of Merkel cell polyomavirus DNA detected in the normal skin of Japanese patients with Merkel cell carcinoma. J Clin Virol. 2016;82:101–7.

    Article  CAS  PubMed  Google Scholar 

  11. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Martel-Jantin C, Filippone C, Tortevoye P, Afonso PV, Betsem E, Descorps-Declere S, Nicol JT, Touzé A, Coursaget P, Crouzat M, Berthet N, Cassar O, Gessain A. Molecular epidemiology of Merkel cell polyomavirus: evidence for geographically related variant genotypes. J Clin Microbiol. 2014;52:1687–90.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Hashida Y, Kamioka M, Tanaka M, Hosokawa S, Murakami M, Nakajima K, Kikuchi H, Fujieda M, Sano S, Daibata M. Ecology of Merkel cell polyomavirus in healthy skin among individuals in an Asian cohort. J Infect Dis. 2016;213:1708–16.

    Article  PubMed  Google Scholar 

  14. Hashida Y, Higuchi T, Matsui K, Shibata Y, Nakajima K, Sano S, Daibata M. Genetic variability of the noncoding control region of cutaneous Merkel cell polyomavirus: identification of geographically related genotypes. J Infect Dis. 2018;217:1601–11.

    Article  CAS  PubMed  Google Scholar 

  15. Oh J, Byrd AL, Park M, NISC Comparative Sequencing Program, Kong HH, Segre JA. Temporal stability of the human skin microbiome. Cell. 2016;165:854–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Rascovan N, Duraisamy R, Desnues C. Metagenomics and the human virome in asymptomatic individuals. Annu Rev Microbiol. 2016;70:125–41.

    Article  CAS  PubMed  Google Scholar 

  17. Hashida Y, Higuchi T, Matsuzaki S, Nakajima K, Sano S, Daibata M. Prevalence and genetic variability of human polyomaviruses 6 and 7 in healthy skin among asymptomatic individuals. J Infect Dis. 2018;217:483–93.

    Article  PubMed  Google Scholar 

  18. Heidari S, Babor TF, De Castro P, Tort S, Curno M. Sex and gender equity in research: rationale for the SAGER guidelines and recommended use. Res Integr Peer Rev. 2016;1:2.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Hashida Y, Imajoh M, Nemoto Y, Kamioka M, Taniguchi A, Taguchi T, Kume M, Orihashi K, Daibata M. Detection of Merkel cell polyomavirus with a tumour-specific signature in non-small cell lung cancer. Br J Cancer. 2013;108:629–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Seven tips for avoiding DNA contamination in your PCR. Accessed 1 Oct 2020.

  21. Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22:4673–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35:1547–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Basic Local Alignment Search Tool. National Center for Biotechnology Information. Accessed 14 Feb 2023.

  24. Kanda Y. Investigation of the freely available easy-to-use software ‘EZR’ for medical statistics. Bone Marrow Transplant. 2013;48:452–8.

    Article  CAS  PubMed  Google Scholar 

  25. Söderlund-Venermo M. Emerging human parvoviruses: the rocky road to fame. Annu Rev Virol. 2019;6:71–91.

    Article  PubMed  Google Scholar 

  26. Prast-Nielsen S, Tobin AM, Adamzik K, Powles A, Hugerth LW, Sweeney C, Kirby B, Engstrand L, Fry L. Investigation of the skin microbiome: swabs vs. biopsies. Br J Dermatol. 2019;181:572–9.

    Article  CAS  PubMed  Google Scholar 

  27. Mohanraj U, Jokinen M, Thapa RR, Paloniemi M, Vesikari T, Lappalainen M, Tarkka E, Nora-Krūkle Z, Vilmane A, Vettenranta K, Mangani C, Oikarinen S, Fan YM, Ashorn P, Väisänen E, Söderlund-Venermo M. Human protoparvovirus DNA and IgG in children and adults with and without respiratory or gastrointestinal infections. Viruses. 2021;13:483.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Ghazawi FM, Alghazawi N, Le M, Netchiporouk E, Glassman SJ, Sasseville D, Litvinov IV. Environmental and other extrinsic risk factors contributing to the pathogenesis of cutaneous T cell lymphoma (CTCL). Front Oncol. 2019;9:300.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Su C, Tang R, Bai HX, Girardi M, Karakousis G, Zhang PJ, Xiao R, Zhang G. Disease site as a prognostic factor for mycosis fungoides: an analysis of 2428 cases from the US National Cancer Database. Br J Haematol. 2019;185:592–5.

    Article  PubMed  Google Scholar 

  30. Bradford PT, Devesa SS, Anderson WF, Toro JR. Cutaneous lymphoma incidence patterns in the United States: a population-based study of 3884 cases. Blood. 2009;113:5064–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Wilson LD, Hinds GA, Yu JB. Age, race, sex, stage, and incidence of cutaneous lymphoma. Clin Lymphoma Myeloma Leuk. 2012;12:291–6.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Söderlund-Venermo M, Hokynar K, Nieminen J, Rautakorpi H, Hedman K. Persistence of human parvovirus B19 in human tissues. Pathol Biol. 2002;50:307–16.

    Article  PubMed  Google Scholar 

  33. Hokynar K, Söderlund-Venermo M, Pesonen M, Ranki A, Kiviluoto O, Partio EK, Hedman K. A new parvovirus genotype persistent in human skin. Virology. 2002;302:224–8.

    Article  CAS  PubMed  Google Scholar 

  34. Vuorinen T, Lammintausta K, Kotilainen P, Nikkari S. Presence of parvovirus B19 DNA in chronic urticaric and healthy human skin. J Clin Virol. 2002;25:217–21.

    Article  CAS  PubMed  Google Scholar 

  35. Santonja C, Santos-Briz A, Palmedo G, Kutzner H, Requena L. Detection of human parvovirus B19 DNA in 22% of 1815 cutaneous biopsies of a wide variety of dermatological conditions suggests viral persistence after primary infection and casts doubts on its pathogenic significance. Br J Dermatol. 2017;177:1060–5.

    Article  CAS  PubMed  Google Scholar 

  36. Servant A, Laperche S, Lallemand F, Marinho V, De Saint Maur G, Meritet JF, Garbarg-Chenon A. Genetic diversity within human erythroviruses: identification of three genotypes. J Virol. 2002;76:9124–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Hübschen JM, Mihneva Z, Mentis AF, Schneider F, Aboudy Y, Grossman Z, Rudich H, Kasymbekova K, Sarv I, Nedeljkovic J, Tahita MC, Tarnagda Z, Ouedraogo JB, Gerasimova AG, Moskaleva TN, Tikhonova NT, Chitadze N, Forbi JC, Faneye AO, Otegbayo JA, Charpentier E, Muller CP. Phylogenetic analysis of human parvovirus B19 sequences from eleven different countries confirms the predominance of genotype 1 and suggests the spread of genotype 3b. J Clin Microbiol. 2009;47:3735–8.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Jain A, Kant R. Genotypes of erythrovirus B19, their geographical distribution & circulation in cases with various clinical manifestations. Indian J Med Res. 2018;147:239–47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Dobos G, Pohrt A, Ram-Wolff C, Lebbé C, Bouaziz JD, Battistella M, Bagot M, de Masson A. Epidemiology of cutaneous T-cell lymphomas: a systematic review and meta-analysis of 16,953 patients. Cancers. 2020;12:2921.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


Not applicable.


This work was supported by the Japan Society for the Promotion of Science [grant number 22K08434].

Author information

Authors and Affiliations



YH was involved in formal analysis, funding acquisition, investigation, methodology, resources, visualization, and writing—original draft. TH performed investigation. MD was involved in conceptualization, funding acquisition, methodology, project administration, resources, supervision, writing—original draft, and writing—review and editing. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Masanori Daibata.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics Committee of Kochi Medical School, Kochi University. Written informed consent was obtained from all participants. In the case of the enrollment of children aged < 18 or elderly persons with impaired decision-making capacity, their parents or families provided written informed consent on their behalf.

Consent for publication

Not applicable.

Competing interests

The authors have no financial or non-financial interests to disclose.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Additional file 1: table S1 Sequences of primers used for PCR analysis.

Additional file 2: figure S1 Phylogenetic analysis of cutavirus based on complete VP2 sequences.


Additional file 3: table S2 Nucleotide identities between the near-full-length cutavirus sequences (4455 bp) identified in the current study.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hashida, Y., Higuchi, T. & Daibata, M. Cutavirus on the skin in an Asian cohort: identification of a novel geographically related genotype. Virol J 20, 69 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: