E5 nucleotide polymorphisms suggest quasispecies occurrence in BPV-1 sub-clinically infected horses

Complete genomic characterization of bovine papillomavirus type 1 and 2 strains infers ongoing cross-species transmission between cattle and horses

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..

New approach for genomic characterisation of equine sarcoid-derived BPV-1/-2 using nanopore-based sequencing

BACKGROUND

Bovine papillomavirus (BPV) types 1 and 2 play a central role in the etiology of the most common neoplasm in horses, the equine sarcoid. The unknown mechanism behind the unique variety in clinical presentation on the one hand and the host dependent clinical outcome of BPV-1 infection on the other hand indicate the involvement of additional factors. Earlier studies have reported the potential functional significance of intratypic sequence variants, along with the existence of sarcoid-sourced BPV variants. Therefore, intratypic sequence variation seems to be an important emerging viral factor. This study aimed to give a broad insight in sarcoid-sourced BPV variation and explore its potential association with disease presentation.

METHODS

In order to do this, a nanopore sequencing approach was successfully optimized for screening a wide spectrum of clinical samples. Specimens of each tumour were initially screened for BPV-1/-2 by quantitative real-time PCR. A custom-designed primer set was used on BPV-positive samples to amplify the complete viral genome in two multiplex PCR reactions, resulting in a set of overlapping amplicons. For phylogenetic analysis, separate alignments were made of all available complete genome sequences for BPV-1/-2. The resulting alignments were used to infer Bayesian phylogenetic trees.

RESULTS

We found substantial genetic variation among sarcoid-derived BPV-1, although this variation could not be linked to disease severity. Several of the BPV-1 genomes had multiple major deletions. Remarkably, the majority of them cluster within the region coding for late viral genes. Together with the extensiveness (up to 603 nucleotides) of the described deletions, this suggests an altered function of L1/L2 in disease pathogenesis.

CONCLUSIONS

By generating a significant amount of complete-length BPV genomes, we succeeded to introduce next-generation sequencing into veterinary research focusing on the equine sarcoid, thus facilitating the first report of both nanopore-based sequencing of complete sarcoid-sourced BPV-1/-2 and the simultaneous nanopore sequencing of multiple complete genomes originating from a single clinical sample.

Genomic characterisation of bovine papillomavirus types 1 and 2 identified in equine sarcoids in Japan

BACKGROUND

Bovine papillomavirus types 1 and 2 (BPV1/2) infection in horses has been associated with the development of equine sarcoids. Previous findings revealed the presence of sarcoid-associated BPV sequence variants that have been proposed as a key factor of cross-species infection in horses. To verify this hypothesis, sarcoid-associated BPV variants should be identified regardless of geographic location.

OBJECTIVES

Sequence analyses of BPV1/2 derived from both horses and cattle were conducted to clarify the sarcoid-associated sequence variants. The aim of this study was to clarify the correlation between BPV phylogeny and the geographic origin/host species.

STUDY DESIGN

Cross-sectional study.

METHODS

Conventional PCR to detect BPV1/2 was performed with genomic DNA extracted from equine sarcoid (n = 10) and bovine papilloma (n = 10) samples collected in Japan. Direct sequencing results were compared between equine and bovine (equine/bovine)-derived BPV to identify sarcoid-associated variants of two early regions (E2, E5), one late region (L1) and the long control region (LCR). Phylogenetic and phylogeny-trait correlation were analysed using Bayesian Markov chain Monte Carlo (MCMC) method and Bayesian tip-association significance testing (BaTS).

RESULTS

Seven BPV1 and three BPV2 were identified from equine sarcoids using PCR and direct sequencing. Sequence analysis of equine/bovine-derived samples showed no sarcoid-associated variants in four regions (E2, E5, L1 and LCR) of either BPV1 or BPV2. The phylogenetic tree of BPV1 E2, L1 and LCR tended to cluster within its geographic origins. BaTS analysis demonstrated that BPV1 sequence variability may be due to the geographic origin rather than host species difference.

MAIN LIMITATIONS

There was a limitation in sample numbers.

CONCLUSIONS

This study supports the geographic-specific hypothesis of sequence variability, suggesting that BPV1 is shared between local equids and bovids. However, more extensively collected sequences worldwide and functional evaluations are needed to verify the geographic-specific sequence variability of BPV1/2 between equine- and bovine-derived sequence.

Origin and evolution of papillomavirus (onco)genes and genomes

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

Molecular Detection of Bovine Papillomavirus DNA in the Placenta and Blood of Healthy Mares and Respective Foals

Despite the characteristic species specificity of Papillomaviruses (PVs), the bovine papillomavirus (BPV) types 1, 2, and-more rarely-13, can cross-infect equids, where they are involved in the pathogenesis of sarcoid neoplasms. Sarcoids are locally invasive fibroblastic skin tumors that represent the most common skin neoplasms in horses worldwide. The transmission mechanism of BPV is still controversial in horses. Thus far, direct and indirect routes have been implicated, while vertical transmission has been suggested after the detection of viral DNA in the semen of healthy stallions. Testing of the blood and placenta of non-sarcoid baring mares and their respective foals revealed that the equine placenta can harbor BPV DNA, leading us to speculate a possible prenatal vertical DNA transmission in equids.

The possible role of Stomoxys calcitrans in equine sarcoid transmission

The association between bovine papillomavirus (BPV) and equine sarcoids is well established, but it is unclear how the virus spreads. Although evidence in support of viral spread through direct animal contact exists, this does not explain sarcoid development in isolated equids. BPV DNA has been detected in flies, which could indicate that these insects serve as a vector. This study aimed to investigate whether BPV-negative stable flies (Stomoxys calcitrans) become positive for BPV DNA after exposure to equine sarcoid or bovine papilloma tissue under experimental conditions and, if so, for how long. A total of 420 stable flies were caught alive and exposed to BPV positive equine sarcoid or bovine papilloma tissue. During the following week, dead flies were collected daily and BPV loads were determined by quantitative PCR. There was a significant rise in BPV load after tissue exposure both in sarcoid and papilloma exposed flies, but the viral load was higher and remained high for a longer time after exposure to papilloma tissue compared to sarcoid tissue. Within days, viral loads decreased again and became indifferent from loads before exposure. The results of these experiments indicate that BPV transmission by S. calcitrans seems possible and is more likely to occur after contact with bovine papillomas than with equine sarcoids. Transmission seems only possible shortly after tissue exposure. Further research could include experimental induction of sarcoids with BPV positive stable flies, or a repeat of the experiment with micro-dissection prior to PCR.

Genomic comparison of bovine papillomavirus 1 isolates from bovine, equine and asinine lesional tissue samples

Several attempts have been made to categorize equid- and bovid-specific bovine papillomavirus 1 (BPV1) isolates based on sequence tags. This study includes newly determined sequence information from 33 BPV1 isolates of equine, asinine and bovine origin and investigates sequence bias due to host species. Twenty of the viral genomes were sequenced over their entire length and a further thirteen were sequenced, including flanking sequences, at two specific sites, the LCR and the E5 ORF. Alignment and analyses of the sequences did not reveal statistically significant site differences between the sequences of bovine and equid origin. None of the proposed sites of divergence noted by other authors demonstrated significant species-specific characteristics. Our results suggest that BPV1 is shared between equine, asinine and bovine host species, and that viral transfer between bovines and equids is a repeated and ongoing phenomenon.