Papillomavirus in healthy skin of Australian animals

Detection and quantification of bovine papillomavirus DNA by digital droplet PCR in sheep blood

Highly pathogenic bovine papillomaviruses (BPVs) were detected and quantified for the first time using digital droplet polymerase chain reaction (ddPCR) by liquid biopsy in 103 clinically healthy sheep. Overall, ddPCR detected BPVs in 68 blood samples (66%). BPV infection by a single genotype was revealed in 61.8% of the blood samples, and BPV coinfection by double, triple or quadruple genotypes was observed in 38.2% of liquid biopsies. The BPV-2 genotype was most frequently seen in sheep, whereas BPV-1 was the least common. Furthermore, ddPCR was very useful for detection and quantification; the BPV-14 genotype was observed for the first time in ovine species, displaying the highest prevalence in some geographical areas (Apulia). In 42 of the positive samples (61.8%), a single BPV infection was observed, 26 of which were caused by BPV-2 (61.9%) and 7 by BPV-13 (16.7%). BPV-14 was responsible for 7 single infections (16.7%) and BPV-1 for 2 single infections (4.7%). Multiple BPV coinfections were observed in the remaining 26 positive samples (38.2%), with dual BPV-2/BPV-13 infection being the most prevalent (84.6%). BPV infection by triple and quadruple genotypes was also observed in 11.5% and 3.8% of cases, respectively. The present study showed that ddPCR, a biotechnological refinement of conventional PCR, is by far the most sensitive and accurate assay for BPV detection compared to conventional qPCR. Therefore, ddPCR displayed an essential diagnostic and epidemiological value very useful for the identification of otherwise undetectable BPV genotypes as well as their geographical distributions and suggesting that animal husbandry practices contribute to cross-species transmission of BPVs.

Squamous cell carcinoma in a digit of the hind limb with systemic metastasis in a 17-year-old female koala

We encountered a case of cutaneous squamous cell carcinoma (SCC) in a 17-year-old female koala at a zoo. A fragile, papillary, elevated mass was found on the third digit of the right hind limb. SCC was identified histopathologically: squamous cell-like polygonal tumor cells showed a nest-like growth pattern with epidermal down growth, central keratinization and necrotic foci, and invaded dermal connective tissues. Metastatic lesions were observed in various organs, including the lung and axillary lymph node: in the lung, multiple metastatic foci similar to the primary lesion, and in the axillary lymph node, individual polygonal tumor cells infiltrated the sinusoids. Immunohistochemistry revealed that the tumor cells were positive for proliferating cell nuclear antigen, which exhibited 32–33% of labeling indices in the tumor cells. To our knowledge, this is the first report of a case of SCC in a digit of a koala.

Leopardus wiedii Papillomavirus type 1, a novel papillomavirus species in the tree ocelot, suggests Felidae Lambdapapillomavirus polyphyletic origin and host-independent evolution

The limited knowledge on Papillomavirus diversity (particularly in wild animal species) influences the accuracy of PVs phylogeny and their evolutionary history, and hinders the comprehension of PVs pathogenicity, especially the mechanism of virus - related cancer progression. This study reports the identification of Leopardus wiedii Papillomavirus type 1 (LwiePV1), the first PV type within Lambdapapillomavirus in a Leopardus host. LwiePV1 full genome sequencing allowed the investigation of its taxonomic position and phylogeny. Based on results, LwiePV1 should be assigned to a novel PV species providing evidence for a polyphyletic origin of feline lambda PVs, and representing an exception to codivergence between feline lambda PVs and their hosts. Results improve our knowledge on PV diversity and pave the way to future studies investigating biological and evolutionary features of animal PVs.

Origin and evolution of papillomavirus (onco)genes and genomes

