Novel papillomavirus isolated from the oral mucosa of a polar bear does not cluster with other papillomaviruses of carnivores

Diverse papillomaviruses identified from Antarctic fur seals, leopard seals and Weddell seals from the Antarctic

Papillomaviruses (family Papillomaviridae) are non-enveloped, circular, double-stranded DNA viruses known to infect squamous and mucosal epithelial cells. In the family Papillomaviridae there are 53 genera and 133 viral species whose members infect a variety of mammalian, avian, reptilian, and fish species. Within the Antarctic context, papillomaviruses (PVs) have been identified in Adélie penguins (Pygoscelis adeliae, 2 PVs), Weddell seals (Leptonychotes weddellii, 7 PVs), and emerald notothen (Trematomus bernacchii, 1 PV) in McMurdo Sound and Ross Island in eastern Antarctica. Here we identified 13 diverse PVs from buccal swabs of Antarctic fur seals (Arctocephalus gazella, 2 PVs) and leopard seal (Hydrurga leptonyx, 3 PVs) in western Antarctica (Antarctic Peninsula), and vaginal and nasal swabs of Weddell seals (8 PVs) in McMurdo Sound. These PV genomes group into four genera representing 11 new papillomavirus types, of which five are from two Antarctic fur seals and a leopard seal and six from Weddell seals.

Canis Familiaris Papillomavirus Type 26: A Novel Papillomavirus of Dogs and the First Canine Papillomavirus within the Omegapapillomavirus Genus

Domestic dogs are currently recognized as being infected by 25 different canine papillomavirus (CPV) types classified into three genera. A short sequence from a novel CPV type was amplified, along with CPV1, from a papilloma (wart) from the mouth of a dog. The entire 7499 bp genome was amplified, and CPV26 contained putative coding regions that were predicted to produce four early proteins and two late ones. The ORF L1 showed less than 62% similarity for all previously sequenced CPV types but over 69% similarity to multiple Omegapapillomavirus types from a variety of Caniform species including the giant panda, Weddel seal, and polar bear. Phylogenetic analysis confirmed CPV26 clusters within the Omegapapillomavirus genus. Specific primers were used to investigate the presence of CPV26 DNA within a series of 37 canine proliferative lesions. CPV26 DNA was amplified from one lesion, a cutaneous papilloma that also contained CPV6. This is the first time a PV type within the Omegapapillomavirus genus has been detected in a non-domestic species and this provides evidence that the omegapapillomaviruses infected a common ancestor of, and then co-evolved with, the Caniform species. Whether CPV26 causes disease is uncertain, but the absence of an E7 protein may suggest low pathogenicity.

Genomic sequencing of a virus representing a novel type within the species Dyopipapillomavirus 1 in an Indian River Lagoon bottlenose dolphin

Fecal samples collected from free-ranging Atlantic bottlenose dolphins (BDs) in the Indian River Lagoon of Florida were processed for viral discovery using a next-generation sequencing (NGS) approach. A 693-bp contig identified in the NGS data was nearly identical to the partial L1 gene sequence of a papillomavirus (PV) previously found in a penile papilloma in a killer whale (Orcinus orca). Based on this partial bottlenose dolphin papillomavirus (BDPV) sequence, a nested inverse PCR and primer-walking strategy was employed to generate the complete genome sequence. The full BDPV genome consisted of 7299 bp and displayed a typical PV genome organization. The BDPV E6 protein contained a PDZ-binding motif, which has been shown to be involved in carcinogenic transformation involving high-risk genital human PVs. Screening of 12 individual fecal samples using a specific endpoint PCR assay revealed that the feces from a single female BD displaying a genital papilloma was positive for the BDPV. Genetic analysis indicated that this BDPV (Tursiops truncatus papillomavirus 8; TtPV8) is a new type of Dyopipapillomavirus 1, previously sequenced from an isolate obtained from a penile papilloma in a harbor porpoise (Phocoena phocoena). Although only a partial L1 sequence has been determined for a PV detected in a killer whale genital papilloma, our finding of a nearly identical sequence in an Atlantic BD may indicate that members of this viral species are capable of host jumping. Future work is needed to determine if this virus is a high-risk PV that is capable of inducing carcinogenic transformation and whether it poses a significant health risk to wild delphinid populations.

Complete Genome Sequences of Papillomavirus Isolates from the Oral Cavity, Skin, and Feces of Wild Rats

