Characterization of a novel close-to-root papillomavirus from a Florida manatee by using multiply primed rolling-circle amplification: Trichechus manatus latirostris papillomavirus type 1

Genetic characterization of Sus scrofa papillomavirus type 1 from domestic pigs in Guangxi Province, China

Sus scrofa papillomatosis (SsP) is a tumour caused by Sus scrofa papillomaviruses (SsPVs). To investigate the presence of SsPVs in China, 354 domestic pig skin samples collected from Guangxi Province were examined for SsPV DNA by PCR. Three SsPV1s (GX12, GX14, and GX18) were identified with a prevalence of 0.847% (3/354). Sequence analysis showed that L1 of SsPV1/GX12 and SsPV1/GX14 had 99.7% and 99.6% nucleotide identify with the reference SsPV1a, respectively. Phylogenetic and evolutionary analyses showed that SsPV1/GX12 and SsPV1/14 clustered into SsPV1a and that SsPV1/GX18 clustered into SsPV1b. Compared with other SsPV L1 and L2 proteins, we found that the SsPV1/GX18 and SsPV1b strains shared the same unique substitutions, and SsPV1/GX12, SsPV1/GX14, and SsPV1a shared almost identical amino acid sequences. This study reports the first detection of SsPV DNA in China based on whole genome information and provides a scientific basis for the development of SsPV pathogenic biology, epidemiology, and prevention, as well as control technology research.

Co-infection of distinct papillomavirus types in a captive North American porcupine

BACKGROUND

Only two cases of papillomavirus infections in North American porcupines (Erethizon dorsatum) have been described thus far, and molecular investigation linked these cases to two distinct papillomavirus species.

METHODS

In this report, we present the clinical, histological and molecular investigation of a third case of a porcupine papillomavirus infection. Papillomatous lesions occurred on the upper and lower lip of an otherwise healthy three-year old female that was kept in captivity. Within one month, the lesions progressed into exophytic black nodules, followed by a temporary stabilization and ultimately spontaneous regression within seven months of their initial observation. PCR-based screening using specific primers for Erethizon dorsatum papillomavirus 1 and 2 revealed the presence of both these virus types, after which nanopore sequencing was used to determine the complete sequences of the two virus genomes.

RESULTS

One of the genomes shares 99.9% similarity with the only known sequence for Erethizon dorsatum papillomavirus 1, while the second represents a distinct lineage of Erethizon dorsatum papillomavirus 2, sharing only 93.3% similarity with the previously discovered strain.

CONCLUSIONS

This report marks the first observation of a papillomavirus co-infection in a North American porcupine, although the individual contribution of the two virus types to the clinical presentation was not assessed.

Characterization of novel papillomavirus from free-ranging Antillean manatee Trichechus manatus manatus with genital papillomatosis

The Antillean manatee Trichechus manatus manatus is an Endangered species living along the Atlantic coasts of the Americas from Florida (USA), throughout the Caribbean, to Brazil. In July 2020, a manatee with multiple wounds due to boat-inflicted trauma was rescued from the coast east of Cayo Mata, Salinas, Puerto Rico. This manatee had neutropenia, leukopenia, and monocytosis associated with immunosuppression and nutritional deficiency anemia, as well as bacteria and fungi within the lesions. The manatee had genital lesions which included papules and linear plaques, microscopically characterized by mucosal hyperplasia with cytopathic changes typical of papillomavirus infection. Superficial epithelial cells had strong nuclear immunolabeling when examined using a monoclonal antibody specific to papillomavirus. The sequencing data of PCR products with papillomavirus-specific degenerative primers indicated that these lesions contained a novel manatee papillomavirus (Trichechus manatus papillomavirus, TmPV). The genomic DNA was amplified using a rolling circle amplification, and fully sequenced to be 7586 bp (GenBank accession no. OK073977). Other TmPVs were previously isolated from Florida manatees T. manatus latirostris. This novel virus was designated TmPV type 5 (TmPV5) based on its genomic characterization and sequence comparison. The TmPV5 genome shared 50.7, 48.9, 69.4, and 62.1% similarities with TmPV1, TmPV2, TmPV3, and TmPV4, respectively. TmPV5 is classified in the genus Rhopapillomavirus together with other manatee papillomaviruses. After 2.5 mo of veterinary treatment and rehabilitation, the manatee recovered and was released. This is the first report of papillomatosis in a free-ranging Antillean manatee.

Detection of Prorocentrum minimum by hyperbranched rolling circle amplification coupled with lateral flow dipstick

A novel method referred to as hyperbranched rolling circle amplification (HRCA) coupled with lateral flow dipstick (LFD) (HRCA-LFD) here was developed for specific, sensitive, rapid, and simple detection of Prorocentrum minimum. HRCA-LFD relies on a padlock probe (PLP) consisting of a common ligation sequence, two terminal sequences that complement the target DNA, and a manually designed detection probe (LFD probe). The two terminal sequences of the PLP were designed against the species-specific sites of the large subunit ribosomal DNA (LSU rDNA) D1-D2 region of P. minimum. The optimum parameters for HRCA were as follows: PLP concentration of 20 pM, ligation time of 30 min, ligation temperature of 59 °C, enzymic digestion time of 105 min, amplification time of 45 min, and amplification temperature of 58 °C. The HRCA-LFD displaying high specificity could accurately distinguish P. minimum from other microalgae. The detection limit of HRCA-LFD was as low as 1.42 × 10^-7 ng μL^-1 for genomic DNA, 1.03 × 10^-7 ng μL^-1 (approximately 27 copies) for recombinant plasmid containing the inserted LSU rDNA D1-D2, and 0.17 cells for crude DNA extract of P. minimum, which was consistently 100 times more sensitive than regular PCR. Interfering test suggested that the performance of HRCA-LFD is stable and would not be affected by other non-target species. The HRCA-LFD results of field samples that are comparable with microscopic examination confirmed that the developed method is competent for detection of target cells in field samples. In conclusion, the developed HRCA-LFD exhibiting stable performance is specific, sensitive, and rapid, which provides a good alternative to traditional microscopic examination for the detection of P. minimum cells in field samples.