Papillomaviruses (PVs) are ancient viruses infecting vertebrates, from fishes to mammals. Although the genomes of PVs are small and show conserved synteny, PVs display large genotypic diversity and ample variation in the phenotypic presentation of the infection. Most PV genomes contain two small early genes E6 and E7. In a bunch of closely related human papillomaviruses (HPVs), the E6 and E7 proteins provide the viruses with oncogenic potential. The recent discoveries of PVs without E6 and E7 in different fish species place a new root on the PV tree, and suggest that ancestral PVs consisted of the minimal PV backbone E1-E2-L2-L1. Bayesian phylogenetic analyses date the most recent common ancestor of the PV backbone to 424 million years ago (Ma). Common ancestry tests on extant E6 and E7 genes indicate that they share a common ancestor dating back to at least 184 Ma. In AlphaPVs infecting Old World monkeys and apes, the appearance of the E5 oncogene 53-58 Ma concurred with (i) a significant increase in substitution rate, (ii) a basal radiation and (iii) key gain of functions in E6 and E7. This series of events was instrumental to construct the extant phenotype of oncogenic HPVs. Our results assemble the current knowledge on PV diversity and present an ancient evolutionary timeline punctuated by evolutionary innovations in the history of this successful viral family. This article is part of the theme issue 'Silent cancer agents: multi-disciplinary modelling of human DNA oncoviruses'.

New Insight Into Avian Papillomavirus Ecology and Evolution From Characterization of Novel Wild Bird Papillomaviruses

Viruses in the family Papillomaviridae have circular dsDNA genomes of approximately 5.7–8.6 kb that are packaged within non-enveloped, icosahedral capsids. The known papillomavirus (PV) representatives infect vertebrates, and there are currently more than 130 recognized PV species in more than 50 genera. We identified 12 novel avian papillomavirus (APV) types in wild birds that could represent five distinct species and two genera. Viruses were detected in paired oropharyngeal/cloacal swabs collected from six bird species, increasing the number of avian species known to harbor PVs by 40%. A new duck PV (DuPV-3) was found in mallard and American black duck (27.6% estimated prevalence) that was monophyletic with other known DuPVs. A single viral type was identified in Atlantic puffin (PuPV-1, 9.8% estimated prevalence), while a higher genetic diversity was found in other Charadriiformes. Specifically, three types [gull PV-1 (GuPV-1), -2, and -3] were identified in two gull species (estimated prevalence of 17% and 2.6% in American herring and great black-backed gull, respectively), and seven types [kittiwake PV-1 (KiPV-1) through -7] were found in black-legged kittiwake (81.3% estimated prevalence). Significantly higher DuPV-3 circulation was observed in spring compared to fall and in adults compared to juveniles. The studied host species’ tendencies to be in crowded environments likely affect infection rates and their migratory behaviors could explain the high viral diversity, illustrating how host behavior can influence viral ecology and distribution. For DuPV-3, GuPV-1, PuPV-1, and KiPV-2, we obtained the complete genomic sequences, which showed the same organization as other known APVs. Phylogenetic analyses showed evidence for virus–host co-divergence at the host taxonomic levels of family, order, and inter-order, but we also observed that host-specificity constraints are relaxed among highly related hosts as we found cross-species transmission within ducks and within gulls. Furthermore, the phylogeny of viruses infecting the Charadriiformes did not match the host phylogeny and gull viruses formed distinct monophyletic clades with kittiwake viruses, possibly reflecting past host-switching events. Considering the vast PV genotype diversity in other hosts and the large number of bird species, many more APVs likely remain to be discovered.

Adaptation of Alpha-Papillomavirus over Millennia

Papillomaviruses (PVs) are a group of small DNA viruses that, with around 350 million years of evolution, acquired the capacity of infecting a broad range of vertebrates, including humans. To date, more than 300 PV types have been isolated. Viruses that have a long common evolutionary history with their host typically cause unapparent infections. However, in some Alpha-PV infections, lesions become apparent and may cause benign proliferative disorders or even malignant proliferative lesions of the cervix, vulva, vagina, anus, penis, and oropharynx. The incongruence observed between the topology of the phylogenetic tree of Alpha-PVs and that of their hosts suggests that virus-host codivergence is not the only evolutionary force that has driven the progression of PVs. The integration of the precursors of E5, E6, and E7 on the genome of the ancestral Alpha-PV was important and made the colonization of new niches and the emergence of carcinogenic types possible.

Oral Squamous Cell Carcinoma in Three Related Kowari (Dasyuroides byrnei)

We report three kowari (Dasyuroides byrnei) with squamous cell carcinoma affecting the gingiva. These cases occurred in rapid succession in a related group of individuals of similar age, suggesting a familial tendency to this condition and a typical age of presentation. Other conditions affecting the oral cavity can mimic the appearance of oral squamous cell carcinoma in this species, and so knowledge of this condition can assist the veterinarian in making rapid decisions regarding prognosis and improving the welfare of these animals.