Six genome sequences of papillomavirus were determined from oral and skin swabs and fecal samples collected from wild rats. Three genomes were 7,722 bp, two genomes were 7,716 bp, and one was 7,730 bp, displaying typical papillomavirus genome organizations. Phylogenetic analysis revealed that these six genomes belonged to two different clusters.

Diverse papillomaviruses identified in Weddell seals

Papillomaviridae is a diverse family of circular, double-stranded DNA (dsDNA) viruses that infect a broad range of mammalian, avian and fish hosts. While papillomaviruses have been characterized most extensively in humans, the study of non-human papillomaviruses has contributed greatly to our understanding of their pathogenicity and evolution. Using high-throughput sequencing approaches, we identified 7 novel papillomaviruses from vaginal swabs collected from 81 adult female Weddell seals (Leptonychotes weddellii) in the Ross Sea of Antarctica between 2014-2017. These seven papillomavirus genomes were amplified from seven individual seals, and six of the seven genomes represented novel species with distinct evolutionary lineages. This highlights the diversity of papillomaviruses among the relatively small number of Weddell seal samples tested. Viruses associated with large vertebrates are poorly studied in Antarctica, and this study adds information about papillomaviruses associated with Weddell seals and contributes to our understanding of the evolutionary history of papillomaviruses.

Genomic characterization of the first oral avian papillomavirus in a colony of breeding canaries (Serinus canaria)

Papillomaviruses are non-enveloped, DNA viruses that infect skin and mucosa of a wide variety of vertebrates, causing neoplasias or simply persisting asymptomatically. Avian papillomaviruses, with six fully sequenced genomes, are the second most studied group after mammalian papillomaviruses. In this study, we describe the first oral avian papillomavirus, detected in the tongue of a dead Yorkshire canary (Serinus canaria) and in oral swabs of the same bird and other two live canaries from an aviary in Madrid, Spain. Its genome is 8,071 bp and presents the canonical papillomavirus architecture with six early (E6, E7, E1, E9, E2, E4) and two late open reading frames (L1 and L2) and a long control region between L1 and E6. This new avian papillomavirus L1 gene shares a 64% pairwise identity with FcPV1 L1, so it has been classified as a new species (ScPV1) within the Ethapapillomavirus genus. Although the canary died after showing breathing problems, there is no evidence that the papillomavirus caused those symptoms so it could be part of the oral microbiota of the birds. Hence, future investigations are needed to evaluate the clinical relevance of the virus.

Association of papillomavirus E6 proteins with either MAML1 or E6AP clusters E6 proteins by structure, function, and evolutionary relatedness

Papillomavirus E6 proteins bind to LXXLL peptide motifs displayed on targeted cellular proteins. Alpha genus HPV E6 proteins associate with the cellular ubiquitin ligase E6AP (UBE3A), by binding to an LXXLL peptide (ELTLQELLGEE) displayed by E6AP, thereby stimulating E6AP ubiquitin ligase activity. Beta, Gamma, and Delta genera E6 proteins bind a similar LXXLL peptide (WMSDLDDLLGS) on the cellular transcriptional co-activator MAML1 and thereby repress Notch signaling. We expressed 45 different animal and human E6 proteins from diverse papillomavirus genera to ascertain the overall preference of E6 proteins for E6AP or MAML1. E6 proteins from all HPV genera except Alpha preferentially interacted with MAML1 over E6AP. Among animal papillomaviruses, E6 proteins from certain ungulate (SsPV1 from pigs) and cetacean (porpoises and dolphins) hosts functionally resembled Alpha genus HPV by binding and targeting the degradation of E6AP. Beta genus HPV E6 proteins functionally clustered with Delta, Pi, Tau, Gamma, Chi, Mu, Lambda, Iota, Dyokappa, Rho, and Dyolambda E6 proteins to bind and repress MAML1. None of the tested E6 proteins physically and functionally interacted with both MAML1 and E6AP, indicating an evolutionary split. Further, interaction of an E6 protein was insufficient to activate degradation of E6AP, indicating that E6 proteins that target E6AP co-evolved to separately acquire both binding and triggering of ubiquitin ligase activation. E6 proteins with similar biological function clustered together in phylogenetic trees and shared structural features. This suggests that the divergence of E6 proteins from either MAML1 or E6AP binding preference is a major event in papillomavirus evolution.

Rodent Papillomaviruses