Bovine Papillomavirus 1 Gets Out of the Flock: Detection in an Ovine Wart in Sicily

A proliferative cauliflower lesion was excised from the udder of a sheep. Histological investigation confirmed the macroscopic classification of the lesion as a papilloma, without any fibroblastic proliferation. PCR revealed the presence of bovine papillomavirus (BPV), which was further confirmed by the identification of a Deltapapillomavirus4 by Next Generation Sequencing analysis. This was subsequently classified as bovine papillomavirus type 1. Negative staining electron microscopy (EM) analyses produced negative test results for papillomavirus particles. RNA in situ hybridization (ISH) confirmed the presence of BPV-1. The results further confirm the ability of BPVs belonging to the Deltapapillomavirus genus to infect distantly related species and to cause lesions that are different from sarcoids.

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.

Genomic characterization of Erethizon dorsatum papillomavirus 2, a new papillomavirus species marked by its exceptional genome size

We report here the complete sequence and genome organization of a new papillomavirus, Erethizon dorsatum papillomavirus 2 (EdPV2), which was isolated from cutaneous lesions observed on the muzzle of a North American porcupine. The complete genome is 8809 nucleotides long and encodes five early (E6-E7-E1-E2-E4) and two late proteins (L2-L1). In addition to the upstream regulatory region, the EdPV2 genome contains an exceptionally large secondary non-coding region with no apparent functional relevance. EdPV2 is strongly divergent from the previously described porcupine papillomavirus EdPV1 and phylogenetic analysis shows EdPV2 clustering near members of the genus Pipapillomavirus, a group of rodent papillomaviruses. Pairwise sequence comparison based on the L1 open reading frame identifies Rattus norvegicus papillomavirus 1 as the closest related virus (59.97 % similarity). Based on its low sequence similarity to other known papillomaviruses, EdPV2 is thought to represent a new genus in the family Papillomaviridae.

Papillomavirus Infection in Humans and Dromedary Camels in Eastern Sudan

Cases of wart-like lesions in humans and dromedary camels occurred in eastern Sudan in 2015 were described. Involvement of papillomavirus (PV) in causing these cases was affirmed by PCR and immunoperoxidase test. Mostly, the lesions were observed on the skin of the chest and forearms in addition to lips and mandible. Sequence analysis revealed Camelus dromedarius PV types 1 and 2 genotypes as the causative genotypes. We also observed cases of wart-like lesions on hands and legs of two herders attending the infected camel herd. Partial genome sequencing revealed human PV type 2 in one of the two human samples providing no indications for interspecies transmission of camel PVs, yet provides, for the first time evidence of active circulation of this virus in eastern Sudan.

Identification of G-quadruplex forming sequences in three manatee papillomaviruses

The Florida manatee (Trichechus manatus latirotris) is a threatened aquatic mammal in United States coastal waters. Over the past decade, the appearance of papillomavirus-induced lesions and viral papillomatosis in manatees has been a concern for those involved in the management and rehabilitation of this species. To date, three manatee papillomaviruses (TmPVs) have been identified in Florida manatees, one forming cutaneous lesions (TmPV1) and two forming genital lesions (TmPV3 and TmPV4). We identified DNA sequences with the potential to form G-quadruplex structures (G4) across the three genomes. G4 were located on both DNA strands and across coding and non-coding regions on all TmPVs, offering multiple targets for viral control. Although G4 have been identified in several viral genomes, including human PVs, most research has focused on canonical structures comprised of three G-tetrads. In contrast, the vast majority of sequences we identified would allow the formation of non-canonical structures with only two G-tetrads. Our biophysical analysis confirmed the formation of G4 with parallel topology in three such sequences from the E2 region. Two of the structures appear comprised of multiple stacked two G-tetrad structures, perhaps serving to increase structural stability. Computational analysis demonstrated enrichment of G4 sequences on all TmPVs on the reverse strand in the E2/E4 region and on both strands in the L2 region. Several G4 sequences occurred at similar regional locations on all PVs, most notably on the reverse strand in the E2 region. In other cases, G4 were identified at similar regional locations only on PVs forming genital lesions. On all TmPVs, G4 sequences were located in the non-coding region near putative E2 binding sites. Together, these findings suggest that G4 are possible regulatory elements in TmPVs.

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.

Papillomaviruses are a large family of viruses with great medical and veterinary importance. This study explores the viral E6 oncoproteins from diverse papillomavirus genera to determine how E6 distinguishes in interaction between cellular proteins. E6 proteins have been previously found to interact with a ubiquitin ligase called E6AP and thereby target particular cellular proteins for degradation, or to interact with MAML family proteins to repress Notch signaling and thereby alter cellular differentiation. It has been unclear if diverse families of papillomavirus E6 proteins interact with only E6AP or MAML (or possibly both), how E6 distinguishes between these interactions, and if interaction of E6 with E6AP is coupled to ubiquitin ligase activation. We find here that none of the tested E6 proteins physically and functionally interacted with both E6AP and MAML1, indicating an evolutionary split that clustered E6 proteins by sequence similarity analysis. Currently, the categorization of papillomaviruses is complex, with thirty-eight genera so far described. This study establishes an early evolutionary split among most papillomavirus genera between those viruses that encode E6 proteins that physically and functionally associate with MAML compared to E6AP. This provides a structural and functional basis for categorizing most currently described papillomaviruses into two major functional groups.

