Unique genome organization of non-mammalian papillomaviruses provides insights into the evolution of viral early proteins
Identifieur interne : 002D41 ( Pmc/Corpus ); précédent : 002D40; suivant : 002D42Unique genome organization of non-mammalian papillomaviruses provides insights into the evolution of viral early proteins
Auteurs : Koenraad Van Doorslaer ; Valeria Ruoppolo ; Annie Schmidt ; Amelie Lescroël ; Dennis Jongsomjit ; Megan Elrod ; Simona Kraberger ; Daisy Stainton ; Katie M. Dugger ; Grant Ballard ; David G. Ainley ; Arvind VarsaniSource :
- Virus Evolution [ 2057-1577 ] ; 2017.
Abstract
The family
Url:
DOI: 10.1093/ve/vex027
PubMed: 29026649
PubMed Central: 5632515
Links to Exploration step
PMC:5632515Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Unique genome organization of non-mammalian papillomaviruses provides insights into the evolution of viral early proteins</title><author><name sortKey="Van Doorslaer, Koenraad" sort="Van Doorslaer, Koenraad" uniqKey="Van Doorslaer K" first="Koenraad" last="Van Doorslaer">Koenraad Van Doorslaer</name><affiliation><nlm:aff id="vex027-aff1">ACBS and Bio5, University of Arizona, 1657 E Helen St, Tucson, AZ 85721, USA</nlm:aff></affiliation></author><author><name sortKey="Ruoppolo, Valeria" sort="Ruoppolo, Valeria" uniqKey="Ruoppolo V" first="Valeria" last="Ruoppolo">Valeria Ruoppolo</name><affiliation><nlm:aff id="vex027-aff2">Laboratório de Patologia Comparada de Animais Selvagens (LAPCOM), Departamento de Patologia, Faculdade de Medicina Veterinária e Zootecnia, Universidade de São Paulo, São Paulo, Brazil</nlm:aff></affiliation></author><author><name sortKey="Schmidt, Annie" sort="Schmidt, Annie" uniqKey="Schmidt A" first="Annie" last="Schmidt">Annie Schmidt</name><affiliation><nlm:aff id="vex027-aff3">Point Blue Conservation Science, Petaluma, CA 94954, USA</nlm:aff></affiliation></author><author><name sortKey="Lescroel, Amelie" sort="Lescroel, Amelie" uniqKey="Lescroel A" first="Amelie" last="Lescroël">Amelie Lescroël</name><affiliation><nlm:aff id="vex027-aff3">Point Blue Conservation Science, Petaluma, CA 94954, USA</nlm:aff></affiliation><affiliation><nlm:aff id="vex027-aff4">Centre d'Ecologie Fonctionnelle et Evolutive – CNRS, UMR 5175, Montpellier, France</nlm:aff></affiliation></author><author><name sortKey="Jongsomjit, Dennis" sort="Jongsomjit, Dennis" uniqKey="Jongsomjit D" first="Dennis" last="Jongsomjit">Dennis Jongsomjit</name><affiliation><nlm:aff id="vex027-aff3">Point Blue Conservation Science, Petaluma, CA 94954, USA</nlm:aff></affiliation></author><author><name sortKey="Elrod, Megan" sort="Elrod, Megan" uniqKey="Elrod M" first="Megan" last="Elrod">Megan Elrod</name><affiliation><nlm:aff id="vex027-aff3">Point Blue Conservation Science, Petaluma, CA 94954, USA</nlm:aff></affiliation></author><author><name sortKey="Kraberger, Simona" sort="Kraberger, Simona" uniqKey="Kraberger S" first="Simona" last="Kraberger">Simona Kraberger</name><affiliation><nlm:aff id="vex027-aff5">The Biodesign Center for Fundamental and Applied Microbiomics, Center for Evolution and Medicine, School of Life sciences, Arizona State University, Tempe, AZ 85287, USA</nlm:aff></affiliation></author><author><name sortKey="Stainton, Daisy" sort="Stainton, Daisy" uniqKey="Stainton D" first="Daisy" last="Stainton">Daisy Stainton</name><affiliation><nlm:aff id="vex027-aff6">School of Biological Sciences, The University of Queensland, Brisbane, QLD, Australia</nlm:aff></affiliation></author><author><name sortKey="Dugger, Katie M" sort="Dugger, Katie M" uniqKey="Dugger K" first="Katie M" last="Dugger">Katie M. Dugger</name><affiliation><nlm:aff id="vex027-aff7">U.S. Geological Survey, Oregon Cooperative Fish and Wildlife Research Unit, Department of Fisheries and Wildlife, Oregon State University, Corvallis, OR 97331, USA</nlm:aff></affiliation></author><author><name sortKey="Ballard, Grant" sort="Ballard, Grant" uniqKey="Ballard G" first="Grant" last="Ballard">Grant Ballard</name><affiliation><nlm:aff id="vex027-aff3">Point Blue Conservation Science, Petaluma, CA 94954, USA</nlm:aff></affiliation></author><author><name sortKey="Ainley, David G" sort="Ainley, David G" uniqKey="Ainley D" first="David G" last="Ainley">David G. Ainley</name><affiliation><nlm:aff id="vex027-aff8">HT Harvey and Associates, Los Gatos, CA 95032, USA</nlm:aff></affiliation></author><author><name sortKey="Varsani, Arvind" sort="Varsani, Arvind" uniqKey="Varsani A" first="Arvind" last="Varsani">Arvind Varsani</name><affiliation><nlm:aff id="vex027-aff5">The Biodesign Center for Fundamental and Applied Microbiomics, Center for Evolution and Medicine, School of Life sciences, Arizona State University, Tempe, AZ 85287, USA</nlm:aff></affiliation><affiliation><nlm:aff id="vex027-aff9">Structural Biology Research Unit, Department of Clinical Laboratory Sciences, University of Cape Town, Observatory, Cape Town, 7925, South Africa</nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">29026649</idno><idno type="pmc">5632515</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC5632515</idno><idno type="RBID">PMC:5632515</idno><idno type="doi">10.1093/ve/vex027</idno><date when="2017">2017</date><idno type="wicri:Area/Pmc/Corpus">002D41</idno><idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">002D41</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Unique genome organization of non-mammalian papillomaviruses provides insights into the evolution of viral early proteins</title><author><name sortKey="Van Doorslaer, Koenraad" sort="Van Doorslaer, Koenraad" uniqKey="Van Doorslaer K" first="Koenraad" last="Van Doorslaer">Koenraad Van Doorslaer</name><affiliation><nlm:aff id="vex027-aff1">ACBS and Bio5, University of Arizona, 1657 E Helen St, Tucson, AZ 85721, USA</nlm:aff></affiliation></author><author><name sortKey="Ruoppolo, Valeria" sort="Ruoppolo, Valeria" uniqKey="Ruoppolo V" first="Valeria" last="Ruoppolo">Valeria Ruoppolo</name><affiliation><nlm:aff id="vex027-aff2">Laboratório de Patologia Comparada de Animais Selvagens (LAPCOM), Departamento de Patologia, Faculdade de Medicina Veterinária e Zootecnia, Universidade de São Paulo, São Paulo, Brazil</nlm:aff></affiliation></author><author><name sortKey="Schmidt, Annie" sort="Schmidt, Annie" uniqKey="Schmidt A" first="Annie" last="Schmidt">Annie Schmidt</name><affiliation><nlm:aff id="vex027-aff3">Point Blue Conservation Science, Petaluma, CA 94954, USA</nlm:aff></affiliation></author><author><name sortKey="Lescroel, Amelie" sort="Lescroel, Amelie" uniqKey="Lescroel A" first="Amelie" last="Lescroël">Amelie Lescroël</name><affiliation><nlm:aff id="vex027-aff3">Point Blue Conservation Science, Petaluma, CA 94954, USA</nlm:aff></affiliation><affiliation><nlm:aff id="vex027-aff4">Centre d'Ecologie Fonctionnelle et Evolutive – CNRS, UMR 5175, Montpellier, France</nlm:aff></affiliation></author><author><name sortKey="Jongsomjit, Dennis" sort="Jongsomjit, Dennis" uniqKey="Jongsomjit D" first="Dennis" last="Jongsomjit">Dennis Jongsomjit</name><affiliation><nlm:aff id="vex027-aff3">Point Blue Conservation Science, Petaluma, CA 94954, USA</nlm:aff></affiliation></author><author><name sortKey="Elrod, Megan" sort="Elrod, Megan" uniqKey="Elrod M" first="Megan" last="Elrod">Megan Elrod</name><affiliation><nlm:aff id="vex027-aff3">Point Blue Conservation Science, Petaluma, CA 94954, USA</nlm:aff></affiliation></author><author><name sortKey="Kraberger, Simona" sort="Kraberger, Simona" uniqKey="Kraberger S" first="Simona" last="Kraberger">Simona Kraberger</name><affiliation><nlm:aff id="vex027-aff5">The Biodesign Center for Fundamental and Applied Microbiomics, Center for Evolution and Medicine, School of Life sciences, Arizona State University, Tempe, AZ 85287, USA</nlm:aff></affiliation></author><author><name sortKey="Stainton, Daisy" sort="Stainton, Daisy" uniqKey="Stainton D" first="Daisy" last="Stainton">Daisy Stainton</name><affiliation><nlm:aff id="vex027-aff6">School of Biological Sciences, The University of Queensland, Brisbane, QLD, Australia</nlm:aff></affiliation></author><author><name sortKey="Dugger, Katie M" sort="Dugger, Katie M" uniqKey="Dugger K" first="Katie M" last="Dugger">Katie M. Dugger</name><affiliation><nlm:aff id="vex027-aff7">U.S. Geological Survey, Oregon Cooperative Fish and Wildlife Research Unit, Department of Fisheries and Wildlife, Oregon State University, Corvallis, OR 97331, USA</nlm:aff></affiliation></author><author><name sortKey="Ballard, Grant" sort="Ballard, Grant" uniqKey="Ballard G" first="Grant" last="Ballard">Grant Ballard</name><affiliation><nlm:aff id="vex027-aff3">Point Blue Conservation Science, Petaluma, CA 94954, USA</nlm:aff></affiliation></author><author><name sortKey="Ainley, David G" sort="Ainley, David G" uniqKey="Ainley D" first="David G" last="Ainley">David G. Ainley</name><affiliation><nlm:aff id="vex027-aff8">HT Harvey and Associates, Los Gatos, CA 95032, USA</nlm:aff></affiliation></author><author><name sortKey="Varsani, Arvind" sort="Varsani, Arvind" uniqKey="Varsani A" first="Arvind" last="Varsani">Arvind Varsani</name><affiliation><nlm:aff id="vex027-aff5">The Biodesign Center for Fundamental and Applied Microbiomics, Center for Evolution and Medicine, School of Life sciences, Arizona State University, Tempe, AZ 85287, USA</nlm:aff></affiliation><affiliation><nlm:aff id="vex027-aff9">Structural Biology Research Unit, Department of Clinical Laboratory Sciences, University of Cape Town, Observatory, Cape Town, 7925, South Africa</nlm:aff></affiliation></author></analytic><series><title level="j">Virus Evolution</title><idno type="eISSN">2057-1577</idno><imprint><date when="2017">2017</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><title>Abstract</title><p>The family <italic>Papillomaviridae</italic> contains more than 320 papillomavirus types, with most having been identified as infecting skin and mucosal epithelium in mammalian hosts. To date, only nine non-mammalian papillomaviruses have been described from birds (<italic>n</italic> = 5), a fish (<italic>n</italic> = 1), a snake (<italic>n</italic> = 1), and turtles (<italic>n</italic> = 2). The identification of papillomaviruses in sauropsids and a sparid fish suggests that early ancestors of papillomaviruses were already infecting the earliest Euteleostomi. The Euteleostomi clade includes more than 90 per cent of the living vertebrate species, and progeny virus could have been passed on to all members of this clade, inhabiting virtually every habitat on the planet. As part of this study, we isolated a novel papillomavirus from a 16-year-old female Adélie penguin (<italic>Pygoscelis adeliae</italic>) from Cape Crozier, Ross Island (Antarctica). The new papillomavirus shares ∼64 per cent genome-wide identity to a previously described Adélie penguin papillomavirus. Phylogenetic analyses show that the non-mammalian viruses (expect the python, <italic>Morelia spilota</italic>, associated papillomavirus) cluster near the base of the papillomavirus evolutionary tree. A papillomavirus isolated from an avian host (Northern fulmar; <italic>Fulmarus glacialis</italic>), like the two turtle papillomaviruses, lacks a putative E9 protein that is found in all other avian papillomaviruses. Furthermore, the Northern fulmar papillomavirus has an E7 more similar to the mammalian viruses than the other avian papillomaviruses. Typical E6 proteins of mammalian papillomaviruses have two Zinc finger motifs, whereas the sauropsid papillomaviruses only have one such motif. Furthermore, this motif is absent in the fish papillomavirus. Thus, it is highly likely that the most recent common ancestor of the mammalian and sauropsid papillomaviruses had a single motif E6. It appears that a motif duplication resulted in mammalian papillomaviruses having a double Zinc finger motif in E6. We estimated the divergence time between Northern fulmar-associated papillomavirus and the other Sauropsid papillomaviruses be to around 250 million years ago, during the Paleozoic-Mesozoic transition and our analysis dates the root of the papillomavirus tree between 400 and 600 million years ago. Our analysis shows evidence for niche adaptation and that these non-mammalian viruses have highly divergent E6 and E7 proteins, providing insights into the evolution of the early viral (onco-)proteins.</p></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Aiewsakun, P" uniqKey="Aiewsakun P">P. Aiewsakun</name></author><author><name sortKey="Katzourakis, A" uniqKey="Katzourakis A">A. Katzourakis</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Aiewsakun, P" uniqKey="Aiewsakun P">P. Aiewsakun</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bell, M A" uniqKey="Bell M">M. A. Bell</name></author><author><name sortKey="Lloyd, G" uniqKey="Lloyd G">G. Lloyd</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bergvall, M" uniqKey="Bergvall M">M. Bergvall</name></author><author><name sortKey="Melendy, T" uniqKey="Melendy T">T. Melendy</name></author><author><name sortKey="Archambault, J" uniqKey="Archambault J">J. Archambault</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bernard, H U" uniqKey="Bernard H">H. U. Bernard</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bernard, H U" uniqKey="Bernard H">H. U. Bernard</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bouckaert, R" uniqKey="Bouckaert R">R. Bouckaert</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Brooks, D R" uniqKey="Brooks D">D. R. Brooks</name></author><author><name sortKey="Ferrao, A L" uniqKey="Ferrao A">A. L. Ferrao</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Buck, C B" uniqKey="Buck C">C. B. Buck</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Dos Reis, M" uniqKey="Dos Reis M">M. dos Reis</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Drummond, A J" uniqKey="Drummond A">A. J. Drummond</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Duchene, S" uniqKey="Duchene S">S. Duchene</name></author><author><name sortKey="Holmes, E C" uniqKey="Holmes E">E. C. Holmes</name></author><author><name sortKey="Ho, S Y" uniqKey="Ho S">S. Y. Ho</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Duchene, S" uniqKey="Duchene S">S. Duchene</name></author><author><name sortKey="Lanfear, R" uniqKey="Lanfear R">R. Lanfear</name></author><author><name sortKey="Ho, S Y" uniqKey="Ho S">S. Y. Ho</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gaynor, A M" uniqKey="Gaynor A">A. M. Gaynor</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gottschling, M" uniqKey="Gottschling M">M. Gottschling</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Guindon, S" uniqKey="Guindon S">S. Guindon</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hedges, S B" uniqKey="Hedges S">S. B. Hedges</name></author><author><name sortKey="Dudley, J" uniqKey="Dudley J">J. Dudley</name></author><author><name sortKey="Kumar, S" uniqKey="Kumar S">S. Kumar</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hedges, S B" uniqKey="Hedges S">S. B. Hedges</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Heled, J" uniqKey="Heled J">J. Heled</name></author><author><name sortKey="Drummond, A J" uniqKey="Drummond A">A. J. Drummond</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Herbst, L H" uniqKey="Herbst L">L. H. Herbst</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Huelsenbeck, J P" uniqKey="Huelsenbeck J">J. P. Huelsenbeck</name></author><author><name sortKey="Ronquist, F" uniqKey="Ronquist F">F. Ronquist</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jeckel, S" uniqKey="Jeckel S">S. Jeckel</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Katoh, K" uniqKey="Katoh K">K. Katoh</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Katoh, K" uniqKey="Katoh K">K. Katoh</name></author><author><name sortKey="Standley, D M" uniqKey="Standley D">D. M. Standley</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kearse, M" uniqKey="Kearse M">M. Kearse</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lanfear, R" uniqKey="Lanfear R">R. Lanfear</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lange, C E" uniqKey="Lange C">C. E. Lange</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Le, S Q" uniqKey="Le S">S. Q. Le</name></author><author><name sortKey="Gascuel, O" uniqKey="Gascuel O">O. Gascuel</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lopez Bueno, A" uniqKey="Lopez Bueno A">A. Lopez-Bueno</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Miller, M A" uniqKey="Miller M">M. A. Miller</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Morin, G" uniqKey="Morin G">G. Morin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Muhire, B M" uniqKey="Muhire B">B. M. Muhire</name></author><author><name sortKey="Varsani, A" uniqKey="Varsani A">A. Varsani</name></author><author><name sortKey="Martin, D P" uniqKey="Martin D">D. P. Martin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nylander, J A" uniqKey="Nylander J">J. A. Nylander</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Paradis, E" uniqKey="Paradis E">E. Paradis</name></author><author><name sortKey="Claude, J" uniqKey="Claude J">J. Claude</name></author><author><name sortKey="Strimmer, K" uniqKey="Strimmer K">K. Strimmer</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Pattengale, N D" uniqKey="Pattengale N">N. D. Pattengale</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Prum, R O" uniqKey="Prum R">R. O. Prum</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Rambaut, A" uniqKey="Rambaut A">A. Rambaut</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Rector, A" uniqKey="Rector A">A. Rector</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Rector, A" uniqKey="Rector A">A. Rector</name></author><author><name sortKey="Van Ranst, M" uniqKey="Van Ranst M">M. Van Ranst</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Rey, K" uniqKey="Rey K">K. Rey</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ripley, B D" uniqKey="Ripley B">B. D. Ripley</name></author><author><name sortKey="Thompson, M" uniqKey="Thompson M">M. Thompson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Roman, A" uniqKey="Roman A">A. Roman</name></author><author><name sortKey="Munger, K" uniqKey="Munger K">K. Munger</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ronquist, F" uniqKey="Ronquist F">F. Ronquist</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ronquist, F" uniqKey="Ronquist F">F. Ronquist</name></author><author><name sortKey="Huelsenbeck, J P" uniqKey="Huelsenbeck J">J. P. Huelsenbeck</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ruboyianes, R" uniqKey="Ruboyianes R">R. Ruboyianes</name></author><author><name sortKey="Worobey, M" uniqKey="Worobey M">M. Worobey</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Schmidt, H A" uniqKey="Schmidt H">H. A. Schmidt</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shah, S D" uniqKey="Shah S">S. D. Shah</name></author><author><name sortKey="Doorbar, J" uniqKey="Doorbar J">J. Doorbar</name></author><author><name sortKey="Goldstein, R A" uniqKey="Goldstein R">R. A. Goldstein</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shimodaira, H" uniqKey="Shimodaira H">H. Shimodaira</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shimodaira, H" uniqKey="Shimodaira H">H. Shimodaira</name></author><author><name sortKey="Hasegawa, M" uniqKey="Hasegawa M">M. Hasegawa</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Stamatakis, A" uniqKey="Stamatakis A">A. Stamatakis</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tachezy, R" uniqKey="Tachezy R">R. Tachezy</name></author><author><name sortKey="Rector, A" uniqKey="Rector A">A. Rector</name></author><author><name sortKey="Havelkova, M" uniqKey="Havelkova M">M. Havelkova</name></author><author><name sortKey="Wollants, E" uniqKey="Wollants E">E. Wollants</name></author><author><name sortKey="Fiten, P" uniqKey="Fiten P">P. Fiten</name></author><author><name sortKey="Opdenakker, G" uniqKey="Opdenakker G">G. Opdenakker</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Terai, M" uniqKey="Terai M">M. Terai</name></author><author><name sortKey="Desalle, R" uniqKey="Desalle R">R. DeSalle</name></author><author><name sortKey="Burk, R D" uniqKey="Burk R">R. D. Burk</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Van Doorslaer, K" uniqKey="Van Doorslaer K">K. Van Doorslaer</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Van Doorslaer, K" uniqKey="Van Doorslaer K">K. Van Doorslaer</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Van Doorslaer, K" uniqKey="Van Doorslaer K">K. Van Doorslaer</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Van Doorslaer, K" uniqKey="Van Doorslaer K">K. Van Doorslaer</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Vande Pol, S B" uniqKey="Vande Pol S">S. B. Vande Pol</name></author><author><name sortKey="Klingelhutz, A J" uniqKey="Klingelhutz A">A. J. Klingelhutz</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Varsani, A" uniqKey="Varsani A">A. Varsani</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Varsani, A" uniqKey="Varsani A">A. Varsani</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zanier, K" uniqKey="Zanier K">K. Zanier</name></author></analytic></biblStruct></listBibl></div1></back></TEI><pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">Virus Evol</journal-id><journal-id journal-id-type="iso-abbrev">Virus Evol</journal-id><journal-id journal-id-type="publisher-id">vevolu</journal-id><journal-title-group><journal-title>Virus Evolution</journal-title></journal-title-group><issn pub-type="epub">2057-1577</issn><publisher><publisher-name>Oxford University Press</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">29026649</article-id><article-id pub-id-type="pmc">5632515</article-id><article-id pub-id-type="doi">10.1093/ve/vex027</article-id><article-id pub-id-type="publisher-id">vex027</article-id><article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject></subj-group></article-categories><title-group><article-title>Unique genome organization of non-mammalian papillomaviruses provides insights into the evolution of viral early proteins</article-title></title-group><contrib-group><contrib contrib-type="author" corresp="yes"><contrib-id contrib-id-type="orcid" authenticated="false">http://orcid.org/0000-0002-2985-0733</contrib-id><name><surname>Van Doorslaer</surname><given-names>Koenraad</given-names></name><xref ref-type="aff" rid="vex027-aff1">1</xref><xref ref-type="corresp" rid="vex027-cor1"></xref><xref ref-type="author-notes" rid="vex027-FM1"></xref><pmc-comment>vandoorslaer@email.arizona.edu </pmc-comment>
</contrib><contrib contrib-type="author"><name><surname>Ruoppolo</surname><given-names>Valeria</given-names></name><xref ref-type="aff" rid="vex027-aff2">2</xref></contrib><contrib contrib-type="author"><name><surname>Schmidt</surname><given-names>Annie</given-names></name><xref ref-type="aff" rid="vex027-aff3">3</xref></contrib><contrib contrib-type="author"><name><surname>Lescroël</surname><given-names>Amelie</given-names></name><xref ref-type="aff" rid="vex027-aff3">3</xref><xref ref-type="aff" rid="vex027-aff4">4</xref></contrib><contrib contrib-type="author"><name><surname>Jongsomjit</surname><given-names>Dennis</given-names></name><xref ref-type="aff" rid="vex027-aff3">3</xref></contrib><contrib contrib-type="author"><name><surname>Elrod</surname><given-names>Megan</given-names></name><xref ref-type="aff" rid="vex027-aff3">3</xref></contrib><contrib contrib-type="author"><name><surname>Kraberger</surname><given-names>Simona</given-names></name><xref ref-type="aff" rid="vex027-aff5">5</xref></contrib><contrib contrib-type="author"><name><surname>Stainton</surname><given-names>Daisy</given-names></name><xref ref-type="aff" rid="vex027-aff6">6</xref></contrib><contrib contrib-type="author"><name><surname>Dugger</surname><given-names>Katie M</given-names></name><xref ref-type="aff" rid="vex027-aff7">7</xref></contrib><contrib contrib-type="author"><name><surname>Ballard</surname><given-names>Grant</given-names></name><xref ref-type="aff" rid="vex027-aff3">3</xref></contrib><contrib contrib-type="author"><name><surname>Ainley</surname><given-names>David G</given-names></name><xref ref-type="aff" rid="vex027-aff8">8</xref></contrib><contrib contrib-type="author" corresp="yes"><contrib-id contrib-id-type="orcid" authenticated="false">http://orcid.org/0000-0003-4111-2415</contrib-id><name><surname>Varsani</surname><given-names>Arvind</given-names></name><xref ref-type="aff" rid="vex027-aff5">5</xref><xref ref-type="aff" rid="vex027-aff9">9</xref><xref ref-type="corresp" rid="vex027-cor1"></xref><xref ref-type="author-notes" rid="vex027-FM2"></xref><pmc-comment>vandoorslaer@email.arizona.edu </pmc-comment>
</contrib></contrib-group><aff id="vex027-aff1"><label>1</label>ACBS and Bio5, University of Arizona, 1657 E Helen St, Tucson, AZ 85721, USA</aff><aff id="vex027-aff2"><label>2</label>Laboratório de Patologia Comparada de Animais Selvagens (LAPCOM), Departamento de Patologia, Faculdade de Medicina Veterinária e Zootecnia, Universidade de São Paulo, São Paulo, Brazil</aff><aff id="vex027-aff3"><label>3</label>Point Blue Conservation Science, Petaluma, CA 94954, USA</aff><aff id="vex027-aff4"><label>4</label>Centre d'Ecologie Fonctionnelle et Evolutive – CNRS, UMR 5175, Montpellier, France</aff><aff id="vex027-aff5"><label>5</label>The Biodesign Center for Fundamental and Applied Microbiomics, Center for Evolution and Medicine, School of Life sciences, Arizona State University, Tempe, AZ 85287, USA</aff><aff id="vex027-aff6"><label>6</label>School of Biological Sciences, The University of Queensland, Brisbane, QLD, Australia</aff><aff id="vex027-aff7"><label>7</label>U.S. Geological Survey, Oregon Cooperative Fish and Wildlife Research Unit, Department of Fisheries and Wildlife, Oregon State University, Corvallis, OR 97331, USA</aff><aff id="vex027-aff8"><label>8</label>HT Harvey and Associates, Los Gatos, CA 95032, USA</aff><aff id="vex027-aff9"><label>9</label>Structural Biology Research Unit, Department of Clinical Laboratory Sciences, University of Cape Town, Observatory, Cape Town, 7925, South Africa</aff><author-notes><corresp id="vex027-cor1"><label>*</label>Corresponding authors: E-mails: <email>vandoorslaer@email.arizona.edu</email> and <email>arvind.varsani@asu.edu</email></corresp><fn id="vex027-FM1"><label>†</label><p><ext-link ext-link-type="uri" xlink:href="http://orcid.org/0000-0002-2985-0733">http://orcid.org/0000-0002-2985-0733</ext-link></p></fn><fn id="vex027-FM2"><label>‡</label><p><ext-link ext-link-type="uri" xlink:href="http://orcid.org/0000-0003-4111-2415">http://orcid.org/0000-0003-4111-2415</ext-link></p></fn></author-notes><pub-date pub-type="collection"><month>7</month><year>2017</year></pub-date><pub-date pub-type="epub" iso-8601-date="2017-10-06"><day>06</day><month>10</month><year>2017</year></pub-date><pub-date pub-type="pmc-release"><day>06</day><month>10</month><year>2017</year></pub-date><pmc-comment> PMC Release delay is 0 months and 0 days and was based on the
<volume>3</volume><issue>2</issue><elocation-id>vex027</elocation-id><permissions><copyright-statement>© The Author 2017. Published by Oxford University Press.</copyright-statement><copyright-year>2017</copyright-year><license xlink:href="http://creativecommons.org/licenses/by/4.0/" license-type="cc-by"><license-p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">http://creativecommons.org/licenses/by/4.0/</ext-link>), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.</license-p></license></permissions><self-uri xlink:href="vex027.pdf"></self-uri><abstract><title>Abstract</title><p>The family <italic>Papillomaviridae</italic> contains more than 320 papillomavirus types, with most having been identified as infecting skin and mucosal epithelium in mammalian hosts. To date, only nine non-mammalian papillomaviruses have been described from birds (<italic>n</italic> = 5), a fish (<italic>n</italic> = 1), a snake (<italic>n</italic> = 1), and turtles (<italic>n</italic> = 2). The identification of papillomaviruses in sauropsids and a sparid fish suggests that early ancestors of papillomaviruses were already infecting the earliest Euteleostomi. The Euteleostomi clade includes more than 90 per cent of the living vertebrate species, and progeny virus could have been passed on to all members of this clade, inhabiting virtually every habitat on the planet. As part of this study, we isolated a novel papillomavirus from a 16-year-old female Adélie penguin (<italic>Pygoscelis adeliae</italic>) from Cape Crozier, Ross Island (Antarctica). The new papillomavirus shares ∼64 per cent genome-wide identity to a previously described Adélie penguin papillomavirus. Phylogenetic analyses show that the non-mammalian viruses (expect the python, <italic>Morelia spilota</italic>, associated papillomavirus) cluster near the base of the papillomavirus evolutionary tree. A papillomavirus isolated from an avian host (Northern fulmar; <italic>Fulmarus glacialis</italic>), like the two turtle papillomaviruses, lacks a putative E9 protein that is found in all other avian papillomaviruses. Furthermore, the Northern fulmar papillomavirus has an E7 more similar to the mammalian viruses than the other avian papillomaviruses. Typical E6 proteins of mammalian papillomaviruses have two Zinc finger motifs, whereas the sauropsid papillomaviruses only have one such motif. Furthermore, this motif is absent in the fish papillomavirus. Thus, it is highly likely that the most recent common ancestor of the mammalian and sauropsid papillomaviruses had a single motif E6. It appears that a motif duplication resulted in mammalian papillomaviruses having a double Zinc finger motif in E6. We estimated the divergence time between Northern fulmar-associated papillomavirus and the other Sauropsid papillomaviruses be to around 250 million years ago, during the Paleozoic-Mesozoic transition and our analysis dates the root of the papillomavirus tree between 400 and 600 million years ago. Our analysis shows evidence for niche adaptation and that these non-mammalian viruses have highly divergent E6 and E7 proteins, providing insights into the evolution of the early viral (onco-)proteins.</p></abstract><kwd-group><kwd>papillomavirus</kwd><kwd>evolution</kwd><kwd>bird</kwd><kwd>reptile</kwd><kwd>avian</kwd></kwd-group><funding-group><award-group award-type="grant"><funding-source><named-content content-type="funder-name">Cooperative State Research, Education, and Extension Service</named-content><named-content content-type="funder-identifier">10.13039/100007014</named-content></funding-source><award-id>NC229</award-id></award-group><award-group award-type="grant"><funding-source><named-content content-type="funder-name">National Science Foundation</named-content><named-content content-type="funder-identifier">10.13039/100000001</named-content></funding-source><award-id>ANT-0944411</award-id></award-group></funding-group><counts><page-count count="12"></page-count></counts></article-meta></front><body><sec><title>1. Introduction</title><p>The family of the <italic>Papillomaviridae</italic> contains more than 320 distinct viral species. In recent years, there has been an exponential increase in the identification of non-human-associated papillomaviruses. The majority of known papillomaviruses have been isolated from humans and other mammalian hosts. To date, nine non-mammalian viruses have been described. These viruses infect birds including, yellow-necked francolin (<italic>Francolinus leucoscepus</italic>), common chaffinch (<italic>Fringilla coelebs</italic>), northern fulmar (<italic>Fulmar glacialis</italic>), African Grey parrot (<italic>Psittacus erithacus</italic>) and Adélie penguin (<italic>Pygoscelis adeliae</italic>) (<xref rid="vex027-B53" ref-type="bibr">Terai et al. 2002</xref>; <xref rid="vex027-B51" ref-type="bibr">Tachezy et al. 2002</xref>; <xref rid="vex027-B55" ref-type="bibr">Van Doorslaer et al. 2009</xref>; <xref rid="vex027-B59" ref-type="bibr">Varsani et al. 2014</xref>; <xref rid="vex027-B14" ref-type="bibr">Gaynor et al. 2015</xref>; <xref rid="vex027-B60" ref-type="bibr">Varsani et al. 2015</xref> ), as well as green sea turtle (<italic>Chelonia mydas</italic>), loggerhead sea turtle (<italic>Caretta caretta</italic>) (<xref rid="vex027-B20" ref-type="bibr">Herbst et al. 2009</xref>), Carpet python (<italic>Morelia spilota</italic>) (<xref rid="vex027-B27" ref-type="bibr">Lange et al. 2011</xref>), and gilt-head bream fish (<italic>Sparus aurata</italic>) (<xref rid="vex027-B29" ref-type="bibr">Lopez-Bueno et al. 2016</xref>).</p><p>Papillomaviruses had been thought to be restricted to amniotes (<xref rid="vex027-B39" ref-type="bibr">Rector and Van Ranst 2013</xref>). However, recent studies that have identified papillomaviruses in fish (<xref rid="vex027-B29" ref-type="bibr">Lopez-Bueno et al. 2016</xref>) reveal that early ancestors of the papillomaviruses were already infecting the earliest members of Euteleostomi. That clade includes more than 90 per cent of the living vertebrate host species, and progeny viruses could have been passed on to all members of this clade, inhabiting virtually every habitat on the planet. Also, the hosts that have been extensively studied (humans, cattle, dogs, etc.), have an extensive repertoire of highly species-specific papillomaviruses. The observation that papillomaviruses cause benign infections unable to cross the hosts' species-barrier has led to the hypothesis that papillomaviruses coevolved alongside their hosts (<xref rid="vex027-B5" ref-type="bibr">Bernard 1994</xref>). However, the papillomavirus phylogeny shows evidence of initial, ancestral niche sorting events followed by virus–host linked speciation (<xref rid="vex027-B47" ref-type="bibr">Shah et al. 2010</xref>; <xref rid="vex027-B15" ref-type="bibr">Gottschling et al. 2011</xref>; <xref rid="vex027-B54" ref-type="bibr">Van Doorslaer 2013</xref>). This suggests that the acquisition of new ecological niches on the host (gain/loss of fur, the evolution of sweat glands, etc.), not host speciation per se may be driving papillomaviruses evolution (<xref rid="vex027-B8" ref-type="bibr">Brooks and Ferrao 2005</xref>).</p><p>As part of this study, we report the identification of a novel papillomavirus from an Adélie penguin from Cape Crozier, Ross Island (Antarctica). We show that this novel virus is most closely related to a previously identified Adélie penguin papillomavirus and that together they cluster with the other non-mammalian viruses near the base of the papillomavirus evolutionary tree. Comprehensive analysis of the ten non-mammalian viral genomes provides evidence for niche adaptation followed by co-evolution of these papillomaviruses. Furthermore, these non-mammalian viruses have highly divergent E6 and E7 proteins, consideration of which allowed for a deeper understanding of the evolutionary history of these papillomaviruses. The E6 and E7 protein’s unique features may enable a better understanding of the evolution of these oncoproteins.</p></sec><sec><title>2. Materials and methods</title><sec><title>2.1 Ethics statement</title><p>The samples were collected under Oregon State University’s Institutional Animal Care and Use Committee approved Animal Care and Use Protocol (ACUP permit #s 4130 and 4622) Antarctic Conservation Act Permit #2011-002 M#7 from the National Foundation of Science (USA) through H.T. Harvey & Associates.</p></sec><sec><title>2.2 Sample preparation, DNA enrichment, and sequencing</title><p>The study was conducted at the Cape Crozier Adélie penguin colony (77°27′S, 169°12′E), Ross Sea, Antarctica. All cloacal samples were taken under the auspices of the long-term research program in place at Cape Crozier since the 1996–97 austral summer which individually marked a sample (500–1,000) of near-fledglings with numbered stainless-steel flipper bands each year.</p><p>Cloacal swabs of twenty-five birds ranging 7–18 years of age (eight females and seventeen males) were taken in December 2014 (during the breeding season) and stored in UTM™ Viral Transport Media (Copan). UTM™ Viral Transport Media (1000 µl) was filtered through 0.2 µm syringe filters, and 200 µl of this was subsequently used to extract viral DNA using the High Pure Viral Nucleic Acid Kit (Roche Diagnostics). Viral DNA (1 µl) from each sample was used to amplify circular DNA molecules using rolling-circle amplification (RCA) with the TempliPhi™ kit (GE Healthcare). The RCA DNA (5 µl aliquot from each sample) was pooled and sequenced on an Illumina 4000 (Illumina) sequencer at Macrogen Inc. (Korea).</p><p>The paired-end reads were <italic>de novo</italic> assembled using ABySS v1.9 (Simpson et al. 2009) and contigs >750 nucleotides were analyzed using BLASTx against a local viral database. A contig of 7,717 nucleotides was identified that encodes proteins that are homologous to those of papillomavirus. Further analysis of this contig revealed a 63 nucleotide terminal repeat and with the repeated region removed this contig represented a unit length papillomavirus. A pair of abutting primers (AdP2F 5'- CCA TGA GAA ACC TTT TAG TCT TCC CGG G-3' and AdP2R 5'-GGT ATG AAA GAG CCA GTG TCA GGC TTA TTG-3') were designed in the minor capsid protein (L2) open reading frame (ORF) to screen and recover the full genome of this papillomavirus. The RCA product (0.5 µl) from the viral DNA from the cloacal swab samples was used as a template for polymerase chain reaction (PCR) to screen individual Adélie penguin samples for the papillomavirus with the AdP2F/R primers with the following thermal cycling protocol with Kapa HiFi Hotstart DNA polymerase: 95 °C for 3 minutes; twenty-five cycles of 98 °C for 20 seconds, 60 °C for 15 seconds, 72 °C for 8 minutes, and a final extension of 72 °C for 8 minutes followed by a cooling step to 4 °C for 10 minutes. The PCR products were resolved on a 0.7 per cent agarose gel stained with SYBR safe, and the ∼7.5 kb amplicon was excised, purified, and cloned into the pJET1.2 plasmid vector (ThermoFisher), and Sanger sequenced by primer walking at Macrogen Inc. (Korea). The Sanger sequencing contigs were assembled using DNA Baser (Heracle BioSoft S.R.L. Romania).</p></sec><sec><title>2.3 Sequence analysis</title><p>Viral ORFs were predicted using a custom, papillomavirus-specific annotation tool as described previously (<xref rid="vex027-B56" ref-type="bibr">Van Doorslaer et al. 2013</xref>). Genome-wide identity analysis was performed using default settings within SDT v1.2 (<xref rid="vex027-B32" ref-type="bibr">Muhire et al. 2014</xref>). Pairwise sequence comparisons with PaPV1 (<xref rid="vex027-B59" ref-type="bibr">Varsani et al. 2014</xref>) were performed within Geneious version 10.2 (<xref rid="vex027-B25" ref-type="bibr">Kearse et al. 2012</xref>). The annotated sequence of PaPV1 was downloaded from the PaVE database (<xref rid="vex027-B56" ref-type="bibr">Van Doorslaer et al. 2013</xref>; <xref rid="vex027-B57" ref-type="bibr">Van Doorslaer et al. 2017</xref>). The L-INS-I algorithm within MAFFT v7.3 (<xref rid="vex027-B23" ref-type="bibr">Katoh et al. 2002</xref>; <xref rid="vex027-B24" ref-type="bibr">Katoh and Standley 2013</xref>) was used to align the translated protein sequence for each individual ORF (using the Blosum62 scoring matrix). Translated protein sequences were back translated to obtain codon-conserving nucleotide alignments. To aid in L1-based classification, the L1 ORF located between nucleotides 6,114 and 7,637 (PaPV1) and 6,116 and 7,654 (PaPV2) was used. The sequence of novel papillomavirus from Adélie penguin is deposited on GenBank under accession number MF168943.</p></sec><sec><title>2.4 Maximum likelihood phylogenetic tree construction</title><p>For genera containing more than five virus types, five viruses were randomly chosen to represent the genus, for other genera all virus sequences were included. Unclassified viruses, i.e. viruses not yet officially assigned to a genus by the ICTV, were excluded from the analysis. The sequences for the viral E1, E2 and L1 proteins of the 120 viruses were downloaded from PaVE (<xref rid="vex027-B56" ref-type="bibr">Van Doorslaer et al. 2013</xref>; <xref rid="vex027-B57" ref-type="bibr">Van Doorslaer et al. 2017</xref>). Individual proteins were aligned using MAFFT v7.3 (<xref rid="vex027-B23" ref-type="bibr">Katoh et al. 2002</xref>; <xref rid="vex027-B24" ref-type="bibr">Katoh and Standley 2013</xref>), E1 and E2 were aligned using the E-INS-I algorithm, G-INS-I was used for the L1 alignment. Aligned sequences were concatenated. A partitioning scheme allowing each protein to evolve according to the LG + I+G + F model of evolution was identified to be the best fit under the corrected Akaike information criterion (AICc) as implemented within PartitionFinder (<xref rid="vex027-B16" ref-type="bibr">Guindon et al. 2010</xref>; <xref rid="vex027-B26" ref-type="bibr">Lanfear et al. 2017</xref>). Maximum likelihood (ML) trees were estimated using RAxML-HPC v.8 on XSEDE (<xref rid="vex027-B30" ref-type="bibr">Miller et al. 2015</xref>), performing a rapid bootstrap analysis followed by a search for the best tree. The number of bootstrap replicates (<italic>n</italic> = 108; Weighted Robinson-Foulds = 2.46) were determined according to the autoMRE criteria implemented in RAxML (<xref rid="vex027-B35" ref-type="bibr">Pattengale et al. 2010</xref>; <xref rid="vex027-B50" ref-type="bibr">Stamatakis 2014</xref>). The resulting phylogenetic tree displayed represents the best ML tree.</p></sec><sec><title>2.5 Evolutionary hypothesis testing</title><p>Four distinct evolutionary hypotheses were compared (phylogenetic trees in Newick format are available on FigShare). The constrained trees were estimated using the same settings as used for the non-constrained ML tree (previous paragraph). The RAxML calculations started with an incomplete constraint, i.e. only the relevant nodes were constrained while the position of the remaining taxa was optimized under ML respecting the given constraints. Site-likelihoods associated with each tree were calculated using Tree-Puzzle (<xref rid="vex027-B46" ref-type="bibr">Schmidt et al. 2002</xref>) and compared using different hypothesis testing methods implemented in CONSEL (<xref rid="vex027-B49" ref-type="bibr">Shimodaira and Hasegawa 2001</xref>). The tests implemented are based on the RELL bootstrapping method. Hence, the parameters and branch lengths estimated by RAxML are used for the analysis. A significance limit of α<italic> </italic>=<italic> </italic>0.05 was used for all analyses.</p></sec><sec><title>2.6 Dating of the FgPV1 divergence</title><p>The E1, E2, and L1 protein sequences were aligned as described above. The divergence times were estimated using an uncorrelated relaxed clock (<xref rid="vex027-B11" ref-type="bibr">Drummond et al. 2006</xref>) with a calibrated Yule speciation model (<xref rid="vex027-B19" ref-type="bibr">Heled and Drummond 2012</xref>) in BEAST 2.4 (<xref rid="vex027-B28" ref-type="bibr">Le and Gascuel 2008</xref>; <xref rid="vex027-B19" ref-type="bibr">Heled and Drummond 2012</xref>; <xref rid="vex027-B7" ref-type="bibr">Bouckaert et al. 2014</xref>) as implemented on CIPRES (<xref rid="vex027-B30" ref-type="bibr">Miller et al. 2015</xref>). The partitioning scheme and evolutionary model were the same as in that used in evolutionary hypothesis testing described above. Site models and phylogenetic trees were estimated for the entire alignment. An uncorrelated relaxed clock model was assigned to each Individual protein. The previously estimated ML tree was used as a starting tree, with monophyletic requirements imposed on the calibration nodes. Lognormal priors were used to calibrate nodes 2, 9, 25, and 26 based on the divergence estimates of the hosts. Calibration means and 95 per cent confidence interval were downloaded from TimeTree (<ext-link ext-link-type="uri" xlink:href="http://www.timetree.org">www.timetree.org</ext-link>; <xref rid="vex027-B17" ref-type="bibr">Hedges et al. 2006</xref>; <xref rid="vex027-B18" ref-type="bibr">Hedges et al. 2015</xref>). The MCMC chains were run for 1 × 10<sup>8</sup> generations. Samples were collected every 5,000th generation. The analysis was repeated for a total of six independent runs. The resulting log files and trees were combined using LogCombiner and resampled (every 10,000th generation) while discarding a 20 per cent burn-in. An estimated sample size of >200 was achieved for most parameters. The reported estimates (<xref ref-type="fig" rid="vex027-F2">Fig. 2</xref>; <xref ref-type="supplementary-material" rid="sup2">Supplementary Table S1</xref>) are based on the combined dataset.</p><p>A posterior distribution of potential phylogenetic trees for the E1–E2–L1 concatenated alignment described above was estimated using MrBayes 3.2.6 (<xref rid="vex027-B21" ref-type="bibr">Huelsenbeck and Ronquist 2001</xref>; <xref rid="vex027-B44" ref-type="bibr">Ronquist and Huelsenbeck 2003</xref>; <xref rid="vex027-B43" ref-type="bibr">Ronquist et al. 2012</xref>; <xref rid="vex027-B30" ref-type="bibr">Miller et al. 2015</xref>). Two chains of MCMC were run for 1,000,000 steps, and parameters were logged every 1,000 steps with the initial 25 per cent discarded as burn in. Chain convergence was confirmed by comparing different parameter’s potential scale reduction factors (∼1.000), by analyzing individual traces in Tracer v1.6 (<ext-link ext-link-type="uri" xlink:href="http://beast.bio.ed.ac.uk/Tracer">http://beast.bio.ed.ac.uk/Tracer</ext-link>) and by using AWTY (<xref rid="vex027-B33" ref-type="bibr">Nylander et al. 2008</xref>). The analysis was repeated three separate times, and post burn-in samples were combined. The time-dependent rate phenomenon can be modeled using a power-law rate-decay model (T = αS<sup>β</sup>). For each of the estimated posterior Bayesian phylogenies (<italic>n</italic> = 5905), a set of patristic distances (S) was calculated starting from either PePV1 or CPV6. Matching divergence time estimates (T) associated with the host species were downloaded from TimeTree (<ext-link ext-link-type="uri" xlink:href="http://www.timetree.org">www.timetree .org</ext-link>) (<xref rid="vex027-B17" ref-type="bibr">Hedges et al. 2006</xref>; <xref rid="vex027-B18" ref-type="bibr">Hedges et al. 2015</xref>). The values of S and T were log-transformed (base 10), and the average ±SD was plotted. Generalized Deming regression (<xref rid="vex027-B41" ref-type="bibr">Ripley and Thompson 1987</xref>) was performed using the Deming function in R 3.2.1 (<xref rid="vex027-B52" ref-type="bibr">Team RC 2014</xref>). The linear model was then extrapolated to estimate the timescale of other nodes from their S estimates. Data, Python, and R scripts to perform these analyses are available through FigShare. </p></sec></sec><sec><title>3. Results</title><sec><title>3.1 Identification of a novel papillomavirus isolated from an Adélie penguin</title><p>One of the twenty-five cloacal swab samples of Adélie penguins from Cape Crozier on Ross Island (16-year-old female) was positive for a novel papillomavirus. This novel papillomavirus shares ∼64 per cent genome-wide sequence identity with Pygoscelis adeliae papillomavirus 1 (PaPV1), which was identified from a fecal sample retrieved at Cape Crozier (<xref rid="vex027-B59" ref-type="bibr">Varsani et al. 2014</xref>). Based on the recommendation by the animal papillomavirus reference center (<xref rid="vex027-B6" ref-type="bibr">Bernard et al. 2010</xref>), the novel papillomavirus is designated Pygoscelis adeliae papillomavirus type 2 (PaPV2).</p><p>The 7,654-bp double-stranded genome of PaPV2 contains seven ORFs ( E6, E7, E1, E9, E2, L2, and L1) on the same coding strand (<xref ref-type="fig" rid="vex027-F1">Fig. 1</xref>). Both pairwise identities and phylogenetic analysis show that PaPV2 is most closely related to PaPV1 (<xref ref-type="fig" rid="vex027-F1">Figs 1 and 2</xref>). The L1 gene of PaVP1 and PaPV2 share 62.5 per cent nucleotide sequence identity thus suggesting that both these Adélie penguin-associated papillomaviruses likely belong to the genus <italic>Treisepsilonpapillomavirus</italic> (<xref ref-type="table" rid="vex027-T1">Table 1</xref>).
<table-wrap id="vex027-T1" orientation="portrait" position="float"><label>Table 1.</label><caption><p> Pairwise sequence comparison between PaPV1 and PaPV2.</p></caption><table frame="hsides" rules="groups"><colgroup span="1"><col valign="top" align="left" span="1"></col><col valign="top" align="char" char="." span="1"></col><col valign="top" align="char" char="." span="1"></col></colgroup><thead><tr><th rowspan="1" colspan="1">PaPV1 vs. PaPV2</th><th rowspan="1" colspan="1">Nucleotide identity<xref ref-type="table-fn" rid="tblfn1"><sup>a</sup></xref></th><th rowspan="1" colspan="1">Amino acid similarity<xref ref-type="table-fn" rid="tblfn2"><sup>b</sup></xref></th></tr></thead><tbody><tr><td rowspan="1" colspan="1">E6</td><td rowspan="1" colspan="1">43.4</td><td rowspan="1" colspan="1">60.0</td></tr><tr><td rowspan="1" colspan="1">E7</td><td rowspan="1" colspan="1">45.3</td><td rowspan="1" colspan="1">62.1</td></tr><tr><td rowspan="1" colspan="1">E1</td><td rowspan="1" colspan="1">63.4</td><td rowspan="1" colspan="1">80.2</td></tr><tr><td rowspan="1" colspan="1">E9</td><td rowspan="1" colspan="1">53.7</td><td rowspan="1" colspan="1">79.5</td></tr><tr><td rowspan="1" colspan="1">E2</td><td rowspan="1" colspan="1">49.8</td><td rowspan="1" colspan="1">68.7</td></tr><tr><td rowspan="1" colspan="1">L1</td><td rowspan="1" colspan="1">62.5</td><td rowspan="1" colspan="1">86.2</td></tr><tr><td rowspan="1" colspan="1">L2</td><td rowspan="1" colspan="1">53.2</td><td rowspan="1" colspan="1">74.5</td></tr></tbody></table><table-wrap-foot><fn id="tblfn1"><p><sup>a</sup>Pairwise, codon conserving aligment.</p></fn><fn id="tblfn2"><p><sup>b</sup>Blosum62 matrix with threshold = 0.</p></fn></table-wrap-foot></table-wrap></p><fig id="vex027-F1" orientation="portrait" position="float"><label>Figure 1.</label><caption><p> Genomic organization of Sauropsid papillomaviruses. Illustration of the genomic organization of the known bird and turtle papillomaviruses, i.e. Caretta caretta papillomavirus 1 (CcPV1), Chelonia mydas papillomavirus 1 (CmPV1), Francolinus leucoscepus papillomavirus 1, (FlPV1), Fringilla coelebs papillomavirus 1 (PcPV1), Psittacus erithacus papillomavirus 1 (PePV1), Fulmar glacialis papillomavirus type 1 (FgPV1), and Pygoscelis adeliae papillomavirus 1 and 2 (PaPV1 and PaPV2). PaPV2, the virus isolated as part of this report is underlined. The open reading frames are color coded. Phylogenetic tree on the left shows the evolutionary relationships between the different viruses. Colored/numbered branches refer to “molecular morphology traits”, see main text for details. The mammalian virus clade contains MsPV1, the virus isolated from a lesion on a Carpet python. The ‘H–Θ–3E’ viruses are indicated.</p></caption><graphic xlink:href="vex027f1"></graphic></fig><fig id="vex027-F2" orientation="portrait" position="float"><label>Figure 2.</label><caption><p> Maximum likelihood phylogenetic tree showing the relationship between the sauropsid and mammalian papillomaviruses. Colored nodes indicate bootstrap support (blue=100%, green >85%). The branches leading to the sauropsid viruses are shown in bold and orange. Despite being isolated from a python, the Morelia spilota papillomavirus type 1 (MsPV1; indicated with an asterisk) clusters with the mammalian viruses. TempEst (<xref rid="vex027-B37" ref-type="bibr">Rambaut et al. 2016</xref>) was used to identify the root of the tree.</p></caption><graphic xlink:href="vex027f2"></graphic></fig></sec><sec><title>3.2 Phylogenetic analysis clusters mammalian and sauropsid papillomaviruses into separate monophyletic clades</title><p>A substantial part of papillomavirus evolutionary history can be explained by virus-host co-speciation (<xref rid="vex027-B54" ref-type="bibr">Van Doorslaer 2013</xref>). The recent identification of fish papillomaviruses (<xref rid="vex027-B29" ref-type="bibr">Lopez-Bueno et al. 2016</xref>) extends the evolutionary history of this group of viruses to at least 400 million years ago (MYA) (Osteichthyes and Amniota shared a last common ancestor around 400 MYA (<xref rid="vex027-B10" ref-type="bibr">dos Reis et al. 2015</xref>). To reconstruct the papillomavirus evolutionary history, we estimated a phylogenetic tree including all known non-mammalian papillomaviruses (including the fish papillomavirus Sparus aurata papillomavirus 1 (SaPV1), bird (<italic>n</italic> = 6), snake (MsPV1; <italic>n</italic> = 1), and testudines (<italic>n</italic> = 2) infecting papillomaviruses), as well as representative papillomaviruses associated with different mammalian hosts (<italic>n</italic> = 109). A ML tree was constructed based on a concatenated protein alignment of E1, E2, and L1 (<xref ref-type="fig" rid="vex027-F2">Fig. 2</xref>). Assuming the existence of a molecular clock, the software program TempEst (<xref rid="vex027-B37" ref-type="bibr">Rambaut et al. 2016</xref>) identified the fish papillomavirus (SaPV1) as the root of the papillomavirus phylogenetic tree. Consistent with the hypothesis of co-speciation, the initial split between the mammalian and non-mammalian viruses matches the equivalent split between their hosts; the papillomavirus tree shows two well-separated, monophyletic clades. The first clade clusters the mammalian viruses, with the turtle and bird viruses clustering in a separate, monophyletic clade. As a notable exception, MsPV1 (the papillomavirus isolated from a <italic>Morelia spilota</italic> (Carpet python) lesion), clusters with the mammalian viruses.</p></sec><sec><title>3.3 No strict co-speciation within the turtle-bird clade</title><p>The phylogenetic tree with the highest likelihood (<xref ref-type="fig" rid="vex027-F2">Fig. 2</xref>) suggests that the avian papillomaviruses (<italic>n</italic> = 5) form a paraphyletic group. This reconstruction places the Fulmar glacialis papillomavirus 1 (FgPV1), a virus infecting the Northern Fulmar, sister to a clade containing the turtle and all other avian papillomaviruses. To formally test this hypothesis, we compared different tree topologies using the approximately unbiased (AU) test (<xref rid="vex027-B48" ref-type="bibr">Shimodaira 2002</xref>) as implemented within CONSEL (<xref rid="vex027-B49" ref-type="bibr">Shimodaira and Hasegawa 2001</xref>). We tested the following alternative hypotheses: 1, FgPV1 speciation predates divergence of bird and turtle viruses (i.e. the ML solution); 2, all avian papillomaviruses cluster together, separate from the turtle papillomaviruses; 3, two monophyletic clades containing Caretta caretta papillomavirus 1 (CcPV1), Chelonia mydas papillomavirus 1 (CmPV1), and avian papillomaviruses (Francolinus leucoscepus papillomavirus 1, FlPV1; Fringilla coelebs papillomavirus 1, PcPV1; Psittacus erithacus papillomavirus 1, PePV1; FgPV1, PaPV1, and PaPV2), respectively; 4, constraint based on the accepted host phylogeny (<xref rid="vex027-B36" ref-type="bibr">Prum et al. 2015</xref>). Using an α of 0.05, only the strict host phylogeny constraint (option 4) was rejected as an alternative topology (<xref ref-type="table" rid="vex027-T2">Table 2</xref>). This suggests that the exact phylogenetic position of FgPV1 cannot be determined. In addition, as is the case in the mammalian clade, strict co-speciation within the non-mammalian papillomavirus clade should be rejected in favor of more complex evolutionary history (<xref rid="vex027-B54" ref-type="bibr">Van Doorslaer 2013</xref>).
<table-wrap id="vex027-T2" orientation="portrait" position="float"><label>Table 2.</label><caption><p> Statistical comparison of four distinct phylogenetic hypotheses.</p></caption><table frame="hsides" rules="groups"><colgroup span="1"><col valign="top" align="left" span="1"></col><col valign="top" align="char" char="." span="1"></col><col valign="top" align="char" char="." span="1"></col><col valign="top" align="char" char="." span="1"></col><col valign="top" align="char" char="." span="1"></col><col valign="top" align="char" char="." span="1"></col><col valign="top" align="char" char="." span="1"></col><col valign="top" align="char" char="." span="1"></col></colgroup><thead><tr><th rowspan="1" colspan="1">Phylogenetic tree (newick)</th><th rowspan="1" colspan="1">Rank</th><th rowspan="1" colspan="1">δ log likelihood to best tree</th><th rowspan="1" colspan="1">AU</th><th rowspan="1" colspan="1">KH</th><th rowspan="1" colspan="1">SH</th><th rowspan="1" colspan="1">wKH</th><th rowspan="1" colspan="1">wSH</th></tr></thead><tbody><tr><td rowspan="1" colspan="1">Unconstrained</td><td rowspan="1" colspan="1">1</td><td rowspan="1" colspan="1">0</td><td rowspan="1" colspan="1">0.756</td><td rowspan="1" colspan="1">0.636</td><td rowspan="1" colspan="1">0.919</td><td rowspan="1" colspan="1">0.636</td><td rowspan="1" colspan="1">0.928</td></tr><tr><td rowspan="1" colspan="1">(SaPV1,(((CcPV1,CmPV1),(<bold>FgPV1</bold>,FlPV1,PePV1,FcPV1,PaPV2,PaPV1)), (mammalian*)));</td><td rowspan="1" colspan="1">2</td><td rowspan="1" colspan="1">6.6</td><td rowspan="1" colspan="1">0.476</td><td rowspan="1" colspan="1">0.364</td><td rowspan="1" colspan="1">0.805</td><td rowspan="1" colspan="1">0.364</td><td rowspan="1" colspan="1">0.795</td></tr><tr><td rowspan="1" colspan="1">(SaPV1,((<bold>FgPV1</bold>,(CcPV1,CmPV1,FlPV1,PePV1,FcPV1,PaPV2,PaPV1)), (mammalian*)));</td><td rowspan="1" colspan="1">3</td><td rowspan="1" colspan="1">17.6</td><td rowspan="1" colspan="1">0.242</td><td rowspan="1" colspan="1">0.205</td><td rowspan="1" colspan="1">0.476</td><td rowspan="1" colspan="1">0.205</td><td rowspan="1" colspan="1">0.394</td></tr><tr><td rowspan="1" colspan="1">(SaPV1,(((<bold>FgPV1</bold>,CcPV1,CmPV1),(FlPV1,PePV1,FcPV1,PaPV2,PaPV1)), (mammalian*)));</td><td rowspan="1" colspan="1">4</td><td rowspan="1" colspan="1">18.4</td><td rowspan="1" colspan="1">0.242</td><td rowspan="1" colspan="1">0.221</td><td rowspan="1" colspan="1">0.488</td><td rowspan="1" colspan="1">0.221</td><td rowspan="1" colspan="1">0.445</td></tr><tr><td rowspan="1" colspan="1">(SaPV1,(((CcPV1,CmPV1),(FlPV1,((FcPV1,PePV1),(<bold>FgPV1</bold>,(PaPV2,PaPV1))))), (mammalian*)));</td><td rowspan="1" colspan="1">5</td><td rowspan="1" colspan="1">119.2</td><td rowspan="1" colspan="1">0.002</td><td rowspan="1" colspan="1">0</td><td rowspan="1" colspan="1">3.00E-04</td><td rowspan="1" colspan="1">0</td><td rowspan="1" colspan="1">0</td></tr></tbody></table><table-wrap-foot><fn id="tblfn3"><p>AU, <italic>P</italic>-value for the approximately unbiased test; KH, <italic>P</italic>-value for the Kishino–Hasegawa test; SH, <italic>P</italic>-value for the Shimodaira–Hasegawa test; wKH, <italic>P</italic>-value for the weighted Kishino–Hasegawa test; wSH, <italic>P</italic>-value for the weighted Shimodaira–Hasegawa test.</p></fn></table-wrap-foot></table-wrap></p></sec><sec><title>3.4 The E9 protein is a synapomorphy of the H−Θ–3E Clade</title><p>Phylogenetic analysis based on three viral proteins did not enable us to determine the exact evolutionary history of FgPV1. However, all tested topologies clustered the other bird viruses (FcPV1, FlPV1, PePV1, PaPV1, and PaPV2 belonging to the <italic>Eta-</italic>, <italic>Theta-</italic>, and <italic>Treisepsilonpapillomavirus</italic> genera) and turtle viruses (CcPV1 and CmPV1; genus <italic>Dyozetapapillomavirus</italic>) into two distinct, monophyletic clades. From here on we refer to the avian viruses in the <italic>Eta-</italic>, <italic>Theta-</italic>, and <italic>Treisepsilonpapillomavirus</italic> monophyletic clade as ‘H−Θ−3E’. To date, H − Θ−3E clusters all known avian papillomaviruses except the fulmar FgPV1 virus.</p><p>Of note, all known members of the H − Θ−3E clade contain a unique, highly conserved ORF located within the E1 + 1 ORF (light blue bars in <xref ref-type="fig" rid="vex027-F3">Fig. 3A</xref>), known as E9. The E9 protein is 144–200 amino acids in length. Importantly, E9 is located within the N-terminal regulatory domain of E1. Unlike the other regions of the E1 protein, this N-terminal regulatory domain is relatively variable (<xref ref-type="fig" rid="vex027-F3">Fig. 3A;</xref><xref rid="vex027-B4" ref-type="bibr">Bergvall et al. 2013</xref>), suggesting that E9 is not simply the result of sequence constraints dictated by the E1 ORF. Importantly, FgPV1 and the other non-H − Θ−3E viruses do not contain an E9-like ORF.
</p><fig id="vex027-F3" orientation="portrait" position="float"><label>Figure 3.</label><caption><p> Genetic plasticity in the sauropsid papillomaviruses. (A) The E1 ORF of the sauropsid papillomaviruses and Bovine papillomavirus type 1 (BPV1) were aligned. In this alignment, thin lines indicate indels, while solid black bars indicate confidently aligned regions. The phylogenetic tree indicates the evolutionary history. The bars above the alignment indicate the different E1 protein domains identified in BPV1 (<xref rid="vex027-B4" ref-type="bibr">Bergvall et al. 2013</xref>). The +1 frame of this alignment was investigated for the presence of a large, conserved ORF known as E9 (light blue bars, <xref rid="vex027-B55" ref-type="bibr">Van Doorslaer et al. 2009</xref>). The avian viruses in the <italic>Eta-</italic>, <italic>Theta-</italic>, and <italic>Treisepsilonpapillomavirus</italic> genera (‘H–Θ–3E ’) contain this ORF. Like all mammalian viruses, the BPV1 E8 exon is located in a similar position (green bar). (B) The E7 protein of the sauropsid papillomaviruses and human papillomavirus type 16 (HPV16) were aligned. In this alignment, thin lines indicate indels, solid gray bars indicate confidently aligned regions, and black bars indicate identical residues. The phylogenetic tree indicates the evolutionary history. The bars above the alignment indicate the different E7 protein domains identified in HPV16 (<xref rid="vex027-B42" ref-type="bibr">Roman and Munger 2013</xref>). The rectangle indicates the part of the alignment that is enlarged in C. (C) Amino-acid level view of part of the alignment shown in B. The alignment is shaded by conservation. The conserved LxCxE and Zinc-binding motifs are indicated (<xref rid="vex027-B42" ref-type="bibr">Roman and Munger 2013</xref>).</p></caption><graphic xlink:href="vex027f3"></graphic></fig></sec><sec><title>3.5 The use of ‘molecular morphology’ provides support for a single topology</title><p>In addition to the presence of a unique E9 within the members of the H–Θ-3E clade, all sauropsid papillomaviruses have several unique features that may help with anchoring the evolutionary history of these viruses. Unlike the prototypical ‘mammalian virus’ E6 protein that consists of two 70-residue zinc-binding repeats (<xref rid="vex027-B58" ref-type="bibr">Vande Pol and Klingelhutz 2013</xref>), the sauropsid papillomaviruses contain an E6 protein that consists of a single zinc-binding domain (<xref rid="vex027-B55" ref-type="bibr">Van Doorslaer et al. 2009</xref>). Since all sauropsid papillomaviruses have this feature, the single-domain E6 was likely present in the most recent common ancestor of the (known) sauropsid papillomaviruses, implying it was present before the acquisition of the E9 protein (blue branch indicated by ‘1’ in <xref ref-type="fig" rid="vex027-F1">Fig. 1</xref>). Interestingly, PePV1 (one of the H − Θ−3E viruses) appears to have lost the E6 ORF (purple branch labeled ‘2’ in <xref ref-type="fig" rid="vex027-F1">Fig. 1</xref>).</p><p>Furthermore, the typical ‘mammalian virus’ E7 consists of a natively unfolded N-terminal region followed by a zinc-binding domain, which folds as an obligate homodimer (<xref rid="vex027-B42" ref-type="bibr">Roman and Munger 2013</xref>). The avian E7 proteins contain an extended unfolded N terminus, while the zinc-binding domain is relatively smaller when compared to the mammalian E7 (<xref rid="vex027-B55" ref-type="bibr">Van Doorslaer et al. 2009</xref>). Remarkably, the turtle E7-like proteins have completely lost the N-terminal region and acquired a shorter C-terminal extension (<xref ref-type="fig" rid="vex027-F3">Fig. 3B;</xref> orange branch labeled ‘3’ in <xref ref-type="fig" rid="vex027-F1">Fig. 1</xref>). Despite the loss of the CR1 and CR2 motifs, the turtle E7 zinc binding motif is more closely related to the H − Θ−3E papillomaviruses than to the FgPV1 E7 (<xref ref-type="fig" rid="vex027-F3">Fig. 3C</xref>).</p><p>The most parsimonious explanation suggests that the last common ancestor of FgPV1 and the sauropsid papillomaviruses contained a single domain E6, and an E7 protein consisting of CR1, CR2, and a relatively short zinc binding motif. Following the divergence of the turtle and the H − Θ−3E papillomaviruses, the avian papillomaviruses gained an E9 ORF, while the turtle papillomaviruses lost the CR1-CR2 motifs typical of E7. In addition, the duplication event giving rise to the prototypical E6 and the expansion of the E7 zinc binding motif must have occurred specifically on the clade giving rise to the mammalian papillomaviruses. This scenario can be traced along the branches of the ML solution tree (<xref ref-type="fig" rid="vex027-F2">Fig. 2</xref>). Despite being statistically equivalent, the other topologies tested above would need to invoke additional events to explain the evolution of the E6, E7, and E9 proteins.</p><p>Combining the inconclusive phylogenetic clustering with this additional molecular morphology data enables us to hypothesize that FgPV1 diverged from the other sauropsid papillomaviruses before the avian-turtle host animal diverged. This conclusion is further supported by the observation that FgPV1 lesions are uniquely associated with cartilage, while all other known papillomaviruses are primarily epithelial trophic (<xref rid="vex027-B14" ref-type="bibr">Gaynor et al. 2015</xref>).</p></sec><sec><title>3.6 Divergence of FgPV1 from the other sauropsid viruses may coincide with the paleozoic–mesozoic transition</title><p>We estimated the divergence time between FgPV1 and the other Sauropsid papillomaviruses. The time-calibrated phylogenetic tree (<xref ref-type="fig" rid="vex027-F4">Fig. 4</xref>) was calibrated using host divergence estimates. The most recent common ancestor (MRCA) of extant amniotes (node 2; 312 MYA (297–326 MYA)) and turtles (node 9; 47.5 MYA (30.5-62.8 MYA)), as well as the divergence estimates for the MRCA of raccoons and dogs (node 25; 54 MYA (52–57 MYA)), and hyenas and dogs (node 26; 40 MYA (33–46 MYA)) were used. Calibration means and 95 per cent confidence interval were downloaded from TimeTree (<ext-link ext-link-type="uri" xlink:href="http://www.timetree.org">www.timetree.org</ext-link>; <xref rid="vex027-B17" ref-type="bibr">Hedges et al. 2006</xref>; <xref rid="vex027-B18" ref-type="bibr">Hedges et al. 2015</xref>). This analysis dates the divergence of FgPV1 to about 262 MYA (95% highest probability density 228–294 MYA). Interestingly, this analysis dates the root of phylogenetic tree around 481 MYA (326–656 MYA). However, while the older divergence dates matched the hosts fairly well, this analysis showed that many of the more recent viral divergence estimates are significantly older than what has been shown for the associated hosts. It has been shown that evolutionary rate estimates depend on the age of the nodes used as calibration points (Duchene et al. <xref rid="vex027-B12" ref-type="bibr">2014a</xref>,<xref rid="vex027-B13" ref-type="bibr">b</xref>; <xref rid="vex027-B1" ref-type="bibr">Aiewsakun and Katzourakis 2015</xref>). It was recently suggested that this time-dependent rate phenomenon (TDRP) could be modeled using a power-law decay function (<xref rid="vex027-B1" ref-type="bibr">Aiewsakun and Katzourakis 2015</xref>, <xref rid="vex027-B2" ref-type="bibr">2017</xref>). Therefore, the relationship between the Patristic Distance (S) and evolutionary timescales (T) should also follow the same power-law function. To correct for TDRP, we constructed a model for mammalian papillomaviruses by tracing the canine CPV6 backward in time (<xref ref-type="fig" rid="vex027-F5">Fig. 5A</xref>). CPV6 and its close relatives belong to the <italic>Gammapapillomaviridae</italic>, a genus that has previously been shown to share a co-speciation history with their mammalian hosts (<xref rid="vex027-B38" ref-type="bibr">Rector et al. 2007</xref>). Similarly, a TDRP model was generated by calculating backward from the parrot virus PePV1 (<xref ref-type="fig" rid="vex027-F5">Fig. 5A</xref>). Each approach provided us with five pairs of S and T estimates we used for model construction (<xref ref-type="fig" rid="vex027-F5">Fig. 5B</xref>). Importantly, both TDRP models had overlapping error, suggesting that they give similar solutions (<xref ref-type="fig" rid="vex027-F5">Fig. 5C</xref>). We then extrapolated both models separately to estimate the timescale of FgPV1 divergence based on its patristic distance estimate (S). By assuming that viruses infecting modern-day mammals and birds share similar TDRP dynamics as their ancestors while incorporating phylogenetic uncertainty, the FgPV1 divergence event was inferred to have occurred around between 230 and 250 MYA (PePV1: 235 MYA (208–264 MYA), CPV6: 252 MYA (227–280 MYA)). Both the BEAST and TDRP derived approaches provide similar estimates (<xref ref-type="table" rid="vex027-T3">Table 3</xref>). Collectively, these approaches estimate the niche adaptation event associated with the divergence of FgPV1 to have occurred around 250 MYA (<xref ref-type="table" rid="vex027-T3">Table 3</xref>), during the Paleozoic-Mesozoic transition. Furthermore, these approaches date the root of the papillomavirus tree between 400 and 600 MYA (<xref ref-type="table" rid="vex027-T3">Table 3</xref>).
<table-wrap id="vex027-T3" orientation="portrait" position="float"><label>Table 3.</label><caption><p> Point estimates and associated 95 per cent highest probability density interval.</p></caption><table frame="hsides" rules="groups"><colgroup span="1"><col valign="top" align="left" span="1"></col><col valign="top" align="char" char="." span="1"></col><col valign="top" align="center" span="1"></col><col valign="top" align="char" char="." span="1"></col><col valign="top" align="center" span="1"></col><col valign="top" align="char" char="." span="1"></col><col valign="top" align="center" span="1"></col></colgroup><thead><tr><th rowspan="1" colspan="1"></th><th colspan="2" rowspan="1">BEAST analysis<hr></hr></th><th colspan="4" rowspan="1">TDRP analysis<hr></hr></th></tr><tr><th rowspan="1" colspan="1"></th><th colspan="2" rowspan="1"></th><th colspan="2" rowspan="1">PePV1<hr></hr></th><th colspan="2" rowspan="1">CPV6<hr></hr></th></tr><tr><th rowspan="1" colspan="1"></th><th rowspan="1" colspan="1">Mean</th><th rowspan="1" colspan="1">95% HPD</th><th rowspan="1" colspan="1">Mean</th><th rowspan="1" colspan="1">95% HPD</th><th rowspan="1" colspan="1">Mean</th><th rowspan="1" colspan="1">95% HPD</th></tr></thead><tbody><tr><td rowspan="1" colspan="1">FgPV1</td><td rowspan="1" colspan="1">262</td><td rowspan="1" colspan="1">228–294</td><td rowspan="1" colspan="1">235</td><td rowspan="1" colspan="1">208–264</td><td rowspan="1" colspan="1">252</td><td rowspan="1" colspan="1">227–280</td></tr><tr><td rowspan="1" colspan="1">root</td><td rowspan="1" colspan="1">481</td><td rowspan="1" colspan="1">326–656</td><td rowspan="1" colspan="1">604</td><td rowspan="1" colspan="1">533–685</td><td rowspan="1" colspan="1">489</td><td rowspan="1" colspan="1">437–547</td></tr></tbody></table><table-wrap-foot><fn id="tblfn4"><p>95% HPD: 95 per cent highest probability density interval.</p></fn></table-wrap-foot></table-wrap></p><fig id="vex027-F4" orientation="portrait" position="float"><label>Figure 4.</label><caption><p> Time calibrated phylogenetic tree of 120 papillomaviruses. Host divergence time estimates were used to time calibrate the papillomavirus tree using BEAST2 (<xref rid="vex027-B7" ref-type="bibr">Bouckaert et al. 2014</xref>). The nodes indicated with the black circle correspond to the calibration nodes. Nodes were calibrated using a lognormal prior (see <xref ref-type="supplementary-material" rid="sup2">Supplementary Table S1</xref> for exact values). The blue bar around node 3 indicates the 95 per cent highest probability density associated with the time estimate for the FgPV1 divergence. Support values and exact estimates (with 95% highest posterior density) for each numbered node are presented in <xref ref-type="supplementary-material" rid="sup2">Supplementary Table S1</xref>. Sauropsid papillomavirus branches are highlighted in orange and bold. The tree was constructed in BEAST2. Tree and geological column were visualized using the APE (<xref rid="vex027-B34" ref-type="bibr">Paradis et al. 2004</xref>), phyloch (available from <ext-link ext-link-type="uri" xlink:href="http://www.christophheibl.de/Rpackages.html">http://www.christophheibl.de/Rpackages.html</ext-link>), and strap (<xref rid="vex027-B3" ref-type="bibr">Bell and Lloyd 2014</xref>) packages within R. The scale bar indicates millions of years before the present.</p></caption><graphic xlink:href="vex027f4"></graphic></fig><fig id="vex027-F5" orientation="portrait" position="float"><label>Figure 5.</label><caption><p> Evolutionary timescale of FgPV1 estimated by using the power-law rate-decay model. (A) The power-law rate-decay model (T = αS<sup>β</sup>) accounting for the time dependent rate phenomenon requires the patristic distance between different viruses (S) and the associated host divergence times (T) (<xref rid="vex027-B2" ref-type="bibr">Aiewsakun and Katzourakis 2017</xref>). To account for phylogenetic uncertainty, the post burn in posterior tree sample of three independent MrBayes (<xref rid="vex027-B21" ref-type="bibr">Huelsenbeck and Ronquist 2001</xref>; <xref rid="vex027-B44" ref-type="bibr">Ronquist and Huelsenbeck 2003</xref>; <xref rid="vex027-B43" ref-type="bibr">Ronquist et al. 2012</xref>) runs was used as the basis to calculate S (n<sub>trees</sub>=5905). S values were derived by back calculating from CPV6 and PePV1. T estimates for host species pairs were downloaded from <ext-link ext-link-type="uri" xlink:href="http://www.timetree.org">www.timetree.org</ext-link> (<xref rid="vex027-B17" ref-type="bibr">Hedges et al. 2006</xref>; <xref rid="vex027-B18" ref-type="bibr">Hedges et al. 2015</xref>). Individual S and T values were log<sub>10</sub> transformed, the mean and SD is shown. (B) S and T mean values with associated standard deviation were plotted. The numbers by the data points refer back to panel (A). The gray bars show the associated error for the FgPV1 and SaPV1 S estimates. Deming regression was used to obtain mean (±SD) values for the alpha and beta parameters. To obtain real time estimates, the reverse log of the solved values is reported (mean with 95% confidence interval). The regression line and associated confidence interval were drawn using regular linear regression, and is for illustrative purposes only. (C) Overlay of the curves in B showing overlap of derived estimates. (D) Posterior distribution of log likelihood values showing that the three separate MrBayes analyses converged on the same solution.</p></caption><graphic xlink:href="vex027f5"></graphic></fig></sec></sec><sec><title>4. Discussion</title><p>We identified a novel papillomavirus (PaPV2) isolated from a 16-year-old female Adélie penguin. This novel virus shares 62.5 per cent nucleotide identity across the L1 ORF when compared to its closest neighbor, PaPV1 (<xref ref-type="table" rid="vex027-T1">Table 1</xref>). This suggests that PaPV1 and PaPV2 (both isolated from Adélie penguin samples at Cape Crozier, Ross Island, Antarctica) likely belong to the genus <italic>Treisepsilonpapillomavirus</italic>. There are roughly twice as many species of birds compared to mammals (9956 versus 5416), yet only six papillomaviruses have been described associated with avian hosts (Adélie penguin, African Grey parrot, common chaffinch, northern fulmar, and yellow-necked francolin). Despite only having isolated a few non-mammalian viruses, these viruses have some unique features not shared with the mammalian viruses, suggesting that it is worthwhile to expand our understanding of non-mammalian papillomavirus genomes.</p><p>Interestingly, FgPV1, a virus isolated from a Northern fulmar does not contain many features shared by the other sauropsid viruses. Indeed, based on phylogenetic hypothesis testing, we were not able to pinpoint the exact evolutionary history of FgPV1. As in most biological systems, strict co-evolution of host and species does not fully explain the papillomavirus evolutionary history. However, a model based on niche adaptation near the root of the phylogenetic tree, followed by host-virus co-speciation has been proposed (<xref rid="vex027-B47" ref-type="bibr">Shah et al. 2010</xref>; <xref rid="vex027-B15" ref-type="bibr">Gottschling et al. 2011</xref>; <xref rid="vex027-B54" ref-type="bibr">Van Doorslaer 2013</xref>). Indeed, a sauropsid papillomavirus phylogeny corresponding to strict co-evolution was rejected (<xref ref-type="table" rid="vex027-T2">Table 2</xref>), suggesting a possible role for niche adaptation. However, since papillomaviruses typically have very similar genomic organization and infect epithelial cells, it has been hard to provide evidence for these niche switches as they are likely very subtle changes. We show that the use of ‘molecular morphological’ information allows us to more precisely determine the phylogenetic position of FgPV1 and more importantly provide evidence for such a niche adaptation event. Also, the unique plasticity of the early genes provides us with insights into the evolutionary history of the viral oncogenes.</p><p>The prototypical (mammalian virus) E6 protein consist of two highly conserved Zinc finger motifs (CxxC-x<sub>29</sub>-CxxC; where C is a cysteine and ‘x’ can be any amino acid) (<xref rid="vex027-B61" ref-type="bibr">Zanier et al. 2013</xref>). However, the E6 proteins of all known sauropsid papillomaviruses contain only one such motif. Since the fish papillomavirus (SaPV1) does not contain any such motif (<xref rid="vex027-B14" ref-type="bibr">Gaynor et al. 2015</xref>), it is likely that the most recent common ancestor of the mammalian and sauropsid papillomaviruses had a single motif E6, which duplicated to give rise to the mammalian papillomavirus double Zinc finger motif in E6 (<xref rid="vex027-B55" ref-type="bibr">Van Doorslaer et al. 2009</xref>). The sauropsid (putative) E7 protein may be even more remarkable. As shown in <xref ref-type="fig" rid="vex027-F3">Fig. 3</xref>, the sauropsid papillomaviruses express three distinct variants of this protein. While FgPV1 contains an E7 that is very similar to the mammalian E7 protein, the other bird viruses’ (‘H–Θ-3E’) E7 is more compact (<xref rid="vex027-B55" ref-type="bibr">Van Doorslaer et al. 2009</xref>). Remarkably, the two papillomaviruses associated with lesions on turtles have a unique E7 protein. Unlike the other sauropsid viruses, the turtle E7 protein completely lacks the non-folded CR1 and CR2 motifs typically observed in E7. For the mammalian viruses, these motifs have been shown to interact with many cellular partners (<xref rid="vex027-B42" ref-type="bibr">Roman and Munger 2013</xref>). For example, it is likely that these E7 proteins will not be able to bind to pocket family proteins (e.g. pRb), a function believed to be important for the (α-) papillomavirus lifecycle (<xref rid="vex027-B42" ref-type="bibr">Roman and Munger 2013</xref>). The E6 and E7 proteins of these sauropsid viruses seem to be highly plastic. Remarkably, compared to the core proteins (E1, E2, L1, and L2), the E6 and E7 proteins are significantly less variable, suggesting that these ‘adaptive’ proteins may be under less functional constraints.</p><p>In addition to encoding these sauropsid-specific E6 and E7 proteins, the ‘H–Θ-3E’ viruses may also contain a unique E9 ORF. This putative protein is encoded within the +1 ORF of E1 and is embedded within the N-terminal regulatory domain, a part of E1 that is relatively less conserved (<xref rid="vex027-B4" ref-type="bibr">Bergvall et al. 2013</xref>). This suggests that the evolutionary pressure to maintain the E9 protein is not (solely) due to constraints on the E1 protein. The mammalian viruses express a spliced transcript that encodes fusion protein that consists of a short peptide from an ORF designated E8 fused to the C-terminal domain of E2 (<xref rid="vex027-B22" ref-type="bibr">Jeckel et al. 2003</xref>). The resulting proteins (designated E8^E2) function as repressors of viral transcription and replication (<xref rid="vex027-B22" ref-type="bibr">Jeckel et al. 2003</xref>). Importantly, while the E8 exon of E8^E2 is located in the same general region of the E1 protein as E9, there are no other obvious similarities between both proteins, suggesting that both proteins have separate origins. It had been previously suggested that the N-terminal region of E1 may be a privileged segment of the protein more malleable to evolutionary innovation (<xref rid="vex027-B31" ref-type="bibr">Morin et al. 2011</xref>). While it is intriguing to note that a subset of bird papillomaviruses carry such an E9 ORF, in the absence of functional data it is hard to speculate about the functionality of this protein.</p><p>Importantly, MsPV1 (associated with a lesion on a python) does not cluster with any of the other known non-mammalian tropic papillomaviruses (<xref ref-type="fig" rid="vex027-F2">Fig. 2</xref>). The phylogenetic position of this virus significantly complicates the co-speciation hypothesis (with or without niche adaptation) between papillomaviruses and their hosts. This dilemma could be explained by invoking a cross-species infection (<xref rid="vex027-B27" ref-type="bibr">Lange et al. 2011</xref>). However, it is believed that papillomaviruses are highly species specific. Until this branch of the phylogenetic tree becomes populated with more reptilian viruses, it is impossible to rule out alternative explanations (e.g. contamination of the python lesion with mammalian virus DNA).</p><p>Inclusion of ‘molecular morphological’ supports the ML solution reached by the evolutionary sequence analysis, i.e. FgPV1 speciated before the split between the turtle and bird viruses. Of note, FgPV1 viral particles were found in the cartilaginous tissue of a Northern fulmar’s foot (<xref rid="vex027-B14" ref-type="bibr">Gaynor et al. 2015</xref>). Since papillomaviruses are typically associated with epithelial tissues, FgPV1 may be the prototype of a new group of papillomaviruses able to replicate within cartilage, a dramatically different host tissue. To estimate the timing of this niche adaptation event, we constructed a traditional host constrained time tree for viral evolution. In addition, due to concerns about TDRP, we estimated the date using a complimentary approach that has been suggested to be less sensitive to TDRP (<xref rid="vex027-B2" ref-type="bibr">Aiewsakun and Katzourakis 2017</xref>). Using both approaches, we estimated divergence of FgPV1 to have occurred around 250 MYA ago. This overlaps with the Paleozoic-Mesozoic transition. Not only was this transition associated with mass extinction but also coincides with the rise of endothermia (<xref rid="vex027-B40" ref-type="bibr">Rey et al. 2017</xref>). The combination of mass extinctions and (molecular) changes associated with the ability regulate internal temperature may have provided the right evolutionary environment for this niche adaptation event to be successful. Based on our analysis, the ancestor of the amniote-associated viruses infected epithelial cells. Along the evolutionary path toward the sauropsid viruses, intra-host divergence (i.e. independent of host speciation) (<xref rid="vex027-B9" ref-type="bibr">Buck et al. 2016</xref>) created an ancestral virus capable of infecting a novel host niche (i.e. chondrocytes). This hypothesis suggests that there should be several (if not hundreds) papillomaviruses with a tropism for cartilage. In addition, this analysis estimates that papillomaviruses have a marine origin and may have evolved along with their jawed vertebrate hosts. This is reminiscent of the ancient marine origin of retroviruses (<xref rid="vex027-B45" ref-type="bibr">Ruboyianes and Worobey 2016</xref>; <xref rid="vex027-B2" ref-type="bibr">Aiewsakun and Katzourakis 2017</xref>;).</p><p>As part of this study we isolated a novel papillomavirus from Adélie penguins. Detailed sequence analysis of this novel virus and other non-mammalian papillomaviruses provided us with novel insights into the evolution of these viruses. We provide evidence for genomic plasticity (specifically within the viral E6 and E7 genes). Furthermore, we generated the testable hypothesis that viral niche adaptation gave rise to a clade of viruses with a tropism for cartilage. A concerted effort to identify additional, novel viruses will allow us to fine tune our hypothesis on genome plasticity and whether niche adaptation events similar to the one seen in the Norther Fulmar are more common.</p></sec><sec sec-type="supplementary-material"><title>Supplementary Material</title><supplementary-material content-type="local-data" id="sup1"><label>Supplementary Text</label><media xlink:href="vex027_Supp.txt"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="sup2"><label>Supplementary Table 1</label><media xlink:href="vex027_Supp_Table1.xlsx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material></sec></body><back><ack><title>Acknowledgments</title><p>The fieldwork for this study was supported by the US National Science Foundation (NSF) under grant ANT-0944411 awarded to D.G.A., K.M.D. and G.B., with logistics supplied by the US Antarctic Program. K.V.D. is supported by a State of Arizona Improving Health TRIF and by the National Institute of Food and Agriculture, U.S. Department of Agriculture, Hatch NC229. The molecular work described in this study is supported by Center of Evolution and Medicine Venture Fund (Center of Evolution and Medicine, Arizona State University, USA) grant awarded to A.V. Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the US Government.</p><sec><title>Data availability</title><p>The annotated genome sequence of PapV2 is available on GenBank (Accession number: MF168943). Associated files and scripts to repeat the analyses are available on FigShare under the following DOIs: 10.6084/m9.figshare.5035202; 10.6084/m9.figshare.5035205; 10.6084/m9.figshare.5035208; 10.6084/m9.figshare.5035211; 10.6084/m9.figshare.5035214; 10.6084/m9.figshare.5035220; 10.6084/m9.figshare.5036153.</p></sec></ack><sec><title>Supplementary data</title><p><xref ref-type="supplementary-material" rid="sup1">Supplementary data</xref> are available at <italic>Virus Evolution</italic> online.</p><p>Conflict of interest: None declared.</p></sec><ref-list><title>References</title><ref id="vex027-B1"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Aiewsakun</surname><given-names>P.</given-names></name>, <name name-style="western"><surname>Katzourakis</surname><given-names>A.</given-names></name></person-group> (<year>2015</year>) ‘<article-title>Time dependency of foamy virus evolutionary rate estimates</article-title>’, <source>BMC Evolutionary Biology</source>, <volume>15</volume>: <fpage>119</fpage>.<pub-id pub-id-type="pmid">26111824</pub-id></mixed-citation></ref><ref id="vex027-B2"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Aiewsakun</surname><given-names>P.</given-names></name></person-group> (<year>2017</year>) ‘<article-title>Marine origin of retroviruses in the early Palaeozoic Era</article-title>’, <source>Nature Communications</source>, <volume>8</volume>: <fpage>13954</fpage>.</mixed-citation></ref><ref id="vex027-B3"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Bell</surname><given-names>M. A.</given-names></name>, <name name-style="western"><surname>Lloyd</surname><given-names>G.</given-names></name></person-group> (<year>2014</year>) ‘<article-title>strap: an R package for plotting phylogenies against stratigraphy and assessing their stratigraphic congruence</article-title>’, <source>Palaentology</source>, <volume>58</volume>/<issue>2</issue>: <fpage>379</fpage>–<lpage>89</lpage></mixed-citation></ref><ref id="vex027-B4"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Bergvall</surname><given-names>M.</given-names></name>, <name name-style="western"><surname>Melendy</surname><given-names>T.</given-names></name>, <name name-style="western"><surname>Archambault</surname><given-names>J.</given-names></name></person-group> (<year>2013</year>) ‘<article-title>The E1 proteins</article-title>’, <source>Virology</source>, <volume>445</volume>/<issue>1-2</issue>: <fpage>35</fpage>–<lpage>56</lpage>.<pub-id pub-id-type="pmid">24029589</pub-id></mixed-citation></ref><ref id="vex027-B5"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Bernard</surname><given-names>H. U.</given-names></name></person-group> (<year>1994</year>) ‘<article-title>Coevolution of papillomaviruses with human populations</article-title>’, <source>Trends Microbiology</source>, <volume>2</volume>/<issue>4</issue>: <fpage>140</fpage>–<lpage>3</lpage>.</mixed-citation></ref><ref id="vex027-B6"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Bernard</surname><given-names>H. U.</given-names></name></person-group><etal></etal> (<year>2010</year>) ‘<article-title>Classification of papillomaviruses (PVs) based on 189 PV types and proposal of taxonomic amendments</article-title>’, <source>Virology</source>, <volume>401</volume>/<issue>1</issue>: <fpage>70</fpage>–<lpage>9</lpage>.<pub-id pub-id-type="pmid">20206957</pub-id></mixed-citation></ref><ref id="vex027-B7"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Bouckaert</surname><given-names>R.</given-names></name></person-group><etal></etal> (<year>2014</year>) ‘<article-title>BEAST 2: a software platform for Bayesian evolutionary analysis</article-title>’, <source>PLOS Computational Biology</source>, <volume>10</volume>/<issue>4</issue>: <fpage>e1003537</fpage>.<pub-id pub-id-type="pmid">24722319</pub-id></mixed-citation></ref><ref id="vex027-B8"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Brooks</surname><given-names>D. R.</given-names></name>, <name name-style="western"><surname>Ferrao</surname><given-names>A. L.</given-names></name></person-group> (<year>2005</year>) ‘<article-title>The historical biogeography of co-evolution: emerging infectious diseases are evolutionary accidents waiting to happen</article-title>’, <source>Journal of Biogeography</source>, <volume>32</volume>: <fpage>1291</fpage>–<lpage>9</lpage>.</mixed-citation></ref><ref id="vex027-B9"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Buck</surname><given-names>C. B.</given-names></name></person-group><etal></etal> (<year>2016</year>) ‘<article-title>The ancient evolutionary history of polyomaviruses</article-title>’, <source>PLoS Pathogens</source>, <volume>12</volume>/<issue>4</issue>: <fpage>e1005574</fpage>.<pub-id pub-id-type="pmid">27093155</pub-id></mixed-citation></ref><ref id="vex027-B10"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>dos Reis</surname><given-names>M.</given-names></name></person-group><etal></etal> (<year>2015</year>) ‘<article-title>Uncertainty in the timing of origin of animals and the limits of precision in molecular timescales</article-title>’, <source>Current Biology</source>, <volume>25</volume>/<issue>22</issue>: <fpage>2939</fpage>–<lpage>50</lpage>.<pub-id pub-id-type="pmid">26603774</pub-id></mixed-citation></ref><ref id="vex027-B11"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Drummond</surname><given-names>A. J.</given-names></name></person-group><etal></etal> (<year>2006</year>) ‘<article-title>Relaxed phylogenetics and dating with confidence</article-title>’, <source>PLoS Biology</source>, <volume>4</volume>/<issue>5</issue>: <fpage>e88</fpage>.<pub-id pub-id-type="pmid">16683862</pub-id></mixed-citation></ref><ref id="vex027-B12"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Duchene</surname><given-names>S.</given-names></name>, <name name-style="western"><surname>Holmes</surname><given-names>E. C.</given-names></name>, <name name-style="western"><surname>Ho</surname><given-names>S. Y.</given-names></name></person-group> (<year>2014a</year>) ‘<article-title>Analyses of evolutionary dynamics in viruses are hindered by a time-dependent bias in rate estimates</article-title>’, <source>Proceedings Biological Sciences</source>, <volume>281</volume>/<issue>1786</issue>.</mixed-citation></ref><ref id="vex027-B13"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Duchene</surname><given-names>S.</given-names></name>, <name name-style="western"><surname>Lanfear</surname><given-names>R.</given-names></name>, <name name-style="western"><surname>Ho</surname><given-names>S. Y.</given-names></name></person-group> (<year>2014b</year>) ‘<article-title>The impact of calibration and clock-model choice on molecular estimates of divergence times</article-title>’, <source>Molecular Phylogenetics and Evolution</source>, <volume>78</volume>: <fpage>277</fpage>–<lpage>89</lpage>.<pub-id pub-id-type="pmid">24910154</pub-id></mixed-citation></ref><ref id="vex027-B14"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Gaynor</surname><given-names>A. M.</given-names></name></person-group><etal></etal> (<year>2015</year>) ‘<article-title>Identification of a novel papillomavirus in a Northern Fulmar (<italic>Fulmarus glacialis</italic>) with viral production in cartilage</article-title>’, <source>Veterinary Pathology</source>, <volume>52</volume>/<issue>3</issue>: <fpage>553</fpage>–<lpage>61</lpage>.<pub-id pub-id-type="pmid">25034110</pub-id></mixed-citation></ref><ref id="vex027-B15"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Gottschling</surname><given-names>M.</given-names></name></person-group><etal></etal> (<year>2011</year>) ‘<article-title>Quantifying the phylodynamic forces driving papillomavirus evolution</article-title>’, <source>Molecular Biology and Evolution</source>, <volume>28</volume>/<issue>7</issue>: <fpage>2101</fpage>–<lpage>13</lpage>.<pub-id pub-id-type="pmid">21285031</pub-id></mixed-citation></ref><ref id="vex027-B16"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Guindon</surname><given-names>S.</given-names></name></person-group><etal></etal> (<year>2010</year>) ‘<article-title>New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0</article-title>’, <source>Systematic Biology</source>, <volume>59</volume>/<issue>3</issue>: <fpage>307</fpage>–<lpage>21</lpage>.<pub-id pub-id-type="pmid">20525638</pub-id></mixed-citation></ref><ref id="vex027-B17"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Hedges</surname><given-names>S. B.</given-names></name>, <name name-style="western"><surname>Dudley</surname><given-names>J.</given-names></name>, <name name-style="western"><surname>Kumar</surname><given-names>S.</given-names></name></person-group> (<year>2006</year>) ‘<article-title>TimeTree: a public knowledge-base of divergence times among organisms</article-title>’, <source>Bioinformatics</source>, <volume>22</volume>/<issue>23</issue>: <fpage>2971</fpage>–<lpage>2</lpage>.<pub-id pub-id-type="pmid">17021158</pub-id></mixed-citation></ref><ref id="vex027-B18"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Hedges</surname><given-names>S. B.</given-names></name></person-group><etal></etal> (<year>2015</year>) ‘<article-title>Tree of life reveals clock-like speciation and diversification</article-title>’, <source>Molecular Biology and Evolution</source>, <volume>32</volume>/<issue>4</issue>: <fpage>835</fpage>–<lpage>45</lpage>.<pub-id pub-id-type="pmid">25739733</pub-id></mixed-citation></ref><ref id="vex027-B19"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Heled</surname><given-names>J.</given-names></name>, <name name-style="western"><surname>Drummond</surname><given-names>A. J.</given-names></name></person-group> (<year>2012</year>) ‘<article-title>Calibrated tree priors for relaxed phylogenetics and divergence time estimation</article-title>’, <source>Systematic Biology</source>, <volume>61</volume>/<issue>1</issue>: <fpage>138</fpage>–<lpage>49</lpage> doi: 10.1093/sysbio/syr087. PubMed PMID: 21856631; PubMed Central PMCID: PMCPMC3243734.<pub-id pub-id-type="pmid">21856631</pub-id></mixed-citation></ref><ref id="vex027-B20"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Herbst</surname><given-names>L. H.</given-names></name></person-group><etal></etal> (<year>2009</year>) ‘<article-title>Genomic characterization of two novel reptilian papillomaviruses, Chelonia mydas papillomavirus 1 and Caretta caretta papillomavirus 1</article-title>’, <source>Virology</source>, <volume>383</volume>/<issue>1</issue>: <fpage>131</fpage>–<lpage>5</lpage>. doi: 10.1016/j.virol.2008.09.022. PubMed PMID: 18973915.<pub-id pub-id-type="pmid">18973915</pub-id></mixed-citation></ref><ref id="vex027-B21"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Huelsenbeck</surname><given-names>J. P.</given-names></name>, <name name-style="western"><surname>Ronquist</surname><given-names>F.</given-names></name></person-group> (<year>2001</year>) ‘<article-title>MRBAYES: Bayesian inference of phylogenetic trees</article-title>’, <source>Bioinformatics</source>, <volume>17</volume>/<issue>8</issue>: <fpage>754</fpage>–<lpage>5</lpage> PubMed PMID: 11524383.<pub-id pub-id-type="pmid">11524383</pub-id></mixed-citation></ref><ref id="vex027-B22"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Jeckel</surname><given-names>S.</given-names></name></person-group><etal></etal> (<year>2003</year>) ‘ <article-title>Identification of the E9/E2C cDNA and functional characterization of the gene product reveal a new repressor of transcription and replication in cottontail rabbit papillomavirus</article-title>’, <source>Journal of Virology</source>, <volume>77</volume>/<issue>16</issue>: <fpage>8736</fpage>–<lpage>44</lpage>.<pub-id pub-id-type="pmid">12885893</pub-id></mixed-citation></ref><ref id="vex027-B23"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Katoh</surname><given-names>K.</given-names></name></person-group><etal></etal> (<year>2002</year>) ‘<article-title>MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform</article-title>’, <source>Nucleic Acids Research</source>, <volume>30</volume>/<issue>14</issue>: <fpage>3059</fpage>–<lpage>66</lpage>.<pub-id pub-id-type="pmid">12136088</pub-id></mixed-citation></ref><ref id="vex027-B24"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Katoh</surname><given-names>K.</given-names></name>, <name name-style="western"><surname>Standley</surname><given-names>D. M.</given-names></name></person-group> (<year>2013</year>) ‘<article-title>MAFFT multiple sequence alignment software version 7: improvements in performance and usability</article-title>’, <source>Molecular Biology and Evolution</source>, <volume>30</volume>/<issue>4</issue>: <fpage>772</fpage>–<lpage>80</lpage>.<pub-id pub-id-type="pmid">23329690</pub-id></mixed-citation></ref><ref id="vex027-B25"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Kearse</surname><given-names>M.</given-names></name></person-group><etal></etal> (<year>2012</year>) ‘<article-title>Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data</article-title>’, <source>Bioinformatics</source>, <volume>28</volume>/<issue>12</issue>: <fpage>1647</fpage>–<lpage>9</lpage>.<pub-id pub-id-type="pmid">22543367</pub-id></mixed-citation></ref><ref id="vex027-B26"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Lanfear</surname><given-names>R.</given-names></name></person-group><etal></etal> (<year>2017</year>) ‘<article-title>PartitionFinder 2: New Methods for Selecting Partitioned Models of Evolution for Molecular and Morphological Phylogenetic Analyses</article-title>’, <source>Molecular Biology and Evolution</source>, <volume>34</volume>/<issue>3</issue>: <fpage>772</fpage>–<lpage>3</lpage>.<pub-id pub-id-type="pmid">28013191</pub-id></mixed-citation></ref><ref id="vex027-B27"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Lange</surname><given-names>C. E.</given-names></name></person-group><etal></etal> (<year>2011</year>) ‘<article-title>Novel snake papillomavirus does not cluster with other non-mammalian papillomaviruses</article-title>’, <source>Virology Journal</source>, <volume>8</volume>: <fpage>436</fpage>.<pub-id pub-id-type="pmid">21910860</pub-id></mixed-citation></ref><ref id="vex027-B28"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Le</surname><given-names>S. Q.</given-names></name>, <name name-style="western"><surname>Gascuel</surname><given-names>O.</given-names></name></person-group> (<year>2008</year>) ‘<article-title>An improved general amino acid replacement matrix</article-title>’, <source>Molecular Biology and Evolution</source>, <volume>25</volume>/<issue>7</issue>: <fpage>1307</fpage>–<lpage>20</lpage>.<pub-id pub-id-type="pmid">18367465</pub-id></mixed-citation></ref><ref id="vex027-B29"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Lopez-Bueno</surname><given-names>A.</given-names></name></person-group><etal></etal> (<year>2016</year>) ‘<article-title>Concurrence of iridovirus, polyomavirus, and a unique member of a new group of fish papillomaviruses in lymphocystis disease-affected gilthead sea bream</article-title>’, <source>Journal of Virology</source>, <volume>90</volume>/<issue>19</issue>: <fpage>8768</fpage>–<lpage>79</lpage>.<pub-id pub-id-type="pmid">27440877</pub-id></mixed-citation></ref><ref id="vex027-B30"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Miller</surname><given-names>M. A.</given-names></name></person-group><etal></etal> (<year>2015</year>) ‘<article-title>A RESTful API for access to phylogenetic tools via the CIPRES Science Gateway</article-title>’, <source>Evolutionary Bioinformatics Online</source>, <volume>11</volume>: <fpage>43</fpage>–<lpage>8</lpage>.<pub-id pub-id-type="pmid">25861210</pub-id></mixed-citation></ref><ref id="vex027-B31"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Morin</surname><given-names>G.</given-names></name></person-group><etal></etal> (<year>2011</year>) ‘<article-title>A conserved amphipathic helix in the N-terminal regulatory region of the papillomavirus E1 helicase is required for efficient viral DNA replication</article-title>’, <source>Journal of Virology</source>, <volume>85</volume>/<issue>11</issue>: <fpage>5287</fpage>–<lpage>300</lpage>.<pub-id pub-id-type="pmid">21450828</pub-id></mixed-citation></ref><ref id="vex027-B32"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Muhire</surname><given-names>B. M.</given-names></name>, <name name-style="western"><surname>Varsani</surname><given-names>A.</given-names></name>, <name name-style="western"><surname>Martin</surname><given-names>D. P.</given-names></name></person-group> (<year>2014</year>) ‘<article-title>SDT: a virus classification tool based on pairwise sequence alignment and identity calculation</article-title>’, <source>PLoS One</source>, <volume>9</volume>/<issue>9</issue>: <fpage>e108277</fpage>.<pub-id pub-id-type="pmid">25259891</pub-id></mixed-citation></ref><ref id="vex027-B33"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Nylander</surname><given-names>J. A.</given-names></name></person-group><etal></etal> (<year>2008</year>) ‘<article-title>AWTY (are we there yet?): a system for graphical exploration of MCMC convergence in Bayesian phylogenetics</article-title>’, <source>Bioinformatics</source>, <volume>24</volume>/<issue>4</issue>: <fpage>581</fpage>–<lpage>3</lpage>.<pub-id pub-id-type="pmid">17766271</pub-id></mixed-citation></ref><ref id="vex027-B34"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Paradis</surname><given-names>E.</given-names></name>, <name name-style="western"><surname>Claude</surname><given-names>J.</given-names></name>, <name name-style="western"><surname>Strimmer</surname><given-names>K.</given-names></name></person-group> (<year>2004</year>) ‘<article-title>APE: analyses of phylogenetics and evolution in R language</article-title>’, <source>Bioinformatics</source>, <volume>20</volume>/<issue>2</issue>: <fpage>289</fpage>–<lpage>90</lpage>.<pub-id pub-id-type="pmid">14734327</pub-id></mixed-citation></ref><ref id="vex027-B35"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Pattengale</surname><given-names>N. D.</given-names></name></person-group><etal></etal> (<year>2010</year>) ‘<article-title>How many bootstrap replicates are necessary?</article-title> ’, <source>Journal of Computational Biology</source>, <volume>17</volume>/<issue>3</issue>: <fpage>337</fpage>–<lpage>54</lpage>.<pub-id pub-id-type="pmid">20377449</pub-id></mixed-citation></ref><ref id="vex027-B36"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Prum</surname><given-names>R. O.</given-names></name></person-group><etal></etal> (<year>2015</year>) ‘<article-title>A comprehensive phylogeny of birds (Aves) using targeted next-generation DNA sequencing</article-title>’, <source>Nature</source>, <volume>526</volume>/<issue>7574</issue>: <fpage>569</fpage>–<lpage>73</lpage>.<pub-id pub-id-type="pmid">26444237</pub-id></mixed-citation></ref><ref id="vex027-B37"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Rambaut</surname><given-names>A.</given-names></name></person-group><etal></etal> (<year>2016</year>) ‘<article-title>Exploring the temporal structure of heterochronous sequences using TempEst (formerly Path-O-Gen)</article-title>’, <source>Virus Evolution</source>, <volume>2</volume>/<issue>1</issue>: <fpage>vew007</fpage>.<pub-id pub-id-type="pmid">27774300</pub-id></mixed-citation></ref><ref id="vex027-B38"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Rector</surname><given-names>A.</given-names></name></person-group><etal></etal> (<year>2007</year>) ‘<article-title>Ancient papillomavirus-host co-speciation in Felidae</article-title>’, <source>Genome Biology</source>, <volume>8</volume>/<issue>4</issue>: <fpage>R57</fpage>.<pub-id pub-id-type="pmid">17430578</pub-id></mixed-citation></ref><ref id="vex027-B39"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Rector</surname><given-names>A.</given-names></name>, <name name-style="western"><surname>Van Ranst</surname><given-names>M.</given-names></name></person-group> (<year>2013</year>) ‘<article-title>Animal papillomaviruses</article-title>’, <source>Virology</source>, <volume>445</volume>/<issue>1-2</issue>: <fpage>213</fpage>–<lpage>23</lpage>.<pub-id pub-id-type="pmid">23711385</pub-id></mixed-citation></ref><ref id="vex027-B40"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Rey</surname><given-names>K.</given-names></name></person-group><etal></etal> (<year>2017</year>) ‘<article-title>Oxygen isotopes suggest elevated thermometabolism within multiple Permo-Triassic therapsid clades</article-title>’, <source>Elife</source>, <volume>6</volume>: <fpage>e28589</fpage>.<pub-id pub-id-type="pmid">28716184</pub-id></mixed-citation></ref><ref id="vex027-B41"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Ripley</surname><given-names>B. D.</given-names></name>, <name name-style="western"><surname>Thompson</surname><given-names>M.</given-names></name></person-group> (<year>1987</year>) ‘<article-title>Regression techniques for the detection of analytical bias</article-title>’, <source>Analyst</source>, <volume>10987112</volume>/<issue>4</issue>: <fpage>377</fpage>–<lpage>83</lpage>.</mixed-citation></ref><ref id="vex027-B42"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Roman</surname><given-names>A.</given-names></name>, <name name-style="western"><surname>Munger</surname><given-names>K.</given-names></name></person-group> (<year>2013</year>) ‘<article-title>The papillomavirus E7 proteins</article-title>’, <source>Virology</source>, <volume>445</volume>/<issue>1-2</issue>: <fpage>138</fpage>–<lpage>68</lpage>.<pub-id pub-id-type="pmid">23731972</pub-id></mixed-citation></ref><ref id="vex027-B43"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Ronquist</surname><given-names>F.</given-names></name></person-group><etal></etal> (<year>2012</year>) ‘<article-title>MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space</article-title>’, <source>Systematic Biology</source>, <volume>61</volume>/<issue>3</issue>: <fpage>539</fpage>–<lpage>42</lpage>.<pub-id pub-id-type="pmid">22357727</pub-id></mixed-citation></ref><ref id="vex027-B44"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Ronquist</surname><given-names>F.</given-names></name>, <name name-style="western"><surname>Huelsenbeck</surname><given-names>J. P.</given-names></name></person-group> (<year>2003</year>) ‘<article-title>MrBayes 3: Bayesian phylogenetic inference under mixed models</article-title>’, <source>Bioinformatics</source>, <volume>19</volume>/<issue>12</issue>: <fpage>1572</fpage>–<lpage>4</lpage>.<pub-id pub-id-type="pmid">12912839</pub-id></mixed-citation></ref><ref id="vex027-B45"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Ruboyianes</surname><given-names>R.</given-names></name>, <name name-style="western"><surname>Worobey</surname><given-names>M.</given-names></name></person-group> (<year>2016</year>) ‘<article-title>Foamy-like endogenous retroviruses are extensive and abundant in teleosts</article-title>’, <source>Virus Evolution</source>, <volume>2</volume>/<issue>2</issue>: <fpage>vew032</fpage>.<pub-id pub-id-type="pmid">28058112</pub-id></mixed-citation></ref><ref id="vex027-B46"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Schmidt</surname><given-names>H. A.</given-names></name></person-group><etal></etal> (<year>2002</year>) ‘<article-title>TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing</article-title>’, <source>Bioinformatics</source>, <volume>18</volume>/<issue>3</issue>: <fpage>502</fpage>–<lpage>4</lpage>.<pub-id pub-id-type="pmid">11934758</pub-id></mixed-citation></ref><ref id="vex027-B47"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Shah</surname><given-names>S. D.</given-names></name>, <name name-style="western"><surname>Doorbar</surname><given-names>J.</given-names></name>, <name name-style="western"><surname>Goldstein</surname><given-names>R. A.</given-names></name></person-group> (<year>2010</year>) ‘<article-title>Analysis of host-parasite incongruence in papillomavirus evolution using importance sampling</article-title>’, <source>Molecular Biology and Evolution</source>, <volume>27</volume>/<issue>6</issue>: <fpage>1301</fpage>–<lpage>14</lpage>.<pub-id pub-id-type="pmid">20093429</pub-id></mixed-citation></ref><ref id="vex027-B48"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Shimodaira</surname><given-names>H.</given-names></name></person-group> (<year>2002</year>) ‘<article-title>An approximately unbiased test of phylogenetic tree selection</article-title>’, <source>Systematic Biology</source>, <volume>51</volume>/<issue>3</issue>: <fpage>492</fpage>–<lpage>508</lpage>.<pub-id pub-id-type="pmid">12079646</pub-id></mixed-citation></ref><ref id="vex027-B49"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Shimodaira</surname><given-names>H.</given-names></name>, <name name-style="western"><surname>Hasegawa</surname><given-names>M.</given-names></name></person-group> (<year>2001</year>) ‘<article-title>CONSEL: for assessing the confidence of phylogenetic tree selection</article-title>’, <source>Bioinformatics</source>, <volume>17</volume>/<issue>12</issue>: <fpage>1246</fpage>–<lpage>7</lpage>.<pub-id pub-id-type="pmid">11751242</pub-id></mixed-citation></ref><ref id="vex027-B50"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Stamatakis</surname><given-names>A.</given-names></name></person-group> (<year>2014</year>) ‘<article-title>RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies</article-title>’, <source>Bioinformatics</source>, <volume>30</volume>/<issue>9</issue>: <fpage>1312</fpage>–<lpage>3</lpage>.<pub-id pub-id-type="pmid">24451623</pub-id></mixed-citation></ref><ref id="vex027-B51"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Tachezy</surname><given-names>R.</given-names></name>, <name name-style="western"><surname>Rector</surname><given-names>A.</given-names></name>, <name name-style="western"><surname>Havelkova</surname><given-names>M.</given-names></name>, <name name-style="western"><surname>Wollants</surname><given-names>E.</given-names></name>, <name name-style="western"><surname>Fiten</surname><given-names>P.</given-names></name>, <name name-style="western"><surname>Opdenakker</surname><given-names>G.</given-names></name></person-group><etal></etal> (<year>2002</year>) ‘<article-title>Avian papillomaviruses: the parrot Psittacus erithacus papillomavirus (PePV) genome has a unique organization of the early protein region and is phylogenetically related to the chaffinch papillomavirus</article-title>’, <source>BMC Microbiology</source>, <volume>2</volume>: <fpage>19</fpage>.<pub-id pub-id-type="pmid">12110158</pub-id></mixed-citation></ref><ref id="vex027-B52"><mixed-citation publication-type="other"><collab>R Core Team</collab>. (<year>2014</year>) <italic>R: A Language and Environment for Statistical Computing</italic>. Vienna, Austria: R Foundation Statistical Computing1:409.</mixed-citation></ref><ref id="vex027-B53"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Terai</surname><given-names>M.</given-names></name>, <name name-style="western"><surname>DeSalle</surname><given-names>R.</given-names></name>, <name name-style="western"><surname>Burk</surname><given-names>R. D.</given-names></name></person-group> (<year>2002</year>) ‘<article-title>Lack of canonical E6 and E7 open reading frames in bird papillomaviruses: Fringilla coelebs papillomavirus and Psittacus erithacus timneh papillomavirus</article-title>’, <source>Journal of Virology</source>, <volume>76</volume>/<issue>19</issue>: <fpage>10020</fpage>–<lpage>3</lpage>.<pub-id pub-id-type="pmid">12208979</pub-id></mixed-citation></ref><ref id="vex027-B54"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Van Doorslaer</surname><given-names>K.</given-names></name></person-group> (<year>2013</year>) ‘<article-title>Evolution of the papillomaviridae</article-title>’, <source>Virology</source>, <volume>445</volume>/<issue>1-2</issue>: <fpage>11</fpage>–<lpage>20</lpage>.<pub-id pub-id-type="pmid">23769415</pub-id></mixed-citation></ref><ref id="vex027-B55"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Van Doorslaer</surname><given-names>K.</given-names></name></person-group><etal></etal> (<year>2009</year>) ‘<article-title>Identification of unusual E6 and E7 proteins within avian papillomaviruses: cellular localization, biophysical characterization, and phylogenetic analysis</article-title>’, <source>Journal of Virology</source>, <volume>83</volume>/<issue>17</issue>: <fpage>8759</fpage>–<lpage>70</lpage>.<pub-id pub-id-type="pmid">19553340</pub-id></mixed-citation></ref><ref id="vex027-B56"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Van Doorslaer</surname><given-names>K.</given-names></name></person-group><etal></etal> (<year>2013</year>) ‘<article-title>The Papillomavirus Episteme: a central resource for papillomavirus sequence data and analysis</article-title>’, <source>Nucleic Acids Research</source>, <volume>41(Database issue)</volume>: <fpage>D571</fpage>–<lpage>8</lpage>.<pub-id pub-id-type="pmid">23093593</pub-id></mixed-citation></ref><ref id="vex027-B57"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Van Doorslaer</surname><given-names>K.</given-names></name></person-group><etal></etal> (<year>2017</year>) ‘<article-title>The papillomavirus episteme: a major update to the papillomavirus sequence database</article-title>’, <source>Nucleic Acids Research</source>, <volume>45</volume>/<issue>D1</issue>: <fpage>D499</fpage>–<lpage>506</lpage>.<pub-id pub-id-type="pmid">28053164</pub-id></mixed-citation></ref><ref id="vex027-B58"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Vande Pol</surname><given-names>S. B.</given-names></name>, <name name-style="western"><surname>Klingelhutz</surname><given-names>A. J.</given-names></name></person-group> (<year>2013</year>) ‘<article-title>Papillomavirus E6 oncoproteins</article-title>’, <source>Virology</source>, <volume>445</volume>/<issue>1-2</issue>: <fpage>115</fpage>–<lpage>37</lpage>.<pub-id pub-id-type="pmid">23711382</pub-id></mixed-citation></ref><ref id="vex027-B59"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Varsani</surname><given-names>A.</given-names></name></person-group><etal></etal> (<year>2014</year>) ‘<article-title>A novel papillomavirus in Adelie penguin (<italic>Pygoscelis adeliae</italic>) faeces sampled at the Cape Crozier colony, Antarctica</article-title>’, <source>Journal of General Virology</source>, <volume>95(Pt 6)</volume>: <fpage>1352</fpage>–<lpage>65</lpage>.<pub-id pub-id-type="pmid">24686913</pub-id></mixed-citation></ref><ref id="vex027-B60"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Varsani</surname><given-names>A.</given-names></name></person-group> (<year>2015</year>) ‘<article-title>Identification of an avian polyomavirus associated with Adelie penguins (<italic>Pygoscelis adeliae</italic>)</article-title>’, <source>Journal of General Virology</source>, <volume>96(Pt 4)</volume>: <fpage>851</fpage>–<lpage>7</lpage>.<pub-id pub-id-type="pmid">25537375</pub-id></mixed-citation></ref><ref id="vex027-B61"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Zanier</surname><given-names>K.</given-names></name></person-group><etal></etal> (<year>2013</year>) ‘<article-title>Structural basis for hijacking of cellular LxxLL motifs by papillomavirus E6 oncoproteins</article-title>’, <source>Science</source>, <volume>339</volume>/<issue>6120</issue>: <fpage>694</fpage>–<lpage>8</lpage>.<pub-id pub-id-type="pmid">23393263</pub-id></mixed-citation></ref></ref-list></back></pmc></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Wicri/Asie/explor/AustralieFrV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 002D41 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 002D41 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Wicri/Asie
|area= AustralieFrV1
|flux= Pmc
|étape= Corpus
|type= RBID
|clé= PMC:5632515
|texte= Unique genome organization of non-mammalian papillomaviruses provides insights into the evolution of viral early proteins
}}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:29026649" \
| HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \
| NlmPubMed2Wicri -a AustralieFrV1
| This area was generated with Dilib version V0.6.33. |  |