Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Cancer, Warts, or Asymptomatic Infections: Clinical Presentation Matches Codon Usage Preferences in Human Papillomaviruses</title>
<author><name sortKey="Felez Sanchez, Marta" sort="Felez Sanchez, Marta" uniqKey="Felez Sanchez M" first="Marta" last="Félez-Sánchez">Marta Félez-Sánchez</name>
<affiliation><nlm:aff id="evv129-AFF1">Infections and Cancer Laboratory, Catalan Institute of Oncology, L'Hospitalet de Llobregat, Barcelona, Spain</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="evv129-AFF2">Virus and Cancer Laboratory. Bellvitge Institute of Biomedical Research (IDIBELL), L'Hospitalet de Llobregat, Barcelona, Spain</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Trosemeier, Jan Hendrik" sort="Trosemeier, Jan Hendrik" uniqKey="Trosemeier J" first="Jan-Hendrik" last="Trösemeier">Jan-Hendrik Trösemeier</name>
<affiliation><nlm:aff id="evv129-AFF3">Molecular Bioinformatics, Institute of Computer Science, Johann Wolfgang Goethe University, Frankfurt am Main, Germany</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="evv129-AFF4">Paul-Ehrlich-Institut, Federal Institute for Vaccines and Biomedicines, Langen, Germany</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Bedhomme, Stephanie" sort="Bedhomme, Stephanie" uniqKey="Bedhomme S" first="Stéphanie" last="Bedhomme">Stéphanie Bedhomme</name>
<affiliation><nlm:aff id="evv129-AFF1">Infections and Cancer Laboratory, Catalan Institute of Oncology, L'Hospitalet de Llobregat, Barcelona, Spain</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="evv129-AFF2">Virus and Cancer Laboratory. Bellvitge Institute of Biomedical Research (IDIBELL), L'Hospitalet de Llobregat, Barcelona, Spain</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="evv129-AFF5">Département d'Ecologie Evolutive Centre d’Ecologie Fonctionnelle et Evolutive, CNRS – UMR 5175, Montpellier, France</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Gonzalez Bravo, Maria Isabel" sort="Gonzalez Bravo, Maria Isabel" uniqKey="Gonzalez Bravo M" first="Maria Isabel" last="González-Bravo">Maria Isabel González-Bravo</name>
<affiliation><nlm:aff id="evv129-AFF6">Facultad de Economía y Empresa, University Salamanca, Salamanca, Spain</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Kamp, Christel" sort="Kamp, Christel" uniqKey="Kamp C" first="Christel" last="Kamp">Christel Kamp</name>
<affiliation><nlm:aff id="evv129-AFF4">Paul-Ehrlich-Institut, Federal Institute for Vaccines and Biomedicines, Langen, Germany</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Bravo, Ignacio G" sort="Bravo, Ignacio G" uniqKey="Bravo I" first="Ignacio G." last="Bravo">Ignacio G. Bravo</name>
<affiliation><nlm:aff id="evv129-AFF1">Infections and Cancer Laboratory, Catalan Institute of Oncology, L'Hospitalet de Llobregat, Barcelona, Spain</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="evv129-AFF2">Virus and Cancer Laboratory. Bellvitge Institute of Biomedical Research (IDIBELL), L'Hospitalet de Llobregat, Barcelona, Spain</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">26139833</idno>
<idno type="pmc">4558848</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4558848</idno>
<idno type="RBID">PMC:4558848</idno>
<idno type="doi">10.1093/gbe/evv129</idno>
<date when="2015">2015</date>
<idno type="wicri:Area/Pmc/Corpus">000005</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000005</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Cancer, Warts, or Asymptomatic Infections: Clinical Presentation Matches Codon Usage Preferences in Human Papillomaviruses</title>
<author><name sortKey="Felez Sanchez, Marta" sort="Felez Sanchez, Marta" uniqKey="Felez Sanchez M" first="Marta" last="Félez-Sánchez">Marta Félez-Sánchez</name>
<affiliation><nlm:aff id="evv129-AFF1">Infections and Cancer Laboratory, Catalan Institute of Oncology, L'Hospitalet de Llobregat, Barcelona, Spain</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="evv129-AFF2">Virus and Cancer Laboratory. Bellvitge Institute of Biomedical Research (IDIBELL), L'Hospitalet de Llobregat, Barcelona, Spain</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Trosemeier, Jan Hendrik" sort="Trosemeier, Jan Hendrik" uniqKey="Trosemeier J" first="Jan-Hendrik" last="Trösemeier">Jan-Hendrik Trösemeier</name>
<affiliation><nlm:aff id="evv129-AFF3">Molecular Bioinformatics, Institute of Computer Science, Johann Wolfgang Goethe University, Frankfurt am Main, Germany</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="evv129-AFF4">Paul-Ehrlich-Institut, Federal Institute for Vaccines and Biomedicines, Langen, Germany</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Bedhomme, Stephanie" sort="Bedhomme, Stephanie" uniqKey="Bedhomme S" first="Stéphanie" last="Bedhomme">Stéphanie Bedhomme</name>
<affiliation><nlm:aff id="evv129-AFF1">Infections and Cancer Laboratory, Catalan Institute of Oncology, L'Hospitalet de Llobregat, Barcelona, Spain</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="evv129-AFF2">Virus and Cancer Laboratory. Bellvitge Institute of Biomedical Research (IDIBELL), L'Hospitalet de Llobregat, Barcelona, Spain</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="evv129-AFF5">Département d'Ecologie Evolutive Centre d’Ecologie Fonctionnelle et Evolutive, CNRS – UMR 5175, Montpellier, France</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Gonzalez Bravo, Maria Isabel" sort="Gonzalez Bravo, Maria Isabel" uniqKey="Gonzalez Bravo M" first="Maria Isabel" last="González-Bravo">Maria Isabel González-Bravo</name>
<affiliation><nlm:aff id="evv129-AFF6">Facultad de Economía y Empresa, University Salamanca, Salamanca, Spain</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Kamp, Christel" sort="Kamp, Christel" uniqKey="Kamp C" first="Christel" last="Kamp">Christel Kamp</name>
<affiliation><nlm:aff id="evv129-AFF4">Paul-Ehrlich-Institut, Federal Institute for Vaccines and Biomedicines, Langen, Germany</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Bravo, Ignacio G" sort="Bravo, Ignacio G" uniqKey="Bravo I" first="Ignacio G." last="Bravo">Ignacio G. Bravo</name>
<affiliation><nlm:aff id="evv129-AFF1">Infections and Cancer Laboratory, Catalan Institute of Oncology, L'Hospitalet de Llobregat, Barcelona, Spain</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="evv129-AFF2">Virus and Cancer Laboratory. Bellvitge Institute of Biomedical Research (IDIBELL), L'Hospitalet de Llobregat, Barcelona, Spain</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">Genome Biology and Evolution</title>
<idno type="eISSN">1759-6653</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>Viruses rely completely on the hosts’ machinery for translation of viral transcripts. However, for most viruses infecting humans, codon usage preferences (CUPrefs) do not match those of the host. Human papillomaviruses (HPVs) are a showcase to tackle this paradox: they present a large genotypic diversity and a broad range of phenotypic presentations, from asymptomatic infections to productive lesions and cancer. By applying phylogenetic inference and dimensionality reduction methods, we demonstrate first that genes in HPVs are poorly adapted to the average human CUPrefs, the only exception being capsid genes in viruses causing productive lesions. Phylogenetic relationships between HPVs explained only a small proportion of CUPrefs variation. Instead, the most important explanatory factor for viral CUPrefs was infection phenotype, as orthologous genes in viruses with similar clinical presentation displayed similar CUPrefs. Moreover, viral genes with similar spatiotemporal expression patterns also showed similar CUPrefs. Our results suggest that CUPrefs in HPVs reflect either variations in the mutation bias or differential selection pressures depending on the clinical presentation and expression timing. We propose that poor viral CUPrefs may be central to a trade-off between strong viral gene expression and the potential for eliciting protective immune response.</p>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Akashi, H" uniqKey="Akashi H">H Akashi</name>
</author>
<author><name sortKey="Eyre Walker, A" uniqKey="Eyre Walker A">A Eyre-Walker</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Aota, S" uniqKey="Aota S">S Aota</name>
</author>
<author><name sortKey="Ikemura, T" uniqKey="Ikemura T">T Ikemura</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Aragones, L" uniqKey="Aragones L">L Aragones</name>
</author>
<author><name sortKey="Bosch, A" uniqKey="Bosch A">A Bosch</name>
</author>
<author><name sortKey="Pinto, Rm" uniqKey="Pinto R">RM Pinto</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Aragones, L" uniqKey="Aragones L">L Aragones</name>
</author>
<author><name sortKey="Guix, S" uniqKey="Guix S">S Guix</name>
</author>
<author><name sortKey="Ribes, E" uniqKey="Ribes E">E Ribes</name>
</author>
<author><name sortKey="Bosch, A" uniqKey="Bosch A">A Bosch</name>
</author>
<author><name sortKey="Pinto, Rm" uniqKey="Pinto R">RM Pinto</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bahir, I" uniqKey="Bahir I">I Bahir</name>
</author>
<author><name sortKey="Fromer, M" uniqKey="Fromer M">M Fromer</name>
</author>
<author><name sortKey="Prat, Y" uniqKey="Prat Y">Y Prat</name>
</author>
<author><name sortKey="Linial, M" uniqKey="Linial M">M Linial</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bedell, Ma" uniqKey="Bedell M">MA Bedell</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Belalov, Is" uniqKey="Belalov I">IS Belalov</name>
</author>
<author><name sortKey="Lukashev, An" uniqKey="Lukashev A">AN Lukashev</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bernard, Hu" uniqKey="Bernard H">HU Bernard</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 Felez-Sanchez</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bravo, Ig" uniqKey="Bravo I">IG Bravo</name>
</author>
<author><name sortKey="Muller, M" uniqKey="Muller M">M Müller</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bulmer, M" uniqKey="Bulmer M">M Bulmer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Castresana, J" uniqKey="Castresana J">J Castresana</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chen, L" uniqKey="Chen L">L Chen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chen, Sl" uniqKey="Chen S">SL Chen</name>
</author>
<author><name sortKey="Lee, W" uniqKey="Lee W">W Lee</name>
</author>
<author><name sortKey="Hottes, Ak" uniqKey="Hottes A">AK Hottes</name>
</author>
<author><name sortKey="Shapiro, L" uniqKey="Shapiro L">L Shapiro</name>
</author>
<author><name sortKey="Mcadams, Hh" uniqKey="Mcadams H">HH McAdams</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chen, Sl" uniqKey="Chen S">SL Chen</name>
</author>
<author><name sortKey="Tsao, Yp" uniqKey="Tsao Y">YP Tsao</name>
</author>
<author><name sortKey="Lee, Jw" uniqKey="Lee J">JW Lee</name>
</author>
<author><name sortKey="Sheu, Wc" uniqKey="Sheu W">WC Sheu</name>
</author>
<author><name sortKey="Liu, Yt" uniqKey="Liu Y">YT Liu</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chen, Z" uniqKey="Chen Z">Z Chen</name>
</author>
<author><name sortKey="Schiffman, M" uniqKey="Schiffman M">M Schiffman</name>
</author>
<author><name sortKey="Herrero, R" uniqKey="Herrero R">R Herrero</name>
</author>
<author><name sortKey="Desalle, R" uniqKey="Desalle R">R Desalle</name>
</author>
<author><name sortKey="Burk, Rd" uniqKey="Burk R">RD Burk</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cid Arregui, A" uniqKey="Cid Arregui A">A Cid-Arregui</name>
</author>
<author><name sortKey="Juarez, V" uniqKey="Juarez V">V Juarez</name>
</author>
<author><name sortKey="Zur Hausen, H" uniqKey="Zur Hausen H">H zur Hausen</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="Cladel, Nm" uniqKey="Cladel N">NM Cladel</name>
</author>
<author><name sortKey="Hu, J" uniqKey="Hu J">J Hu</name>
</author>
<author><name sortKey="Balogh, Kk" uniqKey="Balogh K">KK Balogh</name>
</author>
<author><name sortKey="Christensen, Nd" uniqKey="Christensen N">ND Christensen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cox, Tf" uniqKey="Cox T">TF Cox</name>
</author>
<author><name sortKey="Cox, Maa" uniqKey="Cox M">MAA Cox</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="De Villiers, Em" uniqKey="De Villiers E">EM de Villiers</name>
</author>
<author><name sortKey="Fauquet, C" uniqKey="Fauquet C">C Fauquet</name>
</author>
<author><name sortKey="Broker, Tr" uniqKey="Broker T">TR Broker</name>
</author>
<author><name sortKey="Bernard, Hu" uniqKey="Bernard H">HU Bernard</name>
</author>
<author><name sortKey="Zur Hausen, H" uniqKey="Zur Hausen H">H zur Hausen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Disbrow, Gl" uniqKey="Disbrow G">GL Disbrow</name>
</author>
<author><name sortKey="Sunitha, I" uniqKey="Sunitha I">I Sunitha</name>
</author>
<author><name sortKey="Baker, Cc" uniqKey="Baker C">CC Baker</name>
</author>
<author><name sortKey="Hanover, J" uniqKey="Hanover J">J Hanover</name>
</author>
<author><name sortKey="Schlegel, R" uniqKey="Schlegel R">R Schlegel</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Doorbar, J" uniqKey="Doorbar J">J Doorbar</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Doorbar, J" uniqKey="Doorbar J">J Doorbar</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Drummond, Da" uniqKey="Drummond D">DA Drummond</name>
</author>
<author><name sortKey="Wilke, Co" uniqKey="Wilke C">CO Wilke</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Duret, L" uniqKey="Duret L">L Duret</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Duret, L" uniqKey="Duret L">L Duret</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Edgar, Rc" uniqKey="Edgar R">RC Edgar</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Eriksson, A" uniqKey="Eriksson A">A Eriksson</name>
</author>
<author><name sortKey="Herron, Jr" uniqKey="Herron J">JR Herron</name>
</author>
<author><name sortKey="Yamada, T" uniqKey="Yamada T">T Yamada</name>
</author>
<author><name sortKey="Wheeler, Cm" uniqKey="Wheeler C">CM Wheeler</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Forslund, O" uniqKey="Forslund O">O Forslund</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Garcia Vallve, S" uniqKey="Garcia Vallve S">S Garcia-Vallve</name>
</author>
<author><name sortKey="Alonso, A" uniqKey="Alonso A">A Alonso</name>
</author>
<author><name sortKey="Bravo, Ig" uniqKey="Bravo I">IG Bravo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gingold, H" uniqKey="Gingold H">H Gingold</name>
</author>
<author><name sortKey="Pilpel, Y" uniqKey="Pilpel Y">Y Pilpel</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gottschling, M" uniqKey="Gottschling M">M Gottschling</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gruener, M" uniqKey="Gruener M">M Gruener</name>
</author>
<author><name sortKey="Bravo, Ig" uniqKey="Bravo I">IG Bravo</name>
</author>
<author><name sortKey="Momburg, F" uniqKey="Momburg F">F Momburg</name>
</author>
<author><name sortKey="Alonso, A" uniqKey="Alonso A">A Alonso</name>
</author>
<author><name sortKey="Tomakidi, P" uniqKey="Tomakidi P">P Tomakidi</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gu, W" uniqKey="Gu W">W Gu</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ham, J" uniqKey="Ham J">J Ham</name>
</author>
<author><name sortKey="Dostatni, N" uniqKey="Dostatni N">N Dostatni</name>
</author>
<author><name sortKey="Gauthier, Jm" uniqKey="Gauthier J">JM Gauthier</name>
</author>
<author><name sortKey="Yaniv, M" uniqKey="Yaniv M">M Yaniv</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hibma, Mh" uniqKey="Hibma M">MH Hibma</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hughes, Al" uniqKey="Hughes A">AL Hughes</name>
</author>
<author><name sortKey="Hughes, Ma" uniqKey="Hughes M">MA Hughes</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Iftner, A" uniqKey="Iftner A">A Iftner</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Iftner, T" uniqKey="Iftner T">T Iftner</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jenkins, Gm" uniqKey="Jenkins G">GM Jenkins</name>
</author>
<author><name sortKey="Holmes, Ec" uniqKey="Holmes E">EC Holmes</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Karlin, S" uniqKey="Karlin S">S Karlin</name>
</author>
<author><name sortKey="Blaisdell, Be" uniqKey="Blaisdell B">BE Blaisdell</name>
</author>
<author><name sortKey="Schachtel, Ga" uniqKey="Schachtel G">GA Schachtel</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kent, P" uniqKey="Kent P">P Kent</name>
</author>
<author><name sortKey="Kongsted, A" uniqKey="Kongsted A">A Kongsted</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kent, Wj" uniqKey="Kent W">WJ Kent</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Li, M" uniqKey="Li M">M Li</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Liu, W" uniqKey="Liu W">W Liu</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Longworth, Ms" uniqKey="Longworth M">MS Longworth</name>
</author>
<author><name sortKey="Laimins, La" uniqKey="Laimins L">LA Laimins</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Marin, A" uniqKey="Marin A">A Marin</name>
</author>
<author><name sortKey="Bertranpetit, J" uniqKey="Bertranpetit J">J Bertranpetit</name>
</author>
<author><name sortKey="Oliver, Jl" uniqKey="Oliver J">JL Oliver</name>
</author>
<author><name sortKey="Medina, Jr" uniqKey="Medina J">JR Medina</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Martincorena, I" uniqKey="Martincorena I">I Martincorena</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Moriyama, En" uniqKey="Moriyama E">EN Moriyama</name>
</author>
<author><name sortKey="Powell, Jr" uniqKey="Powell J">JR Powell</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mossadegh, N" uniqKey="Mossadegh N">N Mossadegh</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Munk, C" uniqKey="Munk C">C Munk</name>
</author>
<author><name sortKey="Willemsen, A" uniqKey="Willemsen A">A Willemsen</name>
</author>
<author><name sortKey="Bravo, Ig" uniqKey="Bravo I">IG Bravo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nakamura, Y" uniqKey="Nakamura Y">Y Nakamura</name>
</author>
<author><name sortKey="Gojobori, T" uniqKey="Gojobori T">T Gojobori</name>
</author>
<author><name sortKey="Ikemura, T" uniqKey="Ikemura T">T Ikemura</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Narechania, A" uniqKey="Narechania A">A Narechania</name>
</author>
<author><name sortKey="Terai, M" uniqKey="Terai M">M Terai</name>
</author>
<author><name sortKey="Burk, Rd" uniqKey="Burk R">RD Burk</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nindl, I" uniqKey="Nindl I">I Nindl</name>
</author>
<author><name sortKey="Gottschling, M" uniqKey="Gottschling M">M Gottschling</name>
</author>
<author><name sortKey="Stockfleth, E" uniqKey="Stockfleth E">E Stockfleth</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nobre, Rj" uniqKey="Nobre R">RJ Nobre</name>
</author>
<author><name sortKey="De Almeida, Lp" uniqKey="De Almeida L">LP de Almeida</name>
</author>
<author><name sortKey="Martins, Tc" uniqKey="Martins T">TC Martins</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Park, P" uniqKey="Park P">P Park</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Peh, Wl" uniqKey="Peh W">WL Peh</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Plotkin, Jb" uniqKey="Plotkin J">JB Plotkin</name>
</author>
<author><name sortKey="Kudla, G" uniqKey="Kudla G">G Kudla</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Protic Sabljic, M" uniqKey="Protic Sabljic M">M Protic-Sabljic</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Rancurel, C" uniqKey="Rancurel C">C Rancurel</name>
</author>
<author><name sortKey="Khosravi, M" uniqKey="Khosravi M">M Khosravi</name>
</author>
<author><name sortKey="Dunker, Ak" uniqKey="Dunker A">AK Dunker</name>
</author>
<author><name sortKey="Romero, Pr" uniqKey="Romero P">PR Romero</name>
</author>
<author><name sortKey="Karlin, D" uniqKey="Karlin D">D Karlin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Rocha, Ep" uniqKey="Rocha E">EP Rocha</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Roychoudhury, S" uniqKey="Roychoudhury S">S Roychoudhury</name>
</author>
<author><name sortKey="Mukherjee, D" uniqKey="Mukherjee D">D Mukherjee</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Samorski, R" uniqKey="Samorski R">R Samorski</name>
</author>
<author><name sortKey="Gissmann, L" uniqKey="Gissmann L">L Gissmann</name>
</author>
<author><name sortKey="Osen, W" uniqKey="Osen W">W Osen</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="Sharp, Pm" uniqKey="Sharp P">PM Sharp</name>
</author>
<author><name sortKey="Averof, M" uniqKey="Averof M">M Averof</name>
</author>
<author><name sortKey="Lloyd, At" uniqKey="Lloyd A">AT Lloyd</name>
</author>
<author><name sortKey="Matassi, G" uniqKey="Matassi G">G Matassi</name>
</author>
<author><name sortKey="Peden, Jf" uniqKey="Peden J">JF Peden</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sharp, Pm" uniqKey="Sharp P">PM Sharp</name>
</author>
<author><name sortKey="Li, Wh" uniqKey="Li W">WH Li</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sharp, Pm" uniqKey="Sharp P">PM Sharp</name>
</author>
<author><name sortKey="Li, Wh" uniqKey="Li W">WH Li</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Shields, Dc" uniqKey="Shields D">DC Shields</name>
</author>
<author><name sortKey="Sharp, Pm" uniqKey="Sharp P">PM Sharp</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sohn, Wj" uniqKey="Sohn W">WJ Sohn</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="Stenico, M" uniqKey="Stenico M">M Stenico</name>
</author>
<author><name sortKey="Lloyd, At" uniqKey="Lloyd A">AT Lloyd</name>
</author>
<author><name sortKey="Sharp, Pm" uniqKey="Sharp P">PM Sharp</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sugiyama, T" uniqKey="Sugiyama T">T Sugiyama</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Suzuki, H" uniqKey="Suzuki H">H Suzuki</name>
</author>
<author><name sortKey="Brown, Cj" uniqKey="Brown C">CJ Brown</name>
</author>
<author><name sortKey="Forney, Lj" uniqKey="Forney L">LJ Forney</name>
</author>
<author><name sortKey="Top, Em" uniqKey="Top E">EM Top</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Temiz, Na" uniqKey="Temiz N">NA Temiz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tindle, Rw" uniqKey="Tindle R">RW Tindle</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Turabelidze, A" uniqKey="Turabelidze A">A Turabelidze</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Vartanian, Jp" uniqKey="Vartanian J">JP Vartanian</name>
</author>
<author><name sortKey="Guetard, D" uniqKey="Guetard D">D Guetard</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="Vartanian, Jp" uniqKey="Vartanian J">JP Vartanian</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zhao, Kn" uniqKey="Zhao K">KN Zhao</name>
</author>
<author><name sortKey="Gu, W" uniqKey="Gu W">W Gu</name>
</author>
<author><name sortKey="Fang, Nx" uniqKey="Fang N">NX Fang</name>
</author>
<author><name sortKey="Saunders, Na" uniqKey="Saunders N">NA Saunders</name>
</author>
<author><name sortKey="Frazer, Ih" uniqKey="Frazer I">IH Frazer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zhao, Kn" uniqKey="Zhao K">KN Zhao</name>
</author>
<author><name sortKey="Liu, Wj" uniqKey="Liu W">WJ Liu</name>
</author>
<author><name sortKey="Frazer, Ih" uniqKey="Frazer I">IH Frazer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zheng, Zm" uniqKey="Zheng Z">ZM Zheng</name>
</author>
<author><name sortKey="Baker, Cc" uniqKey="Baker C">CC Baker</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zhou, J" uniqKey="Zhou J">J Zhou</name>
</author>
<author><name sortKey="Liu, Wj" uniqKey="Liu W">WJ Liu</name>
</author>
<author><name sortKey="Peng, Sw" uniqKey="Peng S">SW Peng</name>
</author>
<author><name sortKey="Sun, Xy" uniqKey="Sun X">XY Sun</name>
</author>
<author><name sortKey="Frazer, I" uniqKey="Frazer I">I Frazer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zhou, M" uniqKey="Zhou M">M Zhou</name>
</author>
<author><name sortKey="Wang, T" uniqKey="Wang T">T Wang</name>
</author>
<author><name sortKey="Fu, J" uniqKey="Fu J">J Fu</name>
</author>
<author><name sortKey="Xiao, G" uniqKey="Xiao G">G Xiao</name>
</author>
<author><name sortKey="Liu, Y" uniqKey="Liu Y">Y Liu</name>
</author>
</analytic>
</biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">Genome Biol Evol</journal-id>
<journal-id journal-id-type="iso-abbrev">Genome Biol Evol</journal-id>
<journal-id journal-id-type="publisher-id">gbe</journal-id>
<journal-id journal-id-type="hwp">gbe</journal-id>
<journal-title-group><journal-title>Genome Biology and Evolution</journal-title>
</journal-title-group>
<issn pub-type="epub">1759-6653</issn>
<publisher><publisher-name>Oxford University Press</publisher-name>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">26139833</article-id>
<article-id pub-id-type="pmc">4558848</article-id>
<article-id pub-id-type="doi">10.1093/gbe/evv129</article-id>
<article-id pub-id-type="publisher-id">evv129</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject>
</subj-group>
</article-categories>
<title-group><article-title>Cancer, Warts, or Asymptomatic Infections: Clinical Presentation Matches Codon Usage Preferences in Human Papillomaviruses</article-title>
</title-group>
<contrib-group><contrib contrib-type="author"><name><surname>Félez-Sánchez</surname>
<given-names>Marta</given-names>
</name>
<xref ref-type="aff" rid="evv129-AFF1"><sup>1</sup>
</xref>
<xref ref-type="aff" rid="evv129-AFF2"><sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Trösemeier</surname>
<given-names>Jan-Hendrik</given-names>
</name>
<xref ref-type="aff" rid="evv129-AFF3"><sup>3</sup>
</xref>
<xref ref-type="aff" rid="evv129-AFF4"><sup>4</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Bedhomme</surname>
<given-names>Stéphanie</given-names>
</name>
<xref ref-type="aff" rid="evv129-AFF1"><sup>1</sup>
</xref>
<xref ref-type="aff" rid="evv129-AFF2"><sup>2</sup>
</xref>
<xref ref-type="aff" rid="evv129-AFF5"><sup>5</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>González-Bravo</surname>
<given-names>Maria Isabel</given-names>
</name>
<xref ref-type="aff" rid="evv129-AFF6"><sup>6</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Kamp</surname>
<given-names>Christel</given-names>
</name>
<xref ref-type="aff" rid="evv129-AFF4"><sup>4</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Bravo</surname>
<given-names>Ignacio G.</given-names>
</name>
<xref ref-type="aff" rid="evv129-AFF1"><sup>1</sup>
</xref>
<xref ref-type="aff" rid="evv129-AFF2"><sup>2</sup>
</xref>
<xref ref-type="corresp" rid="evv129-COR1">*</xref>
</contrib>
<aff id="evv129-AFF1"><sup>1</sup>
Infections and Cancer Laboratory, Catalan Institute of Oncology, L'Hospitalet de Llobregat, Barcelona, Spain</aff>
<aff id="evv129-AFF2"><sup>2</sup>
Virus and Cancer Laboratory. Bellvitge Institute of Biomedical Research (IDIBELL), L'Hospitalet de Llobregat, Barcelona, Spain</aff>
<aff id="evv129-AFF3"><sup>3</sup>
Molecular Bioinformatics, Institute of Computer Science, Johann Wolfgang Goethe University, Frankfurt am Main, Germany</aff>
<aff id="evv129-AFF4"><sup>4</sup>
Paul-Ehrlich-Institut, Federal Institute for Vaccines and Biomedicines, Langen, Germany</aff>
<aff id="evv129-AFF5"><sup>5</sup>
Département d'Ecologie Evolutive Centre d’Ecologie Fonctionnelle et Evolutive, CNRS – UMR 5175, Montpellier, France</aff>
<aff id="evv129-AFF6"><sup>6</sup>
Facultad de Economía y Empresa, University Salamanca, Salamanca, Spain</aff>
</contrib-group>
<author-notes><corresp id="evv129-COR1">*Corresponding author: E-mail: <email>igbravo@iconcologia.net</email>
.</corresp>
<fn id="evv129-FN11"><p><bold>Associate editor:</bold>
Purificación López-García</p>
</fn>
</author-notes>
<pub-date pub-type="collection"><month>8</month>
<year>2015</year>
</pub-date>
<pub-date pub-type="epub"><day>01</day>
<month>7</month>
<year>2015</year>
</pub-date>
<pub-date pub-type="pmc-release"><day>01</day>
<month>7</month>
<year>2015</year>
</pub-date>
<pmc-comment> PMC Release delay is 0 months and 0 days and was based on the
. </pmc-comment>
<volume>7</volume>
<issue>8</issue>
<fpage>2117</fpage>
<lpage>2135</lpage>
<history><date date-type="accepted"><day>28</day>
<month>6</month>
<year>2015</year>
</date>
</history>
<permissions><copyright-statement>© The Author(s) 2015. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.</copyright-statement>
<copyright-year>2015</copyright-year>
<license license-type="creative-commons" xlink:href="http://creativecommons.org/licenses/by-nc/4.0/"><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 journals.permissions@oup.com</license-p>
</license>
</permissions>
<abstract><p>Viruses rely completely on the hosts’ machinery for translation of viral transcripts. However, for most viruses infecting humans, codon usage preferences (CUPrefs) do not match those of the host. Human papillomaviruses (HPVs) are a showcase to tackle this paradox: they present a large genotypic diversity and a broad range of phenotypic presentations, from asymptomatic infections to productive lesions and cancer. By applying phylogenetic inference and dimensionality reduction methods, we demonstrate first that genes in HPVs are poorly adapted to the average human CUPrefs, the only exception being capsid genes in viruses causing productive lesions. Phylogenetic relationships between HPVs explained only a small proportion of CUPrefs variation. Instead, the most important explanatory factor for viral CUPrefs was infection phenotype, as orthologous genes in viruses with similar clinical presentation displayed similar CUPrefs. Moreover, viral genes with similar spatiotemporal expression patterns also showed similar CUPrefs. Our results suggest that CUPrefs in HPVs reflect either variations in the mutation bias or differential selection pressures depending on the clinical presentation and expression timing. We propose that poor viral CUPrefs may be central to a trade-off between strong viral gene expression and the potential for eliciting protective immune response.</p>
</abstract>
<kwd-group><kwd>human viruses</kwd>
<kwd>codon usage preferences</kwd>
<kwd>mutation</kwd>
<kwd>translational selection</kwd>
<kwd>immune system</kwd>
<kwd>clinical presentation</kwd>
<kwd>genotype–phenotype</kwd>
<kwd>warts</kwd>
<kwd>cancer</kwd>
<kwd>chronic infection</kwd>
<kwd>acute infection</kwd>
</kwd-group>
<counts><page-count count="19"></page-count>
</counts>
</article-meta>
</front>
<body><sec sec-type="intro"><title>Introduction</title>
<p>Synonymous codons are not used at random (<xref rid="evv129-B2" ref-type="bibr">Aota and Ikemura 1986</xref>
; <xref rid="evv129-B68" ref-type="bibr">Shields and Sharp 1987</xref>
). Codon usage preferences (CUPrefs) vary between species, and between genes within the same genome (<xref rid="evv129-B47" ref-type="bibr">Marin et al. 1989</xref>
). CUPrefs have arisen from a complex interplay between several evolutionary processes, essentially mutation and selection (<xref rid="evv129-B11" ref-type="bibr">Bulmer 1991</xref>
). The mutational model postulates that the main factor influencing average codon usage is nucleotide composition (Guanine-Cytosine [GC] content) in the genome (<xref rid="evv129-B14" ref-type="bibr">Chen et al. 2004</xref>
). This model considers changes in synonymous codon usage neutral: It assumes that no fitness effect is associated with the preferential use of a given synonymous codon (<xref rid="evv129-B58" ref-type="bibr">Plotkin and Kudla 2011</xref>
). The selection-related model postulates coadaptation between synonymous codon usage and the translation machinery (e.g., differential transfer RNA [tRNA] abundance) to optimize translational speed and enhance translational accuracy (<xref rid="evv129-B65" ref-type="bibr">Sharp et al. 1995</xref>
; <xref rid="evv129-B26" ref-type="bibr">Duret 2000</xref>
; <xref rid="evv129-B61" ref-type="bibr">Rocha 2004</xref>
). Hence, the selection model claims that synonymous mutations can indeed influence the fitness of an organism (<xref rid="evv129-B58" ref-type="bibr">Plotkin and Kudla 2011</xref>
). The mutation model and the selection model are not mutually exclusive. In fast-growing organisms with large population sizes, such as <italic>Escherichia coli</italic>
or <italic>Saccharomyces cerevisiae</italic>
, experimental evidence supports the idea that translation selection is the main factor conditioning CUPrefs (<xref rid="evv129-B71" ref-type="bibr">Stenico et al. 1994</xref>
; <xref rid="evv129-B49" ref-type="bibr">Moriyama and Powell 1997</xref>
). In contrast, in slowly growing organisms with small population sizes, such as mammals, natural selection may be inefficient to strongly pattern CUPrefs, and its effect on codon usage remains controversial (<xref rid="evv129-B27" ref-type="bibr">Duret 2002</xref>
). Besides selection for translational efficiency and accuracy, the choice of synonymous codons may also be under the selective pressure for optimal translation kinetics, to ensure correct messenger RNA (mRNA) structure and protein folding (<xref rid="evv129-B58" ref-type="bibr">Plotkin and Kudla 2011</xref>
). Inappropriate CUPrefs may impair translation kinetics, thus leading to fitness costs associated to low quantity of functional protein, but also to waste of cellular resources incurred through accumulation of erroneous and misfolded protein, increased toxicity, and cleaning costs (<xref rid="evv129-B31" ref-type="bibr">Gingold and Pilpel 2011</xref>
).</p>
<p>Codon usage in viruses seems to be shaped both by selection and mutation. On one side, all viruses depend on host translational machinery, in particular viruses that do not encode their own tRNAs (as it is the case for human viruses), and CUPrefs in viral genes tend to match protein-specific requirements (<xref rid="evv129-B1" ref-type="bibr">Akashi and Eyre-Walker 1998</xref>
): proteins required in large amount are usually encoded by genes optimized to the host CUPrefs, while maladaptation of CUPrefs results in reduced protein production (<xref rid="evv129-B5" ref-type="bibr">Bahir et al. 2009</xref>
). On the other side, genomic GC content is often a strong predictive variable for codon usage in viruses (<xref rid="evv129-B66" ref-type="bibr">Sharp and Li 1986</xref>
; <xref rid="evv129-B41" ref-type="bibr">Karlin et al. 1990</xref>
), revealing that genome-wide mutational pressures play an important role in patterning viral CUPrefs (<xref rid="evv129-B64" ref-type="bibr">Shackelton et al. 2006</xref>
). Other studies suggest that additional selective factors such as fine-tuning selection on translation kinetics and escape from antiviral cellular responses may also underlie viral CUPrefs (<xref rid="evv129-B72" ref-type="bibr">Sugiyama et al. 2005</xref>
; <xref rid="evv129-B3" ref-type="bibr">Aragones et al. 2008</xref>
, <xref rid="evv129-B4" ref-type="bibr">2010</xref>
).</p>
<p>Papillomaviruses (PVs) are nonenveloped, double-stranded DNA viruses with a circular genome of approximately 8 kbp. PVs infect epithelia in a wide spectrum of vertebrates, at cutaneous and mucosal sites (<xref rid="evv129-B9" ref-type="bibr">Bravo and Felez-Sanchez 2015</xref>
). The PV life cycle depends on keratinocyte differentiation (<xref rid="evv129-B6" ref-type="bibr">Bedell et al. 1991</xref>
). Viral genomes are primarily present as nuclear episomes, which replicate in parallel to cell division. As the daughter cell migrates upwards and undergoes differentiation, the viral DNA is amplified and the viral expression pattern is modified, eventually leading to nonlytic virion release (<xref rid="evv129-B46" ref-type="bibr">Longworth and Laimins 2004</xref>
). PV genomes typically contain eight well-defined open reading frames (ORFs) classified as early (<italic>E1</italic>
, <italic>E2</italic>
, <italic>E4</italic>
, <italic>E5</italic>
, <italic>E6</italic>
, and <italic>E7</italic>
) and late (<italic>L1</italic>
and <italic>L2</italic>
) based on their temporal expression during viral life cycle. The early <italic>E5</italic>
, <italic>E6</italic>
, and <italic>E7</italic>
oncoproteins induce cell immortalization and transformation; <italic>E1</italic>
and <italic>E2</italic>
are associated with viral genome replication; and the <italic>E4</italic>
protein is associated with cytokeratin filament collapse. The late <italic>L1</italic>
and <italic>L2</italic>
genes encode for capsid proteins that are strongly immunogenic (<xref rid="evv129-B81" ref-type="bibr">Zheng and Baker 2006</xref>
).</p>
<p>Human PVs belong to five monophyletic genera (<ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">supplementary fig. S1</ext-link>
, <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">Supplementary Material</ext-link>
online): <italic>Alphapapillomaviruses</italic>
(AlphaPVs), <italic>Betapapillomaviruses</italic>
(BetaPVs), <italic>Gammapapillomaviruses</italic>
(GammaPVs), <italic>Mupapillomaviruses</italic>
(MuPVs), and <italic>Nupapillomaviruses</italic>
(NuPVs) (<xref rid="evv129-B8" ref-type="bibr">Bernard et al. 2010</xref>
). In PV taxonomy, two PV genomes sharing more than 60% nucleotide identity in the <italic>L1</italic>
gene belong to the same PV genera (<xref rid="evv129-B21" ref-type="bibr">de Villiers et al. 2004</xref>
). Different human PV lineages have adapted to specific epithelial niches, with different types showing differences in cell tropism, natural history of the infection, prevalence, and association with disease (<xref rid="evv129-B9" ref-type="bibr">Bravo and Felez-Sanchez 2015</xref>
). The large majority of the hitherto known human PVs, essentially Beta- and GammaPVs, cause asymptomatic infections and can be detected in healthy skin swabs or, for a reduced number of GammaPVs, also in mucosal rinses (<xref rid="evv129-B54" ref-type="bibr">Nindl et al. 2007</xref>
). Mu- and NuPVs cause conspicuous, productive cutaneous lesions usually at palmar and plantar epithelial sites (<xref rid="evv129-B8" ref-type="bibr">Bernard et al. 2010</xref>
). Finally, AlphaPVs are very diverse in terms of tropism and clinical manifestation of the disease. They include viruses with cutaneous tropism causing warts in the hands, lips, or eyelids; viruses with a very defined tropism and causing sexually transmitted warts and condylomas; and viruses causing less productive, long-lasting infections associated with certain human cancers, such as cervical cancer, other anogenital cancers, and a fraction of head and neck tumors (<xref rid="evv129-B8" ref-type="bibr">Bernard et al. 2010</xref>
). Thus, there is no sharp correspondence between clinical manifestation and phylogenetic relationships for human PVs, as members from different genera could produce similar clinical presentations (e.g., Beta- and GammaPVs essentially causing asymptomatic mucocutaneous lesions [MucCutAsym]), but certain viruses belonging to the same genera could cause different clinical presentations (e.g., AlphaPVs can cause cutaneous warts [CutW], genital warts [GenW], or other mucosal lesions [MucL]).</p>
<p>PVs do not encode for any element of the transcription or the translation machinery. They rely on the host cellular apparatus for gene expression and it would be expected that PV CUPrefs match those of the host. It is thus puzzling, however, that CUPrefs in human PVs are different from CUPrefs in humans (<xref rid="evv129-B80" ref-type="bibr">Zhao et al. 2003</xref>
; <xref rid="evv129-B10" ref-type="bibr">Bravo and Müller 2005</xref>
): the preferred codons in human PV genes are the less-preferred codons in the average human genes (<xref rid="evv129-B10" ref-type="bibr">Bravo and Müller 2005</xref>
) and show a strong bias toward codons ending with Adenine/Thymine [A/T] (<xref rid="evv129-B10" ref-type="bibr">Bravo and Müller 2005</xref>
; <xref rid="evv129-B18" ref-type="bibr">Cladel et al. 2010</xref>
). These compositional differences possibly reflect a bias in the mutation/selection evolutionary processes that still needs to be understood. Hitherto, two adaptive explanations for the biased CUPrefs in PV genes have been proposed. First, it has been suggested that PV CUPrefs have been selected for because they decrease viral protein synthesis, thereby lowering immune exposure (<xref rid="evv129-B75" ref-type="bibr">Tindle 2002</xref>
; <xref rid="evv129-B17" ref-type="bibr">Cid-Arregui et al. 2003</xref>
). Second, it has been postulated that PV CUPrefs may have evolved to differentially match the varying tRNA profile of the host cell in which viral protein actually occurs: the differentiating keratinocyte (<xref rid="evv129-B82" ref-type="bibr">Zhou et al. 1999</xref>
; <xref rid="evv129-B34" ref-type="bibr">Gu et al. 2004</xref>
; <xref rid="evv129-B4" ref-type="bibr">Aragones et al. 2010</xref>
; <xref rid="evv129-B18" ref-type="bibr">Cladel et al. 2010</xref>
).</p>
<p>We have analyzed here the CUPrefs for 156 human PVs from five distinct phylogenetic groups to determine whether variations in CUPrefs could be explained by differences in tissue tropism, association with disease, and/or timing of gene expression. Due to the high dimensionality of CUPrefs data, dimensionality reduction techniques were applied: Multidimensional scaling (MDS), correspondence analysis (CA), and cluster analyses.</p>
</sec>
<sec sec-type="materials"><title>Materials and Methods</title>
<sec><title>Human PVs Gene Sequences</title>
<p>The ORFs of all human PVs available at the Papillomavirus Episteme Database (<ext-link ext-link-type="uri" xlink:href="http://pave.niaid.nih.gov">http://pave.niaid.nih.gov</ext-link>
) were collected between March and April 2013. Using an in-house PERL script, the ORFs were examined by checking the start codon, stop codon, and internal stop codons to guarantee that only true ORFs were used. The final data set included 156 HPV types: 63 AlphaPVs, 45 BetaPVs, 45 GammaPVs, 2 MuPVs, and 1 NuPV. Names, accession numbers, and other information are detailed in <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">supplementary table S1</ext-link>
, <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">Supplementary Material</ext-link>
online.</p>
</sec>
<sec><title>Clinical Manifestations of Human PV Infections</title>
<p>Human PVs were classified according to their phenotypic clinical presentation characteristics. This classification took into account both the nature of the infection and its tropism. As for the nature of the infection, most human PVs are recovered from healthy skin and healthy mucosa, and generate unapparent, nonproductive infections. Other PVs cause highly productive infections that cause self-limited benign proliferative lesions, chiefly warts. Finally, a few human PVs cause long-lasting, low productive infections that can lead to the development of malignant proliferative lesions, essentially anogenital cancers (<xref rid="evv129-B24" ref-type="bibr">Doorbar et al. 2012</xref>
). As for tropism, human PV infections are either mucosal or cutaneous. The following four groups were defined (<ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">supplementary fig. S1</ext-link>
, <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">Supplementary Material</ext-link>
online, and <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">supplementary table S1</ext-link>
, <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">Supplementary Material</ext-link>
online): Mucocutaneous asymptomatic (MucCutAsym), including Beta and GammaPVs that cause unapparent infections; GenW group, including AlphaPVs causing proliferative lesions at mucosal sites; MucL group contained AlphaPVs causing other lesions at mucosal sites and with potential for malignisation; finally, the CutW group included Alpha-, Mu-, and NuPVs causing proliferative cutaneous lesions.</p>
</sec>
<sec><title>Codon Usage Preferences Data</title>
<p>Detailed codon composition for each genus is provided in <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">supplementary table S4</ext-link>
, <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">Supplementary Material</ext-link>
online. Patterns of synonymous codon usage were analyzed in the <italic>E1</italic>
, <italic>E2</italic>
, <italic>E4</italic>
, <italic>E6</italic>
, <italic>E7</italic>
, <italic>L1</italic>
, and <italic>L2</italic>
genes. The <italic>E5</italic>
gene was excluded from the analysis because it is absent in most human PVs. The relative frequency (RF) distribution of 59 codons (excluding Met, Trp, and stop codons) was calculated using an in-house PERL script. The abundance of each codon in a gene was calculated and pondered by a factor corresponding to the sum of all synonymous codons:
<disp-formula><mml:math id="MM1"><mml:mrow><mml:mi>R</mml:mi>
<mml:msub><mml:mi>F</mml:mi>
<mml:mrow><mml:mi>a</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac><mml:mrow><mml:msub><mml:mi>n</mml:mi>
<mml:mrow><mml:mi>a</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo>
<mml:mrow><mml:mi>c</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow><mml:msub><mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
</mml:mrow>
</mml:munderover>
<mml:mrow><mml:msub><mml:mi>n</mml:mi>
<mml:mrow><mml:mi>a</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</disp-formula>
where <italic>n</italic>
<sub>ac</sub>
is the number of events in which the c-th codon for the a-th amino acid is used, and <italic>t</italic>
<sub>a</sub>
the total number of synonymous codons that encode the ath amino acid. The final representation of the codon usage data for each gene is thus a vector of 59 positions with values between 0 and 1.</p>
<p>We calculated the pairwise CUPrefs distances as the Euclidean distances between the RF vectors of the corresponding human PVs.</p>
</sec>
<sec><title>Codon Adaptation Index</title>
<p>To analyze the relationship of CUPrefs between human PVs and humans, we employed the codon adaptation index (CAI) (<xref rid="evv129-B67" ref-type="bibr">Sharp and Li 1987</xref>
). This index evaluates the match between the CUPrefs of a particular gene and those in a reference set. In our case, the reference values were the average CUPrefs in the human genome, as retrieved from the Kazusa codon usage database, under <ext-link ext-link-type="uri" xlink:href="http://www.kazusa.or.jp/codon/cgi-bin/showcodon.cgi?species=9606">http://www.kazusa.or.jp/codon/cgi-bin/showcodon.cgi?species=9606</ext-link>
, last accessed April 2013) (<xref rid="evv129-B52" ref-type="bibr">Nakamura et al. 2000</xref>
). CAI values were calculated for all human PV genes. CAI values were also calculated for the subset of human genes differentially expressed in the epithelium. Epithelial genes were retrieved from the UCSC browser in May 2015 (<xref rid="evv129-B43" ref-type="bibr">Kent et al. 2002</xref>
) filtering skin genes by their expression (log<sub>2</sub>
(tissue/reference) with a maximum value of 5 and minimum of −5). No differences in CUPrefs were identified in genes underexpressed or overexpressed in skin compared with the average CUPrefs. We therefore chose the CAI values for the subset of epithelial genes as a reference (<ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">supplementary fig. S2</ext-link>
, <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">Supplementary Material</ext-link>
online). CAI was calculated using in-house PERL Scripts. For each gene, the output of the CAI calculation was a single value, between 0 and 1, with higher values reflecting a higher similarity in CUPrefs to the reference. The maximum value of 1 is only achieved if, for each synonymous codons set, all amino acids in the considered gene are encoded by the most used codon in the reference set.</p>
</sec>
<sec sec-type="methods"><title>Dimensionality Reduction Methods: Multidimensional Scaling and Correspondence Analysis</title>
<p>Our final data set for analysis was a matrix in which the rows correspond to the genes of one human PV genome and the columns to the 59 codons, such that each row has the codon usage information for a specific gene, in terms of relative frequencies. This data set was subjected to dimensionality reduction techniques to analyze similarities among codon usage of human PV genes, by applying MDS and CA.</p>
<p>MDS refers to a broad class of procedures that create low-dimensional representations of complex data with preservation of the similarities between data points (<xref rid="evv129-B20" ref-type="bibr">Cox TF and Cox MAA 1994</xref>
). In an <italic>n</italic>
-dimensional representation, samples with very similar codon usage profiles are displayed close together. We performed a nonmetric MDS with column wise <italic>Z</italic>
-transformation of the variables using SPSS Statistics Version 17.0 (IBM, Chicago IL). For MDS analysis, the matrix based on codon RF values was used in order to avoid biases linked to amino acid composition. In order to determine the appropriate dimensionality in which data should be scaled, we used scree plots (data not shown), which display stress as a function of dimensionality. Based on the stress values and on interpretability, we chose two dimensions as the “best” solution to plot our data.</p>
<p>CA consists in a multivariate statistical method widely used to summarize the lack of independence between objects represented through rows and columns of a matrix (here genes and codons, respectively) as a small number of derived variables, called axes. By definition, the axes are ordered according to the amount of variance in the data explained by them. Data were plotted on the first two axes with the information on the amount of variance explained in these two-dimensional representations.</p>
<p>CA of codon usage data is a widely used method in sequence analysis, which can be refined through internal correspondence analysis (ICA) to account for the variability in amino acid composition between proteins as a confounding factor when one wants to analyze synonymous codon usage variability. ICA, which is basically a double within-between-CA, has been found to be the best method in generating axes that reflect variations in synonymous codon usage (<xref rid="evv129-B73" ref-type="bibr">Suzuki et al. 2008</xref>
). This method further allows distinguishing within and between group variability with respect to genotype, gene, or clinical manifestations. ICA was run in R-3.0.1 with package ade4-1.5.2 and cross-checked with the implementations in vegan-2.0.9 and FactoMineR-1.25. We additionally used seqinr-3.0.7 for computing codon usage counts and boot1.3-11 for bootstrapping to assess the variance in each codon’s contribution to the first principal axis under sampling from a sequence population. Here, we sampled with replacement from the population of all PVs, performed the ICA, and obtained the projections of codons on the first axis of the ICA. Repeating this 1,000 times rendered a distribution of projections that informed us on how reliable the ICA analysis is with regards to sampling from a population.</p>
</sec>
<sec><title>Cluster Analysis</title>
<p>Statistical clustering was used to determine the optimal number of natural conglomerates within the data, and to classify each individual gene into one of the identified conglomerates (<xref rid="evv129-B42" ref-type="bibr">Kent and Kongsted 2012</xref>
). We applied the statistical clustering technique implemented as two-step cluster analysis (SPSS Statistics; IBM). The first step groups the cases into many small subclusters. The second step groups the subclusters into the final, optimal number of clusters, estimated using the Bayesian information criterion (BIC).</p>
</sec>
<sec><title>Phylogenetic Analysis</title>
<p>Amino acid sequences were aligned with MUSCLE (<xref rid="evv129-B84" ref-type="bibr">Edgar et al. 2004</xref>
) and back-translated into codon-aligned nucleotide sequences. Informative positions were filtered with GBLOCKS under nonstringent conditions (<xref rid="evv129-B12" ref-type="bibr">Castresana 2000</xref>
). Phylogenetic relationships were inferred in a maximum likelihood framework using RAxML v.8 (<ext-link ext-link-type="uri" xlink:href="http://www.exelixis-lab.org/">http://www.exelixis-lab.org/</ext-link>
) (<xref rid="evv129-B70" ref-type="bibr">Stamatakis 2014</xref>
) at the nucleotide level. We used the GTR+ Γ4 model, considering three partitions (one per codon position). The number of required bootstrap cycles was determined with the –autoMRE command (<xref rid="evv129-B70" ref-type="bibr">Stamatakis 2014</xref>
). Pairwise evolutionary distances between terminal taxa were estimated on the best-known maximum likelihood tree using RAxML.</p>
</sec>
<sec><title>Statistical Analysis</title>
<p>The Huber’s M-estimator calculated with R was used as robust central estimator for CAI values. The significance of the differences in CAI values between PV genes and human genes was tested with a Wilcoxon and Mann–Whitney (WMW) test, implemented in R package stats. The effect of clinical presentation and gene (and their interaction) on CAI values was tested through a two-way analysis of variance (ANOVA) followed by Tukey post hoc test (SPSS Statistics; IBM). Pearson’s correlations (implemented in R) were used to determine if there was a relationship between evolutionary and codon usage-based distances. Moreover, in order to identify major sources of variation among human PVs on the axes generated by ICA, we tested for the correlation between projections to the first principal axis with CAI and average GC content at the third position in 4-fold family codons and also for the 4-fold component in the 6-fold family codons, for each gene separately. Pearson’s correlation (<italic>r</italic>
) was calculated. The square of <italic>r</italic>
indicated the percentage of the variance in the first axis that is explained by the variance in the gene feature values.</p>
</sec>
</sec>
<sec sec-type="results"><title>Results</title>
<sec><title>Genes in Human PVs Do Not Follow Human Codon Usage Preferences</title>
<p>We first evaluated the global adaptation of CUPrefs in human PV genes to the human CUPrefs, by calculating the CAI (<xref rid="evv129-B67" ref-type="bibr">Sharp and Li 1987</xref>
). This index evaluates the match between the CUPrefs of a particular gene and those in a reference set, in our case the human average CUPrefs. CAI values were calculated for each gene in the human PV data set. Because human PVs infect epithelial cells, CAI values were also calculated for the subset of human genes expressed in the epithelium. <xref ref-type="fig" rid="evv129-F1">Figure 1</xref>
shows the cumulative frequency of CAI for human epithelial genes and for each gene in the human PVs. For all human PV genes, frequency distributions are shifted to lower CAI values compared with the reference human epithelial genes. WMW tests showed that for each of the seven human PV genes, the median CAI value is significantly lower than the median CAI value for the human epithelial genes (<italic>P</italic>
< 0.001). A Kolmogorov–Smirnov (KS) test demonstrated further that for all human PV genes the distribution of the CAI values were significantly different from that of human genes (<italic>P</italic>
< 0.001). All these results point toward a low level of adaptation of the human PV genes to the CUPrefs of their host.
<fig id="evv129-F1" orientation="portrait" position="float"><label>F<sc>ig</sc>
. 1.—</label>
<caption><p>Cumulative frequency plot showing the distribution of CAI values of HPVs and human epithelial genes. CAI values for HPV genes were calculated for every gene separately, using the human mean codon usage as a reference set. CAI values of HPV genes are represented in black. CAI values for human genes expressed at epithelial cells are represented in grey. Statistical differences in CAI values between human and HPV genes were assessed by a WMW test.</p>
</caption>
<graphic xlink:href="evv129f1p"></graphic>
</fig>
</p>
</sec>
<sec><title>Capsid Genes in Human PVs Causing Cutaneous Warts Present a Greater Adaptation to Human Codon Usage</title>
<p>In order to understand the link between CUPrefs and virus clinical presentation, we analyzed CAI values stratifying human PVs by their clinical presentation, essentially in terms of productivity of the lesion and tropism. Four groups were defined (<ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">supplementary fig. S1</ext-link>
, <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">Supplementary Material</ext-link>
online, and <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">supplementary table S1</ext-link>
, <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">Supplementary Material</ext-link>
online): MucCutAsym PVs including mucocutaneous human PVs typically associated with unapparent infections belonging to BetaPVs and GammaPVs; GenW group including AlphaPVs causing proliferative lesions at mucosal sites; MucL group including other AlphaPVs causing other lesions at mucosal sites and with potential for malignization; and finally, the cutaneous warts group (CutW) including Alpha-, Mu-, and NuPVs causing proliferative cutaneous lesions.</p>
<p>We observed that CAI values to the human average for all viral genes and for all viral clinical presentations were statistically lower than CAI values for human epithelial genes (WMW test, <italic>P</italic>
< 0.001; KS test, <italic>P</italic>
< 0.001). The only exception was the <italic>L1</italic>
capsid gene of human PVs causing cutaneous warts (WMW test, <italic>P</italic>
> 0.5). CAI cumulative distributions for each gene and clinical presentation are represented in <xref ref-type="fig" rid="evv129-F2">figure 2</xref>
central values are given in <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">supplementary table S2</ext-link>
, <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">Supplementary Material</ext-link>
online. All genes in human PVs are thus poorly adapted to the average human CUPrefs, independently of their clinical presentation, to the exception of one of the capsid genes in the very productive group of viruses causing cutaneous warts.
<fig id="evv129-F2" orientation="portrait" position="float"><label>F<sc>ig</sc>
. 2.—</label>
<caption><p>Cumulative frequency plot of CAI values for human and HPV genes stratified by clinical manifestation. CAI values of each HPV gene were computed on the basis of CUPrefs in human genes. In green are represented CAI values of mucocutaneous asymptomatic HPVs (both Beta- and GammaPVs). In red, HPVs that cause mucosal lesions. In yellow, HPVs that cause genital warts. In blue, HPVs causing cutaneous warts. In black, human epithelial genes.</p>
</caption>
<graphic xlink:href="evv129f2p"></graphic>
</fig>
</p>
<p>A two-way ANOVA with clinical presentation, gene, and their interaction as factors was performed to analyze its effect on CAI values (<xref ref-type="table" rid="evv129-T1">table 1</xref>
). This analysis revealed that all factors had a significant effect. The significant “clinical presentation × gene” interaction indicated that different genes from human PVs with different clinical presentation had different degrees of adaptation to the average human CUPrefs. Then a one-way ANOVA followed by Tukey test was performed for each gene with clinical presentation as factor (<ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">supplementary table S2</ext-link>
, <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">Supplementary Material</ext-link>
online). For the <italic>E1</italic>
, <italic>E6</italic>
, <italic>L1</italic>
, and <italic>L2</italic>
genes, human PVs causing cutaneous warts showed significantly greater adaptation to human CUPrefs (<italic>P</italic>
< 0.005, Tukey post hoc test) than those with other clinical presentations. The lowest level of adaptation in <italic>L1</italic>
and <italic>L2</italic>
was found in the MucCutAsym group (<italic>P</italic>
< 0.005, Tukey post hoc test). In contrast, for the <italic>E1</italic>
and <italic>E6</italic>
genes, the lowest level of adaptation was seen in mucosal PVs (<italic>P</italic>
< 0.005, Tukey post hoc test). For <italic>E2</italic>
, the highest level of adaptation was found in the MucCutAsym group, and the lowest in PVs causing MucLs (<italic>P</italic>
< 0.005, Tukey post hoc test) (<xref ref-type="fig" rid="evv129-F2">fig. 2</xref>
<ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">supplementary table S2</ext-link>
, <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">Supplementary Material</ext-link>
online). Finally, for <italic>E7</italic>
and <italic>E4</italic>
, all cutaneous PVs—MucCutAsym and CutW—showed higher adaptation to the average human CUPrefs than mucosal PVs—GenW and MucL (<xref ref-type="fig" rid="evv129-F2">fig. 2</xref>
<ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">supplementary table S2</ext-link>
, <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">Supplementary Material</ext-link>
online).
<table-wrap id="evv129-T1" orientation="portrait" position="float"><label>Table 1</label>
<caption><p>Effect of Clinical Manifestation and Gene on CAI</p>
</caption>
<table frame="hsides" rules="groups"><thead align="left"><tr><th rowspan="1" colspan="1">Factor</th>
<th rowspan="1" colspan="1">df</th>
<th rowspan="1" colspan="1"><italic>F</italic>
ratio</th>
<th rowspan="1" colspan="1"><italic>P</italic>
value</th>
</tr>
</thead>
<tbody align="left"><tr><td rowspan="1" colspan="1">Clinical manifestation</td>
<td rowspan="1" colspan="1">3</td>
<td rowspan="1" colspan="1">74.52</td>
<td rowspan="1" colspan="1"><0.001</td>
</tr>
<tr><td rowspan="1" colspan="1">Gene</td>
<td rowspan="1" colspan="1">6</td>
<td rowspan="1" colspan="1">33.16</td>
<td rowspan="1" colspan="1"><0.001</td>
</tr>
<tr><td rowspan="1" colspan="1">Clinical manifestation × gene</td>
<td rowspan="1" colspan="1">18</td>
<td rowspan="1" colspan="1">16.69</td>
<td rowspan="1" colspan="1"><0.001</td>
</tr>
</tbody>
</table>
<table-wrap-foot><fn id="evv129-TF1"><p>N<sc>ote</sc>
.—A two-way ANOVA was performed to assess the influence of clinical manifestation (mucocutaneous asymptomatic, mucosal lesions, genital warts, and cutaneous warts) and gene (<italic>E6</italic>
, <italic>E7</italic>
, <italic>E1</italic>
, <italic>E2</italic>
, <italic>E4</italic>
, <italic>L1</italic>
, and <italic>L2</italic>
) on CAI. Both factors and their interaction show a significant effect on CAI. df, degree of freedom.</p>
</fn>
</table-wrap-foot>
</table-wrap>
</p>
<p>We also performed a one-way ANOVA followed by Tukey test for each clinical manifestation with gene as a factor (<ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">fig. S2</ext-link>
and <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">supplementary table S2</ext-link>
, <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">Supplementary Material</ext-link>
online). For Human PVs causing CutW, GenW, and MucL, the <italic>L1</italic>
gene showed the highest level of adaptation to human CUPrefs (<italic>P</italic>
< 0.005, Tukey post hoc test). In contrast, for MucCutAsym human PVs, late genes (<italic>L1</italic>
and <italic>L2</italic>
) showed the lowest level of adaptation. The highest level of adaptation in MucCutAsym was found in <italic>E4</italic>
, whereas for GenW and MucL, the <italic>E4</italic>
gene exhibited the lowest level of adaptation (<italic>P</italic>
< 0.005, Tukey post hoc test) (<ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">supplementary fig. S3</ext-link>
, <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">Supplementary Material</ext-link>
online, and <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">supplementary table S2</ext-link>
, <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">Supplementary Material</ext-link>
online).</p>
</sec>
<sec><title>The E2 Hinge Region Shows Higher Adaptation to Human CUPrefs than the Overlapping E4 Gene</title>
<p>The <italic>E4</italic>
ORF is nested within the <italic>E2</italic>
sequence in a different reading frame, with the <italic>E4</italic>
coding sequence overlapping the so-called hinge region of the <italic>E2</italic>
gene (<xref ref-type="fig" rid="evv129-F3">fig. 3</xref>
<italic>A</italic>
). To assess the influence of the presence of an overlapping region on CUPrefs, we calculated the CAI values for the <italic>E2</italic>
hinge region, containing <italic>E4</italic>
, as well as for the nonoverlapping region (<xref ref-type="fig" rid="evv129-F3">fig. 3</xref>
<italic>B</italic>
). We performed a WMW paired test in order to compare the adaptation of these two regions for each human PV. We found that the CAI values for the hinge region were significantly higher than those for the nonoverlapping region. Moreover, the <italic>E2</italic>
hinge region also showed significantly higher degree of adaptation compared with <italic>E4</italic>
(<italic>P</italic>
< 0.005, WMW paired test). We performed the same analysis but stratifying each human PV by their clinical manifestation (<ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">supplementary fig. S4</ext-link>
, <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">Supplementary Material</ext-link>
online). The results of the WMW paired test indicated again that for all clinical manifestations the CAI values for the <italic>E2</italic>
hinge region were significantly higher than those for the nonoverlapping region (<ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">supplementary fig. S4</ext-link>
, <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">Supplementary Material</ext-link>
online). Also, for all clinical manifestations, the CAI values for <italic>E2</italic>
hinge region were significantly higher than for <italic>E4</italic>
. Finally, we found that for MucL and GenW, <italic>E4</italic>
genes were significantly less adapted than <italic>E2</italic>
genes. On the contrary, for MucCutAsym and CutW, we found the opposite pattern, with <italic>E2</italic>
being significantly less adapted than <italic>E4</italic>
(<ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">supplementary fig. S4</ext-link>
, <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">Supplementary Material</ext-link>
online).
<fig id="evv129-F3" orientation="portrait" position="float"><label>F<sc>ig</sc>
. 3.—</label>
<caption><p>(<italic>A</italic>
) Schematic map of the <italic>E2</italic>
and <italic>E4</italic>
genes. The <italic>E4</italic>
ORF overlaps the hinge region of the <italic>E2</italic>
gene. This region is a flexible, essentially disordered connector between the functionally conserved transcription activation domain and the DNA-binding domain. (<italic>B</italic>
) Cumulative frequency plot of <italic>E2</italic>
and <italic>E4</italic>
genes, and the overlapping and nonoverlapping region of both genes. CAI values of each HPV gene were computed on the basis of CUPrefs in human genes. Code for genes: Orange, <italic>E2</italic>
; brown, <italic>E4</italic>
; violet, <italic>E2</italic>
hinge region; and green, <italic>E2</italic>
nonoverlapping region with <italic>E4</italic>
.</p>
</caption>
<graphic xlink:href="evv129f3p"></graphic>
</fig>
</p>
</sec>
<sec><title>Evolutionary Distances Explain only a Low Proportion of Codon Usage Differences</title>
<p>In order to address the relationship between CUPrefs in human PVs and their evolutionary history, we analyzed the correlation between pairwise evolutionary distances and pairwise CUPrefs distances for each of the seven genes studied (<xref ref-type="fig" rid="evv129-F4">fig. 4</xref>
). Pairwise evolutionary distances were calculated from the best-known tree reconstructed by maximum likelihood techniques (see Methods). Pairwise CUPrefs distances were calculated as the Euclidean distances between relative frequencies vectors of synonymous codon usage. They represent the similarities in terms of codon usage between pairs of PVs. The correlation between both variables was significant and positive for all genes, with <italic>r</italic>
<sup>2</sup>
ranging from 0.01 (<italic>E4</italic>
) to 0.41 (<italic>E2</italic>
) with a median value of 0.19. This means that on average, evolutionary distances explained only a fifth of the variance in codon-based distance. No clear trend was obvious depending on gene expression timing, as early genes displayed both the highest and the lowest <italic>r</italic>
<sup>2</sup>
, for <italic>E2</italic>
and <italic>E4</italic>
, respectively.
<fig id="evv129-F4" orientation="portrait" position="float"><label>F<sc>ig</sc>
. 4.—</label>
<caption><p>Scatter plot of the evolutionary distances and the CUPref-based distances in the complete data set for each HPV gene. For all plots, <italic>x</italic>
-axis represents evolutionary pairwise distances and <italic>y</italic>
-axis represents CUPref-based pairwise distances. Phylogenetic distances were obtained with RAxML for every gene separately. The correlation coefficient and <italic>P</italic>
value obtained after Pearson’s analysis are reported for each bivariate analysis.</p>
</caption>
<graphic xlink:href="evv129f4p"></graphic>
</fig>
</p>
</sec>
<sec><title>Orthologous Genes in Viruses with Similar Clinical Presentation Display Similar Codon Usage Preferences</title>
<p>In order to elucidate common patterns on CUPrefs among human PVs, MDS in combination with cluster analysis was performed for each gene independently. The 59 RF variables in the CUPrefs matrix were reduced to two dimensions using MDS procedure (<xref ref-type="fig" rid="evv129-F5">fig. 5</xref>
). The closer to one another two points lie in this two-dimensional space, the more similar the two corresponding taxa are in terms of codon usage. However, MDS does not classify individuals into clusters, and it is impossible to capture all the variability of the multidimensional data in a lower dimensional display. Hence, in parallel and independently of MDS, we performed a tag-free two-step cluster analysis on the 59 relative frequencies variable matrix for each gene independently. The optimal number of clusters was also inferred blindly using the BIC. The results of the cluster analysis were incorporated into the MDS plot in <xref ref-type="fig" rid="evv129-F5">figure 5</xref>
(see also clustering analysis results in <xref ref-type="table" rid="evv129-T2">table 2</xref>
).
<fig id="evv129-F5" orientation="portrait" position="float"><label>F<sc>ig</sc>
. 5.—</label>
<caption><p>MDS plot of codon usage differences among HPVs. The different HPVs are classified by an unsupervised two-step clustering algorithm, and visualized by MDS dimension reduction. The clusters are colored yellow (cluster1), blue (cluster 2), and green (cluster 3), respectively. Cluster analysis has been conducted for each gene separately; for this reason, “cluster 1” represents a totally independent set of PVs for each gene. Each symbol represents one individual HPV, and the distance between points is proportional to the overall dissimilarity of CUPrefs. Codes for phenotypic presentation of the infection: Squares, mucocutaneous asymptomatic; dots, mucosal lesions; inverse triangles, genital warts; triangles, cutaneous warts.</p>
</caption>
<graphic xlink:href="evv129f5p"></graphic>
</fig>
<table-wrap id="evv129-T2" orientation="portrait" position="float"><label>Table 2</label>
<caption><p>Distribution of HPVs by their Clinical Presentation in Each Conglomerate Obtained by Two-Step Cluster Analysis</p>
</caption>
<table frame="hsides" rules="groups"><thead align="left"><tr><th rowspan="1" colspan="1"></th>
<th colspan="3" align="center" rowspan="1">E6<xref ref-type="table-fn" rid="evv129-TF3"><sup>a</sup>
</xref>
, %<hr></hr>
</th>
<th colspan="2" align="center" rowspan="1">E7, %<hr></hr>
</th>
<th rowspan="1" colspan="1"></th>
<th rowspan="1" colspan="1"></th>
</tr>
<tr><th rowspan="1" colspan="1"></th>
<th rowspan="1" colspan="1">Cluster1</th>
<th rowspan="1" colspan="1">Cluster2</th>
<th rowspan="1" colspan="1">Cluster3</th>
<th rowspan="1" colspan="1">Cluster1</th>
<th rowspan="1" colspan="1">Cluster2</th>
<th rowspan="1" colspan="1"></th>
<th rowspan="1" colspan="1"></th>
</tr>
</thead>
<tbody align="left"><tr><td rowspan="1" colspan="1">MucCutAsym (<italic>n</italic>
= 90)</td>
<td rowspan="1" colspan="1">97.8</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">2.2</td>
<td rowspan="1" colspan="1">94.4</td>
<td rowspan="1" colspan="1">5.6</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td rowspan="1" colspan="1">MucL (<italic>n</italic>
= 25)</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td rowspan="1" colspan="1">GenW (<italic>n</italic>
= 13)</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td rowspan="1" colspan="1">CutW (<italic>n</italic>
= 28)</td>
<td rowspan="1" colspan="1">3.6</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">96.4</td>
<td rowspan="1" colspan="1">3.6</td>
<td rowspan="1" colspan="1">96.4</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
</tr>
</tbody>
</table>
<table frame="hsides" rules="groups"><thead align="left"><tr><th rowspan="2" colspan="1"></th>
<th colspan="3" align="center" rowspan="1">E1, %<hr></hr>
</th>
<th colspan="3" align="center" rowspan="1">E2, %<hr></hr>
</th>
<th colspan="2" align="center" rowspan="1">E4, %<hr></hr>
</th>
</tr>
<tr><th rowspan="1" colspan="1">Cluster1</th>
<th rowspan="1" colspan="1">Cluster2</th>
<th rowspan="1" colspan="1">Cluster3</th>
<th rowspan="1" colspan="1">Cluster1</th>
<th rowspan="1" colspan="1">Cluster2</th>
<th rowspan="1" colspan="1">Cluster3</th>
<th rowspan="1" colspan="1">Cluster1</th>
<th rowspan="1" colspan="1">Cluster2</th>
</tr>
</thead>
<tbody align="left"><tr><td rowspan="1" colspan="1">MucCutAsym (<italic>n</italic>
= 90)</td>
<td rowspan="1" colspan="1">94.4</td>
<td rowspan="1" colspan="1">3.4</td>
<td rowspan="1" colspan="1">94.4</td>
<td rowspan="1" colspan="1">94.4</td>
<td rowspan="1" colspan="1">3.4</td>
<td rowspan="1" colspan="1">2.2</td>
<td rowspan="1" colspan="1">94.4</td>
<td rowspan="1" colspan="1">5.6</td>
</tr>
<tr><td rowspan="1" colspan="1">MucL (<italic>n</italic>
= 25)</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
</tr>
<tr><td rowspan="1" colspan="1">GenW (<italic>n</italic>
= 13)</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
</tr>
<tr><td rowspan="1" colspan="1">CutW (<italic>n</italic>
= 28)</td>
<td rowspan="1" colspan="1">3.6</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">3.6</td>
<td rowspan="1" colspan="1">3.6</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">96.4</td>
<td rowspan="1" colspan="1">3.6</td>
<td rowspan="1" colspan="1">96.4</td>
</tr>
</tbody>
</table>
<table frame="hsides" rules="groups"><thead align="left"><tr><th rowspan="1" colspan="1"></th>
<th colspan="3" align="center" rowspan="1">L2, %<hr></hr>
</th>
<th colspan="3" align="center" rowspan="1">L1, %<hr></hr>
</th>
<th rowspan="1" colspan="1"></th>
<th rowspan="1" colspan="1"></th>
</tr>
<tr><th rowspan="1" colspan="1"></th>
<th rowspan="1" colspan="1">Cluster1</th>
<th rowspan="1" colspan="1">Cluster2</th>
<th rowspan="1" colspan="1">Cluster3</th>
<th rowspan="1" colspan="1">Cluster1</th>
<th rowspan="1" colspan="1">Cluster2</th>
<th rowspan="1" colspan="1">Cluster3</th>
<th rowspan="1" colspan="1"></th>
<th rowspan="1" colspan="1"></th>
</tr>
</thead>
<tbody align="left"><tr><td rowspan="1" colspan="1">MucCutAsym (n = 90)</td>
<td rowspan="1" colspan="1">94.4</td>
<td rowspan="1" colspan="1">3.4</td>
<td rowspan="1" colspan="1">2.2</td>
<td rowspan="1" colspan="1">94.4</td>
<td rowspan="1" colspan="1">3.4</td>
<td rowspan="1" colspan="1">2.2</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td rowspan="1" colspan="1">MucL (n = 25)</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td rowspan="1" colspan="1">GenW (n = 13)</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">92.3</td>
<td rowspan="1" colspan="1">7.7</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">84.6</td>
<td rowspan="1" colspan="1">15.4</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td rowspan="1" colspan="1">CutW (n = 28)</td>
<td rowspan="1" colspan="1">3.6</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">96.4</td>
<td rowspan="1" colspan="1">3.6</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">96.4</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
</tr>
</tbody>
</table>
<table-wrap-foot><fn id="evv129-TF2"><p>Table should be read as follows, taking E7 as an example: For E7, two-step cluster analysis identifies two different clusters: Cluster 1 spans 94.4% (<italic>n</italic>
= 85) of MucCutAsym HPVs and 3.6% (<italic>n</italic>
= 1) of HPVs causing CutW. Cluster 2 spans 100% (<italic>n</italic>
= 25) of HPV causing mucosal lesions, 100% (<italic>n</italic>
= 13) of HPVs causing GenW, 96.4% (<italic>n</italic>
= 27) of HPVs causing CutW and 5.6% (<italic>n</italic>
= 5) of HPVs causing MucCutAsym infections.</p>
</fn>
<fn id="evv129-TF3"><p><sup>a</sup>
CutA in E6, <italic>n</italic>
= 87.</p>
</fn>
</table-wrap-foot>
</table-wrap>
</p>
<p>In line with our previous findings (<xref ref-type="fig" rid="evv129-F4">fig. 4</xref>
), common ancestry did not explain similarity in CUPrefs, as the cluster analysis did not identify the three main genera, Alpha-, Beta-, and GammaPVs (<ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">supplementary table S3</ext-link>
, <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">Supplementary Material</ext-link>
online), for any gene. Instead, the main factor driving codon usage-based grouping was virus clinical presentation, as MucCutAsym including Beta- and GammaPVs appeared together for all genes, both for cluster and MDS analysis. Human PVs in the MuPV genus, which cause cutaneous warts, grouped together with very distant PVs in the AlphaPV genus and associated with similar clinical manifestations. Furthermore, HPV41, the only member in the NuPV genus, which has been associated with cutaneous warts, clustered together with the MucCutAsym group. HPV4 and HPV65, classified as GammaPVs but associated to cutaneous warts, also appeared together with CutW group. The analyses also revealed that for all genes (except for <italic>E6</italic>
, due to the absence of this gene in certain GammaPVs) three GammaPVs (in species Gamma-6) grouped with the phylogenetically distant GenW group. Finally, we found that HPV32 and HPV54, belonging to the GenW group, appeared together with CutW group in <italic>L1</italic>
and <italic>L2</italic>
for HPV32 and <italic>L1</italic>
for HPV54.</p>
<p>Viruses in the AlphaPVs genus (encompassing viruses with different clinical presentations, namely GenW, CutW, and MucL) clustered together only for the <italic>E7</italic>
and <italic>E4</italic>
genes, but grouped separately for other genes, according to their clinical presentation (<xref ref-type="fig" rid="evv129-F5">fig. 5</xref>
). For <italic>E6</italic>
, the main driving factor for CUPrefs was productivity of the infections, as CutW and GenW groups clustered together, separate from MucL. On the other hand, for <italic>E1</italic>
, <italic>E2</italic>
, <italic>L2</italic>
, and <italic>L1</italic>
, the main driving factor seemed to be tropism, as the GenW and MucL groups (both with mucosal tropism) clustered separately from viruses in the CutW group (cutaneous tropism).</p>
<p>We also performed ICA, a different method of dimensionality reduction, in order to cross-validate the MDS results (<xref ref-type="fig" rid="evv129-F6">fig. 6</xref>
). ICA is a powerful CA method for the analysis of synonymous codon usage (<xref rid="evv129-B73" ref-type="bibr">Suzuki et al. 2008</xref>
). It is able to dissociate the effects of amino acid composition from the effects directly related to synonymous codon usage. Both approaches, ICA and MDS, rendered globally similar results, further sustaining the claim that clinical presentation was the main driving factor of CUPrefs in human PVs.
<fig id="evv129-F6" orientation="portrait" position="float"><label>F<sc>ig</sc>
. 6.—</label>
<caption><p>ICA plots of synonymous CUPrefs among HPVs. The different HPVs are classified by an unsupervised two-step clustering algorithm. Cluster analysis has been conducted for each gene separately; for this reason, “cluster 1” represents a totally independent set of PVs in for each gene. Each symbol represents an individual HPV, and the distance between points is proportional to the overall dissimilarity of CUPrefs. Codes for phenotypic presentation of the infection: Squares, mucocutaneous asymptomatic; dots, mucosal lesions; inverse triangles, genital warts; triangles, cutaneous warts.</p>
</caption>
<graphic xlink:href="evv129f6p"></graphic>
</fig>
</p>
</sec>
<sec><title>Differences in GC Content in the Third Position and in Codon Adaptation Index Partly Account for Codon Usage Preferences</title>
<p>In order to identify major sources of variation among genes on the axes generated by ICA of codon usage data, we assessed separately for each gene the linear correlation between projections on the first principal axis and values of CAI and of average GC content in the third codon position of 4-fold degenerate codon families and in the 4-fold component of 6-fold degenerate codon families (GC3_4) (<xref ref-type="fig" rid="evv129-F7">fig. 7</xref>
<ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">supplementary figs. S5–S11</ext-link>
, <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">Supplementary Material</ext-link>
online). Both CAI and average GC content appeared as important sources of variation in L1 and L2 genes. For <italic>E6</italic>
and <italic>E1</italic>
genes, only GC content correlated significantly with the first ICA axis. For <italic>E7</italic>
, <italic>E2</italic>
, and <italic>E4</italic>
, none of these gene features correlated with the first axis of the ICA. These results demonstrate that both gene features contribute to variations in synonymous codon usage among genes, but their contributions vary among different genes as also may vary the correlations between GC3_4 content and CAI.
<fig id="evv129-F7" orientation="portrait" position="float"><label>F<sc>ig</sc>
. 7.—</label>
<caption><p>Projections of the first principal axis in combination with either CAI or GC content at the third position of 4-fold degenerated codon families and at the 4-fold component of the 6-fold degenerated codon families (GC3_4). Codes for phenotypic presentation of the infection: Squares, mucocutaneous asymptomatic; dots, mucosal lesions; inverse triangles, genital warts; triangles, cutaneous warts.</p>
</caption>
<graphic xlink:href="evv129f7p"></graphic>
</fig>
</p>
</sec>
<sec><title>Differences in Gene Temporal Expression Largely Determine Differential Codon Usage Preferences</title>
<p>Finally, we approached the connection between CUPrefs of the different genes and the corresponding gene expression patterns, analyzing separately the three main genera, that is, Alpha-, Beta-, and GammaPVs, by means of unsupervised two-step clustering. For each genus, the optimal number of clusters was automatically determined using the BIC criterion. For every genus, gene temporal expression was the main factor driving data clustering (<xref ref-type="fig" rid="evv129-F8">fig. 8</xref>
<xref ref-type="table" rid="evv129-T3">table 3</xref>
), and genes expressed simultaneously during virus life cycle showed similar CUPrefs and clustered together. For AlphaPVs, the late genes largely differed in CUPrefs from the early genes, which showed a larger variability. Two-step clustering also clustered late stage genes <italic>L1</italic>
and <italic>L2</italic>
together for BetaPVs and GammaPVs in combination with some early stage genes (<xref ref-type="fig" rid="evv129-F8">fig. 8</xref>
). We confirmed the results obtained by MDS using ICA, which however shows subtle differences in particular for the E2 gene (<xref ref-type="fig" rid="evv129-F9">fig. 9</xref>
).
<fig id="evv129-F8" orientation="portrait" position="float"><label>F<sc>ig</sc>
. 8.—</label>
<caption><p>MDS plot of codon usage differences for HPV genes within each genus. HPVs are classified by an unsupervised two-step clustering algorithm and visualized by MDS dimension reduction. The clusters are colored yellow (cluster1), blue (cluster 2), and green (cluster 3), respectively. Codes for genes: Squares, <italic>E6</italic>
; diamonds, <italic>E7</italic>
; triangles, <italic>E1</italic>
; inverse triangles, <italic>E2</italic>
; crosses, <italic>E4</italic>
; stars, <italic>L2</italic>
; circles, <italic>L1</italic>
.</p>
</caption>
<graphic xlink:href="evv129f8p"></graphic>
</fig>
<fig id="evv129-F9" orientation="portrait" position="float"><label>F<sc>ig</sc>
. 9.—</label>
<caption><p>ICA showing variance of synonymous CUPrefs between genes of HPV genera. The difference between genes in terms of synonymous codon usage is shown for Alpha-, Beta-, and GammaPVs (from left to right). Codes for genes: Squares, E6; diamonds, E7; triangles, E1; inverse triangles, E2; crosses, E4; stars, L2; circles, L1.</p>
</caption>
<graphic xlink:href="evv129f9p"></graphic>
</fig>
<table-wrap id="evv129-T3" orientation="portrait" position="float"><label>Table 3</label>
<caption><p>Distribution of Genes Stratified by Genera in Each Conglomerate of the Two-Step Cluster Analysis</p>
</caption>
<table frame="hsides" rules="groups"><thead align="left"><tr><th rowspan="2" colspan="1"></th>
<th colspan="3" align="center" rowspan="1">AlphaPVs, % (<italic>n</italic>
= 63)<hr></hr>
</th>
<th colspan="3" align="center" rowspan="1">BetaPVs, % (<italic>n</italic>
= 45)<hr></hr>
</th>
<th colspan="3" align="center" rowspan="1">GammaPVs<xref ref-type="table-fn" rid="evv129-TF5"><sup>a</sup>
</xref>
, % (<italic>n</italic>
= 45)<hr></hr>
</th>
</tr>
<tr><th rowspan="1" colspan="1">Cluster 1</th>
<th rowspan="1" colspan="1">Cluster 2</th>
<th rowspan="1" colspan="1">Cluster 3</th>
<th rowspan="1" colspan="1">Cluster 1</th>
<th rowspan="1" colspan="1">Cluster 2</th>
<th rowspan="1" colspan="1">Cluster 3</th>
<th rowspan="1" colspan="1">Cluster 1</th>
<th rowspan="1" colspan="1">Cluster 2</th>
<th rowspan="1" colspan="1">Cluster 3</th>
</tr>
</thead>
<tbody align="left"><tr><td rowspan="1" colspan="1">E6</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">—</td>
</tr>
<tr><td rowspan="1" colspan="1">E7</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">—</td>
</tr>
<tr><td rowspan="1" colspan="1">E1</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
</tr>
<tr><td rowspan="1" colspan="1">E2</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
</tr>
<tr><td rowspan="1" colspan="1">E4</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
</tr>
<tr><td rowspan="1" colspan="1">L2</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
</tr>
<tr><td rowspan="1" colspan="1">L1</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">—</td>
<td rowspan="1" colspan="1">100</td>
<td rowspan="1" colspan="1">—</td>
</tr>
</tbody>
</table>
<table-wrap-foot><fn id="evv129-TF4"><p>N<sc>ote</sc>
.—Table should be read as follows taking AlphaPVs as an example: For AlphaPVs, two-step cluster analysis distinguishes three different clusters. Cluster 1 groups <italic>E6</italic>
, <italic>E7</italic>
, and <italic>E4</italic>
genes. Cluster 2 groups <italic>E1</italic>
and <italic>E2</italic>
genes. Finally, cluster 3 groups <italic>L1</italic>
and <italic>L2</italic>
genes.</p>
</fn>
<fn id="evv129-TF5"><p><sup>a</sup>
For E6, <italic>n</italic>
= 42.</p>
</fn>
</table-wrap-foot>
</table-wrap>
</p>
</sec>
</sec>
<sec sec-type="discussion"><title>Discussion</title>
<p>We have sought to understand the differential contribution of the different evolutionary forces shaping human PV CUPrefs. We first assessed adaptation to the host’s CUPrefs, and showed that all genes in human PVs display poor adaptation to human CUPrefs. Our null, most parsimonious, hypothesis was that CUPrefs in human PVs should be close to CUPrefs in humans, because: 1) PVs do not encode for any element involved in translation and rely completely on their host’s cell machinery to translate their genes and 2) the relationship of these viruses with their amniote hosts lineage is ancient (<xref rid="evv129-B30" ref-type="bibr">Garcia-Vallve et al. 2005</xref>
; <xref rid="evv129-B32" ref-type="bibr">Gottschling et al. 2011</xref>
; <xref rid="evv129-B9" ref-type="bibr">Bravo and Felez-Sanchez 2015</xref>
). Deviation from host’s CUPrefs may lead to inefficient translation, in terms of decreased amount of translated protein, decreased quality of the synthesized protein, and/or decreased amount of properly folded protein (<xref rid="evv129-B10" ref-type="bibr">Bravo and Müller 2005</xref>
; <xref rid="evv129-B25" ref-type="bibr">Drummond and Wilke 2008</xref>
).</p>
<p>We tried then to evaluate whether closely related PVs displayed closely related CUPrefs. Our analyses revealed however that evolutionary distances between PVs explain only a small fraction of CUPref variation among human PVs (<xref ref-type="fig" rid="evv129-F4">fig. 4</xref>
). Similar analyses in Herpesviruses had also shown that codon usage cannot be explicitly tied to species evolution (<xref rid="evv129-B62" ref-type="bibr">Roychoudhury and Mukherjee 2009</xref>
). Furthermore, we found that Euclidean distances within AlphaPVs based on CUPrefs were larger than within BetaPVs or within GammaPVs distances (<ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">supplementary fig. S12</ext-link>
, <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">Supplementary Material</ext-link>
online). We interpret that large Euclidean CUPrefs-based distances between AlphaPVs reflect the broad diversity in clinical presentations for these viruses, spanning viruses with cutaneous and mucosal tropisms. Our results suggest that clinical variables related to tropism, productivity, and immune exposure could be responsible for variation in CUPrefs in PVs. Indeed, detailed analyses by means of different methods of dimensionality reduction (MDS and ICA) showed that orthologous genes in viruses with similar clinical presentation display similar patterns of CUPrefs.</p>
<p>A major finding obtained through unsupervised cluster analysis of viral CUPrefs was an almost perfect match between groups of viruses showing similar CUPrefs and the clinical, phenotypic presentation of the infection in terms of tropism, productivity, and potential for malignization (<xref ref-type="fig" rid="evv129-F5">figs. 5</xref>
<xref ref-type="fig" rid="evv129-F6">6</xref>
). A reduced number of human PVs were initially possible exceptions to this pattern, as they presented a tissue tropism different from the rest of their sister taxa. However, a closer analysis showed that indeed viruses with exceptional tropism also displayed exceptional CUPrefs compared with their genetically close counterparts: 1) Gamma-6 PVs phylogenetically group with asymptomatic cutaneous GammaPVs, but were isolated from cervical lesions (<xref rid="evv129-B16" ref-type="bibr">Chen et al. 2007</xref>
; <xref rid="evv129-B55" ref-type="bibr">Nobre et al. 2008</xref>
), and CUPrefs of Gamma-6 PVs are similar in all six genes to those of AlphaPVs causing GenW, and not to GammaPVs; 2) two members of Gamma-1 PVs (HPV2 and HPV65) exhibited similar CUPrefs to those of CutWarts, consistent with their association with skin warts (<xref rid="evv129-B15" ref-type="bibr">Chen et al. 1993</xref>
; <xref rid="evv129-B38" ref-type="bibr">Iftner et al. 2003</xref>
); and 3) MuPVs, causing cutaneous warts, shared common CUPrefs in all genes with AlphaPVs also causing cutaneous warts, although these viruses show large phylogenetic distances. Remarkably, HPV41, the only member of NuPV genus, showed similar CUPrefs to viruses causing MucCutAsym infections. Although this virus has been associated to cutaneous warts, it has also been associated with squamous cell carcinoma lesions as it has also been the case for some BetaPVs (<xref rid="evv129-B29" ref-type="bibr">Forslund et al. 2007</xref>
).</p>
<p>We tried finally to evaluate whether gene expression pattern during the virus life cycle could also influence CUPrefs. We found a strikingly sharp pattern with genes expressed at similar stages of the viral infection cycle sharing similar CUPrefs (<xref ref-type="fig" rid="evv129-F7">figs. 7</xref>
<xref ref-type="fig" rid="evv129-F8">8</xref>
). In AlphaPVs, genes clustered separately into three groups depending on their CUPrefs. Early genes involved in replication exhibited similar patterns of codon usage, which differ from codon usage patterns of oncogenes (<italic>E4</italic>
, <italic>E6</italic>
, and <italic>E7</italic>
). Finally, structural genes (<italic>L1</italic>
and <italic>L2</italic>
) expressed in differentiated keratinocytes also shared patterns of CUPrefs. It had been described that early genes (<italic>E1</italic>
and <italic>E2</italic>
) in AlphaPVs differed in codon usage from late genes (<italic>L1</italic>
and <italic>L2</italic>
) (<xref rid="evv129-B18" ref-type="bibr">Cladel et al. 2010</xref>
). Our results thus suggest that in AlphaPVs, differential CUPrefs match differences in temporal expression. Such differences may reflect changes in tRNA availability, as it has been reported that keratinocytes express different tRNA profiles as they differentiate (<xref rid="evv129-B79" ref-type="bibr">Zhao et al. 2005</xref>
). In Beta- and GammaPVs, the gene clustering based on CUPrefs was different from that found in AlphaPVs. Life cycle in Beta- and GammaPVs has not been studied in detail and little information is available about the temporal expression of the genes from these viruses (<xref rid="evv129-B24" ref-type="bibr">Doorbar et al. 2012</xref>
). Their transcription patterns, life cycle, and other molecular characteristics are often inferred by homology with those of the best-described AlphaPVs. This might lead to an overgeneralization of the HPV biology, overlooking the fact that different PVs within a given genus are not genetically homogeneous and that different PVs infecting the same host (here human PVs within Alpha- and BetaPVs) might indeed present different biologies (<xref rid="evv129-B18" ref-type="bibr">Cladel et al. 2010</xref>
; <xref rid="evv129-B9" ref-type="bibr">Bravo and Felez-Sanchez 2015</xref>
). We propose that studying CUPrefs of PVs infecting the same host could allow inferring life-cycle characteristics that could then be experimentally explored.</p>
<p>The <italic>E4</italic>
gene is expressed following a splice event including a small number of amino acids from the <italic>E1</italic>
gene, while most of the coding sequence overlaps the hinge region of the <italic>E2</italic>
gene in a different reading frame. The <italic>E2</italic>
hinge region is a flexible connector between the functionally conserved transcription activation domain and the DNA-binding domain (<xref ref-type="fig" rid="evv129-F3">fig. 3</xref>
<italic>A</italic>
) (<xref rid="evv129-B35" ref-type="bibr">Ham et al. 1991</xref>
; <xref rid="evv129-B28" ref-type="bibr">Eriksson et al. 1999</xref>
). We have evaluated the influence of the overlapping nature of the <italic>E2</italic>
and <italic>E4</italic>
reading frames on their CUPrefs (<xref ref-type="fig" rid="evv129-F3">fig. 3</xref>
). Our analyses revealed first that for all clinical manifestations the <italic>E4</italic>
gene was less adapted to human CUPrefs than the E2 hinge region; and second that the <italic>E2</italic>
nonoverlapping region was also significantly less adapted than the hinge region (fig. 3B). This is a striking result for two main reasons. On the one hand, because previous studies on the differential evolution of these overlapping ORFs had shown that the <italic>E4</italic>
region presents an excess of synonymous mutations, compared with and excess of nonsynonymous mutation in the <italic>E2</italic>
hinge region (<xref rid="evv129-B37" ref-type="bibr">Hughes AL and Hughes MA 2005</xref>
; <xref rid="evv129-B53" ref-type="bibr">Narechania et al. 2005</xref>
). Hence, diversifying selection in the hinge region of <italic>E2</italic>
coexists with purifying selection on the overlapping region of <italic>E4</italic>
(<xref rid="evv129-B53" ref-type="bibr">Narechania et al. 2005</xref>
). In this context, our results suggest that purifying selection in <italic>E4</italic>
concurs with deadaptation to human CUPrefs, thus implying that lowering CAI values to the human average could indeed be adaptive. On the other hand because the <italic>E2</italic>
hinge region is essentially disordered but displays the highest CAI values to the human average (<xref rid="evv129-B60" ref-type="bibr">Rancurel et al. 2009</xref>
)This is in sharp contrast with previous findings showing that disordered proteins tend to be encoded by genes with poor CUPrefs (<xref rid="evv129-B83" ref-type="bibr">Zhou et al. 2015</xref>
). Overall our results for the <italic>E2</italic>
/<italic>E4</italic>
gene pair suggest that we are still far from understanding the evolutionary interplay between overlapping genes, and that CUPrefs may play a major role that may not be captured by the standard d<italic>N</italic>
/d<italic>S</italic>
analyses.</p>
<p>The evolutionary forces shaping CUPrefs are basically mutation and/or selection, and they are not always easy to disentangle. Mutational pressures are produced by differential probability of different nucleotide exchange types, leading to nucleotide composition biases and eventually shaping CUPrefs (<xref rid="evv129-B64" ref-type="bibr">Shackelton et al. 2006</xref>
; <xref rid="evv129-B7" ref-type="bibr">Belalov and Lukashev 2013</xref>
). We have shown that nucleotide composition of human PV genes differs from that of human genes, in line with previous results (<xref rid="evv129-B10" ref-type="bibr">Bravo and Müller 2005</xref>
; <xref rid="evv129-B18" ref-type="bibr">Cladel et al. 2010</xref>
). Our analysis also revealed a correlation between GC composition and CUPrefs of different viral clinical presentations for the <italic>E1</italic>
, <italic>E6</italic>
, <italic>L1</italic>
, and <italic>L2</italic>
genes (<xref ref-type="fig" rid="evv129-F7">fig. 7</xref>
). As some authors assume that GC3_4 composition reflects to some extent mutational pressure, this correlation would suggest that mutational pressures have a role patterning CUPrefs of human PVs genes, as it has been proposed for other DNA viruses (<xref rid="evv129-B64" ref-type="bibr">Shackelton et al. 2006</xref>
). However, PVs do not encode for any element of the genome replication machinery and rely instead in cellular polymerases for replication (<xref rid="evv129-B56" ref-type="bibr">Park et al. 1994</xref>
). Mutational biases associated to viral DNA replication should thus <italic>a priori</italic>
be similar to those associated with cellular DNA replication, but it does not seem to be the case. Several mutation-based explanations between viral and cellular replication may account for the observed compositional differences. First, nucleotide composition biases in the human genome reflect mutational biases during replication in the germinal line, while nucleotide composition biases in PVs reflect mutational biases during replication in somatic cells (<xref rid="evv129-B48" ref-type="bibr">Martincorena et al. 2015</xref>
). Analogously, the spectrum of mutations accumulated in human cancers is different from that arising through common ancestry (<xref rid="evv129-B74" ref-type="bibr">Temiz et al. 2015</xref>
). Second, PV DNA is the target of the APOBEC3 internal mutators, a family of cellular cytidine deaminases that introduce directional C > T substitutions (<xref rid="evv129-B77" ref-type="bibr">Vartanian et al. 2008</xref>
). The APOBEC3 locus has been under strong selection in the primate lineage, possibly reflecting an arms race between virus and host (<xref rid="evv129-B51" ref-type="bibr">Munk et al. 2012</xref>
). In this sense, the decreased GC content in human PVs, but also in most human DNA viruses, could reflect the sustained edition pressure by the APOBEC3 proteins (<xref rid="evv129-B78" ref-type="bibr">Vartanian et al. 2010</xref>
). Additionally, some experimental studies (<xref rid="evv129-B13" ref-type="bibr">Chen et al. 2010</xref>
; <xref rid="evv129-B76" ref-type="bibr">Turabelidze et al. 2014</xref>
) have shown differences in APOBEC3 expression patterns between oral mucosa and cutaneous skin, thus allowing for differential edition of PV genomes depending on their tropism. Third, viral genome replication in PVs infecting sun-exposed keratinocytes will be subject to a higher mutation rate via the error-prone mechanisms of resolution and repair of cyclobutane pyrimidine dimers (<xref rid="evv129-B59" ref-type="bibr">Protic-Sabljic et al. 1986</xref>
). Thus, nucleotide bias in human PV genes with respect to human genes may be accounted for to some extent by differences in biochemical environment and in gene expression context in infected cells compared with the germinal line.</p>
<p>Regarding adaptation, CUPrefs can increase viral fitness by directly modulating protein production and by indirectly modulating immune exposure and expression timing. The counterintuitive result of selection on PV CUPrefs is that expression of PV genes in human cell culture from the wild-type gene sequence leads to very low protein amounts, independently of the promoter and the cell line used. Instead, high gene expression levels are only achieved when the gene sequence has been “humanized” by modification of CUPrefs (<xref rid="evv129-B45" ref-type="bibr">Liu et al. 2002</xref>
; <xref rid="evv129-B17" ref-type="bibr">Cid-Arregui et al. 2003</xref>
; <xref rid="evv129-B22" ref-type="bibr">Disbrow et al. 2003</xref>
; <xref rid="evv129-B50" ref-type="bibr">Mossadegh et al. 2004</xref>
; <xref rid="evv129-B63" ref-type="bibr">Samorski et al. 2006</xref>
; <xref rid="evv129-B33" ref-type="bibr">Gruener et al. 2007</xref>
; <xref rid="evv129-B19" ref-type="bibr">Cladel et al. 2008</xref>
). We communicate here a correlation between the adaptation to human codon usage and CUPrefs of different viral clinical presentations in <italic>E6</italic>
, <italic>L1</italic>
, and <italic>L2</italic>
(<xref ref-type="fig" rid="evv129-F7">fig. 7</xref>
), suggesting an important role of adaptation in shaping CUPrefs of at least these genes in human PVs. Our analyses reveal that the only exceptions to the systematic CUPrefs maladaptation are the capsid genes (<italic>L1</italic>
and <italic>L2</italic>
) of human PVs causing cutaneous warts, which show the highest similarity to human CUPrefs. We interpret that the observed variation in mismatch between human PV genes and human CUPrefs is related to differential virus clinical presentation. Cutaneous warts are very productive lesions, and the infected keratinocytes contain a large amount of virions that are released when dead cells shed off from skin surface (<xref rid="evv129-B24" ref-type="bibr">Doorbar et al. 2012</xref>
). Productive lesions require large amounts of the capsid proteins and we propose that the increased similarity with human CUPrefs in late genes of PVs causing cutaneous warts enhances gene expression by facilitating capsid protein synthesis. For viruses in other viral families, the highest levels of adaptation of viral proteins to host’s CUPrefs are also observed for proteins appearing abundantly in the virion (<xref rid="evv129-B41" ref-type="bibr">Karlin et al. 1990</xref>
; <xref rid="evv129-B5" ref-type="bibr">Bahir et al. 2009</xref>
). Although GenW are also productive lesions, capsid genes in human PVs responsible for these infections do not show a higher adaptation to human CUPrefs. We suggest that differences on their human codon usage adaptation may arise from differences in productivity between the two wart types, as it has been demonstrated that human PVs causing GenW produce fewer virions compared with those producing cutaneous warts (<xref rid="evv129-B57" ref-type="bibr">Peh et al. 2002</xref>
). In contrast, Beta- and GammaPVs also have a mucocutaneous tropism but different clinical presentation (asymptomatic infections) and showed lower adaptation to human CUPrefs in the capsid components of the virion (<italic>L1</italic>
and <italic>L2</italic>
). CUPref maladaptation could be adaptive if decreased protein synthesis lead to a less intense host immune response (<xref rid="evv129-B75" ref-type="bibr">Tindle 2002</xref>
). Because <italic>L1</italic>
and <italic>L2</italic>
are the most immunogenic proteins (<xref rid="evv129-B36" ref-type="bibr">Hibma 2012</xref>
), we propose that limiting the expression of structural proteins by means of codon usage maladaptation allows these viruses to better escape immune surveillance for a prolonged period of time without compromising their life cycle. These findings are supported by previous reports showing that seroprevalence against Beta- and GammaPVs exhibits a delayed but long-lasting antibody response: the immune response was low in children and increased continuously with age (<xref rid="evv129-B39" ref-type="bibr">Iftner et al. 2010</xref>
). In human PVs causing MucCutAsym infections, the highest level of adaptation to human CUPrefs was found in the <italic>E4</italic>
gene. In some cutaneous infections, <italic>E4</italic>
can be expressed at higher levels than the virion coat proteins, and can account for as much as 30% of the total protein content (<xref rid="evv129-B23" ref-type="bibr">Doorbar 2013</xref>
).</p>
<p>We hypothesize that an evolutionary trade-off exists in virus clinical presentation between a potential for strong gene expression and a potential for eliciting strong immune responses. Modulation of viral CUPrefs with respect to the host’s CUPrefs may help push the equilibrium in one direction or another. For viruses associated to chronic infections, such as human PVs, the adaptive strategy could thus be to sacrifice virion productivity to avoid the generation of strong, protective immune responses, resulting in long-lasting infections and allowing for reinfection of a previously infected host. For viruses associated to acute infections, on the contrary, large virion production is accompanied by induction of a strong immune response that may eventually render the infected host nonsusceptible to subsequent reinfections by closely related viruses. Experimental evidence of innate immune activity of the schlafen 11 gene against viral infections further sustains our hypothesis. Expression of schlafen 11 is triggered by cellular exposure to interferon, as a response to viral infections (<xref rid="evv129-B69" ref-type="bibr">Sohn et al. 2007</xref>
). The activity of schlafen 11 protein is to selectively inhibit translation from mRNAs enriched in A/T-ending codons (<xref rid="evv129-B44" ref-type="bibr">Li et al. 2012</xref>
), and many viruses infecting mammals are enriched precisely in A/T-ending codons (<xref rid="evv129-B40" ref-type="bibr">Jenkins and Holmes 2003</xref>
; <xref rid="evv129-B64" ref-type="bibr">Shackelton et al. 2006</xref>
). The result is that, as a response against a viral infection, the cellular machinery shuts down specifically the translation of transcripts that are potentially of viral origin, using CUPrefs as a guide for pinpointing viral transcripts.</p>
<p>In summary, we have presented here a thorough analysis of CUPrefs in human PV genes, connecting codon preferences with virus infection clinical presentation and with gene expression patterns. We have shown that, for viruses with a well-characterized infection cycle, genes expressed simultaneously tend to show similar CUPrefs. Furthermore, closely related viruses did not necessarily display closely related CUPrefs, while orthologous genes in distantly related viruses but with similar tropism tend to show similar CUPrefs. Finally, we propose that modulation of viral CUPrefs, as a result of differential mutation and/or selection pressures, may have an adaptive value, as they may strongly condition expression efficiency, virion production, immune exposure, and propensity toward chronic/acute virus lifestyle. Comparative research, with insight into the different life-history traits of virus lifestyle and not remaining merely on descriptions of preferences, will be required to elucidate the role and the evolutionary forces fuelling the evolution of viral CUPrefs.</p>
</sec>
<sec sec-type="supplementary-material"><title>Supplementary Material</title>
<p><ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">Supplementary tables S1–S4</ext-link>
and <ext-link ext-link-type="uri" xlink:href="http://gbe.oxfordjournals.org/lookup/suppl/doi:10.1093/gbe/evv129/-/DC1">files S1–S12</ext-link>
are available at <italic>Genome Biology and Evolution</italic>
online (<ext-link ext-link-type="uri" xlink:href="http://www.gbe.oxfordjournals.org/">http://www.gbe.oxfordjournals.org/</ext-link>
).</p>
<supplementary-material id="PMC_1" content-type="local-data"><caption><title>Supplementary Data</title>
</caption>
<media mimetype="text" mime-subtype="html" xlink:href="supp_7_8_2117__index.html"></media>
<media xlink:role="associated-file" mimetype="application" mime-subtype="pdf" xlink:href="supp_evv129_Draft_v29_GBE_codon_usage_HPV_supplementary.pdf"></media>
</supplementary-material>
</sec>
</body>
<back><ack><title>Acknowledgments</title>
<p>This work was supported by the <funding-source>disappeared Spanish Ministry for Science and Innovation</funding-source>
(grant <award-id>CGL2010-16713</award-id>
to I.G.B.), and by grants by the <funding-source>IDIBELL</funding-source>
(PhD fellowship to M.F.S.), by the Dexeus Foundation for Women's Health (to I.G.B.) and by the Adolf Messer Stiftung (personal grant to J.H.T.).</p>
</ack>
<ref-list><title>Literature Cited</title>
<ref id="evv129-B1"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Akashi</surname>
<given-names>H</given-names>
</name>
<name><surname>Eyre-Walker</surname>
<given-names>A</given-names>
</name>
</person-group>
<year>1998</year>
<article-title>Translational selection and molecular evolution</article-title>
. <source>Curr Opin Genet Dev.</source>
<volume>8</volume>
:<fpage>688</fpage>
–<lpage>693</lpage>
.<pub-id pub-id-type="pmid">9914211</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B2"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Aota</surname>
<given-names>S</given-names>
</name>
<name><surname>Ikemura</surname>
<given-names>T</given-names>
</name>
</person-group>
<year>1986</year>
<article-title>Diversity in G + C content at the third position of codons in vertebrate genes and its cause</article-title>
. <source>Nucleic Acids Res.</source>
<volume>14</volume>
:<fpage>6345</fpage>
–<lpage>6355</lpage>
.<pub-id pub-id-type="pmid">3748815</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B3"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Aragones</surname>
<given-names>L</given-names>
</name>
<name><surname>Bosch</surname>
<given-names>A</given-names>
</name>
<name><surname>Pinto</surname>
<given-names>RM</given-names>
</name>
</person-group>
<year>2008</year>
<article-title>Hepatitis A virus mutant spectra under the selective pressure of monoclonal antibodies: codon usage constraints limit capsid variability</article-title>
. <source>J Virol.</source>
<volume>82</volume>
:<fpage>1688</fpage>
–<lpage>1700</lpage>
.<pub-id pub-id-type="pmid">18057242</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B4"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Aragones</surname>
<given-names>L</given-names>
</name>
<name><surname>Guix</surname>
<given-names>S</given-names>
</name>
<name><surname>Ribes</surname>
<given-names>E</given-names>
</name>
<name><surname>Bosch</surname>
<given-names>A</given-names>
</name>
<name><surname>Pinto</surname>
<given-names>RM</given-names>
</name>
</person-group>
<year>2010</year>
<article-title>Fine-tuning translation kinetics selection as the driving force of codon usage bias in the hepatitis A virus capsid</article-title>
. <source>PLoS Pathog.</source>
<volume>6</volume>
:<fpage>e1000797</fpage>
.<pub-id pub-id-type="pmid">20221432</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B5"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bahir</surname>
<given-names>I</given-names>
</name>
<name><surname>Fromer</surname>
<given-names>M</given-names>
</name>
<name><surname>Prat</surname>
<given-names>Y</given-names>
</name>
<name><surname>Linial</surname>
<given-names>M</given-names>
</name>
</person-group>
<year>2009</year>
<article-title>Viral adaptation to host: a proteome-based analysis of codon usage and amino acid preferences</article-title>
. <source>Mol Syst Biol.</source>
<volume>5</volume>
:<fpage>311</fpage>
.<pub-id pub-id-type="pmid">19888206</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B6"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bedell</surname>
<given-names>MA</given-names>
</name>
<etal></etal>
</person-group>
<year>1991</year>
<article-title>Amplification of human papillomavirus genomes in vitro is dependent on epithelial differentiation</article-title>
. <source>J Virol.</source>
<volume>65</volume>
:<fpage>2254</fpage>
–<lpage>2260</lpage>
.<pub-id pub-id-type="pmid">1850010</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B7"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Belalov</surname>
<given-names>IS</given-names>
</name>
<name><surname>Lukashev</surname>
<given-names>AN</given-names>
</name>
</person-group>
<year>2013</year>
<article-title>Causes and implications of codon usage bias in RNA viruses</article-title>
. <source>PLoS One</source>
<volume>8</volume>
:<fpage>e56642</fpage>
.<pub-id pub-id-type="pmid">23451064</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B8"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bernard</surname>
<given-names>HU</given-names>
</name>
<etal></etal>
</person-group>
<year>2010</year>
<article-title>Classification of papillomaviruses (PVs) based on 189 PV types and proposal of taxonomic amendments</article-title>
. <source>Virology</source>
<volume>401</volume>
:<fpage>70</fpage>
–<lpage>79</lpage>
.<pub-id pub-id-type="pmid">20206957</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B9"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bravo</surname>
<given-names>IG</given-names>
</name>
<name><surname>Felez-Sanchez</surname>
<given-names>M</given-names>
</name>
</person-group>
<year>2015</year>
<article-title>Papillomaviruses: viral evolution, cancer and evolutionary medicine</article-title>
. <source>Evol Med Public Health.</source>
<volume>2015</volume>
:<fpage>32</fpage>
–<lpage>51</lpage>
.<pub-id pub-id-type="pmid">25634317</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B10"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bravo</surname>
<given-names>IG</given-names>
</name>
<name><surname>Müller</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 Rep.</source>
<volume>16</volume>
:<fpage>63</fpage>
–<lpage>72</lpage>
.</mixed-citation>
</ref>
<ref id="evv129-B11"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bulmer</surname>
<given-names>M</given-names>
</name>
</person-group>
<year>1991</year>
<article-title>The selection-mutation-drift theory of synonymous codon usage</article-title>
. <source>Genetics</source>
<volume>129</volume>
:<fpage>897</fpage>
–<lpage>907</lpage>
.<pub-id pub-id-type="pmid">1752426</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B12"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Castresana</surname>
<given-names>J</given-names>
</name>
</person-group>
<year>2000</year>
<article-title>Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis</article-title>
. <source>Mol Biol Evol.</source>
<volume>17</volume>
:<fpage>540</fpage>
–<lpage>552</lpage>
.<pub-id pub-id-type="pmid">10742046</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B13"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname>
<given-names>L</given-names>
</name>
<etal></etal>
</person-group>
<year>2010</year>
<article-title>Positional differences in the wound transcriptome of skin and oral mucosa</article-title>
. <source>BMC Genomics</source>
<volume>11</volume>
:<fpage>471</fpage>
.<pub-id pub-id-type="pmid">20704739</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B14"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname>
<given-names>SL</given-names>
</name>
<name><surname>Lee</surname>
<given-names>W</given-names>
</name>
<name><surname>Hottes</surname>
<given-names>AK</given-names>
</name>
<name><surname>Shapiro</surname>
<given-names>L</given-names>
</name>
<name><surname>McAdams</surname>
<given-names>HH</given-names>
</name>
</person-group>
<year>2004</year>
<article-title>Codon usage between genomes is constrained by genome-wide mutational processes</article-title>
. <source>Proc Natl Acad Sci U S A.</source>
<volume>101</volume>
:<fpage>3480</fpage>
–<lpage>3485</lpage>
.<pub-id pub-id-type="pmid">14990797</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B15"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname>
<given-names>SL</given-names>
</name>
<name><surname>Tsao</surname>
<given-names>YP</given-names>
</name>
<name><surname>Lee</surname>
<given-names>JW</given-names>
</name>
<name><surname>Sheu</surname>
<given-names>WC</given-names>
</name>
<name><surname>Liu</surname>
<given-names>YT</given-names>
</name>
</person-group>
<year>1993</year>
<article-title>Characterization and analysis of human papillomaviruses of skin warts</article-title>
. <source>Arch Dermatol Res.</source>
<volume>285</volume>
:<fpage>460</fpage>
–<lpage>465</lpage>
.<pub-id pub-id-type="pmid">8274034</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B16"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname>
<given-names>Z</given-names>
</name>
<name><surname>Schiffman</surname>
<given-names>M</given-names>
</name>
<name><surname>Herrero</surname>
<given-names>R</given-names>
</name>
<name><surname>Desalle</surname>
<given-names>R</given-names>
</name>
<name><surname>Burk</surname>
<given-names>RD</given-names>
</name>
</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>453</lpage>
.<pub-id pub-id-type="pmid">17125811</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B17"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cid-Arregui</surname>
<given-names>A</given-names>
</name>
<name><surname>Juarez</surname>
<given-names>V</given-names>
</name>
<name><surname>zur Hausen</surname>
<given-names>H</given-names>
</name>
</person-group>
<year>2003</year>
<article-title>A synthetic E7 gene of human papillomavirus type 16 that yields enhanced expression of the protein in mammalian cells and is useful for DNA immunization studies</article-title>
. <source>J Virol.</source>
<volume>77</volume>
:<fpage>4928</fpage>
–<lpage>4937</lpage>
.<pub-id pub-id-type="pmid">12663798</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B18"><mixed-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>340</lpage>
.<pub-id pub-id-type="pmid">20157772</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B19"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cladel</surname>
<given-names>NM</given-names>
</name>
<name><surname>Hu</surname>
<given-names>J</given-names>
</name>
<name><surname>Balogh</surname>
<given-names>KK</given-names>
</name>
<name><surname>Christensen</surname>
<given-names>ND</given-names>
</name>
</person-group>
<year>2008</year>
<article-title>CRPV genomes with synonymous codon optimizations in the CRPV E7 gene show phenotypic differences in growth and altered immunity upon E7 vaccination</article-title>
. <source>PLoS One</source>
<volume>3</volume>
:<fpage>e2947</fpage>
.<pub-id pub-id-type="pmid">18698362</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B20"><mixed-citation publication-type="book"><person-group person-group-type="author"><name><surname>Cox</surname>
<given-names>TF</given-names>
</name>
<name><surname>Cox</surname>
<given-names>MAA</given-names>
</name>
</person-group>
<year>1994</year>
<source>Multidimensional scaling</source>
. <publisher-loc>London</publisher-loc>
: <publisher-name>Chapman & Hall</publisher-name>
.</mixed-citation>
</ref>
<ref id="evv129-B21"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>de Villiers</surname>
<given-names>EM</given-names>
</name>
<name><surname>Fauquet</surname>
<given-names>C</given-names>
</name>
<name><surname>Broker</surname>
<given-names>TR</given-names>
</name>
<name><surname>Bernard</surname>
<given-names>HU</given-names>
</name>
<name><surname>zur Hausen</surname>
<given-names>H</given-names>
</name>
</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>
</mixed-citation>
</ref>
<ref id="evv129-B22"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Disbrow</surname>
<given-names>GL</given-names>
</name>
<name><surname>Sunitha</surname>
<given-names>I</given-names>
</name>
<name><surname>Baker</surname>
<given-names>CC</given-names>
</name>
<name><surname>Hanover</surname>
<given-names>J</given-names>
</name>
<name><surname>Schlegel</surname>
<given-names>R</given-names>
</name>
</person-group>
<year>2003</year>
<article-title>Codon optimization of the HPV-16 E5 gene enhances protein expression</article-title>
. <source>Virology</source>
<volume>311</volume>
:<fpage>105</fpage>
–<lpage>114</lpage>
.<pub-id pub-id-type="pmid">12832208</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B23"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Doorbar</surname>
<given-names>J</given-names>
</name>
</person-group>
<year>2013</year>
<article-title>The E4 protein; structure, function and patterns of expression</article-title>
. <source>Virology</source>
<volume>445</volume>
:<fpage>80</fpage>
–<lpage>98</lpage>
.<pub-id pub-id-type="pmid">24016539</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B24"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Doorbar</surname>
<given-names>J</given-names>
</name>
<etal></etal>
</person-group>
<year>2012</year>
<article-title>The biology and life-cycle of human papillomaviruses</article-title>
. <source>Vaccine</source>
<volume>30</volume>
(<issue>Suppl 5</issue>
):<fpage>F55</fpage>
–<lpage>F70</lpage>
.<pub-id pub-id-type="pmid">23199966</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B25"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Drummond</surname>
<given-names>DA</given-names>
</name>
<name><surname>Wilke</surname>
<given-names>CO</given-names>
</name>
</person-group>
<year>2008</year>
<article-title>Mistranslation-induced protein misfolding as a dominant constraint on coding-sequence evolution</article-title>
. <source>Cell</source>
<volume>134</volume>
:<fpage>341</fpage>
–<lpage>352</lpage>
.<pub-id pub-id-type="pmid">18662548</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B26"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Duret</surname>
<given-names>L</given-names>
</name>
</person-group>
<year>2000</year>
<article-title>tRNA gene number and codon usage in the <italic>C. elegans</italic>
genome are co-adapted for optimal translation of highly expressed genes</article-title>
. <source>Trends Genet.</source>
<volume>16</volume>
:<fpage>287</fpage>
–<lpage>289</lpage>
.<pub-id pub-id-type="pmid">10858656</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B27"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Duret</surname>
<given-names>L</given-names>
</name>
</person-group>
<year>2002</year>
<article-title>Detecting genomic features under weak selective pressure: the example of codon usage in animals and plants</article-title>
. <source>Bioinformatics</source>
<volume>18</volume>
(<issue>Suppl 2</issue>
):<fpage>S91</fpage>
.<pub-id pub-id-type="pmid">12385989</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B84"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Edgar</surname>
<given-names>RC</given-names>
</name>
</person-group>
<year>2004</year>
<article-title>MUSCLE: multiple sequence alignment with high accuracy and high throughput</article-title>
. <source>Nucleic Acids Res.</source>
<volume>32</volume>
(<issue>5</issue>
):<fpage>1792</fpage>
–<lpage>1797</lpage>
.<pub-id pub-id-type="pmid">15034147</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B28"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Eriksson</surname>
<given-names>A</given-names>
</name>
<name><surname>Herron</surname>
<given-names>JR</given-names>
</name>
<name><surname>Yamada</surname>
<given-names>T</given-names>
</name>
<name><surname>Wheeler</surname>
<given-names>CM</given-names>
</name>
</person-group>
<year>1999</year>
<article-title>Human papillomavirus type 16 variant lineages characterized by nucleotide sequence analysis of the E5 coding segment and the E2 hinge region</article-title>
. <source>J Gen Virol.</source>
<volume>80</volume>
(<issue>Pt 3</issue>
):<fpage>595</fpage>
–<lpage>600</lpage>
.<pub-id pub-id-type="pmid">10091997</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B29"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Forslund</surname>
<given-names>O</given-names>
</name>
<etal></etal>
</person-group>
<year>2007</year>
<article-title>Cutaneous human papillomaviruses found in sun-exposed skin: beta-papillomavirus species 2 predominates in squamous cell carcinoma</article-title>
. <source>J Infect Dis.</source>
<volume>196</volume>
:<fpage>876</fpage>
–<lpage>883</lpage>
.<pub-id pub-id-type="pmid">17703418</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B30"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Garcia-Vallve</surname>
<given-names>S</given-names>
</name>
<name><surname>Alonso</surname>
<given-names>A</given-names>
</name>
<name><surname>Bravo</surname>
<given-names>IG</given-names>
</name>
</person-group>
<year>2005</year>
<article-title>Papillomaviruses: different genes have different histories</article-title>
. <source>Trends Microbiol.</source>
<volume>13</volume>
:<fpage>514</fpage>
–<lpage>521</lpage>
.<pub-id pub-id-type="pmid">16181783</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B31"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gingold</surname>
<given-names>H</given-names>
</name>
<name><surname>Pilpel</surname>
<given-names>Y</given-names>
</name>
</person-group>
<year>2011</year>
<article-title>Determinants of translation efficiency and accuracy</article-title>
. <source>Mol Syst Biol.</source>
<volume>7</volume>
:<fpage>481</fpage>
.<pub-id pub-id-type="pmid">21487400</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B32"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gottschling</surname>
<given-names>M</given-names>
</name>
<etal></etal>
</person-group>
<year>2011</year>
<article-title>Quantifying the phylodynamic forces driving papillomavirus evolution</article-title>
. <source>Mol Biol Evol.</source>
<volume>28</volume>
:<fpage>2101</fpage>
–<lpage>2113</lpage>
.<pub-id pub-id-type="pmid">21285031</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B33"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gruener</surname>
<given-names>M</given-names>
</name>
<name><surname>Bravo</surname>
<given-names>IG</given-names>
</name>
<name><surname>Momburg</surname>
<given-names>F</given-names>
</name>
<name><surname>Alonso</surname>
<given-names>A</given-names>
</name>
<name><surname>Tomakidi</surname>
<given-names>P</given-names>
</name>
</person-group>
<year>2007</year>
<article-title>The E5 protein of the human papillomavirus type 16 down-regulates HLA-I surface expression in calnexin-expressing but not in calnexin-deficient cells</article-title>
. <source>Virol J.</source>
<volume>4</volume>
:<fpage>116</fpage>
.<pub-id pub-id-type="pmid">17971213</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B34"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gu</surname>
<given-names>W</given-names>
</name>
<etal></etal>
</person-group>
<year>2004</year>
<article-title>tRNASer(CGA) differentially regulates expression of wild-type and codon-modified papillomavirus L1 genes</article-title>
. <source>Nucleic Acids Res.</source>
<volume>32</volume>
:<fpage>4448</fpage>
–<lpage>4461</lpage>
.<pub-id pub-id-type="pmid">15319446</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B35"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ham</surname>
<given-names>J</given-names>
</name>
<name><surname>Dostatni</surname>
<given-names>N</given-names>
</name>
<name><surname>Gauthier</surname>
<given-names>JM</given-names>
</name>
<name><surname>Yaniv</surname>
<given-names>M</given-names>
</name>
</person-group>
<year>1991</year>
<article-title>The papillomavirus E2 protein: a factor with many talents</article-title>
. <source>Trends Biochem Sci.</source>
<volume>16</volume>
:<fpage>440</fpage>
–<lpage>444</lpage>
.<pub-id pub-id-type="pmid">1663669</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B36"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hibma</surname>
<given-names>MH</given-names>
</name>
</person-group>
<year>2012</year>
<article-title>The immune response to papillomavirus during infection persistence and regression</article-title>
. <source>Open Virol J.</source>
<volume>6</volume>
:<fpage>241</fpage>
–<lpage>248</lpage>
.<pub-id pub-id-type="pmid">23341859</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B37"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hughes</surname>
<given-names>AL</given-names>
</name>
<name><surname>Hughes</surname>
<given-names>MA</given-names>
</name>
</person-group>
<year>2005</year>
<article-title>Patterns of nucleotide difference in overlapping and non-overlapping reading frames of papillomavirus genomes</article-title>
. <source>Virus Res.</source>
<volume>113</volume>
:<fpage>81</fpage>
–<lpage>88</lpage>
.<pub-id pub-id-type="pmid">15913825</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B38"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Iftner</surname>
<given-names>A</given-names>
</name>
<etal></etal>
</person-group>
<year>2003</year>
<article-title>The prevalence of human papillomavirus genotypes in nonmelanoma skin cancers of nonimmunosuppressed individuals identifies high-risk genital types as possible risk factors</article-title>
. <source>Cancer Res.</source>
<volume>63</volume>
:<fpage>7515</fpage>
–<lpage>7519</lpage>
.<pub-id pub-id-type="pmid">14612553</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B39"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Iftner</surname>
<given-names>T</given-names>
</name>
<etal></etal>
</person-group>
<year>2010</year>
<article-title>Prevalence of low-risk and high-risk types of human papillomavirus and other risk factors for HPV infection in Germany within different age groups in women up to 30 years of age: an epidemiological observational study</article-title>
. <source>J Med Virol.</source>
<volume>82</volume>
:<fpage>1928</fpage>
–<lpage>1939</lpage>
.<pub-id pub-id-type="pmid">20872721</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B40"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jenkins</surname>
<given-names>GM</given-names>
</name>
<name><surname>Holmes</surname>
<given-names>EC</given-names>
</name>
</person-group>
<year>2003</year>
<article-title>The extent of codon usage bias in human RNA viruses and its evolutionary origin</article-title>
. <source>Virus Res.</source>
<volume>92</volume>
:<fpage>1</fpage>
–<lpage>7</lpage>
.<pub-id pub-id-type="pmid">12606071</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B41"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Karlin</surname>
<given-names>S</given-names>
</name>
<name><surname>Blaisdell</surname>
<given-names>BE</given-names>
</name>
<name><surname>Schachtel</surname>
<given-names>GA</given-names>
</name>
</person-group>
<year>1990</year>
<article-title>Contrasts in codon usage of latent versus productive genes of Epstein-Barr virus: data and hypotheses</article-title>
. <source>J Virol.</source>
<volume>64</volume>
:<fpage>4264</fpage>
–<lpage>4273</lpage>
.<pub-id pub-id-type="pmid">2166815</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B42"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kent</surname>
<given-names>P</given-names>
</name>
<name><surname>Kongsted</surname>
<given-names>A</given-names>
</name>
</person-group>
<year>2012</year>
<article-title>Identifying clinical course patterns in SMS data using cluster analysis</article-title>
. <source>Chiropr Man Therap.</source>
<volume>20</volume>
:<fpage>20</fpage>
.</mixed-citation>
</ref>
<ref id="evv129-B43"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kent</surname>
<given-names>WJ</given-names>
</name>
<etal></etal>
</person-group>
<year>2002</year>
<article-title>The human genome browser at UCSC</article-title>
. <source>Genome Res.</source>
<volume>12</volume>
:<fpage>996</fpage>
–<lpage>1006</lpage>
.<pub-id pub-id-type="pmid">12045153</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B44"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname>
<given-names>M</given-names>
</name>
<etal></etal>
</person-group>
<year>2012</year>
<article-title>Codon-usage-based inhibition of HIV protein synthesis by human schlafen 11</article-title>
. <source>Nature</source>
<volume>491</volume>
:<fpage>125</fpage>
–<lpage>128</lpage>
.<pub-id pub-id-type="pmid">23000900</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B45"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname>
<given-names>W</given-names>
</name>
<etal></etal>
</person-group>
<year>2002</year>
<article-title>Codon modified human papillomavirus type 16 E7 DNA vaccine enhances cytotoxic T-lymphocyte induction and anti-tumour activity</article-title>
. <source>Virology</source>
<volume>301</volume>
:<fpage>43</fpage>
–<lpage>52</lpage>
.<pub-id pub-id-type="pmid">12359445</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B46"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Longworth</surname>
<given-names>MS</given-names>
</name>
<name><surname>Laimins</surname>
<given-names>LA</given-names>
</name>
</person-group>
<year>2004</year>
<article-title>Pathogenesis of human papillomaviruses in differentiating epithelia</article-title>
. <source>Microbiol Mol Biol Rev.</source>
<volume>68</volume>
:<fpage>362</fpage>
–<lpage>372</lpage>
.<pub-id pub-id-type="pmid">15187189</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B47"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Marin</surname>
<given-names>A</given-names>
</name>
<name><surname>Bertranpetit</surname>
<given-names>J</given-names>
</name>
<name><surname>Oliver</surname>
<given-names>JL</given-names>
</name>
<name><surname>Medina</surname>
<given-names>JR</given-names>
</name>
</person-group>
<year>1989</year>
<article-title>Variation in G + C-content and codon choice: differences among synonymous codon groups in vertebrate genes</article-title>
. <source>Nucleic Acids Res.</source>
<volume>17</volume>
:<fpage>6181</fpage>
–<lpage>6189</lpage>
.<pub-id pub-id-type="pmid">2570402</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B48"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Martincorena</surname>
<given-names>I</given-names>
</name>
<etal></etal>
</person-group>
<year>2015</year>
<article-title>Tumor evolution. High burden and pervasive positive selection of somatic mutations in normal human skin</article-title>
. <source>Science</source>
<volume>348</volume>
:<fpage>880</fpage>
–<lpage>886</lpage>
.<pub-id pub-id-type="pmid">25999502</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B49"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Moriyama</surname>
<given-names>EN</given-names>
</name>
<name><surname>Powell</surname>
<given-names>JR</given-names>
</name>
</person-group>
<year>1997</year>
<article-title>Codon usage bias and tRNA abundance in <italic>Drosophila</italic>
</article-title>
. <source>J Mol Evol.</source>
<volume>45</volume>
:<fpage>514</fpage>
–<lpage>523</lpage>
.<pub-id pub-id-type="pmid">9342399</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B50"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mossadegh</surname>
<given-names>N</given-names>
</name>
<etal></etal>
</person-group>
<year>2004</year>
<article-title>Codon optimization of the human papillomavirus 11 (HPV 11) L1 gene leads to increased gene expression and formation of virus-like particles in mammalian epithelial cells</article-title>
. <source>Virology</source>
<volume>326</volume>
:<fpage>57</fpage>
–<lpage>66</lpage>
.<pub-id pub-id-type="pmid">15262495</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B51"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Munk</surname>
<given-names>C</given-names>
</name>
<name><surname>Willemsen</surname>
<given-names>A</given-names>
</name>
<name><surname>Bravo</surname>
<given-names>IG</given-names>
</name>
</person-group>
<year>2012</year>
<article-title>An ancient history of gene duplications, fusions and losses in the evolution of APOBEC3 mutators in mammals</article-title>
. <source>BMC Evol Biol.</source>
<volume>12</volume>
:<fpage>71</fpage>
.<pub-id pub-id-type="pmid">22640020</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B52"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nakamura</surname>
<given-names>Y</given-names>
</name>
<name><surname>Gojobori</surname>
<given-names>T</given-names>
</name>
<name><surname>Ikemura</surname>
<given-names>T</given-names>
</name>
</person-group>
<year>2000</year>
<article-title>Codon usage tabulated from international DNA sequence databases: status for the year 2000</article-title>
. <source>Nucleic Acids Res.</source>
<volume>28</volume>
:<fpage>292</fpage>
.<pub-id pub-id-type="pmid">10592250</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B53"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Narechania</surname>
<given-names>A</given-names>
</name>
<name><surname>Terai</surname>
<given-names>M</given-names>
</name>
<name><surname>Burk</surname>
<given-names>RD</given-names>
</name>
</person-group>
<year>2005</year>
<article-title>Overlapping reading frames in closely related human papillomaviruses result in modular rates of selection within E2</article-title>
. <source>J Gen Virol.</source>
<volume>86</volume>
:<fpage>1307</fpage>
–<lpage>1313</lpage>
.<pub-id pub-id-type="pmid">15831941</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B54"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nindl</surname>
<given-names>I</given-names>
</name>
<name><surname>Gottschling</surname>
<given-names>M</given-names>
</name>
<name><surname>Stockfleth</surname>
<given-names>E</given-names>
</name>
</person-group>
<year>2007</year>
<article-title>Human papillomaviruses and non-melanoma skin cancer: basic virology and clinical manifestations</article-title>
. <source>Dis Markers.</source>
<volume>23</volume>
:<fpage>247</fpage>
–<lpage>259</lpage>
.<pub-id pub-id-type="pmid">17627060</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B55"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nobre</surname>
<given-names>RJ</given-names>
</name>
<name><surname>de Almeida</surname>
<given-names>LP</given-names>
</name>
<name><surname>Martins</surname>
<given-names>TC</given-names>
</name>
</person-group>
<year>2008</year>
<article-title>Complete genotyping of mucosal human papillomavirus using a restriction fragment length polymorphism analysis and an original typing algorithm</article-title>
. <source>J Clin Virol.</source>
<volume>42</volume>
:<fpage>13</fpage>
–<lpage>21</lpage>
.<pub-id pub-id-type="pmid">18304866</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B56"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Park</surname>
<given-names>P</given-names>
</name>
<etal></etal>
</person-group>
<year>1994</year>
<article-title>The cellular DNA polymerase alpha-primase is required for papillomavirus DNA replication and associates with the viral E1 helicase</article-title>
. <source>Proc Natl Acad Sci U S A.</source>
<volume>91</volume>
:<fpage>8700</fpage>
–<lpage>8704</lpage>
.<pub-id pub-id-type="pmid">8078945</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B57"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Peh</surname>
<given-names>WL</given-names>
</name>
<etal></etal>
</person-group>
<year>2002</year>
<article-title>Life cycle heterogeneity in animal models of human papillomavirus-associated disease</article-title>
. <source>J Virol.</source>
<volume>76</volume>
:<fpage>10401</fpage>
–<lpage>10416</lpage>
.<pub-id pub-id-type="pmid">12239317</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B58"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Plotkin</surname>
<given-names>JB</given-names>
</name>
<name><surname>Kudla</surname>
<given-names>G</given-names>
</name>
</person-group>
<year>2011</year>
<article-title>Synonymous but not the same: the causes and consequences of codon bias</article-title>
. <source>Nat Rev Genet.</source>
<volume>12</volume>
:<fpage>32</fpage>
–<lpage>42</lpage>
.<pub-id pub-id-type="pmid">21102527</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B59"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Protic-Sabljic</surname>
<given-names>M</given-names>
</name>
<etal></etal>
</person-group>
<year>1986</year>
<article-title>UV light-induced cyclobutane pyrimidine dimers are mutagenic in mammalian cells</article-title>
. <source>Mol Cell Biol.</source>
<volume>6</volume>
:<fpage>3349</fpage>
–<lpage>3356</lpage>
.<pub-id pub-id-type="pmid">3540589</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B60"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rancurel</surname>
<given-names>C</given-names>
</name>
<name><surname>Khosravi</surname>
<given-names>M</given-names>
</name>
<name><surname>Dunker</surname>
<given-names>AK</given-names>
</name>
<name><surname>Romero</surname>
<given-names>PR</given-names>
</name>
<name><surname>Karlin</surname>
<given-names>D</given-names>
</name>
</person-group>
<year>2009</year>
<article-title>Overlapping genes produce proteins with unusual sequence properties and offer insight into de novo protein creation</article-title>
. <source>J Virol.</source>
<volume>83</volume>
:<fpage>10719</fpage>
–<lpage>10736</lpage>
.<pub-id pub-id-type="pmid">19640978</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B61"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rocha</surname>
<given-names>EP</given-names>
</name>
</person-group>
<year>2004</year>
<article-title>Codon usage bias from tRNA’s point of view: redundancy, specialization, and efficient decoding for translation optimization</article-title>
. <source>Genome Res.</source>
<volume>14</volume>
:<fpage>2279</fpage>
–<lpage>2286</lpage>
.<pub-id pub-id-type="pmid">15479947</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B62"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Roychoudhury</surname>
<given-names>S</given-names>
</name>
<name><surname>Mukherjee</surname>
<given-names>D</given-names>
</name>
</person-group>
<year>2009</year>
<article-title>A detailed comparative analysis on the overall codon usage pattern in herpesviruses</article-title>
. <source>Virus Res.</source>
<volume>148</volume>
:<fpage>31</fpage>
–<lpage>43</lpage>
.<pub-id pub-id-type="pmid">19969032</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B63"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Samorski</surname>
<given-names>R</given-names>
</name>
<name><surname>Gissmann</surname>
<given-names>L</given-names>
</name>
<name><surname>Osen</surname>
<given-names>W</given-names>
</name>
</person-group>
<year>2006</year>
<article-title>Codon optimized expression of HPV 16 E6 renders target cells susceptible to E6-specific CTL recognition</article-title>
. <source>Immunol Lett.</source>
<volume>107</volume>
:<fpage>41</fpage>
–<lpage>49</lpage>
.<pub-id pub-id-type="pmid">16949679</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B64"><mixed-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>J Mol Evol.</source>
<volume>62</volume>
:<fpage>551</fpage>
–<lpage>563</lpage>
.<pub-id pub-id-type="pmid">16557338</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B65"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sharp</surname>
<given-names>PM</given-names>
</name>
<name><surname>Averof</surname>
<given-names>M</given-names>
</name>
<name><surname>Lloyd</surname>
<given-names>AT</given-names>
</name>
<name><surname>Matassi</surname>
<given-names>G</given-names>
</name>
<name><surname>Peden</surname>
<given-names>JF</given-names>
</name>
</person-group>
<year>1995</year>
<article-title>DNA sequence evolution: the sounds of silence</article-title>
. <source>Philos Trans R Soc Lond B Biol Sci.</source>
<volume>349</volume>
:<fpage>241</fpage>
–<lpage>247</lpage>
.<pub-id pub-id-type="pmid">8577834</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B66"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sharp</surname>
<given-names>PM</given-names>
</name>
<name><surname>Li</surname>
<given-names>WH</given-names>
</name>
</person-group>
<year>1986</year>
<article-title>Codon usage in regulatory genes in <italic>Escherichia coli</italic>
does not reflect selection for ‘rare’ codons</article-title>
. <source>Nucleic Acids Res.</source>
<volume>14</volume>
:<fpage>7737</fpage>
–<lpage>7749</lpage>
.<pub-id pub-id-type="pmid">3534792</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B67"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sharp</surname>
<given-names>PM</given-names>
</name>
<name><surname>Li</surname>
<given-names>WH</given-names>
</name>
</person-group>
<year>1987</year>
<article-title>The codon Adaptation Index—a measure of directional synonymous codon usage bias, and its potential applications</article-title>
. <source>Nucleic Acids Res.</source>
<volume>15</volume>
:<fpage>1281</fpage>
–<lpage>1295</lpage>
.<pub-id pub-id-type="pmid">3547335</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B68"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Shields</surname>
<given-names>DC</given-names>
</name>
<name><surname>Sharp</surname>
<given-names>PM</given-names>
</name>
</person-group>
<year>1987</year>
<article-title>Synonymous codon usage in <italic>Bacillus subtilis</italic>
reflects both translational selection and mutational biases</article-title>
. <source>Nucleic Acids Res.</source>
<volume>15</volume>
:<fpage>8023</fpage>
–<lpage>8040</lpage>
.<pub-id pub-id-type="pmid">3118331</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B69"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sohn</surname>
<given-names>WJ</given-names>
</name>
<etal></etal>
</person-group>
<year>2007</year>
<article-title>Novel transcriptional regulation of the schlafen-2 gene in macrophages in response to TLR-triggered stimulation</article-title>
. <source>Mol Immunol.</source>
<volume>44</volume>
:<fpage>3273</fpage>
–<lpage>3282</lpage>
.<pub-id pub-id-type="pmid">17434208</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B70"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Stamatakis</surname>
<given-names>A</given-names>
</name>
</person-group>
<year>2014</year>
<article-title>RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies</article-title>
. <source>Bioinformatics</source>
<volume>30</volume>
:<fpage>1312</fpage>
–<lpage>1313</lpage>
.<pub-id pub-id-type="pmid">24451623</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B71"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Stenico</surname>
<given-names>M</given-names>
</name>
<name><surname>Lloyd</surname>
<given-names>AT</given-names>
</name>
<name><surname>Sharp</surname>
<given-names>PM</given-names>
</name>
</person-group>
<year>1994</year>
<article-title>Codon usage in <italic>Caenorhabditis elegans</italic>
: delineation of translational selection and mutational biases</article-title>
. <source>Nucleic Acids Res.</source>
<volume>22</volume>
:<fpage>2437</fpage>
–<lpage>2446</lpage>
.<pub-id pub-id-type="pmid">8041603</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B72"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sugiyama</surname>
<given-names>T</given-names>
</name>
<etal></etal>
</person-group>
<year>2005</year>
<article-title>CpG RNA: identification of novel single-stranded RNA that stimulates human CD14+CD11c+ monocytes</article-title>
. <source>J Immunol.</source>
<volume>174</volume>
:<fpage>2273</fpage>
–<lpage>2279</lpage>
.<pub-id pub-id-type="pmid">15699162</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B73"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Suzuki</surname>
<given-names>H</given-names>
</name>
<name><surname>Brown</surname>
<given-names>CJ</given-names>
</name>
<name><surname>Forney</surname>
<given-names>LJ</given-names>
</name>
<name><surname>Top</surname>
<given-names>EM</given-names>
</name>
</person-group>
<year>2008</year>
<article-title>Comparison of correspondence analysis methods for synonymous codon usage in bacteria</article-title>
. <source>DNA Res.</source>
<volume>15</volume>
:<fpage>357</fpage>
–<lpage>365</lpage>
.<pub-id pub-id-type="pmid">18940873</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B74"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Temiz</surname>
<given-names>NA</given-names>
</name>
<etal></etal>
</person-group>
<year>2015</year>
<article-title>The somatic autosomal mutation matrix in cancer genomes</article-title>
. <source>Hum Genet.</source>
<volume>134</volume>
(8):<fpage>851</fpage>
–<lpage>864</lpage>
.<pub-id pub-id-type="pmid">26001532</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B75"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tindle</surname>
<given-names>RW</given-names>
</name>
</person-group>
<year>2002</year>
<article-title>Immune evasion in human papillomavirus-associated cervical cancer</article-title>
. <source>Nat Rev Cancer.</source>
<volume>2</volume>
:<fpage>59</fpage>
–<lpage>65</lpage>
.<pub-id pub-id-type="pmid">11902586</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B76"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Turabelidze</surname>
<given-names>A</given-names>
</name>
<etal></etal>
</person-group>
<year>2014</year>
<article-title>Intrinsic differences between oral and skin keratinocytes</article-title>
. <source>PLoS One</source>
<volume>9</volume>
:<fpage>e101480</fpage>
.<pub-id pub-id-type="pmid">25198578</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B77"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Vartanian</surname>
<given-names>JP</given-names>
</name>
<name><surname>Guetard</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>233</lpage>
.<pub-id pub-id-type="pmid">18403710</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B78"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Vartanian</surname>
<given-names>JP</given-names>
</name>
<etal></etal>
</person-group>
<year>2010</year>
<article-title>Massive APOBEC3 editing of hepatitis B viral DNA in cirrhosis</article-title>
. <source>PLoS Pathog.</source>
<volume>6</volume>
:<fpage>e1000928</fpage>
.<pub-id pub-id-type="pmid">20523896</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B79"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhao</surname>
<given-names>KN</given-names>
</name>
<name><surname>Gu</surname>
<given-names>W</given-names>
</name>
<name><surname>Fang</surname>
<given-names>NX</given-names>
</name>
<name><surname>Saunders</surname>
<given-names>NA</given-names>
</name>
<name><surname>Frazer</surname>
<given-names>IH</given-names>
</name>
</person-group>
<year>2005</year>
<article-title>Gene codon composition determines differentiation-dependent expression of a viral capsid gene in keratinocytes in vitro and in vivo</article-title>
. <source>Mol Cell Biol.</source>
<volume>25</volume>
:<fpage>8643</fpage>
–<lpage>8655</lpage>
.<pub-id pub-id-type="pmid">16166644</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B80"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhao</surname>
<given-names>KN</given-names>
</name>
<name><surname>Liu</surname>
<given-names>WJ</given-names>
</name>
<name><surname>Frazer</surname>
<given-names>IH</given-names>
</name>
</person-group>
<year>2003</year>
<article-title>Codon usage bias and A+T content variation in human papillomavirus genomes</article-title>
. <source>Virus Res.</source>
<volume>98</volume>
:<fpage>95</fpage>
–<lpage>104</lpage>
.<pub-id pub-id-type="pmid">14659556</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B81"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zheng</surname>
<given-names>ZM</given-names>
</name>
<name><surname>Baker</surname>
<given-names>CC</given-names>
</name>
</person-group>
<year>2006</year>
<article-title>Papillomavirus genome structure, expression, and post-transcriptional regulation</article-title>
. <source>Front Biosci.</source>
<volume>11</volume>
:<fpage>2286</fpage>
–<lpage>2302</lpage>
.<pub-id pub-id-type="pmid">16720315</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B82"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhou</surname>
<given-names>J</given-names>
</name>
<name><surname>Liu</surname>
<given-names>WJ</given-names>
</name>
<name><surname>Peng</surname>
<given-names>SW</given-names>
</name>
<name><surname>Sun</surname>
<given-names>XY</given-names>
</name>
<name><surname>Frazer</surname>
<given-names>I</given-names>
</name>
</person-group>
<year>1999</year>
<article-title>Papillomavirus capsid protein expression level depends on the match between codon usage and tRNA availability</article-title>
. <source>J Virol.</source>
<volume>73</volume>
:<fpage>4972</fpage>
–<lpage>4982</lpage>
.<pub-id pub-id-type="pmid">10233959</pub-id>
</mixed-citation>
</ref>
<ref id="evv129-B83"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhou</surname>
<given-names>M</given-names>
</name>
<name><surname>Wang</surname>
<given-names>T</given-names>
</name>
<name><surname>Fu</surname>
<given-names>J</given-names>
</name>
<name><surname>Xiao</surname>
<given-names>G</given-names>
</name>
<name><surname>Liu</surname>
<given-names>Y</given-names>
</name>
</person-group>
<year>2015</year>
<article-title>Non-optimal codon usage influences protein structure in intrinsically disordered regions</article-title>
. <source>Mol Microbiol.</source>
<comment>Advance Access published June 25, 2015; doi: 10.1111/mmi.13079.</comment>
</mixed-citation>
</ref>
</ref-list>
</back>
</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 0000050 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 0000050 | 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é= |texte= }}
![]() | This area was generated with Dilib version V0.6.31. | ![]() |