Viral metagenomics of six bat species in close contact with humans in southern China

Accumulating studies have shown that bats could harbor various important pathogenic viruses that could be transmitted to humans and other animals. Extensive metagenomic studies of different organs/tissues from bats have revealed a large number of novel or divergent viruses. To elucidate viral diversity and epidemiological and phylogenetic characteristics, six pooled fecal samples from bats were generated (based on bat species and geographic regions characteristic for virome analysis). These contained 500 fecal samples from six bat species, collected in four geographic regions. Metagenomic analysis revealed a plethora of divergent viruses originally found in bats. Multiple contigs from influenza A virus and coronaviruses in bats shared high identity with those from humans, suggesting possible cross-species transmission, whereas a number of contigs, whose sequences were taxonomically classifiable within Alphapapillomavirus, Betaretrovirus, Alpharetrovirus, Varicellovirus, Cyprinivirus, Chlorovirus and Cucumovirus had low identity to viruses in existing databases, which indicated possible evolution of novel viral species. None of the established caliciviruses and picornaviruses were found in the 500 fecal specimens. Papillomaviruses with high amino acid identity were found in Scotophilus kuhlii and Rhinolophus blythi, challenging the hypotheses regarding the strict host specificity and co-evolution of papillomaviruses. Phylogenetic analysis showed that four bat rotavirus A strains might be tentative G3 strains, according to the Rotavirus Classification Working Group classification.

Investigation on papillomavirus infection in dromedary camels in Al-Ahsa, Saudi Arabia

We investigated two outbreaks of papillomatosis between 2013 and 2015 in Al Ahsa region of eastern Saudi Arabia involving fourteen dromedary camels. The disease affected both young and adult animals and occurred in coincidence with demodectic mange infestation. Diagnosis was made based on gross and histopathological characteristics of the wart lesion and was confirmed by PCR. Rolling circle amplification followed by degenerate primer PCR and sequencing of the amplicons revealed the presence of both Camelus dromedarius papillomavirus types 1 and 2, previously identified in infected dromedaries in Sudan.

How many papillomavirus species can go undetected in papilloma lesions?

A co-infection comprising to at least seven papillomavirus (PV) types was detected by next generation sequencing (NGS) of randomly primed rolling circle amplification (RCA) products of a bovine (Bos taurus) papilloma lesion from the Brazilian Amazon region. Six putative new PV types that could not be detected by commonly used PCR protocols were identified. Their overall L1 nucleotide identities were less than 90% compared to described PV species and types. L1 nucleotide BLAST sequence hits showed that each new type was related to Beta, Gamma, Dyokappa, Dyoeta, and Xipapillomavirus, as well as two likely new unclassified genera. Our results show that the employment of NGS is relevant to the detection and characterization of distantly related PV and is of major importance in co-infection studies. This knowledge will help us understand the biology and pathogenesis of PV, as well as contribute to disease control. Moreover, we can also conclude that there are many unknown circulating PVs.

A single bat species in Cameroon harbors multiple highly divergent papillomaviruses in stool identified by metagenomics analysis

A number of PVs have been described in bats but to the best of our knowledge not from feces. Using a previously described NetoVIR protocol, Eidolon helvum pooled fecal samples (Eh) were treated and sequenced by Illumina next generation sequencing technology. Two complete genomes of novel PVs (EhPV2 and EhPV3) and 3 partial sequences (BATPV61, BATPV890a and BATPV890b) were obtained and analysis showed that the EhPV2 and EhPV3 major capsid proteins cluster with and share 60-64% nucleotide identity with that of Rousettus aegyptiacus PV1, thus representing new species of PVs within the genus Psipapillomavirus. The other PVs clustered in different branches of our phylogenetic tree and may potentially represent novel species and/or genera. This points to the vast diversity of PVs in bats and in Eidolon helvum bats in particular, therefore adding support to the current concept that PV evolution is more complex than merely strict PV-host co-evolution.

Molecular characterization of novel mucosotropic papillomaviruses from a Florida manatee (Trichechus manatus latirostris)

We isolated two new manatee papillomavirus (PV) types, TmPV3 and TmPV4, from a Florida manatee (Trichechus manatus latirostris). Two PV types were previously isolated from this species. TmPV1 is widely dispersed amongst manatees and a close-to-root PV; not much is known about TmPV2. The genomes of TmPV3 and TmPV4 were 7622 and 7771 bp in size, respectively. Both PVs had a genomic organization characteristic of all PVs, with one non-coding region and seven ORFs, including the E7 ORF that is absent in other cetacean PVs. Although these PVs were isolated from separate genital lesions of the same manatee, an enlarged E2/E4 ORF was found only in the TmPV4 genome. The full genome and L1 sequence similarities between TmPV3 and TmPV4 were 63.2 and 70.3 %, respectively. These genomes shared only 49.1 and 50.2 % similarity with TmPV1. The pairwise alignment of L1 nucleotide sequences indicated that the two new PVs nested in a monophyletic group of the genus Rhopapillomavirus, together with the cutaneotropic TmPV1 and TmPV2.

