Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Identification of a Novel Human Papillomavirus by Metagenomic Analysis of Samples from Patients with Febrile Respiratory Illness</title><author><name sortKey="Mokili, John L" sort="Mokili, John L" uniqKey="Mokili J" first="John L." last="Mokili">John L. Mokili</name><affiliation><nlm:aff id="aff1"><addr-line>Department of Biology, San Diego State University, San Diego, California, United States of America</addr-line></nlm:aff></affiliation><affiliation><nlm:aff id="aff2"><addr-line>Global Viral Forecasting (now known as Metabiota), San Francisco, California, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Dutilh, Bas E" sort="Dutilh, Bas E" uniqKey="Dutilh B" first="Bas E." last="Dutilh">Bas E. Dutilh</name><affiliation><nlm:aff id="aff3"><addr-line>Centre for Molecular and Biomolecular Informatics/Nijmegen Centre for Molecular Life Sciences, Radboud University Medical Centre, Nijmegen, The Netherlands</addr-line></nlm:aff></affiliation></author><author><name sortKey="Lim, Yan Wei" sort="Lim, Yan Wei" uniqKey="Lim Y" first="Yan Wei" last="Lim">Yan Wei Lim</name><affiliation><nlm:aff id="aff1"><addr-line>Department of Biology, San Diego State University, San Diego, California, United States of America</addr-line></nlm:aff></affiliation><affiliation><nlm:aff id="aff2"><addr-line>Global Viral Forecasting (now known as Metabiota), San Francisco, California, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Schneider, Bradley S" sort="Schneider, Bradley S" uniqKey="Schneider B" first="Bradley S." last="Schneider">Bradley S. Schneider</name><affiliation><nlm:aff id="aff2"><addr-line>Global Viral Forecasting (now known as Metabiota), San Francisco, California, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Taylor, Travis" sort="Taylor, Travis" uniqKey="Taylor T" first="Travis" last="Taylor">Travis Taylor</name><affiliation><nlm:aff id="aff2"><addr-line>Global Viral Forecasting (now known as Metabiota), San Francisco, California, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Haynes, Matthew R" sort="Haynes, Matthew R" uniqKey="Haynes M" first="Matthew R." last="Haynes">Matthew R. Haynes</name><affiliation><nlm:aff id="aff1"><addr-line>Department of Biology, San Diego State University, San Diego, California, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Metzgar, David" sort="Metzgar, David" uniqKey="Metzgar D" first="David" last="Metzgar">David Metzgar</name><affiliation><nlm:aff id="aff4"><addr-line>Naval Health Research Center, U.S. Navy, San Diego, California, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Myers, Christopher A" sort="Myers, Christopher A" uniqKey="Myers C" first="Christopher A." last="Myers">Christopher A. Myers</name><affiliation><nlm:aff id="aff4"><addr-line>Naval Health Research Center, U.S. Navy, San Diego, California, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Blair, Patrick J" sort="Blair, Patrick J" uniqKey="Blair P" first="Patrick J." last="Blair">Patrick J. Blair</name><affiliation><nlm:aff id="aff4"><addr-line>Naval Health Research Center, U.S. Navy, San Diego, California, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Nosrat, Bahador" sort="Nosrat, Bahador" uniqKey="Nosrat B" first="Bahador" last="Nosrat">Bahador Nosrat</name><affiliation><nlm:aff id="aff1"><addr-line>Department of Biology, San Diego State University, San Diego, California, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Wolfe, Nathan D" sort="Wolfe, Nathan D" uniqKey="Wolfe N" first="Nathan D." last="Wolfe">Nathan D. Wolfe</name><affiliation><nlm:aff id="aff2"><addr-line>Global Viral Forecasting (now known as Metabiota), San Francisco, California, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Rohwer, Forest" sort="Rohwer, Forest" uniqKey="Rohwer F" first="Forest" last="Rohwer">Forest Rohwer</name><affiliation><nlm:aff id="aff1"><addr-line>Department of Biology, San Diego State University, San Diego, California, United States of America</addr-line></nlm:aff></affiliation><affiliation><nlm:aff id="aff5"><addr-line>Center for Microbial Sciences, San Diego State University, San Diego, California, United States of America</addr-line></nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">23554892</idno><idno type="pmc">3600855</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3600855</idno><idno type="RBID">PMC:3600855</idno><idno type="doi">10.1371/journal.pone.0058404</idno><date when="2013">2013</date><idno type="wicri:Area/Pmc/Corpus">000129</idno><idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000129</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Identification of a Novel Human Papillomavirus by Metagenomic Analysis of Samples from Patients with Febrile Respiratory Illness</title><author><name sortKey="Mokili, John L" sort="Mokili, John L" uniqKey="Mokili J" first="John L." last="Mokili">John L. Mokili</name><affiliation><nlm:aff id="aff1"><addr-line>Department of Biology, San Diego State University, San Diego, California, United States of America</addr-line></nlm:aff></affiliation><affiliation><nlm:aff id="aff2"><addr-line>Global Viral Forecasting (now known as Metabiota), San Francisco, California, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Dutilh, Bas E" sort="Dutilh, Bas E" uniqKey="Dutilh B" first="Bas E." last="Dutilh">Bas E. Dutilh</name><affiliation><nlm:aff id="aff3"><addr-line>Centre for Molecular and Biomolecular Informatics/Nijmegen Centre for Molecular Life Sciences, Radboud University Medical Centre, Nijmegen, The Netherlands</addr-line></nlm:aff></affiliation></author><author><name sortKey="Lim, Yan Wei" sort="Lim, Yan Wei" uniqKey="Lim Y" first="Yan Wei" last="Lim">Yan Wei Lim</name><affiliation><nlm:aff id="aff1"><addr-line>Department of Biology, San Diego State University, San Diego, California, United States of America</addr-line></nlm:aff></affiliation><affiliation><nlm:aff id="aff2"><addr-line>Global Viral Forecasting (now known as Metabiota), San Francisco, California, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Schneider, Bradley S" sort="Schneider, Bradley S" uniqKey="Schneider B" first="Bradley S." last="Schneider">Bradley S. Schneider</name><affiliation><nlm:aff id="aff2"><addr-line>Global Viral Forecasting (now known as Metabiota), San Francisco, California, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Taylor, Travis" sort="Taylor, Travis" uniqKey="Taylor T" first="Travis" last="Taylor">Travis Taylor</name><affiliation><nlm:aff id="aff2"><addr-line>Global Viral Forecasting (now known as Metabiota), San Francisco, California, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Haynes, Matthew R" sort="Haynes, Matthew R" uniqKey="Haynes M" first="Matthew R." last="Haynes">Matthew R. Haynes</name><affiliation><nlm:aff id="aff1"><addr-line>Department of Biology, San Diego State University, San Diego, California, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Metzgar, David" sort="Metzgar, David" uniqKey="Metzgar D" first="David" last="Metzgar">David Metzgar</name><affiliation><nlm:aff id="aff4"><addr-line>Naval Health Research Center, U.S. Navy, San Diego, California, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Myers, Christopher A" sort="Myers, Christopher A" uniqKey="Myers C" first="Christopher A." last="Myers">Christopher A. Myers</name><affiliation><nlm:aff id="aff4"><addr-line>Naval Health Research Center, U.S. Navy, San Diego, California, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Blair, Patrick J" sort="Blair, Patrick J" uniqKey="Blair P" first="Patrick J." last="Blair">Patrick J. Blair</name><affiliation><nlm:aff id="aff4"><addr-line>Naval Health Research Center, U.S. Navy, San Diego, California, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Nosrat, Bahador" sort="Nosrat, Bahador" uniqKey="Nosrat B" first="Bahador" last="Nosrat">Bahador Nosrat</name><affiliation><nlm:aff id="aff1"><addr-line>Department of Biology, San Diego State University, San Diego, California, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Wolfe, Nathan D" sort="Wolfe, Nathan D" uniqKey="Wolfe N" first="Nathan D." last="Wolfe">Nathan D. Wolfe</name><affiliation><nlm:aff id="aff2"><addr-line>Global Viral Forecasting (now known as Metabiota), San Francisco, California, United States of America</addr-line></nlm:aff></affiliation></author><author><name sortKey="Rohwer, Forest" sort="Rohwer, Forest" uniqKey="Rohwer F" first="Forest" last="Rohwer">Forest Rohwer</name><affiliation><nlm:aff id="aff1"><addr-line>Department of Biology, San Diego State University, San Diego, California, United States of America</addr-line></nlm:aff></affiliation><affiliation><nlm:aff id="aff5"><addr-line>Center for Microbial Sciences, San Diego State University, San Diego, California, United States of America</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>As part of a virus discovery investigation using a metagenomic approach, a highly divergent novel <italic>Human papillomavirus</italic> type was identified in pooled convenience nasal/oropharyngeal swab samples collected from patients with febrile respiratory illness. Phylogenetic analysis of the whole genome and the L1 gene reveals that the new HPV identified in this study clusters with previously described gamma papillomaviruses, sharing only 61.1% (whole genome) and 63.1% (L1) sequence identity with its closest relative in the Papillomavirus episteme (PAVE) database. This new virus was named HPV_SD2 pending official classification. The complete genome of HPV-SD2 is 7,299 bp long (36.3% G/C) and contains 7 open reading frames (L2, L1, E6, E7, E1, E2 and E4) and a non-coding long control region (LCR) between L1 and E6. The metagenomic procedures, coupled with the bioinformatic methods described herein are well suited to detect small circular genomes such as those of human papillomaviruses.</p></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Bernard, Hu" uniqKey="Bernard H">HU Bernard</name></author><author><name sortKey="Burk, Rd" uniqKey="Burk R">RD Burk</name></author><author><name sortKey="Chen, Z" uniqKey="Chen Z">Z Chen</name></author><author><name sortKey="Van Doorslaer, K" uniqKey="Van Doorslaer K">K van Doorslaer</name></author><author><name sortKey="Hausen, H" uniqKey="Hausen H">H Hausen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="De Villiers, Em" uniqKey="De Villiers E">EM de Villiers</name></author><author><name sortKey="Fauquet, C" uniqKey="Fauquet C">C Fauquet</name></author><author><name sortKey="Broker, Tr" uniqKey="Broker T">TR Broker</name></author><author><name sortKey="Bernard, Hu" uniqKey="Bernard H">HU Bernard</name></author><author><name sortKey="Zur Hausen, H" uniqKey="Zur Hausen H">H zur Hausen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ramamoorthy, S" uniqKey="Ramamoorthy S">S Ramamoorthy</name></author><author><name sortKey="Liu, Yt" uniqKey="Liu Y">YT Liu</name></author><author><name sortKey="Luo, L" uniqKey="Luo L">L Luo</name></author><author><name sortKey="Miyai, K" uniqKey="Miyai K">K Miyai</name></author><author><name sortKey="Lu, Q" uniqKey="Lu Q">Q Lu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bartholomew, Da" uniqKey="Bartholomew D">DA Bartholomew</name></author><author><name sortKey="Luff, Rd" uniqKey="Luff R">RD Luff</name></author><author><name sortKey="Quigley, Nb" uniqKey="Quigley N">NB Quigley</name></author><author><name sortKey="Curtis, M" uniqKey="Curtis M">M Curtis</name></author><author><name sortKey="Olson, Mc" uniqKey="Olson M">MC Olson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chaturvedi, Ak" uniqKey="Chaturvedi A">AK Chaturvedi</name></author><author><name sortKey="Engels, Ea" uniqKey="Engels E">EA Engels</name></author><author><name sortKey="Pfeiffer, Rm" uniqKey="Pfeiffer R">RM Pfeiffer</name></author><author><name sortKey="Hernandez, By" uniqKey="Hernandez B">BY Hernandez</name></author><author><name sortKey="Xiao, W" uniqKey="Xiao W">W Xiao</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jemal, A" uniqKey="Jemal A">A Jemal</name></author><author><name sortKey="Bray, F" uniqKey="Bray F">F Bray</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tse, H" uniqKey="Tse H">H Tse</name></author><author><name sortKey="Tsang, Ak" uniqKey="Tsang A">AK Tsang</name></author><author><name sortKey="Tsoi, Hw" uniqKey="Tsoi H">HW Tsoi</name></author><author><name sortKey="Leung, As" uniqKey="Leung A">AS Leung</name></author><author><name sortKey="Ho, Cc" uniqKey="Ho C">CC Ho</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Forslund, O" uniqKey="Forslund O">O Forslund</name></author><author><name sortKey="Ly, H" uniqKey="Ly H">H Ly</name></author><author><name sortKey="Higgins, G" uniqKey="Higgins G">G Higgins</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Riesenfeld, Cs" uniqKey="Riesenfeld C">CS Riesenfeld</name></author><author><name sortKey="Schloss, Pd" uniqKey="Schloss P">PD Schloss</name></author><author><name sortKey="Handelsman, J" uniqKey="Handelsman J">J Handelsman</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Edwards, Ra" uniqKey="Edwards R">RA Edwards</name></author><author><name sortKey="Rohwer, F" uniqKey="Rohwer F">F Rohwer</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Delwart, El" uniqKey="Delwart E">EL Delwart</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Mokili, Jl" uniqKey="Mokili J">JL Mokili</name></author><author><name sortKey="Rohwer, F" uniqKey="Rohwer F">F Rohwer</name></author><author><name sortKey="Dutilh, Be" uniqKey="Dutilh B">BE Dutilh</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Dong, X" uniqKey="Dong X">X Dong</name></author><author><name sortKey="Stothard, P" uniqKey="Stothard P">P Stothard</name></author><author><name sortKey="Forsythe, Ij" uniqKey="Forsythe I">IJ Forsythe</name></author><author><name sortKey="Wishart, Ds" uniqKey="Wishart D">DS Wishart</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Deutsch, Mj" uniqKey="Deutsch M">MJ Deutsch</name></author><author><name sortKey="Ott, E" uniqKey="Ott E">E Ott</name></author><author><name sortKey="Papior, P" uniqKey="Papior P">P Papior</name></author><author><name sortKey="Schepers, A" uniqKey="Schepers A">A Schepers</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ott, E" uniqKey="Ott E">E Ott</name></author><author><name sortKey="Norio, P" uniqKey="Norio P">P Norio</name></author><author><name sortKey="Ritzi, M" uniqKey="Ritzi M">M Ritzi</name></author><author><name sortKey="Schildkraut, C" uniqKey="Schildkraut C">C Schildkraut</name></author><author><name sortKey="Schepers, A" uniqKey="Schepers A">A Schepers</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sarkari, F" uniqKey="Sarkari F">F Sarkari</name></author><author><name sortKey="Sanchez Alcaraz, T" uniqKey="Sanchez Alcaraz T">T Sanchez-Alcaraz</name></author><author><name sortKey="Wang, S" uniqKey="Wang S">S Wang</name></author><author><name sortKey="Holowaty, Mn" uniqKey="Holowaty M">MN Holowaty</name></author><author><name sortKey="Sheng, Y" uniqKey="Sheng Y">Y Sheng</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Krubke, J" uniqKey="Krubke J">J Krubke</name></author><author><name sortKey="Kraus, J" uniqKey="Kraus J">J Kraus</name></author><author><name sortKey="Delius, H" uniqKey="Delius H">H Delius</name></author><author><name sortKey="Chow, L" uniqKey="Chow L">L Chow</name></author><author><name sortKey="Broker, T" uniqKey="Broker T">T Broker</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Androphy, Ej" uniqKey="Androphy E">EJ Androphy</name></author><author><name sortKey="Lowy, Dr" uniqKey="Lowy D">DR Lowy</name></author><author><name sortKey="Schiller, Jt" uniqKey="Schiller J">JT Schiller</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Spalholz, Ba" uniqKey="Spalholz B">BA Spalholz</name></author><author><name sortKey="Yang, Yc" uniqKey="Yang Y">YC Yang</name></author><author><name sortKey="Howley, Pm" uniqKey="Howley P">PM Howley</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Figge, J" uniqKey="Figge J">J Figge</name></author><author><name sortKey="Webster, T" uniqKey="Webster T">T Webster</name></author><author><name sortKey="Smith, Tf" uniqKey="Smith T">TF Smith</name></author><author><name sortKey="Paucha, E" uniqKey="Paucha E">E Paucha</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Castresana, J" uniqKey="Castresana J">J Castresana</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Dean, Fb" uniqKey="Dean F">FB Dean</name></author><author><name sortKey="Nelson, Jr" uniqKey="Nelson J">JR Nelson</name></author><author><name sortKey="Giesler, Tl" uniqKey="Giesler T">TL Giesler</name></author><author><name sortKey="Lasken, Rs" uniqKey="Lasken R">RS Lasken</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wegley, L" uniqKey="Wegley L">L Wegley</name></author><author><name sortKey="Mosier Boss, P" uniqKey="Mosier Boss P">P Mosier-Boss</name></author><author><name sortKey="Lieberman, S" uniqKey="Lieberman S">S Lieberman</name></author><author><name sortKey="Andrews, J" uniqKey="Andrews J">J Andrews</name></author><author><name sortKey="Graff Baker, A" uniqKey="Graff Baker A">A Graff-Baker</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Margulies, M" uniqKey="Margulies M">M Margulies</name></author><author><name sortKey="Egholm, M" uniqKey="Egholm M">M Egholm</name></author><author><name sortKey="Altman, We" uniqKey="Altman W">WE Altman</name></author><author><name sortKey="Attiya, S" uniqKey="Attiya S">S Attiya</name></author><author><name sortKey="Bader, Js" uniqKey="Bader J">JS Bader</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Schmieder, R" uniqKey="Schmieder R">R Schmieder</name></author><author><name sortKey="Edwards, R" uniqKey="Edwards R">R Edwards</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Schmieder, R" uniqKey="Schmieder R">R Schmieder</name></author><author><name sortKey="Lim, Yw" uniqKey="Lim Y">YW Lim</name></author><author><name sortKey="Rohwer, F" uniqKey="Rohwer F">F Rohwer</name></author><author><name sortKey="Edwards, R" uniqKey="Edwards R">R Edwards</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tomlins, Sa" uniqKey="Tomlins S">SA Tomlins</name></author><author><name sortKey="Mehra, R" uniqKey="Mehra R">R Mehra</name></author><author><name sortKey="Rhodes, Dr" uniqKey="Rhodes D">DR Rhodes</name></author><author><name sortKey="Shah, Rb" uniqKey="Shah R">RB Shah</name></author><author><name sortKey="Rubin, Ma" uniqKey="Rubin M">MA Rubin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Altschul, Sf" uniqKey="Altschul S">SF Altschul</name></author><author><name sortKey="Gish, W" uniqKey="Gish W">W Gish</name></author><author><name sortKey="Miller, W" uniqKey="Miller W">W Miller</name></author><author><name sortKey="Myers, Ew" uniqKey="Myers E">EW Myers</name></author><author><name sortKey="Lipman, Dj" uniqKey="Lipman D">DJ Lipman</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Trindade Silva, Ae" uniqKey="Trindade Silva A">AE Trindade-Silva</name></author><author><name sortKey="Rua, C" uniqKey="Rua C">C Rua</name></author><author><name sortKey="Silva, Gg" uniqKey="Silva G">GG Silva</name></author><author><name sortKey="Dutilh, Be" uniqKey="Dutilh B">BE Dutilh</name></author><author><name sortKey="Moreira, Ap" uniqKey="Moreira A">AP Moreira</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Guindon, S" uniqKey="Guindon S">S Guindon</name></author><author><name sortKey="Gascuel, O" uniqKey="Gascuel O">O Gascuel</name></author></analytic></biblStruct><biblStruct></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">23554892</article-id><article-id pub-id-type="pmc">3600855</article-id><article-id pub-id-type="publisher-id">PONE-D-12-33549</article-id><article-id pub-id-type="doi">10.1371/journal.pone.0058404</article-id><article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Biology</subject><subj-group><subject>Genomics</subject><subj-group><subject>Metagenomics</subject></subj-group></subj-group><subj-group><subject>Microbiology</subject><subj-group><subject>Virology</subject></subj-group></subj-group></subj-group></article-categories><title-group><article-title>Identification of a Novel Human Papillomavirus by Metagenomic Analysis of Samples from Patients with Febrile Respiratory Illness</article-title><alt-title alt-title-type="running-head">Discovery of a Novel HPV by Metagenomic Approach</alt-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Mokili</surname><given-names>John L.</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="aff" rid="aff2"><sup>2</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author"><name><surname>Dutilh</surname><given-names>Bas E.</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref></contrib><contrib contrib-type="author"><name><surname>Lim</surname><given-names>Yan Wei</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author"><name><surname>Schneider</surname><given-names>Bradley S.</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author"><name><surname>Taylor</surname><given-names>Travis</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author"><name><surname>Haynes</surname><given-names>Matthew R.</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="author-notes" rid="fn2"><sup>¤b</sup></xref></contrib><contrib contrib-type="author"><name><surname>Metzgar</surname><given-names>David</given-names></name><xref ref-type="aff" rid="aff4"><sup>4</sup></xref><xref ref-type="author-notes" rid="fn1"><sup>¤a</sup></xref></contrib><contrib contrib-type="author"><name><surname>Myers</surname><given-names>Christopher A.</given-names></name><xref ref-type="aff" rid="aff4"><sup>4</sup></xref></contrib><contrib contrib-type="author"><name><surname>Blair</surname><given-names>Patrick J.</given-names></name><xref ref-type="aff" rid="aff4"><sup>4</sup></xref></contrib><contrib contrib-type="author"><name><surname>Nosrat</surname><given-names>Bahador</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="author-notes" rid="fn2"><sup>¤b</sup></xref></contrib><contrib contrib-type="author"><name><surname>Wolfe</surname><given-names>Nathan D.</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author"><name><surname>Rohwer</surname><given-names>Forest</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="aff" rid="aff5"><sup>5</sup></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>Department of Biology, San Diego State University, San Diego, California, United States of America</addr-line></aff><aff id="aff2"><label>2</label><addr-line>Global Viral Forecasting (now known as Metabiota), San Francisco, California, United States of America</addr-line></aff><aff id="aff3"><label>3</label><addr-line>Centre for Molecular and Biomolecular Informatics/Nijmegen Centre for Molecular Life Sciences, Radboud University Medical Centre, Nijmegen, The Netherlands</addr-line></aff><aff id="aff4"><label>4</label><addr-line>Naval Health Research Center, U.S. Navy, San Diego, California, United States of America</addr-line></aff><aff id="aff5"><label>5</label><addr-line>Center for Microbial Sciences, San Diego State University, San Diego, California, United States of America</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>jmokili@gmail.com</email></corresp><fn fn-type="conflict"><p><bold>Competing Interests: </bold>JLM, YWL, BSS and TT were associated with the Global Viral Forecasting Initiative, now known as Metabiota, during the time of the execution of the study. These authors have declared that no competing interests exist. DM was associated with the Naval Health Research Center, U.S. Navy, San Diego, CA, USA, during the study period. DM has moved to Ibis Biosciences, Carlsbad, CA. DM has declared that no competing interests exist. These authors fully adhere to all the PLOS ONE policies on sharing data and materials.</p></fn><fn fn-type="con"><p>Collected samples, obtained IRB clearance: PJB DM CAM. Read, edited and approved the manuscript: JLM BED YWL BSS TT MRH DM CAM PJB BN NDW FLR. Conceived and designed the experiments: JLM FLR NDW PJB DM. Performed the experiments: JLM YWL MRH. Analyzed the data: JLM BED TT BN. Contributed reagents/materials/analysis tools: NDW BSS FLR. Wrote the paper: JLM.</p></fn><fn id="fn1" fn-type="current-aff"><label>¤a</label><p>Current address: Ibis Biosciences, Carlsbad, California, United States of America</p></fn><fn id="fn2" fn-type="current-aff"><label>¤b</label><p>Current address: Department of Energy, Joint Genome Institute, Walnut Creek, California, United States of America</p></fn></author-notes><pub-date pub-type="collection"><year>2013</year></pub-date><pub-date pub-type="epub"><day>15</day><month>3</month><year>2013</year></pub-date><volume>8</volume><issue>3</issue><elocation-id>e58404</elocation-id><history><date date-type="received"><day>27</day><month>10</month><year>2012</year></date><date date-type="accepted"><day>4</day><month>2</month><year>2013</year></date></history><permissions><copyright-year>2013</copyright-year><copyright-holder>Mokili 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>As part of a virus discovery investigation using a metagenomic approach, a highly divergent novel <italic>Human papillomavirus</italic> type was identified in pooled convenience nasal/oropharyngeal swab samples collected from patients with febrile respiratory illness. Phylogenetic analysis of the whole genome and the L1 gene reveals that the new HPV identified in this study clusters with previously described gamma papillomaviruses, sharing only 61.1% (whole genome) and 63.1% (L1) sequence identity with its closest relative in the Papillomavirus episteme (PAVE) database. This new virus was named HPV_SD2 pending official classification. The complete genome of HPV-SD2 is 7,299 bp long (36.3% G/C) and contains 7 open reading frames (L2, L1, E6, E7, E1, E2 and E4) and a non-coding long control region (LCR) between L1 and E6. The metagenomic procedures, coupled with the bioinformatic methods described herein are well suited to detect small circular genomes such as those of human papillomaviruses.</p></abstract><funding-group><funding-statement>BED is supported by the Dutch Science Foundation (NWO) Veni grant (016.111.075). Global Viral Forecasting (JLM, BSS and NDW) was graciously supported by the U.S. Department of Defense Global Emerging Infections, Surveillance and Response Systems (DoD GEIS) and the Defense Threat Reduction Agency (DTRA), Cooperative Biological Engagement Program (CBEP), Google.org, the Skoll Foundation, and the U.S. Agency for International Development (USAID) Emerging Pandemic Threats Program, PREDICT project, under the terms of Cooperative Agreement Number GHN-A-OO-09-00010-00. The Surveillance conducted by Naval Health Research Center (NHRC) was supported by the Global Emerging Infections System division of the Armed Forces Health Surveillance Center. The views expressed in this work are those of the authors and do not reflect the official policy of the Department of the Navy, Department of Defense, or the US Government. Approved for public releases; distribution is unlimited. This research has been conducted in compliance with all applicable federal regulations governing the protection of human subjects in research. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</funding-statement></funding-group><counts><page-count count="12"></page-count></counts></article-meta></front><body><sec id="s1"><title>Introduction</title><p>Papillomaviruses (PV) are small, circular, double-stranded DNA viruses belonging to <italic>Papillomaviridae</italic>, a large family with over 180 viruses phylogenetically classified in groups, genera, and species <xref ref-type="bibr" rid="pone.0058404-Bernard1">[1]</xref>, <xref ref-type="bibr" rid="pone.0058404-deVilliers1">[2]</xref>. The PVs infecting human, the human PV (HPVs), encompass highly pathogenic viruses causing a wide range of diseases including cervical, vulval, penile, oropharyngeal and anal warts and cancers <xref ref-type="bibr" rid="pone.0058404-Ramamoorthy1">[3]</xref>–<xref ref-type="bibr" rid="pone.0058404-Jemal1">[6]</xref>. Over 150 genomes of HPVs have been fully sequenced, characterized and catalogued by the Papillomavirus Episteme (PAVE) database. It is likely that more unknown HPV types are yet to be discovered, and may be circulating between humans, animals and their respective environments. It is necessary to determine the full picture of the genetic diversity, the host-range and the clinical relevance of the vast majority of HPVs, which thus far remain unclear. One obstacle to discovering novel HPVs viruses has been the lack of suitable laboratory techniques and computational methods to survey environments and biological samples.</p><p>While various laboratory techniques for the detection of known HPVs are readily available, they are unpredictable for find highly divergent viruses in biological and ecological environments. Although generally viewed as the gold standard of novel virus discovery, tissue culture techniques are inefficient to propagate HPVs <xref ref-type="bibr" rid="pone.0058404-deVilliers1">[2]</xref>. Therefore, the detection of HPVs primarily relies on molecular techniques such as PCR using consensus primers <xref ref-type="bibr" rid="pone.0058404-Tse1">[7]</xref>, <xref ref-type="bibr" rid="pone.0058404-Forslund1">[8]</xref>, and the sequencing of a subset of the genome. The subsequent genotyping of the partial genome of HPV is generally carried out by pair-wise comparison of the sequences with the corresponding genes of known viruses. This approach can lead to analytical difficulties when non-overlapping genomic regions of different viruses are compared. Moreover, PCR using consensus primers is inadequate for highly distant and novel HPVs that can only detected by metagenomics <xref ref-type="bibr" rid="pone.0058404-Tse1">[7]</xref>. In order to assign taxonomical classification, the PV Working Group of the International Committee on Taxonomy of Viruses recommends that researchers analyze at least the complete gene encoding the major capsid protein, L1 <xref ref-type="bibr" rid="pone.0058404-Bernard1">[1]</xref>.</p><p>With the advent of metagenomics using high-throughput sequencing technology <xref ref-type="bibr" rid="pone.0058404-Riesenfeld1">[9]</xref>–<xref ref-type="bibr" rid="pone.0058404-Mokili1">[12]</xref>, it is possible to generate the full genome of phylogenetically HPVs at relatively low cost and unprecedented speed, regardless of genetic divergence from known viruses. Metagenomics is a culture-independent technique and does not depend on <italic>a priori</italic> knowledge of the genetic information of the organism to be detected <xref ref-type="bibr" rid="pone.0058404-Riesenfeld1">[9]</xref>–<xref ref-type="bibr" rid="pone.0058404-Mokili1">[12]</xref>.</p><p>As part of a virus discovery consortium between San Diego State University, Global Viral Forecasting and the Naval Health Research Center, we applied the metagenomic approach to pooled convenience samples obtained from patients with respiratory illness, and identified a novel gamma papillomavirus that is highly divergent compared to previously known HPVs.</p></sec><sec id="s2"><title>Results and Discussion</title><sec id="s2a"><title>Overall Taxonomic Classification of Metagenomic Reads</title><p>Sequencing reads generated by high-throughput sequencing were assembled with a threshold of 95% identity and at least 45 nt overlap and subjected to taxonomical classification using BLASTn against the Genbank database employing a stringent cut-off E-value of ≤10<sup>−5</sup>.</p><p>The proportion of sequence reads with significant hits to viruses, bacteria and eukaryotes are summarized in <xref ref-type="fig" rid="pone-0058404-g001">Figure 1</xref> (details provided as Document S1). Overall, of the 2,027,673 reads generated in this study, 120,801 (6%) were classified as bacterial, 30,259 (1.5%) mapped to known viruses and 6,373 (0.3%) were classified as originating from eukaryotic cells. Details of viruses identified in this study by BLASTn analysis are presented in <xref ref-type="table" rid="pone-0058404-t001">Table 1</xref>, showing that both DNA and RNA viruses were identified in nasopharyngeal samples. These included human enterovirus, Torque teno virus, human rhinovirus, SEN virus and HPV. The HPV identified by BLAST consists of 5 contigs with significant hit to HPV_49. We combined these 5 contigs into a scaffold of length 6,356 bp, and will further refer to this scaffold as HPV_SD1 in the further analyses (below). The remaining 1,870,240 reads (92%) could not be classified as belonging to previously defined taxa. Such sequences are referred to as ‘unknowns’ throughout this paper.</p><fig id="pone-0058404-g001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0058404.g001</object-id><label>Figure 1</label><caption><title>Taxonomical assignment of sequence reads convenient respiratory waste samples.</title><p>A. Overall proportions of sequences with homolog in Genbank (the known, 7.8%) compared to sequences with no homolog in Genbank (Unknown, 92.2%). Unknown/Divergent indicates the proportion of highly divergent and/or novel sequences with no homology to NCBI. B. Proportion of sequences classified as eukaryotes, bacteria and viruses. Sequences were classified using BLASTn search against all non-redundant nucleotide sequences in the NCBI nt database with an E-value cutoff of 10<sup>−5</sup>.</p></caption><graphic xlink:href="pone.0058404.g001"></graphic></fig><table-wrap id="pone-0058404-t001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0058404.t001</object-id><label>Table 1</label><caption><title>Summary of known viruses detected in nasopharyngeal samples.</title></caption><alternatives><graphic id="pone-0058404-t001-1" xlink:href="pone.0058404.t001"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Viral species</td><td align="left" rowspan="1" colspan="1">DNA metagenome</td><td align="left" rowspan="1" colspan="1">RNA metagenome</td><td align="left" rowspan="1" colspan="1">Total Read counts</td><td align="left" rowspan="1" colspan="1">Percentage total read</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">17,313</td><td align="left" rowspan="1" colspan="1">12,946</td><td align="left" rowspan="1" colspan="1">30,259</td><td align="left" rowspan="1" colspan="1">100%</td></tr><tr><td align="left" rowspan="1" colspan="1">Human enterovirus C</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">9,926</td><td align="left" rowspan="1" colspan="1">9,926</td><td align="left" rowspan="1" colspan="1">36.32</td></tr><tr><td align="left" rowspan="1" colspan="1">Torque teno virus</td><td align="left" rowspan="1" colspan="1">8,090</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">8,090</td><td align="left" rowspan="1" colspan="1">29.61</td></tr><tr><td align="left" rowspan="1" colspan="1">Torque teno mini virus</td><td align="left" rowspan="1" colspan="1">3,722</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">3,722</td><td align="left" rowspan="1" colspan="1">13.62</td></tr><tr><td align="left" rowspan="1" colspan="1">Human rhinovirus A</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">1,531</td><td align="left" rowspan="1" colspan="1">1,531</td><td align="left" rowspan="1" colspan="1">5.60</td></tr><tr><td align="left" rowspan="1" colspan="1">SEN virus</td><td align="left" rowspan="1" colspan="1">1,424</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">1,424</td><td align="left" rowspan="1" colspan="1">5.21</td></tr><tr><td align="left" rowspan="1" colspan="1">Human rhinovirus C</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">589</td><td align="left" rowspan="1" colspan="1">589</td><td align="left" rowspan="1" colspan="1">2.16</td></tr><tr><td align="left" rowspan="1" colspan="1">Human papillomavirus</td><td align="left" rowspan="1" colspan="1">571</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">571</td><td align="left" rowspan="1" colspan="1">2.09</td></tr><tr><td align="left" rowspan="1" colspan="1">Lactobacillus phage mv4</td><td align="left" rowspan="1" colspan="1">429</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">429</td><td align="left" rowspan="1" colspan="1">1.57</td></tr><tr><td align="left" rowspan="1" colspan="1">Small anellovirus</td><td align="left" rowspan="1" colspan="1">254</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">254</td><td align="left" rowspan="1" colspan="1">0.93</td></tr><tr><td align="left" rowspan="1" colspan="1">Human rhinovirus B</td><td align="left" rowspan="1" colspan="1">245</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">245</td><td align="left" rowspan="1" colspan="1">0.90</td></tr><tr><td align="left" rowspan="1" colspan="1">Human rhinovirus sp.</td><td align="left" rowspan="1" colspan="1">202</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">202</td><td align="left" rowspan="1" colspan="1">0.74</td></tr><tr><td align="left" rowspan="1" colspan="1">Cyclovirus Chimp11</td><td align="left" rowspan="1" colspan="1">168</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">168</td><td align="left" rowspan="1" colspan="1">0.61</td></tr><tr><td align="left" rowspan="1" colspan="1">Human papillomavirus 49</td><td align="left" rowspan="1" colspan="1">76</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">76</td><td align="left" rowspan="1" colspan="1">0.28</td></tr><tr><td align="left" rowspan="1" colspan="1">TTV-like mini virus</td><td align="left" rowspan="1" colspan="1">67</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">67</td><td align="left" rowspan="1" colspan="1">0.25</td></tr><tr><td align="left" rowspan="1" colspan="1">Torque teno mini virus 5</td><td align="left" rowspan="1" colspan="1">32</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">32</td><td align="left" rowspan="1" colspan="1">0.12</td></tr><tr><td align="left" rowspan="1" colspan="1">Streptococcus phage 5093</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">0.01</td></tr></tbody></table></alternatives><table-wrap-foot><fn id="nt101"><label></label><p>BLASTn results. The new HPV-SD2 was not detected by the BLASTn analysis using a cut off E value of 10<sup>−5</sup>.</p></fn></table-wrap-foot></table-wrap></sec><sec id="s2b"><title>Deciphering the Unknowns</title><p>Further analyses were performed to determine the origin of the unknown sequences. First, we developed a script (circular_from_ACE.pl, Document S2) to determine whether the contigs were circular. Second, any circular genomes identified were subjected to PlasMapper analysis <xref ref-type="bibr" rid="pone.0058404-Dong1">[13]</xref> to determine if they were bacteria or plasmids. Third, all unknown sequences were tested with tBLASTx using a very relaxed E value cutoff. Finally, sequences with hits to known viruses with relaxed E value were confirmed by phylogenetic analysis.</p><p>The circular_from_ACE.pl script used the ACE files generated by the <italic>De Novo</italic> Assembler and determined whether exact copies of overhanging sequences at both 5′ and 3′ ends of the assembly were found within the contigs (<xref ref-type="fig" rid="pone-0058404-g002">Figure 2</xref>.A). <xref ref-type="fig" rid="pone-0058404-g002">Figure 2B</xref> depicts the striking similarity between overhanging fragments at the 5′ and the 3′ ends, and sequences in their respective opposite ends of the contig. Of 2,075 contigs analyzed, 193 (9.3%) met the criteria suggestive of having a circular DNA structure.</p><fig id="pone-0058404-g002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0058404.g002</object-id><label>Figure 2</label><caption><title>Coverage plot of the circular structure of HPV-SD2.</title><p>A. Screenshot of the assembly of sequences using the <italic>De Novo Assembler</italic> set for 98% similarity and 45 bp overlap. The coverage plot is shown in green. Overhanging sequences at the 5′ and 3′ ends are partially assembled. B. Box at 5′ end and box at 3′ end are highlighted to show the similarity between unassembled overhang reads at 5′ and 3′ ends and corresponding ends of the consensus sequence.</p></caption><graphic xlink:href="pone.0058404.g002"></graphic></fig><p>Many viral genomes including HPVs and all bacterial plasmids are characteristically circular. However, plasmids have characteristic features, including the origin of replication, which can distinguish them from viruses. The feature is seldom found among viruses, with rare exceptions such as Epstein-Barr virus, which uses an origin of replication to remain in latency in human cells <xref ref-type="bibr" rid="pone.0058404-Deutsch1">[14]</xref>–<xref ref-type="bibr" rid="pone.0058404-Sarkari1">[16]</xref>. Thus we investigated whether the circular genomic elements were plasmids or viruses by screening the apparently circular contigs for plasmid-specific motifs, including the origin of replication, using PlasMapper <xref ref-type="bibr" rid="pone.0058404-Dong1">[13]</xref>. We found no evidence for these signals. The method for identification of circular genomes described here permitted us to characterize 9.3% of the unknown sequences including the putative novel papillomavirus identified in this study. Additional studies are needed to fully characterize these putative novel viral genomic structures.</p><p>All contigs were further analyzed with tBLASTx against NCBI sequences to assign taxonomy. This analysis revealed that a contig constructed with 599 sequencing reads had a significant similarity to HPV sequences. The contig was named HPV_SD2.</p></sec><sec id="s2c"><title>Genomic Organization of and Phylogenetic Analysis of HPV_SD2</title><p>The complete genome of HPV_SD2 (<xref ref-type="fig" rid="pone-0058404-g003">Figure 3</xref>) is 7,299 bp long with a GC content of 36.3%. The genomic organization of HPV_SD2 includes 7 open reading frames (L2, L1, E6, E7, E1, E2 and E4) characteristic of other previously described HPVs. In addition, there was a 530 bp non-coding long control region (LCR), also known as an upstream regulatory region (URR) located between the L1 and E6 genes, at position 3,244 through 3,773 of the genome (<xref ref-type="fig" rid="pone-0058404-g003">Figure 3B</xref>). This LCR contains a TATA box (TATAAA, nucleotide positions 3,735–3,740), a polyadenylated site (AATAAA, nucleotide positions 3,336–3,341), 3 palindrome sites (ACCG-N<sub>4</sub>-CGGT; nucleotide positions 3,483–3,494; 3,650–3,661 and 3,720–3,731) and 1 degenerate palindrome (ACC-N<sub>6</sub>-GGT, positions 3,524–3,535). The palindrome motif ACCG-N<sub>4</sub>-CGGT is known to play a crucial role in binding of the E2 protein <xref ref-type="bibr" rid="pone.0058404-Krubke1">[17]</xref>, <xref ref-type="bibr" rid="pone.0058404-Androphy1">[18]</xref>, which mediates viral transcription in <italic>trans</italic><xref ref-type="bibr" rid="pone.0058404-Spalholz1">[19]</xref>. The E6 protein contained 2 zinc-binding domains (CxxC-N<sub>29</sub>-CxxC) at positions 3,853–3,964 and 4,072–4,182.</p><fig id="pone-0058404-g003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0058404.g003</object-id><label>Figure 3</label><caption><title>Genomic organization of the HPV_SD2 virus.</title><p>Open reading frames (L2, L1, E6, E7, E1, E2, E4) and a 530 bp non-coding long control region (LCR) are shown. B. Details of the LCR region showing the TATA box (TATAAA, positions 3735–3740), a polyadenylated site (AATAAA, positions 3336–3341), 3 palindrome sites (ACCG-N<sub>4</sub>-CGGT; positions 3483–3494, 3650–3661, 3720–3731) and 1 degenerate palindrome (ACC-N<sub>6</sub>-GGT, positions 3524–3535). C. Metal-binding domains in deduced E6 and E7 proteins.</p></caption><graphic xlink:href="pone.0058404.g003"></graphic></fig><p>The E7 protein is considered as an oncogene of high-risk HPVs because it shares structural and functional similarities with known DNA viruses such as SV40 and adenovirus <xref ref-type="bibr" rid="pone.0058404-Krubke1">[17]</xref>, <xref ref-type="bibr" rid="pone.0058404-Figge1">[20]</xref>. A sequence coding for a metal-binding domain, including two Cys-X-X-Cys motifs (CGTC and CSGC) separated by 29 amino acids, was found within the deduced E7 gene at positions 4,344 through 4,447.</p><p>The HPV_SD2 contig from this study were aligned with 155 complete HPV genomes in the PAVE database. The alignment was carried out using ClustalW 2.1 (default parameters) and all ambiguities and gaps were removed with GBlocks 0.91 b <xref ref-type="bibr" rid="pone.0058404-Castresana1">[21]</xref>. We also included another genome from this study, HPV_SD1, identified by BLAST analysis. The phylogenetic analysis of both the complete genome and the L1 gene showed that while HPV_SD1 was closely related to HPV_49, HPV_SD2 was a distant relative of gamma papillomaviruses. Indeed, HPV_SD2 clustered strongly with the Gamma-human papillomavirus genus, with a bootstrap value of 100% (<xref ref-type="fig" rid="pone-0058404-g004">Figure 4</xref>). A pair-wise distance analysis performed with all 155 previously described complete genome sequences and the L1 gene of papillomavirus showed that the identity between HPV_SD2 and the closest taxon (HPV128) was only 61.1% and 63.1%, respectively (Document S3). The average pair-wise identity between HPV_SD2 and all previously described complete HPV genomes depicted in <xref ref-type="fig" rid="pone-0058404-g004">Figure 4</xref> was 49.4% (standard deviation 4.2%, Document S3).</p><fig id="pone-0058404-g004" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0058404.g004</object-id><label>Figure 4</label><caption><title>Maximum likelihood tree showing the clustering of HPV_SD2 with previously documented full-length genomes of human papillomaviruses.</title><p>The tree was visualized with the Interactive Tree of Life (iTOL) <xref ref-type="bibr" rid="pone.0058404-Letunic1">[32]</xref>. Bootstrap values with at least 50%, 75% or 100% (100 re-samplings) are indicated. Note, another HPV_SD1 from this study was added in this analysis clustered with reference HPV_49.</p></caption><graphic xlink:href="pone.0058404.g004"></graphic></fig></sec><sec id="s2d"><title>Verification of the Genome Structure of HPV_SD2 by PCR and Sanger Sequencing</title><p>We used PCR as an alternative detection method to determine if the HPV_SD2 contig generated by 454 sequencing could be found in the sample pool as a double stranded DNA. Three PCR amplifications (<xref ref-type="fig" rid="pone-0058404-g005">Figure 5A</xref>) using primer-pairs 671F-892R, 6838F-7076R and 6838F-6972R were designed to yield amplicons of 222, 135 and 239 base pairs, respectively. Because the <italic>Taq</italic> DNA polymerase was capable to amplify the templates beyond the primer sites, the primer sets were also expected to generate larger amplicons containing one or more full genome of HPV_SD2 (7,299 bp) plus a fragment of the size corresponding to the small amplicon product of the same primer pair (222, 135 and 239 bp). As depicted in <xref ref-type="fig" rid="pone-0058404-g005">Figure 5B</xref>. the PCR reactions with each primer pair generated the predicted amplicons of the projected sizes. The presence of these bands indicates that the HPV_SD2 is indeed a double stranded circular genome.</p><fig id="pone-0058404-g005" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0058404.g005</object-id><label>Figure 5</label><caption><title>Confirmatory PCRs to verify HPV_SD2 is a circular double-stranded genome.</title><p>A. Graphical representation of the binding sites of primers 671F-892R on the putative circular structure of HPV_SD2, and the predicted PCR product sizes (I: 222 bp, II: 7,591 bp, III: 14,890 bp). The predicted short band (I) indicates the amplification of the proximal region between primers 671F and 892R. The large band (II) indicates that <italic>Taq</italic> DNA polymerase would amplify the region between the primers 671F and 892R by making the full circle of the HPV genome. The large band (III) indicates the <italic>Taq</italic> DNA polymerase would make 2 full circle around the HPV genome. PCRs were also performed using primer sets 6838F-6972R, 6838F-7076R. B. Agarose gel (0.5%) showing the amplified HPV_SD2. Primer sets used are shown: primer sets 6838F-6972R, 6838F-7076R. For each PCR, the same sample pool was tested at different concentrations [1∶1 (lanes 1, 4, 7), 1∶10 (lanes 2, 5, 8) and 1∶50 (lanes (3, 6, 9)] using 1 µl per reaction. L: DNA ladder. Amplicons can be seen at the expected sizes.</p></caption><graphic xlink:href="pone.0058404.g005"></graphic></fig><p>To confirm the PCR products were HPV_SD2, the large amplicons containing one full genome shown in <xref ref-type="fig" rid="pone-0058404-g005">Figure 5B</xref> were purified then sequenced using the Sanger method. The resulting 13 sequences were aligned along the original contig generated by 454 sequencing (<xref ref-type="fig" rid="pone-0058404-g006">Figure 6A</xref>). The newly generated sequences covered 4,622 nucleotides including the complete L1 region of HPV_SD2 (<xref ref-type="fig" rid="pone-0058404-g006">Figure 6B</xref>). The total length of these new sequences represents 63.3% of the HPV_SD2 genome. A pair-wise comparison between the new sequences and the original contig (HPV_SD2) showed that the sequences newly generated sequences were identical to the contig generated with the 454 sequencing, except for 1 position, which had an insert. Twelve of 13 sequences were identical to the corresponding region of HPV_SD2 contig generated by 454 sequencing. Due to the presence of an insertion, one sequence generated by Sanger method had 99.8% identical position to HPV_SD2. Close look at this sequence suggested the 1 base insert was due to quality call in Sanger sequencing. Indeed, the additional insert would cause a frame-shift in the open reading frame of HPV_SD2, suggesting that the sequencing by 454 was correct. Overall, the confirmatory Sanger sequencing results were in agreement with our data by 99.9% (<xref ref-type="fig" rid="pone-0058404-g006">Figure 6C</xref>.). Therefore, the PCR amplification and sequencing confirm the discovery of a novel circular genome of human papillomavirus HPV_SD2.</p><fig id="pone-0058404-g006" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0058404.g006</object-id><label>Figure 6</label><caption><title>Confirmatory sequencing of full-genome amplicons generated by PCR.</title><p>A. For each amplicons, the sequencing was carried out with primers other than those used for full-genome amplification. The linear sequencing template (black line), the PCR primers are shown at both ends of the amplicons. Bleu lines (amplicons 671F-892R), green lines (amplicons 6838F-6972R), red lines (amplicons 6838F-7076R) represent the sequenced regions attached to the corresponding sequencing primers. Numbers included in primer IDs represent the position relative to HPV_SD2 genome. B. Alignment of sequences generated by Sanger method along the HPV_SD2 generated by 454 sequencing. Bars in blue, green and red correspond to sequences shown in <xref ref-type="fig" rid="pone-0058404-g006">Figure 6A</xref>. The newly generated sequences covered 4,622 nucleotides representing 63.3% of the total HPV_SD2 genome. C. Comparison of the sequences generated by Sanger method and HPV_SD2 contig generated by 454 sequencing. Size of each Sanger sequence and percentage identity are shown.</p></caption><graphic xlink:href="pone.0058404.g006"></graphic></fig></sec><sec id="s2e"><title>Conclusion</title><p>This study was carried out to develop a standard operating procedure for virus discovery using a metagenomic approach. To increase our sampling range, samples were tested in pools. This limited our ability to trace back any identified virus to an individual patient. Thus, we cannot link HPV_SD2 or any of the other viruses detected to clinical symptoms. Nevertheless, the metagenomic approach has been shown to be a powerful technique to detect and characterize new viruses that could have been missed by culture-dependent approaches or by sequence-dependent detection using probes. The metagenomic procedures and bioinformatic methods described in this study are suited for the detection of novel circular genomes of viruses, including unknown human papillomaviruses.</p><p>More HPV types are likely to be discovered as newer sequencing capabilities and bioinformatics procedures are being developed. Based on the L1 ORF, different papillomaviruses of the same genera and species share at least 60% and 90% pairwise sequence identity, respectively <xref ref-type="bibr" rid="pone.0058404-deVilliers1">[2]</xref>. According to the definition, the observed low degree of similarity and high genetic distance to previously described HPVs (E-M.de Villiers personal communication), we propose that the HPV_SD2 genome described in this study represents a novel human papillomavirus type of the Gamma genus.</p></sec></sec><sec sec-type="materials|methods" id="s3"><title>Materials and Methods</title><p>The materials and methods used in this study are outlined in <xref ref-type="fig" rid="pone-0058404-g007">Figure 7</xref>, and are detailed below.</p><fig id="pone-0058404-g007" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0058404.g007</object-id><label>Figure 7</label><caption><title>Outline of methods.</title></caption><graphic xlink:href="pone.0058404.g007"></graphic></fig><sec id="s3a"><title>Ethics Statement</title><p>Convenience samples analyzed in this study were originally collected as part of a febrile respiratory illness study conducted by the US Naval Health Research Center (NHRC). The study cohort included military recruits from various training facilities throughout the United States, residents of the US/Mexico border in and near San Diego, and military dependents. That study was approved by the NHRC Institutional review board and a signed consent was obtained from all patients prior to enrollment. Samples were transferred to SDSU as part of virus discovery collaboration between NHRC, Global Viral Forecasting, Inc. and SDSU, and included patients with pneumonia, febrile respiratory illness, and healthy controls. The metagenomic study described here was approved by SDSU’s Health Service Biological Use Authorization (BUA 06-02-062R).</p></sec><sec id="s3b"><title>Sample Preparation</title><p>Samples were collected with nasopharyngeal and oropharyngeal swabs and stored at −80°C in transport media (Remel Microtest™ M4 or Remel Microtest™ M4RT, Lenexa, KS, USA). Testing was performed in pools of 10, 35 or 75 samples using 40 µl, 100 µl and 100 µl, respectively (<xref ref-type="table" rid="pone-0058404-t002">Table 2</xref>). Pooled samples were mixed with an equal volume of sodium magnesium solution (35.8 g/L, NaCl; 2 g/L, MgSO<sub>4</sub> 7H<sub>2</sub>O; 12.5 mM TrisCl, pH 7.4) and filtered through 0.45 µm filters (Whatman, USA) to remove eukaryotic cells, bacteria and large debris. The filtrates containing viral particles, bacteria and cells that could not be removed by filtration, were treated with chloroform (20% of the total volume) to lyse any eukaryotic host cells and denature lipid bi-layers of enveloped viruses and bacteria. Samples were then treated with DNase I (0.7 units/µl of sample) to remove any free DNA. DNA and RNA were extracted concurrently using the QIAamp MinElute virus spin kit (Qiagen, cat. 51104, USA). The extracted materials containing both DNA and RNA was treated with RNase-free DNase I, then the RNA was reverse transcribed to cDNA using the Retroscript kit (Ambion, USA) (Libraries: LIB3RNA, LIB8RNA and LIB10RNA). Alternatively, the extraction of RNA for 2 libraries (LIB24RNA and LIB25RNA) was carried out using the Trizol-LS method according to manufacturer instructions, followed by cDNA synthesis using the TransPlex Complete Whole Transcriptome Amplification Kit (Sigma-Aldrich, Cat. WTA2, St. Louis, MO, USA). Four aliquots of each DNA and cDNA library were amplified by a phi29 DNA polymerase full-genome amplification method <xref ref-type="bibr" rid="pone.0058404-Dean1">[22]</xref> then pooled together to minimize preferential amplification of small circular viruses. Amplified DNA and cDNA were then purified using the NucleoSpin® Tissue Kit (Clontech, Mountain View, CA 94043, USA). Epifluorescence microscopy was used to ascertain the presence of microbes and viruses before and after chloroform treatment <xref ref-type="bibr" rid="pone.0058404-Wegley1">[23]</xref>.</p><table-wrap id="pone-0058404-t002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0058404.t002</object-id><label>Table 2</label><caption><title>Composition of sample pools and methods for DNA and RNA extraction.</title></caption><alternatives><graphic id="pone-0058404-t002-2" xlink:href="pone.0058404.t002"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Sample Pools</td><td align="left" rowspan="1" colspan="1">Samples (n)</td><td align="left" rowspan="1" colspan="1">Library ID(Type)</td><td align="left" rowspan="1" colspan="1">Extraction method</td><td align="left" rowspan="1" colspan="1">Nucleic acid amplification method</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Pool 1</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">LIB1 (DNA)</td><td align="left" rowspan="1" colspan="1">Qiamp Minelute virus spin</td><td align="left" rowspan="1" colspan="1">phi29</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">LIB3 (RNA)<xref ref-type="table-fn" rid="nt104">*</xref></td><td align="left" rowspan="1" colspan="1">Qiamp Minelute virus spin</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">Pool 2</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">LIB7 (DNA)</td><td align="left" rowspan="1" colspan="1">Qiamp Minelute virus spin</td><td align="left" rowspan="1" colspan="1">phi29</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">LIB8 (RNA)<xref ref-type="table-fn" rid="nt104">*</xref></td><td align="left" rowspan="1" colspan="1">Qiamp Minelute virus spin</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">Pool 3</td><td align="left" rowspan="1" colspan="1">75</td><td align="left" rowspan="1" colspan="1">LIB9 (DNA)</td><td align="left" rowspan="1" colspan="1">Qiamp Minelute virus spin</td><td align="left" rowspan="1" colspan="1">phi29</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">LIB10 (RNA)<xref ref-type="table-fn" rid="nt104">*</xref></td><td align="left" rowspan="1" colspan="1">Qiamp Minelute virus spin</td><td align="left" rowspan="1" colspan="1"></td></tr><tr><td align="left" rowspan="1" colspan="1">Pool 4</td><td align="left" rowspan="1" colspan="1">35</td><td align="left" rowspan="1" colspan="1">LIB24 (RNA)</td><td align="left" rowspan="1" colspan="1">Trisol-LS</td><td align="left" rowspan="1" colspan="1">WTA2</td></tr><tr><td align="left" rowspan="1" colspan="1">Pool 5</td><td align="left" rowspan="1" colspan="1">35</td><td align="left" rowspan="1" colspan="1">LIB25 (RNA)</td><td align="left" rowspan="1" colspan="1">Trisol-LS</td><td align="left" rowspan="1" colspan="1">WTA2</td></tr></tbody></table></alternatives><table-wrap-foot><fn id="nt102"><label></label><p>Phi29: Phi29 DNA polymerase full-genome amplification method <xref ref-type="bibr" rid="pone.0058404-Dean1">[22]</xref>.</p></fn><fn id="nt103"><label></label><p>WTA2: TransPlex Complete Whole Transcriptome Amplification Kit (Sigma-Aldrich, Cat. WTA2, St. Louis, MO, USA).</p></fn><fn id="nt104"><label>*</label><p>cDNA synthesis was carried out using Retroscript kit (Ambion, USA) prior to sequencing.</p></fn></table-wrap-foot></table-wrap><p>Sequencing was carried out using the Roche 454 FLX titanium according to the manufacturer’s protocol (454 Life Sciences Corporation, Roche, Brandford, CT, USA). Approximately 0.5–5 µg of purified DNA was fragmented by applying a 30 psi nitrogen pressure for 1 minute. Fragmentation was not required for the cDNA amplicons, as their sizes were within the required range for optimal pyrosequencing (150–800 bp). The fragmented DNA and cDNA ligated to 2 different oligo adaptors, each specific to 3′ or 5′ ends. The DNA/cDNA-adaptor complexes were attached to beads, which were encapsulated into individual micro-reactor micelles in which the fragments were amplified by emulsion PCR (emPCR) <xref ref-type="bibr" rid="pone.0058404-Margulies1">[24]</xref>. The clonally amplified DNA and cDNA were then loaded into a picotiter plate device in the 454 FLX sequencer to generate DNA sequences <xref ref-type="bibr" rid="pone.0058404-Margulies1">[24]</xref>.</p></sec><sec id="s3c"><title>Bioinformatics and Taxonomic Classification</title><p>Prior to data analysis, all metagenomic data was pre-processed through a pipeline of software including, TagCleaner, PRINSEQ and DeconSeq <xref ref-type="bibr" rid="pone.0058404-Schmieder1">[25]</xref>, <xref ref-type="bibr" rid="pone.0058404-Schmieder2">[26]</xref>. TagCleaner was used to remove primers added during cDNA synthesis with WTA Kit <xref ref-type="bibr" rid="pone.0058404-Tomlins1">[27]</xref>. PrinSeq was used to remove low-quality sequences and artifacts. Human DNA was filtered out of each metagenome using the DeconSeq software. Sequence reads in each metagenome were analyzed either unassembled or as assembled contigs to determine the viral and microbial relative abundance and diversity of each sample pool. The taxonomic assignment of each sequence or assembled contig was carried out using BLASTn (version 2.2.18) by comparing metagenomic sequences with the GenBank non-redundant nucleotide database <xref ref-type="bibr" rid="pone.0058404-Altschul1">[28]</xref> using a threshold E-value of ≤10<sup>−5</sup>. In the case where several related taxa yielded an equally high scoring top hit, reads were assigned to most recent common ancestor. Scores for these ancestral clades were calculated as the sum of the scores in all daughter clades <xref ref-type="bibr" rid="pone.0058404-TrindadeSilva1">[29]</xref>. Sequences with no homology in Genbank are referred to as unknowns.</p></sec><sec id="s3d"><title>Linearity and Circularity Analysis</title><p>Sequences from each metagenome with no significant hit to Genbank were assembled using the GS <italic>De Novo</italic> Assembler (454 Life Sciences Corporation, Roche, Banford, CT, USA) with requirements set at 95% identity and at least 45 bp overlaps. The assembled sequence reads and their related contigs were aligned using Geneious version 5.4 <xref ref-type="bibr" rid="pone.0058404-Drummond1">[30]</xref>. A Perl script (Document S2: circular_from_ACE.pl) was developed to determine the linearity or circularity of the unknown contigs. The analysis was performed on the ACE files generated by the <italic>De Novo</italic> Assembler.</p></sec><sec id="s3e"><title>Screening for Plasmid-specific Features</title><p>The PlasMapper <xref ref-type="bibr" rid="pone.0058404-Dong1">[13]</xref>, a web server interface that automatically generates and annotates circular plasmid maps, was used to determine if circular genomes found in this study had plasmid-specific patterns including promoters, terminators and replication origins. In addition to plasmid-specific features, PlasMapper <xref ref-type="bibr" rid="pone.0058404-Dong1">[13]</xref> and Genious Pro 5.4.3 <xref ref-type="bibr" rid="pone.0058404-Drummond1">[30]</xref> were used to identify open reading frames in circular contigs.</p></sec><sec id="s3f"><title>Human Papillomavirus Phylogeny</title><p>Taxonomical classification was performed initially with BLASTn and tBLASTx, and contigs with significant hits to human papillomavirus were aligned with complete genomes of human papilloma virus (HPV) reference sequences obtained from the Papillomavirus Episteme (PAVE) database. HPV_SD1 and HPV_SD2 from this study were aligned with PAVE sequences using ClustalW 2.1 (default parameters) and all ambiguities and gaps were removed with GBlocks 0.91 b <xref ref-type="bibr" rid="pone.0058404-Castresana1">[21]</xref>. A phylogenetic tree was constructed using PhyML 3.0.1 <xref ref-type="bibr" rid="pone.0058404-Guindon1">[31]</xref> with the following parameters: NNI tree topology search, HKY85 model of nucleotide substitution, discrete gamma model with 4 categories, estimated proportion of invariant sites and transition/transversion ratio. Bootstrap analysis using 100 resampling iterations was performed to validate the phylogenetic assignment of previously described HPV full-length genomes and to provide quality assurance of the classification of sequences. The tree was visualized with the Interactive Tree of Life (iTOL) <xref ref-type="bibr" rid="pone.0058404-Letunic1">[32]</xref>.</p></sec><sec id="s3g"><title>Confirmation by Complete Genome Amplification and Sanger Sequencing</title><p>Three polymerase chain reaction assays were performed to verify if the sequence reads assembled by the metagenomic approach derived from a complete circular double-stranded DNA in the swab sample. We designed 5 PCR primers 671F: 5′-ACCA GGCC CAAC TCCC CCAA A-3′; 892R: 5′-GACG GTCC CGCC TTTT CTTG A-3′; 7076R: <named-content content-type="gene">5′-CCGT TTGG CACT GGGG GAGG-3′</named-content>; 6838F: <named-content content-type="gene">5′-CCGC CTGC GACT CCGA AGAA-3′</named-content> and 6972R: <named-content content-type="gene">5′-GCGG CGAT GCGG TGGT AGTT-3′</named-content>, each identical to HPV_SD2 at positions 671–691, 872–892, 7057–7076, 6838–6857, 6953–7072. Three PCR reactions were designed to amplify the HPV_SD2 using primer sets 671F-892R, 6838F-7076R and 6838F-6972R (<xref ref-type="fig" rid="pone-0058404-g005">Figure 5A</xref>). The PCR amplification was performed using the Thermo Scientific Extensor Hi-Fidelity PCR Master Mix following manufacturer instruction (Thermo Fisher Scientific, Inc.). Given the large size of the amplicons and different melting temperatures of primers, we employed a touchdown approach that allowed us to amplify all the PCR products using the same conditions. The cycling conditions were as follows: 1) Denaturation, 94°C for 2 minutes; 2) denaturation, 94°C for 10 seconds; 3) Annealing, 58°C for 30 seconds; Touch-down, −0.5°C per cycle; 4) Elongation, 68°C for 8 minutes; 5) Repeat steps 2–4 28 times; 6) Elongation, 68°C for 7 minutes; store sample at 4°C until use. A 10 µl aliquot of the PCR reaction was run in a 0.5% agarose gel and any amplified amplicon was viewed by ultra violet trans-luminescence. The sizes of the predicted amplicon for these primer sets were 222, 135 and 239 bp, respectively. Assuming that HPV_SD2 is circular (see <xref ref-type="fig" rid="pone-0058404-g005">Figure 5A</xref>), each primer set (671F-892R, 6838F-7076R and 6838F-6972R) would also generate at least one larger fragment representing one or more complete genome of HPV_SD2 (7,299 bp) plus the corresponding small amplicon (222, 135, 239 base pairs, respectively) generated by the same primer set.</p><p>The resulting full-genome linear PCR amplicons were excised from the agarose gels and then purified using Zymoclean™ Gel DNA Recovery kit according to the manufacturer’s instructions. Purified amplicons were shipped on dry ice to Eton Bioscience (5820 Oberlin Drive, Suite 100; San Diego, CA 92121) for sequencing. The sequencing was carried out using the Sanger method with primers other than those used for the full-genome amplification (<xref ref-type="fig" rid="pone-0058404-g005">Figure 5</xref>).</p></sec><sec id="s3h"><title>Sequence Data</title><p>The HPV_SD2 was linearized inside the E2 gene at nucleotide 1087 of this gene (i.e. 140 nt before the end of the E2 gene) and deposited to Genbank under the accession id: KC113191. The raw metagenomic pyrosequencing reads were submitted to Sequencing Read Archive (SRA) database under the accession id: SRA051429.</p></sec></sec><sec sec-type="supplementary-material" id="s4"><title>Supporting Information</title><supplementary-material content-type="local-data" id="pone.0058404.s001"><label>Document S1</label><caption><p><bold>Results of the BLASTn analysis of all sequence reads against the Genbank database.</bold></p><p>(XLSX)</p></caption><media xlink:href="pone.0058404.s001.xlsx" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0058404.s002"><label>Document S2</label><caption><p><bold>This is a Perl script developed to determine whether a contig was circular or non-circular.</bold></p><p>(PL)</p></caption><media xlink:href="pone.0058404.s002.pl" mimetype="text" mime-subtype="plain"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0058404.s003"><label>Document S3</label><caption><p><bold>Identity matrix.</bold> Pair-wise comparison between 155 HPV full-genomes and L1region in PAVE database and 2 HPVs sequences from this study. Values shown represent percentage identity.</p><p>(XLSX)</p></caption><media xlink:href="pone.0058404.s003.xlsx" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material></sec></body><back><ack><p>The authors would like to thank Dr. Ethel-Michele de Villiers from the Reference Center for Papillomaviruses for useful advice.</p></ack><ref-list><title>References</title><ref id="pone.0058404-Bernard1"><label>1</label><mixed-citation publication-type="journal"><name><surname>Bernard</surname><given-names>HU</given-names></name>, <name><surname>Burk</surname><given-names>RD</given-names></name>, <name><surname>Chen</surname><given-names>Z</given-names></name>, <name><surname>van Doorslaer</surname><given-names>K</given-names></name>, <name><surname>Hausen</surname><given-names>H</given-names></name>, <etal>et al</etal> (<year>2010</year>) <article-title>Classification of papillomaviruses (PVs) based on 189 PV types and proposal of taxonomic amendments</article-title>. <source>Virology</source><volume>401</volume>: <fpage>70</fpage>–<lpage>79</lpage><pub-id pub-id-type="pmid">20206957</pub-id></mixed-citation></ref><ref id="pone.0058404-deVilliers1"><label>2</label><mixed-citation publication-type="journal"><name><surname>de Villiers</surname><given-names>EM</given-names></name>, <name><surname>Fauquet</surname><given-names>C</given-names></name>, <name><surname>Broker</surname><given-names>TR</given-names></name>, <name><surname>Bernard</surname><given-names>HU</given-names></name>, <name><surname>zur Hausen</surname><given-names>H</given-names></name> (<year>2004</year>) <article-title>Classification of papillomaviruses</article-title>. <source>Virology</source><volume>324</volume>: <fpage>17</fpage>–<lpage>27</lpage><pub-id pub-id-type="pmid">15183049</pub-id></mixed-citation></ref><ref id="pone.0058404-Ramamoorthy1"><label>3</label><mixed-citation publication-type="journal"><name><surname>Ramamoorthy</surname><given-names>S</given-names></name>, <name><surname>Liu</surname><given-names>YT</given-names></name>, <name><surname>Luo</surname><given-names>L</given-names></name>, <name><surname>Miyai</surname><given-names>K</given-names></name>, <name><surname>Lu</surname><given-names>Q</given-names></name>, <etal>et al</etal> (<year>2010</year>) <article-title>Detection of multiple human papillomavirus genotypes in anal carcinoma</article-title>. <source>Infect Agent Cancer</source><volume>5</volume>: <fpage>17</fpage><pub-id pub-id-type="pmid">20939896</pub-id></mixed-citation></ref><ref id="pone.0058404-Bartholomew1"><label>4</label><mixed-citation publication-type="journal"><name><surname>Bartholomew</surname><given-names>DA</given-names></name>, <name><surname>Luff</surname><given-names>RD</given-names></name>, <name><surname>Quigley</surname><given-names>NB</given-names></name>, <name><surname>Curtis</surname><given-names>M</given-names></name>, <name><surname>Olson</surname><given-names>MC</given-names></name> (<year>2011</year>) <article-title>Analytical performance of Cervista HPV 16/18 genotyping test for cervical cytology samples</article-title>. <source>J Clin Virol</source><volume>51</volume>: <fpage>38</fpage>–<lpage>43</lpage><pub-id pub-id-type="pmid">21376660</pub-id></mixed-citation></ref><ref id="pone.0058404-Chaturvedi1"><label>5</label><mixed-citation publication-type="journal"><name><surname>Chaturvedi</surname><given-names>AK</given-names></name>, <name><surname>Engels</surname><given-names>EA</given-names></name>, <name><surname>Pfeiffer</surname><given-names>RM</given-names></name>, <name><surname>Hernandez</surname><given-names>BY</given-names></name>, <name><surname>Xiao</surname><given-names>W</given-names></name>, <etal>et al</etal> (<year>2011</year>) <article-title>Human papillomavirus and rising oropharyngeal cancer incidence in the United States</article-title>. <source>J Clin Oncol</source><volume>29</volume>: <fpage>4294</fpage>–<lpage>4301</lpage><pub-id pub-id-type="pmid">21969503</pub-id></mixed-citation></ref><ref id="pone.0058404-Jemal1"><label>6</label><mixed-citation publication-type="journal"><name><surname>Jemal</surname><given-names>A</given-names></name>, <name><surname>Bray</surname><given-names>F</given-names></name> (<year>2011</year>) <collab>Center MM, Ferlay J, Ward E, et al</collab> (<year>2011</year>) <article-title>Global cancer statistics</article-title>. <source>CA Cancer J Clin</source><volume>61</volume>: <fpage>69</fpage>–<lpage>90</lpage><pub-id pub-id-type="pmid">21296855</pub-id></mixed-citation></ref><ref id="pone.0058404-Tse1"><label>7</label><mixed-citation publication-type="journal"><name><surname>Tse</surname><given-names>H</given-names></name>, <name><surname>Tsang</surname><given-names>AK</given-names></name>, <name><surname>Tsoi</surname><given-names>HW</given-names></name>, <name><surname>Leung</surname><given-names>AS</given-names></name>, <name><surname>Ho</surname><given-names>CC</given-names></name>, <etal>et al</etal> (<year>2012</year>) <article-title>Identification of a novel bat papillomavirus by metagenomics</article-title>. <source>PLoS ONE</source><volume>7</volume>: <fpage>e43986</fpage><pub-id pub-id-type="pmid">22937142</pub-id></mixed-citation></ref><ref id="pone.0058404-Forslund1"><label>8</label><mixed-citation publication-type="journal"><name><surname>Forslund</surname><given-names>O</given-names></name>, <name><surname>Ly</surname><given-names>H</given-names></name>, <name><surname>Higgins</surname><given-names>G</given-names></name> (<year>2003</year>) <article-title>Improved detection of cutaneous human papillomavirus DNA by single tube nested 'hanging droplet' PCR</article-title>. <source>J Virol Methods</source><volume>110</volume>: <fpage>129</fpage>–<lpage>136</lpage><pub-id pub-id-type="pmid">12798239</pub-id></mixed-citation></ref><ref id="pone.0058404-Riesenfeld1"><label>9</label><mixed-citation publication-type="journal"><name><surname>Riesenfeld</surname><given-names>CS</given-names></name>, <name><surname>Schloss</surname><given-names>PD</given-names></name>, <name><surname>Handelsman</surname><given-names>J</given-names></name> (<year>2004</year>) <article-title>Metagenomics: genomic analysis of microbial communities</article-title>. <source>Annu Rev Genet</source><volume>38</volume>: <fpage>525</fpage>–<lpage>552</lpage><pub-id pub-id-type="pmid">15568985</pub-id></mixed-citation></ref><ref id="pone.0058404-Edwards1"><label>10</label><mixed-citation publication-type="journal"><name><surname>Edwards</surname><given-names>RA</given-names></name>, <name><surname>Rohwer</surname><given-names>F</given-names></name> (<year>2005</year>) <article-title>Viral metagenomics</article-title>. <source>Nat Rev Microbiol</source><volume>3</volume>: <fpage>504</fpage>–<lpage>510</lpage><pub-id pub-id-type="pmid">15886693</pub-id></mixed-citation></ref><ref id="pone.0058404-Delwart1"><label>11</label><mixed-citation publication-type="journal"><name><surname>Delwart</surname><given-names>EL</given-names></name> (<year>2007</year>) <article-title>Viral metagenomics</article-title>. <source>Rev Med Virol</source><volume>17</volume>: <fpage>115</fpage>–<lpage>131</lpage><pub-id pub-id-type="pmid">17295196</pub-id></mixed-citation></ref><ref id="pone.0058404-Mokili1"><label>12</label><mixed-citation publication-type="journal"><name><surname>Mokili</surname><given-names>JL</given-names></name>, <name><surname>Rohwer</surname><given-names>F</given-names></name>, <name><surname>Dutilh</surname><given-names>BE</given-names></name> (<year>2012</year>) <article-title>Metagenomics and future perspectives in virus discovery</article-title>. <source>Curr Opin Virology</source><volume>2</volume>: <fpage>1</fpage>–<lpage>15</lpage></mixed-citation></ref><ref id="pone.0058404-Dong1"><label>13</label><mixed-citation publication-type="journal"><name><surname>Dong</surname><given-names>X</given-names></name>, <name><surname>Stothard</surname><given-names>P</given-names></name>, <name><surname>Forsythe</surname><given-names>IJ</given-names></name>, <name><surname>Wishart</surname><given-names>DS</given-names></name> (<year>2004</year>) <article-title>PlasMapper: a web server for drawing and auto-annotating plasmid maps</article-title>. <source>Nucleic Acids Res</source><volume>32</volume>: <fpage>W660</fpage>–<lpage>664</lpage><pub-id pub-id-type="pmid">15215471</pub-id></mixed-citation></ref><ref id="pone.0058404-Deutsch1"><label>14</label><mixed-citation publication-type="journal"><name><surname>Deutsch</surname><given-names>MJ</given-names></name>, <name><surname>Ott</surname><given-names>E</given-names></name>, <name><surname>Papior</surname><given-names>P</given-names></name>, <name><surname>Schepers</surname><given-names>A</given-names></name> (<year>2010</year>) <article-title>The latent origin of replication of Epstein-Barr virus directs viral genomes to active regions of the nucleus</article-title>. <source>J Virol</source><volume>84</volume>: <fpage>2533</fpage>–<lpage>2546</lpage><pub-id pub-id-type="pmid">20032186</pub-id></mixed-citation></ref><ref id="pone.0058404-Ott1"><label>15</label><mixed-citation publication-type="journal"><name><surname>Ott</surname><given-names>E</given-names></name>, <name><surname>Norio</surname><given-names>P</given-names></name>, <name><surname>Ritzi</surname><given-names>M</given-names></name>, <name><surname>Schildkraut</surname><given-names>C</given-names></name>, <name><surname>Schepers</surname><given-names>A</given-names></name> (<year>2011</year>) <article-title>The dyad symmetry element of Epstein-Barr virus is a dominant but dispensable replication origin</article-title>. <source>PLoS One</source><volume>6</volume>: <fpage>e18609</fpage><pub-id pub-id-type="pmid">21603652</pub-id></mixed-citation></ref><ref id="pone.0058404-Sarkari1"><label>16</label><mixed-citation publication-type="journal"><name><surname>Sarkari</surname><given-names>F</given-names></name>, <name><surname>Sanchez-Alcaraz</surname><given-names>T</given-names></name>, <name><surname>Wang</surname><given-names>S</given-names></name>, <name><surname>Holowaty</surname><given-names>MN</given-names></name>, <name><surname>Sheng</surname><given-names>Y</given-names></name>, <etal>et al</etal> (<year>2009</year>) <article-title>EBNA1-mediated recruitment of a histone H2B deubiquitylating complex to the Epstein-Barr virus latent origin of DNA replication</article-title>. <source>PLoS Pathog</source><volume>5</volume>: <fpage>e1000624</fpage><pub-id pub-id-type="pmid">19834552</pub-id></mixed-citation></ref><ref id="pone.0058404-Krubke1"><label>17</label><mixed-citation publication-type="journal"><name><surname>Krubke</surname><given-names>J</given-names></name>, <name><surname>Kraus</surname><given-names>J</given-names></name>, <name><surname>Delius</surname><given-names>H</given-names></name>, <name><surname>Chow</surname><given-names>L</given-names></name>, <name><surname>Broker</surname><given-names>T</given-names></name>, <etal>et al</etal> (<year>1987</year>) <article-title>Genetic relationship among human papillomaviruses associated with benign and malignant tumours of patients with epidermodysplasia verruciformis</article-title>. <source>J Gen Virol 68 (Pt</source><volume>12)</volume>: <fpage>3091</fpage>–<lpage>3103</lpage></mixed-citation></ref><ref id="pone.0058404-Androphy1"><label>18</label><mixed-citation publication-type="journal"><name><surname>Androphy</surname><given-names>EJ</given-names></name>, <name><surname>Lowy</surname><given-names>DR</given-names></name>, <name><surname>Schiller</surname><given-names>JT</given-names></name> (<year>1987</year>) <article-title>Bovine papillomavirus E2 trans-activating gene product binds to specific sites in papillomavirus DNA</article-title>. <source>Nature</source><volume>325</volume>: <fpage>70</fpage>–<lpage>73</lpage><pub-id pub-id-type="pmid">3025749</pub-id></mixed-citation></ref><ref id="pone.0058404-Spalholz1"><label>19</label><mixed-citation publication-type="journal"><name><surname>Spalholz</surname><given-names>BA</given-names></name>, <name><surname>Yang</surname><given-names>YC</given-names></name>, <name><surname>Howley</surname><given-names>PM</given-names></name> (<year>1985</year>) <article-title>Transactivation of a bovine papilloma virus transcriptional regulatory element by the E2 gene product</article-title>. <source>Cell</source><volume>42</volume>: <fpage>183</fpage>–<lpage>191</lpage><pub-id pub-id-type="pmid">2990724</pub-id></mixed-citation></ref><ref id="pone.0058404-Figge1"><label>20</label><mixed-citation publication-type="journal"><name><surname>Figge</surname><given-names>J</given-names></name>, <name><surname>Webster</surname><given-names>T</given-names></name>, <name><surname>Smith</surname><given-names>TF</given-names></name>, <name><surname>Paucha</surname><given-names>E</given-names></name> (<year>1988</year>) <article-title>Prediction of similar transforming regions in simian virus 40 large T, adenovirus E1A, and myc oncoproteins</article-title>. <source>J Virol</source><volume>62</volume>: <fpage>1814</fpage>–<lpage>1818</lpage><pub-id pub-id-type="pmid">2965766</pub-id></mixed-citation></ref><ref id="pone.0058404-Castresana1"><label>21</label><mixed-citation publication-type="journal"><name><surname>Castresana</surname><given-names>J</given-names></name> (<year>2000</year>) <article-title>Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis</article-title>. <source>Mol Biol Evol</source><volume>17</volume>: <fpage>540</fpage>–<lpage>552</lpage><pub-id pub-id-type="pmid">10742046</pub-id></mixed-citation></ref><ref id="pone.0058404-Dean1"><label>22</label><mixed-citation publication-type="journal"><name><surname>Dean</surname><given-names>FB</given-names></name>, <name><surname>Nelson</surname><given-names>JR</given-names></name>, <name><surname>Giesler</surname><given-names>TL</given-names></name>, <name><surname>Lasken</surname><given-names>RS</given-names></name> (<year>2001</year>) <article-title>Rapid amplification of plasmid and phage DNA using Phi 29 DNA polymerase and multiply-primed rolling circle amplification</article-title>. <source>Genome Res</source><volume>11</volume>: <fpage>1095</fpage>–<lpage>1099</lpage><pub-id pub-id-type="pmid">11381035</pub-id></mixed-citation></ref><ref id="pone.0058404-Wegley1"><label>23</label><mixed-citation publication-type="journal"><name><surname>Wegley</surname><given-names>L</given-names></name>, <name><surname>Mosier-Boss</surname><given-names>P</given-names></name>, <name><surname>Lieberman</surname><given-names>S</given-names></name>, <name><surname>Andrews</surname><given-names>J</given-names></name>, <name><surname>Graff-Baker</surname><given-names>A</given-names></name>, <etal>et al</etal> (<year>2006</year>) <article-title>Rapid estimation of microbial numbers in water using bulk fluorescence</article-title>. <source>Environ Microbiol</source><volume>8</volume>: <fpage>1775</fpage>–<lpage>1782</lpage><pub-id pub-id-type="pmid">16958758</pub-id></mixed-citation></ref><ref id="pone.0058404-Margulies1"><label>24</label><mixed-citation publication-type="journal"><name><surname>Margulies</surname><given-names>M</given-names></name>, <name><surname>Egholm</surname><given-names>M</given-names></name>, <name><surname>Altman</surname><given-names>WE</given-names></name>, <name><surname>Attiya</surname><given-names>S</given-names></name>, <name><surname>Bader</surname><given-names>JS</given-names></name>, <etal>et al</etal> (<year>2005</year>) <article-title>Genome sequencing in microfabricated high-density picolitre reactors</article-title>. <source>Nature</source><volume>437</volume>: <fpage>376</fpage>–<lpage>380</lpage><pub-id pub-id-type="pmid">16056220</pub-id></mixed-citation></ref><ref id="pone.0058404-Schmieder1"><label>25</label><mixed-citation publication-type="journal"><name><surname>Schmieder</surname><given-names>R</given-names></name>, <name><surname>Edwards</surname><given-names>R</given-names></name> (<year>2011</year>) <article-title>Quality control and preprocessing of metagenomic datasets</article-title>. <source>Bioinformatics</source><volume>27</volume>: <fpage>863</fpage>–<lpage>864</lpage><pub-id pub-id-type="pmid">21278185</pub-id></mixed-citation></ref><ref id="pone.0058404-Schmieder2"><label>26</label><mixed-citation publication-type="journal"><name><surname>Schmieder</surname><given-names>R</given-names></name>, <name><surname>Lim</surname><given-names>YW</given-names></name>, <name><surname>Rohwer</surname><given-names>F</given-names></name>, <name><surname>Edwards</surname><given-names>R</given-names></name> (<year>2010</year>) <article-title>TagCleaner: Identification and removal of tag sequences from genomic and metagenomic datasets</article-title>. <source>BMC Bioinformatics</source><volume>11</volume>: <fpage>341</fpage><pub-id pub-id-type="pmid">20573248</pub-id></mixed-citation></ref><ref id="pone.0058404-Tomlins1"><label>27</label><mixed-citation publication-type="journal"><name><surname>Tomlins</surname><given-names>SA</given-names></name>, <name><surname>Mehra</surname><given-names>R</given-names></name>, <name><surname>Rhodes</surname><given-names>DR</given-names></name>, <name><surname>Shah</surname><given-names>RB</given-names></name>, <name><surname>Rubin</surname><given-names>MA</given-names></name>, <etal>et al</etal> (<year>2006</year>) <article-title>Whole transcriptome amplification for gene expression profiling and development of molecular archives</article-title>. <source>Neoplasia</source><volume>8</volume>: <fpage>153</fpage>–<lpage>162</lpage><pub-id pub-id-type="pmid">16611408</pub-id></mixed-citation></ref><ref id="pone.0058404-Altschul1"><label>28</label><mixed-citation publication-type="journal"><name><surname>Altschul</surname><given-names>SF</given-names></name>, <name><surname>Gish</surname><given-names>W</given-names></name>, <name><surname>Miller</surname><given-names>W</given-names></name>, <name><surname>Myers</surname><given-names>EW</given-names></name>, <name><surname>Lipman</surname><given-names>DJ</given-names></name> (<year>1990</year>) <article-title>Basic local alignment search tool</article-title>. <source>J Mol Biol</source><volume>215</volume>: <fpage>403</fpage>–<lpage>410</lpage><pub-id pub-id-type="pmid">2231712</pub-id></mixed-citation></ref><ref id="pone.0058404-TrindadeSilva1"><label>29</label><mixed-citation publication-type="journal"><name><surname>Trindade-Silva</surname><given-names>AE</given-names></name>, <name><surname>Rua</surname><given-names>C</given-names></name>, <name><surname>Silva</surname><given-names>GG</given-names></name>, <name><surname>Dutilh</surname><given-names>BE</given-names></name>, <name><surname>Moreira</surname><given-names>AP</given-names></name>, <etal>et al</etal> (<year>2012</year>) <article-title>Taxonomic and functional microbial signatures of the endemic marine sponge Arenosclera brasiliensis</article-title>. <source>PLoS ONE</source><volume>7</volume>: <fpage>e39905</fpage><pub-id pub-id-type="pmid">22768320</pub-id></mixed-citation></ref><ref id="pone.0058404-Drummond1"><label>30</label><mixed-citation publication-type="other">Drummond AJ, Ashton B, Buxton S, Cheung M, Cooper A, <etal>et al</etal>. (2011) Geneious v5.4.</mixed-citation></ref><ref id="pone.0058404-Guindon1"><label>31</label><mixed-citation publication-type="journal"><name><surname>Guindon</surname><given-names>S</given-names></name>, <name><surname>Gascuel</surname><given-names>O</given-names></name> (<year>2003</year>) <article-title>A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood</article-title>. <source>Syst Biol</source><volume>52</volume>: <fpage>696</fpage>–<lpage>704</lpage><pub-id pub-id-type="pmid">14530136</pub-id></mixed-citation></ref><ref id="pone.0058404-Letunic1"><label>32</label><mixed-citation publication-type="other">Letunic I, Bork P (2011) Interactive Tree Of Life v2: online annotation and display of phylogenetic trees made easy. Nucleic Acids Res.</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 000129 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000129 | 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. |  |