Divergent bufavirus harboured in megabats represents a new lineage of parvoviruses

Characterization of a mammalian orthoreovirus isolated from the large flying fox, Pteropus vampyrus, in Indonesia

Fruit bats serve as an important reservoir for many zoonotic pathogens, including Nipah virus, Hendra virus, Marburg virus and Lyssavirus. To gain a deeper insight into the virological characteristics, pathogenicity and zoonotic potential of bat-borne viruses, recovery of infectious viruses from field samples is important. Here, we report the isolation and characterization of a mammalian orthoreovirus (MRV) from a large flying fox (Pteropus vampyrus) in Indonesia, which is the first detection of MRV in Southeast Asia. MRV was recovered from faecal samples of three different P. vampyrus in Central Java. Nucleotide sequence analysis revealed that the genome of the three MRV isolates shared more than 99% nucleotide sequence identity. We tentatively named one isolated strain as MRV12-52 for further analysis and characterization. Among 10 genome segments, MRV12-52 S1 and S4, which encode the cell-attachment protein and outer capsid protein, had 93.6 and 95.1% nucleotide sequence identities with known MRV strains, respectively. Meanwhile, the remaining genome segments of MRV12-52 were divergent with 72.9-80.7 % nucleotide sequence identities. Based on the nucleotide sequence of the S1 segment, MRV12-52 was grouped into serotype 2, and phylogenetic analysis demonstrated evidence of past reassortment events. In vitro characterization of MRV12-52 showed that the virus efficiently replicated in BHK-21, HEK293T and A549 cells. In addition, experimental infection of laboratory mice with MRV12-52 caused severe pneumonia with 75% mortality. This study highlights the presence of pathogenic MRV in Indonesia, which could serve as a potential animal and public health concern.

Structural Characterization of Human Bufavirus 1: Receptor Binding and Endosomal pH-Induced Changes

Bufaviruses (BuV) are members of the Parvoviridae of the Protoparvovirus genus. They are non-enveloped, T = 1 icosahedral ssDNA viruses isolated from patients exhibiting acute diarrhea. The lack of treatment options and a limited understanding of their disease mechanisms require studying these viruses on a molecular and structural level. In the present study, we utilize glycan arrays and cell binding assays to demonstrate that BuV1 capsid binds terminal sialic acid (SIA) glycans. Furthermore, using cryo-electron microscopy (cryo-EM), SIA is shown to bind on the 2/5-fold wall of the capsid surface. Interestingly, the capsid residues stabilizing SIA binding are conserved in all human BuVs identified to date. Additionally, biophysical assays illustrate BuV1 capsid stabilization during endo-lysosomal (pH 7.4-pH 4) trafficking and capsid destabilization at pH 3 and less, which correspond to the pH of the stomach. Hence, we determined the cryo-EM structures of BuV1 capsids at pH 7.4, 4.0, and 2.6 to 2.8 Å, 3.2 Å, and 2.7 Å, respectively. These structures reveal capsid structural rearrangements during endo-lysosomal escape and provide a potential mechanism for this process. The structural insights gained from this study will add to the general knowledge of human pathogenic parvoviruses. Furthermore, the identification of the conserved SIA receptor binding site among BuVs provides a possible targetable surface-accessible pocket for the design of small molecules to be developed as anti-virals for these viruses.

Viral metagenomic analysis reveals diverse viruses and a novel bocaparvovirus in the enteric virome of snow leopard (Panthera uncia)

The enteric virome, comprising a complex community of viruses inhabiting the gastrointestinal tract, plays a significant role in health and disease dynamics. In this study, the fecal sample of a wild snow leopard was subjected to viral metagenomic analysis using a double barcode Illumina MiSeq platform. The resulting reads were de novo assembled into contigs with SOAPdenovo2 version r240. Additional bioinformatic analysis of the assembled genome and genome annotation was done using the Geneious prime software (version 2022.0.2). Following viral metagenomic analysis and bioinformatic analysis, a total of 7 viral families and a novel specie of bocaparvovirus tentatively named Panthera uncia bocaparvovirus (PuBOV) with GenBank accession number OQ627713 were identified. The complete genome of PuBOV was predicted to contain 3 open reading frames (ORFs), contains 5433 nucleotides and has a G + C content of 47.40 %. BLASTx analysis and pairwise sequence comparison indicated the novel virus genome was a new species in the genus Bocaparvovirus based on the species demarcation criteria of the International Committee on the Taxonomy of Viruses. This study provides valuable insights into the diversity and composition of the enteric virome in wild endangered snow leopards. The identification and characterization of viruses in wildlife is crucial for developing effective strategies to manage and mitigate potential zoonotic and other viral disease threats to human and animal health.

