Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Genetic Variation of Human Papillomavirus Type 16 in Individual Clinical Specimens Revealed by Deep Sequencing</title>
<author><name sortKey="Kukimoto, Iwao" sort="Kukimoto, Iwao" uniqKey="Kukimoto I" first="Iwao" last="Kukimoto">Iwao Kukimoto</name>
<affiliation><nlm:aff id="aff1"><addr-line>Pathogen Genomics Center, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Maehama, Tomohiko" sort="Maehama, Tomohiko" uniqKey="Maehama T" first="Tomohiko" last="Maehama">Tomohiko Maehama</name>
<affiliation><nlm:aff id="aff2"><addr-line>Department of Biochemistry and Cell Biology, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Sekizuka, Tsuyoshi" sort="Sekizuka, Tsuyoshi" uniqKey="Sekizuka T" first="Tsuyoshi" last="Sekizuka">Tsuyoshi Sekizuka</name>
<affiliation><nlm:aff id="aff1"><addr-line>Pathogen Genomics Center, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Ogasawara, Yumiko" sort="Ogasawara, Yumiko" uniqKey="Ogasawara Y" first="Yumiko" last="Ogasawara">Yumiko Ogasawara</name>
<affiliation><nlm:aff id="aff1"><addr-line>Pathogen Genomics Center, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Kondo, Kazunari" sort="Kondo, Kazunari" uniqKey="Kondo K" first="Kazunari" last="Kondo">Kazunari Kondo</name>
<affiliation><nlm:aff id="aff3"><addr-line>NTT Medical Center Tokyo, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Kusumoto Matsuo, Rika" sort="Kusumoto Matsuo, Rika" uniqKey="Kusumoto Matsuo R" first="Rika" last="Kusumoto-Matsuo">Rika Kusumoto-Matsuo</name>
<affiliation><nlm:aff id="aff1"><addr-line>Pathogen Genomics Center, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Mori, Seiichiro" sort="Mori, Seiichiro" uniqKey="Mori S" first="Seiichiro" last="Mori">Seiichiro Mori</name>
<affiliation><nlm:aff id="aff1"><addr-line>Pathogen Genomics Center, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Ishii, Yoshiyuki" sort="Ishii, Yoshiyuki" uniqKey="Ishii Y" first="Yoshiyuki" last="Ishii">Yoshiyuki Ishii</name>
<affiliation><nlm:aff id="aff1"><addr-line>Pathogen Genomics Center, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Takeuchi, Takamasa" sort="Takeuchi, Takamasa" uniqKey="Takeuchi T" first="Takamasa" last="Takeuchi">Takamasa Takeuchi</name>
<affiliation><nlm:aff id="aff1"><addr-line>Pathogen Genomics Center, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Yamaji, Toshiyuki" sort="Yamaji, Toshiyuki" uniqKey="Yamaji T" first="Toshiyuki" last="Yamaji">Toshiyuki Yamaji</name>
<affiliation><nlm:aff id="aff2"><addr-line>Department of Biochemistry and Cell Biology, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Takeuchi, Fumihiko" sort="Takeuchi, Fumihiko" uniqKey="Takeuchi F" first="Fumihiko" last="Takeuchi">Fumihiko Takeuchi</name>
<affiliation><nlm:aff id="aff1"><addr-line>Pathogen Genomics Center, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Hanada, Kentaro" sort="Hanada, Kentaro" uniqKey="Hanada K" first="Kentaro" last="Hanada">Kentaro Hanada</name>
<affiliation><nlm:aff id="aff2"><addr-line>Department of Biochemistry and Cell Biology, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Kuroda, Makoto" sort="Kuroda, Makoto" uniqKey="Kuroda M" first="Makoto" last="Kuroda">Makoto Kuroda</name>
<affiliation><nlm:aff id="aff1"><addr-line>Pathogen Genomics Center, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">24236186</idno>
<idno type="pmc">3827439</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3827439</idno>
<idno type="RBID">PMC:3827439</idno>
<idno type="doi">10.1371/journal.pone.0080583</idno>
<date when="2013">2013</date>
<idno type="wicri:Area/Pmc/Corpus">000131</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000131</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Genetic Variation of Human Papillomavirus Type 16 in Individual Clinical Specimens Revealed by Deep Sequencing</title>
<author><name sortKey="Kukimoto, Iwao" sort="Kukimoto, Iwao" uniqKey="Kukimoto I" first="Iwao" last="Kukimoto">Iwao Kukimoto</name>
<affiliation><nlm:aff id="aff1"><addr-line>Pathogen Genomics Center, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Maehama, Tomohiko" sort="Maehama, Tomohiko" uniqKey="Maehama T" first="Tomohiko" last="Maehama">Tomohiko Maehama</name>
<affiliation><nlm:aff id="aff2"><addr-line>Department of Biochemistry and Cell Biology, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Sekizuka, Tsuyoshi" sort="Sekizuka, Tsuyoshi" uniqKey="Sekizuka T" first="Tsuyoshi" last="Sekizuka">Tsuyoshi Sekizuka</name>
<affiliation><nlm:aff id="aff1"><addr-line>Pathogen Genomics Center, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Ogasawara, Yumiko" sort="Ogasawara, Yumiko" uniqKey="Ogasawara Y" first="Yumiko" last="Ogasawara">Yumiko Ogasawara</name>
<affiliation><nlm:aff id="aff1"><addr-line>Pathogen Genomics Center, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Kondo, Kazunari" sort="Kondo, Kazunari" uniqKey="Kondo K" first="Kazunari" last="Kondo">Kazunari Kondo</name>
<affiliation><nlm:aff id="aff3"><addr-line>NTT Medical Center Tokyo, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Kusumoto Matsuo, Rika" sort="Kusumoto Matsuo, Rika" uniqKey="Kusumoto Matsuo R" first="Rika" last="Kusumoto-Matsuo">Rika Kusumoto-Matsuo</name>
<affiliation><nlm:aff id="aff1"><addr-line>Pathogen Genomics Center, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Mori, Seiichiro" sort="Mori, Seiichiro" uniqKey="Mori S" first="Seiichiro" last="Mori">Seiichiro Mori</name>
<affiliation><nlm:aff id="aff1"><addr-line>Pathogen Genomics Center, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Ishii, Yoshiyuki" sort="Ishii, Yoshiyuki" uniqKey="Ishii Y" first="Yoshiyuki" last="Ishii">Yoshiyuki Ishii</name>
<affiliation><nlm:aff id="aff1"><addr-line>Pathogen Genomics Center, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Takeuchi, Takamasa" sort="Takeuchi, Takamasa" uniqKey="Takeuchi T" first="Takamasa" last="Takeuchi">Takamasa Takeuchi</name>
<affiliation><nlm:aff id="aff1"><addr-line>Pathogen Genomics Center, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Yamaji, Toshiyuki" sort="Yamaji, Toshiyuki" uniqKey="Yamaji T" first="Toshiyuki" last="Yamaji">Toshiyuki Yamaji</name>
<affiliation><nlm:aff id="aff2"><addr-line>Department of Biochemistry and Cell Biology, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Takeuchi, Fumihiko" sort="Takeuchi, Fumihiko" uniqKey="Takeuchi F" first="Fumihiko" last="Takeuchi">Fumihiko Takeuchi</name>
<affiliation><nlm:aff id="aff1"><addr-line>Pathogen Genomics Center, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Hanada, Kentaro" sort="Hanada, Kentaro" uniqKey="Hanada K" first="Kentaro" last="Hanada">Kentaro Hanada</name>
<affiliation><nlm:aff id="aff2"><addr-line>Department of Biochemistry and Cell Biology, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Kuroda, Makoto" sort="Kuroda, Makoto" uniqKey="Kuroda M" first="Makoto" last="Kuroda">Makoto Kuroda</name>
<affiliation><nlm:aff id="aff1"><addr-line>Pathogen Genomics Center, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">PLoS ONE</title>
<idno type="eISSN">1932-6203</idno>
<imprint><date when="2013">2013</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><p>Viral genetic diversity within infected cells or tissues, called viral quasispecies, has been mostly studied for RNA viruses, but has also been described among DNA viruses, including human papillomavirus type 16 (HPV16) present in cervical precancerous lesions. However, the extent of HPV genetic variation in cervical specimens, and its involvement in HPV-induced carcinogenesis, remains unclear. Here, we employ deep sequencing to comprehensively analyze genetic variation in the HPV16 genome isolated from individual clinical specimens. Through overlapping full-circle PCR, approximately 8-kb DNA fragments covering the whole HPV16 genome were amplified from HPV16-positive cervical exfoliated cells collected from patients with either low-grade squamous intraepithelial lesion (LSIL) or invasive cervical cancer (ICC). Deep sequencing of the amplified HPV16 DNA enabled <italic>de novo</italic>
assembly of the full-length HPV16 genome sequence for each of 7 specimens (5 LSIL and 2 ICC samples). Subsequent alignment of read sequences to the assembled HPV16 sequence revealed that 2 LSILs and 1 ICC contained nucleotide variations within E6, E1 and the non-coding region between E5 and L2 with mutation frequencies of 0.60% to 5.42%. In transient replication assays, a novel E1 mutant found in ICC, E1 Q381E, showed reduced ability to support HPV16 origin-dependent replication. In addition, partially deleted E2 genes were detected in 1 LSIL sample in a mixed state with the intact E2 gene. Thus, the methods used in this study provide a fundamental framework for investigating the influence of HPV somatic genetic variation on cervical carcinogenesis. </p>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Zur Hausen, H" uniqKey="Zur Hausen H">H zur Hausen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="De Sanjose, S" uniqKey="De Sanjose S">S de Sanjose</name>
</author>
<author><name sortKey="Quint, Wg" uniqKey="Quint W">WG Quint</name>
</author>
<author><name sortKey="Alemany, L" uniqKey="Alemany L">L Alemany</name>
</author>
<author><name sortKey="Geraets, Dt" uniqKey="Geraets D">DT Geraets</name>
</author>
<author><name sortKey="Klaustermeier, Je" uniqKey="Klaustermeier J">JE Klaustermeier</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kajitani, N" uniqKey="Kajitani N">N Kajitani</name>
</author>
<author><name sortKey="Satsuka, A" uniqKey="Satsuka A">A Satsuka</name>
</author>
<author><name sortKey="Kawate, A" uniqKey="Kawate A">A Kawate</name>
</author>
<author><name sortKey="Sakai, H" uniqKey="Sakai H">H Sakai</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pett, M" uniqKey="Pett M">M Pett</name>
</author>
<author><name sortKey="Coleman, N" uniqKey="Coleman N">N Coleman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Xu, B" uniqKey="Xu B">B Xu</name>
</author>
<author><name sortKey="Chotewutmontri, S" uniqKey="Chotewutmontri S">S Chotewutmontri</name>
</author>
<author><name sortKey="Wolf, S" uniqKey="Wolf S">S Wolf</name>
</author>
<author><name sortKey="Klos, U" uniqKey="Klos U">U Klos</name>
</author>
<author><name sortKey="Schmitz, M" uniqKey="Schmitz M">M Schmitz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hebner, Cm" uniqKey="Hebner C">CM Hebner</name>
</author>
<author><name sortKey="Laimins, La" uniqKey="Laimins L">LA Laimins</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bodily, J" uniqKey="Bodily J">J Bodily</name>
</author>
<author><name sortKey="Laimins, La" uniqKey="Laimins L">LA Laimins</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Beerenwinkel, N" uniqKey="Beerenwinkel N">N Beerenwinkel</name>
</author>
<author><name sortKey="Gunthard, Hf" uniqKey="Gunthard H">HF Gunthard</name>
</author>
<author><name sortKey="Roth, V" uniqKey="Roth V">V Roth</name>
</author>
<author><name sortKey="Metzner, Kj" uniqKey="Metzner K">KJ Metzner</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ekstrom, J" uniqKey="Ekstrom J">J Ekström</name>
</author>
<author><name sortKey="Bzhalava, D" uniqKey="Bzhalava D">D Bzhalava</name>
</author>
<author><name sortKey="Svenback, D" uniqKey="Svenback D">D Svenback</name>
</author>
<author><name sortKey="Forslund, O" uniqKey="Forslund O">O Forslund</name>
</author>
<author><name sortKey="Dillner, J" uniqKey="Dillner J">J Dillner</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Meiring, Tl" uniqKey="Meiring T">TL Meiring</name>
</author>
<author><name sortKey="Salimo, At" uniqKey="Salimo A">AT Salimo</name>
</author>
<author><name sortKey="Coetzee, B" uniqKey="Coetzee B">B Coetzee</name>
</author>
<author><name sortKey="Maree, Hj" uniqKey="Maree H">HJ Maree</name>
</author>
<author><name sortKey="Moodley, J" uniqKey="Moodley J">J Moodley</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Foulongne, V" uniqKey="Foulongne V">V Foulongne</name>
</author>
<author><name sortKey="Sauvage, V" uniqKey="Sauvage V">V Sauvage</name>
</author>
<author><name sortKey="Hebert, C" uniqKey="Hebert C">C Hebert</name>
</author>
<author><name sortKey="Dereure, O" uniqKey="Dereure O">O Dereure</name>
</author>
<author><name sortKey="Cheval, J" uniqKey="Cheval J">J Cheval</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lauring, As" uniqKey="Lauring A">AS Lauring</name>
</author>
<author><name sortKey="Andino, R" uniqKey="Andino R">R Andino</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Domingo, E" uniqKey="Domingo E">E Domingo</name>
</author>
<author><name sortKey="Sheldon, J" uniqKey="Sheldon J">J Sheldon</name>
</author>
<author><name sortKey="Perales, C" uniqKey="Perales C">C Perales</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Van Velzen, M" uniqKey="Van Velzen M">M van Velzen</name>
</author>
<author><name sortKey="Van Loenen, Fb" uniqKey="Van Loenen F">FB van Loenen</name>
</author>
<author><name sortKey="Meesters, Rj" uniqKey="Meesters R">RJ Meesters</name>
</author>
<author><name sortKey="De Graaf, M" uniqKey="De Graaf M">M de Graaf</name>
</author>
<author><name sortKey="Remeijer, L" uniqKey="Remeijer L">L Remeijer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kwok, H" uniqKey="Kwok H">H Kwok</name>
</author>
<author><name sortKey="Tong, Ah" uniqKey="Tong A">AH Tong</name>
</author>
<author><name sortKey="Lin, Ch" uniqKey="Lin C">CH Lin</name>
</author>
<author><name sortKey="Lok, S" uniqKey="Lok S">S Lok</name>
</author>
<author><name sortKey="Farrell, Pj" uniqKey="Farrell P">PJ Farrell</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 Guétard</name>
</author>
<author><name sortKey="Henry, M" uniqKey="Henry M">M Henry</name>
</author>
<author><name sortKey="Wain Hobson, S" uniqKey="Wain Hobson S">S Wain-Hobson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jeon, S" uniqKey="Jeon S">S Jeon</name>
</author>
<author><name sortKey="Allen Hoffmann, Bl" uniqKey="Allen Hoffmann B">BL Allen-Hoffmann</name>
</author>
<author><name sortKey="Lambert, Pf" uniqKey="Lambert P">PF Lambert</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Stanley, Ma" uniqKey="Stanley M">MA Stanley</name>
</author>
<author><name sortKey="Browne, Hm" uniqKey="Browne H">HM Browne</name>
</author>
<author><name sortKey="Appleby, M" uniqKey="Appleby M">M Appleby</name>
</author>
<author><name sortKey="Minson, Ac" uniqKey="Minson A">AC Minson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Arad, U" uniqKey="Arad U">U Arad</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Estrade, C" uniqKey="Estrade C">C Estrade</name>
</author>
<author><name sortKey="Menoud, Pa" uniqKey="Menoud P">PA Menoud</name>
</author>
<author><name sortKey="Nardelli Haefliger, D" uniqKey="Nardelli Haefliger D">D Nardelli-Haefliger</name>
</author>
<author><name sortKey="Sahli, R" uniqKey="Sahli R">R Sahli</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kanda, T" uniqKey="Kanda T">T Kanda</name>
</author>
<author><name sortKey="Watanabe, S" uniqKey="Watanabe S">S Watanabe</name>
</author>
<author><name sortKey="Yoshiike, K" uniqKey="Yoshiike K">K Yoshiike</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Simpson, Jt" uniqKey="Simpson J">JT Simpson</name>
</author>
<author><name sortKey="Wong, K" uniqKey="Wong K">K Wong</name>
</author>
<author><name sortKey="Jackman, Sd" uniqKey="Jackman S">SD Jackman</name>
</author>
<author><name sortKey="Schein, Je" uniqKey="Schein J">JE Schein</name>
</author>
<author><name sortKey="Jones, Sj" uniqKey="Jones S">SJ Jones</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gordon, D" uniqKey="Gordon D">D Gordon</name>
</author>
<author><name sortKey="Abajian, C" uniqKey="Abajian C">C Abajian</name>
</author>
<author><name sortKey="Green, P" uniqKey="Green P">P Green</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Li, H" uniqKey="Li H">H Li</name>
</author>
<author><name sortKey="Durbin, R" uniqKey="Durbin R">R Durbin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Futschik, A" uniqKey="Futschik A">A Futschik</name>
</author>
<author><name sortKey="Schlotterer, C" uniqKey="Schlotterer C">C Schlötterer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Thorvaldsd Ttir, H" uniqKey="Thorvaldsd Ttir H">H Thorvaldsdóttir</name>
</author>
<author><name sortKey="Robinson, Jt" uniqKey="Robinson J">JT Robinson</name>
</author>
<author><name sortKey="Mesirov, Jp" uniqKey="Mesirov J">JP Mesirov</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Milne, I" uniqKey="Milne I">I Milne</name>
</author>
<author><name sortKey="Bayer, M" uniqKey="Bayer M">M Bayer</name>
</author>
<author><name sortKey="Cardle, L" uniqKey="Cardle L">L Cardle</name>
</author>
<author><name sortKey="Shaw, P" uniqKey="Shaw P">P Shaw</name>
</author>
<author><name sortKey="Stephen, G" uniqKey="Stephen G">G Stephen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Li, H" uniqKey="Li H">H Li</name>
</author>
<author><name sortKey="Handsaker, B" uniqKey="Handsaker B">B Handsaker</name>
</author>
<author><name sortKey="Wysoker, A" uniqKey="Wysoker A">A Wysoker</name>
</author>
<author><name sortKey="Fennell, T" uniqKey="Fennell T">T Fennell</name>
</author>
<author><name sortKey="Ruan, J" uniqKey="Ruan J">J Ruan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Koboldt, Dc" uniqKey="Koboldt D">DC Koboldt</name>
</author>
<author><name sortKey="Chen, K" uniqKey="Chen K">K Chen</name>
</author>
<author><name sortKey="Wylie, T" uniqKey="Wylie T">T Wylie</name>
</author>
<author><name sortKey="Larson, De" uniqKey="Larson D">DE Larson</name>
</author>
<author><name sortKey="Mclellan, Md" uniqKey="Mclellan M">MD McLellan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Fradet Turcotte, A" uniqKey="Fradet Turcotte A">A Fradet-Turcotte</name>
</author>
<author><name sortKey="Morin, G" uniqKey="Morin G">G Morin</name>
</author>
<author><name sortKey="Lehoux, M" uniqKey="Lehoux M">M Lehoux</name>
</author>
<author><name sortKey="Bullock, Pa" uniqKey="Bullock P">PA Bullock</name>
</author>
<author><name sortKey="Archambault, J" uniqKey="Archambault J">J Archambault</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Yee, C" uniqKey="Yee C">C Yee</name>
</author>
<author><name sortKey="Krishnan Hewlett, I" uniqKey="Krishnan Hewlett I">I Krishnan-Hewlett</name>
</author>
<author><name sortKey="Baker, Cc" uniqKey="Baker C">CC Baker</name>
</author>
<author><name sortKey="Schlegel, R" uniqKey="Schlegel R">R Schlegel</name>
</author>
<author><name sortKey="Howley, Pm" uniqKey="Howley P">PM Howley</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Meissner, Jd" uniqKey="Meissner J">JD Meissner</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cornet, I" uniqKey="Cornet I">I Cornet</name>
</author>
<author><name sortKey="Gheit, T" uniqKey="Gheit T">T Gheit</name>
</author>
<author><name sortKey="Franceschi, S" uniqKey="Franceschi S">S Franceschi</name>
</author>
<author><name sortKey="Vignat, J" uniqKey="Vignat J">J Vignat</name>
</author>
<author><name sortKey="Burk, Rd" uniqKey="Burk R">RD Burk</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lin, Qq" uniqKey="Lin Q">QQ Lin</name>
</author>
<author><name sortKey="Yu, Sz" uniqKey="Yu S">SZ Yu</name>
</author>
<author><name sortKey="Qu, W" uniqKey="Qu W">W Qu</name>
</author>
<author><name sortKey="Cruz, Y" uniqKey="Cruz Y">Y Cruz</name>
</author>
<author><name sortKey="Burk, Rd" uniqKey="Burk R">RD Burk</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Miura, S" uniqKey="Miura S">S Miura</name>
</author>
<author><name sortKey="Matsumoto, K" uniqKey="Matsumoto K">K Matsumoto</name>
</author>
<author><name sortKey="Oki, A" uniqKey="Oki A">A Oki</name>
</author>
<author><name sortKey="Satoh, T" uniqKey="Satoh T">T Satoh</name>
</author>
<author><name sortKey="Tsunoda, H" uniqKey="Tsunoda H">H Tsunoda</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kondo, K" uniqKey="Kondo K">K Kondo</name>
</author>
<author><name sortKey="Uenoyama, A" uniqKey="Uenoyama A">A Uenoyama</name>
</author>
<author><name sortKey="Kitagawa, R" uniqKey="Kitagawa R">R Kitagawa</name>
</author>
<author><name sortKey="Tsunoda, H" uniqKey="Tsunoda H">H Tsunoda</name>
</author>
<author><name sortKey="Kusumoto Matsuo, R" uniqKey="Kusumoto Matsuo R">R Kusumoto-Matsuo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Titolo, S" uniqKey="Titolo S">S Titolo</name>
</author>
<author><name sortKey="Pelletier, A" uniqKey="Pelletier A">A Pelletier</name>
</author>
<author><name sortKey="Pulichino, Am" uniqKey="Pulichino A">AM Pulichino</name>
</author>
<author><name sortKey="Brault, K" uniqKey="Brault K">K Brault</name>
</author>
<author><name sortKey="Wardrop, E" uniqKey="Wardrop E">E Wardrop</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Melendy, T" uniqKey="Melendy T">T Melendy</name>
</author>
<author><name sortKey="Sedman, J" uniqKey="Sedman J">J Sedman</name>
</author>
<author><name sortKey="Stenlund, A" uniqKey="Stenlund A">A Stenlund</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Smith, B" uniqKey="Smith B">B Smith</name>
</author>
<author><name sortKey="Chen, Z" uniqKey="Chen Z">Z Chen</name>
</author>
<author><name sortKey="Reimers, L" uniqKey="Reimers L">L Reimers</name>
</author>
<author><name sortKey="Van Doorslaer, K" uniqKey="Van Doorslaer K">K van Doorslaer</name>
</author>
<author><name sortKey="Schiffman, M" uniqKey="Schiffman M">M Schiffman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jackson, Sp" uniqKey="Jackson S">SP Jackson</name>
</author>
<author><name sortKey="Bartek, J" uniqKey="Bartek J">J Bartek</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Amador Molina, A" uniqKey="Amador Molina A">A Amador-Molina</name>
</author>
<author><name sortKey="Gonzalez Montoya, Jl" uniqKey="Gonzalez Montoya J">JL González-Montoya</name>
</author>
<author><name sortKey="Garcia Carranca, A" uniqKey="Garcia Carranca A">A García-Carrancá</name>
</author>
<author><name sortKey="Mohar, A" uniqKey="Mohar A">A Mohar</name>
</author>
<author><name sortKey="Lizano, M" uniqKey="Lizano M">M Lizano</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Fu, L" uniqKey="Fu L">L Fu</name>
</author>
<author><name sortKey="Terai, M" uniqKey="Terai M">M Terai</name>
</author>
<author><name sortKey="Matsukura, T" uniqKey="Matsukura T">T Matsukura</name>
</author>
<author><name sortKey="Herrero, R" uniqKey="Herrero R">R Herrero</name>
</author>
<author><name sortKey="Burk, Rd" uniqKey="Burk R">RD Burk</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mcbride, Aa" uniqKey="Mcbride A">AA McBride</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zou, N" uniqKey="Zou N">N Zou</name>
</author>
<author><name sortKey="Lin, By" uniqKey="Lin B">BY Lin</name>
</author>
<author><name sortKey="Duan, F" uniqKey="Duan F">F Duan</name>
</author>
<author><name sortKey="Lee, Ky" uniqKey="Lee K">KY Lee</name>
</author>
<author><name sortKey="Jin, G" uniqKey="Jin G">G Jin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Arias Pulido, H" uniqKey="Arias Pulido H">H Arias-Pulido</name>
</author>
<author><name sortKey="Peyton, Cl" uniqKey="Peyton C">CL Peyton</name>
</author>
<author><name sortKey="Joste, Ne" uniqKey="Joste N">NE Joste</name>
</author>
<author><name sortKey="Vargas, H" uniqKey="Vargas H">H Vargas</name>
</author>
<author><name sortKey="Wheeler, Cm" uniqKey="Wheeler C">CM Wheeler</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Matsumoto, K" uniqKey="Matsumoto K">K Matsumoto</name>
</author>
<author><name sortKey="Yoshikawa, H" uniqKey="Yoshikawa H">H Yoshikawa</name>
</author>
<author><name sortKey="Nakagawa, S" uniqKey="Nakagawa S">S Nakagawa</name>
</author>
<author><name sortKey="Tang, X" uniqKey="Tang X">X Tang</name>
</author>
<author><name sortKey="Yasugi, T" uniqKey="Yasugi T">T Yasugi</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pande, S" uniqKey="Pande S">S Pande</name>
</author>
<author><name sortKey="Jain, N" uniqKey="Jain N">N Jain</name>
</author>
<author><name sortKey="Prusty, Bk" uniqKey="Prusty B">BK Prusty</name>
</author>
<author><name sortKey="Bhambhani, S" uniqKey="Bhambhani S">S Bhambhani</name>
</author>
<author><name sortKey="Gupta, S" uniqKey="Gupta S">S Gupta</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sichero, L" uniqKey="Sichero L">L Sichero</name>
</author>
<author><name sortKey="Sobrinho, Js" uniqKey="Sobrinho J">JS Sobrinho</name>
</author>
<author><name sortKey="Villa, Ll" uniqKey="Villa L">LL Villa</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cornet, I" uniqKey="Cornet I">I Cornet</name>
</author>
<author><name sortKey="Gheit, T" uniqKey="Gheit T">T Gheit</name>
</author>
<author><name sortKey="Iannacone, Mr" uniqKey="Iannacone M">MR Iannacone</name>
</author>
<author><name sortKey="Vignat, J" uniqKey="Vignat J">J Vignat</name>
</author>
<author><name sortKey="Sylla, Bs" uniqKey="Sylla B">BS Sylla</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Burk, Rd" uniqKey="Burk R">RD Burk</name>
</author>
<author><name sortKey="Harari, A" uniqKey="Harari A">A Harari</name>
</author>
<author><name sortKey="Chen, Z" uniqKey="Chen Z">Z Chen</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">PLoS One</journal-id>
<journal-id journal-id-type="iso-abbrev">PLoS ONE</journal-id>
<journal-id journal-id-type="publisher-id">plos</journal-id>
<journal-id journal-id-type="pmc">plosone</journal-id>
<journal-title-group><journal-title>PLoS ONE</journal-title>
</journal-title-group>
<issn pub-type="epub">1932-6203</issn>
<publisher><publisher-name>Public Library of Science</publisher-name>
<publisher-loc>San Francisco, USA</publisher-loc>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">24236186</article-id>
<article-id pub-id-type="pmc">3827439</article-id>
<article-id pub-id-type="publisher-id">PONE-D-13-21385</article-id>
<article-id pub-id-type="doi">10.1371/journal.pone.0080583</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject>
</subj-group>
</article-categories>
<title-group><article-title>Genetic Variation of Human Papillomavirus Type 16 in Individual Clinical Specimens Revealed by Deep Sequencing</article-title>
<alt-title alt-title-type="running-head">Variation Analysis of HPV16 by Deep Sequencing</alt-title>
</title-group>
<contrib-group><contrib contrib-type="author" equal-contrib="yes"><name><surname>Kukimoto</surname>
<given-names>Iwao</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
<xref ref-type="corresp" rid="cor1"><sup>*</sup>
</xref>
</contrib>
<contrib contrib-type="author" equal-contrib="yes"><name><surname>Maehama</surname>
<given-names>Tomohiko</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Sekizuka</surname>
<given-names>Tsuyoshi</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Ogasawara</surname>
<given-names>Yumiko</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Kondo</surname>
<given-names>Kazunari</given-names>
</name>
<xref ref-type="aff" rid="aff3"><sup>3</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Kusumoto-Matsuo</surname>
<given-names>Rika</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Mori</surname>
<given-names>Seiichiro</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Ishii</surname>
<given-names>Yoshiyuki</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Takeuchi</surname>
<given-names>Takamasa</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Yamaji</surname>
<given-names>Toshiyuki</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Takeuchi</surname>
<given-names>Fumihiko</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Hanada</surname>
<given-names>Kentaro</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Kuroda</surname>
<given-names>Makoto</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1"><label>1</label>
<addr-line>Pathogen Genomics Center, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</aff>
<aff id="aff2"><label>2</label>
<addr-line>Department of Biochemistry and Cell Biology, National Institute of Infectious Diseases, Tokyo, Japan</addr-line>
</aff>
<aff id="aff3"><label>3</label>
<addr-line>NTT Medical Center Tokyo, Tokyo, Japan</addr-line>
</aff>
<contrib-group><contrib contrib-type="editor"><name><surname>Burk</surname>
<given-names>Robert D.</given-names>
</name>
<role>Editor</role>
<xref ref-type="aff" rid="edit1"></xref>
</contrib>
</contrib-group>
<aff id="edit1"><addr-line>Albert Einstein College of Medicine, United States of America</addr-line>
</aff>
<author-notes><corresp id="cor1">* E-mail: <email>ikuki@niid.go.jp</email>
</corresp>
<fn fn-type="conflict"><p><bold>Competing Interests: </bold>
The authors have declared that no competing interests exist.</p>
</fn>
<fn fn-type="con"><p>Conceived and designed the experiments: IK TM. Performed the experiments: IK TM TS YO RK SM YI TT TY. Analyzed the data: IK TM TS FT KH MK. Contributed reagents/materials/analysis tools: TS KK. Wrote the manuscript: IK TM. </p>
</fn>
</author-notes>
<pub-date pub-type="collection"><year>2013</year>
</pub-date>
<pub-date pub-type="epub"><day>13</day>
<month>11</month>
<year>2013</year>
</pub-date>
<volume>8</volume>
<issue>11</issue>
<elocation-id>e80583</elocation-id>
<history><date date-type="received"><day>16</day>
<month>5</month>
<year>2013</year>
</date>
<date date-type="accepted"><day>4</day>
<month>10</month>
<year>2013</year>
</date>
</history>
<permissions><copyright-year>2013</copyright-year>
<copyright-holder>Kukimoto et al</copyright-holder>
<license><license-p>This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</license-p>
</license>
</permissions>
<abstract><p>Viral genetic diversity within infected cells or tissues, called viral quasispecies, has been mostly studied for RNA viruses, but has also been described among DNA viruses, including human papillomavirus type 16 (HPV16) present in cervical precancerous lesions. However, the extent of HPV genetic variation in cervical specimens, and its involvement in HPV-induced carcinogenesis, remains unclear. Here, we employ deep sequencing to comprehensively analyze genetic variation in the HPV16 genome isolated from individual clinical specimens. Through overlapping full-circle PCR, approximately 8-kb DNA fragments covering the whole HPV16 genome were amplified from HPV16-positive cervical exfoliated cells collected from patients with either low-grade squamous intraepithelial lesion (LSIL) or invasive cervical cancer (ICC). Deep sequencing of the amplified HPV16 DNA enabled <italic>de novo</italic>
assembly of the full-length HPV16 genome sequence for each of 7 specimens (5 LSIL and 2 ICC samples). Subsequent alignment of read sequences to the assembled HPV16 sequence revealed that 2 LSILs and 1 ICC contained nucleotide variations within E6, E1 and the non-coding region between E5 and L2 with mutation frequencies of 0.60% to 5.42%. In transient replication assays, a novel E1 mutant found in ICC, E1 Q381E, showed reduced ability to support HPV16 origin-dependent replication. In addition, partially deleted E2 genes were detected in 1 LSIL sample in a mixed state with the intact E2 gene. Thus, the methods used in this study provide a fundamental framework for investigating the influence of HPV somatic genetic variation on cervical carcinogenesis. </p>
</abstract>
<funding-group><funding-statement>This work was supported by Grants-in-Aid for Scientific Research from the Japanese Ministry of Education, Culture, Sports, and Science; by Grants-in-Aid for Scientific Research from the Ministry of Health, Labor, Welfare, and Technology; and by a grant from the Mitsubishi Foundation. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</funding-statement>
</funding-group>
</article-meta>
</front>
<body><sec id="s1"><title>Introduction</title>
<p>Human papillomaviruses (HPVs) are small DNA viruses having a circular double-stranded DNA genome of approximately 8-kb, some of which can induce benign and malignant hyper-proliferative lesions in the skin or mucosa [<xref ref-type="bibr" rid="B1">1</xref>
]. A subset of mucosal HPVs, referred to as high-risk HPVs, are causally related to the development of cervical cancer, among which HPV type 16 (HPV16) accounts for at least 50% of cervical cancer cases worldwide [<xref ref-type="bibr" rid="B2">2</xref>
]. Mucosal HPVs infect basal cells in the cervix and establish their genomes as extrachromosomal episomes. When infected cells migrate from the basal layer and begin the process of differentiation, viral genome amplification is initiated, followed by the induction of viral capsid expression and the assembly of infectious virions in the upper layers of epithelia [<xref ref-type="bibr" rid="B3">3</xref>
]. This differentiation-dependent propagation of HPVs can manifest clinically as low-grade squamous intraepithelial lesion (LSIL). In contrast to the episomal state in their normal life-cycle, DNAs of high-risk HPVs are frequently found integrated into the host genome in invasive cervical cancer (ICC) [<xref ref-type="bibr" rid="B4">4</xref>
,<xref ref-type="bibr" rid="B5">5</xref>
]. Viral DNA integration leads to enhanced expression of two viral transforming genes, E6 and E7, which play critical roles in carcinogenesis [<xref ref-type="bibr" rid="B6">6</xref>
,<xref ref-type="bibr" rid="B7">7</xref>
]. </p>
<p>In recent years next generation sequencing technologies have been increasingly incorporated into viral genomic research with the aim of comprehensively defining viral genomic sequences in clinical samples [<xref ref-type="bibr" rid="B8">8</xref>
]. HPV genomic research is no exception; the application of ultra high-throughput sequencing has been applied in the identification of known and unknown HPV genotypes in both cervical and cutaneous lesions [<xref ref-type="bibr" rid="B9">9</xref>
-<xref ref-type="bibr" rid="B11">11</xref>
]. In these studies, amplification of full-length HPV genome sequences from clinical specimens has been performed using unbiased amplification techniques such as rolling circle amplification or multiple displacement amplification with phi29 DNA polymerase, followed by deep sequencing analyses or next generation sequencing (NGS) [<xref ref-type="bibr" rid="B10">10</xref>
,<xref ref-type="bibr" rid="B11">11</xref>
]. </p>
<p>In addition, the genetic diversity of viral quasispecies in infected cells or tissues, which is a typical characteristic of RNA viruses [<xref ref-type="bibr" rid="B12">12</xref>
,<xref ref-type="bibr" rid="B13">13</xref>
], has been extensively examined by employing deep sequencing techniques. RNA viral populations consist of ‘mutant clouds’, rather than genomes with the same nucleotide sequence, due to the fact that they utilize low-fidelity viral RNA polymerases for their genome replication. Viral quasispecies play major roles in many of the biological activities of RNA viruses. For example, their capacity to change cell tropism or host range or to overcome internal and external selective constraints, such as immune responses and antiviral agents, depends on the presence of variants in mutant clouds. Viral genomic diversity in individual clinical specimens has also been described among DNA viruses, including herpes simplex virus type 1 [<xref ref-type="bibr" rid="B14">14</xref>
] and Epstein-Barr virus (EBV) [<xref ref-type="bibr" rid="B15">15</xref>
]. Intriguingly, the detection of hypermutated HPV1a and HPV16 genome sequences in plantar warts and precancerous cervical biopsies, respectively, has also been reported [<xref ref-type="bibr" rid="B16">16</xref>
]. However, the mutational landscape in the full-length HPV genome present in a clinical specimen and its physiological significance in viral carcinogenesis have not yet been investigated in a comprehensive manner. </p>
<p>In this study, we aimed to utilize NGS to analyze genetic variation in the entire HPV genome isolated from varying individual clinical samples. To this end, we have used a high-fidelity long-range PCR enzyme to amplify full-length HPV16 genomes from cultured cells and clinical samples. We report that PCR primers targeting a short region of the HPV16 genome yielded amplicons covering the entire viral genome sequence in a single reaction. By analyzing the amplified HPV16 DNA on the Illumina next generation sequencer, we have established a series of procedures for reconstructing the full-length HPV16 genome sequence and extracting genetic variations, such as nucleotide substitutions and insertions/deletions, from millions of short read sequences. This work could form the basis for a comprehensive analysis of HPV genetic variation present in individual clinical specimens. </p>
</sec>
<sec sec-type="materials|methods" id="s2"><title>Materials and Methods</title>
<sec id="s2.1"><title>Ethics Statement</title>
<p>Human cervical exfoliated cells were collected from patients at NTT Medical Center Tokyo with written informed consent. The study adhered to the declaration of Helsinki and was approved by the Ethics Committee at NTT Medical Center Tokyo. </p>
</sec>
<sec id="s2.2"><title>Cell culture</title>
<p>W12 cell clones 20850 (5 to 50 copies of HPV16 episomes per cell) and 20863 (500 to 1,000 copies of HPV16 episomes per cell) were cultured in the presence of mitomycin C-treated 3T3m feeder cells as described previously [<xref ref-type="bibr" rid="B17">17</xref>
,<xref ref-type="bibr" rid="B18">18</xref>
]. CaSki and human embryonic kidney 293 (HEK293) cells were cultured at 37°C with 5% CO<sub>2</sub>
in DMEM supplemented with 10% fetal bovine serum, 50 U/mL of penicillin, and 50 µg/mL of streptomycin.</p>
</sec>
<sec id="s2.3"><title>Extraction of episomal and total DNA</title>
<p>To extract episomal HPV DNA from W12 cells, a modified Hirt method was used [<xref ref-type="bibr" rid="B19">19</xref>
]. Briefly, subconfluent cells on a 6-cm dish were digested with 500 μL of lysis buffer consisting of 25 mM Tris-HCl (pH 7.5), 5 mM EDTA, 0.6% SDS, and 50 μg/mL of RNase A at room temperature for 5 min. After the addition of 350 μL of precipitation buffer consisting of 3 M CsCl, 1 M potassium acetate, and 0.67 M acetic acid, the solution was incubated on ice for 15 min, and precipitated by centrifugation (14,000 x g, 15 min, 4°C). Cleared supernatant was applied onto a Miniprep column (Favorgen, Ping-Tung, Taiwan) and passed through the column by centrifugation (14,000 x g, 1 min). The column was then washed with 750 μL of wash buffer consisting of 10 mM Tris-HCl (pH 7.5), 80 mM potassium acetate, 40 μM EDTA, and 60% ethanol. DNA was eluted with 60 μL of elution buffer consisting of 10 mM Tris-HCl (pH 7.5) and 0.1 mM EDTA. The resultant DNA fraction mainly contains extrachromosomal DNAs. </p>
<p>Total DNA was isolated from cervical exfoliated cells using the QIAamp DNA blood kit (QIAGEN, Hilden, Germany). The DNA samples were genotyped by PGMY-reverse blot hybridization assay [<xref ref-type="bibr" rid="B20">20</xref>
,<xref ref-type="bibr" rid="B21">21</xref>
], and HPV16-positive samples without other detectable HPVs from LSIL and ICC specimens were used for subsequent PCR. Total DNA was extracted from CaSki cells as described above for clinical samples. The concentrations of extracted DNAs were about 80 ng/μL from the cultured cells and 2 to 4 ng/μL from clinical samples.</p>
</sec>
<sec id="s2.4"><title>Full-circle PCR</title>
<p>PCR was performed in 10 μL consisting of 1× PrimeSTAR<sup>®</sup>
GXL Buffer (Takara, Ohtsu, Japan), 200 μM dNTP mixture, 0.25 μM forward primer, 0.25 μM reverse primer, 0.25 U of PrimeSTAR<sup>®</sup>
GXL DNA polymerase (Takara), and 2 to 4 ng of template DNA. Primers used are as follows: 1742F, <named-content content-type="gene">5’-tgc tgt cta aac tat tat gtg tgt ctc-3’</named-content>
; 1873R, <named-content content-type="gene">5’-gcg tgt ctc cat aca ctt ca-3’</named-content>
; 1869F, <named-content content-type="gene">5’-cac gcc aga atg gat aca aa-3’</named-content>
; 1746R, <named-content content-type="gene">5’-cag caa ttt ttc aat tgt ttc tc-3’</named-content>
. Positions of 5’-end nucleotide of primers in the reference HPV16 genome (HPV16REF in Papillomavirus <italic>Episteme</italic>
, <ext-link ext-link-type="uri" xlink:href="http://pave.niaid.nih.gov"><underline>http://pave.niaid.nih.gov</underline>
</ext-link>
) are used as primer names with either forward (F) or reverse (R) orientation. The plasmid containing the full-length HPV16 genome, HPV16/pUC19 [<xref ref-type="bibr" rid="B22">22</xref>
], was used as a template for positive PCR. The PCR protocol was 98°C for 30 s, followed by 30 cycles of 98°C for 10 s, 60°C for 15 s, 68°C for 2 min. The PCR products were separated on a 0.7% agarose gel in 1 x TAE buffer and visualized with ethidium bromide. </p>
</sec>
<sec id="s2.5"><title>Deep sequencing</title>
<p>Approximately 8-kb amplicons generated by full-circle PCR with the primer-pair 1742F/1873R were separated on a 0.7% agarose gel and purified using the Wizard gel purification kit (Promega, Madison, WI). To prepare DNA libraries compatible with Illumina NGS, the purified DNA (16 ng) was fragmented and adaptor-ligated using the Nextra DNA sample prep kit (Epicentre, Madison, WI), followed by DNA purification with the QIAquick PCR purification kit (QIAGEN). The resultant DNA libraries were amplified using a minimal PCR amplification step with Phusion High-Fidelity DNA polymerase (New England Biolab, Ipswich, MA) as follows: 72°C for 3 min, 95°C for 30 s, followed by 13 cycles of 95°C for 10 s, 62°C for 30 s, 72°C for 3 min. Adaptor-ligated DNAs in the range of 500 to 700 bp were then size-selected by 1% agarose gel electrophoresis and purified using the Wizard kit. Five to 8 libraries from different clinical samples were multiplexed with adaptors containing different index sequences and simultaneously subjected to deep sequencing. Cluster generation and sequencing were performed for 151 cycles on a MiSeq sequencer (Illumina, San Diego, CA) with the MiSeq Reagent kit (300 cycle) (Illumina) according to the manufacturer’s instructions. Fluorescent images were analyzed using the Illumina MCS1.1/RTA1.13.56 base-calling pipeline to obtain FASTQ-formatted sequence data. To construct contiguous DNA sequences (contigs), total sequences of paired-end 151-mer reads were assembled using ABySS-pe v1.3.3 [<xref ref-type="bibr" rid="B23">23</xref>
] with the following parameters: k80, n1200, c1000, t10, e10, q20, and s160; or k80, n120, c100, t10, e10, q20, and s160. The assembly generated 2 to 20 contigs, followed by further assembly of these sequences into the complete circular sequence by Phrap [<xref ref-type="bibr" rid="B24">24</xref>
]. </p>
<p>The fidelity of <italic>de novo</italic>
assembly was examined by alignment of paired-end reads to the assembled sequence using Burrows-Wheeler Aligner (BWA) v0.6.1 [<xref ref-type="bibr" rid="B25">25</xref>
]. Prior to read mapping, read sequences of Phred quality score <underline>></underline>
20 were selected by PoPoolation v1.2.2 [<xref ref-type="bibr" rid="B26">26</xref>
], and trimmed for 20 bases at 5’-terminus and 51 bases at 3’-terminus because 5’ and 3’-terminus of reads generate relatively low fidelity sequences. The processed 80-base sequences were mapped to the assembled sequence as a paired-end mode by BWA. The mapping results were visualized by Integrative Genomics Viewer (IGV) v2.0.34 [<xref ref-type="bibr" rid="B27">27</xref>
] or Tablet v1.12.12.05 [<xref ref-type="bibr" rid="B28">28</xref>
]. The numbers of paired-end reads obtained from each sample and used for <italic>de novo</italic>
assembly are described in <xref ref-type="table" rid="pone-0080583-t001">Table 1</xref>
, and the average depth of the HPV16 genome ranged from x 7,624 for ICC sample 6 to x 48,221 for W12 sample. Nucleotide mismatches compared to the assembled reference genome and positions of SNPs in each sample were identified by using SAMtools v0.1.18 [<xref ref-type="bibr" rid="B29">29</xref>
] with in-house Perl scripts (available upon request) or VarScan v2.3.2 [<xref ref-type="bibr" rid="B30">30</xref>
]. Based on a quality score confidence threshold of Phred quality score > 30 (error probability <underline><</underline>
0.001) that was used to extract variation positions in the read sequences, we defined a position as homogeneous if the mutation frequency is < 0.5% and a position to be heterogeneous if the mutation frequency is > 0.5%. The presence of nucleotide substitutions and insertions/deletions was finally verified by visual inspection of mismatched read sequences using IGV and Tablet.</p>
<table-wrap id="pone-0080583-t001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0080583.t001</object-id>
<label>Table 1</label>
<caption><title><italic>De novo</italic>
assembly of complete HPV16 genome sequences from short-read sequence data.</title>
</caption>
<table frame="hsides" rules="groups" cellpadding="1"><colgroup span="1"><col span="1"></col>
<col span="1"></col>
<col span="1"></col>
<col span="1"></col>
<col span="1"></col>
<col span="1"></col>
<col span="1"></col>
</colgroup>
<thead><tr><th rowspan="1" colspan="1">Sample</th>
<th rowspan="1" colspan="1">Cytology</th>
<th rowspan="1" colspan="1">Read number</th>
<th rowspan="1" colspan="1">Length (bp)</th>
<th rowspan="1" colspan="1">Variant lineage<xref ref-type="table-fn" rid="ngtab1.1">*</xref>
</th>
<th rowspan="1" colspan="1">Average depth (per nucleotide)</th>
<th rowspan="1" colspan="1">DDBJ accession</th>
</tr>
</thead>
<tbody><tr><td rowspan="1" colspan="1">W12</td>
<td rowspan="1" colspan="1">LSIL</td>
<td rowspan="1" colspan="1">2,539,664</td>
<td rowspan="1" colspan="1">7,904</td>
<td rowspan="1" colspan="1">EUR (A1)</td>
<td rowspan="1" colspan="1">48,221</td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td rowspan="1" colspan="1">#1</td>
<td rowspan="1" colspan="1">LSIL</td>
<td rowspan="1" colspan="1">2,057,386</td>
<td rowspan="1" colspan="1">7,903</td>
<td rowspan="1" colspan="1">EUR (A3)</td>
<td rowspan="1" colspan="1">39,064</td>
<td rowspan="1" colspan="1">AB818687</td>
</tr>
<tr><td rowspan="1" colspan="1">#2</td>
<td rowspan="1" colspan="1">LSIL</td>
<td rowspan="1" colspan="1">1,739,070</td>
<td rowspan="1" colspan="1">7,906</td>
<td rowspan="1" colspan="1">EUR (A3)</td>
<td rowspan="1" colspan="1">33,020</td>
<td rowspan="1" colspan="1">AB818688</td>
</tr>
<tr><td rowspan="1" colspan="1">#3</td>
<td rowspan="1" colspan="1">LSIL</td>
<td rowspan="1" colspan="1">3,577,966</td>
<td rowspan="1" colspan="1">7,907</td>
<td rowspan="1" colspan="1">NA (D1)</td>
<td rowspan="1" colspan="1">36,233</td>
<td rowspan="1" colspan="1">AB818689</td>
</tr>
<tr><td rowspan="1" colspan="1">#4</td>
<td rowspan="1" colspan="1">LSIL</td>
<td rowspan="1" colspan="1">921,386</td>
<td rowspan="1" colspan="1">7,903</td>
<td rowspan="1" colspan="1">AFR2 (C)</td>
<td rowspan="1" colspan="1">17,495</td>
<td rowspan="1" colspan="1">AB818690</td>
</tr>
<tr><td rowspan="1" colspan="1">#5</td>
<td rowspan="1" colspan="1">LSIL</td>
<td rowspan="1" colspan="1">1,913,428</td>
<td rowspan="1" colspan="1">7,906</td>
<td rowspan="1" colspan="1">As (A4)</td>
<td rowspan="1" colspan="1">19,376</td>
<td rowspan="1" colspan="1">AB818691</td>
</tr>
<tr><td rowspan="1" colspan="1">#6</td>
<td rowspan="1" colspan="1">ICC</td>
<td rowspan="1" colspan="1">401,546</td>
<td rowspan="1" colspan="1">7,905</td>
<td rowspan="1" colspan="1">As (A4)</td>
<td rowspan="1" colspan="1">7,624</td>
<td rowspan="1" colspan="1">AB818692</td>
</tr>
<tr><td rowspan="1" colspan="1">#7</td>
<td rowspan="1" colspan="1">ICC</td>
<td rowspan="1" colspan="1">1,349,444</td>
<td rowspan="1" colspan="1">7,905</td>
<td rowspan="1" colspan="1">As (A4)</td>
<td rowspan="1" colspan="1">25,622</td>
<td rowspan="1" colspan="1">AB818693</td>
</tr>
<tr><td rowspan="1" colspan="1">CaSki</td>
<td rowspan="1" colspan="1">ICC</td>
<td rowspan="1" colspan="1">1,394,918</td>
<td rowspan="1" colspan="1">7,905</td>
<td rowspan="1" colspan="1">EUR (A1)</td>
<td rowspan="1" colspan="1">26,486</td>
<td rowspan="1" colspan="1"></td>
</tr>
</tbody>
</table>
<table-wrap-foot><fn id="ngtab1.1"><label>*</label>
<p>Sublineage designations based on recent classification by Burk et al. [<xref ref-type="bibr" rid="B51">51</xref>
] are in parenthesis.</p>
</fn>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="s2.6"><title>Transient HPV16 replication assay</title>
<p>An HPV16 origin-positive plasmid was constructed by insertion of a <italic>Pst</italic>
I fragment of the HPV16 genome (from 7,005 to 7,906 and 1 to 880) into pGL4.50 (Promega). Expression plasmids for N-terminal FLAG-tagged HPV16 E1 and N-terminal FLAG-tagged HPV16 E2 were constructed by cloning of the codon-optimized cDNA of E1 and E2 into pCMV-β (Clontech, Palo Alto, CA). Expression plasmids for E1 mutants were generated by using the QuickChange Lightning Multi Site-Directed Mutagenesis Kit (Agilent Technologies, La Jolla, CA). Transient replication assays were performed as described by Fradet-Turcotte et al. [<xref ref-type="bibr" rid="B31">31</xref>
] with some modifications. Briefly, HEK293 cells were plated 20 h before transfection in a 24-well plate at a density of 32,000 cells/well, and transfected with 2 ng of the origin-positive plasmid and 2 ng of pGL4.75 (Promega), together with 20 ng of the E2 expression plasmid and the indicated amounts of E1 expression plasmid using the FuGENE6 reagent (Promega). The total quantity of transfected plasmid DNA was adjusted to 100 ng with the empty plasmid p3xFLAG-CMV10 (Sigma-Aldrich, St. Louis, MO) as carrier DNA. At 72 h after transfection, firefly and <italic>Renilla</italic>
luciferase activities were measured using the Dual-Glo Luciferase assay system (Promega) on an ARVO MX luminescence counter (PerkinElmer, Waltham, MA), and the level of replication was quantified as the ratio of the two luciferase activities.</p>
</sec>
<sec id="s2.7"><title>Nucleotide sequence accession number</title>
<p>All assembled HPV genome sequences have been deposited in the DNA Data Bank of Japan (DDBJ) under accession numbers AB818687 to AB818693 for HPV16, AB819272 to AB819274 for HPV52, and AB819275 to AB819279 for HPV58. The sequence reads of the HPV genomes are available from the DDBJ Sequence Read Archive under accession number DRA001083.</p>
</sec>
</sec>
<sec sec-type="results" id="s3"><title>Results</title>
<sec id="s3.1"><title>Full-circle PCR amplifying full-length HPV16 genome sequences</title>
<p>Using PrimeSTAR<sup>®</sup>
GXL polymerase, which is specifically designed for the fast and accurate amplification of long DNA sequences, and the HPV16-specific primers 1742F and 1873R, we generated two amplicons from HPV16/pUC19 plasmid DNA: a short fragment of ~130 bp and a large fragment consisting of the entire HPV16/pUC19 sequence of approximately 11 kb (with ~130-bp overlap) (<xref ref-type="fig" rid="pone-0080583-g001">Figure 1A</xref>
, lanes 2 to 7, and see <xref ref-type="fig" rid="pone-0080583-g001">Figure 1B</xref>
). This PCR reaction is slightly less efficient at creating the large amplicon than the same reaction using the non-overlapping primers 1869F and 1746R, which create only one amplicon (<xref ref-type="fig" rid="pone-0080583-g001">Figure 1A</xref>
, lanes 9 to 14). The reactions with PrimeSTAR<sup>®</sup>
GXL polymerase were both specific (only 1 or 2 PCR amplicons, depending on the primers used) and sensitive (we successfully amplified as little as 0.01 pg of input template DNA, which corresponds to about 900 copies of the plasmid). The introduction of cellular DNA did not inhibit this reaction, even at the sub-picogram level of template DNA (<xref ref-type="supplementary-material" rid="pone.0080583.s002">Figure S1</xref>
). </p>
<fig id="pone-0080583-g001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0080583.g001</object-id>
<label>Figure 1</label>
<caption><title>Amplification of full-length HPV16 genomes by full-circle PCR.</title>
<p>(A) PCR was performed with PrimeSTAR<sup>®</sup>
GXL DNA polymerase and HPV16-specific primer-pairs as indicated. The amounts of HPV16/pUC19 used for the PCR template are also indicated above. The PCR products were analyzed by agarose gel electrophoresis. M, DNA size markers. (B) Scheme for full-circle PCR. PrimeSTAR<sup>®</sup>
GXL DNA polymerase generates short and long DNA products with primer-pair 1742F/1873R. (C) Full-circle PCR with DNA extracted from W12 cells, clone 20863 (high-copy HPV16 episomes) (lane 2) and clone 20850 (low-copy HPV16 episomes) (lane 3). M, DNA size marker (lanes 1) (D) Full-circle PCR using DNA isolated from 7 clinical specimens: 5 LSIL (lanes 2 to 6), and 2 ICC (lanes 7 and 8). M, DNA size marker (lanes 1).</p>
</caption>
<graphic xlink:href="pone.0080583.g001"></graphic>
</fig>
<p>Using the primer-pair 1742F/1873R and PrimeSTAR<sup>®</sup>
GXL polymerase, we also amplified two DNA fragments from extrachromosomal HPV16 DNA isolated from W12 cells that were originally isolated from an LSIL biopsy (<xref ref-type="fig" rid="pone-0080583-g001">Figure 1C</xref>
). The short and long amplicons matched the predicted sizes of 132 bp and ~8 kb, respectively (<xref ref-type="fig" rid="pone-0080583-g001">Figure 1B</xref>
), and the long amplicon covered the complete HPV16 genome, as confirmed by Sanger sequencing (results not shown). Hence we refer to this overlapping PCR hereafter as full-circle PCR. To further test the utility of full-circle PCR with PrimeSTAR<sup>®</sup>
GXL polymerase, DNA prepared from CaSki cells, known to contain 60 to 600 copies of integrated tandem repeated HPV16 genomes [<xref ref-type="bibr" rid="B32">32</xref>
,<xref ref-type="bibr" rid="B33">33</xref>
], was used as the PCR template. As expected, a mixture of 8 and 6.5-kb amplicons were obtained (<xref ref-type="supplementary-material" rid="pone.0080583.s003">Figure S2</xref>
). </p>
<p>Importantly, full-circle PCR was able to generate approximately 8-kb amplicons using DNA isolated from HPV16-positive cervical exfoliated cells in patients with either LSIL (<xref ref-type="fig" rid="pone-0080583-g001">Figure 1D</xref>
, lanes 2 to 6) or ICC (lanes 7 and 8). This observation prompted us to use these long PCR products for genetic variation analyses of full-length HPV16 genomes derived from clinical samples.</p>
</sec>
<sec id="s3.2"><title><italic>De novo</italic>
assembly of full-length HPV16 genome sequences by deep sequencing</title>
<p>In order to comprehensively analyze HPV genetic variation through deep sequencing, we sought a methodology capable of determining the full-length reference sequence of the HPV16 genome from high-throughput short-read sequence data. First, to develop the system, the 8-kb DNA amplified from W12 DNA by full-circle PCR was fragmented and adaptor-ligated to yield a DNA library compatible with Illumina next generation sequencers. Deep sequencing of the resultant library generated sequence data from more than a million paired-end reads. <italic>De novo</italic>
assembly of those paired-end reads by ABySS-pe and Phrap generated one contiguous sequence of 7,904 bp length, which was a perfect match for the full-length W12 HPV16 genome sequence in GenBank (accession No. AF125673). Next, the long PCR products from CaSki DNA were subjected to deep sequencing analysis. The <italic>de novo</italic>
assembly procedure once again resulted in construction of one contiguous sequence of 7,905-bp length, identical to the full-length CaSki HPV16 genome sequence in GenBank (accession No. U89348). Mapping of the total read sequences to the CaSki reference sequence demonstrated the presence of two large deletions; from nucleotide positions 1 to 470 and 6,907 to 7,905 (<xref ref-type="supplementary-material" rid="pone.0080583.s004">Figure S3</xref>
), which perfectly matched a deleted region previously reported in CaSki cells [<xref ref-type="bibr" rid="B33">33</xref>
]. Thus, it is likely that the 6.5-kb DNA products obtained from CaSki DNA contain this 1,469-bp deletion. Overall, these results demonstrate the accuracy of long PCR, deep sequencing and bioinformatics to obtain full-length HPV16 genome sequences. </p>
<p>To determine the complete genome sequence of HPV16 contained in clinical specimens, we analyzed the 8-kb full-circle PCR products from 7 clinical samples by deep sequencing and assembled the complete viral genome sequence as described above. <xref ref-type="table" rid="pone-0080583-t001">Table 1</xref>
presents the results of this analysis and demonstrates that <italic>de novo</italic>
assembly yielded one contiguous sequence for each sample. The length of each sequence ranged from 7,903 to 7,907 bp due to the differences in the length of the non-coding region between E5 and L2 (<xref ref-type="supplementary-material" rid="pone.0080583.s005">Figure S4</xref>
). </p>
<p>The HPV16 genome is classified into 4 major variant-lineages based on a unique combination of single nucleotide polymorphisms (SNPs): European-Asian, including the sublineages European (EUR) and Asian (As), African 1 (AFR1), African 2 (AFR2), and North-American/Asian-American (NA/AA) [<xref ref-type="bibr" rid="B34">34</xref>
]. As presented in <xref ref-type="table" rid="pone-0080583-t001">Table 1</xref>
, SNP analysis of the assembled HPV16 genomes assigned the following variants to each clinical sample: 2 LSIL samples, EUR; 1 LSIL sample, NA; 1 LSIL sample, AFR2; and 1 LSIL and 2 ICC samples, As. Furthermore, phylogenetic tree analysis with the assembled complete genome sequences of HPV16 confirmed the assignments in variant-lineages (<xref ref-type="supplementary-material" rid="pone.0080583.s006">Figure S5</xref>
). These results indicate that one HPV16 variant predominates in each clinical specimen. </p>
</sec>
<sec id="s3.3"><title>Analysis of nucleotide variations in HPV16 genomes</title>
<p> Before analyzing the sequence variation present in clinical samples, we first validated the efficacy and error rates of deep sequencing on our platform. For this purpose, we used a linearized full-length HPV16 DNA that was excised from HPV16/pUC19 by <italic>Bam</italic>
HI digestion. The read sequences obtained were aligned to the prototype HPV16 sequence (GenBank accession No. K02718), and the number of reads containing nucleotide mismatches was determined at each position of the viral genome. This analysis enabled us to calculate mutation frequencies at each nucleotide position (the average depth of reads count per nucleotide position was 7,035). As shown in <xref ref-type="fig" rid="pone-0080583-g002">Figure 2</xref>
, an overall profile of the mutation/error frequency demonstrated an error rate below 0.3% distributed throughout the whole sequence, which reflects intrinsic errors generated in the deep sequencing analysis. The average error rate per nucleotide position throughout the whole sequence was determined to be 0.0070%. Visualization and manual inspection of aligned read sequences revealed that the two minor peaks observed at nt positions 6,242 and 6,434 (see <xref ref-type="fig" rid="pone-0080583-g002">Figure 2</xref>
) were derived from mis-mapping of reads to ends of the linear HPV16 DNA (<italic>Bam</italic>
HI site, at nt 6,154). </p>
<fig id="pone-0080583-g002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0080583.g002</object-id>
<label>Figure 2</label>
<caption><title>Mutation frequency profile of full-length HPV16 genomes prepared from a plasmid.</title>
<p>The read sequences obtained with full-length HPV16 genomes prepared from a cloned plasmid were aligned to the reference HPV16 sequence, and mutation/error frequencies at each nucleotide position relative to the total coverage (number of reads that encompass each nucleotide, up to the maximum x 8,000) are presented in the landscape of the full-length HPV16 genome. A threshold line for a reliable mutation frequency (0.5%) is indicated with the red dotted line. The genome organization of HPV16 is indicated above: p97, the early promoter; p670, the late promoter; polyA(early) and polyA(late), the early and late polyadenylation signals, respectively. </p>
</caption>
<graphic xlink:href="pone.0080583.g002"></graphic>
</fig>
<p>Our deep sequencing analysis of HPV16 DNA amplified by PrimeSTAR<sup>®</sup>
GXL polymerase demonstrated that the average error rate was 0.0076% (results not shown), which is comparable to the error rate reported by the manufacturer (0.0062%), indicating high fidelity DNA amplification by this polymerase. Based upon this low error rate and our analytical method’s ability to pick up nucleotide variation with a minimum confidence score threshold of Phred quality score 30 (error probability <underline><</underline>
0.001), we decided to use 0.5% as a threshold value of mutation frequency to distinguish real mutations from putative errors. </p>
<p>We applied this variation analysis method to HPV16 genome sequences obtained from W12 cells. As shown in <xref ref-type="fig" rid="pone-0080583-g003">Figure 3</xref>
, alignment of the paired-end reads to the W12 reference sequence revealed no major nucleotide substitutions with a frequency >0.5% in the read coverage , suggesting that HPV16 genome sequences are highly homogeneous in W12 cells. </p>
<fig id="pone-0080583-g003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0080583.g003</object-id>
<label>Figure 3</label>
<caption><title>Mutation frequency profile of full-length HPV16 genomes in W12 cells.</title>
<p>The read sequences obtained with full-length HPV16 genomes prepared from W12 cells were aligned to the reference HPV16 sequence (AF125673), and mutation/error frequencies at each nucleotide position are presented in the landscape of the full-length HPV16 genome. A threshold line for a reliable mutation frequency (0.5%) is indicated with the red dotted line. The genome organization of HPV16 is indicated above.</p>
</caption>
<graphic xlink:href="pone.0080583.g003"></graphic>
</fig>
<p> Next we employed this established analysis procedure to examine HPV genomic variation in clinical samples. We report two types of mutation profiles found in 7 clinical samples: a highly homogeneous profile as observed with W12 cells, and a profile with 1 to 3 specific positions of nucleotide substitution. Four samples (#1, #4, #5, and #7) showed an almost homogeneous profile without specific nucleotide substitution (two representative profiles are presented in <xref ref-type="supplementary-material" rid="pone.0080583.s007">Figure S6</xref>
) and 3 samples (#2, #3, and #6) demonstrated a mutation profile with several positions of nucleotide substitution (<xref ref-type="fig" rid="pone-0080583-g004">Figure 4</xref>
). The nucleotide substitutions were located in the early region of the HPV16 genome, E6 and E1, and the non-coding region between E5 and L2, with frequencies from 0.60 to 5.42% (<xref ref-type="table" rid="pone-0080583-t002">Table 2</xref>
). Intriguingly, the nucleotide substitutions observed in the E1 gene of sample 6 cause amino-acid (aa) changes in the E1 protein: methionine at aa position 326 to isoleucine, and glutamine at aa position 381 to glutamic acid. The nucleotide substitution in the E6 gene detected in sample 6 also leads to an aa change in the E6 protein, glutamic acid to aspartic acid at aa position 25.</p>
<fig id="pone-0080583-g004" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0080583.g004</object-id>
<label>Figure 4</label>
<caption><title>Mutation frequency profile of full-length HPV16 genomes in clinical specimens.</title>
<p>The read sequences obtained with full-length HPV16 genomes prepared from clinical specimens (2 LSIL samples, #2 and #3; 1 ICC sample, #6) were aligned to their <italic>de</italic>
<italic>novo</italic>
assembled complete genome sequences, and mutation/error frequencies at each nucleotide position are presented in the landscape of the full-length HPV16 genome. A threshold line for a reliable mutation frequency (0.5%) is indicated with the red dotted line. Peaks above 0.5% are indicated with red asterisks. The genome organization of HPV16 is indicated above.</p>
</caption>
<graphic xlink:href="pone.0080583.g004"></graphic>
</fig>
<table-wrap id="pone-0080583-t002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0080583.t002</object-id>
<label>Table 2</label>
<caption><title>Nucleotide and amino acid changes detected in individual clinical specimens.</title>
</caption>
<table frame="hsides" rules="groups"><colgroup span="1"><col span="1"></col>
<col span="1"></col>
<col span="1"></col>
<col span="1"></col>
<col span="1"></col>
<col span="1"></col>
<col span="1"></col>
</colgroup>
<thead><tr><th rowspan="1" colspan="1">Sample</th>
<th rowspan="1" colspan="1">Nt position</th>
<th rowspan="1" colspan="1">Nt change</th>
<th rowspan="1" colspan="1">*Frequency (%)</th>
<th rowspan="1" colspan="1">Aa position</th>
<th rowspan="1" colspan="1">Aa change</th>
<th rowspan="1" colspan="1">Gene</th>
</tr>
</thead>
<tbody><tr><td rowspan="1" colspan="1">#2</td>
<td rowspan="1" colspan="1">256</td>
<td rowspan="1" colspan="1">C to T</td>
<td rowspan="1" colspan="1">0.60</td>
<td rowspan="1" colspan="1">51</td>
<td rowspan="1" colspan="1">no change</td>
<td rowspan="1" colspan="1">E6</td>
</tr>
<tr><td rowspan="1" colspan="1">#3</td>
<td rowspan="1" colspan="1">4158</td>
<td rowspan="1" colspan="1">A to G</td>
<td rowspan="1" colspan="1">0.70</td>
<td rowspan="1" colspan="1">-</td>
<td rowspan="1" colspan="1">no change</td>
<td rowspan="1" colspan="1">-</td>
</tr>
<tr><td rowspan="1" colspan="1">#6</td>
<td rowspan="1" colspan="1">178</td>
<td rowspan="1" colspan="1">G to T</td>
<td rowspan="1" colspan="1">0.72</td>
<td rowspan="1" colspan="1">25</td>
<td rowspan="1" colspan="1">Glu to Asp</td>
<td rowspan="1" colspan="1">E6</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">1842</td>
<td rowspan="1" colspan="1">G to A</td>
<td rowspan="1" colspan="1">1.44</td>
<td rowspan="1" colspan="1">326</td>
<td rowspan="1" colspan="1">Met to Ile</td>
<td rowspan="1" colspan="1">E1</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">2005</td>
<td rowspan="1" colspan="1">C to G</td>
<td rowspan="1" colspan="1">5.42</td>
<td rowspan="1" colspan="1">381</td>
<td rowspan="1" colspan="1">Gln to Glu</td>
<td rowspan="1" colspan="1">E1</td>
</tr>
</tbody>
</table>
<table-wrap-foot><fn id="ngtab2.0"><p>* Mutation frequencies were calculated by using SAMtools with in-house Perl scripts.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p> We were also able to employ our full-circle PCR and deep sequencing approach to analyze clinical samples that contained HPV52 and HPV58 DNA (<xref ref-type="supplementary-material" rid="pone.0080583.s008">Figure S7</xref>
and <xref ref-type="supplementary-material" rid="pone.0080583.s001">Table S1</xref>
), suggesting that this method is not specific to the analysis of HPV16 and that it is also applicable to the oncogenic HPV types prevalent in East Asian countries including China and Japan [<xref ref-type="bibr" rid="B35">35</xref>
-<xref ref-type="bibr" rid="B37">37</xref>
]. </p>
</sec>
<sec id="s3.4"><title>Effect of E1 amino-acid changes on its replication function</title>
<p> The <italic>de novo</italic>
assembled HPV16 genomes from clinical specimens revealed aa polymorphisms of the E1 protein in its DNA-binding domain at aa positions 294 and 326: a replacement of leucine (prototype) at position 294 with methionine (L294M) in 2 LSIL specimens (#3 and #4) and that of isoleucine (prototype) at position 326 with methionine (I326M) in all the specimens (<xref ref-type="fig" rid="pone-0080583-g005">Figure 5A</xref>
). Further, nucleotide variation analyses showed a novel aa variation at position 381 from glutamine (prototype) to glutamic acid (Q381E) in sample 6. To explore the biological significance of these E1 variations, we mutated the corresponding aa residues in the European prototype E1 protein and tested the mutated E1’s ability to support HPV16 origin-dependent DNA replication. The E1 mutant Y379F (in which tyrosine at 379 is replaced with phenylalanine), which is expected to lose replication activity as reported previously [<xref ref-type="bibr" rid="B38">38</xref>
], was also included as a representative of a replication-deficient E1 protein. </p>
<fig id="pone-0080583-g005" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0080583.g005</object-id>
<label>Figure 5</label>
<caption><title>Characterization of E1 variants by transient HPV16 replication assay.</title>
<p>(A) Left panel: schematic representation of the HPV16 E1 protein. The prototype E1 protein has leucine at position 294 and isoleucine at position 326, both of which are located in the DNA-binding domain, and glutamine at position 381 located in the oligomerization domain. Right panel: summary of aa variations found in the E1 protein from the clinical samples (#1 to #7). Proto, the E1 protein of European prototype. (B) E1-dependent replication of HPV16 origin-positive plasmid. The HPV16 origin-positive plasmid expressing Firefly luciferase (Fluc) and the origin-deficient plasmid expressing <italic>Renilla</italic>
luciferase (Rluc) were transfected into HEK293 cells with increasing amounts of the prototype E1 expression plasmid together with the E2 expression plasmid. At 72 h after transfection, the ratio of the two luciferase activities (Fluc/Rluc) was measured as levels of the origin-dependent replication. The origin-deficient plasmid expressing Firefly luciferase (pGL4.50) was used as a negative control for HPV replication. Statistically significant differences (Welch’s t-test, <italic>p</italic>
<0.01) are indicated with *. (C) Replication activity of E1 variants. Increasing amounts of expression plasmids for FLAG-tagged prototype or variant E1 proteins were transfected into HEK293 cells, and the levels of replication were measured at 72 h after transfection. Error bars represent the standard deviation of triplicate transfections. Statistically significant differences (Welch’s t-test, <italic>p</italic>
<0.01) are indicated with *. The data are representative of three independent experiments. (D) Western blot analysis of E1 variants. FLAG-tagged prototype or variant E1 proteins expressed in HEK293 cells were detected with anti-FLAG antibody (Sigma-Aldrich) and the ECL prime Western blotting detection reagent (GE Healthcare, Buckinghamshire, England). Tubulin was visualized with anti-α/β-tubulin antibody (Cell Signaling, Danvers, MA) as loading control. </p>
</caption>
<graphic xlink:href="pone.0080583.g005"></graphic>
</fig>
<p>In the transient HPV16 replication assay, two reporter plasmids were transfected into HEK293 cells together with expression vectors for FLAG-tagged E1 and E2. Firefly luciferase (Fluc) is expressed from the HPV16 origin-positive plasmid, whereas <italic>Renilla</italic>
luciferase (Rluc) is expressed from the origin-deficient plasmid, so the ratio of the two luciferase activities (Fluc/Rluc) in the presence of the E1/E2 expression vectors represents levels of HPV origin-dependent replication. As previously reported [<xref ref-type="bibr" rid="B31">31</xref>
], the Fluc/Rluc ratio obtained with the origin-positive plasmid increased in a manner dependent on the amount of the E1 expression plasmid (<xref ref-type="fig" rid="pone-0080583-g005">Figure 5B</xref>
). Under these assay conditions, two E1 variants, L294M and I326M, showed replication activities almost comparable to the parental E1 (<xref ref-type="fig" rid="pone-0080583-g005">Figure 5C</xref>
). In contrast, Q381E exhibited a reduced activity for supporting HPV16 replication, and Y379F completely lost the replication activity as expected. Western blot analyses revealed that the expression levels of I326M, Q381E, and Y379F were comparable to that of the prototype E1 (<xref ref-type="fig" rid="pone-0080583-g005">Figure 5D</xref>
), indicating that the reduced replication activity of Q381E is not due to reduced protein level. </p>
</sec>
<sec id="s3.5"><title>Detection of deletions in HPV16 genomes</title>
<p>Because our deep sequencing analysis was able to precisely define the deletion in the HPV16 genome of CaSki cells, we examined HPV16 genome deletions in the 7 clinical samples by read mapping to the assembled reference sequence. Of the 5 LSIL samples, 4 LSIL samples did not show any deletions and 1 LSIL sample (#2) demonstrated the discontinuous distribution of read depth indicative of genome deletions (<xref ref-type="supplementary-material" rid="pone.0080583.s004">Figure S3</xref>
, arrows). Further detailed visualization of the mapped read sequences revealed that 27 bases were deleted from nucleotide positions 3,476 to 3,502 in 78.9% of total read sequences covering this region (<xref ref-type="fig" rid="pone-0080583-g006">Figure 6A</xref>
). Independent full-circle PCR for sample 2 and subsequent deep sequencing analysis yielded a comparable deletion frequency (86.6%, results not shown), demonstrating the reproducibility of the analysis overall. PCR amplification of the E2 region (<xref ref-type="supplementary-material" rid="pone.0080583.s009">Figure S8</xref>
) and subsequent Sanger sequencing confirmed the presence of the deletion (results not shown). In addition, read sequences lacking 1 to 3 bases in the region between E5 and L2 were also detected in sample 2 in 89.0% of the total read coverage of that region (<xref ref-type="fig" rid="pone-0080583-g006">Figure 6B</xref>
). Although the read depth in the E5 region was low (<xref ref-type="supplementary-material" rid="pone.0080583.s004">Figure S3</xref>
), no deletion was found in this region. These data suggest that deleted viral genomes are present as a mixed state with intact full-length HPV16 genomes in this clinical specimen. In contrast to HPV16 genomes in CaSki cells, no deletion was detected in HPV16 genome sequences obtained from the 2 ICC samples.</p>
<fig id="pone-0080583-g006" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0080583.g006</object-id>
<label>Figure 6</label>
<caption><title>HPV16 genome deletions in a clinical specimen revealed by read mapping.</title>
<p>Mapping of paired-end read sequences to <italic>de</italic>
<italic>novo</italic>
assembled full-length HPV16 genome sequence for LSIL sample 2 is visualized by Tablet. Deletions in the E2/E4 region (A) and in the region between E5 and L2 (B) are shown. (C) The genome organization of HPV16. The deleted regions in E2/E4 and between E5 and L2 are indicated with red and blue dotted-line boxes, respectively.</p>
</caption>
<graphic xlink:href="pone.0080583.g006"></graphic>
</fig>
</sec>
</sec>
<sec sec-type="discussion" id="s4"><title>Discussion</title>
<p>In this study, utilizing a recently developed long PCR enzyme, PrimeSTAR<sup>®</sup>
GXL polymerase, we have succeeded in amplifying full-length HPV16, HPV52, and HPV58 genome sequences, using DNA isolated from clinical specimens, in a single reaction. The use of this highly-processive DNA polymerase is critical, since traditional PCR enzymes such as <italic>Taq</italic>
DNA polymerase are unable to amplify such long DNAs. Furthermore, by applying the NGS technology to the amplified HPV DNA, we have developed a methodology capable of analyzing the genetic variation present in individual clinical samples. Although library construction is still laborious, an improved Nextera DNA sample preparation kit (Illumina) facilitates rapid and simple library preparation for up to 96 samples with different index primers, thus allowing each sample to be analyzed in a single NGS run and further reducing the overall cost. </p>
<p>The <italic>de novo</italic>
assembly procedure described in this study provides a simple and reliable methodology to construct a reference sequence of the full-length HPV genome from millions of short-read sequences. Construction of a contiguous reference sequence from short-read sequence data is critical for both determining the complete genome sequence of known or unknown HPVs and examining sequence variations in HPV genomes compared to the reference sequence. In this regard, it is worth noting that the overlapping sequence at both ends of the long PCR products (see <xref ref-type="fig" rid="pone-0080583-g001">Figure 1B</xref>
) contributes to an accurate and faithful reconstruction of the circular HPV genome. </p>
<p>By aligning short-read sequences to the assembled reference sequence and extracting nucleotides that do not match the reference sequence, we have presented a comprehensive profile of HPV genetic variation in the landscape of the whole viral genome. Overall, HPV16 genome sequences in individual cervical specimens showed a highly homogeneous profile. Among the 7 clinical samples, 4 samples did not contain nucleotide substitution with a frequency >0.5% anywhere in their genomes, suggesting that HPV16 genome sequences in infected cells of the cervix are highly static. However, the remaining 3 samples demonstrated nucleotide substitutions (1 to 3 positions) with a frequency >0.5% in the total read coverage (<xref ref-type="fig" rid="pone-0080583-g004">Figure 4</xref>
). The mutations were located in E6, E1 or the non-coding region between E5 and L2, and the levels of the mutation frequency ranged from 0.60 to 5.42%, which is considerably lower than those observed for other viruses, such as HCV and HBV [<xref ref-type="bibr" rid="B13">13</xref>
]. The low frequency of mutation in the HPV genome likely reflects its strategy for cellular replication that relies completely on host DNA polymerases [<xref ref-type="bibr" rid="B39">39</xref>
]. Human DNA polymerases have proof-reading activities that enable faithful replication of the host genome, thus also contributing to the accurate replication of HPV genomes. Nevertheless, we cannot completely exclude the possibility of co-infection by highly related HPV16 variants. In particular, since A at nucleotide position 1842 is a diagnostic SNP to assign the EUR sublineage [<xref ref-type="bibr" rid="B40">40</xref>
], the G to A variation observed at this position in sample 6 (<xref ref-type="table" rid="pone-0080583-t002">Table 2</xref>
) may suggest the coexistence of two closely related variants, As and EUR, in a single specimen. However, no read sequences containing other diagnostic SNPs for EUR (C at nucleotide position 3159, G at nucleotide position 3249, and C at nucleotide position 3787) were detected in this sample.</p>
<p>One possible mechanism for the introduction of low-frequency mutations into the HPV16 genome is an impairment of the p53-dependent DNA damage response (DDR) by the E6 protein. To preserve genetic information when DNA is damaged, eukaryotic cells evoke the DDR involving detection of the DNA damage by sensor proteins, transduction of the signal from DNA lesions to inhibit cell cycle progression and recruitment of repair factors [<xref ref-type="bibr" rid="B41">41</xref>
]. In the course of the DDR, p53 plays important roles in promoting cell survival by activating cell-cycle arrest and DNA repair, or cell death by directing cells to undergo apoptosis thus eliminating extensively damaged cells. To overcome the functions of p53 and create a cellular environment that favors viral propagation, the E6 protein of high-risk HPVs binds to and inactivates p53 via proteasomal degradation [<xref ref-type="bibr" rid="B7">7</xref>
]. In cervical basal cells retaining HPV genomes, or differentiating cells supporting HPV genome amplification, E6 may facilitate the introduction of mutations into both cellular and HPV genomes through inhibition of the p53-dependent DDR. </p>
<p>The mutations found in the E1 gene are non-synonymous substitutions and cause aa changes in the E1 protein, a DNA helicase that plays essential roles in the initiation and progression of viral genome replication. While the amino-acid replacements at positions 294 and 326 located in the DNA-binding domain (see <xref ref-type="fig" rid="pone-0080583-g005">Figure 5A</xref>
), and seen in multiple HPV16 lineages [<xref ref-type="bibr" rid="B40">40</xref>
], did not exert a major effect on the replication activity of E1, the previously unreported mutation at position 381 found in 1 ICC sample led to reduced HPV replication. Since glutamine at 381 is located in the oligomerization domain of E1 (<xref ref-type="fig" rid="pone-0080583-g005">Figure 5A</xref>
) and the nearby Y379F mutant completely lost replication activity (<xref ref-type="fig" rid="pone-0080583-g005">Figure 5C</xref>
), possibly due to a failure to form E1-E1 oligomers [<xref ref-type="bibr" rid="B38">38</xref>
], we suggest that the Q381E mutation may also affect oligomerization of E1. Although amino-acid polymorphisms in E1 have been previously reported for HPV18 and are thought to be responsible for functional differences between HPV18 E1 intra-type variants [<xref ref-type="bibr" rid="B42">42</xref>
], our present study suggests the presence of a mixture of E1 variants in a single clinical specimen as previously reported for HPV62 [<xref ref-type="bibr" rid="B43">43</xref>
]. We have not yet determined whether the HPV16 genomes containing the E1 mutations were episomal or integrated, because full-circle PCR is expected to amplify tandem multiple-copies of the HPV16 genome, that are occasionally found in cervical cancer [<xref ref-type="bibr" rid="B4">4</xref>
], and as we observed with CaSki cells (<xref ref-type="supplementary-material" rid="pone.0080583.s003">Figure S2</xref>
). Frame-shift mutations in the E1 open-reading frame have also been reported for several HPV genotypes including HPV16 [<xref ref-type="bibr" rid="B43">43</xref>
], and we suggest that any loss of E1 activity due to E1-disrupted genomes can be compensated by E1 proteins produced from E1-intact genomes. Although the physiological impact of the mixed state of multiple E1 proteins on HPV life cycle is not clear, somehow the attenuated replication phenotype of the E1 protein, as observed with Q381E, may facilitate the oncogenic progression of HPV-infected lesions, since E1 functions are frequently silenced or lost during progression to cervical cancer. </p>
<p> The 27-bp deletion found in the E2/E4 region in 1 LSIL specimen causes in-frame deletions in both the E2 and E4 proteins, which may have an impact on their functions. The E2 protein is a site-specific DNA-binding protein required for HPV replication, transcription and genome segregation, and consists of an amino-terminal <italic>trans</italic>
-acting domain, a central hinge domain, and a carboxyl-terminal protein dimerization and DNA-binding domain [<xref ref-type="bibr" rid="B44">44</xref>
]. Because the deletion is located in the E2 hinge domain that plays a role in nuclear localization and nuclear matrix association [<xref ref-type="bibr" rid="B45">45</xref>
], it may impair these activities. Intriguingly, deletions in the E2 hinge region have previously been reported in about 80% of cervical carcinoma in situ and invasive cervical cancer specimens [<xref ref-type="bibr" rid="B46">46</xref>
]. Given that the E2 protein may be a transcriptional repressor of the HPV16 promoter responsible for E6/E7 expression, and its functional loss by integration into the host genome may contribute to up-regulation of the E6/E7 expression and development of cervical cancer [<xref ref-type="bibr" rid="B1">1</xref>
], further studies with larger sample sizes will be needed to assess whether the deletion of the hinge domain has any general consequences for the progression of precancerous lesions to ICC.</p>
<p>Although a previous study demonstrated the presence of A/T-rich hypermutation in the LCR of HPV16 in precancerous lesions [<xref ref-type="bibr" rid="B16">16</xref>
], we could not detect such hypermutated sequences with a frequency >0.5% in the total coverage at any positions in the viral genome. This discrepancy is likely due to the extremely high sensitivity of differential DNA denaturation-PCR (3D-PCR) utilized in the previous study to detect hypermutated A/T-rich DNA. To denature target DNA, 3D-PCR uses lower melting temperatures than that used in conventional PCR, which favors selective amplification of A/T-rich DNA and results in enrichment of A/T-rich mutated DNA in the final PCR product. Thus 3D-PCR suggests the presence of hypermutation in the HPV genome, but provides no information about the overall frequency of such hypermutation events. In contrast, we have used an unbiased PCR method to recover the full-length HPV16 genome sequence in a single reaction: a method that is more suitable for capturing the overall frequency of mutation throughout the viral genome. Our analysis of mutation frequency suggests that if such A/T-rich hypermutation were present the frequency is below 0.5%, suggesting that it constitutes only a minor fraction of the viral population. The exact frequency and physiological significance of A/T-rich hypermutation requires further investigation.</p>
<p>Lastly, as shown in this study, deep sequencing of the full-length HPV genome allows comprehensive analyses of those genetic variations, including base substitutions and insertions/deletions, which occur during persistent HPV infection and that may have a link to disease progression. Further investigation of full-length HPV genome sequences in individual clinical samples is essential for a better understanding of the genetic and epidemiological basis for HPV-associated carcinogenesis because some E6 and long control region variants of HPV16 are known to possess enhanced oncogenicity [<xref ref-type="bibr" rid="B47">47</xref>
-<xref ref-type="bibr" rid="B49">49</xref>
] and HPV16 variant lineages show different distributions in cervical cancer cases by country and region [<xref ref-type="bibr" rid="B50">50</xref>
]. Although the level of genetic variation in HPV genomes isolated from cervical specimens observed in this study is seemingly low, its influence on carcinogenesis, diagnostics, and clinical management will undoubtedly be an important topic for future HPV research.</p>
</sec>
<sec sec-type="supplementary-material"><title>Supporting Information</title>
<supplementary-material content-type="local-data" id="pone.0080583.s001"><label>Table S1</label>
<caption><p><bold><italic>De</italic>
<italic>novo</italic>
assembly of complete HPV52/58 genome sequences from short-read sequence data.</bold>
</p>
<p>(DOC)</p>
</caption>
<media xlink:href="pone.0080583.s001.doc"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0080583.s002"><label>Figure S1</label>
<caption><p><bold>Sensitivity and specificity of full-circle PCR.</bold>
Full-circle PCR was performed with PrimeSTAR<sup>®</sup>
GXL DNA polymerase and the primer-pair 1742F/1873R in the presence (lanes 2 to 6) or absence (lanes 8 to 12) of cellular DNA (10 ng per reaction). The amounts of HPV16/pUC19 used for the PCR template are also indicated. M, DNA size markers.</p>
<p>(PDF)</p>
</caption>
<media xlink:href="pone.0080583.s002.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0080583.s003"><label>Figure S2</label>
<caption><p><bold>Amplification of full-length HPV16 genomes from CaSki cells.</bold>
(A) Scheme for PCR amplification from integrated HPV16 genomes in cervical cancer. Tandem multiple-copies of integrated HPV16 DNA can be a template for full-circle PCR as well as episomal HPV16 genomes, while single-copy integrated HPV16 DNA cannot. (B) PCR was performed with primer-pair 1742F/1873R and increasing amounts of total DNA extracted from CaSki cells (lanes 2 and 3). M, DNA size marker (lane 1). </p>
<p>(PDF)</p>
</caption>
<media xlink:href="pone.0080583.s003.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0080583.s004"><label>Figure S3</label>
<caption><p><bold>Read-depth profile obtained by deep sequencing of long PCR products.</bold>
(A) Paired-end read sequences were aligned using BWA to each <italic>de</italic>
<italic>novo</italic>
assembled complete HPV16 genome sequence, and resultant read-depth profiles for CaSki and LSIL sample 2 are shown. Arrows indicate positions of the discontinuous distribution of read depth. The maximum read-depths were 30,173 for CaSki and 23,901 for sample 2. (B) Scheme for the genomic organization of HPV16. </p>
<p>(PDF)</p>
</caption>
<media xlink:href="pone.0080583.s004.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0080583.s005"><label>Figure S4</label>
<caption><p><bold>Alignment of <italic>de</italic>
<italic>novo</italic>
assembled HPV16 genome sequences.</bold>
Eight <italic>de</italic>
<italic>novo</italic>
assembled complete HPV16 genome sequences (samples 1 to 7, and W12) were aligned against each other by MAFFT. The sequence alignment around the non-coding region between E5 and L2 is presented. The stop codon of E5 is indicated with dotted-line box. Numbering of nucleotide positions is based on the sequence of #1.</p>
<p>(PDF)</p>
</caption>
<media xlink:href="pone.0080583.s005.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0080583.s006"><label>Figure S5</label>
<caption><p><bold>Phylogenetic analysis of <italic>de</italic>
<italic>novo</italic>
assembled complete HPV16</bold>
<bold>genome sequences</bold>
. Sixty-two complete HPV16 genome sequences previously reported (Smith et al., PLoS One: 6, e21375, 2011) and 7 complete HPV16 genome sequences <italic>de</italic>
<italic>novo</italic>
assembled in this study were aligned against each other by mafft. Maximum likelihood phylogenetic tree was constructed using the aligned sequences with 1,000-fold bootstrapping. The scale indicates a branch length of 0.01 that shows 1% difference between the nucleotide sequences at the beginning and end of the branch. The number at each branch node represents the bootstrapping value. Vertical bars indicate major clades of HPV16 variants: AFR1, AFR2, NA, AA, As, and EUR.</p>
<p>(PDF)</p>
</caption>
<media xlink:href="pone.0080583.s006.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0080583.s007"><label>Figure S6</label>
<caption><p><bold>Mutation frequency profile of full-length HPV16 genomes in clinical specimens.</bold>
The read sequences obtained with full-length HPV16 genomes prepared from clinical specimens (#1, LSIL; #7, ICC) were aligned to their <italic>de</italic>
<italic>novo</italic>
assembled complete genome sequences, and mutation/error frequencies at each nucleotide position are presented in the landscape of the full-length HPV16 genome. A threshold line for a reliable mutation frequency (0.5%) is indicated with the red dotted line. The genome organization of HPV16 is indicated above.</p>
<p>(PDF)</p>
</caption>
<media xlink:href="pone.0080583.s007.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0080583.s008"><label>Figure S7</label>
<caption><p><bold>Amplification of full-length genomes of HPV52 and HPV58 from clinical specimens.</bold>
Full-circle PCR was performed with HPV52 or HPV58-specific primers and DNA isolated from 8 HPV52-positive LSIL (lanes 1 to 8) and 8 HPV58-positive LSIL (lanes 10 to 17) specimens. Primers are as follows: HPV52-F, <named-content content-type="gene">5’-ACC AGA AAC ACA TAT GGT AAT AGA ACC-3’</named-content>
; HPV52-R, <named-content content-type="gene">5’-GTA ATA CTG TTT GTT GTT CTA TCC ATT C-3’</named-content>
; HPV58-F, <named-content content-type="gene">5’-TAC TAT CAA TTC CTG AAA CAT GTA TGA-3’</named-content>
; HPV58-R, <named-content content-type="gene">5’-AAT CTA TCT ATC CAT TCT GGT GTT G-3’</named-content>
. M, DNA size marker (lane 9).</p>
<p>(PDF)</p>
</caption>
<media xlink:href="pone.0080583.s008.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0080583.s009"><label>Figure S8</label>
<caption><p><bold>Amplification of deleted E2 sequence from a clinical specimen.</bold>
Conventional PCR using AmpliTaq Gold was performed with HPV16-specific primers in the E2 gene and DNA extracted from LSIL specimen 2 (lane 2) and HPV16/pUC19 (lane 3). M, DNA size marker (lane 1).</p>
<p>(PDF)</p>
</caption>
<media xlink:href="pone.0080583.s009.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
</sec>
</body>
<back><ack><p>We thank Drs. Tomomi Nakahara and Paul Lambert for providing us with the W12 cell line. </p>
</ack>
<ref-list><title>References</title>
<ref id="B1"><label>1</label>
<mixed-citation publication-type="journal"><name><surname>zur Hausen</surname>
<given-names>H</given-names>
</name>
(<year>2002</year>
) <article-title>Papillomaviruses and cancer: from basic studies to clinical application</article-title>
. <source>Nat Rev Cancer</source>
<volume>2</volume>
: <fpage>342</fpage>
-<lpage>350</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1038/nrc798">10.1038/nrc798</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/12044010">12044010</ext-link>
<pub-id pub-id-type="pmid">12044010</pub-id>
</mixed-citation>
</ref>
<ref id="B2"><label>2</label>
<mixed-citation publication-type="journal"><name><surname>de Sanjose</surname>
<given-names>S</given-names>
</name>
, <name><surname>Quint</surname>
<given-names>WG</given-names>
</name>
, <name><surname>Alemany</surname>
<given-names>L</given-names>
</name>
, <name><surname>Geraets</surname>
<given-names>DT</given-names>
</name>
, <name><surname>Klaustermeier</surname>
<given-names>JE</given-names>
</name>
<etal>et al.</etal>
(<year>2010</year>
) <article-title>Human papillomavirus genotype attribution in invasive cervical cancer: a retrospective cross-sectional worldwide study</article-title>
. <source>Lancet Oncol</source>
<volume>11</volume>
: <fpage>1048</fpage>
-<lpage>1056</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1016/S1470-2045(10)70230-8">10.1016/S1470-2045(10)70230-8</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/20952254">20952254</ext-link>
<pub-id pub-id-type="pmid">20952254</pub-id>
</mixed-citation>
</ref>
<ref id="B3"><label>3</label>
<mixed-citation publication-type="journal"><name><surname>Kajitani</surname>
<given-names>N</given-names>
</name>
, <name><surname>Satsuka</surname>
<given-names>A</given-names>
</name>
, <name><surname>Kawate</surname>
<given-names>A</given-names>
</name>
, <name><surname>Sakai</surname>
<given-names>H</given-names>
</name>
(<year>2012</year>
) <article-title>Productive Lifecycle of Human Papillomaviruses that Depends Upon Squamous Epithelial Differentiation</article-title>
. <source>Front Microbiol</source>
<volume>3</volume>
: <fpage>152</fpage>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/22536200">22536200</ext-link>
<pub-id pub-id-type="pmid">22536200</pub-id>
</mixed-citation>
</ref>
<ref id="B4"><label>4</label>
<mixed-citation publication-type="journal"><name><surname>Pett</surname>
<given-names>M</given-names>
</name>
, <name><surname>Coleman</surname>
<given-names>N</given-names>
</name>
(<year>2007</year>
) <article-title>Integration of high-risk human papillomavirus: a key event in cervical carcinogenesis?</article-title>
<source>J Pathol</source>
<volume>212</volume>
: <fpage>356</fpage>
-<lpage>367</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1002/path.2192">10.1002/path.2192</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/17573670">17573670</ext-link>
<pub-id pub-id-type="pmid">17573670</pub-id>
</mixed-citation>
</ref>
<ref id="B5"><label>5</label>
<mixed-citation publication-type="journal"><name><surname>Xu</surname>
<given-names>B</given-names>
</name>
, <name><surname>Chotewutmontri</surname>
<given-names>S</given-names>
</name>
, <name><surname>Wolf</surname>
<given-names>S</given-names>
</name>
, <name><surname>Klos</surname>
<given-names>U</given-names>
</name>
, <name><surname>Schmitz</surname>
<given-names>M</given-names>
</name>
<etal>et al.</etal>
(<year>2013</year>
) <article-title>Multiplex Identification of Human Papillomavirus 16 DNA Integration Sites in Cervical Carcinomas</article-title>
. <source>PLOS ONE</source>
<volume>8</volume>
: <fpage>e66693</fpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1371/journal.pone.0066693">10.1371/journal.pone.0066693</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/23824673">23824673</ext-link>
<pub-id pub-id-type="pmid">23824673</pub-id>
</mixed-citation>
</ref>
<ref id="B6"><label>6</label>
<mixed-citation publication-type="journal"><name><surname>Hebner</surname>
<given-names>CM</given-names>
</name>
, <name><surname>Laimins</surname>
<given-names>LA</given-names>
</name>
(<year>2006</year>
) <article-title>Human papillomaviruses: basic mechanisms of pathogenesis and oncogenicity</article-title>
. <source>Rev Med Virol</source>
<volume>16</volume>
: <fpage>83</fpage>
-<lpage>97</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1002/rmv.488">10.1002/rmv.488</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/16287204">16287204</ext-link>
<pub-id pub-id-type="pmid">16287204</pub-id>
</mixed-citation>
</ref>
<ref id="B7"><label>7</label>
<mixed-citation publication-type="journal"><name><surname>Bodily</surname>
<given-names>J</given-names>
</name>
, <name><surname>Laimins</surname>
<given-names>LA</given-names>
</name>
(<year>2011</year>
) <article-title>Persistence of human papillomavirus infection: keys to malignant progression</article-title>
. <source>Trends Microbiol</source>
<volume>19</volume>
: <fpage>33</fpage>
-<lpage>39</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1016/j.tim.2010.10.002">10.1016/j.tim.2010.10.002</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/21050765">21050765</ext-link>
<pub-id pub-id-type="pmid">21050765</pub-id>
</mixed-citation>
</ref>
<ref id="B8"><label>8</label>
<mixed-citation publication-type="journal"><name><surname>Beerenwinkel</surname>
<given-names>N</given-names>
</name>
, <name><surname>Gunthard</surname>
<given-names>HF</given-names>
</name>
, <name><surname>Roth</surname>
<given-names>V</given-names>
</name>
, <name><surname>Metzner</surname>
<given-names>KJ</given-names>
</name>
(<year>2012</year>
) <article-title>Challenges and opportunities in estimating viral genetic diversity from next-generation sequencing data</article-title>
. <source>Front Microbiol</source>
<volume>3</volume>
: <fpage>329</fpage>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/22973268">22973268</ext-link>
<pub-id pub-id-type="pmid">22973268</pub-id>
</mixed-citation>
</ref>
<ref id="B9"><label>9</label>
<mixed-citation publication-type="journal"><name><surname>Ekström</surname>
<given-names>J</given-names>
</name>
, <name><surname>Bzhalava</surname>
<given-names>D</given-names>
</name>
, <name><surname>Svenback</surname>
<given-names>D</given-names>
</name>
, <name><surname>Forslund</surname>
<given-names>O</given-names>
</name>
, <name><surname>Dillner</surname>
<given-names>J</given-names>
</name>
(<year>2011</year>
) <article-title>High throughput sequencing reveals diversity of Human Papillomaviruses in cutaneous lesions</article-title>
. <source>Int J Cancer</source>
<volume>129</volume>
: <fpage>2643</fpage>
-<lpage>2650</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1002/ijc.26204">10.1002/ijc.26204</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/21630257">21630257</ext-link>
<pub-id pub-id-type="pmid">21630257</pub-id>
</mixed-citation>
</ref>
<ref id="B10"><label>10</label>
<mixed-citation publication-type="journal"><name><surname>Meiring</surname>
<given-names>TL</given-names>
</name>
, <name><surname>Salimo</surname>
<given-names>AT</given-names>
</name>
, <name><surname>Coetzee</surname>
<given-names>B</given-names>
</name>
, <name><surname>Maree</surname>
<given-names>HJ</given-names>
</name>
, <name><surname>Moodley</surname>
<given-names>J</given-names>
</name>
<etal>et al.</etal>
(<year>2012</year>
) <article-title>Next-generation sequencing of cervical DNA detects human papillomavirus types not detected by commercial kits</article-title>
. <source>Virol J</source>
<volume>9</volume>
: <fpage>164</fpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1186/1743-422X-9-164">10.1186/1743-422X-9-164</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/22897914">22897914</ext-link>
<pub-id pub-id-type="pmid">22897914</pub-id>
</mixed-citation>
</ref>
<ref id="B11"><label>11</label>
<mixed-citation publication-type="journal"><name><surname>Foulongne</surname>
<given-names>V</given-names>
</name>
, <name><surname>Sauvage</surname>
<given-names>V</given-names>
</name>
, <name><surname>Hebert</surname>
<given-names>C</given-names>
</name>
, <name><surname>Dereure</surname>
<given-names>O</given-names>
</name>
, <name><surname>Cheval</surname>
<given-names>J</given-names>
</name>
<etal>et al.</etal>
(<year>2012</year>
) <article-title>Human skin microbiota: high diversity of DNA viruses identified on the human skin by high throughput sequencing</article-title>
. <source>PLOS ONE</source>
<volume>7</volume>
: <fpage>e38499</fpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1371/journal.pone.0038499">10.1371/journal.pone.0038499</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/22723863">22723863</ext-link>
<pub-id pub-id-type="pmid">22723863</pub-id>
</mixed-citation>
</ref>
<ref id="B12"><label>12</label>
<mixed-citation publication-type="journal"><name><surname>Lauring</surname>
<given-names>AS</given-names>
</name>
, <name><surname>Andino</surname>
<given-names>R</given-names>
</name>
(<year>2010</year>
) <article-title>Quasispecies theory and the behavior of RNA viruses</article-title>
. <source>PLOS Pathog</source>
<volume>6</volume>
: <fpage>e1001005</fpage>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/20661479">20661479</ext-link>
<pub-id pub-id-type="pmid">20661479</pub-id>
</mixed-citation>
</ref>
<ref id="B13"><label>13</label>
<mixed-citation publication-type="journal"><name><surname>Domingo</surname>
<given-names>E</given-names>
</name>
, <name><surname>Sheldon</surname>
<given-names>J</given-names>
</name>
, <name><surname>Perales</surname>
<given-names>C</given-names>
</name>
(<year>2012</year>
) <article-title>Viral quasispecies evolution</article-title>
. <source>Microbiol Mol Biol Rev</source>
<volume>76</volume>
: <fpage>159</fpage>
-<lpage>216</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1128/MMBR.05023-11">10.1128/MMBR.05023-11</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/22688811">22688811</ext-link>
<pub-id pub-id-type="pmid">22688811</pub-id>
</mixed-citation>
</ref>
<ref id="B14"><label>14</label>
<mixed-citation publication-type="journal"><name><surname>van Velzen</surname>
<given-names>M</given-names>
</name>
, <name><surname>van Loenen</surname>
<given-names>FB</given-names>
</name>
, <name><surname>Meesters</surname>
<given-names>RJ</given-names>
</name>
, <name><surname>de Graaf</surname>
<given-names>M</given-names>
</name>
, <name><surname>Remeijer</surname>
<given-names>L</given-names>
</name>
<etal>et al.</etal>
(<year>2012</year>
) <article-title>Latent acyclovir-resistant herpes simplex virus type 1 in trigeminal ganglia of immunocompetent individuals</article-title>
. <source>J Infect Dis</source>
<volume>205</volume>
: <fpage>1539</fpage>
-<lpage>1543</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1093/infdis/jis237">10.1093/infdis/jis237</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/22457282">22457282</ext-link>
<pub-id pub-id-type="pmid">22457282</pub-id>
</mixed-citation>
</ref>
<ref id="B15"><label>15</label>
<mixed-citation publication-type="journal"><name><surname>Kwok</surname>
<given-names>H</given-names>
</name>
, <name><surname>Tong</surname>
<given-names>AH</given-names>
</name>
, <name><surname>Lin</surname>
<given-names>CH</given-names>
</name>
, <name><surname>Lok</surname>
<given-names>S</given-names>
</name>
, <name><surname>Farrell</surname>
<given-names>PJ</given-names>
</name>
<etal>et al.</etal>
(<year>2012</year>
) <article-title>Genomic sequencing and comparative analysis of Epstein-Barr virus genome isolated from primary nasopharyngeal carcinoma biopsy</article-title>
. <source>PLOS ONE</source>
<volume>7</volume>
: <fpage>e36939</fpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1371/journal.pone.0036939">10.1371/journal.pone.0036939</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/22590638">22590638</ext-link>
<pub-id pub-id-type="pmid">22590638</pub-id>
</mixed-citation>
</ref>
<ref id="B16"><label>16</label>
<mixed-citation publication-type="journal"><name><surname>Vartanian</surname>
<given-names>JP</given-names>
</name>
, <name><surname>Guétard</surname>
<given-names>D</given-names>
</name>
, <name><surname>Henry</surname>
<given-names>M</given-names>
</name>
, <name><surname>Wain-Hobson</surname>
<given-names>S</given-names>
</name>
(<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>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1126/science.1153201">10.1126/science.1153201</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/18403710">18403710</ext-link>
<pub-id pub-id-type="pmid">18403710</pub-id>
</mixed-citation>
</ref>
<ref id="B17"><label>17</label>
<mixed-citation publication-type="journal"><name><surname>Jeon</surname>
<given-names>S</given-names>
</name>
, <name><surname>Allen-Hoffmann</surname>
<given-names>BL</given-names>
</name>
, <name><surname>Lambert</surname>
<given-names>PF</given-names>
</name>
(<year>1995</year>
) <article-title>Integration of human papillomavirus type 16 into the human genome correlates with a selective growth advantage of cells</article-title>
. <source>J Virol</source>
<volume>69</volume>
: <fpage>2989</fpage>
-<lpage>2997</lpage>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/7707525">7707525</ext-link>
<pub-id pub-id-type="pmid">7707525</pub-id>
</mixed-citation>
</ref>
<ref id="B18"><label>18</label>
<mixed-citation publication-type="journal"><name><surname>Stanley</surname>
<given-names>MA</given-names>
</name>
, <name><surname>Browne</surname>
<given-names>HM</given-names>
</name>
, <name><surname>Appleby</surname>
<given-names>M</given-names>
</name>
, <name><surname>Minson</surname>
<given-names>AC</given-names>
</name>
(<year>1989</year>
) <article-title>Properties of a non-tumorigenic human cervical keratinocyte cell line</article-title>
. <source>Int J Cancer</source>
<volume>43</volume>
: <fpage>672</fpage>
-<lpage>676</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1002/ijc.2910430422">10.1002/ijc.2910430422</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/2467886">2467886</ext-link>
<pub-id pub-id-type="pmid">2467886</pub-id>
</mixed-citation>
</ref>
<ref id="B19"><label>19</label>
<mixed-citation publication-type="journal"><name><surname>Arad</surname>
<given-names>U</given-names>
</name>
(<year>1998</year>
) <article-title>Modified Hirt procedure for rapid purification of extrachromosomal DNA from mammalian cells</article-title>
. <source>BioTechniques</source>
<volume>24</volume>
: <fpage>760</fpage>
-<lpage>762</lpage>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/9591124">9591124</ext-link>
<pub-id pub-id-type="pmid">9591124</pub-id>
</mixed-citation>
</ref>
<ref id="B20"><label>20</label>
<mixed-citation publication-type="book"><collab>WHO</collab>
(<year>2009</year>
) <source>HUMAN PAPILLOMAVIRUS LABORATORY MANUAL</source>
, <edition>1st</edition>
edition. <comment>(WNO/IVB/10.12) WHO, Geneva (Switzerland), WHO website. Available: </comment>
<ext-link ext-link-type="uri" xlink:href="http://www.who.int/immunization/documents/who_ivb_10.12/en/index.html">http://www.who.int/immunization/documents/WHO_IVB_10.12/en/index.html</ext-link>
<comment>. Accessed 2013 Oct 16.</comment>
</mixed-citation>
</ref>
<ref id="B21"><label>21</label>
<mixed-citation publication-type="journal"><name><surname>Estrade</surname>
<given-names>C</given-names>
</name>
, <name><surname>Menoud</surname>
<given-names>PA</given-names>
</name>
, <name><surname>Nardelli-Haefliger</surname>
<given-names>D</given-names>
</name>
, <name><surname>Sahli</surname>
<given-names>R</given-names>
</name>
(<year>2011</year>
) <article-title>Validation of a low-cost human papillomavirus genotyping assay based on PGMY PCR and reverse blotting hybridization with reusable membranes</article-title>
. <source>J Clin Microbiol</source>
<volume>49</volume>
: <fpage>3474</fpage>
-<lpage>3481</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1128/JCM.05039-11">10.1128/JCM.05039-11</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/21832011">21832011</ext-link>
<pub-id pub-id-type="pmid">21832011</pub-id>
</mixed-citation>
</ref>
<ref id="B22"><label>22</label>
<mixed-citation publication-type="journal"><name><surname>Kanda</surname>
<given-names>T</given-names>
</name>
, <name><surname>Watanabe</surname>
<given-names>S</given-names>
</name>
, <name><surname>Yoshiike</surname>
<given-names>K</given-names>
</name>
(<year>1987</year>
) <article-title>Human papillomavirus type 16 transformation of rat 3Y1 cells</article-title>
. <source>Jpn J Cancer Res</source>
<volume>78</volume>
: <fpage>103</fpage>
-<lpage>108</lpage>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/3030985">3030985</ext-link>
<pub-id pub-id-type="pmid">3030985</pub-id>
</mixed-citation>
</ref>
<ref id="B23"><label>23</label>
<mixed-citation publication-type="journal"><name><surname>Simpson</surname>
<given-names>JT</given-names>
</name>
, <name><surname>Wong</surname>
<given-names>K</given-names>
</name>
, <name><surname>Jackman</surname>
<given-names>SD</given-names>
</name>
, <name><surname>Schein</surname>
<given-names>JE</given-names>
</name>
, <name><surname>Jones</surname>
<given-names>SJ</given-names>
</name>
<etal>et al.</etal>
(<year>2009</year>
) <article-title>ABySS: a parallel assembler for short read sequence data</article-title>
. <source>Genome Res</source>
<volume>19</volume>
: <fpage>1117</fpage>
-<lpage>1123</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1101/gr.089532.108">10.1101/gr.089532.108</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/19251739">19251739</ext-link>
<pub-id pub-id-type="pmid">19251739</pub-id>
</mixed-citation>
</ref>
<ref id="B24"><label>24</label>
<mixed-citation publication-type="journal"><name><surname>Gordon</surname>
<given-names>D</given-names>
</name>
, <name><surname>Abajian</surname>
<given-names>C</given-names>
</name>
, <name><surname>Green</surname>
<given-names>P</given-names>
</name>
(<year>1998</year>
) <article-title>Consed: a graphical tool for sequence finishing</article-title>
. <source>Genome Res</source>
<volume>8</volume>
: <fpage>195</fpage>
-<lpage>202</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1101/gr.8.3.195">10.1101/gr.8.3.195</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/9521923">9521923</ext-link>
<pub-id pub-id-type="pmid">9521923</pub-id>
</mixed-citation>
</ref>
<ref id="B25"><label>25</label>
<mixed-citation publication-type="journal"><name><surname>Li</surname>
<given-names>H</given-names>
</name>
, <name><surname>Durbin</surname>
<given-names>R</given-names>
</name>
(<year>2010</year>
) <article-title>Fast and accurate long-read alignment with Burrows-Wheeler transform</article-title>
. <source>Bioinformatics</source>
<volume>26</volume>
: <fpage>589</fpage>
-<lpage>595</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1093/bioinformatics/btp698">10.1093/bioinformatics/btp698</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/20080505">20080505</ext-link>
<pub-id pub-id-type="pmid">20080505</pub-id>
</mixed-citation>
</ref>
<ref id="B26"><label>26</label>
<mixed-citation publication-type="journal"><name><surname>Futschik</surname>
<given-names>A</given-names>
</name>
, <name><surname>Schlötterer</surname>
<given-names>C</given-names>
</name>
(<year>2010</year>
) <article-title>The next generation of molecular markers from massively parallel sequencing of pooled DNA samples</article-title>
. <source>Genetics</source>
<volume>186</volume>
: <fpage>207</fpage>
-<lpage>218</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1534/genetics.110.114397">10.1534/genetics.110.114397</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/20457880">20457880</ext-link>
<pub-id pub-id-type="pmid">20457880</pub-id>
</mixed-citation>
</ref>
<ref id="B27"><label>27</label>
<mixed-citation publication-type="journal"><name><surname>Thorvaldsdóttir</surname>
<given-names>H</given-names>
</name>
, <name><surname>Robinson</surname>
<given-names>JT</given-names>
</name>
, <name><surname>Mesirov</surname>
<given-names>JP</given-names>
</name>
(<year>2013</year>
) <article-title>Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration</article-title>
. <source>Brief Bioinform</source>
<volume>14</volume>
: <fpage>178</fpage>
-<lpage>192</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1093/bib/bbs017">10.1093/bib/bbs017</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/22517427">22517427</ext-link>
<pub-id pub-id-type="pmid">22517427</pub-id>
</mixed-citation>
</ref>
<ref id="B28"><label>28</label>
<mixed-citation publication-type="journal"><name><surname>Milne</surname>
<given-names>I</given-names>
</name>
, <name><surname>Bayer</surname>
<given-names>M</given-names>
</name>
, <name><surname>Cardle</surname>
<given-names>L</given-names>
</name>
, <name><surname>Shaw</surname>
<given-names>P</given-names>
</name>
, <name><surname>Stephen</surname>
<given-names>G</given-names>
</name>
<etal>et al.</etal>
(<year>2010</year>
) <article-title>Tablet--next generation sequence assembly visualization</article-title>
. <source>Bioinformatics</source>
<volume>26</volume>
: <fpage>401</fpage>
-<lpage>402</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1093/bioinformatics/btp666">10.1093/bioinformatics/btp666</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/19965881">19965881</ext-link>
<pub-id pub-id-type="pmid">19965881</pub-id>
</mixed-citation>
</ref>
<ref id="B29"><label>29</label>
<mixed-citation publication-type="journal"><name><surname>Li</surname>
<given-names>H</given-names>
</name>
, <name><surname>Handsaker</surname>
<given-names>B</given-names>
</name>
, <name><surname>Wysoker</surname>
<given-names>A</given-names>
</name>
, <name><surname>Fennell</surname>
<given-names>T</given-names>
</name>
, <name><surname>Ruan</surname>
<given-names>J</given-names>
</name>
<etal>et al.</etal>
(<year>2009</year>
) <article-title>The Sequence Alignment/Map format and SAMtools</article-title>
. <source>Bioinformatics</source>
<volume>25</volume>
: <fpage>2078</fpage>
-<lpage>2079</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1093/bioinformatics/btp352">10.1093/bioinformatics/btp352</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/19505943">19505943</ext-link>
<pub-id pub-id-type="pmid">19505943</pub-id>
</mixed-citation>
</ref>
<ref id="B30"><label>30</label>
<mixed-citation publication-type="journal"><name><surname>Koboldt</surname>
<given-names>DC</given-names>
</name>
, <name><surname>Chen</surname>
<given-names>K</given-names>
</name>
, <name><surname>Wylie</surname>
<given-names>T</given-names>
</name>
, <name><surname>Larson</surname>
<given-names>DE</given-names>
</name>
, <name><surname>McLellan</surname>
<given-names>MD</given-names>
</name>
<etal>et al.</etal>
(<year>2009</year>
) <article-title>VarScan: variant detection in massively parallel sequencing of individual and pooled samples</article-title>
. <source>Bioinformatics</source>
<volume>25</volume>
: <fpage>2283</fpage>
-<lpage>2285</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1093/bioinformatics/btp373">10.1093/bioinformatics/btp373</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/19542151">19542151</ext-link>
<pub-id pub-id-type="pmid">19542151</pub-id>
</mixed-citation>
</ref>
<ref id="B31"><label>31</label>
<mixed-citation publication-type="journal"><name><surname>Fradet-Turcotte</surname>
<given-names>A</given-names>
</name>
, <name><surname>Morin</surname>
<given-names>G</given-names>
</name>
, <name><surname>Lehoux</surname>
<given-names>M</given-names>
</name>
, <name><surname>Bullock</surname>
<given-names>PA</given-names>
</name>
, <name><surname>Archambault</surname>
<given-names>J</given-names>
</name>
(<year>2010</year>
) <article-title>Development of quantitative and high-throughput assays of polyomavirus and papillomavirus DNA replication</article-title>
. <source>Virology</source>
<volume>399</volume>
: <fpage>65</fpage>
-<lpage>76</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1016/j.virol.2009.12.026">10.1016/j.virol.2009.12.026</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/20079917">20079917</ext-link>
<pub-id pub-id-type="pmid">20079917</pub-id>
</mixed-citation>
</ref>
<ref id="B32"><label>32</label>
<mixed-citation publication-type="journal"><name><surname>Yee</surname>
<given-names>C</given-names>
</name>
, <name><surname>Krishnan-Hewlett</surname>
<given-names>I</given-names>
</name>
, <name><surname>Baker</surname>
<given-names>CC</given-names>
</name>
, <name><surname>Schlegel</surname>
<given-names>R</given-names>
</name>
, <name><surname>Howley</surname>
<given-names>PM</given-names>
</name>
(<year>1985</year>
) <article-title>Presence and expression of human papillomavirus sequences in human cervical carcinoma cell lines</article-title>
. <source>Am J Pathol</source>
<volume>119</volume>
: <fpage>361</fpage>
-<lpage>366</lpage>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/2990217">2990217</ext-link>
<pub-id pub-id-type="pmid">2990217</pub-id>
</mixed-citation>
</ref>
<ref id="B33"><label>33</label>
<mixed-citation publication-type="journal"><name><surname>Meissner</surname>
<given-names>JD</given-names>
</name>
(<year>1999</year>
) <article-title>Nucleotide sequences and further characterization of human papillomavirus DNA present in the CaSki, SiHa and HeLa cervical carcinoma cell lines</article-title>
. <source>J Gen Virol</source>
<volume>80</volume>
( <issue-id>7</issue-id>
): <fpage>1725</fpage>
-<lpage>1733</lpage>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/10423141">10423141</ext-link>
<pub-id pub-id-type="pmid">10423141</pub-id>
</mixed-citation>
</ref>
<ref id="B34"><label>34</label>
<mixed-citation publication-type="journal"><name><surname>Cornet</surname>
<given-names>I</given-names>
</name>
, <name><surname>Gheit</surname>
<given-names>T</given-names>
</name>
, <name><surname>Franceschi</surname>
<given-names>S</given-names>
</name>
, <name><surname>Vignat</surname>
<given-names>J</given-names>
</name>
, <name><surname>Burk</surname>
<given-names>RD</given-names>
</name>
<etal>et al.</etal>
(<year>2012</year>
) <article-title>Human Papillomavirus Type 16 Genetic Variants: Phylogeny and Classification Based on E6 and LCR</article-title>
. <source>J Virol</source>
<volume>86</volume>
: <fpage>6855</fpage>
-<lpage>6861</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1128/JVI.00483-12">10.1128/JVI.00483-12</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/22491459">22491459</ext-link>
<pub-id pub-id-type="pmid">22491459</pub-id>
</mixed-citation>
</ref>
<ref id="B35"><label>35</label>
<mixed-citation publication-type="journal"><name><surname>Lin</surname>
<given-names>QQ</given-names>
</name>
, <name><surname>Yu</surname>
<given-names>SZ</given-names>
</name>
, <name><surname>Qu</surname>
<given-names>W</given-names>
</name>
, <name><surname>Cruz</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Burk</surname>
<given-names>RD</given-names>
</name>
(<year>1998</year>
) <article-title>Human papillomavirus types 52 and 58</article-title>
. <source>Int J Cancer</source>
<volume>75</volume>
: <fpage>484</fpage>
-<lpage>485</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1002/(SICI)1097-0215(19980130)75:3">10.1002/(SICI)1097-0215(19980130)75:3</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/9455813">9455813</ext-link>
<pub-id pub-id-type="pmid">9455813</pub-id>
</mixed-citation>
</ref>
<ref id="B36"><label>36</label>
<mixed-citation publication-type="journal"><name><surname>Miura</surname>
<given-names>S</given-names>
</name>
, <name><surname>Matsumoto</surname>
<given-names>K</given-names>
</name>
, <name><surname>Oki</surname>
<given-names>A</given-names>
</name>
, <name><surname>Satoh</surname>
<given-names>T</given-names>
</name>
, <name><surname>Tsunoda</surname>
<given-names>H</given-names>
</name>
<etal>et al.</etal>
(<year>2006</year>
) <article-title>Do we need a different strategy for HPV screening and vaccination in East Asia?</article-title>
<source>Int J Cancer</source>
<volume>119</volume>
: <fpage>2713</fpage>
-<lpage>2715</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1002/ijc.22195">10.1002/ijc.22195</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/16929495">16929495</ext-link>
<pub-id pub-id-type="pmid">16929495</pub-id>
</mixed-citation>
</ref>
<ref id="B37"><label>37</label>
<mixed-citation publication-type="journal"><name><surname>Kondo</surname>
<given-names>K</given-names>
</name>
, <name><surname>Uenoyama</surname>
<given-names>A</given-names>
</name>
, <name><surname>Kitagawa</surname>
<given-names>R</given-names>
</name>
, <name><surname>Tsunoda</surname>
<given-names>H</given-names>
</name>
, <name><surname>Kusumoto-Matsuo</surname>
<given-names>R</given-names>
</name>
<etal>et al.</etal>
(<year>2012</year>
) <article-title>Genotype distribution of human papillomaviruses in Japanese women with abnormal cervical cytology. Open</article-title>
. <source>Virol J</source>
<volume>6</volume>
: <fpage>277</fpage>
-<lpage>283</lpage>
</mixed-citation>
</ref>
<ref id="B38"><label>38</label>
<mixed-citation publication-type="journal"><name><surname>Titolo</surname>
<given-names>S</given-names>
</name>
, <name><surname>Pelletier</surname>
<given-names>A</given-names>
</name>
, <name><surname>Pulichino</surname>
<given-names>AM</given-names>
</name>
, <name><surname>Brault</surname>
<given-names>K</given-names>
</name>
, <name><surname>Wardrop</surname>
<given-names>E</given-names>
</name>
<etal>et al.</etal>
(<year>2000</year>
) <article-title>Identification of domains of the human papillomavirus type 11 E1 helicase involved in oligomerization and binding to the viral origin</article-title>
. <source>J Virol</source>
<volume>74</volume>
: <fpage>7349</fpage>
-<lpage>7361</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1128/JVI.74.16.7349-7361.2000">10.1128/JVI.74.16.7349-7361.2000</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/10906188">10906188</ext-link>
<pub-id pub-id-type="pmid">10906188</pub-id>
</mixed-citation>
</ref>
<ref id="B39"><label>39</label>
<mixed-citation publication-type="journal"><name><surname>Melendy</surname>
<given-names>T</given-names>
</name>
, <name><surname>Sedman</surname>
<given-names>J</given-names>
</name>
, <name><surname>Stenlund</surname>
<given-names>A</given-names>
</name>
(<year>1995</year>
) <article-title>Cellular factors required for papillomavirus DNA replication</article-title>
. <source>J Virol</source>
<volume>69</volume>
: <fpage>7857</fpage>
-<lpage>7867</lpage>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/7494298">7494298</ext-link>
<pub-id pub-id-type="pmid">7494298</pub-id>
</mixed-citation>
</ref>
<ref id="B40"><label>40</label>
<mixed-citation publication-type="journal"><name><surname>Smith</surname>
<given-names>B</given-names>
</name>
, <name><surname>Chen</surname>
<given-names>Z</given-names>
</name>
, <name><surname>Reimers</surname>
<given-names>L</given-names>
</name>
, <name><surname>van Doorslaer</surname>
<given-names>K</given-names>
</name>
, <name><surname>Schiffman</surname>
<given-names>M</given-names>
</name>
<etal>et al.</etal>
(<year>2011</year>
) <article-title>Sequence imputation of HPV16 genomes for genetic association studies</article-title>
. <source>PLOS ONE</source>
<volume>6</volume>
: <fpage>e21375</fpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1371/journal.pone.0021375">10.1371/journal.pone.0021375</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/21731721">21731721</ext-link>
<pub-id pub-id-type="pmid">21731721</pub-id>
</mixed-citation>
</ref>
<ref id="B41"><label>41</label>
<mixed-citation publication-type="journal"><name><surname>Jackson</surname>
<given-names>SP</given-names>
</name>
, <name><surname>Bartek</surname>
<given-names>J</given-names>
</name>
(<year>2009</year>
) <article-title>The DNA-damage response in human biology and disease</article-title>
. <source>Nature</source>
<volume>461</volume>
: <fpage>1071</fpage>
-<lpage>1078</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1038/nature08467">10.1038/nature08467</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/19847258">19847258</ext-link>
<pub-id pub-id-type="pmid">19847258</pub-id>
</mixed-citation>
</ref>
<ref id="B42"><label>42</label>
<mixed-citation publication-type="journal"><name><surname>Amador-Molina</surname>
<given-names>A</given-names>
</name>
, <name><surname>González-Montoya</surname>
<given-names>JL</given-names>
</name>
, <name><surname>García-Carrancá</surname>
<given-names>A</given-names>
</name>
, <name><surname>Mohar</surname>
<given-names>A</given-names>
</name>
, <name><surname>Lizano</surname>
<given-names>M</given-names>
</name>
(<year>2013</year>
) <article-title>Intratypic changes of the E1 gene and the long control region affect ori function of human papillomavirus type 18 variants</article-title>
. <source>J Gen Virol</source>
<volume>94</volume>
: <fpage>393</fpage>
-<lpage>402</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1099/vir.0.045807-0">10.1099/vir.0.045807-0</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/23100366">23100366</ext-link>
<pub-id pub-id-type="pmid">23100366</pub-id>
</mixed-citation>
</ref>
<ref id="B43"><label>43</label>
<mixed-citation publication-type="journal"><name><surname>Fu</surname>
<given-names>L</given-names>
</name>
, <name><surname>Terai</surname>
<given-names>M</given-names>
</name>
, <name><surname>Matsukura</surname>
<given-names>T</given-names>
</name>
, <name><surname>Herrero</surname>
<given-names>R</given-names>
</name>
, <name><surname>Burk</surname>
<given-names>RD</given-names>
</name>
(<year>2004</year>
) <article-title>Codetection of a mixed population of candHPV62 containing wild-type and disrupted E1 open-reading frame in a 45-year-old woman with normal cytology</article-title>
. <source>J Infect Dis</source>
<volume>190</volume>
: <fpage>1303</fpage>
-<lpage>1309</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1086/423855">10.1086/423855</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/15346342">15346342</ext-link>
<pub-id pub-id-type="pmid">15346342</pub-id>
</mixed-citation>
</ref>
<ref id="B44"><label>44</label>
<mixed-citation publication-type="journal"><name><surname>McBride</surname>
<given-names>AA</given-names>
</name>
(<year>2008</year>
) <article-title>Replication and partitioning of papillomavirus genomes</article-title>
. <source>Adv Virus Res</source>
<volume>72</volume>
: <fpage>155</fpage>
-<lpage>205</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1016/S0065-3527(08)00404-1">10.1016/S0065-3527(08)00404-1</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/19081491">19081491</ext-link>
<pub-id pub-id-type="pmid">19081491</pub-id>
</mixed-citation>
</ref>
<ref id="B45"><label>45</label>
<mixed-citation publication-type="journal"><name><surname>Zou</surname>
<given-names>N</given-names>
</name>
, <name><surname>Lin</surname>
<given-names>BY</given-names>
</name>
, <name><surname>Duan</surname>
<given-names>F</given-names>
</name>
, <name><surname>Lee</surname>
<given-names>KY</given-names>
</name>
, <name><surname>Jin</surname>
<given-names>G</given-names>
</name>
<etal>et al.</etal>
(<year>2000</year>
) <article-title>The hinge of the human papillomavirus type 11 E2 protein contains major determinants for nuclear localization and nuclear matrix association</article-title>
. <source>J Virol</source>
<volume>74</volume>
: <fpage>3761</fpage>
-<lpage>3770</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1128/JVI.74.8.3761-3770.2000">10.1128/JVI.74.8.3761-3770.2000</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/10729151">10729151</ext-link>
<pub-id pub-id-type="pmid">10729151</pub-id>
</mixed-citation>
</ref>
<ref id="B46"><label>46</label>
<mixed-citation publication-type="journal"><name><surname>Arias-Pulido</surname>
<given-names>H</given-names>
</name>
, <name><surname>Peyton</surname>
<given-names>CL</given-names>
</name>
, <name><surname>Joste</surname>
<given-names>NE</given-names>
</name>
, <name><surname>Vargas</surname>
<given-names>H</given-names>
</name>
, <name><surname>Wheeler</surname>
<given-names>CM</given-names>
</name>
(<year>2006</year>
) <article-title>Human papillomavirus type 16 integration in cervical carcinoma in situ and in invasive cervical cancer</article-title>
. <source>J Clin Microbiol</source>
<volume>44</volume>
: <fpage>1755</fpage>
-<lpage>1762</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1128/JCM.44.5.1755-1762.2006">10.1128/JCM.44.5.1755-1762.2006</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/16672403">16672403</ext-link>
<pub-id pub-id-type="pmid">16672403</pub-id>
</mixed-citation>
</ref>
<ref id="B47"><label>47</label>
<mixed-citation publication-type="journal"><name><surname>Matsumoto</surname>
<given-names>K</given-names>
</name>
, <name><surname>Yoshikawa</surname>
<given-names>H</given-names>
</name>
, <name><surname>Nakagawa</surname>
<given-names>S</given-names>
</name>
, <name><surname>Tang</surname>
<given-names>X</given-names>
</name>
, <name><surname>Yasugi</surname>
<given-names>T</given-names>
</name>
<etal>et al.</etal>
(<year>2000</year>
) <article-title>Enhanced oncogenicity of human papillomavirus type 16 (HPV16) variants in Japanese population</article-title>
. <source>Cancer Lett</source>
<volume>156</volume>
: <fpage>159</fpage>
-<lpage>165</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1016/S0304-3835(00)00457-2">10.1016/S0304-3835(00)00457-2</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/10880765">10880765</ext-link>
<pub-id pub-id-type="pmid">10880765</pub-id>
</mixed-citation>
</ref>
<ref id="B48"><label>48</label>
<mixed-citation publication-type="journal"><name><surname>Pande</surname>
<given-names>S</given-names>
</name>
, <name><surname>Jain</surname>
<given-names>N</given-names>
</name>
, <name><surname>Prusty</surname>
<given-names>BK</given-names>
</name>
, <name><surname>Bhambhani</surname>
<given-names>S</given-names>
</name>
, <name><surname>Gupta</surname>
<given-names>S</given-names>
</name>
<etal>et al.</etal>
(<year>2008</year>
) <article-title>Human papillomavirus type 16 variant analysis of E6, E7, and L1 genes and long control region in biopsy samples from cervical cancer patients in north India</article-title>
. <source>J Clin Microbiol</source>
<volume>46</volume>
: <fpage>1060</fpage>
-<lpage>1066</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1128/JCM.02202-07">10.1128/JCM.02202-07</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/18199779">18199779</ext-link>
<pub-id pub-id-type="pmid">18199779</pub-id>
</mixed-citation>
</ref>
<ref id="B49"><label>49</label>
<mixed-citation publication-type="journal"><name><surname>Sichero</surname>
<given-names>L</given-names>
</name>
, <name><surname>Sobrinho</surname>
<given-names>JS</given-names>
</name>
, <name><surname>Villa</surname>
<given-names>LL</given-names>
</name>
(<year>2012</year>
) <article-title>Oncogenic potential diverge among human papillomavirus type 16 natural variants</article-title>
. <source>Virology</source>
<volume>432</volume>
: <fpage>127</fpage>
-<lpage>132</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1016/j.virol.2012.06.011">10.1016/j.virol.2012.06.011</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/22748253">22748253</ext-link>
<pub-id pub-id-type="pmid">22748253</pub-id>
</mixed-citation>
</ref>
<ref id="B50"><label>50</label>
<mixed-citation publication-type="journal"><name><surname>Cornet</surname>
<given-names>I</given-names>
</name>
, <name><surname>Gheit</surname>
<given-names>T</given-names>
</name>
, <name><surname>Iannacone</surname>
<given-names>MR</given-names>
</name>
, <name><surname>Vignat</surname>
<given-names>J</given-names>
</name>
, <name><surname>Sylla</surname>
<given-names>BS</given-names>
</name>
<etal>et al.</etal>
(<year>2013</year>
) <article-title>HPV16 genetic variation and the development of cervical cancer worldwide</article-title>
. <source>Br J Cancer</source>
<volume>108</volume>
: <fpage>240</fpage>
-<lpage>244</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1038/bjc.2012.508">10.1038/bjc.2012.508</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/23169278">23169278</ext-link>
<pub-id pub-id-type="pmid">23169278</pub-id>
</mixed-citation>
</ref>
<ref id="B51"><label>51</label>
<mixed-citation publication-type="journal"><name><surname>Burk</surname>
<given-names>RD</given-names>
</name>
, <name><surname>Harari</surname>
<given-names>A</given-names>
</name>
, <name><surname>Chen</surname>
<given-names>Z</given-names>
</name>
(<year>2013</year>
) <article-title>Human papillomavirus genome variants</article-title>
. <source>Virology</source>
<volume>445</volume>
: <fpage>232</fpage>
-<lpage>243</lpage>
.<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1016/j.virol.2013.07.018">10.1016/j.virol.2013.07.018</ext-link>
PubMed: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/pubmed/23998342">23998342</ext-link>
<pub-id pub-id-type="pmid">23998342</pub-id>
</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 000131 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000131 | 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. | ![]() |