ABSTRACT
Infection with bovine papillomavirus (BPV) types 1 and 2 results in the most common skin tumor of horses, termed equine sarcoid. The persistent and recurrent nature of this tumor stands in contrast to the regressive nature of BPV-1/- 2 induced cutaneous papillomas in cattle. The circulation of horse-specific BPV-1/- 2 variants within equine populations has been suggested as a possible explanation for the difference in clinical presentation of BPV-1/- 2 infection between horses and cattle. In order to investigate this hypothesis, we identified 98 complete BPV-1/- 2 genomes using a Nanopore sequencing approach. Separate BPV-1/- 2 alignments were used to infer Bayesian phylogenetic trees. Phylogeny-trait association concerning host species was investigated using Bayesian Tip-association Significance software (BaTS) Overall, 179 unique BPV-1 and 128 BPV-2 substitutions were found. The E2 coding region in the viral genome exhibited an exceptionally high rate of non-synonymous mutations (81 %, n = 13/16). Interestingly, extensive deletions in the L1/L2 region (up to 1.5 kb) were found exclusively in horse-derived samples. Nevertheless, the most frequently detected single nucleotide polymorphisms were shared between equine and bovine hosts, which is in agreement with BaTS results indicating no phylogeny-host correlation. We found indications that horse-specific mutations might exist in subpopulations of equine derived BPV-1/- 2, but these did not result in horse-adapted genetic variants. Based on these observations, cross-species transmission from cattle to horses seems to be an ongoing process, rather than an ancient occurrence that has been followed by the circulation of horse-adapted BPV variants in the horse population..