Molecular Surveillance for Bocaparvoviruses and Bufaviruses in the European Hedgehog (Erinaceus europaeus)

The presence of bocaparvoviruses (BoVs) and bufaviruses (BuVs) in the European hedgehog (Erinaceus europaeus) was investigated by screening duodenal and liver samples collected from 183 carcasses, delivered to wildlife rescue centers located in northwestern Italy. BoV DNA was detected in 15 animals (8.2%), with prevalences of 7.1% (13/183) and 2.7% (5/183) in intestine and liver samples, respectively. Upon the sequence analyses of the NS1 gene, two highly divergent BoVs (65.5-67.8% nt identities) were identified. Fourteen strains showed the highest identity (98.3-99.4% nt) to the hedgehog BoV strains recently detected in China in Amur hedgehogs (Erinaceus amurensis), whilst four strains were genetically related (98.9-99.4% nt identities) to the porcine BoVs identified in pigs and classified in the species Bocaparvovirus ungulate 4, which included related viruses also found in rats, minks, shrews, and mice. BuV DNA was detected in the duodenal samples of two hedgehogs, with a prevalence rate of 1.1%. The nearly full-length genome of two BuV strains, Hedgehog/331DU-2022/ITA and Hedgehog/1278DU/2019/ITA, was reconstructed. Upon phylogenetic analysis based on the NS and VP aa sequences, the Italian hedgehog BuVs tightly clustered with the BuVs recently identified in the Chinese Amur hedgehogs, within a potential novel candidate species of the genus Protoparvovirus.

Isolation and characterization of an orthoreovirus from Indonesian fruit bats

Nelson Bay orthoreovirus (NBV) is an emerging bat-borne virus and causes respiratory tract infections in humans sporadically. Over the last two decades, several strains genetically related to NBV were isolated from humans and various bat species, predominantly in Southeast Asia (SEA), suggesting a high prevalence of the NBV species in this region. In this study, an orthoreovirus (ORV) belonging to the NBV species was isolated from Indonesian fruit bats' feces, tentatively named Paguyaman orthoreovirus (PgORV). Serological studies revealed that 81.2% (108/133) of Indonesian fruit bats sera had neutralizing antibodies against PgORV. Whole-genome sequencing and phylogenetic analysis of PgORV suggested the occurrence of past reassortments with other NBV strains isolated in SEA, indicating the dispersal and circulation of NBV species among bats in this region. Intranasal PgORV inoculation of laboratory mice caused severe pneumonia. Our study characterized PgORV's unique genetic background and highlighted the potential risk of PgORV-related diseases in Indonesia.

First report of canine bufavirus in India

Canine bufavirus (CBuV), a novel protoparvovirus of dogs that is associated with enteric and respiratory symptoms, has been reported only in Italy and China. The enteric prevalence of CBuV in India was investigated, and the nearly complete genome sequence (4292 bp) was amplified and reconstructed for one strain. A nucleotide sequence alignment indicated 93.42–98.81% identity to the other available CBuV sequences and 70.88–73.39% and 54.4–54.8% identity to human bufavirus and canine parvovirus 2 (CPV-2), respectively. The current strain is most closely related to Chinese CBuV strains, which together form an Asian lineage. This first report of the prevalence of CBuV in India emphasizes the need for further epidemiological surveillance.

[Investigation of viruses harbored by wild animals: toward pre-emptive measures against future zoonotic diseases]

Zoonoses are caused by pathogens transmitted from animals. To prepare mitigating measures against emerging zoonoses, it is imperative to identify animal reservoirs that carry potential pathogens and also elucidate the transmission routes of these pathogens. Under the continuous collaboration with counterparts from Zambia and Indonesia, we have so far identified various viruses in wild animals. Some of the identified viruses were phylogenetically distinct from known virus species and this finding led to approved new virus species by the International Committee on Taxonomy of Viruses (ICTV). Our studies provided new insights into the divergence, natural hosts and lifecycle of viruses. Through the exploration and characterization of viruses in animals, we will endeavor to contribute to the existing knowledge on viral pathogens in wild animals. This is cardinal for evidence-based preemptive measures against future zoonoses.