Bovine papillomavirus type 7 in Italy: complete genomes and sequence variants

Two novel bovine papillomavirus type 7 (BPV-7) variants have been identified in teat cutaneous papillomas affecting dairy cows in northern Italy. The entire genome sequences of two BPV-7 Italian variants showed major sequence differences in the long control region (LCR) and in the L2 gene compared to the Japanese reference strain. In order to define the stability of these genetic variants, the L2 and LCR sequences of seven further BPV-7 positive isolates were characterized. An insertion of six amino acids in the L2 structural protein has been detected in all samples while different genetic variants have been identified for the LCR. These findings provide new insights on intra-type variability of BPVs and represent a starting point for future studies aimed at establishing the biological role of the different BPV genomic regions and investigating the pathogenic potential of papillomavirus variants.

Evidence of zoonotic Poxviridae coinfections in clinically diagnosed papillomas using a newly developed mini-array test

Our study describes a newly developed mini-array test for the rapid detection of poxviruses in animals and humans. The method is based on detection that combines target nucleic acid amplification by polymerase chain reaction and specific hybridization, using enzyme-linked antibodies, allowing identification of zoonotic orthopoxviruses and parapoxviruses in animal and human biological samples. With 100% specificity, the test rules out the possibility of cross-reactions with viral agents causing look-alike diseases. The assay was employed in the field to investigate the causes of several outbreaks of a malignant proliferative skin disease that affected domestic ruminants in Sicily during 2011-2014. Due to specific aspects of the lesions, the animals were clinically diagnosed with papillomatosis. The mini-array test allowed the identification of coinfections caused by more than 1 viral species belonging to the Parapoxvirus and Orthopoxvirus genera, either in goats or in cattle. Our study suggests that the so-called "papillomatosis" can be the result of multiple infections with epitheliotropic viruses, including zoonotic poxviruses that cannot be properly identified with classical diagnostic techniques.

Emergence of infectious malignant thrombocytopenia in Japanese macaques (Macaca fuscata) by SRV-4 after transmission to a novel host

We discovered a lethal hemorrhagic syndrome arising from severe thrombocytopenia in Japanese macaques kept at the Primate Research Institute, Kyoto University. Extensive investigation identified that simian retrovirus type 4 (SRV-4) was the causative agent of the disease. SRV-4 had previously been isolated only from cynomolgus macaques in which it is usually asymptomatic. We consider that the SRV-4 crossed the so-called species barrier between cynomolgus and Japanese macaques, leading to extremely severe acute symptoms in the latter. Infectious agents that cross the species barrier occasionally amplify in virulence, which is not observed in the original hosts. In such cases, the new hosts are usually distantly related to the original hosts. However, Japanese macaques are closely related to cynomolgus macaques, and can even hybridize when given the opportunity. This lethal outbreak of a novel pathogen in Japanese macaques highlights the need to modify our expectations about virulence with regards crossing species barriers.

Oral papillomatosis caused by Enhydra lutris papillomavirus 1 (ElPV-1) in southern sea otters (Enhydra lutris nereis) in California, USA

The southern sea otter (Enhydra lutris nereis) is a threatened marine sentinel. During postmortem investigations of stranded sea otters from 2004 to 2013 in California, US, papillomas were detected in the oral cavity of at least seven otters via necropsy and histopathology. Next-generation sequencing of viral particles purified from a single papilloma revealed a novel papillomavirus, Enhydra lutris papillomavirus 1 (ElPV-1). The genome of ElPV-1 was obtained, representing the first fully sequenced viral genome from southern sea otters. Phylogenetic analysis of the entire L1 gene, as well as a concatenated protein identities plot of all papillomaviral genes revealed that ElPV-1 is a λ-papillomavirus, related to a raccoon papillomavirus (Procyon lotor papillomavirus type 1) and a canine oral papillomavirus. Immunohistochemical staining, using a cross-reactive bovine papillomavirus antibody, suggested that ElPV-1 is present in intranuclear inclusions and intracytoplasmic keratin granules. Virus-infected cells were scattered throughout the stratum granulosum and stratum spinosum of the gingival and buccal papillomas. Using ElPV-1-specific PCR, we confirmed viral DNA in oral papillomas from all seven stranded sea otters, with identical L1 sequences. This virus is associated with the development of oral papillomatosis in southern sea otters.

Development of rolling circle amplification based surface-enhanced Raman spectroscopy method for 35S promoter gene detection

In this study, we developed the genetically modified organism detection method by using the combination of rolling circle amplification (RCA) and surface-enhanced Raman spectroscopy (SERS). An oligonucleotide probe which is specific for 35S DNA promoter target was immobilised onto the gold slide and a RCA reaction was performed. A self-assembled monolayer was formed on gold nanorods using 5,5'-dithiobis(2-nitrobenzoic acid) (DTNB) and the second probe of the 35S DNA promoter target was immobilised on the activated gold coated slide surfaces. Probes on the nanoparticles were hybridised with the target oligonucleotide. Quantification of the target concentration was performed via SERS spectra of DTNB on the nanorods. SERS spectra of target molecules were enhanced through the RCA reaction and the detection limit was found to be 6.3fM. The sensitivity of the developed RCA-SERS method was compared with another method which had been performed without using RCA reaction, and the detection limit was found to be 0.1pM. The correlation between the target concentration and the SERS signal was found to be linear, within the range of 1pM to 10nM for the traditional assay and 100fM to 100nM for the RCA assay. For the developed RCA-SERS assay, the specificity tests were performed using the 35S promoter of Bt-176 maize gene. It was found out that the developed RCA-SERS sandwich assay method is quite sensitive, selective and specific for target sequences in model and real systems.