Preclinical infection model systems are extremely valuable tools to aid in our understanding of Human Papillomavirus (HPV) biology, disease progression, prevention, and treatments. In this context, rodent papillomaviruses and their respective infection models are useful tools but remain underutilized resources in the field of papillomavirus biology. Two rodent papillomaviruses, MnPV1, which infects the Mastomys species of multimammate rats, and MmuPV1, which infects laboratory mice, are currently the most studied rodent PVs. Both of these viruses cause malignancy in the skin and can provide attractive infection models to study the lesser understood cutaneous papillomaviruses that have been frequently associated with HPV-related skin cancers. Of these, MmuPV1 is the first reported rodent papillomavirus that can naturally infect the laboratory strain of mice. MmuPV1 is an attractive model virus to study papillomavirus pathogenesis because of the ubiquitous availability of lab mice and the fact that this mouse species is genetically modifiable. In this review, we have summarized the knowledge we have gained about PV biology from the study of rodent papillomaviruses and point out the remaining gaps that can provide new research opportunities.

Virome comparisons in wild-diseased and healthy captive giant pandas

BACKGROUND

The giant panda (Ailuropoda melanoleuca) is a vulnerable mammal herbivore living wild in central China. Viral infections have become a potential threat to the health of these endangered animals, but limited information related to these infections is available.

METHODS

Using a viral metagenomic approach, we surveyed viruses in the feces, nasopharyngeal secretions, blood, and different tissues from a wild giant panda that died from an unknown disease, a healthy wild giant panda, and 46 healthy captive animals.

RESULTS

The previously uncharacterized complete or near complete genomes of four viruses from three genera in Papillomaviridae family, six viruses in a proposed new Picornaviridae genus (Aimelvirus), two unclassified viruses related to posaviruses in Picornavirales order, 19 anelloviruses in four different clades of Anelloviridae family, four putative circoviruses, and 15 viruses belonging to the recently described Genomoviridae family were sequenced. Reflecting the diet of giant pandas, numerous insect virus sequences related to the families Iflaviridae, Dicistroviridae, Iridoviridae, Baculoviridae, Polydnaviridae, and subfamily Densovirinae and plant viruses sequences related to the families Tombusviridae, Partitiviridae, Secoviridae, Geminiviridae, Luteoviridae, Virgaviridae, and Rhabdoviridae; genus Umbravirus, Alphaflexiviridae, and Phycodnaviridae were also detected in fecal samples. A small number of insect virus sequences were also detected in the nasopharyngeal secretions of healthy giant pandas and lung tissues from the dead wild giant panda. Although the viral families present in the sick giant panda were also detected in the healthy ones, a higher proportion of papillomaviruses, picornaviruses, and anelloviruses reads were detected in the diseased panda.

CONCLUSION

This viral survey increases our understanding of eukaryotic viruses in giant pandas and provides a baseline for comparison to viruses detected in future infectious disease outbreaks. The similar viral families detected in sick and healthy giant pandas indicate that these viruses result in commensal infections in most immuno-competent animals.

Assessing parallel gene histories in viral genomes

BACKGROUND

The increasing abundance of sequence data has exacerbated a long known problem: gene trees and species trees for the same terminal taxa are often incongruent. Indeed, genes within a genome have not all followed the same evolutionary path due to events such as incomplete lineage sorting, horizontal gene transfer, gene duplication and deletion, or recombination. Considering conflicts between gene trees as an obstacle, numerous methods have been developed to deal with these incongruences and to reconstruct consensus evolutionary histories of species despite the heterogeneity in the history of their genes. However, inconsistencies can also be seen as a source of information about the specific evolutionary processes that have shaped genomes.

RESULTS

The goal of the approach here proposed is to exploit this conflicting information: we have compiled eleven variables describing phylogenetic relationships and evolutionary pressures and submitted them to dimensionality reduction techniques to identify genes with similar evolutionary histories. To illustrate the applicability of the method, we have chosen two viral datasets, namely papillomaviruses and Turnip mosaic virus (TuMV) isolates, largely dissimilar in genome, evolutionary distance and biology. Our method pinpoints viral genes with common evolutionary patterns. In the case of papillomaviruses, gene clusters match well our knowledge on viral biology and life cycle, illustrating the potential of our approach. For the less known TuMV, our results trigger new hypotheses about viral evolution and gene interaction.

CONCLUSIONS