Small but mighty: old and new parvoviruses of veterinary significance

In line with the Latin expression "sed parva forti" meaning "small but mighty," the family Parvoviridae contains many of the smallest known viruses, some of which result in fatal or debilitating infections. In recent years, advances in metagenomic viral discovery techniques have dramatically increased the identification of novel parvoviruses in both diseased and healthy individuals. While some of these discoveries have solved etiologic mysteries of well-described diseases in animals, many of the newly discovered parvoviruses appear to cause mild or no disease, or disease associations remain to be established. With the increased use of animal parvoviruses as vectors for gene therapy and oncolytic treatments in humans, it becomes all the more important to understand the diversity, pathogenic potential, and evolution of this diverse family of viruses. In this review, we discuss parvoviruses infecting vertebrate animals, with a special focus on pathogens of veterinary significance and viruses discovered within the last four years.

Genetic and phylogenetic analysis of canine bufavirus from Anhui Province, Eastern China

Bufavirus is a novel virus associated with canine gastroenteritis. Three strains of bufavirus were first detected in dog feces collected from Anhui province in Eastern China. The near-complete genome sequences were amplified. Sequence alignment showed 98.3-99.5% homology between the three bufavirus strains and reference strains. Phylogenetic analysis showed the distributed viruses forming a cluster of close relationships. Selective pressure analysis of the VP2 region indicated that the canine bufavirus (CBuV) was mainly subject to negative selection during evolution. The negative selection site was located on the residue of B-cell epitopes, indicating minimal change to the virus's immunogenicity. Since this is the first report of CBuV circulating in Anhui Province, this study will provide further understanding of the phylogenetic and molecular characteristics of CBuV and serve as a reference for prevention and vaccine development.

Genetic heterogeneity of canine bufaviruses

Canine bufavirus (CBuV) is a protoparvovirus, genetically related to human and non-human primate bufaviruses and distantly related to canine parvovirus type 2 (CPV-2). CBuV was initially identified from young dogs with respiratory signs but subsequent studies revealed that this virus is also a common component of the canine enteric virome. In this survey, by assessing archival and recent collections of dogs faecal samples, CBuV DNA was detected with a higher prevalence rate (8.8%) in animals with enteritis than in control animals (5.0%), although this difference was not statistically significant. The rate of co-infections with other enteric viruses in diarrhoeic dogs was high (84.6%), mostly in association with canine parvovirus CPV-2 (90.1%). The complete ORF2 gene was determined in five samples, and the nearly full-length genome was reconstructed for three strains, 62/2017/ITA, 9AS/2005/ITA and 35/2018/ITA. Upon sequence comparison, the viruses appeared highly conserved in the NS1 (97.2%-97.9% nt and 97.5%-98.1% aa identities). In the complete VP2 coding region, three strains were similar to the prototype viruses (99.7-99.8 nt and 99.6%-99.8% aa) whilst strains 9AS/2005/ITA and 35/2016/ITA were distantly related (87.6%-89.3% nt and 93.9%-95.1% aa identities). Interestingly, genetic diversification occurred downstream conserved regions such as the VP1/VP2 splicing signals and/or the G-rich motif in the N terminus of the VP2, suggesting a potential recombination nature. Upon phylogenetic analysis, the two divergent CBuV strains formed a distinct cluster/genotype.

First report and genetic diversity of porcine bufavirus in China

BACKGROUND

Bufavirus is a newly discovered zoonotic virus reported in numerous mammals and humans. However, the epidemiological and genetic characteristics of porcine bufaviruses (PBuVs) in China remain unclear.

METHODS

To detect PBuVs in China, 384 samples (92 fecal and 292 serum specimens) were collected from 2017 to 2018, covering six provinces in China, and were evaluated by nested PCR. Further, the positive samples from different provinces were selected to obtain the complete genome of Chinese PBuVs.

RESULTS

The prevalence rate of PBuV was 16.7% in Chinese domestic pigs in the Guangdong, Guangxi, Fujian, Jiangxi, Anhui, and Henan provinces. Additionally, the positive rate of fecal specimens was higher than that of the serum samples. Next, we sequenced nine near-complete genomes of Chinese field PBuV strains from different provinces. Homology and phylogenetic analyses indicated that Chinese PBuVs have high genetic variation (93.3-99.2%), showed higher nucleotide identity with an Austrian PBuV strain (KU867071.1), and developed into different branches within the same cluster.

CONCLUSION

