Role of the host restriction factor APOBEC3 on papillomavirus evolution
Identifieur interne : 000141 ( Pmc/Corpus ); précédent : 000140; suivant : 000142Role of the host restriction factor APOBEC3 on papillomavirus evolution
Auteurs : Cody J. Warren ; Koenraad Van Doorslaer ; Ahwan Pandey ; Joaquin M. Espinosa ; Dohun PyeonSource :
- Virus evolution [ 2057-1577 ] ; 2015.
Abstract
More than 270 different types of papillomaviruses have been discovered in a wide array of animal species. Despite the great diversity of papillomaviruses, little is known about the evolutionary processes that drive host tropism and the emergence of oncogenic genotypes. Although host defense mechanisms have evolved to interfere with various aspects of a virus life cycle, viruses have also coevolved copious strategies to avoid host antiviral restriction. Our and other studies have shown that the cytidine deaminase APOBEC3 family members edit HPV genomes and restrict virus infectivity. Thus, we hypothesized that host restriction by APOBEC3 served as selective pressure during papillomavirus evolution. To test this hypothesis, we analyzed the relative abundance of all dinucleotide sequences in full-length genomes of 274 papillomavirus types documented in the Papillomavirus Episteme database (PaVE). Here, we report that TC dinucleotides, the preferred target sequence of several human APOBEC3 proteins (hA3A, hA3B, hA3F, and hA3H), are highly depleted in papillomavirus genomes. Given that HPV infection is highly tissue-specific, the expression levels of APOBEC3 family members were analyzed. The basal expression levels of all APOBEC3 isoforms, excluding hA3B, are significantly higher in mucosal skin compared with cutaneous skin. Interestingly, we reveal that
Url:
DOI: 10.1093/ve/vev015
PubMed: NONE
PubMed Central: 4999249
Links to Exploration step
PMC:4999249Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Role of the host restriction factor APOBEC3 on papillomavirus evolution</title><author><name sortKey="Warren, Cody J" sort="Warren, Cody J" uniqKey="Warren C" first="Cody J." last="Warren">Cody J. Warren</name><affiliation><nlm:aff id="A1">Department of Immunology and Microbiology, University of Colorado School of Medicine, Aurora, CO, USA</nlm:aff></affiliation></author><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="A2">DNA Tumor Virus Section, Laboratory of Viral Diseases, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, MD, USA</nlm:aff></affiliation></author><author><name sortKey="Pandey, Ahwan" sort="Pandey, Ahwan" uniqKey="Pandey A" first="Ahwan" last="Pandey">Ahwan Pandey</name><affiliation><nlm:aff id="A3">Department of Pharmacology, University of Colorado School of Medicine, Aurora, CO, USA</nlm:aff></affiliation><affiliation><nlm:aff id="A4">Linda Crnic Institute for Down Syndrome, University of Colorado School of Medicine, Aurora, CO, USA</nlm:aff></affiliation></author><author><name sortKey="Espinosa, Joaquin M" sort="Espinosa, Joaquin M" uniqKey="Espinosa J" first="Joaquin M." last="Espinosa">Joaquin M. Espinosa</name><affiliation><nlm:aff id="A3">Department of Pharmacology, University of Colorado School of Medicine, Aurora, CO, USA</nlm:aff></affiliation><affiliation><nlm:aff id="A4">Linda Crnic Institute for Down Syndrome, University of Colorado School of Medicine, Aurora, CO, USA</nlm:aff></affiliation></author><author><name sortKey="Pyeon, Dohun" sort="Pyeon, Dohun" uniqKey="Pyeon D" first="Dohun" last="Pyeon">Dohun Pyeon</name><affiliation><nlm:aff id="A1">Department of Immunology and Microbiology, University of Colorado School of Medicine, Aurora, CO, USA</nlm:aff></affiliation><affiliation><nlm:aff id="A5">Division of Infectious Diseases, Department of Medicine, University of Colorado School of Medicine, Aurora, CO, USA</nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmc">4999249</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4999249</idno><idno type="RBID">PMC:4999249</idno><idno type="doi">10.1093/ve/vev015</idno><idno type="pmid">NONE</idno><date when="2015">2015</date><idno type="wicri:Area/Pmc/Corpus">000141</idno><idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000141</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Role of the host restriction factor APOBEC3 on papillomavirus evolution</title><author><name sortKey="Warren, Cody J" sort="Warren, Cody J" uniqKey="Warren C" first="Cody J." last="Warren">Cody J. Warren</name><affiliation><nlm:aff id="A1">Department of Immunology and Microbiology, University of Colorado School of Medicine, Aurora, CO, USA</nlm:aff></affiliation></author><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="A2">DNA Tumor Virus Section, Laboratory of Viral Diseases, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, MD, USA</nlm:aff></affiliation></author><author><name sortKey="Pandey, Ahwan" sort="Pandey, Ahwan" uniqKey="Pandey A" first="Ahwan" last="Pandey">Ahwan Pandey</name><affiliation><nlm:aff id="A3">Department of Pharmacology, University of Colorado School of Medicine, Aurora, CO, USA</nlm:aff></affiliation><affiliation><nlm:aff id="A4">Linda Crnic Institute for Down Syndrome, University of Colorado School of Medicine, Aurora, CO, USA</nlm:aff></affiliation></author><author><name sortKey="Espinosa, Joaquin M" sort="Espinosa, Joaquin M" uniqKey="Espinosa J" first="Joaquin M." last="Espinosa">Joaquin M. Espinosa</name><affiliation><nlm:aff id="A3">Department of Pharmacology, University of Colorado School of Medicine, Aurora, CO, USA</nlm:aff></affiliation><affiliation><nlm:aff id="A4">Linda Crnic Institute for Down Syndrome, University of Colorado School of Medicine, Aurora, CO, USA</nlm:aff></affiliation></author><author><name sortKey="Pyeon, Dohun" sort="Pyeon, Dohun" uniqKey="Pyeon D" first="Dohun" last="Pyeon">Dohun Pyeon</name><affiliation><nlm:aff id="A1">Department of Immunology and Microbiology, University of Colorado School of Medicine, Aurora, CO, USA</nlm:aff></affiliation><affiliation><nlm:aff id="A5">Division of Infectious Diseases, Department of Medicine, University of Colorado School of Medicine, Aurora, CO, USA</nlm:aff></affiliation></author></analytic><series><title level="j">Virus evolution</title><idno type="eISSN">2057-1577</idno><imprint><date when="2015">2015</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><p id="P1">More than 270 different types of papillomaviruses have been discovered in a wide array of animal species. Despite the great diversity of papillomaviruses, little is known about the evolutionary processes that drive host tropism and the emergence of oncogenic genotypes. Although host defense mechanisms have evolved to interfere with various aspects of a virus life cycle, viruses have also coevolved copious strategies to avoid host antiviral restriction. Our and other studies have shown that the cytidine deaminase APOBEC3 family members edit HPV genomes and restrict virus infectivity. Thus, we hypothesized that host restriction by APOBEC3 served as selective pressure during papillomavirus evolution. To test this hypothesis, we analyzed the relative abundance of all dinucleotide sequences in full-length genomes of 274 papillomavirus types documented in the Papillomavirus Episteme database (PaVE). Here, we report that TC dinucleotides, the preferred target sequence of several human APOBEC3 proteins (hA3A, hA3B, hA3F, and hA3H), are highly depleted in papillomavirus genomes. Given that HPV infection is highly tissue-specific, the expression levels of APOBEC3 family members were analyzed. The basal expression levels of all APOBEC3 isoforms, excluding hA3B, are significantly higher in mucosal skin compared with cutaneous skin. Interestingly, we reveal that <italic>Alphapapillomaviruses</italic> (alpha-PVs), a majority of which infects anogenital mucosa, display the most dramatic reduction in TC dinucleotide content. Computer modeling and reconstruction of ancestral alpha-PV genomes suggest that TC depletion occurred after the alpha-PVs diverged from their most recent common ancestor. In addition, we found that TC depletion in alpha-PVs is greatly affected by protein coding potential. Taken together, our results suggest that PVs replicating in tissues with high APOBEC3 levels may have evolved to evade restriction by selecting for variants that contain reduced APOBEC3 target sites in their genomes.</p></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Ahasan, Mm" uniqKey="Ahasan M">MM Ahasan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Barker, D" uniqKey="Barker D">D Barker</name></author><author><name sortKey="Meade, A" uniqKey="Meade A">A Meade</name></author><author><name sortKey="Pagel, M" uniqKey="Pagel M">M Pagel</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Betts, Mj" uniqKey="Betts M">MJ Betts</name></author><author><name sortKey="Russell, Rb" uniqKey="Russell R">RB Russell</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bogerd, Hp" uniqKey="Bogerd H">HP Bogerd</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bravo, Ig" uniqKey="Bravo I">IG Bravo</name></author><author><name sortKey="Felez Sanchez, M" uniqKey="Felez Sanchez M">M Félez-Sánchez</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bravo, Ig" uniqKey="Bravo I">IG Bravo</name></author><author><name sortKey="Felez Sanchez, M" uniqKey="Felez Sanchez M">M Félez-Sánchez</name></author><author><name sortKey="Muller, M" uniqKey="Muller M">M Muller</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Burns, Mb" uniqKey="Burns M">MB Burns</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Burns, Mb" uniqKey="Burns M">MB Burns</name></author><author><name sortKey="Temiz, Na" uniqKey="Temiz N">NA Temiz</name></author><author><name sortKey="Harris, Rs" uniqKey="Harris R">RS Harris</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chen, H" uniqKey="Chen H">H Chen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chen, Z" uniqKey="Chen Z">Z Chen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chiu, Y L" uniqKey="Chiu Y">Y-L Chiu</name></author><author><name sortKey="Greene, Wc" uniqKey="Greene W">WC Greene</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chiu, Y L" uniqKey="Chiu Y">Y-L Chiu</name></author><author><name sortKey="Greene, Wc" uniqKey="Greene W">WC Greene</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Cladel, Nm" uniqKey="Cladel N">NM Cladel</name></author><author><name sortKey="Bertotto, A" uniqKey="Bertotto A">A Bertotto</name></author><author><name sortKey="Christensen, Nd" uniqKey="Christensen N">ND Christensen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Cullen, Br" uniqKey="Cullen B">BR Cullen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Darriba, D" uniqKey="Darriba D">D Darriba</name></author><author><name sortKey="Taboada, Gl" uniqKey="Taboada G">GL Taboada</name></author><author><name sortKey="Doallo, R" uniqKey="Doallo R">R Doallo</name></author><author><name sortKey="Posada, D" uniqKey="Posada D">D Posada</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="De Villiers, E M" uniqKey="De Villiers E">E-M de Villiers</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Duggal, Nk" uniqKey="Duggal N">NK Duggal</name></author><author><name sortKey="Emerman, M" uniqKey="Emerman M">M Emerman</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Esnault, Cc" uniqKey="Esnault C">CC Esnault</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Felez Sanchez, M" uniqKey="Felez Sanchez M">M Félez-Sánchez</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Friedman, Ji" uniqKey="Friedman J">JI Friedman</name></author><author><name sortKey="Stivers, Jt" uniqKey="Stivers J">JT Stivers</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Harris, Rs" uniqKey="Harris R">RS Harris</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Harris, Rs" uniqKey="Harris R">RS Harris</name></author><author><name sortKey="Liddament, Mt" uniqKey="Liddament M">MT Liddament</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hasan, Ua" uniqKey="Hasan U">UA Hasan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Henderson, S" uniqKey="Henderson S">S Henderson</name></author><author><name sortKey="Chakravarthy, A" uniqKey="Chakravarthy A">A Chakravarthy</name></author><author><name sortKey="Fenton, T" uniqKey="Fenton T">T Fenton</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hoelzer, K" uniqKey="Hoelzer K">K Hoelzer</name></author><author><name sortKey="Shackelton, La" uniqKey="Shackelton L">LA Shackelton</name></author><author><name sortKey="Parrish, Cr" uniqKey="Parrish C">CR Parrish</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hultquist, Jf" uniqKey="Hultquist J">JF Hultquist</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jarmuz, A" uniqKey="Jarmuz A">A Jarmuz</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jern, P" uniqKey="Jern P">P Jern</name></author><author><name sortKey="Russell, Ra" uniqKey="Russell R">RA Russell</name></author><author><name sortKey="Pathak, Vk" uniqKey="Pathak V">VK Pathak</name></author><author><name sortKey="Coffin, Jm" uniqKey="Coffin J">JM Coffin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Johnson, We" uniqKey="Johnson W">WE Johnson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Karlin, S" uniqKey="Karlin S">S Karlin</name></author><author><name sortKey="Doerfler, W" uniqKey="Doerfler W">W Doerfler</name></author><author><name sortKey="Cardon, Lr" uniqKey="Cardon L">LR Cardon</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Katoh, K" uniqKey="Katoh K">K Katoh</name></author><author><name sortKey="Standley, Dm" uniqKey="Standley D">DM Standley</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kim, E Y" uniqKey="Kim E">E-Y Kim</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lackey, L" uniqKey="Lackey L">L Lackey</name></author><author><name sortKey="Law, Ek" uniqKey="Law E">EK Law</name></author><author><name sortKey="Brown, Wl" uniqKey="Brown W">WL Brown</name></author><author><name sortKey="Harris, Rs" uniqKey="Harris R">RS Harris</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Larue, Rs" uniqKey="Larue R">RS LaRue</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Liang, G" uniqKey="Liang G">G Liang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Liddament, Mt" uniqKey="Liddament M">MT Liddament</name></author><author><name sortKey="Brown, Wl" uniqKey="Brown W">WL Brown</name></author><author><name sortKey="Schumacher, Aj" uniqKey="Schumacher A">AJ Schumacher</name></author><author><name sortKey="Harris, Rs" uniqKey="Harris R">RS Harris</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lonsdale, J" uniqKey="Lonsdale J">J Lonsdale</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Love, Mi" uniqKey="Love M">MI Love</name></author><author><name sortKey="Huber, W" uniqKey="Huber W">W Huber</name></author><author><name sortKey="Anders, S" uniqKey="Anders S">S Anders</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Miller, M" uniqKey="Miller M">M Miller</name></author><author><name sortKey="Pfeiffer, W" uniqKey="Pfeiffer W">W Pfeiffer</name></author><author><name sortKey="Schwartz, T" uniqKey="Schwartz T">T Schwartz</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Minkah, N" uniqKey="Minkah N">N Minkah</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Mulder, Lcf" uniqKey="Mulder L">LCF Mulder</name></author><author><name sortKey="Harari, A" uniqKey="Harari A">A Harari</name></author><author><name sortKey="Simon, V" uniqKey="Simon V">V Simon</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Munk, C" uniqKey="Munk C">C Münk</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Neil, Sjd" uniqKey="Neil S">SJD Neil</name></author><author><name sortKey="Zang, T" uniqKey="Zang T">T Zang</name></author><author><name sortKey="Bieniasz, Pd" uniqKey="Bieniasz P">PD Bieniasz</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nik Zainal, S" uniqKey="Nik Zainal S">S Nik-Zainal</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nobre, Rj" uniqKey="Nobre R">RJ Nobre</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nylander, Jaa" uniqKey="Nylander J">JAA Nylander</name></author><author><name sortKey="Wilgenbusch, Jc" uniqKey="Wilgenbusch J">JC Wilgenbusch</name></author><author><name sortKey="Warren, Dl" uniqKey="Warren D">DL Warren</name></author><author><name sortKey="Swofford, Dl" uniqKey="Swofford D">DL Swofford</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ong, Ck" uniqKey="Ong C">CK Ong</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ooms, M" uniqKey="Ooms M">M Ooms</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Pagel, M" uniqKey="Pagel M">M Pagel</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Pagel, M" uniqKey="Pagel M">M Pagel</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Pham, P" uniqKey="Pham P">P Pham</name></author><author><name sortKey="Bransteitter, R" uniqKey="Bransteitter R">R Bransteitter</name></author><author><name sortKey="Petruska, J" uniqKey="Petruska J">J Petruska</name></author><author><name sortKey="Goodman, Mf" uniqKey="Goodman M">MF Goodman</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="Shackelton, La" uniqKey="Shackelton L">LA Shackelton</name></author><author><name sortKey="Parrish, Cr" uniqKey="Parrish C">CR Parrish</name></author><author><name sortKey="Holmes, Ec" uniqKey="Holmes E">EC Holmes</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sheehy, Am" uniqKey="Sheehy A">AM Sheehy</name></author><author><name sortKey="Gaddis, Nc" uniqKey="Gaddis N">NC Gaddis</name></author><author><name sortKey="Choi, Jd" uniqKey="Choi J">JD Choi</name></author><author><name sortKey="Malim, Mh" uniqKey="Malim M">MH Malim</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="Suspene, R" uniqKey="Suspene R">R Suspène</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Suspene, R" uniqKey="Suspene R">R Suspène</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Turelli, P" uniqKey="Turelli P">P Turelli</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Upadhyay, M" uniqKey="Upadhyay M">M Upadhyay</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="Vartanian, J P" uniqKey="Vartanian J">J-P Vartanian</name></author><author><name sortKey="Guetard, D" uniqKey="Guetard D">D Guétard</name></author><author><name sortKey="Henry, M" uniqKey="Henry M">M Henry</name></author><author><name sortKey="Wain Hobson, S" uniqKey="Wain Hobson S">S Wain-Hobson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wang, Z" uniqKey="Wang Z">Z Wang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Warren, Cj" uniqKey="Warren C">CJ Warren</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wood, N" uniqKey="Wood N">N Wood</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yang, B" uniqKey="Yang B">B Yang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yu, Q" uniqKey="Yu Q">Q Yu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yu, Q" uniqKey="Yu Q">Q Yu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhang, H" uniqKey="Zhang H">H Zhang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zheng, Y H" uniqKey="Zheng Y">Y-H Zheng</name></author></analytic></biblStruct></listBibl></div1></back></TEI><pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<pmc-dir>properties manuscript</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-journal-id">101664675</journal-id><journal-id journal-id-type="pubmed-jr-id">44168</journal-id><journal-id journal-id-type="nlm-ta">Virus Evol</journal-id><journal-id journal-id-type="iso-abbrev">Virus Evol</journal-id><journal-title-group><journal-title>Virus evolution</journal-title></journal-title-group><issn pub-type="epub">2057-1577</issn></journal-meta><article-meta><article-id pub-id-type="pmc">4999249</article-id><article-id pub-id-type="doi">10.1093/ve/vev015</article-id><article-id pub-id-type="manuscript">NIHMS778642</article-id><article-categories><subj-group subj-group-type="heading"><subject>Article</subject></subj-group></article-categories><title-group><article-title>Role of the host restriction factor APOBEC3 on papillomavirus evolution</article-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Warren</surname><given-names>Cody J.</given-names></name><xref ref-type="aff" rid="A1">1</xref><xref rid="FN2" ref-type="author-notes">*</xref></contrib><contrib contrib-type="author"><name><surname>Van Doorslaer</surname><given-names>Koenraad</given-names></name><xref ref-type="aff" rid="A2">2</xref><xref rid="FN2" ref-type="author-notes">*</xref><xref rid="FN3" ref-type="author-notes">†</xref></contrib><contrib contrib-type="author"><name><surname>Pandey</surname><given-names>Ahwan</given-names></name><xref ref-type="aff" rid="A3">3</xref><xref ref-type="aff" rid="A4">4</xref></contrib><contrib contrib-type="author"><name><surname>Espinosa</surname><given-names>Joaquin M.</given-names></name><xref ref-type="aff" rid="A3">3</xref><xref ref-type="aff" rid="A4">4</xref></contrib><contrib contrib-type="author"><name><surname>Pyeon</surname><given-names>Dohun</given-names></name><xref ref-type="aff" rid="A1">1</xref><xref ref-type="aff" rid="A5">5</xref><xref rid="FN1" ref-type="author-notes">‡</xref></contrib></contrib-group><aff id="A1"><label>1</label>Department of Immunology and Microbiology, University of Colorado School of Medicine, Aurora, CO, USA</aff><aff id="A2"><label>2</label>DNA Tumor Virus Section, Laboratory of Viral Diseases, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, MD, USA</aff><aff id="A3"><label>3</label>Department of Pharmacology, University of Colorado School of Medicine, Aurora, CO, USA</aff><aff id="A4"><label>4</label>Linda Crnic Institute for Down Syndrome, University of Colorado School of Medicine, Aurora, CO, USA</aff><aff id="A5"><label>5</label>Division of Infectious Diseases, Department of Medicine, University of Colorado School of Medicine, Aurora, CO, USA</aff><author-notes><corresp id="FN1"><label>‡</label>Corresponding author: <email>Dohun.Pyeon@ucdenver.edu</email></corresp><fn id="FN2" fn-type="equal"><label>*</label><p>These authors contributed equally to this work.</p></fn><fn id="FN3"><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></author-notes><pub-date pub-type="nihms-submitted"><day>20</day><month>4</month><year>2016</year></pub-date><pub-date pub-type="epub"><day>1</day><month>1</month><year>2015</year></pub-date><pub-date pub-type="ppub"><month>1</month><year>2015</year></pub-date><pub-date pub-type="pmc-release"><day>25</day><month>8</month><year>2016</year></pub-date><volume>1</volume><issue>1</issue><elocation-id>vev015</elocation-id><permissions><license license-type="open-access"><license-p>This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by-nc/4.0/">http://creativecommons.org/licenses/by-nc/4.0/</ext-link>), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited. For commercial re-use, please contact <email>journals.permissions@oup.com</email></license-p></license></permissions><self-uri xlink:href="http://ve.oxfordjournals.org/content/1/1/vev015"></self-uri><abstract><p id="P1">More than 270 different types of papillomaviruses have been discovered in a wide array of animal species. Despite the great diversity of papillomaviruses, little is known about the evolutionary processes that drive host tropism and the emergence of oncogenic genotypes. Although host defense mechanisms have evolved to interfere with various aspects of a virus life cycle, viruses have also coevolved copious strategies to avoid host antiviral restriction. Our and other studies have shown that the cytidine deaminase APOBEC3 family members edit HPV genomes and restrict virus infectivity. Thus, we hypothesized that host restriction by APOBEC3 served as selective pressure during papillomavirus evolution. To test this hypothesis, we analyzed the relative abundance of all dinucleotide sequences in full-length genomes of 274 papillomavirus types documented in the Papillomavirus Episteme database (PaVE). Here, we report that TC dinucleotides, the preferred target sequence of several human APOBEC3 proteins (hA3A, hA3B, hA3F, and hA3H), are highly depleted in papillomavirus genomes. Given that HPV infection is highly tissue-specific, the expression levels of APOBEC3 family members were analyzed. The basal expression levels of all APOBEC3 isoforms, excluding hA3B, are significantly higher in mucosal skin compared with cutaneous skin. Interestingly, we reveal that <italic>Alphapapillomaviruses</italic> (alpha-PVs), a majority of which infects anogenital mucosa, display the most dramatic reduction in TC dinucleotide content. Computer modeling and reconstruction of ancestral alpha-PV genomes suggest that TC depletion occurred after the alpha-PVs diverged from their most recent common ancestor. In addition, we found that TC depletion in alpha-PVs is greatly affected by protein coding potential. Taken together, our results suggest that PVs replicating in tissues with high APOBEC3 levels may have evolved to evade restriction by selecting for variants that contain reduced APOBEC3 target sites in their genomes.</p></abstract><kwd-group><kwd>papillomavirus</kwd><kwd>APOBEC3</kwd><kwd>restriction factor</kwd><kwd>coevolution</kwd></kwd-group></article-meta></front><body><sec id="S1"><title>1. Introduction</title><p id="P2">Papillomaviruses (PVs) are non-enveloped, double-stranded DNA viruses that infect keratinocytes and lead to benign and malignant tumors. Over 270 types of PVs have been discovered in various host species including mammals, reptiles, and birds (<xref rid="R61" ref-type="bibr">Van Doorslaer et al., 2013</xref>). It is believed that PVs have evolved for hundreds of millions of years and successfully established the largest and most diversified family of vertebrate viruses.</p><p id="P3">Of about 180 different types of human papillomaviruses (HPVs) identified so far, high-risk genotypes such as HPV16 and HPV18 are causally associated with multiple human cancers including cervical and head and neck cancers. All high-risk HPV types are classified in a single genus, <italic>Alphapapillomaviridae</italic>, which diverged from a common ancestor about 75 million years ago (<xref rid="R52" ref-type="bibr">Rector et al. 2007</xref>; <xref rid="R61" ref-type="bibr">Van Doorslaer 2013</xref>). This implies that the <italic>Alphapapillomaviruses</italic> (alpha-PVs) diverged well before the emergence of Homo sapiens in the primate lineage (<xref rid="R61" ref-type="bibr">Van Doorslaer 2013</xref>; <xref rid="R47" ref-type="bibr">Ong et al., 1993</xref>). However, little is known about factors contributing to PV evolution (<xref rid="R61" ref-type="bibr">Van Doorslaer 2013</xref>; <xref rid="R5" ref-type="bibr">Bravo and Félez-Sánchez 2015</xref>).</p><p id="P4">The apolipoprotein B mRNA-editing, enzyme-catalytic, polypeptide-like 3 (APOBEC3) family of proteins have recently been discovered as host restriction factors against various viruses and endogenous retroelements (<xref rid="R54" ref-type="bibr">Sheehy et al. 2002</xref>; <xref rid="R70" ref-type="bibr">Zheng et al. 2004</xref>; <xref rid="R4" ref-type="bibr">Bogerd et al. 2006</xref>; <xref rid="R26" ref-type="bibr">Hultquist et al. 2011</xref>; <xref rid="R56" ref-type="bibr">Suspène et al. 2011</xref>; <xref rid="R48" ref-type="bibr">Ooms et al. 2012</xref>; <xref rid="R35" ref-type="bibr">Liang et al. 2013</xref>; <xref rid="R40" ref-type="bibr">Minkah et al. 2014</xref>). APOBEC3 proteins deaminate cytosine (C) residues to uracil (U) in single-stranded DNA (<xref rid="R67" ref-type="bibr">Yu et al. 2004</xref>; <xref rid="R12" ref-type="bibr">Chiu et al. 2005</xref>). These mutations can lead to defects in gene transcription, genome replication and chromosomal integration thereby restricting the viral lifecycle. Of the multiple potential targets of human APOBEC3 (hA3), retroviruses, endogenous retroelements, DNA viruses and surprisingly, host genomic DNA, have been identified as common substrates (<xref rid="R22" ref-type="bibr">Harris and Liddament 2004</xref>; <xref rid="R58" ref-type="bibr">Turelli et al. 2004</xref>; <xref rid="R18" ref-type="bibr">Esnault et al. 2005</xref>; <xref rid="R4" ref-type="bibr">Bogerd et al. 2006</xref>; <xref rid="R9" ref-type="bibr">Chen et al. 2006</xref>; <xref rid="R8" ref-type="bibr">Burns et al. 2013</xref>; <xref rid="R8" ref-type="bibr">Burns, Temiz, and Harris 2013</xref>).</p><p id="P5">The human and non-human primate APOBEC3 family consists of seven members (A3A, A3B, A3C, A3D, A3F, A3G, and A3H), while other mammals encode fewer isoforms (<xref rid="R27" ref-type="bibr">Jarmuz et al. 2002</xref>; <xref rid="R42" ref-type="bibr">Münk et al. 2008</xref>; <xref rid="R34" ref-type="bibr">LaRue et al. 2009</xref>). Substrate specificity varies amongst the seven hA3 family members with each having slightly different target sequences based on the dinucleotide context of the cytosine residue (<xref rid="R70" ref-type="bibr">Zheng et al. 2004</xref>; <xref rid="R26" ref-type="bibr">Hultquist et al. 2011</xref>; <xref rid="R51" ref-type="bibr">Pham et al. 2003</xref>; <xref rid="R36" ref-type="bibr">Liddament et al. 2004</xref>; <xref rid="R67" ref-type="bibr">Yu et al. 2004</xref>; <xref rid="R28" ref-type="bibr">Jern et al. 2009</xref>; <xref rid="R44" ref-type="bibr">Nik-Zainal et al. 2012</xref>). In addition to differences in dinucleotide preferences, hA3s differ in their subcellular distribution during interphase and mitosis (<xref rid="R33" ref-type="bibr">Lackeyet al. 2013</xref>). These differences in localization impact the class of virus targeted by a particular hA3. For example, the predominantly cytoplasmic hA3F and hA3G are important inhibitors of productive replication of viruses whose replication cycle proceeds through obligatory reverse transcribed cDNA intermediates (<xref rid="R58" ref-type="bibr">Turelli et al. 2004</xref>; <xref rid="R12" ref-type="bibr">Chiu et al. 2005</xref>; <xref rid="R14" ref-type="bibr">Cullen 2006</xref>). In contrast, genetic editing of DNA viruses by nuclear hA3s (hA3A, hA3B, and hA3C) has been reported (<xref rid="R9" ref-type="bibr">Chen et al. 2006</xref>; <xref rid="R62" ref-type="bibr">Vartanian et al. 2008</xref>; <xref rid="R56" ref-type="bibr">Suspène et al. 2011</xref>; <xref rid="R33" ref-type="bibr">Lackeyet al. 2013</xref>; <xref rid="R63" ref-type="bibr">Wang et al. 2014</xref>).</p><p id="P6">It has been demonstrated that episomal HPV genomes are edited by hA3A and hA3B in cervical lesions and keratinocytes (<xref rid="R62" ref-type="bibr">Vartanian et al. 2008</xref>; <xref rid="R63" ref-type="bibr">Wang et al. 2014</xref>). Recent work by us and others revealed that HPV induced hA3A expression restricts HPV16 infectivity in human keratinocytes (<xref rid="R1" ref-type="bibr">Ahasan et al. 2015</xref>; <xref rid="R64" ref-type="bibr">Warren et al. 2015</xref>). If deamination by hA3s represents a bona fide host restriction mechanism, PVs may have evolved to limit APOBEC3 recognition motifs in their genomes through mutation and selection. Thus, we hypothesized that the host restriction function of APOBEC3 proteins contributed to PV genome evolution by selecting for depletion of specific dinucleotide sequences within PV genomes. To test this hypothesis, the relative abundance of dinucleotide sequence patterns in the full genome sequences of all known PVs was analyzed.</p><p id="P7">In this study, we investigated the link between APOBEC3 expression levels, tissue tropism of various HPV genotypes, and the extent of dinucleotide depletion in PV genomes. We propose a novel model suggesting that PV types replicating in tissues with high APOBEC3 levels have evolved to evade host restriction by selecting for variants containing reduced APOBEC3 target sites in their genomes. Our results provide important insight into the role of host restriction factors during PV evolution.</p></sec><sec id="S2"><title>2 Materials and methods</title><sec id="S3"><title>2.1 Papillomavirus genomes</title><p id="P8">We retrieved full-length reference sequences of all 274 PV types archived in the PV database (PaVE; pave.niaid.nih.gov, accessed on 1 August 2014).</p></sec><sec id="S4"><title>2.2 Calculation of dinucleotide frequencies</title><p id="P9">Dinucleotide frequencies were determined using three complementary approaches. First, in order to determine a single observed vs. expected (O/E) dinucleotide ratio across the entire viral genome, a custom python script was used (available from a public computer code repository [<ext-link ext-link-type="uri" xlink:href="http://www.github.com/Van-Doorslaer/Warren-et-al">http://www.github.com/Van-Doorslaer/Warren-et-al</ext-link>.]). The script uses the CompSeq program from Emboss. The expected frequencies of dinucleotide ‘words’ were estimated based on the observed frequency of single bases in the sequences. Alternatively, to assess local differences in dinucleotide content across the length of the viral genome, the frequency of each motif was counted within a 1 kb sliding window with a 100 bp overlap (<xref rid="R40" ref-type="bibr">Minkah et al. 2014</xref>). In order to generate a null model, the word occurrence was enumerated in 1,000 randomly shuffled versions of each 1 kb window. If the actual frequency fell within the lowest or highest fifth percentile, the word was considered to be under- or over-represented, respectively. The data are presented as percent 1 kb windows that fall in each category. Finally, the above-mentioned CompSeq based approach does not control for potential effects of coding restraints. To analyze this, the E1, E2, L1, and L2 ORFs were concatenated and the occurrence of each dinucleotide was counted in relationship to the position in a triplet (for coding sequences this triplet equals a codon). As a control, the non-coding upstream regulatory region (URR) was analyzed. In order to obtain the expected distribution, the sequences were randomly shuffled and dinucleotides were counted as earlier. This process was repeated 1,000 times.</p></sec><sec id="S5"><title>2.3 Tissue-specific gene expression analysis</title><p id="P10">Correlation of tissue-specific APOBEC3 gene expression analysis was performed using publicly available RNAseq data from the Genotype-Tissue Expression (GTEx) project (<ext-link ext-link-type="uri" xlink:href="http://gtexportal.org">gtexportal.org</ext-link>) (<xref rid="R37" ref-type="bibr">Lonsdale et al. 2013</xref>). DESeq2 (v 1.6.3) in R (v 3.1.0) was used for analysis of differentially expressed genes using DESeq2’s negative binomial distribution analysis model (<xref rid="R38" ref-type="bibr">Love, Huber, and Anders 2014</xref>).</p></sec><sec id="S6"><title>2.4 Statistical analysis</title><p id="P11">Student’s t test and one- or two-way analysis of variance (ANOVA) were used where appropriate. Data are presented as box-and-whisker plots with Tukey’s method for outliers noted as distinct data points. All graphs were generated using Prism software. Results were considered statistically significant at a <italic>P</italic>-value of <0.05.</p></sec><sec id="S7"><title>2.5 Phylogenetic analysis</title><p id="P12">The E1, E2, L1, and L2 sequences from 274 PVs were downloaded from PaVE. The nucleotide sequences were translated into amino acids and aligned using the MAFFT algorithm (<xref rid="R31" ref-type="bibr">Katoh and Standley 2013</xref>). The resulting alignment was back translated to nucleotides. The individual gene alignments were concatenated into a single supermatrix. A maximum likelihood tree was constructed using RaxML (<xref rid="R55" ref-type="bibr">Stamatakis 2006</xref>) implementing a gamma model allowing for among-site rate variation and variable substitution rates GTR<bold>+</bold>I<bold>+</bold>G; model selected using jModeltest (<xref rid="R15" ref-type="bibr">Darriba et al. 2012</xref>). The phylogenetic tree shown in <xref rid="F3" ref-type="fig">Figure 3</xref> was annotated using EvolView (<xref rid="R69" ref-type="bibr">Zhang et al. 2012</xref>). All computational analyses were performed on the CIPRES Science Gateway (<xref rid="R39" ref-type="bibr">Miller, Pfeiffer, and Schwartz 2010</xref>).</p></sec><sec id="S8"><title>2.6 Ancestral state reconstruction</title><p id="P13">The E1, E2, L1, and L2 nucleotide sequences belonging to the alpha, omega, upsilon, omicron, and dyodelta-PV genera were downloaded from PaVE. The nucleotide sequences were translated into amino acids and aligned using the MAFFT algorithm (<xref rid="R31" ref-type="bibr">Katoh and Standley 2013</xref>). The resulting alignment was reverse translated to nucleotides. Aligned sequences were concatenated and the resulting supermatrix was used to estimate a Bayesian phylogenetic tree using MrBayes (version 3; Cipres science gateway) implementing a GTR<bold>+</bold>I<bold>+</bold>G model of evolution. The MCMC-chain length was set to 1 × 10<sup>7</sup> while the chain was sampled every 1,000th generation. Chain convergence was analyzed using AWTY (<xref rid="R46" ref-type="bibr">Nylander et al. 2008</xref>). Bovine papillomavirus type 1 (BPV1), a member of the <italic>Deltapapillomavirs</italic> genus, was used to root the phylogenetic tree. <xref rid="F5" ref-type="fig">Figure 5</xref> shows the majority rule tree. Ancestral states were estimated using a Bayesian approach within BayesTraits (<xref rid="R2" ref-type="bibr">Barker, Meade, and Pagel 2007</xref>). The applied General Least squares approach removes the need for phylogenetic correction when comparing continuous traits between evolutionarily related species. Furthermore, in order to account for phylogenetic uncertainty during the tree building process, 500 trees were randomly sampled from the post-burn-in posterior sample of MrBayes trees and used for further analysis. Within BayesTraits, the MCMC analysis was run for 1 × 10<sup>9</sup> iterations. After a burn-in of 1 × 10<sup>6</sup> iterations, the chains were sampled every 10,000th iteration. The alpha (estimated root of the tree) and beta (directional change parameter of the directional model) parameters were estimated based on a uniform prior distribution. BayesTraits keeps a running tally of the logarithm of the harmonic mean of the likelihoods. When the analysis is completed, the final harmonic means for both models are compared using Bayes Factors (logBF = 2 (log[harmonic mean(complex model)]−log[harmonic mean(simple model)]). Kappa, lambda, and delta scaling parameters were estimated under both the ‘random walk’ and ‘directional walk’ models (<xref rid="R50" ref-type="bibr">Pagel 1999</xref>). These parameters test the tempo, mode, and phylogenetic associations of trait evolution, respectively. Following selection of appropriate model parameters, these parameters were used to estimate ancestral TC O/E ratios at indicated nodes. To confirm the robustness of the results, the analysis was performed three independent times.</p></sec></sec><sec id="S9"><title>3 Results</title><sec id="S10"><title>3.1 TC dinucleotides are significantly depleted in PV genomes</title><p id="P14">To determine whether there are preferred dinucleotide sequence patterns in PV genomes, we analyzed the relative abundance of all possible dinucleotide combinations in the full-length genomes of all 274 PVs deposited in the Papillomavirus Episteme database (<xref rid="R61" ref-type="bibr">Van Doorslaer et al. 2013</xref>) (PaVE; <ext-link ext-link-type="uri" xlink:href="http://pave.niaid.nih.gov">pave.niaid.nih.gov</ext-link>; accessed on 1 August 2014). The ratio of observed vs. expected (O/E) counts of all dinucleotides were calculated and this O/E ratio was used to determine whether dinucleotides were over-represented (ratio > 1), neutral (ratio = 1), or underrepresented (ratio < 1) in these genome sequences. Among all dinucleotide combinations, CG and TC dinucleotides are the most underrepresented dinucleotides in PV genomes (<xref rid="F1" ref-type="fig">Fig. 1A</xref> and <xref ref-type="supplementary-material" rid="SD2">Supplementary Table S1</xref>; depletion of the CG or TC dinucleotide compared with any other dinucleotide, <italic>P</italic><0.0001 by two-way ANOVA with Tukey correction). CG depletion has been previously noted for several small DNA viruses (<xref rid="R30" ref-type="bibr">Karlin, Doerfler, and Cardon 1994</xref>; <xref rid="R25" ref-type="bibr">Hoelzer, Shackelton, and Parrish 2008</xref>; <xref rid="R59" ref-type="bibr">Upadhyay et al. 2013</xref>), including PVs (<xref rid="R53" ref-type="bibr">Shackelton, Parrish, and Holmes 2006</xref>). Given that T<underline>C</underline> dinucleotides are preferred or specific targets of APOBEC3 DNA editing, these results imply that the PV genomes may be greatly affected by the host restriction functions of APOBEC3.</p><p id="P15">Each APOBEC3 family member preferably targets different dinucleotide sequences for cytidine deamination. For example, the T<underline>C</underline> dinucleotide is preferably targeted by A3A, A3B, A3C, A3F, and A3H (<xref rid="R36" ref-type="bibr">Liddament et al. 2004</xref>; <xref rid="R67" ref-type="bibr">Yu et al. 2004</xref>; <xref rid="R70" ref-type="bibr">Zheng et al. 2004</xref>; <xref rid="R26" ref-type="bibr">Hultquist et al. 2011</xref>; <xref rid="R44" ref-type="bibr">Nik-Zainal et al. 2012</xref>), while C<underline>C</underline> and R<underline>C</underline> (where R represents any pyrimidine) dinucleotides are preferably deaminated by A3G (<xref rid="R28" ref-type="bibr">Jern et al. 2009</xref>) and AID (<xref rid="R51" ref-type="bibr">Pham et al. 2003</xref>), respectively. Interestingly, among cytidine-containing dinucleotides, only CG and TC residues are significantly depleted in all papillomavirus genomes, including HPVs (<xref rid="F1" ref-type="fig">Fig. 1A–C</xref> and <xref ref-type="supplementary-material" rid="SD1">Supplementary Fig. S1</xref>; <italic>P</italic><0.0001 by two-way ANOVA with Tukey correction). Among the APOBEC3 family members preferably editing TC dinucleotides, hA3A, hA3B, and hA3H are known to restrict DNA viruses including HPVs and herpesviruses (<xref rid="R62" ref-type="bibr">Vartanian et al. 2008</xref>; <xref rid="R56" ref-type="bibr">Suspène et al. 2011</xref>; <xref rid="R40" ref-type="bibr">Minkah et al. 2014</xref>; <xref rid="R64" ref-type="bibr">Warren et al. 2015</xref>). Therefore, these results suggest that TC depletion may be influenced by host restriction factors, A3A, A3B, and/or A3H.</p></sec><sec id="S11"><title>3.2 APOBEC3 family members are differentially expressed at mucosal and cutaneous sites</title><p id="P16">HPVs can roughly be categorized into two groups based on tissue tropism, those that infect cutaneous skin, which mainly include beta- and gamma-PVs, and those that infect the mucosa, which predominantly include alpha-PVs (<xref rid="R16" ref-type="bibr">de Villiers et al. 2004</xref>). Because TC depletion is significant amongst all HPV types, we sought to determine whether the expression levels of all human APOBEC3 family members differed in mucosal vs. cutaneous skin. RNA sequencing data from the GTEx project (<xref rid="R37" ref-type="bibr">Lonsdale et al. 2013</xref>) was used to compare the expression levels of each human APOBEC3 gene in cervix and vagina to cutaneous skin using DESeq2 (<xref rid="R38" ref-type="bibr">Love, Huber, and Anders 2014</xref>). Of note, no strong evidence of potential batch effects was observed from PCA plots and hierarchical clustering. High APOBEC3 expression was observed for all isoforms, except A3B, when mucosal tissues (cervix and vagina) were compared with cutaneous skin (<xref rid="F2" ref-type="fig">Fig. 2A and B</xref>).</p></sec><sec id="S12"><title>3.3 TC content is significantly lower in alpha-PVs, compared with beta- and gamma-PVs</title><p id="P17">Gene expression analysis revealed increased basal expression of APOBEC3 family members at mucosal sites when compared with cutaneous skin. Because increased basal APOBEC3 expression may pose an intrinsic barrier to viral infection, it is likely that PVs infecting these tissues would evolve to evade APOBEC3 restriction pressure. Thus, we hypothesized that the TC content in mucosotropic PVs would be significantly lower than PVs with a cutaneous tropism. The TC dinucleotide ratio in the genomes of human alpha- (<italic>n</italic> = 64), beta- (<italic>n</italic> = 44), and gamma- (<italic>n</italic> = 50) PVs was calculated as described earlier. Notably, TC dinucleotides are significantly depleted in the genomes of alpha-PVs (mean O/E ratio <0.6), compared with the TC contents of beta- and gamma-PVs (mean O/E ratio >0.8) (<xref rid="F3" ref-type="fig">Fig. 3A</xref>). In contrast, all three genera contain similar levels of CG dinucleotides (<xref rid="F3" ref-type="fig">Fig. 3B</xref>). The distinct pattern of TC depletion in alpha-PVs becomes apparent when mapped on a PV phylogenetic tree (<xref rid="F3" ref-type="fig">Fig. 3C</xref>). These results suggest that high basal expression of APOBEC3s at mucosal sites may have affected the evolutionary trajectory of the alpha-PV clade.</p><p id="P18">Although alpha-PVs are primarily mucosotropic, a subset is known to infect cutaneous skin (<xref rid="R16" ref-type="bibr">de Villiers et al. 2004</xref>). If higher APOBEC3 levels affect the genomic TC content, mucosotropic alpha-PVs should have decreased TC ratios when compared with cutaneotropic alpha-PVs. When HPVs were stratified by tissue tropism, significantly lower TC contents were observed in mucosal alpha-HPVs compared with cutaneous alpha-HPVs (<xref rid="F4" ref-type="fig">Fig. 4A</xref>; <xref ref-type="supplementary-material" rid="SD2">Supplementary Table S2</xref>). A similar analysis of primarily cutaneous gamma-HPVs and their nearest neighbors (pi, tau, and dyoxi-PVs) revealed no differences in TC O/E ratios suggesting a specific effect of mucosal tropism on TC depletion (<xref rid="F4" ref-type="fig">Fig. 4B</xref>). Interestingly, a few gamma-PVs (including HPV101, 103, and 108) that have been isolated from cervical lesions also show significant depletion of TC dinucleotides (<xref rid="F3" ref-type="fig">Fig. 3C</xref>). This further supports an effect of tropism on TC dinucleotide content (<xref rid="R10" ref-type="bibr">Chen et al. 2007</xref>; <xref rid="R45" ref-type="bibr">Nobre et al. 2009</xref>). This analysis was based on available clinical data (<xref rid="R16" ref-type="bibr">de Villiers et al. 2004</xref>) and the tropism of many PVs is not yet clear. Taken together, our results suggest that one or more APOBEC3 family members that are highly expressed in mucosal tissues may have driven genome evolution of alpha-PVs.</p></sec><sec id="S13"><title>3.4 TC depletion represents an apomorphy of alpha-PVs</title><p id="P19">Our results point towards a role for the APOBEC3 proteins in the evolution of the alpha-PV genus. However, it remains possible that low TC content was a characteristic of the most recent common ancestor (MRCA) of alpha-PVs. In this case, a founder effect may be responsible for the TC contents we observe in extant members, rather than specific APOBEC3 effects. To address this issue, thorough evolutionary trait analysis was performed to model the evolution of the TC O/E ratios and reconstruct the ancestral alpha-PV genomes.</p><p id="P20">The alpha-PV genus forms a monophyletic clade with the omega, upsilon, omicron, and dyodelta-PV genera (<xref rid="F3" ref-type="fig">Fig. 3C</xref>) (<xref rid="R61" ref-type="bibr">Van Doorslaer 2013</xref>). <xref rid="F5" ref-type="fig">Figure 5</xref> shows the evolutionary relationship within this monophyletic clade. The computer program Continuous, as implemented within BayesTraits, was used to model the evolution of the TC O/E ratio. This approach allows for the characterization of trait evolution while automatically correcting for the influence of phylogenetic signal (<xref rid="R49" ref-type="bibr">Pagel 1997</xref>, <xref rid="R50" ref-type="bibr">1999</xref>). These analyses indicate that the TC O/E ratio evolved according to a directional-random walk model (Bayes Factor >12). Furthermore, the Delta parameter was estimated to be significantly <1, indicating that adaptive radiation best explains the extant O/E ratios. This suggests that the initial TC depletion occurred rapidly, followed by slower rates of change among closely related species. Additionally, this approach allows for the estimation of ancestral states at important nodes of the phylogenetic tree (<xref rid="F5" ref-type="fig">Fig. 5A</xref>). The ancestral state reconstruction suggests that the MRCA of the alpha, omega, upsilon, omicron, and dyodelta-PV clades had a TC O/E ratio of 0.78 (<xref rid="F5" ref-type="fig">Fig. 5B</xref>). Likewise, the last common ancestor of the alpha-PV genus had a TC O/E ratio of 0.79. Therefore, the MRCA of the alpha-PVs is predicted to have a TC O/E ratio that is significantly higher than the extant members of this genus. Thus, the predicted rapid TC depletion must have occurred after the alpha-PVs began to diverge. Importantly, this provides strong evidence that TC depletion is unique to alpha-PVs and begins to illuminate the evolutionary path responsible for TC depletion.</p></sec><sec id="S14"><title>3.5 APOBEC3 mediated editing directly affects Alphapapillomavirus evolution</title><p id="P21">It appears that the TC O/E ratios in alpha-PV genomes are dramatically lower when compared with other viral genera. Next, we determined whether APOBEC3 actively edited viral genomes. APOBEC3 deamination of TC residues can result in C-to-T transitions or C-to-G transversions that, when fixed in the viral genome, would increase the overall TT or TG content, respectively (<xref rid="R67" ref-type="bibr">Yu et al. 2004</xref>; <xref rid="R11" ref-type="bibr">Chiu and Greene 2008</xref>; <xref rid="R21" ref-type="bibr">Harris 2013</xref>). Although many mutations will be deleterious to the virus (<xref rid="R66" ref-type="bibr">Yang et al. 2007</xref>; <xref rid="R20" ref-type="bibr">Friedman and Stivers 2010</xref>), some of these mutations will be neutral or may even increase viral fitness (<xref rid="R41" ref-type="bibr">Mulder, Harari, and Simon, 2008</xref>; <xref rid="R28" ref-type="bibr">Jern et al. 2009</xref>; <xref rid="R65" ref-type="bibr">Wood et al. 2009</xref>; <xref rid="R32" ref-type="bibr">Kim et al. 2014</xref>). If APOBEC3 driven deamination events became fixed in the viral population, this should be reflected in a relative gain of the O/E ratios of TT or TG dinucleotides. Alternatively, if mutation followed by selection was responsible, any mutation disrupting the TC recognition motif would be selected for. Although alpha-HPV genomes are not enriched for TT dinucleotides (O/E ratio ~1), TG dinucleotides (and CA on the opposite strand) are the most over-represented dinucleotides (<xref rid="F6" ref-type="fig">Fig. 6A</xref>; increase of TG or CA dinucleotide compared with increase of any other dinucleotide, <italic>P</italic> < 0.0001 by two-way ANOVA with Tukey correction). No significant differences were noted between TG and CA using this analysis (<italic>P</italic> = 0.87). Additionally, GA dinucleotides (TC counterpart on negative strand) are also depleted in alpha-PVs thus suggesting that APOBEC3 effects are observed on both strands. These data suggest that TC depletion might in part be driven by APOBEC3-mediated deamination.</p></sec><sec id="S15"><title>3.6 Effects of coding potential on TC depletion</title><p id="P22">Roughly 85% of the PV genome encodes for viral proteins. The previous analyses considered dinucleotide ratios for the whole viral genome, regardless of coding potential. If APOBEC3 mediated editing affected PV evolution, it is likely that TC content is correlated with the codon position of the cytidine (NNT <underline>C</underline>NN vs. T<underline>C</underline>N vs. NT<underline>C</underline> with N representing any nucleotide). Strikingly, it appears that TC reduction is most pronounced when the cytidine is in the third codon position (NT<underline>C</underline> O/E <0.5; <italic>P</italic> < 0.0001 by two-way ANOVA with Tukey correction) (<xref rid="F6" ref-type="fig">Fig. 6B</xref>). Importantly, this NT<underline>C</underline> bias is not seen in the non-coding URR (<xref rid="F6" ref-type="fig">Fig. 6C</xref>). Because cytidine is at the third position of these codons, substitution of cytidine in any NT<underline>C</underline> codon does not incur any significant change in the corresponding amino acid during protein synthesis (<xref rid="R3" ref-type="bibr">Betts and Russell 2003</xref>). These results suggest that APOBEC3 editing primarily resulted in silent and/or conserved mutations.</p></sec><sec id="S16"><title>3.6 TC dinucleotides are dramatically depleted in the early gene regions for alpha-PVs</title><p id="P23">To examine if there was any bias for TC dinucleotide reduction in particular regions within alpha-PV genomes, the TC O/E ratio was calculated using a sliding window approach (a 1,000 bp window with a 100 bp overlap between two windows). <xref rid="F7" ref-type="fig">Figure 7A</xref> compares genomic TC content across all alpha-and beta-PVs. The alpha-PV genomes have a lower TC ratio across the early gene regions of the viral genome (<xref rid="F7" ref-type="fig">Fig. 7B</xref>).</p></sec></sec><sec id="S17"><title>4 Discussion</title><p id="P24">PV infections have been described in most amniotic hosts. Interestingly, these infections appear to be tissue specific and show a highly restricted host range. Based on these observations, it has been suggested that PVs co-speciate with their hosts (<xref rid="R61" ref-type="bibr">Van Doorslaer 2013</xref>). To date, limited evidence exists supporting positive selection of PV genomes (<xref rid="R61" ref-type="bibr">Van Doorslaer 2013</xref>), suggesting that genetic drift is the primary driving force of PV evolution. However, most studies have investigated the role of evolution in shaping individual viral proteins, not taking into account possible evolutionary selection working at the primary nucleotide level. This study was undertaken to survey extant genomes of all known PVs and to investigate a potential role for APOBEC3-mediated restriction in the evolution of these viruses.</p><p id="P25">Viral genomes are frequently targeted by host restriction factors such as restriction nucleases, RNAi-mediated mechanisms, and DNA editing enzymes such as APOBEC3 and ADAR (<xref rid="R29" ref-type="bibr">Johnson 2013</xref>). In order to achieve productive infection, viruses need to evade these antiviral host mechanisms. Primate lentiviruses, such as human immunodeficiency virus type 1 (HIV-1), encode multiple viral proteins that antagonize host restriction mechanisms. For example, the HIV-1 accessory proteins Vif and Vpu counteract host restriction by hA3G and tetherin, respectively (<xref rid="R54" ref-type="bibr">Sheehy et al. 2002</xref>; <xref rid="R43" ref-type="bibr">Neil, Zang, and Bieniasz 2008</xref>). Likewise, host genomes evolve new restriction mechanisms to detect and eliminate virus infections. Well exemplified by the APOBEC3 genetic locus, host genome evolution alongside constant viral exposure may have selected for beneficial duplication events that increased the number of APOBEC3 genes. A recent study suggests that the APOBEC3 genes might have originated from one ancestral APOBEC3 gene (<xref rid="R42" ref-type="bibr">Münk et al. 2008</xref>). During mammalian evolution, APOBEC3 genes evolved by gene duplication, fusions, and losses, probably to adapt to various virus infections. It is likely that this ‘arms race’ has been indispensible in shaping both viral and host genomes (<xref rid="R17" ref-type="bibr">Duggal and Emerman 2012</xref>).</p><p id="P26">Although viral regulatory proteins may be purposed to antagonize host restriction factors, viruses may also develop more subtle ways to subvert host recognition and restriction mechanisms. We and others have previously shown that HPV infection is restricted by hA3A and that HPV genomes show signs of APOBEC3 hypermutation (<xref rid="R62" ref-type="bibr">Vartanian et al. 2008</xref>; <xref rid="R63" ref-type="bibr">Wang et al. 2014</xref>; <xref rid="R64" ref-type="bibr">Warren et al. 2015</xref>). Given that APOBEC3s would pose a significant barrier to HPV infection, we sought to investigate how PVs may have evolved to evade APOBEC3 restriction.</p><p id="P27">Here, our study provides evidence that PV genomes are significantly depleted in TC dinucleotides, the preferred target sites of several APOBEC3 proteins. Although all PV genomes have decreased TC content, this depletion is strongest in mucosotropic alpha-PVs (<xref rid="F3" ref-type="fig">Figs. 3A</xref> and <xref rid="F4" ref-type="fig">4A</xref>). Importantly, gene expression data clearly shows significantly higher expression levels of all APOBEC3 isoforms (except A3B) in mucosal skin of cervix and vagina, compared with cutaneous skin (<xref rid="F2" ref-type="fig">Fig. 2A and B</xref>). This suggests that increased APOBEC3 expression in these tissues may have affected the evolution of a subset of viral genomes.</p><p id="P28">APOBEC3-mediated TC deamination is correlated with TT and/or TG mutational signatures in tumor tissues (<xref rid="R44" ref-type="bibr">Nik-Zainal et al. 2012</xref>). These mutational signatures arise from distinct pathways, thus providing clues into the mechanism of deaminated cytosine processing. Although C-to-T transitions may result from error prone DNA replication across uracil lesions, C-to-T transversions largely result from errors incurred during base excision repair (BER) pathways (<xref rid="R44" ref-type="bibr">Nik-Zainal et al. 2012</xref>; <xref rid="R24" ref-type="bibr">Henderson, Chakravarthy, Fenton, 2014</xref>). TC depletion led to an increase in TG (and CA on the opposite strand), but not TT dinucleotides in alpha-PVs. This suggests that APOBEC3-mediated deamination, coupled with BER, may have partially contributed to alpha-PV genome evolution (<xref rid="F6" ref-type="fig">Fig. 6A</xref>).</p><p id="P29">Although alpha-PVs appear to have minimized the TC O/E ratio across their genome (<xref rid="F7" ref-type="fig">Fig. 7</xref>), further analysis of coding regions revealed that TC depletion is dependent on the position of cytidine within a codon (<xref rid="F6" ref-type="fig">Fig. 6B and C</xref>). This analysis suggested that while alpha-PV genomes are able to tolerate NT<underline>C</underline> mutations, there was no decrease in T<underline>C</underline>N or NNT <underline>C</underline>NN sites. This is likely due to the degenerate nature of the genetic code (<xref rid="R3" ref-type="bibr">Betts and Russell 2003</xref>). Interestingly, it has been demonstrated that HPV genomes have a pronounced codon usage bias when compared with host genes (<xref rid="R6" ref-type="bibr">Bravo and Muller 2005</xref>; <xref rid="R13" ref-type="bibr">Cladel, Bertotto, and Christensen, 2010</xref>; <xref rid="R19" ref-type="bibr">Félez-Sánchez et al. 2015</xref>). This study may begin to shed light unto why papillomavirus genomes have evolved to use sub-optimal codons. Importantly, the observation that the third codon position is preferentially edited, suggests that these genomes are expected to still be vulnerable to APOBEC3 restriction because mutations would lead to changes in the amino acid sequence.</p><p id="P30">In addition to TC depletion, we also observed a significant reduction in CG dinucleotides in all PVs surveyed. This finding is consistent with others that have reported a significant depletion in CG residues across many small (<30 kb) DNA virus genomes (<xref rid="R30" ref-type="bibr">Karlin, Doerfler, Cardon 1994</xref>; <xref rid="R25" ref-type="bibr">Hoelzer, Shackelton, and Parrish 2008</xref>; <xref rid="R59" ref-type="bibr">Upadhyay et al. 2013</xref>). Previous work suggests that 5-methylcytosine modifications can negatively influence gene expression and lead to mutation via spontaneous deamination, both of which would negatively impact the course of a viral infection (<xref rid="R30" ref-type="bibr">Karlin, Doerfler, Cardon 1994</xref>). Additionally, TLR9 recognition of CG motifs in viral genomes can trigger antiviral immune responses such as production of type I interferons. High-risk HPV16 has previously been shown to inhibit TLR9 expression as a mechanism to avoid innate immune recognition (<xref rid="R23" ref-type="bibr">Hasan et al. 2007</xref>). Collectively, depletion of CG dinucleotides represents an evolutionarily conserved mechanism to evade innate immune responses.</p><p id="P31">Although further studies are necessary, our results consistently support a model by which the mucosal niche was colonized by an ancestral virus. Higher APOBEC3 levels in cervical and vaginal tissues resulted in increased editing of viral genomes and selection for those viruses with reduced TC contents. With about 85% of the genome coding for viral proteins, TC depletion appears to have been limited by coding potential. This is predicted by the evolutionary model suggesting that after the TC O/E ratio was minimized, initial rapid evolution was followed by slow drift (<xref rid="F5" ref-type="fig">Fig. 5</xref>).</p><p id="P32">In conclusion, our data suggest that viral evasion strategies that select against TC dinucleotides may be an important driver of PV evolution. Furthermore, our findings reinforce that tissue tropism and host range across these, and possibly other DNA viruses, may similarly be affected by interactions with APOBEC3.</p></sec><sec sec-type="supplementary-material" id="S18"><title>Supplementary Material</title><supplementary-material content-type="local-data" id="SD1"><label>Supplementary Figure</label><media xlink:href="NIHMS778642-supplement-Supplementary_Figure.pdf" orientation="portrait" xlink:type="simple" id="d37e859" position="anchor"></media></supplementary-material><supplementary-material content-type="local-data" id="SD2"><label>Supplementary Tables</label><media xlink:href="NIHMS778642-supplement-Supplementary_Tables.xlsx" orientation="portrait" xlink:type="simple" id="d37e863" position="anchor"></media></supplementary-material></sec></body><back><ack id="S19"><p>This work was supported by National Institutes of Health (R01 AI091968) to D.P. and (R01 CA117907-07) to J.M.E., the Intramural Research Program of the National Institute of Allergy and Infectious Diseases at the National Institutes of Health to K.V.D., and the Cancer League of Colorado to D.P. We thank Katerina Kechris for consultation on statistical analyses. We also thank Paul Lambert, Paul Ahlquist, Mario Santiago, Matt Daugherty, Rob DeSalle, and members of our lab for useful comments and suggestions.</p></ack><fn-group><fn id="FN4"><p><bold>Supplementary data</bold></p><p><xref ref-type="supplementary-material" rid="SD1">Supplementary data</xref> is available at Virus Evolution online.</p></fn><fn fn-type="COI-statement" id="FN5"><p>Conflict of interest: None declared</p></fn></fn-group><ref-list><ref id="R1"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ahasan</surname><given-names>MM</given-names></name><etal></etal></person-group><year>2015</year><article-title>APOBEC3A and C Decrease Human Papillomavirus 16 Pseudovirion Infectivity</article-title><source>Biochemical and Biophysical Research Communications</source><volume>457</volume><fpage>295</fpage><lpage>9</lpage><pub-id pub-id-type="pmid">25576866</pub-id></element-citation></ref><ref id="R2"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Barker</surname><given-names>D</given-names></name><name><surname>Meade</surname><given-names>A</given-names></name><name><surname>Pagel</surname><given-names>M</given-names></name></person-group><year>2007</year><article-title>Constrained Models of Evolution Lead to Improved Prediction of Functional Linkage from Correlated Gain and Loss of Genes</article-title><source>Bioinformatics</source><volume>23</volume><fpage>14</fpage><lpage>20</lpage><pub-id pub-id-type="pmid">17090580</pub-id></element-citation></ref><ref id="R3"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Betts</surname><given-names>MJ</given-names></name><name><surname>Russell</surname><given-names>RB</given-names></name></person-group><year>2003</year><article-title>Amino Acid Properties and Consequences of Substitutions</article-title><source>Bioinformatics for Geneticists</source><volume>317</volume><fpage>289</fpage></element-citation></ref><ref id="R4"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bogerd</surname><given-names>HP</given-names></name><etal></etal></person-group><year>2006</year><article-title>APOBEC3A and APOBEC3B are Potent Inhibitors of LTR-Retrotransposon Function in Human Cells</article-title><source>Nucleic Acids Research</source><volume>34</volume><fpage>89</fpage><lpage>95</lpage><pub-id pub-id-type="pmid">16407327</pub-id></element-citation></ref><ref id="R5"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bravo</surname><given-names>IG</given-names></name><name><surname>Félez-Sánchez</surname><given-names>M</given-names></name></person-group><year>2015</year><article-title>Papillomaviruses: Viral Evolution, Cancer and Evolutionary Medicine</article-title><source>Evolution, Medicine and Public Health</source><volume>2015</volume><fpage>32</fpage><lpage>51</lpage></element-citation></ref><ref id="R6"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bravo</surname><given-names>IG</given-names></name><name><surname>Félez-Sánchez</surname><given-names>M</given-names></name><name><surname>Muller</surname><given-names>M</given-names></name></person-group><year>2005</year><article-title>Codon Usage in Papillomavirus Genes: Practical and Functional Aspects</article-title><source>Papillomavirus Report</source><volume>16</volume><fpage>63</fpage><lpage>72</lpage></element-citation></ref><ref id="R7"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Burns</surname><given-names>MB</given-names></name><etal></etal></person-group><year>2013</year><article-title>APOBEC3B is an Enzymatic Source of Mutation in Breast Cancer</article-title><source>Nature</source><volume>494</volume><fpage>366</fpage><lpage>70</lpage><pub-id pub-id-type="pmid">23389445</pub-id></element-citation></ref><ref id="R8"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Burns</surname><given-names>MB</given-names></name><name><surname>Temiz</surname><given-names>NA</given-names></name><name><surname>Harris</surname><given-names>RS</given-names></name></person-group><year>2013</year><article-title>Evidence for APOBEC3B Mutagenesis in Multiple Human Cancers</article-title><source>Nature Genetics</source><volume>45</volume><fpage>977</fpage><lpage>83</lpage><pub-id pub-id-type="pmid">23852168</pub-id></element-citation></ref><ref id="R9"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname><given-names>H</given-names></name><etal></etal></person-group><year>2006</year><article-title>APOBEC3A is a Potent Inhibitor of Adeno-Associated Virus and Retrotransposons</article-title><source>Current Biology</source><volume>16</volume><fpage>480</fpage><lpage>5</lpage><pub-id pub-id-type="pmid">16527742</pub-id></element-citation></ref><ref id="R10"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname><given-names>Z</given-names></name><etal></etal></person-group><year>2007</year><article-title>Human Papillomavirus (HPV) Types 101 and 103 Isolated from Cervicovaginal Cells Lack an E6 Open Reading Frame (ORF) and are Related to Gamma-Papillomaviruses</article-title><source>Virology</source><volume>360</volume><fpage>447</fpage><lpage>53</lpage><pub-id pub-id-type="pmid">17125811</pub-id></element-citation></ref><ref id="R11"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chiu</surname><given-names>Y-L</given-names></name><name><surname>Greene</surname><given-names>WC</given-names></name></person-group><year>2008</year><article-title>The APOBEC3 Cytidine Deaminases: An Innate Defensive Network Opposing Exogenous Retroviruses and Endogenous Retroelements</article-title><source>Annual Review of Immunology</source><volume>26</volume><fpage>317</fpage><lpage>53</lpage></element-citation></ref><ref id="R12"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chiu</surname><given-names>Y-L</given-names></name><name><surname>Greene</surname><given-names>WC</given-names></name><etal></etal></person-group><year>2005</year><article-title>Cellular APOBEC3G Restricts HIV-1 Infection in Resting CD4+ T Cells</article-title><source>Nature</source><volume>435</volume><fpage>108</fpage><lpage>14</lpage><pub-id pub-id-type="pmid">15829920</pub-id></element-citation></ref><ref id="R13"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cladel</surname><given-names>NM</given-names></name><name><surname>Bertotto</surname><given-names>A</given-names></name><name><surname>Christensen</surname><given-names>ND</given-names></name></person-group><year>2010</year><article-title>Human Alpha and Beta Papillomaviruses Use Different Synonymous Codon Profiles</article-title><source>Virus Genes</source><volume>40</volume><fpage>329</fpage><lpage>40</lpage><pub-id pub-id-type="pmid">20157772</pub-id></element-citation></ref><ref id="R14"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cullen</surname><given-names>BR</given-names></name></person-group><year>2006</year><article-title>Role and Mechanism of Action of the APOBEC3 Family of Antiretroviral Resistance Factors</article-title><source>Journal of Virology</source><volume>80</volume><fpage>1067</fpage><lpage>76</lpage><pub-id pub-id-type="pmid">16414984</pub-id></element-citation></ref><ref id="R15"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Darriba</surname><given-names>D</given-names></name><name><surname>Taboada</surname><given-names>GL</given-names></name><name><surname>Doallo</surname><given-names>R</given-names></name><name><surname>Posada</surname><given-names>D</given-names></name></person-group><year>2012</year><article-title>jModelTest 2: More Models, New Heuristics and Parallel Computing</article-title><source>Nature Methods</source><volume>9</volume><fpage>772</fpage><pub-id pub-id-type="pmid">22847109</pub-id></element-citation></ref><ref id="R16"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>de Villiers</surname><given-names>E-M</given-names></name><etal></etal></person-group><year>2004</year><article-title>Classification of Papillomaviruses</article-title><source>Virology</source><volume>324</volume><fpage>17</fpage><lpage>27</lpage><pub-id pub-id-type="pmid">15183049</pub-id></element-citation></ref><ref id="R17"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Duggal</surname><given-names>NK</given-names></name><name><surname>Emerman</surname><given-names>M</given-names></name></person-group><year>2012</year><article-title>Evolutionary Conflicts Between Viruses and Restriction Factors Shape Immunity</article-title><source>Nature Review of Immunology</source><volume>12</volume><fpage>687</fpage><lpage>95</lpage></element-citation></ref><ref id="R18"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Esnault</surname><given-names>CC</given-names></name><etal></etal></person-group><year>2005</year><article-title>APOBEC3G Cytidine Deaminase Inhibits Retrotransposition of Endogenous Retroviruses</article-title><source>Nature</source><volume>433</volume><fpage>430</fpage><lpage>3</lpage><pub-id pub-id-type="pmid">15674295</pub-id></element-citation></ref><ref id="R19"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Félez-Sánchez</surname><given-names>M</given-names></name><etal></etal></person-group><year>2015</year><article-title>Cancer, Warts, or Asymptomatic Infections: Clinical Presentation Matches Codon Usage Preferences in Human Papillomaviruses</article-title><source>Genome Biological Evolution</source><volume>7</volume><fpage>2117</fpage><lpage>35</lpage></element-citation></ref><ref id="R20"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Friedman</surname><given-names>JI</given-names></name><name><surname>Stivers</surname><given-names>JT</given-names></name></person-group><year>2010</year><article-title>Detection of Damaged DNA Bases by DNA Glycosylase Enzymes</article-title><source>Biochemistry</source><volume>49</volume><fpage>4957</fpage><lpage>67</lpage><pub-id pub-id-type="pmid">20469926</pub-id></element-citation></ref><ref id="R21"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Harris</surname><given-names>RS</given-names></name></person-group><year>2013</year><article-title>Cancer Mutation Signatures, DNA Damage Mechanisms, and Potential Clinical Implications</article-title><source>Genome Medicine</source><volume>5</volume><fpage>87</fpage><pub-id pub-id-type="pmid">24073723</pub-id></element-citation></ref><ref id="R22"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Harris</surname><given-names>RS</given-names></name><name><surname>Liddament</surname><given-names>MT</given-names></name></person-group><year>2004</year><article-title>Retroviral Restriction by APOBEC Proteins</article-title><source>Nature Reviews Immunology</source><volume>4</volume><fpage>868</fpage><lpage>77</lpage></element-citation></ref><ref id="R23"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hasan</surname><given-names>UA</given-names></name><etal></etal></person-group><year>2007</year><article-title>TLR9 Expression and Function is Abolished by the Cervical Cancer-Associated Human Papillomavirus Type 16</article-title><source>The Journal of Immunology</source><volume>178</volume><fpage>3186</fpage><lpage>97</lpage><pub-id pub-id-type="pmid">17312167</pub-id></element-citation></ref><ref id="R24"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Henderson</surname><given-names>S</given-names></name><name><surname>Chakravarthy</surname><given-names>A</given-names></name><name><surname>Fenton</surname><given-names>T</given-names></name></person-group><year>2014</year><article-title>When Defense Turns into Attack</article-title><source>Molecular and Cellular Oncology</source><volume>1</volume><fpage>e29914</fpage><pub-id pub-id-type="pmid">27308321</pub-id></element-citation></ref><ref id="R25"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hoelzer</surname><given-names>K</given-names></name><name><surname>Shackelton</surname><given-names>LA</given-names></name><name><surname>Parrish</surname><given-names>CR</given-names></name></person-group><year>2008</year><article-title>Presence and Role of Cytosine Methylation in DNA Viruses of Animals</article-title><source>Nucleic Acids Research</source><volume>36</volume><fpage>2825</fpage><lpage>37</lpage><pub-id pub-id-type="pmid">18367473</pub-id></element-citation></ref><ref id="R26"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hultquist</surname><given-names>JF</given-names></name><etal></etal></person-group><year>2011</year><article-title>Human and Rhesus APOBEC3D, APOBEC3F, APOBEC3G, and APOBEC3H Demonstrate a Conserved Capacity to Restrict Vif-Deficient HIV-1</article-title><source>Journal of Virology</source><volume>85</volume><fpage>11220</fpage><lpage>34</lpage><pub-id pub-id-type="pmid">21835787</pub-id></element-citation></ref><ref id="R27"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jarmuz</surname><given-names>A</given-names></name><etal></etal></person-group><year>2002</year><article-title>An Anthropoid-Specific Locus of Orphan C to U RNA-Editing Enzymes on Chromosome 22</article-title><source>Genomics</source><volume>79</volume><fpage>285</fpage><lpage>96</lpage><pub-id pub-id-type="pmid">11863358</pub-id></element-citation></ref><ref id="R28"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jern</surname><given-names>P</given-names></name><name><surname>Russell</surname><given-names>RA</given-names></name><name><surname>Pathak</surname><given-names>VK</given-names></name><name><surname>Coffin</surname><given-names>JM</given-names></name></person-group><year>2009</year><article-title>Likely Role of APOBEC3G-Mediated G-to-A Mutations in HIV-1 Evolution and Drug Resistance</article-title><source>PLoS Pathogens</source><volume>5</volume><fpage>e1000367</fpage><pub-id pub-id-type="pmid">19343218</pub-id></element-citation></ref><ref id="R29"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Johnson</surname><given-names>WE</given-names></name></person-group><year>2013</year><article-title>Rapid Adversarial Co-Evolution of Viruses and Cellular Restriction Factors</article-title><source>Current Topics in Microbiology and Immunology</source><volume>371</volume><fpage>123</fpage><lpage>51</lpage><pub-id pub-id-type="pmid">23686234</pub-id></element-citation></ref><ref id="R30"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Karlin</surname><given-names>S</given-names></name><name><surname>Doerfler</surname><given-names>W</given-names></name><name><surname>Cardon</surname><given-names>LR</given-names></name></person-group><year>1994</year><article-title>Why is CpG Suppressed in the Genomes of Virtually all Small Eukaryotic Viruses but not in those of Large Eukaryotic Viruses?</article-title><source>Journal of Virology</source><volume>68</volume><fpage>2889</fpage><lpage>97</lpage><pub-id pub-id-type="pmid">8151759</pub-id></element-citation></ref><ref id="R31"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Katoh</surname><given-names>K</given-names></name><name><surname>Standley</surname><given-names>DM</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><fpage>772</fpage><lpage>80</lpage><pub-id pub-id-type="pmid">23329690</pub-id></element-citation></ref><ref id="R32"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kim</surname><given-names>E-Y</given-names></name><etal></etal></person-group><year>2014</year><article-title>Human APOBEC3 Induced Mutation of Human Immunodeficiency Virus Type-1 Contributes to Adaptation and Evolution in Natural Infection</article-title><source>PLoS Pathogens</source><volume>10</volume><fpage>e1004281</fpage><pub-id pub-id-type="pmid">25080100</pub-id></element-citation></ref><ref id="R33"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lackey</surname><given-names>L</given-names></name><name><surname>Law</surname><given-names>EK</given-names></name><name><surname>Brown</surname><given-names>WL</given-names></name><name><surname>Harris</surname><given-names>RS</given-names></name></person-group><year>2013</year><article-title>Subcellular Localization of the APOBEC3 Proteins During Mitosis and Implications for Genomic DNA Deamination</article-title><source>Cell Cycle</source><volume>12</volume><fpage>762</fpage><lpage>72</lpage><pub-id pub-id-type="pmid">23388464</pub-id></element-citation></ref><ref id="R34"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>LaRue</surname><given-names>RS</given-names></name><etal></etal></person-group><year>2009</year><article-title>Guidelines for Naming Nonprimate APOBEC3 Genes and Proteins</article-title><source>Journal of Virology</source><volume>83</volume><fpage>494</fpage><lpage>7</lpage><pub-id pub-id-type="pmid">18987154</pub-id></element-citation></ref><ref id="R35"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Liang</surname><given-names>G</given-names></name><etal></etal></person-group><year>2013</year><article-title>RNA Editing of Hepatitis B Virus Transcripts by Activation-Induced Cytidine Deaminase</article-title><source>Proceedings of the National Academy of Sciences</source><volume>110</volume><fpage>2246</fpage><lpage>51</lpage></element-citation></ref><ref id="R36"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Liddament</surname><given-names>MT</given-names></name><name><surname>Brown</surname><given-names>WL</given-names></name><name><surname>Schumacher</surname><given-names>AJ</given-names></name><name><surname>Harris</surname><given-names>RS</given-names></name></person-group><year>2004</year><article-title>APOBEC3F Properties and Hypermutation Preferences Indicate Activity against HIV-1 in vivo</article-title><source>Current Biology</source><volume>14</volume><fpage>1385</fpage><lpage>91</lpage><pub-id pub-id-type="pmid">15296757</pub-id></element-citation></ref><ref id="R37"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lonsdale</surname><given-names>J</given-names></name><etal></etal></person-group><year>2013</year><article-title>The Genotype-Tissue Expression (GTEx) Project</article-title><source>Nature Genetics</source><volume>45</volume><fpage>580</fpage><lpage>5</lpage><pub-id pub-id-type="pmid">23715323</pub-id></element-citation></ref><ref id="R38"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Love</surname><given-names>MI</given-names></name><name><surname>Huber</surname><given-names>W</given-names></name><name><surname>Anders</surname><given-names>S</given-names></name></person-group><year>2014</year><article-title>Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data with DESeq2</article-title><source>Genome Biology</source><volume>15</volume><fpage>550</fpage><pub-id pub-id-type="pmid">25516281</pub-id></element-citation></ref><ref id="R39"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Miller</surname><given-names>M</given-names></name><name><surname>Pfeiffer</surname><given-names>W</given-names></name><name><surname>Schwartz</surname><given-names>T</given-names></name></person-group><year>2010</year><article-title>Creating the CIPRES Science Gateway for Inference of Large Phylogenetic Trees</article-title><source>Gateway Computing Environments Workshop (GCE)</source><fpage>1</fpage><lpage>8</lpage></element-citation></ref><ref id="R40"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Minkah</surname><given-names>N</given-names></name><etal></etal></person-group><year>2014</year><article-title>Host Restriction of Murine Gammaherpesvirus 68 Replication by Human APOBEC3 Cytidine Deaminases but not Murine APOBEC3</article-title><source>Virology</source><fpage>454</fpage><lpage>455</lpage><fpage>215</fpage><lpage>26</lpage><pub-id pub-id-type="pmid">25243334</pub-id></element-citation></ref><ref id="R41"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mulder</surname><given-names>LCF</given-names></name><name><surname>Harari</surname><given-names>A</given-names></name><name><surname>Simon</surname><given-names>V</given-names></name></person-group><year>2008</year><article-title>Cytidine Deamination Induced HIV-1 Drug Resistance</article-title><source>Proceedings of the National Academy of Sciences</source><volume>105</volume><fpage>5501</fpage><lpage>6</lpage></element-citation></ref><ref id="R42"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Münk</surname><given-names>C</given-names></name><etal></etal></person-group><year>2008</year><article-title>Functions, Structure, and Read-Through Alternative Splicing of Feline APOBEC3 Genes</article-title><source>Genome Biology</source><volume>9</volume><fpage>R48</fpage><pub-id pub-id-type="pmid">18315870</pub-id></element-citation></ref><ref id="R43"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Neil</surname><given-names>SJD</given-names></name><name><surname>Zang</surname><given-names>T</given-names></name><name><surname>Bieniasz</surname><given-names>PD</given-names></name></person-group><year>2008</year><article-title>Tetherin Inhibits Retrovirus Release and is Antagonized by HIV-1 Vpu</article-title><source>Nature</source><volume>451</volume><fpage>425</fpage><lpage>30</lpage><pub-id pub-id-type="pmid">18200009</pub-id></element-citation></ref><ref id="R44"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nik-Zainal</surname><given-names>S</given-names></name><etal></etal><collab>Breast Cancer Working Group of the International Cancer Genome Consortium</collab></person-group><year>2012</year><article-title>Mutational Processes Molding the Genomes of 21 Breast Cancers</article-title><source>Cell</source><volume>149</volume><fpage>979</fpage><lpage>93</lpage><pub-id pub-id-type="pmid">22608084</pub-id></element-citation></ref><ref id="R45"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nobre</surname><given-names>RJ</given-names></name><etal></etal></person-group><year>2009</year><article-title>E7 Oncoprotein of Novel Human Papillomavirus Type 108 Lacking the E6 Gene Induces Dysplasia in Organotypic Keratinocyte Cultures</article-title><source>Journal of Virology</source><volume>83</volume><fpage>2907</fpage><lpage>16</lpage><pub-id pub-id-type="pmid">19153227</pub-id></element-citation></ref><ref id="R46"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nylander</surname><given-names>JAA</given-names></name><name><surname>Wilgenbusch</surname><given-names>JC</given-names></name><name><surname>Warren</surname><given-names>DL</given-names></name><name><surname>Swofford</surname><given-names>DL</given-names></name></person-group><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><fpage>581</fpage><lpage>3</lpage><pub-id pub-id-type="pmid">17766271</pub-id></element-citation></ref><ref id="R47"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ong</surname><given-names>CK</given-names></name><etal></etal></person-group><year>1993</year><article-title>Evolution of Human Papillomavirus Type 18: An Ancient Phylogenetic Root in Africa and Intratype Diversity Reflect Coevolution with Human Ethnic Groups</article-title><source>Journal of Virology</source><volume>67</volume><fpage>6424</fpage><lpage>31</lpage><pub-id pub-id-type="pmid">8411344</pub-id></element-citation></ref><ref id="R48"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ooms</surname><given-names>M</given-names></name><etal></etal></person-group><year>2012</year><article-title>APOBEC3A, APOBEC3B, and APOBEC3H Haplotype 2 Restrict Human T-Lymphotropic Virus Type 1</article-title><source>Journal of Virology</source><volume>86</volume><fpage>6097</fpage><lpage>108</lpage><pub-id pub-id-type="pmid">22457529</pub-id></element-citation></ref><ref id="R49"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pagel</surname><given-names>M</given-names></name></person-group><year>1997</year><article-title>Inferring Evolutionary Processes from Phylogenies</article-title><source>Zoologica Scripta</source><volume>26</volume><fpage>331</fpage><lpage>48</lpage></element-citation></ref><ref id="R50"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pagel</surname><given-names>M</given-names></name></person-group><year>1999</year><article-title>Inferring the Historical Patterns of Biological Evolution</article-title><source>Nature</source><volume>401</volume><fpage>877</fpage><lpage>84</lpage><pub-id pub-id-type="pmid">10553904</pub-id></element-citation></ref><ref id="R51"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pham</surname><given-names>P</given-names></name><name><surname>Bransteitter</surname><given-names>R</given-names></name><name><surname>Petruska</surname><given-names>J</given-names></name><name><surname>Goodman</surname><given-names>MF</given-names></name></person-group><year>2003</year><article-title>Processive AID-Catalysed Cytosine Deamination on Single-Stranded DNA Simulates Somatic Hypermutation</article-title><source>Nature</source><volume>424</volume><fpage>103</fpage><lpage>7</lpage><pub-id pub-id-type="pmid">12819663</pub-id></element-citation></ref><ref id="R52"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rector</surname><given-names>A</given-names></name><etal></etal></person-group><year>2007</year><article-title>Ancient Papillomavirus-Host Co-Speciation in Felidae</article-title><source>Genome Biology</source><volume>8</volume><fpage>R57</fpage><pub-id pub-id-type="pmid">17430578</pub-id></element-citation></ref><ref id="R53"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Shackelton</surname><given-names>LA</given-names></name><name><surname>Parrish</surname><given-names>CR</given-names></name><name><surname>Holmes</surname><given-names>EC</given-names></name></person-group><year>2006</year><article-title>Evolutionary Basis of Codon Usage and Nucleotide Composition Bias in Vertebrate DNA Viruses</article-title><source>Journal of Molecular Evolution</source><volume>62</volume><fpage>551</fpage><lpage>63</lpage><pub-id pub-id-type="pmid">16557338</pub-id></element-citation></ref><ref id="R54"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sheehy</surname><given-names>AM</given-names></name><name><surname>Gaddis</surname><given-names>NC</given-names></name><name><surname>Choi</surname><given-names>JD</given-names></name><name><surname>Malim</surname><given-names>MH</given-names></name></person-group><year>2002</year><article-title>Isolation of a Human Gene that Inhibits HIV-1 Infection and is Suppressed by the Viral Vif Protein</article-title><source>Nature</source><volume>418</volume><fpage>646</fpage><lpage>50</lpage><pub-id pub-id-type="pmid">12167863</pub-id></element-citation></ref><ref id="R55"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Stamatakis</surname><given-names>A</given-names></name></person-group><year>2006</year><article-title>RAxML-VI-HPC: Maximum Likelihood-Based Phylogenetic Analyses with Thousands of Taxa and Mixed Models</article-title><source>Bioinformatics</source><volume>22</volume><fpage>2688</fpage><lpage>90</lpage><pub-id pub-id-type="pmid">16928733</pub-id></element-citation></ref><ref id="R56"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Suspène</surname><given-names>R</given-names></name><etal></etal></person-group><year>2011</year><article-title>Genetic Editing of Herpes Simplex Virus 1 and Epstein-Barr Herpesvirus Genomes by Human APOBEC3 Cytidine Deaminases in Culture and in vivo</article-title><source>Journal of Virology</source><volume>85</volume><fpage>7594</fpage><lpage>602</lpage><pub-id pub-id-type="pmid">21632763</pub-id></element-citation></ref><ref id="R57"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Suspène</surname><given-names>R</given-names></name><etal></etal></person-group><year>2011</year><article-title>Genetic Editing of Herpes Simplex Virus 1 and Epstein-Barr Herpesvirus Genomes by Human APOBEC3 Cytidine Deaminases in Culture and in vivo</article-title><source>Journal of Virology</source><volume>85</volume><fpage>7594</fpage><lpage>602</lpage><pub-id pub-id-type="pmid">21632763</pub-id></element-citation></ref><ref id="R58"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Turelli</surname><given-names>P</given-names></name><etal></etal></person-group><year>2004</year><article-title>Inhibition of Hepatitis B Virus Replication by APOBEC3G</article-title><source>Science</source><volume>303</volume><fpage>1829</fpage><pub-id pub-id-type="pmid">15031497</pub-id></element-citation></ref><ref id="R59"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Upadhyay</surname><given-names>M</given-names></name><etal></etal></person-group><year>2013</year><article-title>CpG Dinucleotide Frequencies Reveal the Role of Host Methylation Capabilities in Parvovirus Evolution</article-title><source>Journal of Virology</source><volume>87</volume><fpage>13816</fpage><lpage>24</lpage><pub-id pub-id-type="pmid">24109231</pub-id></element-citation></ref><ref id="R60"><element-citation publication-type="journal"><person-group person-group-type="author"><name><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><fpage>11</fpage><lpage>20</lpage><pub-id pub-id-type="pmid">23769415</pub-id></element-citation></ref><ref id="R61"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Van Doorslaer</surname><given-names>K</given-names></name><etal></etal></person-group><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</volume><fpage>D571</fpage><lpage>8</lpage><pub-id pub-id-type="pmid">23093593</pub-id></element-citation></ref><ref id="R62"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Vartanian</surname><given-names>J-P</given-names></name><name><surname>Guétard</surname><given-names>D</given-names></name><name><surname>Henry</surname><given-names>M</given-names></name><name><surname>Wain-Hobson</surname><given-names>S</given-names></name></person-group><year>2008</year><article-title>Evidence for Editing of Human Papillomavirus DNA by APOBEC3 in Benign and Precancerous Lesions</article-title><source>Science</source><volume>320</volume><fpage>230</fpage><lpage>3</lpage><pub-id pub-id-type="pmid">18403710</pub-id></element-citation></ref><ref id="R63"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname><given-names>Z</given-names></name><etal></etal></person-group><year>2014</year><article-title>APOBEC3 Deaminases Induce Hypermutation in Human Papillomavirus 16 DNA Upon Beta Interferon Stimulation</article-title><source>Journal of Virology</source><volume>88</volume><fpage>1308</fpage><lpage>17</lpage><pub-id pub-id-type="pmid">24227842</pub-id></element-citation></ref><ref id="R64"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Warren</surname><given-names>CJ</given-names></name><etal></etal></person-group><year>2015</year><article-title>APOBEC3A Functions as a Restriction Factor of Human Papillomavirus</article-title><source>Journal of Virology</source><volume>89</volume><fpage>688</fpage><lpage>702</lpage><pub-id pub-id-type="pmid">25355878</pub-id></element-citation></ref><ref id="R65"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wood</surname><given-names>N</given-names></name><etal></etal></person-group><year>2009</year><article-title>HIV Evolution in Early Infection: Selection Pressures, Patterns of Insertion and Deletion, and the Impact of APOBEC</article-title><source>PLoS Pathogens</source><volume>5</volume><fpage>e1000414</fpage><pub-id pub-id-type="pmid">19424423</pub-id></element-citation></ref><ref id="R66"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yang</surname><given-names>B</given-names></name><etal></etal></person-group><year>2007</year><article-title>Virion-Associated Uracil DNA Glycosylase-2 and Apurinic/Apyrimidinic Endonuclease are Involved in the Degradation of APOBEC3G-Edited Nascent HIV-1 DNA</article-title><source>Journal of Biological Chemistry</source><volume>282</volume><fpage>11667</fpage><lpage>75</lpage><pub-id pub-id-type="pmid">17272283</pub-id></element-citation></ref><ref id="R67"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yu</surname><given-names>Q</given-names></name><etal></etal></person-group><year>2004</year><article-title>APOBEC3B and APOBEC3C are Potent Inhibitors of Simian Immunodeficiency Virus Replication</article-title><source>Journal of Biological Chemistry</source><volume>279</volume><fpage>53379</fpage><lpage>86</lpage><pub-id pub-id-type="pmid">15466872</pub-id></element-citation></ref><ref id="R68"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yu</surname><given-names>Q</given-names></name><etal></etal></person-group><year>2004</year><article-title>Single-Strand Specificity of APOBEC3G Accounts for Minus-Strand Deamination of the HIV Genome</article-title><source>Nature Structural and Molecular Biology</source><volume>11</volume><fpage>435</fpage><lpage>42</lpage></element-citation></ref><ref id="R69"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname><given-names>H</given-names></name><etal></etal></person-group><year>2012</year><article-title>EvolView, An Online Tool for Visualizing, Annotating and Managing Phylogenetic Trees</article-title><source>Nucleic Acids Research</source><volume>40</volume><fpage>W569</fpage><lpage>72</lpage><pub-id pub-id-type="pmid">22695796</pub-id></element-citation></ref><ref id="R70"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zheng</surname><given-names>Y-H</given-names></name><etal></etal></person-group><year>2004</year><article-title>Human APOBEC3F is Another Host Factor that Blocks Human Immunodeficiency Virus Type 1 Replication</article-title><source>Journal of Virology</source><volume>78</volume><fpage>6073</fpage><lpage>6</lpage><pub-id pub-id-type="pmid">15141007</pub-id></element-citation></ref></ref-list></back><floats-group><fig id="F1" orientation="portrait" position="float"><label>Figure 1</label><caption><p>TC and CG dinucleotide sequences are significantly depleted in PV genomes. (A) The observed vs. expected (O/E) ratios of each dinucleotide in PV genome sequences from all PV genotypes (<italic>n</italic> = 274) were calculated using a custom wrapper around the CompSeq program from the EMBOSS software suite. (B) The O/E ratios of TC, CC, AC, and GC dinucleotides, preferred editing sites of A3A, A3G, and AID, respectively, were analyzed with the genome sequences of all PVs (A; <italic>n</italic> = 274) and HPVs (B; <italic>n</italic> = 161). Box-and-whisker plots are shown with the outliers (black triangles) identified using Tukey’s method. CG and TC dinucleotide are the most depleted cytidine containing dinucleotides in all PVs and HPVs (depletion of CG or TC compared with depletion of any other dinucleotide, <italic>P</italic><0.0001 by two-way ANOVA with Tukey correction).</p></caption><graphic xlink:href="nihms778642f1"></graphic></fig><fig id="F2" orientation="portrait" position="float"><label>Figure 2</label><caption><p>The expression levels of APOBEC3s are different between cervicovaginal and cutaneous skin. (A) Gene expression levels of APOBEC3s were analyzed using RNA sequencing data retrieved from the Genotype-Tissue Expression project database (<ext-link ext-link-type="uri" xlink:href="http://gtexportal.org">gtexportal.org</ext-link>) (<xref rid="R37" ref-type="bibr">Lonsdale et al. 2013</xref>) and analyzed using DESeq2. Data are presented as fold change in cervix (<italic>n</italic> = 6) or vagina (<italic>n</italic> = 34) compared with cutaneous skin (<italic>n</italic> = 41). (B) Table showing the fold change and <italic>P</italic>-values adjusted for repeated measurements.</p></caption><graphic xlink:href="nihms778642f2"></graphic></fig><fig id="F3" orientation="portrait" position="float"><label>Figure 3</label><caption><p>TC content is significantly lower in alpha-PVs, compared with beta- and gamma-PVs. The O/E ratios of TC (A) and CG (B) dinucleotides were compared between alpha- (<italic>n</italic> = 64), beta- (<italic>n</italic> = 44), and gamma- (<italic>n</italic> = 50) HPVs. Box-and-whisker plots are shown with the outliers (black triangles) identified using Tukey’s method. <italic>P</italic>-values were determined by one-way ANOVA analysis using Prism software. (C) A maximum likelihood phylogenetic tree is shown comparing the O/E ratios of TC dinucleotides from all PV genomes (<italic>n</italic> = 274). The alpha-, beta-, and gamma-PV genera are indicated. Bars indicate percent of 1 kb windows in which TC is underrepresented (green), overrepresented (light green), and neutral (pink). Details are available in ‘Materials and Methods’ section.</p></caption><graphic xlink:href="nihms778642f3"></graphic></fig><fig id="F4" orientation="portrait" position="float"><label>Figure 4</label><caption><p>Extensive TC dinucleotide depletion in mucosal alpha-HPVs. The O/E ratios of TC dinucleotides in PV genomes were analyzed with genome sequences of (A) nearest genera to alpha-HPVs (one omega, seven upsilon, four omikron, two dyodelta, and one unclassified PVs; <italic>n</italic> = 15), cutaneous alpha-HPVs (<italic>n</italic> = 13), and mucosal alpha-HPVs (<italic>n</italic> = 30); and (B) nearest genera to gamma-HPVs (six pi, one dyoxi, three tau, and two unclassified PVs; <italic>n</italic> = 12) and cutaneous gamma-HPVs (<italic>n</italic> = 7). Box-and-whisker plots are shown. <italic>P</italic>-values were calculated by Student’s unpaired, two-tailed <italic>t</italic>-test. *<italic>P</italic><0.0001, ns, not significant.</p></caption><graphic xlink:href="nihms778642f4"></graphic></fig><fig id="F5" orientation="portrait" position="float"><label>Figure 5</label><caption><p>TC depletion represents an apomorphy of the alpha-PVs. (A) The Bayesian phylogenetic tree is constructed based on a concatenated E1-E2-L2-L1 nucleotide sequence alignment. The majority-rule consensus phylogenetic tree shows the evolutionary relationship between the alpha-PVs and the members of the omega, upsilon, omicron, and dyodelta-PV genera. BPV1, a delta-PV was used to root the tree. The solid branches indicate the optimized branch lengths. Dotted lines were added to the branches to facilitate visual inspection of the tree. The alpha-PV clade is highlighted in red. Posterior probability values are indicated using symbols at the nodes (circle = 1, diamond > 0.90, square > 0.80). The nodes for which the ancestral states were estimated are indicated in the figure. (B) Table showing the average (and SD) of the estimates for the ancestral TC ratios at the indicated nodes. The reconstruction process was repeated three times to ensure repeatability.</p></caption><graphic xlink:href="nihms778642f5"></graphic></fig><fig id="F6" orientation="portrait" position="float"><label>Figure 6</label><caption><p>Depletion of TC dinucleotides is associated with amino acid coding potential. (A) The O/E ratios of each dinucleotide in alpha-HPV (<italic>n</italic> = 64) genome sequences was calculated as described in <xref rid="F1" ref-type="fig">Figure 1A</xref>. (B) E1, E2, L1, and L2 sequences from alpha-HPVs were concatenated as described in ‘Materials and Methods’ section. The O/E ratios of TC dinucleotides at various codon positions were calculated as described in ‘Materials and Methods’ section. N represents any of T, C, A, or G. (C) The TC O/E ratios at NNT <underline>C</underline>NN, T<underline>C</underline>N, and NT<underline>C</underline> in the URR were calculated as a noncoding control. Box-and-whisker plots are shown with the outliers (black triangles) identified using Tukey’s method. CA and TG dinucleotides are the most over-represented dinucleotides, while NT<underline>C</underline> cytidines are significantly depleted in alpha-PVs. <italic>P</italic><0.0001 by two-way ANOVA analysis with Tukey correction.</p></caption><graphic xlink:href="nihms778642f6"></graphic></fig><fig id="F7" orientation="portrait" position="float"><label>Figure 7</label><caption><p>The E1 and E2 gene regions of alpha-PV genomes are the most significantly depleted sites for TC dinucleotides. (A) TC O/E ratios were calculated in 1 kb scanning windows using the CompSeq program as described in ‘Materials and Methods’ section. All alpha- (<italic>n</italic> = 77) and beta-PVs (<italic>n</italic> = 49) genomes were included. Each window represents 1,000 kb with a 100 bp step between two windows. The viral genomes were linearized in such a way that the first nucleotide corresponds to the first nucleotide following the L1 stop codon. (B) Genome maps of HPV16 and HPV5, adapted from PaVE, show the various ORFs and URR for representative alpha- and beta-PVs, respectively.</p></caption><graphic xlink:href="nihms778642f7"></graphic></fig></floats-group></pmc></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Ticri/CIDE/explor/EpistemeV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000141 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000141 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Ticri/CIDE
|area= EpistemeV1
|flux= Pmc
|étape= Corpus
|type= RBID
|clé= PMC:4999249
|texte= Role of the host restriction factor APOBEC3 on papillomavirus evolution
}}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:NONE" \
| HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \
| NlmPubMed2Wicri -a EpistemeV1
| This area was generated with Dilib version V0.6.31. |  |