Identification of a novel papillomavirus in a Northern Fulmar (Fulmarus glacialis) with viral production in cartilage

We report the identification of a novel papillomavirus, Fulmarus glacialis papillomavirus 1 (FgPV1), present within an interdigital foot mass of a Northern Fulmar (Fulmarus glacialis). The mass of interest was composed of normal stratified and keratinized epithelium and dense mesenchymal cells with central cartilaginous islands. Within the nuclei of many chondrocytes were loose aggregates or paracrystalline arrays of virions approximately 50 nm in size. Degenerate polymerase chain reaction was used to identify the virus as a putative papillomavirus, and the entire viral genome of 8132 base pairs was subsequently amplified and sequenced. Analysis revealed canonical papillomavirus architecture, including the early open reading frames E6, E7, E1, and E2 and the 2 late proteins L1 and L2. FgPV1 is most closely related to a cluster of avian and reptilian papillomaviruses as visualized by phylogenetic trees. This observation suggests that papillomavirus virion production can occur in mesenchymal cells.

Equus asinus papillomavirus (EaPV1) provides new insights into equine papillomavirus diversity

We detected a novel papillomavirus (EaPV1) from healthy skin and from sun associated cutaneous lesions of an Asinara (Sardinia, Italy) white donkey reared in captivity in a wildlife recovery centre. The entire genome of EaPV1 was cloned, sequenced, and characterised. Genome is 7467 bp long, and shows some characteristic elements of horse papillomaviruses, including a small untranslated region between the early and late regions and the lack of the retinoblastoma tumour suppressor binding domain LXCXE in E7. Additionally, a typical E6 ORF is missing. EaPV1 DNA was detected in low copies in normal skin of white and grey donkeys of the Asinara Island, and does not transform rodent fibroblasts in standard transformation assays. Pairwise nucleotide alignments and phylogenetic analyses based on concatenated E1-E2-L1 amino acid sequences revealed the highest similarity with the Equine papillomavirus type 1. The discovery of EaPV1, the prototype of a novel genus and the first papillomavirus isolated in donkeys, confirms a broad diversity in Equidae papillomaviruses. Taken together, data suggest that EaPV1 is a non-malignant papillomavirus adapted to healthy skin of donkeys.

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.

Bovine papillomavirus type 13 DNA in equine sarcoids

Equine sarcoids are locally aggressive fibroblastic neoplasms considered to be the most common skin tumors of horses worldwide. Bovine papillomavirus types 1 and 2 have typically been associated with sarcoids in equids. Investigations aiming to identify papillomavirus strains, aside from bovine papillomaviruses 1 and 2, which might be associated with sarcoid lesions, have been lacking. The aim of this article is to report the identification of a third bovine papillomavirus type, bovine papillomavirus 13, associated with equine sarcoids. Six sarcoid lesions were collected from diverse anatomical sites on two horses from southern Brazil. To detect a broad spectrum of papillomavirus strains, eight degenerate primer pairs designed to detect conserved regions on the L1 and E1 genes were tested on the DNA samples. Direct sequencing was then performed on the obtained amplicons, and sequence identities were compared with sequences from all bovine papillomavirus types. The FAP59/FAP64, MY09/MY11, and AR-E1F2/AR-E1R4 sequences generated from the sarcoids were shown to present 99 to 100% identity with bovine papillomavirus 13, a new bovine papillomavirus type previously described in cattle. The results from this study suggest that there is a need to identify bovine papillomavirus type 13 and other papillomavirus strains that might be associated with sarcoids in diverse geographical areas; such investigations might establish the frequency of occurrence of this viral type in these common tumors of equids.

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.

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.

Multiply-primed rolling circle amplification of human papillomavirus using sequence-specific primers

Multiply-primed rolling circle amplification (RCA) is a suitable technique for amplification of circular templates and has been used to identify novel human papillomaviruses (HPV). In this study we develop an efficient RCA for whole genome amplification of HPV using HPV-specific primers in clinical samples and establish a protocol for whole genome sequencing using the Sanger method. Amplification of cloned HPV-genomes by RCA was compared using specific primers against random hexamers. Using HPV-specific primers increased the effectiveness on average 15.2 times and the enrichment of HPV relative to human gDNA on average 62.2 times, as compared to using random hexamer. RCA products were sequenced without need for cloning, even when using low-input amounts. The technique was successfully used on 4 patient samples from FTA cards, to generate whole HPV-genome sequences. Degenerated HPV-specific primers for RCA produce DNA of sufficient quality and quantity suitable for sequencing and other potential downstream analyses.

Viral metagenomics as an emerging and powerful tool in veterinary medicine

New diseases continue to emerge in both human and animal populations, and the importance of animals, as reservoirs for viruses that can cause zoonoses are evident. Thus, an increased knowledge of the viral flora in animals, both in healthy and diseased individuals, is important both for animal and human health. Viral metagenomics is a culture-independent approach that is used to investigate the complete viral genetic populations of a sample. This review describes and discusses the different possible steps of a viral metagenomic study utilizing sequence-independent amplification, high-throughput sequencing, and bioinformatics to identify viruses. With this technology, multiple viruses can be detected simultaneously and novel and highly divergent viruses can be discovered and genetically characterized for the first time. This review also briefly discusses the applications of viral metagenomics in veterinary science and lists some of the viruses discovered within this field.