The approach here presented allows turning phylogenetic inconsistencies into evolutionary information, detecting gene assemblies with similar histories, and could be a powerful tool for comparative pathogenomics.

Characterization of a novel polyomavirus isolated from a fibroma on the trunk of an African elephant (Loxodonta africana)

Viruses of the family Polyomaviridae infect a wide variety of avian and mammalian hosts with a broad spectrum of outcomes including asymptomatic infection, acute systemic disease, and tumor induction. In this study a novel polyomavirus, the African elephant polyomavirus 1 (AelPyV-1) found in a protruding hyperplastic fibrous lesion on the trunk of an African elephant (Loxodonta africana) was characterized. The AelPyV-1 genome is 5722 bp in size and is one of the largest polyomaviruses characterized to date. Analysis of the AelPyV-1 genome reveals five putative open-reading frames coding for the classic small and large T antigens in the early region, and the VP1, VP2 and VP3 capsid proteins in the late region. In the area preceding the VP2 start codon three putative open-reading frames, possibly coding for an agnoprotein, could be localized. A regulatory, non-coding region separates the 2 coding regions. Unique for polyomaviruses is the presence of a second 854 bp long non-coding region between the end of the early region and the end of the late region. Based on maximum likelihood phylogenetic analyses of the large T antigen of the AelPyV-1 and 61 other polyomavirus sequences, AelPyV-1 clusters within a heterogeneous group of polyomaviruses that have been isolated from bats, new world primates and rodents.

The papillomavirus E7 proteins

E7 is an accessory protein that is not encoded by all papillomaviruses. The E7 amino terminus contains two regions of similarity to conserved regions 1 and 2 of the adenovirus E1A protein, which are also conserved in the simian vacuolating virus 40 large tumor antigen. The E7 carboxyl terminus consists of a zinc-binding motif, which is related to similar motifs in E6 proteins. E7 proteins play a central role in the human papillomavirus life cycle, reprogramming the cellular environment to be conducive to viral replication. E7 proteins encoded by the cancer-associated alpha human papillomaviruses have potent transforming activities, which together with E6, are necessary but not sufficient to render their host squamous epithelial cell tumorigenic. This article strives to provide a comprehensive summary of the published research studies on human papillomavirus E7 proteins.

Animal papillomaviruses

We provide an overview of the host range, taxonomic classification and genomic diversity of animal papillomaviruses. The complete genomes of 112 non-human papillomavirus types, recovered from 54 different host species, are currently available in GenBank. The recent characterizations of reptilian papillomaviruses extend the host range of the Papillomaviridae to include all amniotes. Although the genetically diverse papillomaviruses have a highly conserved genomic lay-out, deviations from this prototypic genome organization are observed in several animal papillomaviruses, and only the core ORFs E1, E2, L2 and L1 are present in all characterized papillomavirus genomes. The discovery of papilloma-polyoma hybrids BPCV1 and BPCV2, containing a papillomaviral late region but an early region encoding typical polyomaviral nonstructural proteins, and the detection of recombination breakpoints between the early and late coding regions of cetacean papillomaviruses, could indicate that early and late gene cassettes of papillomaviruses are relatively independent entities that can be interchanged by recombination.

Papillomavirus E6 oncoproteins

Papillomaviruses induce benign and malignant epithelial tumors, and the viral E6 oncoprotein is essential for full transformation. E6 contributes to transformation by associating with cellular proteins, docking on specific acidic LXXLL peptide motifs found on these proteins. This review examines insights from recent studies of human and animal E6 proteins that determine the three-dimensional structure of E6 when bound to acidic LXXLL peptides. The structure of E6 is related to recent advances in the purification and identification of E6 associated protein complexes. These E6 protein-complexes, together with other proteins that bind to E6, alter a broad array of biological outcomes including modulation of cell survival, cellular transcription, host cell differentiation, growth factor dependence, DNA damage responses, and cell cycle progression.

Multiple evolutionary origins of bat papillomaviruses