Searching beyond the usual papillomavirus suspects in squamous carcinomas of the vulva, penis and head and neck

Human Papillomaviruses (HPVs) are involved in the etiology of anogenital and head and neck cancers. The HPV DNA prevalence greatly differs by anatomical site. Indeed, the high rates of viral DNA prevalence in anal and cervical carcinomas contrast with the lower fraction of cancer cases attributable to HPVs in other anatomical sites, chiefly the vulva, the penis and head and neck. Here we analyzed 2635 Formalin Fixed Paraffin Embedded surgical samples that had previously tested negative for the presence of HPVs DNA using the SPF10/DEIA procedure, in order to identify the presence of other PVs not explicitly targeted by standard molecular epidemiologic approaches. All samples were reanalyzed using five broad-PV PCR primer sets (CP1/2, FAP6064/FAP64, SKF/SKR, MY9/MY11, MFI/MFII) targeting the main PV main clades. In head and neck carcinoma samples (n=1141), we recovered DNA from two BetaHPVs, namely HPV20 and HPV21, and from three cutaneous AlphaPVs, namely HPV2, HPV57 and HPV61. In vulvar squamous cell carcinoma samples (n=902), we found one of the samples containing DNA of one cutaneous HPV, namely HPV2, and 29 samples contained DNA from essentially mucosal HPVs. In penile squamous cell carcinoma samples (n=592), we retrieved the DNA of HPV16 in 16 samples. Our results show first that the SPF10/DEIA is very sensitive, as we recovered only 2.1% (55/2635) false negative results; second, that although the DNA of cutaneous HPVs can be detected in cancer samples, their relative contribution remains anyway minor (0.23%; 6/2635) and may be neglected for screening and vaccination purposes; and third, their contribution to malignancy is not necessarily warranted and needs to be elucidated.

Geno- and seroprevalence of Felis domesticus Papillomavirus type 2 (FdPV2) in dermatologically healthy cats

Background

Papillomaviruses can cause proliferative skin lesions ranging from benign hyperplasia to squamous cell carcinoma (SCC). However, asymptomatic infection is also possible. Several groups have detected Felis domesticus Papillomavirus type 2 (FdPV2) DNA in association with feline Bowenoid in situ carcinoma (BISC). Therefore, a causative connection has been suggested. However, the knowledge about FdPV2 epidemiology is limited. The aim of this study was to describe the genoprevalence and seroprevalence of FdPV2 in healthy cats.

For this purpose an FdPV2-specific quantitative (q)PCR assay was developed and used to analyse Cytobrush samples collected from 100 dermatologically healthy cats. Moreover, an ELISA was established to test the sera obtained from the same cats for antibodies against the major capsid protein (L1) of FdPV2.

Results

The genoprevalence of FdPV2 was to 98 %. Surprisingly, the quantities of viral DNA detected in some samples from the healthy cats exceeded the amounts detected in control samples from feline BISC lesions. The seroprevalence was much lower, amounting to 22 %. The concentrations of antibodies against FdPV2 were relatively low in healthy cats, whereas they were very high in control cats with BISC.

Conclusion

These observations suggest that FdPV2 is highly prevalent, even among healthy cats. However, cats that carry it on their skin mount in most instances no antibody response. It might be hypothesized that FdPV2 is only rarely productively replicating or its replication is only rarely exposed to the immune system.

Electronic supplementary material

The online version of this article (doi:10.1186/s12917-016-0776-7) contains supplementary material, which is available to authorized users.

Novel Papillomaviral Sequence Detected within Epidermal Plaques in a Wolf (Canis lupus)

We describe numerous pale plaques affecting the inguinal skin of a grey wolf (Canis lupus). Histologically, these were consistent with papillomaviral plaques. Immunohistochemistry confirmed papillomavirus antigens, and partial sequencing of the L1 gene suggests this is a novel papillomavirus most-closely related to Canis familiaris Papillomavirus 5.

Human papillomavirus molecular biology and disease association