Marine mammals as sentinel species for oceans and human health

The long-term consequences of climate change and potential environmental degradation are likely to include aspects of disease emergence in marine plants and animals. In turn, these emerging diseases may have epizootic potential, zoonotic implications, and a complex pathogenesis involving other cofactors such as anthropogenic contaminant burden, genetics, and immunologic dysfunction. The concept of marine sentinel organisms provides one approach to evaluating aquatic ecosystem health. Such sentinels are barometers for current or potential negative impacts on individual- and population-level animal health. In turn, using marine sentinels permits better characterization and management of impacts that ultimately affect animal and human health associated with the oceans. Marine mammals are prime sentinel species because many species have long life spans, are long-term coastal residents, feed at a high trophic level, and have unique fat stores that can serve as depots for anthropogenic toxins. Marine mammals may be exposed to environmental stressors such as chemical pollutants, harmful algal biotoxins, and emerging or resurging pathogens. Since many marine mammal species share the coastal environment with humans and consume the same food, they also may serve as effective sentinels for public health problems. Finally, marine mammals are charismatic megafauna that typically stimulate an exaggerated human behavioral response and are thus more likely to be observed.

Metagenomics and the molecular identification of novel viruses

There have been rapid recent developments in establishing methods for identifying and characterising viruses associated with animal and human diseases. These methodologies, commonly based on hybridisation or PCR techniques, are combined with advanced sequencing techniques termed ‘next generation sequencing’. Allied advances in data analysis, including the use of computational transcriptome subtraction, have also impacted the field of viral pathogen discovery. This review details these molecular detection techniques, discusses their application in viral discovery, and provides an overview of some of the novel viruses discovered. The problems encountered in attributing disease causality to a newly identified virus are also considered.

Genomic analysis of the first laboratory-mouse papillomavirus

A papillomavirus (PV) that naturally infects laboratory mice will provide an extremely valuable tool for PV research. We describe here the isolation, cloning and molecular analysis of the first novel laboratory-mouse PV, designated MusPV. This agent, recently identified in the tissues from florid and asymmetrical papillomas on the face of nude mice (NMRI-Foxn1(nu)/Foxn1(nu)), was demonstrated to be transmissible to immunocompetent mice (Ingle et al., 2010). The MusPV genome is 7510 bp in length, is organized similarly to those of other PVs and has at least seven ORFs (E1, E2, E4, E6, E7, L1 and L2). Phylogenetic analysis indicates that MusPV belongs to the π genus together with four other rodent PVs (McPV2, MaPV1, MmiPV and RnPV1). Of the rodent PVs, MusPV appears most closely related to Mastomys coucha PV (McPV2), with 65 % genomic homogeneity and 80 % L1 amino acid similarity. Rodent PVs, except for MnPV1, do not contain any identifiable retinoblastoma protein (RB) binding sites. MusPV has one putative RB-binding site on the E6 protein but not on the E7 protein. Non-coding regions (NCRs) of PVs maintain multiple binding sites for transcription factors (TFs). The NCR of MusPV has numerous sites for TF binding, of which at least 13 TFs are common to all PVs in the π genus. MusPV provides a potentially valuable, novel mouse model to study mechanisms of infection, oncology and novel preventive and therapeutic approaches in mice that can be translated to diseases caused by human PVs.

The canine papillomavirus and gamma HPV E7 proteins use an alternative domain to bind and destabilize the retinoblastoma protein

The high-risk HPV E6 and E7 proteins cooperate to immortalize primary human cervical cells and the E7 protein can independently transform fibroblasts in vitro, primarily due to its ability to associate with and degrade the retinoblastoma tumor suppressor protein, pRb. The binding of E7 to pRb is mediated by a conserved Leu-X-Cys-X-Glu (LXCXE) motif in the conserved region 2 (CR2) of E7 and this domain is both necessary and sufficient for E7/pRb association. In the current study, we report that the E7 protein of the malignancy-associated canine papillomavirus type 2 encodes an E7 protein that has serine substituted for cysteine in the LXCXE motif. In HPV, this substitution in E7 abrogates pRb binding and degradation. However, despite variation at this critical site, the canine papillomavirus E7 protein still bound and degraded pRb. Even complete deletion of the LXSXE domain of canine E7 failed to interfere with binding to pRb in vitro and in vivo. Rather, the dominant binding site for pRb mapped to the C-terminal domain of canine E7. Finally, while the CR1 and CR2 domains of HPV E7 are sufficient for degradation of pRb, the C-terminal region of canine E7 was also required for pRb degradation. Screening of HPV genome sequences revealed that the LXSXE motif of the canine E7 protein was also present in the gamma HPVs and we demonstrate that the gamma HPV-4 E7 protein also binds pRb in a similar way. It appears, therefore, that the type 2 canine PV and gamma-type HPVs not only share similar properties with respect to tissue specificity and association with immunosuppression, but also the mechanism by which their E7 proteins interact with pRb.

Human papillomaviruses (HPVs) are estimated to cause the most common sexually transmitted infection in the world, and these infections are recognized as the major cause of cervical cancer. One of the papillomavirus oncoproteins, E7, plays a major role in both the viral life cycle and progression to cancer. In cells E7 associates and inactivates pRb, a tumor suppressor protein. For the vast majority of papillomaviruses, E7 binds to pRb using a small amino acid sequence, LXCXE. However, we have now identified a papillomavirus E7 protein that lacks the LXCXE domain yet still binds and degrades pRb. This E7 protein, derived from a carcinogenic canine virus, uses its C-terminal domain to bind pRb. In addition, we discovered that a family of papillomaviruses, the gamma type HPVs, also lacks the LXCXE domain and binds pRb using a similar mechanism.

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.

