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.

The discovery of novel viruses has often been accomplished by using hybridization-based methods that necessitate the availability of a previously characterized virus genome probe or knowledge of the viral nucleotide sequence to construct consensus or degenerate PCR primers. In their natural replication cycle, certain viruses employ a rolling-circle mechanism to propagate their circular genomes, and multiply primed rolling-circle amplification (RCA) with φ29 DNA polymerase has recently been applied in the amplification of circular plasmid vectors used in cloning. We employed an isothermal RCA protocol that uses random hexamer primers to amplify the complete genomes of papillomaviruses without the need for prior knowledge of their DNA sequences. We optimized this RCA technique with extracted human papillomavirus type 16 (HPV-16) DNA from W12 cells, using a real-time quantitative PCR assay to determine amplification efficiency, and obtained a 2.4 × 10 ⁴ -fold increase in HPV-16 DNA concentration. We were able to clone the complete HPV-16 genome from this multiply primed RCA product. The optimized protocol was subsequently applied to a bovine fibropapillomatous wart tissue sample. Whereas no papillomavirus DNA could be detected by restriction enzyme digestion of the original sample, multiply primed RCA enabled us to obtain a sufficient amount of papillomavirus DNA for restriction enzyme analysis, cloning, and subsequent sequencing of a novel variant of bovine papillomavirus type 1. The multiply primed RCA method allows the discovery of previously unknown papillomaviruses, and possibly also other circular DNA viruses, without a priori sequence information.

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