Human papillomaviruses (HPVs) have evolved over millions of years to propagate themselves in a range of different animal species including humans. Viruses that have co‐evolved slowly in this way typically cause chronic inapparent infections, with virion production in the absence of apparent disease. This is the case for many Beta and Gamma HPV types. The Alpha papillomavirus types have however evolved immunoevasion strategies that allow them to cause persistent visible papillomas. These viruses activate the cell cycle as the infected epithelial cell differentiates in order to create a replication competent environment that allows viral genome amplification and packaging into infectious particles. This is mediated by the viral E6, E7, and E5 proteins. High‐risk E6 and E7 proteins differ from their low‐risk counterparts however in being able to drive cell cycle entry in the upper epithelial layers and also to stimulate cell proliferation in the basal and parabasal layers. Deregulated expression of these cell cycle regulators underlies neoplasia and the eventual progression to cancer in individuals who cannot resolve high‐risk HPV infection. Most work to date has focused on the study of high‐risk HPV types such as HPV 16 and 18, which has led to an understanding of the molecular pathways subverted by these viruses. Such approaches will lead to the development of better strategies for disease treatment, including targeted antivirals and immunotherapeutics. Priorities are now focused toward understanding HPV neoplasias at sites other than the cervix (e.g. tonsils, other transformation zones) and toward understanding the mechanisms by which low‐risk HPV types can sometimes give rise to papillomatosis and under certain situations even cancers. Copyright © 2015 John Wiley & Sons, Ltd.

Novel animal papillomavirus sequences and accurate phylogenetic placement

All amniotes are probably infected by specific papillomaviruses (PVs), but knowledge about PV diversity remains sparse. An insufficient taxon sampling, and a focus on humans as hosts, may perturb phylogenetic analyses leading to wrong conclusions about PV evolution. We performed a systematic approach to explore the diversity of PVs combining rolling circle amplification with the use of "universal" primers to search for the presence of novel PV sequences in animal samples. We communicate 12 sequences putatively corresponding to novel PVs gained from 10 host species in eight mammal families: Bovidae, Canidae, Cervidae, Equidae, Hominidae, Phocoenidae, Procyonidae and Pteropodidae. The phylogenetic position of the new sequences was inferred with an evolutionary placement algorithm under a Maximum Likelihood framework using a pre-computed, well-resolved tree constructed with the E1-E2-L1 gene sequences as a backbone. The new sequences were phylogenetically diverse and could be respectively placed with confidence within all four PV crown groups. The prevailing presence of sequences from the crown groups Alpha+Omikron-PVs and Beta+Xi-PVs may correspond to an increased viral diversity in these taxa, or rather reflect a combination of anthropocentric bias and preferential amplification from commonly used "universal" primers. Our results combined with literature data support the view that the number and diversity of animal PVs is overwhelmingly large.

Evidence for an endogenous papillomavirus-like element in the platypus genome

Papillomaviruses (PVs) infect a wide range of vertebrates and have diversified into multiple genetic types, some of which have serious consequences for human health. Although PVs have to date only been characterized as exogenous viral forms, here we report the observation of an endogenous viral element (EPVLoa) in the genome of the platypus (Ornithorhynchus anatinus) that is related to PVs. Further data mining for endogenous PV-like elements is therefore warranted.

Evidence of recombination and positive selection in cetacean papillomaviruses

Papillomaviruses (PVs) are small DNA viruses that have been associated with increased epithelial proliferation. Over one hundred PV types have been identified in humans; however, only three have been identified in bottlenose dolphins (Tursiops truncatus) to date. Using rolling circle amplification and degenerate PCR, we identified four novel PV genomes of bottlenose dolphins. TtPV4, TtPV5 and TtPV6 were identified in genital lesions while TtPV7 was identified in normal genital mucosa. Bayesian analysis of the full-length L1 genes found that TtPV4 and TtPV7 group within the Upsilonpapillomavirus genus while TtPV5 and TtPV6 group with Omikronpapillomavirus. However, analysis of the E1 gene did not distinguish these genera, implying that these genes may not share a common history, consistent with recombination. Recombination analyses identified several probable events. Signals of positive selection were found mostly in the E1 and E2 genes. Recombination and diversifying selection pressures constitute important driving forces of cetacean PV evolution.

Characterization of a novel papillomavirus species (ZcPV1) from two California sea lions (Zalophus californianus)