Multiply-primed rolling-circle amplification (MPRCA) of PCV2 genomes: applications on detection, sequencing and virus isolation

Multiply-primed rolling-circle amplification (MPRCA) was used to amplify porcine circovirus type 2 (PCV2) genomes isolated from tissues of pigs with signs of post-weaning multisystemic wasting syndrome (PMWS). Two of the amplified PCV2 genomes were cloned in prokaryotic plasmids and sequenced. Both were nearly identical (1767 nt) except for one silent substitution in the region coding for the capsid protein (ORF2). In addition, they showed high nucleotide sequence similarity with PCV2 isolates from others countries (93-99%). To investigate whether the MPRCA amplified PCV2 genomes could be used to produce infectious virus, the cloned genomes were isolated from the plasmids, recircularized and used for transfection in PK-15 cells. This procedure led to the production of infectious virus to titres up to 10(5.55) TCID(50)/mL. It was concluded that MPRCA is a useful tool to amplify PCV2 genomes aiming at sequencing and virus isolation strategies, where particularly useful is the fact that it allows straightforward construction of PCV2 infectious clones from amplified genomes. However, it was less sensitive than PCR for diagnostic purposes.

Macaca fascicularis papillomavirus type 1: a non-human primate betapapillomavirus causing rapidly progressive hand and foot papillomatosis

Papillomaviruses (PVs) are a group of small, non-enveloped DNA viruses that cause mucosal or cutaneous neoplasia in a variety of animals. Whilst most papillomas will regress spontaneously, some may persist or undergo malignant transformation. In this study, aggressive, persistent and extensive warts were observed on the hands and feet of a cynomolgus macaque (Macaca fascicularis). The presence of PV in the wart biopsies was identified by immunohistochemistry and PCR amplification of PV DNA. The genomic DNA of this PV was cloned and sequenced, and the PV was designated M. fascicularis papillomavirus type 1 (MfPV-1). Its genome was 7588 bp in length and the organization of its putative open reading frames (E1, E2, E6, E7, L1, L2 and E4) was similar to that of other PVs. MfPV-1 had a short non-coding region (NCR) of 412 bp. Molecular analysis of MfPV-1 genomic DNA classified it into the genus Betapapillomavirus, to which all epidermodysplasia verruciformis (EV)-type PVs belong. Diseases caused by PVs of the genus Betapapillomavirus are usually associated with natural or iatrogenic immunosuppression. The genomic characterization performed in this study showed that MfPV-1 clustered within the genus Betapapillomavirus and also contained EV-type-specific motifs in its NCR. Further characterization of this virus and its host interactions may allow us to develop a non-human primate model for human betapapillomaviruses, a genus populated by human PV types causing EV.

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.

Genomic characterization of novel dolphin papillomaviruses provides indications for recombination within the Papillomaviridae

Phylogenetic analysis of novel dolphin (Tursiops truncatus) papillomavirus sequences, TtPV1, -2, and -3, indicates that the early and late protein coding regions of their genomes differ in evolutionary history. Sliding window bootscan analysis showed a significant a change in phylogenetic clustering, in which the grouped sequences of TtPV1 and -3 move from a cluster with the Phocoena spinipinnis PsPV1 in the early region to a cluster with TtPV2 in the late region. This provides indications for a possible recombination event near the end of E2/beginning of L2. A second possible recombination site could be located near the end of L1, in the upstream regulatory region. Selection analysis by using maximum likelihood models of codon substitutions ruled out the possibility of intense selective pressure, acting asymmetrically on the viral genomes, as an alternative explanation for the observed difference in evolutionary history between the early and late genomic regions of these cetacean papillomaviruses.

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

Papillomatosis has been documented in several carnivores, and papillomavirus (PV) types have been characterized from lesions in a number of carnivore species: the canine oral PV (COPV), the Felis domesticus PV type 1 (FdPV-1) isolated from a Persian cat, the Procyon lotor PV type 1 (PlPV-1) isolated from a raccoon, the canine PV type 2 (CPV-2) from a dog's foot pad lesion and the canine PV type 3 (CPV-3) associated with a canine epidermodysplasia verruciformis - like disease. A tissue sample was taken from a papillomatous lesion on the oral mucosa of a polar bear (Ursus maritimus). Extracted DNA was used as a template for multiply primed rolling-circle amplification (RCA), and restriction enzyme analysis of the RCA product indicated the presence of papillomaviral DNA. The genome of this PV was cloned and the complete genomic sequence was determined. The Ursus maritimus PV type 1 (UmPV-1) genome counts 7582 basepairs and is smaller than that of other papillomaviruses from carnivore species. UmPV-1 contains the typical noncoding region NCR1, but unlike the carnivore PVs of the Lambda genus, UmPV-1 does not possess a second noncoding region NCR2. Phylogenetic analysis based on a nucleotide sequence alignment of the L1 ORF of UmPV-1 and 51 other PV types indicates that UmPV-1 does not cluster with any of the other carnivore PVs, but branches off near the root of the common branch of the genus Alphapapillomavirus.

Ancient papillomavirus-host co-speciation in Felidae

The evolutionary rate of feline papillomaviruses is inferred from the phylogenetic analysis of their hosts, providing evidence for long-term virus-host co-speciation