Infection by papillomaviruses (PVs) has been linked to different types of neoplasias, in both human and non-human hosts. Knowledge about PV diversity is essential to reliably infer the evolutionary history of these pathogens and to elucidate the link between infection and disease. We cloned and sequenced the complete genome of a novel PV, EhelPV1, isolated from hair bulbs from a captive straw-colored fruit bat Eidolon helvum (Pteropodidae, Chiroptera). We also retrieved partial sequences of the E1 and L1 genes from hair bulbs from a captive Indian flying fox Pteropus giganteus (Pteropodidae, Chiroptera). The detected virus (PgigPV1) presumably corresponded to a novel type as well. Maximum likelihood phylogenetic analyses were conducted using a representative collection of 132 PVs. EhelPV1 belonged to the Lambda+Mu-PV crown group and was most closely related to another bat PV, MschPV2. Both fragments of PgigPV1 were placed alongside with EhelPV1. The novel PVs were phylogenetically distant from other previously described bat PVs, namely MrPV1, MschPV1 and RaPV1. We have further characterized the sequence patterns of the E2-binding sites occurring in the upstream regulatory region of Lambda+Mu-PVs. Common fingerprints within this region are shared by certain PVs. However, there is not a sharp correspondence between the repertoire of transcription factor binding sites in the viral regulatory region and host range, tissue tropism or viral life style. Our results reinforce the hypothesis that PVs have undergone an initial radiation prior to the divergence of the mammalian hosts, giving rise to the present-day PV crown groups.

Virome analysis for identification of novel mammalian viruses in bat species from Chinese provinces

Bats are natural hosts for a large variety of zoonotic viruses. This study aimed to describe the range of bat viromes, including viruses from mammals, insects, fungi, plants, and phages, in 11 insectivorous bat species (216 bats in total) common in six provinces of China. To analyze viromes, we used sequence-independent PCR amplification and next-generation sequencing technology (Solexa Genome Analyzer II; Illumina). The viromes were identified by sequence similarity comparisons to known viruses. The mammalian viruses included those of the Adenoviridae, Herpesviridae, Papillomaviridae, Retroviridae, Circoviridae, Rhabdoviridae, Astroviridae, Flaviridae, Coronaviridae, Picornaviridae, and Parvovirinae; insect viruses included those of the Baculoviridae, Iflaviridae, Dicistroviridae, Tetraviridae, and Densovirinae; fungal viruses included those of the Chrysoviridae, Hypoviridae, Partitiviridae, and Totiviridae; and phages included those of the Caudovirales, Inoviridae, and Microviridae and unclassified phages. In addition to the viruses and phages associated with the insects, plants, and bacterial flora related to the diet and habitation of bats, we identified the complete or partial genome sequences of 13 novel mammalian viruses. These included herpesviruses, papillomaviruses, a circovirus, a bocavirus, picornaviruses, a pestivirus, and a foamy virus. Pairwise alignments and phylogenetic analyses indicated that these novel viruses showed little genetic similarity with previously reported viruses. This study also revealed a high prevalence and diversity of bat astroviruses and coronaviruses in some provinces. These findings have expanded our understanding of the viromes of bats in China and hinted at the presence of a large variety of unknown mammalian viruses in many common bat species of mainland China.

Novel snake papillomavirus does not cluster with other non-mammalian papillomaviruses

Papillomaviruses (PVs) are associated with the development of neoplasias and have been found in several different species, most of them in humans and other mammals. We identified, cloned and sequenced PV DNA from pigmented papilloma-like lesions of a diamond python (Morelia spilota spilota). This represents the first complete PV genome discovered in a Squamata host (MsPV1). It consists of 7048 nt and contains the characteristic open reading (ORF) frames E6, E7, E1, E2, L1 and L2. The L1 ORF sequence showed the highest percentage of sequence identities to human PV5 (57.9%) and Caribbean manatee (Trichechus manatus) PV1 (55.4%), thus, establishing a new clade. According to phylogenetic analysis, the MsPV1 genome clusters with PVs of mammalian rather than sauropsid hosts.

Discovery of a new human polyomavirus associated with trichodysplasia spinulosa in an immunocompromized patient