To our knowledge, this is the first report on PBuV in China, expanding the geographic boundaries of PBuV circulation. Our data demonstrate that PBuVs are widely distributed in the six Chinese provinces. Moreover, these Chinese PBuVs exhibit genetic variation and continuous evolution characteristics. Taken together, our findings provide a foundation for future studies on bufaviruses.

Genomic characterization of human adenovirus type 4 strains isolated worldwide since 1953 identifies two separable phylogroups evolving at different rates from their most recent common ancestor

Species Human mastadenovirus E (HAdV-E) comprises several simian types and a single human type: HAdV-E4, a respiratory and ocular pathogen. RFLP analysis for the characterization of intratypic genetic variability has previously distinguished two HAdV-E4 clusters: prototype (p)-like and a-like. Our analysis of whole genome sequences confirmed two distinct lineages, which we refer to as phylogroups (PGs). PGs I and II comprise the p- and a-like genomes, respectively, and differ significantly in their G + C content (57.7% ± 0.013 vs 56.3% ± 0.015). Sequence differences distinguishing the two clades map to several regions of the genome including E3 and ITR. Bayesian analyses showed that the two phylogroups diverged approximately 602 years before the present. A relatively faster evolutionary rate was identified for PG II. Our data provide a rationale for the incorporation of phylogroup identity to HAdV-E4 strain designation to reflect the identified unique genetic characteristics that distinguish PGs I and II.

The spleen morphophysiology of fruit bats

Spleen is one of the important lymphoid organs with wide variations of morphological and physiological functions according to species. Morphology and function of the spleen in bats, which are hosts to several viral strains without exhibiting clinical symptoms, remain to be fully elucidated. This study aims to examine the spleen morphology of fruit bats associated with their physiological functions. Spleen histological observations were performed in three fruit bats species: Cynopterus titthaecheilus (n = 9), Rousettus leschenaultii (n = 3) and Pteropus vampyrus (n = 3). The spleens of these fruit bats were surrounded by a thin capsule. Red pulp consisted of splenic cord and wide vascular space filled with blood. Ellipsoids in all three studied species were found numerously and adjacent to one another forming macrophages aggregates. White pulp consisted of periarteriolar lymphoid sheaths (PALS), lymphoid follicles and marginal zone. The lymphoid follicle contained a germinal centre and a tingible body macrophage that might reflect an active immune system. The marginal zone was prominent and well developed. This study reports some differences in spleen structure of fruit bats compared to other bat species previously reported and discusses possible physiological implications of the spleen based on its morphology.

Novel Parvovirus Related to Primate Bufaviruses in Dogs

A novel protoparvovirus species, related genetically to human bufaviruses, was identified in dogs with respiratory signs. The canine bufavirus was distantly related to the well-known canine protoparvovirus, canine parvovirus type 2, sharing low amino acid identities in the nonstructural protein 1 (40.6%) and in the capsid protein 1 (33.4%). By screening collections of fecal, nasal, and oropharyngeal samples obtained from juvenile dogs (<1 year of age), canine bufavirus DNA appeared as a common component of canine virome. The virus was common in the stool samples of dogs with or without enteric disease and in the nasal and oropharyngeal swab samples of dogs with respiratory signs. However, the virus was not detected in nasal and oropharyngeal swab samples from animals without clinical signs.

Global Distribution of Human Protoparvoviruses

Development of next-generation sequencing and metagenomics has revolutionized detection of novel viruses. Among these viruses are 3 human protoparvoviruses: bufavirus, tusavirus, and cutavirus. These viruses have been detected in feces of children with diarrhea. In addition, cutavirus has been detected in skin biopsy specimens of cutaneous T-cell lymphoma patients in France and in 1 melanoma patient in Denmark. We studied seroprevalences of IgG against bufavirus, tusavirus, and cutavirus in various populations (n = 840), and found a striking geographic difference in prevalence of bufavirus IgG. Although prevalence was low in adult populations in Finland (1.9%) and the United States (3.6%), bufavirus IgG was highly prevalent in populations in Iraq (84.8%), Iran (56.1%), and Kenya (72.3%). Conversely, cutavirus IgG showed evenly low prevalences (0%–5.6%) in all cohorts, and tusavirus IgG was not detected. These results provide new insights on the global distribution and endemic areas of protoparvoviruses.

First identification of a novel parvovirus distantly related to human bufavirus from diarrheal dogs in China

