key: cord-0010994-wm0nk0k3
authors: Murahwa, Alltalents T.; Tshabalala, Mqondisi; Williamson, Anna-Lise
title: Recombination Between High-Risk Human Papillomaviruses and Non-Human Primate Papillomaviruses: Evidence of Ancient Host Switching Among Alphapapillomaviruses
date: 2020-05-08
journal: J Mol Evol
DOI: 10.1007/s00239-020-09946-0
sha: 700624f1fe8380edda093ea9b0895c96404ada18
doc_id: 10994
cord_uid: wm0nk0k3

We use all the currently known 405 Papillomavirus (PV) sequences, 343 curated PV sequences from both humans and animals from the PAVE data base, to analyse the recombination dynamics of these viruses at the whole genome levels. After showing some evidence of human and non-human primate PV recombination, we report a comprehensive recombination analysis of all currently known 82 Alphapapillomaviruses (Alpha-PVs). We carried out an exploratory study and found novel recombination events between High-Risk HPV Types and Macaca fascicularis PV1 (MfPV1), Macaca Fuscata PV2 (MfuPV2) and Pan Paniscus PV1 (PpPV1) Papillomaviruses. This is the first evidence of interactions between PVs from different hosts and hence postulates the likelihood of ancient host switching among Alpha-PVs. Notwithstanding these results should be interpreted with caution because the major and minor parents indicated by RDP4 program are simply the sequences in the alignment that most closely resemble the actual parents. We found statistically significant differences between the phylogenies of the PV sequences with recombination regions and PV sequences without recombination regions using the Shimodaira–Hasegawa phylogenetic incongruence testing. We show that not more than 76MYA Alpha-PVs were in the same biological niche, a pre-requisite for recombination, and as the hosts evolved and diversified, the viruses adapted to specific host niches which eventually led to coevolution with specific hosts before speciation of primate species. Thus providing evidence that in ancient times no earlier than the Cretaceous period of the Mesozoic age, Alpha-PVs recombined and switched hosts, but whether this host switching is occurring currently is unknown. However, a clearer picture of the PVs evolutionary landscape can only be achieved with the incremental discovery of PV sequences, especially from the animal kingdom. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1007/s00239-020-09946-0) contains supplementary material, which is available to authorized users.

The Papillomaviridae (PV) family was once part of the larger family of Papovaviridae which was split into Polyomaviridae and Papillomaviridae by the International Committee on Taxonomy of Viruses (Van Regenmortel 2000) .

Handling Editor: Keith Crandall.