The Polyomaviridae constitute a family of small DNA viruses infecting a variety of hosts. In humans, polyomaviruses can cause infections of the central nervous system, urinary tract, skin, and possibly the respiratory tract. Here we report the identification of a new human polyomavirus in plucked facial spines of a heart transplant patient with trichodysplasia spinulosa, a rare skin disease exclusively seen in immunocompromized patients. The trichodysplasia spinulosa-associated polyomavirus (TSV) genome was amplified through rolling-circle amplification and consists of a 5232-nucleotide circular DNA organized similarly to known polyomaviruses. Two putative “early” (small and large T antigen) and three putative “late” (VP1, VP2, VP3) genes were identified. The TSV large T antigen contains several domains (e.g. J-domain) and motifs (e.g. HPDKGG, pRb family-binding, zinc finger) described for other polyomaviruses and potentially involved in cellular transformation. Phylogenetic analysis revealed a close relationship of TSV with the Bornean orangutan polyomavirus and, more distantly, the Merkel cell polyomavirus that is found integrated in Merkel cell carcinomas of the skin. The presence of TSV in the affected patient's skin was confirmed by newly designed quantitative TSV-specific PCR, indicative of a viral load of 10⁵ copies per cell. After topical cidofovir treatment, the lesions largely resolved coinciding with a reduction in TSV load. PCR screening demonstrated a 4% prevalence of TSV in an unrelated group of immunosuppressed transplant recipients without apparent disease. In conclusion, a new human polyomavirus was discovered and identified as the possible cause of trichodysplasia spinulosa in immunocompromized patients. The presence of TSV also in clinically unaffected individuals suggests frequent virus transmission causing subclinical, probably latent infections. Further studies have to reveal the impact of TSV infection in relation to other populations and diseases.

Diseases that occur exclusively in immunocompromized patients are often of an infectious nature. Trichodysplasia spinulosa (TS) is such a disease characterized by development of papules, spines and alopecia in the face. Fortunately this disease is rare, because facial features can change dramatically, as in the case of an adolescent TS patient who was on immunosuppressive drugs because of heart-transplantation. A viral cause of TS was suspected already for some time because virus particles had been seen in TS lesions. In pursuit of this unknown virus, we isolated DNA from collected TS spines and could detect a unique small circular DNA suggestive of a polyomavirus genome. Additional experiments confirmed the presence in these samples of a new polyomavirus that we tentatively called TS-associated polyomavirus (TSPyV or TSV). TSV shares several properties with other polyomaviruses, such as genome organization and proteome composition, association with disease in immunosuppressed patients and occurence in individuals without overt disease. The latter indicates that TSV circulates in the human population. Future studies have to show how this newly identified polyomavirus spreads, how it causes disease and if it is related to other (skin) conditions as well.

Classification of papillomaviruses (PVs) based on 189 PV types and proposal of taxonomic amendments

We present an expansion of the classification of the family Papillomaviridae, which now contains 29 genera formed by 189 papillomavirus (PV) types isolated from humans (120 types), non-human mammals, birds and reptiles (64, 3 and 2 types, respectively). To accommodate the number of PV genera exceeding the Greek alphabet, the prefix "dyo" is used, continuing after the Omega-PVs with Dyodelta-PVs. The current set of human PVs is contained within five genera, whereas mammalian, avian and reptile PVs are contained within 20, 3 and 1 genera, respectively. We propose standardizations to the names of a number of animal PVs. As prerequisite for a coherent nomenclature of animal PVs, we propose founding a reference center for animal PVs. We discuss that based on emerging species concepts derived from genome sequences, PV types could be promoted to the taxonomic level of species, but we do not recommend implementing this change at the current time.

Isolation and cloning of two variant papillomaviruses from domestic pigs: Sus scrofa papillomaviruses type 1 variants a and b

The healthy skin of two female domestic pigs (Sus scrofa domestica) was sampled with cotton-tipped swabs. Total genomic DNA was extracted from the samples and subjected to PCR with degenerate papillomavirus (PV)-specific primers. Similarity searches performed with blastn showed that partial E1 and L1 sequences of two novel PVs were amplified. Subsequently, the complete genomes of these Sus scrofa papillomaviruses (SsPVs) were amplified by long-template PCR, cloned and sequenced using a transposon insertion method. They contained the typical PV open reading frames (ORFs) E1, E2, E4, E6, L1 and L2, but the E7 ORF was absent in both viruses. Pairwise nucleotide sequence alignment of the L1 ORFs of the SsPVs showed 98.5 % similarity, classifying these viruses as SsPV type 1 'variants' (SsPV-1a and -1b). Based on a concatenated alignment of the E1, E2, L1 and L2 ORFs of SsPV-1 variants a and b, and 81 other human and animal PV type species, a neighbour-joining phylogenetic tree was constructed. This phylogenetic analysis showed that the SsPV-1a and -1b variants did not cluster with the other PVs of artiodactyls (cloven-hoofed) host species, but clustered on the edge of the genus Alphapapillomavirus, very near to the root of this genus.