A seven-year old California sea lion (Zalophus californianus) presented with focally extensive, bilaterally symmetric, proliferative axillary skin lesions and preputial lesions. A second California sea lion in the same population presented with similar proliferative lesions on the underside of the tail. Histopathology revealed epidermal hyperplasia with severe hyperkeratosis, with proliferating keratinocytes forming broad, branching pegs that extended into the dermis. Pan-papillomaviral consensus PCR was used to obtain initial E1 sequence template and the complete genome was determined using a combination of rolling circle amplification and specific-primer PCR. Analysis revealed a novel papillomavirus, Zalophus californianus papillomavirus 1 (ZcPV1), with seven open reading frames encoding five early proteins (E6, E7, E1, E2 and E4) and two late proteins (L1 and L2). Phylogenetic analysis revealed that (ZcPV1) is most closely related to Equine papillomavirus 1 (EcPV1) in the genus Zetapapillomavirus, and Canine papillomaviruses 3 and 4 (CPV3, CPV4) in the genus Chipapillomavirus. The lesions regressed without intervention over a period of several months.

Clinically healthy skin of dogs is a potential reservoir for canine papillomaviruses

Papillomaviruses have been linked to several skin disorders in the dog. In order to have a suitable diagnostic tool for canine papillomavirus detection, eight PCRs with published primer combinations were evaluated. The most sensitive PCR was used to demonstrate that papillomavirus DNA can be detected on nonlesional skin of dogs.

The first complete papillomavirus genome characterized from a marsupial host: a novel isolate from Bettongia penicillata

The first fully sequenced papillomavirus (PV) of marsupials, tentatively named Bettongia penicillata papillomavirus type 1 (BpPV1), was detected in papillomas from a woylie (Bettongia penicillata ogilbyi). The circular, double-stranded DNA genome contains 7,737 bp and encodes 7 open reading frames (ORFs), E6, E7, E1, E2, E4, L2, and L1, in typical PV conformation. BpPV1 is a close-to-root PV with L1 and L2 ORFs most similar to European hedgehog PV and bandicoot papillomatosis carcinomatosis virus types 1 and 2 (BPCV1 and -2). It appears that the BPCVs arose by recombination between an ancient PV and an ancient polyomavirus more than 10 million years ago.

Analysis of host-parasite incongruence in papillomavirus evolution using importance sampling

The papillomaviruses (PVs) are a family of viruses infecting several mammalian and nonmammalian species that cause cervical cancer in humans. The evolutionary history of the PVs as it associated with a wide range of host species is not well understood. Incongruities between the phylogenetic trees of various viral genes as well as between these genes and the host phylogenies suggest historical viral recombination as well as violations of strict virus–host cospeciation. The extent of recombination events among PVs is uncertain, however, and there is little evidence to support a theory of PV spread via recent host transfers. We have investigated incongruence between PV genes and hence, the possibility of recombination, using Bayesian phylogenetic methods. We find significant evidence for phylogenetic incongruence among the six PV genes E1, E2, E6, E7, L1, and L2, indicating substantial recombination. Analysis of E1 and L1 phylogenies suggests ancestral recombination events. We also describe a new method for examining alternative host–parasite association mechanisms by applying importance sampling to Bayesian divergence time estimation. This new approach is not restricted by a fixed viral tree topology or knowledge of viral divergence times, multiple parasite taxa per host may be included, and it can distinguish between prior divergence of the virus before host speciation and host transfer of the virus following speciation. Using this method, we find prior divergence of PV lineages associated with the ancestral mammalian host resulting in at least 6 PV lineages prior to speciation of this host. These PV lineages have then followed paths of prior divergence and cospeciation to eventually become associated with the extant host species. Only one significant instance of host transfer is supported, the transfer of the ancestral L1 gene between a Primate and Hystricognathi host based on the divergence times between the υ human type 41 and porcupine PVs.

In situ hybridization to detect bandicoot papillomatosis carcinomatosis virus type 1 in biopsies from endangered western barred bandicoots (Perameles bougainville)

The western barred bandicoot (Perameles bougainville) is an endangered Australian marsupial species in which a papillomatosis and carcinomatosis syndrome occurs. Bandicoot papillomatosis carcinomatosis virus type 1 (BPCV1) is associated with the lesions of this progressively debilitating syndrome. Five digoxigenin-labelled DNA probes were generated for in situ hybridization (ISH) and the technique was optimized and performed on formalin-fixed paraffin-embedded (FFPE) biopsies. Staining of keratinocyte and sebocyte nuclei within lesions was achieved with all five probes. The sensitivity of ISH (76.9%) surpassed that of PCR (30.8%) for FFPE samples. The sensitivity of ISH varied from 81% (papillomas) and 70% (carcinoma in situ) to 29% (squamous cell carcinomas). The specificity of the test was confirmed using an irrelevant probe and papillomas from other species. These results strengthen the association between BPCV1 and the western barred bandicoot papillomatosis and carcinomatosis syndrome and give insight into the biology of the virus-host interaction.