Bufaviruses are small, nonenveloped, single-stranded DNA viruses belonging to the subfamily Parvovirinae. Human bufaviruses were first identified in 2012 in fecal samples from children with diarrhea. A new parvovirus of canines that was first detected in various samples from dogs with enteric and respiratory symptoms in Italy between 2014 and 2018 is closely related to the newly described human bufavirus. To explore the prevalence and genetic diversity of CBuV in Chinese dogs, 540 canine parvovirus (CPV)-positive serum and diarrhea samples were collected in Guangxi Province between 2016 and 2018. Among the samples, 6.25% (5/80) of rectal swabs and 2.5% (5/200) of CPV PCR-positive samples were positive for CBuV. However, the virus was not detected in CPV PCR-negative samples or nasal swabs. Two CBuV isolates were identified from CPV-positive fecal and serum samples by complete sequence analysis, with 99.8%-99.9% NS1 and VP2 protein identity to each another. Sequence analysis indicated that the CBuV GXNN01-2018 isolate VP2 protein shares 99.6% identity with the Italian CBuV ITA/2015/297 isolate and 62.3%-65.5% identity with human bufavirus. Phylogenetic analysis showed that CBuV was significantly distinct from other known bufaviruses and was most closely related to CBuV ITA/2015/297. This is the first report of the existence of CBuV in China, and our findings will strengthen the understanding of the epidemiology of bufaviruses in different animals.

Detection of adenovirus, papillomavirus and parvovirus in Brazilian bats of the species Artibeus lituratus and Sturnira lilium

Bats play a significant role in maintaining their ecosystems through pollination, dispersal of seeds, and control of insect populations, but they are also known to host many microorganisms and have been described as natural reservoirs for viruses with zoonotic potential. The diversity of viruses in these animals remains largely unknown, however, because studies are limited by species, location, virus target, or sample type. Therefore, the aim of this study was to detect fragments of viral genomes in bat samples. We performed high-throughput sequencing analysis and specific PCR and RT-PCR on pools of anal and oropharyngeal swabs from Artibeus lituratus and Sturnira lilium collected in southern Brazil. As a result, a member of the family Adenoviridae related to human adenovirus C was detected in anal swabs from S. lilium. In addition, we detected a papillomavirus in an anal swab from A. lituratus. Our analyses also allowed the detection of adenoviruses and parvoviruses in oropharyngeal swabs collected from A. lituratus. These results increase our knowledge about viral diversity and illustrate the importance of conducting virus surveillance in bats.

Novel Parvoviruses from Wild and Domestic Animals in Brazil Provide New Insights into Parvovirus Distribution and Diversity

Parvoviruses (family Parvoviridae) are small, single-stranded DNA viruses. Many parvoviral pathogens of medical, veterinary and ecological importance have been identified. In this study, we used high-throughput sequencing (HTS) to investigate the diversity of parvoviruses infecting wild and domestic animals in Brazil. We identified 21 parvovirus sequences (including twelve nearly complete genomes and nine partial genomes) in samples derived from rodents, bats, opossums, birds and cattle in Pernambuco, São Paulo, Paraná and Rio Grande do Sul states. These sequences were investigated using phylogenetic and distance-based approaches and were thereby classified into eight parvovirus species (six of which have not been described previously), representing six distinct genera in the subfamily Parvovirinae. Our findings extend the known biogeographic range of previously characterized parvovirus species and the known host range of three parvovirus genera (Dependovirus, Aveparvovirus and Tetraparvovirus). Moreover, our investigation provides a window into the ecological dynamics of parvovirus infections in vertebrates, revealing that many parvovirus genera contain well-defined sub-lineages that circulate widely throughout the world within particular taxonomic groups of hosts.

Bats host diverse parvoviruses as possible origin of mammalian dependoparvoviruses and source for bat-swine interspecies transmission