The online version of this article (https ://doi.org/10.1007/s0023 9-020-09946 -0) contains supplementary material, which is available to authorized users.

According to the most recent ICTV classification, the PV family includes two subfamilies Firstpapillomavirinae with 52 genera and Secondpapillomavirinae with one genus and one species (Van Doorslaer et al. 2018) . Genera are named according to the Greek alphabet from alpha to omega, and following exhaustion of the alphabet the term dyo and treis (Greek for second time and third, respectively) coined to accommodate the extra genera e.g. dyo-deltapapillomaviruses (Bernard et al. 2010) . There are about 405 known reference PVs uptodate as listed in the curated PAVE database (https ://pave.niaid .nih.gov/), 226 of these are human papillomaviruses (HPVs) and 179 from different animal species (Van Doorslaer and Dillner 2019). HPVs are the most studied and are distributed over 5 genera (Alpha, Beta, Gamma, Mu and Nu) . The other PV genera are from other mammals, birds and reptiles (Bernard et al. 2010) . Below the genus level are species and below the species level are PV types (de Villiers et al. 2004 ). Different genera have less than 60% similarity within the L1 gene, while species share between 60 and 70% similarity and types share between 71 and 89% similarity. The ICTV is responsible for nomenclature of viruses down to species level, and below species level, the International Human Papillomavirus (HPV) Reference Centre (Karonliska Institute, Sweden) and the Animal Papillomavirus Reference Centre (University of Arizona, USA) assigns unique PV-type numbers after the complete genome has been sequenced, cloned and confirmed by the Centre (Mühr et al. 2018) .

Alphapapillomaviruses (Alpha-PVs) genus is the most widely studied of all the PV genera. According to the PAVE (https ://pave.niaid .nih.gov/) database, there are currently 82 Alpha-PV types classified into 14 species (Alpha-1 to . Of the 82 known types, 17 are animal Alpha-PVs but only from non-human primates and 65 from humans.

Traditionally, PVs have been thought to evolve slowly because they replicate by co-opting high-fidelity host cellular DNA polymerases that have an error rate of about 4.3 × 10 − 5 substitutions per nucleotide site per year (Korona et al. 2011) . The general assumption is that PVs have coevolved with their hosts (Chen et al. 2009a; Bravo and Felez-Sanchez 2015; Dube Mandishora et al. 2018 ). Selection pressures due to host immune responses differ among PV genes and result in these genes displaying different substitution rates. Further, the cellular polymerases of different host species may differ in their degree of fidelity such that virus lineages infecting different hosts might display different substitution rates (de Oliveira et al. 2015) . The popular views about PV evolution have been (1) they are static and slow evolving viruses with very low mutation rates (2) host linked co-divergence (3) strict tissue tropism (4) host specific and (5) lack of recombination (Bravo et al. 2010) . However, the understanding of PV evolution has advanced over the years, with various alternative mechanisms such as inter-species transmission, adaptive radiation and recombination (Angulo and Carvajal-Rodriguez 2007; Carvajal-Rodriguez 2008; Gottschling et al. 2011; Robles-Sikisaka et al. 2012) .

The strict definition of recombination incorporates reciprocation, meaning the recipient of a genome portion acts as a donor of the replaced portion in the source, which is not the case with PV recombination (Pérez-Losada et al. 2015) . PV recombination could be more appropriately named, gene conversion, but it has been so widely used that changing it would introduce confusion to the whole subject. Recombination has a potential major impact on PV evolution, pharmacogenomics and vaccine development. In other viruses, recombination has been associated with emergence of novel viruses, increase in virulence and pathogenesis, changes in tissue tropisms and expansion of viral host ranges (Martin et al. 2011; Simon-Loriere and Holmes 2011) .

The biological plausibility of PV recombination is occasioned by the genetic multiplicity of PVs and the high frequencies of observed HPV co-infections (Angulo and Carvajal-Rodriguez 2007) . However, the study of PV recombination has been hampered by technical difficulties associated with the accurate alignment of highly diverse PV gene sequences (Posada and Crandall 2001) . One of the most commonly used approaches to recombination detection is the use of the various recombination analysis tools implemented within the RDP4 software package (Martin and Rybicki 2000) . During recombination detection, RDP4 rigorously tests the quality of sequence alignments to guard against the detection of false-positive recombination signals that arise due to sequence misalignment (Varsani et al. 2006) .

We report here the use of all 405 known PV sequences, 343 curated PV sequences from both humans and animals, to analyse the recombination dynamics of these viruses at the whole genome levels. Specifically, we use these sequences to identify recombination and to determine whether there is intra-genus, inter-species and inter-host species recombination. After stumbling upon evidence of human and nonhuman primate PV recombination, we report also a comprehensive recombination analysis of all currently known 82 Alpha-PVs. Consequently, we carried out the final test for recombination, phylogenetic incongruence, to prove a statistically significant difference between the phylogenies of the PV sequences used.

This was an exploratory study to investigate the likelihood of ancient host-switching events among PVs as previously postulated in literature (Chen et al. , 2019 .

All the 405 currently known PV whole-genome reference sequences from the PAVE database were downloaded in Fasta format (https ://pave.niaid .nih.gov/#searc h/searc h_datab ase). See also Supplementary Material 1a for the fasta format of the sequences. Additionally, Fasta format of the concatenated E1-E2-L1-L2 alignment of 343 sequences was downloaded from the ICTV website (https ://talk.ictvo nline .org/ictv-repor ts/ictv_onlin e_repor t/dsdna -virus es/w/ papil lomav irida e/1214/papil lomav irida e---v2019 11), (see also Supplementary Material 1b for the fasta alignment of the sequences) and all the 82 currently known Alpha-PVs whole-genome sequences were obtained from PAVE database see also Supplementary Material 1c for the fasta format of the sequences.

This was a hierarchical approach in which we started with the whole genomes of all the currently known PVs (405)-R1, then we used a set of all curated genomes 343 from the ICTV database where the alignment was actually downloaded as mentioned in Table 1 , R2, with concatenated gene regions in order to improve the alignment. Then, lastly, we use all the 82 known Alpha-PVs, R3, based on the fact that R1 and R2 had shown recombination among Alpha-PVs.

All the 405 currently known whole-genome PV reference sequences from the PAVE database were included in the initial recombination analysis (R1). We then constructed an alignment of the 405 PV sequences using MUSCLE and the CLUSTALW output was used for the R1 recombination analysis. This alignment was analysed using RDP v4.95 (Martin and Rybicki 2000) (with default settings) which implements analysis of recombination using a suite of 7 recombination detection methods or algorithms: RDP (Martin and Rybicki 2000) , BOOTSCAN (Martin et al. 2005a, b) , CHIMAERA (Martin et al. 2005a, b) , GENE-CONV (Padidam et al. 1999) , MAXIMUM X 2 (Smith 1992 ) and SISCAN (Gibbs et al. 2000) . Only recombination events that were identified by at least four methods and no showed warning flags (messages generated by the software to warn the reader of the possible reasons for each outcome, and to interpret them with caution) in the RDP software were considered.

FASTA format of the concatenated E1-E2-L1-L2 alignment of the 343 sequences (downloaded from ICTV, https ://talk.ictvo nline .org/ictv-repor ts/ictv_onlin e_repor t/ dsdna -virus es/w/papil lomav irida e/1214/papil lomav irida e---v2019 11) were used in the recombination analysis R2. This alignment was also analysed using RDP v4.95 (Martin and Rybicki 2000) (with default settings) as was done for R1 above, and only recombination events that were identified by at least four methods and no showed warning flags in the RDP software were considered.

Eighty-two (82) Alpha-PVs complete genomes (65 Alpha-HPVs and 17 non-human primate Alpha-PVs), which are of all the currently known Alpha-PVs, were obtained from the PAVE database. We then constructed an alignment containing the 82 Alpha-PVs using MUSCLE, the CLUSTALW output was used for the R3 recombination analysis as a followup investigation to findings from R1 and R2. This alignment was also analysed using RDP v4.95 (Martin and Rybicki 2000) (with default settings) as was done for R1 and R2 above, and only recombination events that were identified by at least four methods and no showed warning flags in the RDP software were considered.

The Shimodaira-Hasegawa (SH) test (Shimodaira and Hasegawa 1999) using W-IQ-TREE (Trifinopoulos et al. 2016) . We used CLUSTALW alignments (Sievers et al. (Trifinopoulos et al. 2016 ). The tool tests tree topology, estimates model parameters such as substitution rates and optimizes tree branch lengths to lessen computational usage. We used default settings of the W-IQ-Tree, including best fit model (Kalyaanamoorthy et al. 2017 ) and ultra-fast bootstrap analysis (1000 alignments) (Minh et al. 2013) to run tree topology analysis including the Kishino-Hasegawa (KH) test (Kishino and Hasegawa 1989) , Shimodaira-Hasegawa (SH) test (Shimodaira and Hasegawa 1999) and approximately unbiased (AU) test (Shimodaira 2002) to test if there is a difference in evolutionary patterns among trees generated after removing the recombinant regions from the original sequences.

All the trees performed were denoted A1 and A2, A1 representing the trees generated from the original sequence alignments and A2 denoting the trees generated after removing the recombinant regions from the original sequences, after the recombination analyses. Alignments without recombinant regions were generated automatically from RDP v4.95 (Martin and Rybicki 2000) after each of the recombination analyses.

The sequences of 82 human and animal Alpha-PVs whole genomes (Alphapapillomavirus genus) were aligned using MUSCLE v7.221 (Edgar 2004) . A maximum-likelihood tree of the nucleotide sequences was generated in PhyML using the optimal model of evolution (GTR+G) as determined within MEGA 7 (Kumar et al. 2016) . The Newick format of the tree was uploaded and modified in iTOL https ://itol. embl.de/tree/.

There were a total of 393 recombination events but only 4 were sufficiently supported by at least 4 algorithms as indicated in the priori criterion set in the methodology and no software flags, see Table 2 below for a summary of the 4 events, showing parent and recombinant PV types, the number of methods/algorithms to support the event and positions of the breaking points were the gene conversions occurred the same is true for Tables 3 and 4 for R2 and R3, respectively. (See also Supplementary Fig. 2) .

Event 1 is between cetacean PVs and involve TtPV1 of the Upsilonpapillomavirus genus as a major parent to recombinants TtPV5 and TtPV6 of the Ominkronpapillomvirus genus. This becomes the first evidence of inter-genus recombination, in this case between Upsilon-and Omikronpapillomaviruses. Event 3 is between low-risk human Alpha-PVs. Event 4 is between an Alpha-HPV 68 of the alpha-7 species as a recombinant and a non-human primate PV PpPV1 (Pan Paniscus PV1 from pygmy chimpanzees commonly known as Bonobos, alpha-10 species) as the major parent. This becomes the first evidence of inter-host PV recombination, yet still inter-genus. All the events were sufficiently supported by at least 5 methods. Event 4 prompted further investigation into other recombination events that can occur across PVs that infect different host species. A wellcurated set of 343 sequences was thus used to explore this possibility.

Recombination signals were detected across both the early regions (E genes) and the late regions (L genes) of the concatenated alignment of the PV genomes. There were a total of 456 events detected by the software but only 7 events were sufficiently supported by at least 4 recombination detection algorithms. See Table 3 below for a summary of the 7 events and Supplementary Fig. 3 . Events 1, 2 and 4 are between cetacean PVs. DdPV1 (Delphinus delphis PV1) from the common short-nose dolphin belongs to the Upsilonpapillomavirus-1 species together with TtPV1 (Tursiops truncatus PV1) from the common bottlenose dolphin. TtPV5, TtPV6, PphPV1 (Phocoena phocoena PV1; harbour porpoise) and PsPV1 (Phocoena spinipinnis PV1; Burmeister porpoise) all belong to the Omikronpapillomavirus-1 species. It is apparent from events 1, 2 and 4 that TtPV1 is a minor parent in all the recombination events. Events 1 and 2 are well supported by 7 recombination algorithms while event 4 is supported by 4 algorithms. Event 4 is also an exception as it spans from the late and a small portion of the early region.

Events 3, 5 and 6 are between HPV types. HPV54 a lowrisk HPV type is a major parent in all the events, and only in event 6 it is a parent to HPV82 which is classified as a high-risk HPV type. All the HPV recombination events are intra-genus.

Event 7 is between HPV70 as the recombinant and a nonhuman primate, MfPV7 (Macaca Fascicularis PV7, commonly known as cynomolgus macaque) as the major parent. MfPV7 is from the alpha-12 species while HPV70 is of the alpha-7 species. Event 7 also spans about 692 bp in the late region, about half the size of either the L1 or L2. Event 4 of R1 and event 7 of R2 all pointed to the fact that recombination occurred only among Alpha-PVs, this prompted a third analysis R3, of all currently known Alpha-PVs.

There were a total of 117 events detected by the software but only 10 events were sufficiently supported by at least 4 recombination detection algorithms. See Table 4 below for a summary of the 10 events and Supplementary Fig. 4 . Events 1, 3, 4, 5, 6, 7 and 10 are solely between HPV types, with events 3, 4, 6, 7 and 10 involving high-risk HPV (HR-HPV) types; HPV82, HPV51, HPV45, HPV53 and HPV59 as recombinants, respectively. Events 2, 8 and 9 involve recombination between HPVs and non-human primate Alpha-PVs. In event 2 and 8, h-HPV types HPV39 and HPV68 of the alpha-7 species are recombinants of one major parent PpPV1 (Pan Paniscus PV1, alpha-10 species). In event 9, h-HPV 66 of the alpha-6 species is a recombinant of LR-HPV54 of the alpha-13 species as a major parent and minor parent MfuPV2 (Macaca Fuscata PV2, from Japanese macaques) of the alpha-12 species. All the recombination events are intra-genus but mostly cross-species. All the HPV44, HPV87/HPV59 6 4.695 × 10 02 -9.099 × 10 08 5610-6523 events but event 8 were supported by at least 5 algorithms in the RDP v4.95 software (Martin and Rybicki 2000) . Generally, based on the analysis from R1, R2 and R3 there are 21 recombination events and of these 16 are among high-risk Alpha-PVs, 5 events are between HPVs and nonprimate PVs and 5 among cetacean PVs. Event 3 in R1 is the same as event 5 in R2 (save that the minor parent is HPV71), also the same as event 1 in R3. This common recombination event shows HPV30 as the recombinant and HPV54 as the major parent. HR-HPV54 is also a major parent in 3 other recombination events, event 3 and 6 in R2 and event 9 in R3. Thus, HPV54 is a major parent in a total of six recombination events in all the three analyses. HPV82 is a recombinant in 2 events (event 6 in R2 and event 3 in R3) and HPV44 is a major parent in 2 events (event 7 and event 10 in R3).

To determine if the phylogenetic trees for different Gamma-PV genes were congruent, we used a more conclusive test, the SH test (Shimodaira and Hasegawa 1999) . The null hypothesis of the SH test states that the difference between trees (branch length, topology or likelihoods) is zero. The observed differences for the 405A1 and A2, 343A1 and A2 trees (deltaL values) are significantly greater than zero and the null hypothesis was rejected, thus declaring that these trees are significantly different i.e. incongruent (p < 0.05). The deltaL values of Alpha A1 and A2 trees were < 30 and hence not significantly greater than zero and we thus failed to reject the null hypothesis, thus declaring the Alpha A1 and A2 trees to be similar. The null hypothesis of the SH test states that the difference between trees (branch length, topology or likelihoods) is zero. The observed differences (deltaL values) were significantly greater than zero and rejected the null hypothesis and declare that the trees are significantly different i.e. incongruent (p < 0.05). Table 5 shows the results of the SH test using W-IQ-Tree, indicating that there is substantial phylogenetic incongruence between the 405A1 and A2, 343A1 and A2 trees but not for the Alpha A1 and A2 trees as shown by the p-values (p-SH).

After demonstrating five recombination events between human and non-human primate Alpha-PVs, we sought to show how the 82 Alpha-PV sequences from both humans and animals cluster together by constructing a phylogenetic tree as described below. Members of the alpha-12 species Table 5 Shimodaira-Hasegawa test for incongruence deltaL logL difference from the maximal logl in the set, bp-RELL bootstrap proportion using RELL method (Kishino et al. 1990) , p-KH p-value of one-sided (Kishino and Hasegawa 1989) , p-SH p-value of Shimodaira-Hasegawa test (Shimodaira and Hasegawa 1999) , p-WKH p-value of weighted KH test, p-WSH p-value of weighted SH test, c-ELW expected likelihood weight (Strimmer and Rambaut 2002) , p-AU p-value of approximately unbiased (AU) test (Shimodaira 2002) ; + 95% confidence sets, − significant exclusion (non-human primate PVs) are closely related to the members of the alpha-9 species (High-risk HPVs) and alpha-11 species as also shown in Fig. 1 . MmPV2 and MmPV6 (Macaca mulatta PVs from the Rhesus macaque monkeys) cluster with HPV54 of the alpha-13 species, while MfuPV2 and MmPV3 cluster with the alpha-2 HPVs. CgPV1 from the old world monkey, Colobus guereza, cluster with alpha-14 HPVs.

We report five novel recombination events between HPVs and non-human primate PVs using 3 different sequence sets. These recombination events were between High-Risk HPV Types and Macaca fascicularis PV1 (MfPV1), Macaca Fuscata PV2 (MfuPV2) and Pan Paniscus PV1 Fig. 1 Alphapapillomavirus genus. Phylogenetic analysis of all the 82 human and animal Alpha-PVs whole genomes. The sequences were aligned using MUSCLE v7.221 (Edgar 2004) . A maximumlikelihood tree of the nucleotide sequences was generated in PhyML using the optimal model of evolution (GTR+G) as determined within MEGA 7 (Kumar et al. 2016) . The Newick format of the tree was uploaded and modified in iTOL https ://itol.embl.de/tree/ (PpPV1) PVs. This observation provides the first evidence of interactions between PVs from different hosts and the likelihood of ancient host switching among Alpha-PVs, thus refuting, host specificity, a central dogma in PV evolution. Chen et al 2009b also postulated that an overlapping set of MfPVs in rhesus and cynomolgus macaques indicates non-human primate Alpha-PVs may not strictly be host specific (Chen et al. 2009b) . The same authors also report that members of the alpha-12 species (non-human primate PVs) are closely related to the members of the alpha-9 species (High-risk HPVs) and alpha-11 species as also shown in Fig. 1 . The most important biological prerequisite for recombination to occur is that two virus sequences be in the same space at the same time. So ancestral non-human primate PVs and HPVs should have been in the same space (biological niche) and time if recombination events had to occur as shown above. This unlikely event of cross-host infection by PVs challenges one of PVs central evolutionary tenets. Our findings provide the first likely evidence of such, but we are not the first authors to hypothesize this. Chen et al 2019 showed that non-human primate PVs share similar evolutionary histories and niche adaptation as the human counterparts (Chen et al. 2019) , using the example of HPV16 and MfPV3 that evolved from a most recent common ancestor containing the determinants of carcinogenicity. In another study, Chen et al 2018 also reported on the possibility of niche adaptation and viral transmission of HPVs from archaic hominins to modern humans, the archaic hominin-host-switch model .

It is important in the mathematical modelling of recombination to estimate the time when this recombination occurred. We predict, based on the time tree (https ://www. timet ree.org/about ) that the recombination between High-Risk HPV Types, Cynomolgus & Japanese Macaques and Bonobos PVs may have occurred not more than 76 million years ago (Smelov et al. 2018) , at the time of the last common ancestor of all primates. Chen et al. demonstrated specific host niche adaptation of primate PVs ensued thereafter by host co-divergence to have occurred at least 40MYA (Chen et al. , 2019 , which concurs with our prediction. Our observation suggests that more than 76MYA these viruses were in the same biological niche (same host). As the hosts evolved and diversified, the viruses adapted to specific host niches, which eventually led to coevolution with specific hosts, and this occurred before the speciation events of all primate host species including humans. Recombination events among Cetacean PVs detected in this study have also been previously reported elsewhere (Robles-Sikisaka et al. 2012), but were not the mainstay of this study.

The biological plausibility of Alpha-PVs recombination events is occasioned by the fact that recombination occurs only during viral replication (Pérez-Losada et al. 2015) . Additionally, the high replication rate of Alpha-PVs as seen from their oncogenicity and prevalence in mucosal sites of both humans and non-human primates promote recombination events (Chen et al. 2009b ). Multiple infections with HR-HPV types have been reported extensively in several studies (Vinodhini et al. 2012; Mbulawa et al. 2013; Murahwa et al. 2015; Teixeira et al. 2018) . These infections create a conducive environment for recombination to occur within the same biological niche or anatomical site through availing or enhancing DNA viral load, which in turn increases the probability of genomes to recombine. Multiple infections have also been reported in Betapapillomaviruses (Murahwa et al. 2015) and Gammapapillomaviruses (Meiring et al. 2017) cutaneous infections, but there are limited data on recombination in these genera.

A1 trees generated from the original sequence alignments and A2 trees generated after removing the recombinant regions from the original sequences were incongruent for the 405 PV sequences and the 343 concatenated sequences (i.e. the trees showed different phylogenies), hence supporting recombination as an essential driving force in PV evolution. However, the Alpha-PVs A1 and A2 trees showed congruence, seemingly rendering recombination an unimportant driving force in their evolution. The congruence can be explained by the low number (82) of Alpha-PV sequences currently available for the analysis to make inferential conclusions. Thus, the low number of sequences imply a lack of power for detection of phylogenetic incongruence. Moreover, the fact that essentially 7 sequences (PpPV1, MfPV7, HPV54, HPV70, HPV68, HPV39 and HPV66) of the 82 Alpha-PVs showed recombination events explains how these could not have changed the phylogeny of the whole Alpha-PV genus, with or without recombination regions from these 7 sequences. Thus, recombination plays a role in their evolution as shown in the dispersion of non-human-primate PVs across the Alpha-PV phylogenetic tree in Fig. 1 .

Our current knowledge of PVs is limited and focused on a few medically important and closely related human PVs associated with anogenital cancers and warts (Bravo and Felez-Sanchez 2015) , while the rest of the other plethora of PVs biology is largely understudied and unknown. Hence, assumptions made from studying these cannot be necessarily generalized and applied to all PVs.

Many PV sequences are still being discovered, and until a threshold number of representative sequences are attained, the PV community of researchers will remain underpowered to make assumptions closest to the reality of what happened in the evolution of this group of viruses. The addition of more PV sequences has a bearing on the understanding of the origin, evolution and clinical outcome prediction of given PV genomes. A recent discovery of the first fish PV, Sparus aurata PV1 from the gilthead bream fish showed the smallest PV genome consisting only of E1-E2-L1-L2 backbone (Lopez-Bueno et al. 2016) . It is prudent to hypothesize at this juncture that the PV ancestor is of marine origin unless more sequences become available for further analysis (Puustusmaa et al. 2017 ). More fish PV sequences are needed to elucidate our understanding of PV evolution.

Recombination, without doubt, constitutes an important driving force in Alpha-PVs evolution. We show that not more than 76MYA Alpha-PVs were in the same biological niche, a pre-requisite for recombination, and as the hosts evolved and diversified, the viruses adapted to specific host niches which eventually led to coevolution with specific hosts. Thus providing evidence that in ancient times no earlier than the Cretaceous period of the Mesozoic age, Alpha-PVs recombined and switched hosts, but whether this host switching is occurring currently is unknown. It is important to fully understand the evolutionary history of different PVs to better inform carcinogenicity and novel vaccine development.

Evidence of recombination within human alpha-papillomavirus

Classification of papillomaviruses (PVs) based on 189 PV types and proposal of taxonomic amendments

The clinical importance of understanding the evolution of papillomaviruses

Papillomaviruses: viral evolution, cancer and evolutionary medicine

Detecting recombination and diversifying selection in human alpha-papillomavirus

Evolutionary dynamics of variant genomes of human papillomavirus types 18, 45, and 97

Genomic diversity and interspecies host infection of alpha12 Macaca fascicularis papillomaviruses (MfPVs)

Niche adaptation and viral transmission of human papillomaviruses from archaic hominins to modern humans

Non-human primate papillomaviruses share similar evolutionary histories and niche adaptation as the human counterparts

High-level of viral genomic diversity in cervical cancers: a Brazilian study on human papillomavirus type 16

Classification of papillomaviruses

Intra-host sequence variability in human papillomavirus

MUSCLE: a multiple sequence alignment method with reduced time and space complexity

Sister-Scanning: a Monte Carlo procedure for assessing signals in recombinant sequences

Quantifying the phylodynamic forces driving papillomavirus evolution

ModelFinder: fast model selection for accurate phylogenetic estimates

Evaluation of the maximum likelihood estimate of the evolutionary tree topologies from DNA sequence data, and the branching order in hominoidea

Maximum likelihood inference of protein phylogeny and the origin of chloroplasts

The high fidelity and unique error signature of human DNA polymerase epsilon

MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets

The EMBL-EBI bioinformatics web and programmatic tools framework

Concurrence of iridovirus, polyomavirus, and a unique member of a new group of fish papillomaviruses in lymphocystis disease-affected Gilthead Sea Bream

RDP: detection of recombination amongst aligned sequences

Recombination in eukaryotic single stranded DNA viruses

A modified bootscan algorithm for automated identification of recombinant sequences and recombination breakpoints

RDP2: recombination detection and analysis from sequence alignments

The impact of human immunodeficiency virus on human papillomavirus transmission in heterosexually active couples

Analysis tool web services from the EMBL-EBI

High diversity of alpha, beta and gamma human papillomaviruses in genital samples from HIV-negative and HIV-positive heterosexual South African men

Ultrafast approximation for phylogenetic bootstrap

Towards quality and order in human papillomavirus research

Frequency of betapapillomavirus infections among HIV infected and uninfected Black Zimbabweans with cutaneous lesions

Possible emergence of new geminiviruses by frequent recombination

Recombination in viruses: mechanisms, methods of study, and evolutionary consequences

Evaluation of methods for detecting recombination from DNA sequences: computer simulations

The enigmatic origin of papillomavirus protein domains

Evidence of recombination and positive selection in cetacean papillomaviruses

An approximately unbiased test of phylogenetic tree selection

Multiple comparisons of loglikelihoods with applications to phylogenetic inference

Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega

Why do RNA viruses recombine?

Beta and gamma human papillomaviruses in anal and genital sites among men: prevalence and determinants

Analyzing the mosaic structure of genes

Inferring confidence sets of possibly misspecified gene trees

High risk human papillomavirus prevalence and genotype distribution among women infected with HIV in Manaus

W-IQ-TREE: a fast online phylogenetic tool for maximum likelihood analysis

ICTV virus taxonomy profile: papillomaviridae

The launch of an international animal papillomavirus reference center

Virus Taxonomy. Seventh Report of the International Committee on Taxonomy of Viruses

Martin DP (2006) Evidence of ancient papillomavirus recombination

Prevalence and risk factors of HPV infection among women from various provinces of the world

Acknowledgements We thank Prof Darren Martin for his help rendered with RDP4 software. Dr Alltalents T. Murahwa is the recipient of a Post Doc fellowship from the South African National Research Foundation. This work is based upon research supported by the South African Research Chairs Initiative of the Department of Science and Technology of South Africa and NRF.

Conflict of interest The authors declare that they have no competing interests.