Proliferative dermatitis in a loggerhead turtle, Caretta caretta, and a green turtle, Chelonia mydas, associated with novel papillomaviruses

A subadult loggerhead turtle, Caretta caretta, presented with generalized small, white, raised lesions over its neck, shoulders, and all four flippers. A juvenile green turtle, Chelonia mydas, recently treated for fibropapillomatosis, presented with four similar localized lesions on one flipper. To diagnose the conditions, biopsies of the lesions were taken for histopathology, electron microscopy, and molecular diagnostics. Histopathologic findings were similar in the two turtles and skin lesions were characterized by multifocal areas of epidermal hyperplasia accompanied by variation and abnormalities in the nuclear morphology of keratinocytes and a few intranuclear inclusions in some cells. Transmission electron microscopy revealed multiple epithelial cells with large intranuclear aggregates of virions consistent in morphology with papillomavirus. Papillomavirus was detected in samples from both turtles by polymerase chain reaction (PCR). Sequence analysis of the partial sequence of the papillomavirus E1 gene revealed two viruses (CcPV and CmPV) that were distinct from each other and from other species in Papillomaviridae, and likely represent two novel species and perhaps a new genus.

A novel virus detected in papillomas and carcinomas of the endangered western barred bandicoot (Perameles bougainville) exhibits genomic features of both the Papillomaviridae and Polyomaviridae

Conservation efforts to prevent the extinction of the endangered western barred bandicoot (Perameles bougainville) are currently hindered by a progressively debilitating cutaneous and mucocutaneous papillomatosis and carcinomatosis syndrome observed in captive and wild populations. In this study, we detected a novel virus, designated the bandicoot papillomatosis carcinomatosis virus type 1 (BPCV1), in lesional tissue from affected western barred bandicoots using multiply primed rolling-circle amplification and PCR with the cutaneotropic papillomavirus primer pairs FAP59/FAP64 and AR-L1F8/AR-L1R9. Sequencing of the BPCV1 genome revealed a novel prototype virus exhibiting genomic properties of both the Papillomaviridae and the Polyomaviridae. Papillomaviral properties included a large genome size ( approximately 7.3 kb) and the presence of open reading frames (ORFs) encoding canonical L1 and L2 structural proteins. The genomic organization in which structural and nonstructural proteins were encoded on different strands of the double-stranded genome and the presence of ORFs encoding the nonstructural proteins large T and small t antigens were, on the other hand, typical polyomaviral features. BPCV1 may represent the first member of a novel virus family, descended from a common ancestor of the papillomaviruses and polyomaviruses recognized today. Alternatively, it may represent the product of ancient recombination between members of these two virus families. The discovery of this virus could have implications for the current taxonomic classification of Papillomaviridae and Polyomaviridae and can provide further insight into the evolution of these ancient virus families.

High prevalence of bovine papillomaviral DNA in the normal skin of equine sarcoid-affected and healthy horses

Bovine papillomavirus (BPV), the causative agent of papillomas in cattle, has been shown to play a major role in the pathogenesis of equine sarcoids in horses. BPV has also been detected occasionally in normal equine skin. In this study, presence and activity of BPV in normal skin and peripheral blood of 4 groups of horses were evaluated: sarcoid-affected horses, horses living in contact with sarcoid-affected horses, horses living in contact with papilloma-affected cattle and control horses. From each horse, 3 samples on 4 locations were collected: a swab of the intact skin surface and both a swab and a biopsy after decontamination. BPV DNA was found in the normal skin of 24 of 42 horses (57%). Mainly sarcoid-affected horses and horses living in contact with cattle were carriers (73%), but BPV DNA was also detected in 50% of the horses living in contact with sarcoid-affected horses and in 30% of the control population. BPV mRNA was detected in 58% of the samples positive for BPV DNA, although in a much lower quantity compared to sarcoids. In most of the BPV DNA positive samples mild acanthosis, slight basophilic cytoplasmic swelling of the epidermal layers and/or thickening of the basal membrane were noticed, but these observations were also present in several BPV DNA negative normal skin samples. BPV DNA could not be detected in peripheral blood. These findings suggest latent infection and a wide-spread occurrence of BPV in the horse population.