Compared to the enormous species diversity of bats, relatively few parvoviruses have been reported. We detected diverse and potentially novel parvoviruses from bats in Hong Kong and mainland China. Parvoviruses belonging to Amdoparvovirus, Bocaparvovirus and Dependoparvovirus were detected in alimentary, liver and spleen samples from 16 different chiropteran species of five families by PCR. Phylogenetic analysis of partial helicase sequences showed that they potentially belonged to 25 bocaparvovirus, three dependoparvovirus and one amdoparvovirus species. Nearly complete genome sequencing confirmed the existence of at least four novel bat bocaparvovirus species (Rp-BtBoV1 and Rp-BtBoV2 from Rhinolophus pusillus, Rs-BtBoV2 from Rhinolophus sinicus and Rol-BtBoV1 from Rousettus leschenaultii) and two novel bat dependoparvovirus species (Rp-BtAAV1 from Rhinolophus pusillus and Rs-BtAAV1 from Rhinolophus sinicus). Rs-BtBoV2 was closely related to Ungulate bocaparvovirus 5 with 93, 72.1 and 78.7 % amino acid identities in the NS1, NP1 and VP1/VP2 genes, respectively. The detection of bat bocaparvoviruses, including Rs-BtBoV2, closely related to porcine bocaparvoviruses, suggests recent interspecies transmission of bocaparvoviruses between bats and swine. Moreover, Rp-BtAAV1 and Rs-BtAAV1 were most closely related to human AAV1 with 48.7 and 57.5 % amino acid identities in the rep gene. The phylogenetic relationship between BtAAVs and other mammalian AAVs suggests bats as the ancestral origin of mammalian AAVs. Furthermore, parvoviruses of the same species were detected from multiple bat species or families, supporting the ability of bat parvoviruses to cross species barriers. The results extend our knowledge on the diversity of bat parvoviruses and the role of bats in parvovirus evolution and emergence in humans and animals.

[Discovery of DNA viruses in wildlife in Zambia and Indonesia]

Zoonoses originate from pathogens harbored in domestic and wild animals and therefore it is likely impossible to completely eradicate zoonotic diseases. For pre-emptive measures to attempt to predict the emergence of zoonosis outbreaks and the prevention of future epidemics and pandemics, it is imperative to identify natural host animals carrying potential pathogens and elucidate the routes of pathogen transmission into the human population. Our research team is conducting epidemiological research studies in Zambia and Indonesia for the control of viral zoonotic diseases. In this review, we present the research findings, including the discovery of orthopoxviruses and polyomaviruses in wildlife in Zambia and the identification of herpesviruses in bats in Indonesia among our activities.

Epidemiology of two human protoparvoviruses, bufavirus and tusavirus

Two human parvoviruses were recently discovered by metagenomics in Africa, bufavirus (BuV) in 2012 and tusavirus (TuV) in 2014. These viruses have been studied exclusively by PCR in stool and detected only in patients with diarrhoea, although at low prevalence. Three genotypes of BuV have been identified. We detected, by in-house EIA, BuV1-3 IgG antibodies in 7/228 children (3.1%) and 10/180 adults (5.6%), whereas TuV IgG was found in one child (0.4%). All children and 91% of the adults were Finnish, yet interestingly 3/6 adults of Indian origin were BuV-IgG positive. By competition EIA, no cross-reactivity between the BuVs was detected, indicating that the BuV genotypes represent distinct serotypes. Furthermore, we analysed by BuV qPCR stool and nasal swab samples from 955 children with gastroenteritis, respiratory illness, or both, and found BuV DNA in three stools (0.3%) and for the first time in a nasal swab (0.1%). This is the first study documenting the presence of BuV and TuV antibodies in humans. Although the seroprevalences of both viruses were low in Finland, our results indicate that BuV infections might be widespread in Asia. The BuV-specific humoral immune responses appeared to be strong and long-lasting, pointing to systemic infection in humans.

A new protoparvovirus in human fecal samples and cutaneous T cell lymphomas (mycosis fungoides)

We genetically characterized seven nearly complete genomes in the protoparvovirus genus from the feces of children with diarrhea. The viruses, provisionally named cutaviruses (CutaV), varied by 1-6% nucleotides and shared ~76% and ~82% amino acid identity with the NS1 and VP1 of human bufaviruses, their closest relatives. Using PCR, cutavirus DNA was found in 1.6% (4/245) and 1% (1/100) of diarrhea samples from Brazil and Botswana respectively. In silico analysis of pre-existing metagenomics datasets then revealed closely related parvovirus genomes in skin biopsies from patients with epidermotropic cutaneous T-cell lymphoma (CTCL or mycosis fungoides). PCR of skin biopsies yielded cutavirus DNA in 4/17 CTCL, 0/10 skin carcinoma, and 0/21 normal or noncancerous skin biopsies. In situ hybridization of CTCL skin biopsies detected viral genome within rare individual cells in regions of neoplastic infiltrations. The influence of cutavirus infection on human enteric functions and possible oncolytic role in CTCL progression remain to be determined.