Phylogeny and amino acid analysis in single and mixed bovine papillomavirus infections in Southern Brazil, 2016-2020

Bovine papillomaviruses (BPVs) exhibit a high degree of genetic variability, and several viral types have been identified based on analysis of the L1 gene. The L1 is the main capsid protein and the main target for neutralizing antibodies. We performed a retrospective study on BPVs circulating in Rio Grande do Sul state, Southern Brazil, in 2016-2020. DNA from 43 bovine papilloma samples were amplified using two degenerate primer sets – FAP59/64 and MY09/11 – targeting the L1 region, and analyzed for phylogeny, mixed BPV infections (coinfections) and amino acid (aa) sequences. We also performed an in silico analysis with 114 BPV L1 sequences from the GenBank database to assess the agreement between the phylogeny obtained based on complete L1 sequences versus that based on the region amplified using the FAP59/64 and MY09/11 primer sets. Considering single and coinfections, we identified 31 BPV-1 (31/43; 72.1%), 27 BPV-2 (27/43; 62.8%) and 4 BPV-6 (4/43; 9.3%). Coinfections with BPV-1 and BPV-2 were observed in 61.3% of the samples. Our results are supported by in silico analyses that demonstrate that the classification using FAP59/64 or MY09/11 matches the complete L1 results, except for BPV-17 and -18, which may be mistakenly classified depending on the primers used. Furthermore, we found unique or rare amino acids in at least one L1 sequence of each BPV type identified in our study, some of which have been identified previously in papillomavirus epitopes, suggesting immune-mediated selection. Finally, our study provides an overview of BPVs circulating in Southern Brazil over the last five years and point to the combined use of primers FAP59/64 and MY09/11 for analysis of BPV coinfections and putative epitopes.