Background

Estimating evolutionary rates for slowly evolving viruses such as papillomaviruses (PVs) is not possible using fossil calibrations directly or sequences sampled over a time-scale of decades. An ability to correlate their divergence with a host species, however, can provide a means to estimate evolutionary rates for these viruses accurately. To determine whether such an approach is feasible, we sequenced complete feline PV genomes, previously available only for the domestic cat (Felis domesticus, FdPV1), from four additional, globally distributed feline species: Lynx rufus PV type 1, Puma concolor PV type 1, Panthera leo persica PV type 1, and Uncia uncia PV type 1.

Results

The feline PVs all belong to the Lambdapapillomavirus genus, and contain an unusual second noncoding region between the early and late protein region, which is only present in members of this genus. Our maximum likelihood and Bayesian phylogenetic analyses demonstrate that the evolutionary relationships between feline PVs perfectly mirror those of their feline hosts, despite a complex and dynamic phylogeographic history. By applying host species divergence times, we provide the first precise estimates for the rate of evolution for each PV gene, with an overall evolutionary rate of 1.95 × 10^-8 (95% confidence interval 1.32 × 10^-8 to 2.47 × 10^-8) nucleotide substitutions per site per year for the viral coding genome.

Conclusion

Our work provides evidence for long-term virus-host co-speciation of feline PVs, indicating that viral diversity in slowly evolving viruses can be used to investigate host species evolution. These findings, however, should not be extrapolated to other viral lineages without prior confirmation of virus-host co-divergence.

Complete genomes and phylogenetic positions of bovine papillomavirus type 8 and a variant type from a European bison

The DNA of bovine papillomavirus (BPV) type 8 was extracted from papillomas on cattle kept in Japan, and DNA of bovine papillomavirus BPV-8-EB was extracted from a European bison (Bison bonasus) born in Italy and released into the wild in Slovakia. The DNA genomes of these BPVs were amplified using multiply primed rolling circle amplification and polymerase chain reaction, then characterized by direct sequencing method. The BPV-8 and BPV-8-EB genomes consisted of 7,791 base pairs (bp) and 7,773 bp, respectively (GenBank accession numbers DQ098913 and DQ098917). The nucleotide sequence similarity of these BPVs indicated that BPV-8-EB was a variant of BPV-8. In the genome of BPV-8-EB, one nucleotide substitution was found in the E2 and E5 open reading frame (ORF) and upstream regulatory region (URR), and a short deletion and addition were found in the URR. The high similarity of sequences between the BPV-8 to BPV-5 in total genome (70%) and L1 ORF (75%) as well as a phylogenetic analysis were the bases for classifying BPV-8 in the genus Epsilon papillomavirus. The BPV-8 E6 and E7 ORFs/proteins also showed some characteristic features of genus Epsilon papillomavirus. However, BPV-8 contained E4 ORF, which was not found in BPV-5. In addition, the secondary structure of E5 proteins of BPV-5 and BPV-8 suggested that these proteins may have cell-transforming ability.

Complete genome and phylogenetic position of bovine papillomavirus type 7

Six bovine papillomavirus (BPV) types and 16 putative BPV types have been reported previously. Here, the complete genome sequence of BAPV6, a novel putative BPV type isolated from cattle in Japan, was determined by using multiple-primed rolling-circle amplification. The genome consisted of 7412 bp (G+C content of 46 mol%) that encoded five early (E1, E2, E4, E6 and E7) and two late (L1 and L2) genes, but did not encode the E5 gene. The E6 protein contained a non-consensus CxxC(x)33CxxC and a consensus CxxC(x)29CxxC zinc-binding domain, and the E7 protein lacked the LxCxE motif. The nucleotide sequence of the L1 open reading frame (ORF) was related most closely (57–58 %) to the L1 ORF of member(s) of the genera Betapapillomavirus, Gammapapillomavirus and Pipapillomavirus. Phylogenetic analysis based on the complete L1 ORF suggests that BAPV6 should be classified in a novel genus in the family Papillomaviridae as BPV-7.

Phylogeny and evolution of papillomaviruses based on the E1 and E2 proteins

Papillomaviridae are a family of small double-stranded DNA viruses that infect stratified squamous epithelia in vertebrates. Members of this family are causative agents of malignant tumours, such as cervical cancer while others are associated with benign proliferative lesions. So far, Papillomaviruses (PVs) are classified according to the sequence identity in the capsid gene L1. However, evidence has accumulated indicating a discontinuity in the evolutionary history of the L1 and L2 genes of many PVs, giving rise to differences in the phylogenetic reconstructions of the early and of the late genes. Neither the oncogenes E5, E6 and E7 nor the upstream regulatory region are suitable for phylogenetic inference due to the poor conservation along the Papillomaviridae family. We have analysed here the evolutionary relationships of the PVs with respect to the E1 and E2 proteins, and the results provide both phylogeny and biologic behaviour of the viruses. The hierarchical taxonomic relationships can be structured as an alternative classification system in which mucosal high-risk viruses, mucosal low-risk viruses and viruses associated with cutaneous lesions are grouped separately and do not appear intermingled. Some important trends are also observed: first, evolution of the PVs has not been homogeneous, even in viruses that infect the same host, and second mucosal human PVs have evolved faster than their cutaneous counterparts. The evolutionary analysis based on the E1 and E2 proteins will allow us to better understand the generation of the diversity of the PVs and the development of malignancy associated with these viruses.