The Role of Viral Population Diversity in Adaptation of Bovine Coronavirus to New Host Environments
Identifieur interne : 001382 ( Pmc/Corpus ); précédent : 001381; suivant : 001383The Role of Viral Population Diversity in Adaptation of Bovine Coronavirus to New Host Environments
Auteurs : Monica K. Borucki ; Jonathan E. Allen ; Haiyin Chen-Harris ; Adam Zemla ; Gilda Vanier ; Shalini Mabery ; Clinton Torres ; Pamela Hullinger ; Tom SlezakSource :
- PLoS ONE [ 1932-6203 ] ; 2013.
Abstract
The high mutation rate of RNA viruses enables a diverse genetic population of viral genotypes to exist within a single infected host. In-host genetic diversity could better position the virus population to respond and adapt to a diverse array of selective pressures such as host-switching events. Multiple new coronaviruses, including SARS, have been identified in human samples just within the last ten years, demonstrating the potential of coronaviruses as emergent human pathogens. Deep sequencing was used to characterize genomic changes in coronavirus quasispecies during simulated host-switching. Three bovine nasal samples infected with bovine coronavirus were used to infect human and bovine macrophage and lung cell lines. The virus reproduced relatively well in macrophages, but the lung cell lines were not infected efficiently enough to allow passage of non lab-adapted samples. Approximately 12 kb of the genome was amplified before and after passage and sequenced at average coverages of nearly 950×(454 sequencing) and 38,000×(Illumina). The consensus sequence of many of the passaged samples had a 12 nucleotide insert in the consensus sequence of the spike gene, and multiple point mutations were associated with the presence of the insert. Deep sequencing revealed that the insert was present but very rare in the unpassaged samples and could quickly shift to dominate the population when placed in a different environment. The insert coded for three arginine residues, occurred in a region associated with fusion entry into host cells, and may allow infection of new cell types via heparin sulfate binding. Analysis of the deep sequencing data indicated that two distinct genotypes circulated at different frequency levels in each sample, and support the hypothesis that the mutations present in passaged strains were “selected” from a pre-existing pool rather than through de novo mutation and subsequent population fixation.
Url:
DOI: 10.1371/journal.pone.0052752
PubMed: 23308119
PubMed Central: 3538757
Links to Exploration step
PMC:3538757Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">The Role of Viral Population Diversity in Adaptation of Bovine Coronavirus to New Host Environments</title>
<author><name sortKey="Borucki, Monica K" sort="Borucki, Monica K" uniqKey="Borucki M" first="Monica K." last="Borucki">Monica K. Borucki</name>
<affiliation><nlm:aff id="aff1"></nlm:aff>
</affiliation>
</author>
<author><name sortKey="Allen, Jonathan E" sort="Allen, Jonathan E" uniqKey="Allen J" first="Jonathan E." last="Allen">Jonathan E. Allen</name>
<affiliation><nlm:aff id="aff1"></nlm:aff>
</affiliation>
</author>
<author><name sortKey="Chen Harris, Haiyin" sort="Chen Harris, Haiyin" uniqKey="Chen Harris H" first="Haiyin" last="Chen-Harris">Haiyin Chen-Harris</name>
<affiliation><nlm:aff id="aff1"></nlm:aff>
</affiliation>
</author>
<author><name sortKey="Zemla, Adam" sort="Zemla, Adam" uniqKey="Zemla A" first="Adam" last="Zemla">Adam Zemla</name>
<affiliation><nlm:aff id="aff1"></nlm:aff>
</affiliation>
</author>
<author><name sortKey="Vanier, Gilda" sort="Vanier, Gilda" uniqKey="Vanier G" first="Gilda" last="Vanier">Gilda Vanier</name>
<affiliation><nlm:aff id="aff1"></nlm:aff>
</affiliation>
</author>
<author><name sortKey="Mabery, Shalini" sort="Mabery, Shalini" uniqKey="Mabery S" first="Shalini" last="Mabery">Shalini Mabery</name>
<affiliation><nlm:aff id="aff1"></nlm:aff>
</affiliation>
</author>
<author><name sortKey="Torres, Clinton" sort="Torres, Clinton" uniqKey="Torres C" first="Clinton" last="Torres">Clinton Torres</name>
<affiliation><nlm:aff id="aff1"></nlm:aff>
</affiliation>
</author>
<author><name sortKey="Hullinger, Pamela" sort="Hullinger, Pamela" uniqKey="Hullinger P" first="Pamela" last="Hullinger">Pamela Hullinger</name>
<affiliation><nlm:aff id="aff1"></nlm:aff>
</affiliation>
</author>
<author><name sortKey="Slezak, Tom" sort="Slezak, Tom" uniqKey="Slezak T" first="Tom" last="Slezak">Tom Slezak</name>
<affiliation><nlm:aff id="aff1"></nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">23308119</idno>
<idno type="pmc">3538757</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3538757</idno>
<idno type="RBID">PMC:3538757</idno>
<idno type="doi">10.1371/journal.pone.0052752</idno>
<date when="2013">2013</date>
<idno type="wicri:Area/Pmc/Corpus">001382</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">001382</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">The Role of Viral Population Diversity in Adaptation of Bovine Coronavirus to New Host Environments</title>
<author><name sortKey="Borucki, Monica K" sort="Borucki, Monica K" uniqKey="Borucki M" first="Monica K." last="Borucki">Monica K. Borucki</name>
<affiliation><nlm:aff id="aff1"></nlm:aff>
</affiliation>
</author>
<author><name sortKey="Allen, Jonathan E" sort="Allen, Jonathan E" uniqKey="Allen J" first="Jonathan E." last="Allen">Jonathan E. Allen</name>
<affiliation><nlm:aff id="aff1"></nlm:aff>
</affiliation>
</author>
<author><name sortKey="Chen Harris, Haiyin" sort="Chen Harris, Haiyin" uniqKey="Chen Harris H" first="Haiyin" last="Chen-Harris">Haiyin Chen-Harris</name>
<affiliation><nlm:aff id="aff1"></nlm:aff>
</affiliation>
</author>
<author><name sortKey="Zemla, Adam" sort="Zemla, Adam" uniqKey="Zemla A" first="Adam" last="Zemla">Adam Zemla</name>
<affiliation><nlm:aff id="aff1"></nlm:aff>
</affiliation>
</author>
<author><name sortKey="Vanier, Gilda" sort="Vanier, Gilda" uniqKey="Vanier G" first="Gilda" last="Vanier">Gilda Vanier</name>
<affiliation><nlm:aff id="aff1"></nlm:aff>
</affiliation>
</author>
<author><name sortKey="Mabery, Shalini" sort="Mabery, Shalini" uniqKey="Mabery S" first="Shalini" last="Mabery">Shalini Mabery</name>
<affiliation><nlm:aff id="aff1"></nlm:aff>
</affiliation>
</author>
<author><name sortKey="Torres, Clinton" sort="Torres, Clinton" uniqKey="Torres C" first="Clinton" last="Torres">Clinton Torres</name>
<affiliation><nlm:aff id="aff1"></nlm:aff>
</affiliation>
</author>
<author><name sortKey="Hullinger, Pamela" sort="Hullinger, Pamela" uniqKey="Hullinger P" first="Pamela" last="Hullinger">Pamela Hullinger</name>
<affiliation><nlm:aff id="aff1"></nlm:aff>
</affiliation>
</author>
<author><name sortKey="Slezak, Tom" sort="Slezak, Tom" uniqKey="Slezak T" first="Tom" last="Slezak">Tom Slezak</name>
<affiliation><nlm:aff id="aff1"></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>The high mutation rate of RNA viruses enables a diverse genetic population of viral genotypes to exist within a single infected host. In-host genetic diversity could better position the virus population to respond and adapt to a diverse array of selective pressures such as host-switching events. Multiple new coronaviruses, including SARS, have been identified in human samples just within the last ten years, demonstrating the potential of coronaviruses as emergent human pathogens. Deep sequencing was used to characterize genomic changes in coronavirus quasispecies during simulated host-switching. Three bovine nasal samples infected with bovine coronavirus were used to infect human and bovine macrophage and lung cell lines. The virus reproduced relatively well in macrophages, but the lung cell lines were not infected efficiently enough to allow passage of non lab-adapted samples. Approximately 12 kb of the genome was amplified before and after passage and sequenced at average coverages of nearly 950×(454 sequencing) and 38,000×(Illumina). The consensus sequence of many of the passaged samples had a 12 nucleotide insert in the consensus sequence of the spike gene, and multiple point mutations were associated with the presence of the insert. Deep sequencing revealed that the insert was present but very rare in the unpassaged samples and could quickly shift to dominate the population when placed in a different environment. The insert coded for three arginine residues, occurred in a region associated with fusion entry into host cells, and may allow infection of new cell types via heparin sulfate binding. Analysis of the deep sequencing data indicated that two distinct genotypes circulated at different frequency levels in each sample, and support the hypothesis that the mutations present in passaged strains were “selected” from a pre-existing pool rather than through de novo mutation and subsequent population fixation.</p>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Woolhouse, Mej" uniqKey="Woolhouse M">MEJ Woolhouse</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Elena, Sf" uniqKey="Elena S">SF Elena</name>
</author>
<author><name sortKey="Carrasco, P" uniqKey="Carrasco P">P Carrasco</name>
</author>
<author><name sortKey="Daros, J A" uniqKey="Daros J">J-A Daros</name>
</author>
<author><name sortKey="Sanjuan, R" uniqKey="Sanjuan R">R Sanjuan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Domingo, E" uniqKey="Domingo E">E Domingo</name>
</author>
<author><name sortKey="Sheldon, J" uniqKey="Sheldon J">J Sheldon</name>
</author>
<author><name sortKey="Perales, C" uniqKey="Perales C">C Perales</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Borderia, Av" uniqKey="Borderia A">AV Borderia</name>
</author>
<author><name sortKey="Stapleford, Ka" uniqKey="Stapleford K">KA Stapleford</name>
</author>
<author><name sortKey="Vignuzzi, M" uniqKey="Vignuzzi M">M Vignuzzi</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Graham, Rl" uniqKey="Graham R">RL Graham</name>
</author>
<author><name sortKey="Baric, Rs" uniqKey="Baric R">RS Baric</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Vijgen, L" uniqKey="Vijgen L">L Vijgen</name>
</author>
<author><name sortKey="Keyaerts, E" uniqKey="Keyaerts E">E Keyaerts</name>
</author>
<author><name sortKey="Lemey, P" uniqKey="Lemey P">P Lemey</name>
</author>
<author><name sortKey="Maes, P" uniqKey="Maes P">P Maes</name>
</author>
<author><name sortKey="Van Reeth, K" uniqKey="Van Reeth K">K Van Reeth</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Vijgen, L" uniqKey="Vijgen L">L Vijgen</name>
</author>
<author><name sortKey="Keyaerts, E" uniqKey="Keyaerts E">E Keyaerts</name>
</author>
<author><name sortKey="Moes, E" uniqKey="Moes E">E Moes</name>
</author>
<author><name sortKey="Thoelen, I" uniqKey="Thoelen I">I Thoelen</name>
</author>
<author><name sortKey="Wollants, E" uniqKey="Wollants E">E Wollants</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gonzalez, Jm" uniqKey="Gonzalez J">JM González</name>
</author>
<author><name sortKey="Gomez Puertas, P" uniqKey="Gomez Puertas P">P Gomez-Puertas</name>
</author>
<author><name sortKey="Cavanagh, D" uniqKey="Cavanagh D">D Cavanagh</name>
</author>
<author><name sortKey="Gorbalenya, Ae" uniqKey="Gorbalenya A">AE Gorbalenya</name>
</author>
<author><name sortKey="Enjuanes, L" uniqKey="Enjuanes L">L Enjuanes</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Han, Mg" uniqKey="Han M">MG Han</name>
</author>
<author><name sortKey="Cheon, D S" uniqKey="Cheon D">D-S Cheon</name>
</author>
<author><name sortKey="Zhang, X" uniqKey="Zhang X">X Zhang</name>
</author>
<author><name sortKey="Saif, Lj" uniqKey="Saif L">LJ Saif</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Schultze, B" uniqKey="Schultze B">B Schultze</name>
</author>
<author><name sortKey="Herrler, G" uniqKey="Herrler G">G Herrler</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Li, W" uniqKey="Li W">W Li</name>
</author>
<author><name sortKey="Wong, S K" uniqKey="Wong S">S-K Wong</name>
</author>
<author><name sortKey="Li, F" uniqKey="Li F">F Li</name>
</author>
<author><name sortKey="Kuhn, Jh" uniqKey="Kuhn J">JH Kuhn</name>
</author>
<author><name sortKey="Huang, Ic" uniqKey="Huang I">IC Huang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Peng, G" uniqKey="Peng G">G Peng</name>
</author>
<author><name sortKey="Sun, D" uniqKey="Sun D">D Sun</name>
</author>
<author><name sortKey="Rajashankar, Kr" uniqKey="Rajashankar K">KR Rajashankar</name>
</author>
<author><name sortKey="Qian, Z" uniqKey="Qian Z">Z Qian</name>
</author>
<author><name sortKey="Holmes, Kv" uniqKey="Holmes K">KV Holmes</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zhang, X" uniqKey="Zhang X">X Zhang</name>
</author>
<author><name sortKey="Hasoksuz, M" uniqKey="Hasoksuz M">M Hasoksuz</name>
</author>
<author><name sortKey="Spiro, D" uniqKey="Spiro D">D Spiro</name>
</author>
<author><name sortKey="Halpin, R" uniqKey="Halpin R">R Halpin</name>
</author>
<author><name sortKey="Wang, S" uniqKey="Wang S">S Wang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chouljenko, Vn" uniqKey="Chouljenko V">VN Chouljenko</name>
</author>
<author><name sortKey="Lin, Xq" uniqKey="Lin X">XQ Lin</name>
</author>
<author><name sortKey="Storz, J" uniqKey="Storz J">J Storz</name>
</author>
<author><name sortKey="Kousoulas, Kg" uniqKey="Kousoulas K">KG Kousoulas</name>
</author>
<author><name sortKey="Gorbalenya, Ae" uniqKey="Gorbalenya A">AE Gorbalenya</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Perlman, S" uniqKey="Perlman S">S Perlman</name>
</author>
<author><name sortKey="Netland, J" uniqKey="Netland J">J Netland</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Devaraj, Sg" uniqKey="Devaraj S">SG Devaraj</name>
</author>
<author><name sortKey="Wang, N" uniqKey="Wang N">N Wang</name>
</author>
<author><name sortKey="Chen, Z" uniqKey="Chen Z">Z Chen</name>
</author>
<author><name sortKey="Chen, Z" uniqKey="Chen Z">Z Chen</name>
</author>
<author><name sortKey="Tseng, M" uniqKey="Tseng M">M Tseng</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sun, L" uniqKey="Sun L">L Sun</name>
</author>
<author><name sortKey="Xing, Y" uniqKey="Xing Y">Y Xing</name>
</author>
<author><name sortKey="Chen, X" uniqKey="Chen X">X Chen</name>
</author>
<author><name sortKey="Zheng, Y" uniqKey="Zheng Y">Y Zheng</name>
</author>
<author><name sortKey="Yang, Y" uniqKey="Yang Y">Y Yang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Almeida, Ms" uniqKey="Almeida M">MS Almeida</name>
</author>
<author><name sortKey="Johnson, Ma" uniqKey="Johnson M">MA Johnson</name>
</author>
<author><name sortKey="Herrmann, T" uniqKey="Herrmann T">T Herrmann</name>
</author>
<author><name sortKey="Geralt, M" uniqKey="Geralt M">M Geralt</name>
</author>
<author><name sortKey="Wuthrich, K" uniqKey="Wuthrich K">K Wüthrich</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Brockway, Sm" uniqKey="Brockway S">SM Brockway</name>
</author>
<author><name sortKey="Lu, Xt" uniqKey="Lu X">XT Lu</name>
</author>
<author><name sortKey="Peters, Tr" uniqKey="Peters T">TR Peters</name>
</author>
<author><name sortKey="Dermody, Ts" uniqKey="Dermody T">TS Dermody</name>
</author>
<author><name sortKey="Denison, Mr" uniqKey="Denison M">MR Denison</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Saif, Lj" uniqKey="Saif L">LJ Saif</name>
</author>
<author><name sortKey="Redman, Dr" uniqKey="Redman D">DR Redman</name>
</author>
<author><name sortKey="Brock, Kv" uniqKey="Brock K">KV Brock</name>
</author>
<author><name sortKey="Kohler, Em" uniqKey="Kohler E">EM Kohler</name>
</author>
<author><name sortKey="Heckert, Ra" uniqKey="Heckert R">RA Heckert</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="St Cyr Coats, Ks" uniqKey="St Cyr Coats K">KS St Cyr-Coats</name>
</author>
<author><name sortKey="Storz, J" uniqKey="Storz J">J Storz</name>
</author>
<author><name sortKey="Hussain, Ka" uniqKey="Hussain K">KA Hussain</name>
</author>
<author><name sortKey="Schnorr, Kl" uniqKey="Schnorr K">KL Schnorr</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Parker, Md" uniqKey="Parker M">MD Parker</name>
</author>
<author><name sortKey="Yoo, D" uniqKey="Yoo D">D Yoo</name>
</author>
<author><name sortKey="Cox, Gj" uniqKey="Cox G">GJ Cox</name>
</author>
<author><name sortKey="Babiuk, La" uniqKey="Babiuk L">LA Babiuk</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Belouzard, S" uniqKey="Belouzard S">S Belouzard</name>
</author>
<author><name sortKey="Chu, Vc" uniqKey="Chu V">VC Chu</name>
</author>
<author><name sortKey="Whittaker, Gr" uniqKey="Whittaker G">GR Whittaker</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zemla, A" uniqKey="Zemla A">A Zemla</name>
</author>
<author><name sortKey="Zhou, Ce" uniqKey="Zhou C">CE Zhou</name>
</author>
<author><name sortKey="Slezak, T" uniqKey="Slezak T">T Slezak</name>
</author>
<author><name sortKey="Kuczmarski, T" uniqKey="Kuczmarski T">T Kuczmarski</name>
</author>
<author><name sortKey="Rama, D" uniqKey="Rama D">D Rama</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zagordi, O" uniqKey="Zagordi O">O Zagordi</name>
</author>
<author><name sortKey="Bhattacharya, A" uniqKey="Bhattacharya A">A Bhattacharya</name>
</author>
<author><name sortKey="Eriksson, N" uniqKey="Eriksson N">N Eriksson</name>
</author>
<author><name sortKey="Beerenwinkel, N" uniqKey="Beerenwinkel N">N Beerenwinkel</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zagordi, O" uniqKey="Zagordi O">O Zagordi</name>
</author>
<author><name sortKey="Klein, R" uniqKey="Klein R">R Klein</name>
</author>
<author><name sortKey="Daumer, M" uniqKey="Daumer M">M Daumer</name>
</author>
<author><name sortKey="Beerenwinkel, N" uniqKey="Beerenwinkel N">N Beerenwinkel</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wright, Cf" uniqKey="Wright C">CF Wright</name>
</author>
<author><name sortKey="Morelli, Mj" uniqKey="Morelli M">MJ Morelli</name>
</author>
<author><name sortKey="Thebaud, G" uniqKey="Thebaud G">G Thébaud</name>
</author>
<author><name sortKey="Knowles, Nj" uniqKey="Knowles N">NJ Knowles</name>
</author>
<author><name sortKey="Herzyk, P" uniqKey="Herzyk P">P Herzyk</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hofmann, Ma" uniqKey="Hofmann M">MA Hofmann</name>
</author>
<author><name sortKey="Sethna, Pb" uniqKey="Sethna P">PB Sethna</name>
</author>
<author><name sortKey="Brian, Da" uniqKey="Brian D">DA Brian</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Akaike, T" uniqKey="Akaike T">T Akaike</name>
</author>
<author><name sortKey="Fujii, S" uniqKey="Fujii S">S Fujii</name>
</author>
<author><name sortKey="Kato, A" uniqKey="Kato A">A Kato</name>
</author>
<author><name sortKey="Yoshitake, J" uniqKey="Yoshitake J">J Yoshitake</name>
</author>
<author><name sortKey="Miyamoto, Y" uniqKey="Miyamoto Y">Y Miyamoto</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Watanabe, R" uniqKey="Watanabe R">R Watanabe</name>
</author>
<author><name sortKey="Matsuyama, S" uniqKey="Matsuyama S">S Matsuyama</name>
</author>
<author><name sortKey="Shirato, K" uniqKey="Shirato K">K Shirato</name>
</author>
<author><name sortKey="Maejima, M" uniqKey="Maejima M">M Maejima</name>
</author>
<author><name sortKey="Fukushi, S" uniqKey="Fukushi S">S Fukushi</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="De Haan, Cam" uniqKey="De Haan C">CAM de Haan</name>
</author>
<author><name sortKey="Haijema, Bj" uniqKey="Haijema B">BJ Haijema</name>
</author>
<author><name sortKey="Schellen, P" uniqKey="Schellen P">P Schellen</name>
</author>
<author><name sortKey="Schreur, Pw" uniqKey="Schreur P">PW Schreur</name>
</author>
<author><name sortKey="Te Lintelo, E" uniqKey="Te Lintelo E">E te Lintelo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bernard, Ka" uniqKey="Bernard K">KA Bernard</name>
</author>
<author><name sortKey="Klimstra, Wb" uniqKey="Klimstra W">WB Klimstra</name>
</author>
<author><name sortKey="Johnston, Re" uniqKey="Johnston R">RE Johnston</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="De Haan, Cam" uniqKey="De Haan C">CAM de Haan</name>
</author>
<author><name sortKey="Li, Z" uniqKey="Li Z">Z Li</name>
</author>
<author><name sortKey="Te Lintelo, E" uniqKey="Te Lintelo E">E te Lintelo</name>
</author>
<author><name sortKey="Bosch, Bj" uniqKey="Bosch B">BJ Bosch</name>
</author>
<author><name sortKey="Haijema, Bj" uniqKey="Haijema B">BJ Haijema</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Spillmann, D" uniqKey="Spillmann D">D Spillmann</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jacques, A" uniqKey="Jacques A">A Jacques</name>
</author>
<author><name sortKey="Bleau, C" uniqKey="Bleau C">C Bleau</name>
</author>
<author><name sortKey="Turbide, C" uniqKey="Turbide C">C Turbide</name>
</author>
<author><name sortKey="Beauchemin, N" uniqKey="Beauchemin N">N Beauchemin</name>
</author>
<author><name sortKey="Lamontagne, L" uniqKey="Lamontagne L">L Lamontagne</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bertram, S" uniqKey="Bertram S">S Bertram</name>
</author>
<author><name sortKey="Glowacka, I" uniqKey="Glowacka I">I Glowacka</name>
</author>
<author><name sortKey="Muller, Ma" uniqKey="Muller M">MA Müller</name>
</author>
<author><name sortKey="Lavender, H" uniqKey="Lavender H">H Lavender</name>
</author>
<author><name sortKey="Gnirss, K" uniqKey="Gnirss K">K Gnirss</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chouljenko, Vn" uniqKey="Chouljenko V">VN Chouljenko</name>
</author>
<author><name sortKey="Kousoulas, Kg" uniqKey="Kousoulas K">KG Kousoulas</name>
</author>
<author><name sortKey="Lin, X" uniqKey="Lin X">X Lin</name>
</author>
<author><name sortKey="Storz, J" uniqKey="Storz J">J Storz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zhang, Xm" uniqKey="Zhang X">XM Zhang</name>
</author>
<author><name sortKey="Kousoulas, Kg" uniqKey="Kousoulas K">KG Kousoulas</name>
</author>
<author><name sortKey="Storz, J" uniqKey="Storz J">J Storz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mcroy, Wc" uniqKey="Mcroy W">WC McRoy</name>
</author>
<author><name sortKey="Baric, Rs" uniqKey="Baric R">RS Baric</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Stabel, Jr" uniqKey="Stabel J">JR Stabel</name>
</author>
<author><name sortKey="Stabel, Tj" uniqKey="Stabel T">TJ Stabel</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cho, Ko" uniqKey="Cho K">KO Cho</name>
</author>
<author><name sortKey="Hasoksuz, M" uniqKey="Hasoksuz M">M Hasoksuz</name>
</author>
<author><name sortKey="Nielsen, Pr" uniqKey="Nielsen P">PR Nielsen</name>
</author>
<author><name sortKey="Chang, Ko" uniqKey="Chang K">KO Chang</name>
</author>
<author><name sortKey="Lathrop, S" uniqKey="Lathrop S">S Lathrop</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tamura, K" uniqKey="Tamura K">K Tamura</name>
</author>
<author><name sortKey="Nei, M" uniqKey="Nei M">M Nei</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tamura, K" uniqKey="Tamura K">K Tamura</name>
</author>
<author><name sortKey="Peterson, D" uniqKey="Peterson D">D Peterson</name>
</author>
<author><name sortKey="Peterson, N" uniqKey="Peterson N">N Peterson</name>
</author>
<author><name sortKey="Stecher, G" uniqKey="Stecher G">G Stecher</name>
</author>
<author><name sortKey="Nei, M" uniqKey="Nei M">M Nei</name>
</author>
</analytic>
</biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">PLoS One</journal-id>
<journal-id journal-id-type="iso-abbrev">PLoS ONE</journal-id>
<journal-id journal-id-type="publisher-id">plos</journal-id>
<journal-id journal-id-type="pmc">plosone</journal-id>
<journal-title-group><journal-title>PLoS ONE</journal-title>
</journal-title-group>
<issn pub-type="epub">1932-6203</issn>
<publisher><publisher-name>Public Library of Science</publisher-name>
<publisher-loc>San Francisco, USA</publisher-loc>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">23308119</article-id>
<article-id pub-id-type="pmc">3538757</article-id>
<article-id pub-id-type="publisher-id">PONE-D-12-29217</article-id>
<article-id pub-id-type="doi">10.1371/journal.pone.0052752</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>Genome Evolution</subject>
<subject>Genome Sequencing</subject>
</subj-group>
</subj-group>
<subj-group><subject>Microbiology</subject>
<subj-group><subject>Virology</subject>
<subj-group><subject>Viral classification</subject>
<subj-group><subject>RNA viruses</subject>
</subj-group>
</subj-group>
<subj-group><subject>Emerging Viral Diseases</subject>
<subject>Viral Evolution</subject>
<subject>Viral Transmission and Infection</subject>
</subj-group>
</subj-group>
<subj-group><subject>Emerging Infectious Diseases</subject>
<subject>Host-Pathogen Interaction</subject>
</subj-group>
</subj-group>
</subj-group>
<subj-group subj-group-type="Discipline-v2"><subject>Computer Science</subject>
<subj-group><subject>Computing Methods</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group><article-title>The Role of Viral Population Diversity in Adaptation of Bovine Coronavirus to New Host Environments</article-title>
<alt-title alt-title-type="running-head">Role of Quasispecies in Coronavirus Evolution</alt-title>
</title-group>
<contrib-group><contrib contrib-type="author"><name><surname>Borucki</surname>
<given-names>Monica K.</given-names>
</name>
<xref ref-type="aff" rid="aff1"></xref>
<xref ref-type="corresp" rid="cor1"><sup>*</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Allen</surname>
<given-names>Jonathan E.</given-names>
</name>
<xref ref-type="aff" rid="aff1"></xref>
</contrib>
<contrib contrib-type="author"><name><surname>Chen-Harris</surname>
<given-names>Haiyin</given-names>
</name>
<xref ref-type="aff" rid="aff1"></xref>
</contrib>
<contrib contrib-type="author"><name><surname>Zemla</surname>
<given-names>Adam</given-names>
</name>
<xref ref-type="aff" rid="aff1"></xref>
</contrib>
<contrib contrib-type="author"><name><surname>Vanier</surname>
<given-names>Gilda</given-names>
</name>
<xref ref-type="aff" rid="aff1"></xref>
</contrib>
<contrib contrib-type="author"><name><surname>Mabery</surname>
<given-names>Shalini</given-names>
</name>
<xref ref-type="aff" rid="aff1"></xref>
</contrib>
<contrib contrib-type="author"><name><surname>Torres</surname>
<given-names>Clinton</given-names>
</name>
<xref ref-type="aff" rid="aff1"></xref>
</contrib>
<contrib contrib-type="author"><name><surname>Hullinger</surname>
<given-names>Pamela</given-names>
</name>
<xref ref-type="aff" rid="aff1"></xref>
</contrib>
<contrib contrib-type="author"><name><surname>Slezak</surname>
<given-names>Tom</given-names>
</name>
<xref ref-type="aff" rid="aff1"></xref>
</contrib>
</contrib-group>
<aff id="aff1"><addr-line>Lawrence Livermore National Laboratory, Livermore, California, United States of America</addr-line>
</aff>
<contrib-group><contrib contrib-type="editor"><name><surname>Coffey</surname>
<given-names>Lark L.</given-names>
</name>
<role>Editor</role>
<xref ref-type="aff" rid="edit1"></xref>
</contrib>
</contrib-group>
<aff id="edit1"><addr-line>Blood Systems Research Institute, United States of America</addr-line>
</aff>
<author-notes><corresp id="cor1">* E-mail: <email>borucki2@llnl.gov</email>
</corresp>
<fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>
The authors have declared that no competing interests exist.</p>
</fn>
<fn fn-type="con"><p>Conceived and designed the experiments: MKB JEA TS. Performed the experiments: GV SM AZ CT. Analyzed the data: MKB JEA HC. Contributed reagents/materials/analysis tools: PH TS. Wrote the paper: MKB HC JEA.</p>
</fn>
</author-notes>
<pub-date pub-type="collection"><year>2013</year>
</pub-date>
<pub-date pub-type="epub"><day>7</day>
<month>1</month>
<year>2013</year>
</pub-date>
<volume>8</volume>
<issue>1</issue>
<elocation-id>e52752</elocation-id>
<history><date date-type="received"><day>25</day>
<month>9</month>
<year>2012</year>
</date>
<date date-type="accepted"><day>21</day>
<month>11</month>
<year>2012</year>
</date>
</history>
<permissions><copyright-year>2013</copyright-year>
<license xlink:href="https://creativecommons.org/publicdomain/zero/1.0/"><license-p>This is an open-access article distributed under the terms of the Creative Commons Public Domain declaration, which stipulates that, once placed in the public domain, this work may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose.</license-p>
</license>
</permissions>
<abstract><p>The high mutation rate of RNA viruses enables a diverse genetic population of viral genotypes to exist within a single infected host. In-host genetic diversity could better position the virus population to respond and adapt to a diverse array of selective pressures such as host-switching events. Multiple new coronaviruses, including SARS, have been identified in human samples just within the last ten years, demonstrating the potential of coronaviruses as emergent human pathogens. Deep sequencing was used to characterize genomic changes in coronavirus quasispecies during simulated host-switching. Three bovine nasal samples infected with bovine coronavirus were used to infect human and bovine macrophage and lung cell lines. The virus reproduced relatively well in macrophages, but the lung cell lines were not infected efficiently enough to allow passage of non lab-adapted samples. Approximately 12 kb of the genome was amplified before and after passage and sequenced at average coverages of nearly 950×(454 sequencing) and 38,000×(Illumina). The consensus sequence of many of the passaged samples had a 12 nucleotide insert in the consensus sequence of the spike gene, and multiple point mutations were associated with the presence of the insert. Deep sequencing revealed that the insert was present but very rare in the unpassaged samples and could quickly shift to dominate the population when placed in a different environment. The insert coded for three arginine residues, occurred in a region associated with fusion entry into host cells, and may allow infection of new cell types via heparin sulfate binding. Analysis of the deep sequencing data indicated that two distinct genotypes circulated at different frequency levels in each sample, and support the hypothesis that the mutations present in passaged strains were “selected” from a pre-existing pool rather than through de novo mutation and subsequent population fixation.</p>
</abstract>
<funding-group><funding-statement>Funding was provided by Lawrence Livermore National Laboratory internal funding (LDRD 10-LW-020) and a grant from the Defense Threat Reduction Agency, Transformational Medical Technologies Initiative. 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="11"></page-count>
</counts>
</article-meta>
</front>
<body><sec id="s1"><title>Introduction</title>
<p>Three quarters of the recently discovered human pathogens are viral, and most of those are RNA viruses <xref rid="pone.0052752-Woolhouse1" ref-type="bibr">[1]</xref>
. Some of these emergent viruses, such as HIV and SARS coronavirus (SARS-CoV), are capable of causing epidemics of human disease. RNA virus populations sustain high genetic diversity due to the low fidelity of their polymerase, short genome, high replication rates and large population size <xref rid="pone.0052752-Elena1" ref-type="bibr">[2]</xref>
. For this reason a single RNA virus population can consist of a multiplicity of slightly different genomes, sometimes referred to as a mutant spectra <xref rid="pone.0052752-Domingo1" ref-type="bibr">[3]</xref>
. The high mutation rate of RNA viruses increases the ability of these viruses to adapt to diverse hosts (interspecies transmission events) and the potential cause new human and zoonotic diseases <xref rid="pone.0052752-Borderia1" ref-type="bibr">[4]</xref>
, however, very little is known about the particular mutations that enable interspecies transmission events to occur.</p>
<p>Coronaviruses are particularly adept at adapting to new hosts due in part to their amazing capacity for genome recombination. Coronaviruses have the largest genome of RNA viruses, consisting of 27–30 kb positive sense single-stranded RNA. Although recombination can lead to an interspecies transmission event, as was believed to be the case with SARS-CoV, accumulation of point mutations may also enable the coronaviruses to adapt to new host species <xref rid="pone.0052752-Graham1" ref-type="bibr">[5]</xref>
–<xref rid="pone.0052752-Vijgen2" ref-type="bibr">[7]</xref>
.</p>
<p>The <italic>Coronaviridae</italic>
subfamily <italic>Coronovirinae</italic>
is composed of three genera based on serologic and genetic characteristics: <italic>Alphacoronavirus</italic>
(formerly Group 1) includes viruses that infect pigs, dogs, cats and humans; <italic>Betacoronavirus</italic>
(formerly Group 2) includes bovine, bat, human, horse, pig, rodent, and bat viruses; and <italic>Gammacoronavirus</italic>
(formerly Group 3) which consists of viruses adapted to birds <xref rid="pone.0052752-Vijgen2" ref-type="bibr">[7]</xref>
, <xref rid="pone.0052752-Gonzlez1" ref-type="bibr">[8]</xref>
. Bovine coronavirus (BCoV) is a betacoronavirus which is related to SARS-CoV and has caused disease in humans on at least one occasion <xref rid="pone.0052752-Han1" ref-type="bibr">[9]</xref>
. BCoV is known to use 9-<italic>O-</italic>
acetylated sialic acid to bind to host cells, although a second receptor may be involved <xref rid="pone.0052752-Schultze1" ref-type="bibr">[10]</xref>
. The spike protein which is present on the surface of the virion, determines host range and tissue tropism of coronvaviruses. The receptor binding domain of BCoV has not been determined <xref rid="pone.0052752-Li1" ref-type="bibr">[11]</xref>
, but a recent study by Peng at al. (2011) indicates that it falls within the N terminal domain <xref rid="pone.0052752-Peng1" ref-type="bibr">[12]</xref>
. BCoV infection may cause acute and severe diarrhea and respiratory symptoms in cattle especially under stressful conditions such as transport, but subclinical infections may occur in healthy cattle. BCoV isolates are generally obtained by inoculating nasal or fecal samples from infected cattle in human rectal tumor cells (HRT-18 cells), however the genetic changes in the viral population that allow this bovine virus to adapt to human cell lines have yet to be defined <xref rid="pone.0052752-Zhang1" ref-type="bibr">[13]</xref>
.</p>
<p>We are interested in understanding the role of natural viral population diversity in the adaptation of BCoV to new host environments; in particular, cell types that may enable the virus to spread via the respiratory route or to cause systemic disease. The mutational dynamics of BCoV was determined by serial passage of BCoV nasal samples in human and bovine lung and macrophage cell lines as well as human enteric cells (A549, EBL, THP-1, Bomac, and HRT-18, respectively). The consensus and subconsensus (variant) nucleotide sequences of cell culture passaged samples were compared to that of the natural “unpassaged” viral populations. The genome regions analyzed included nsp1and nsp3 (genes involved with evasion of the host innate immune response and other functions such viral replication), nsp14 (designated nsp11 in NCBI BCoV Reference Sequence: NC_003045.1; involved with polymerase fidelity), and the spike protein gene (determinant of host range) <xref rid="pone.0052752-Vijgen2" ref-type="bibr">[7]</xref>
, <xref rid="pone.0052752-Chouljenko1" ref-type="bibr">[14]</xref>
–<xref rid="pone.0052752-Brockway1" ref-type="bibr">[19]</xref>
. Deep sequencing was used to obtain an accurate picture of viral population diversity before and after passage in cell culture.</p>
</sec>
<sec id="s2"><title>Results</title>
<sec id="s2a"><title>Virus Passage Experiments</title>
<p>Nasal samples were obtained from calves and tested for BCoV RNA via RT-PCR. Naturally-infected BCoV samples are often difficult to adapt to growth in cell culture <xref rid="pone.0052752-Saif1" ref-type="bibr">[20]</xref>
, however, it is essential to use naturally-infected samples rather than laboratory adapted virus strains because isolation in cell culture is likely to result in adaptive mutations not present in the naturally circulating virus. Therefore, every step of the virus passage protocol was optimized to allow for the greatest chance of virus particle recovery before and during cell culture passage including incubation time of the virus on cells, RNA extraction protocol, and RT-PCR protocol. Although trypsin has been shown to enhance infection of HRT-18 cells <xref rid="pone.0052752-Saif1" ref-type="bibr">[20]</xref>
, <xref rid="pone.0052752-StCyrCoats1" ref-type="bibr">[21]</xref>
, it was not used in the infection protocol in order to keep protocols identical for infection of enteric, lung and macrophage cell lines. The three nasal samples (#’s 1, 27, and 59) that had the highest titer as determined by TaqMan RT-PCR were tested using the optimized protocols. A laboratory adapted strain, Nebraska (NEB) was included as a positive control for cell line infectivity.</p>
<p>Infection by the natural BCoV nasal samples was productive in only three of the five cell lines: THP-1, Bomac and HRT-18, as determined by RT-PCR, compared to all five cell lines productively infected by the lab adapted strain, NEB, as determined by RT-PCR. Viral growth was sustained through five passages in the THP-1 and HRT-18 cells and for four passages in the Bomacs (no viral RNA was detected by Taqman assay from passage 5 in Bomacs for any of the viruses except for the NEB strain). For each sample, viral RNA was extracted from samples prior to any passage (“unpassaged virus”) and from the first and the last passage experiments – passage 4 for the Bomacs and passage 5 for the THP-1 and HRT-18 cells.</p>
</sec>
<sec id="s2b"><title>Viral Genome Amplification</title>
<p>Approximately 12 kilobase (kb) of the BCoV 30 kb genome was amplified using 16 primer sets designed to conserved regions of the genome (<xref ref-type="supplementary-material" rid="pone.0052752.s002">Table S1</xref>
). The primer sets were used to amplify unpassaged and passaged viral RNA, however, in some cases not all primer sets yielded PCR product suggesting that changes occurred in the primer binding region during passage.</p>
</sec>
<sec id="s2c"><title>Analysis of Consensus Sequence Data</title>
<p>The nsp1, nsp3, nsp14, and spike protein gene consensus sequences of passaged and unpassaged viruses were compared to identify any changes characteristic of serial passage. Phylogenetic analysis of the consensus sequence showed that two of the naturally-infected nasal samples (#’s 27 and 59) changed much more gradually compared to that of the third such sample (#1) which rapidly changed to resemble that of the laboratory adapted NEB strain after just one passage in cell culture, regardless of cell type (<xref ref-type="fig" rid="pone-0052752-g001">Figure 1</xref>
). The samples formed two major clusters, one which included all the unpassaged samples (“UP” group) and one that included all the NEB reference strain samples and all the passaged samples derived from sample 1 (Passaged group, “P”). Regions of the genome varied in the number of observed mutations in the consensus sequence with the spike and nsp1 genes showing the most mutations due to passage. Interestingly, the consensus sequence of the viruses passaged in Bomac cells changed more rapidly as compared to that of the genomes of viruses passaged in other cell lines (<xref ref-type="fig" rid="pone-0052752-g001">Figure 1</xref>
).</p>
<fig id="pone-0052752-g001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0052752.g001</object-id>
<label>Figure 1</label>
<caption><title>Phylograms generated using consensus sequence for the coding regions of Nsp1, Nsp3, Nsp14, and Spike.</title>
<p>(A–D, respectively.) Samples are labeled using sample sequence identifier (B##), samples #, cell type, and passage number. Phylogeny was inferred by using the Maximum Likelihood method based on the Tamura-Nei model <xref rid="pone.0052752-Tamura1" ref-type="bibr">[42]</xref>
. The tree with the highest log likelihood is shown. The percentage of trees in which the associated taxa clustered together is shown next to the branches. The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. Codon positions included were 1st+2nd+3rd+Noncoding. All positions containing gaps and missing data were eliminated. There were a total of 1671 (nsp1), 851 (nsp3), 1561 (nsp14), and 1457 (spike) positions in the final dataset. The genomic regions included in the analysis span nts 2764–8460 (nsp1), 10858–11714 (nsp3), 17910–19472 (nsp14), and 23641–27732 (spike), as numbered according to accession AB354579. Analyses were conducted using MEGA5 <xref rid="pone.0052752-Tamura2" ref-type="bibr">[43]</xref>
. Samples may be absent from analysis due to inadequate sequence data generated for a particular coding region, however in cases where fewer positions gave the same phylogram as that generated using more positions, fewer positions were used to allow inclusion of more samples in the phylogram.</p>
</caption>
<graphic xlink:href="pone.0052752.g001"></graphic>
</fig>
<p>Most surprisingly, after passage in cell culture, the consensus sequences of many viral samples acquired a 12-nucleotide (nt) insert encoding 4 amino acids (Ser, Arg, Arg, Arg) located at nt 2737 of the spike gene, whereas none of the unpassaged samples contained this insert at the consensus level. This is illustrated in the phylogram of the spike protein (<xref ref-type="fig" rid="pone-0052752-g001">Figure 1D</xref>
) where samples with the insert generally group together in the upper clade of the tree, although other point mutations differed between the groups as well. The four-amino acid insert was located immediately adjacent to the proteolytic cleavage site of the spike protein, and in the S2 subunit of the spike protein which is responsible for fusion entry of the virus into the host cell <xref rid="pone.0052752-Parker1" ref-type="bibr">[22]</xref>
. The addition of multiple basic residues to this region creates a furin cleavage site that may serve as an alternative to trypsin-dependent cleavage of the spike protein <xref rid="pone.0052752-Belouzard1" ref-type="bibr">[23]</xref>
. Data from the 454 sequencing (mean 931×coverage) showed the insert was also present at subconsensus level in some of the cell-passaged samples (<xref ref-type="table" rid="pone-0052752-t001">Table 1</xref>
). An in-depth analysis using Illumina ultra-deep sequencing (mean 37,901×coverage) showed that the insert was present at very low numbers in unpassaged samples but was sometimes enriched during passage in cell culture indicating a possible role in cell culture adaptation.</p>
<table-wrap id="pone-0052752-t001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0052752.t001</object-id>
<label>Table 1</label>
<caption><title>Percentage of 454 reads containing the multibasic insert (Illumina data in parentheses).</title>
</caption>
<alternatives><graphic id="pone-0052752-t001-1" xlink:href="pone.0052752.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">Sample</td>
<td align="left" rowspan="1" colspan="1">27</td>
<td align="left" rowspan="1" colspan="1">59</td>
<td align="left" rowspan="1" colspan="1">1</td>
<td align="left" rowspan="1" colspan="1">NEB</td>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1">UP</td>
<td align="left" rowspan="1" colspan="1">0.1(0.05)</td>
<td align="left" rowspan="1" colspan="1">0(0.01)</td>
<td align="left" rowspan="1" colspan="1">3(0.6)</td>
<td align="left" rowspan="1" colspan="1">na</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">BO.P1</td>
<td align="left" rowspan="1" colspan="1">87</td>
<td align="left" rowspan="1" colspan="1">37</td>
<td align="left" rowspan="1" colspan="1">92</td>
<td align="left" rowspan="1" colspan="1">88</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">BO.P4</td>
<td align="left" rowspan="1" colspan="1">32</td>
<td align="left" rowspan="1" colspan="1">(94)</td>
<td align="left" rowspan="1" colspan="1">94</td>
<td align="left" rowspan="1" colspan="1">95</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">THP.P1</td>
<td align="left" rowspan="1" colspan="1">0.07</td>
<td align="left" rowspan="1" colspan="1">13</td>
<td align="left" rowspan="1" colspan="1">93</td>
<td align="left" rowspan="1" colspan="1">88</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">THP.P5</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">20</td>
<td align="left" rowspan="1" colspan="1">88</td>
<td align="left" rowspan="1" colspan="1">94</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">HRT.P1</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">0.6</td>
<td align="left" rowspan="1" colspan="1">86</td>
<td align="left" rowspan="1" colspan="1">86</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">HRT.P5</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">21</td>
<td align="left" rowspan="1" colspan="1">88</td>
<td align="left" rowspan="1" colspan="1">88</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">EBL.P1</td>
<td align="left" rowspan="1" colspan="1">na</td>
<td align="left" rowspan="1" colspan="1">na</td>
<td align="left" rowspan="1" colspan="1">85</td>
<td align="left" rowspan="1" colspan="1">87</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">A549.P1</td>
<td align="left" rowspan="1" colspan="1">na</td>
<td align="left" rowspan="1" colspan="1">na</td>
<td align="left" rowspan="1" colspan="1">na</td>
<td align="left" rowspan="1" colspan="1">87</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot><fn id="nt101"><label></label>
<p>Columns are for the different samples: 27, 59, 1, and Nebraska strain. Samples are labeled using cell type, and passage number. na: sample data not available.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>To determine if the insert sequence was present in the NEB sample that we obtained from NVSL before it was propagated in our laboratory, a second sample of this strain was requested from NVSL. Unfortunately this strain was no longer available due to a gap in production testing so a sample of the RNA from the NVSL NEB stock was obtained for testing purposes. This RNA sample was tested using a Taqman PCR assay designed to detect the insert and data showed that the insert was present at approximately the consensus level (>50%) in the RNA from NVSL (<xref ref-type="supplementary-material" rid="pone.0052752.s003">Table S2</xref>
).</p>
<p>The insert reached consensus in all passages of nasal sample #1, passage 1 of sample #27 and passage 4 of sample #59 in Bomacs, but in all other sample-passage combinations it was present at subconsensus level between 0.01 and 37% (<xref ref-type="table" rid="pone-0052752-t001">Table 1</xref>
). This 12 nt insert was not present in any BCoV sequences available in GenBank (the NEB strain sequence is not available in GenBank). Because GenBank data typically reflect the consensus sequence present in a sample, it is possible that the insert is present at subconsensus level in other common lab-adapted strains of BCoV. To address this question, BCoV strain Mebus was obtained from BEI Resources and tested for presence of the insert using the Taqman PCR assay. No insert was detected in diluted or undiluted Mebus RNA, even though Mebus polymerase was present at approximately 13,000 copies per ng RNA.</p>
<p>Four additional bovine nasal samples that tested positive for the presence of BCoV polymerase by Taqman assay were tested for the presence of the insert using the Taqman assay. The polymerase gene was detected at very low levels in all of the samples (range 13–367 copies per 5 µl cDNA) and the insert was detected at very low levels in one of these unpassaged samples. Additionally, Illumina and 454 deep sequencing data had been generated from an unpassaged BCoV nasal sample, #74, for a previous study and these data were examined and found to have the insert present at a low level.</p>
<p>To verify the reproducibility of the rapid enrichment of the insert to consensus level in nasal sample #1, we used the Taqman insert assay to test for the presence of the insert in a second sample of #1 that was passed once in Bomac cells and was available from a previous experiment. The Taqman assay detected the insert in the sample in approximately the same numbers as the polymerase gene, thus indicating that the enrichment process is reproducible.</p>
</sec>
<sec id="s2d"><title>Summary of Subconsensus Point Mutations</title>
<p>The number of viral genomic variants (nucleotide polymorphism occurring at subconsensus level) varied widely between samples (<xref ref-type="supplementary-material" rid="pone.0052752.s004">Table S3</xref>
). <xref ref-type="supplementary-material" rid="pone.0052752.s001">Figure S1</xref>
shows that there is little correlation between variant count and coverage, suggesting that the wide differences in variant count could not be explained by variation in coverage by 454 sequencing.</p>
<p>To increase the sensitivity of variant detection, the three unpassaged samples were sequenced a second time using Illumina GA IIx sequencing (sample 59.BO.4 was also resequenced as its 454 data disagreed with Taqman results). Compared to 454 data, Illumina sequencing data detected substantially more variants in nasal sample #27 (81 to 226) and nasal sample #59 (10 to 122) (<xref ref-type="supplementary-material" rid="pone.0052752.s004">Table S3</xref>
). Illumina data for nasal sample #1 showed only a modest increase of variant counts (from 207 to 260) compared to the 454 data. This is consistent with the fact that the majority of variants in nasal sample #1 occurred at higher frequency – high enough that most were detected with 454 sequencing.</p>
<p>To measure the abundance of a variant at a given position in the genome within a sample, we defined the frequency of a variant as the percentage of mapped reads that overlap the query position and contain the target variant divided by the total number of mapped reads that overlap the query position and make a base call of any type. The median variant frequency for samples #27 and #59 was 0.5% but 6% for sample #1. Thus, the increased variant count obtained from Illumina in samples #27 and #59 reflected the presence of a larger pool of variants present at lower frequency levels. On average, the fraction of the sequenced regions of the genome containing a variant was 0.5% for 454 data and 1.5% for the Illumina data (analysis restricted to bases where a consensus nucleotide is called). This is consistent with the higher coverage and sensitivity associated with our Illumina data.</p>
</sec>
<sec id="s2e"><title>Variants in Unpassaged Samples become Consensus in Passaged Samples</title>
<p>We explored the possibility that subconsensus variants support virus adaptation to new cell types by comparing the consensus sequence for each unpassaged parent sample to its passaged descendants. Phylograms generated using the consensus sequences (<xref ref-type="fig" rid="pone-0052752-g001">Figure 1A–D</xref>
for each of 4 proteins studied) highlight those passaged samples that differed substantially from their unpassaged parents (lower clades), as indicated by their clustering to the NEB strain (upper clades). <xref ref-type="fig" rid="pone-0052752-g002">Figure 2</xref>
shows for each sample that clustered away from the unpassaged parent, the percentage of passaged consensus SNPs (single nucleotide polymorphisms) that were present as variants in the unpassaged sample (see Material and Methods for a detailed description of the methods used to rigorously differentiate mutations attributable to sequencing error). The majority of consensus passaged SNPs were present as variants in the respective unpassaged parental samples in all cases except 59.BO1 and 59.BO4, which were impacted by the low variant count in 59 UP (<xref ref-type="supplementary-material" rid="pone.0052752.s004">Table S3</xref>
). Thus results from samples 1 and 27 supports the hypothesis that the mutations present in passage strains were “selected” from a pre-existing pool rather than through <italic>de novo</italic>
mutation and subsequent population fixation.</p>
<fig id="pone-0052752-g002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0052752.g002</object-id>
<label>Figure 2</label>
<caption><title>Percent of consensus SNPs that occur as subconsensus variants in unpassaged samples.</title>
<p>The consensus sequence for each unpassaged parent sample was compared to its passaged descendants. For each sample that clustered away from the unpassaged parent (<xref ref-type="fig" rid="pone-0052752-g001">Figure 1</xref>
), the percentage of passaged consensus SNPs that were present as variants in the unpassaged sample is shown on the y- axis. The x-axis shows the name of the passaged descendant identified by sample (1, 27, or 59), host cell type (BO, THP or HRT) and passage number (1, 4 or 5).</p>
</caption>
<graphic xlink:href="pone.0052752.g002"></graphic>
</fig>
<p>We hypothesized that there are mutation signatures in the genomes of laboratory passaged samples that reflect cell passage adaptation. These mutations are either nonexistent or exist at the subconsensus level in the unpassaged samples, increase in abundance with cell passage, and may eventually reach consensus in highly passaged samples. We identified 186 positions in the BCoV genome where all three unpassaged samples had the same consensus SNP that differed from the consensus SNP found in all eight NEB samples. We defined UP-SNPs and P-SNPs as the consensus SNPs found in the unpassaged parental samples and highly passaged NEB samples respectively, and used the sequence from NEB THP-1 pass 5 (NEB.THP.5) as a highly-passaged representative. Of the 186 P-SNPs, 78%, 24%, and 99% were found to already exist as variants in the unpassaged samples for nasal samples #27, #59 and #1, respectively (using Illumina data), consistent with the results found on a per-sample basis shown in <xref ref-type="fig" rid="pone-0052752-g002">Figure 2</xref>
.</p>
<p>We further made pair-wise comparisons of early vs. late passaged samples derived from the same nasal sample through the same cell line. For example, 27.THP.1 was compared to 27.THP.5, and 59.THP.1 was compared to 59.THP.5. There were 9 such pairs available for this analysis, 3 cell types (Bomac, THP-1, and HRT-18) for each nasal sample (#27, #59 and #1). Of the 186 positions, 42 positions showed an enrichment of the P-SNPs from the early to the late passaged sample in a majority of the cases (at least 5 of the 9 pairs of the passaged samples). Two of these 42 positions are located on the nsp1 gene, the other 40 are all located on the spike gene (<xref ref-type="supplementary-material" rid="pone.0052752.s005">Table S4</xref>
). The point mutations –base change from UP-SNP to P-SNP – at 17 of the 42 locations are non-synonymous (<xref ref-type="table" rid="pone-0052752-t002">Table 2</xref>
).</p>
<table-wrap id="pone-0052752-t002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0052752.t002</object-id>
<label>Table 2</label>
<caption><title>Non-synonymous mutations associated with cell passage.</title>
</caption>
<alternatives><graphic id="pone-0052752-t002-2" xlink:href="pone.0052752.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>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead><tr><td align="left" rowspan="1" colspan="1">Position</td>
<td align="left" rowspan="1" colspan="1">Codon Position</td>
<td align="left" rowspan="1" colspan="1">UP-Codon</td>
<td align="left" rowspan="1" colspan="1">P-Codon</td>
<td align="left" rowspan="1" colspan="1">UP-AA</td>
<td align="left" rowspan="1" colspan="1">P-AA</td>
<td align="left" rowspan="1" colspan="1">Position-AA (spike)</td>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1">6461</td>
<td align="left" rowspan="1" colspan="1">1 (ts)</td>
<td align="left" rowspan="1" colspan="1">TAT</td>
<td align="left" rowspan="1" colspan="1">TGT</td>
<td align="left" rowspan="1" colspan="1">Y</td>
<td align="left" rowspan="1" colspan="1">C</td>
<td align="left" rowspan="1" colspan="1">–</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">23672</td>
<td align="left" rowspan="1" colspan="1">1 (ts)</td>
<td align="left" rowspan="1" colspan="1">ACG</td>
<td align="left" rowspan="1" colspan="1">ATG</td>
<td align="left" rowspan="1" colspan="1">T</td>
<td align="left" rowspan="1" colspan="1">M</td>
<td align="left" rowspan="1" colspan="1">11</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"><bold><italic>23759</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>1 (ts)</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>ACT</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>ATT</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>T</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>I</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>40</italic>
</bold>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"><bold><italic>24083</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>1 (ts)</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>GGT</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>GAT</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>G</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>D</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>148</italic>
</bold>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"><bold><italic>24145</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>0 (tv)</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>AAT</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>CAT</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>N</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>H</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>169</italic>
</bold>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"><bold><italic>24157</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>0 (tv)</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>AAT</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>CAT</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>N</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>H</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>173</italic>
</bold>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"><bold><italic>24160</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>0 (ts)</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>TCT</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>CCT</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>S</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>P</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>174</italic>
</bold>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"><bold><italic>24382</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>0 (tv)</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>ATG</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>CTG</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>M</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>L</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>248</italic>
</bold>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"><bold><italic>24397</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>0 (tv)</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>AAT</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>TCT</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>N</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>S</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>253</italic>
</bold>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"><bold><italic>24398</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>1 (tv)</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>AAT</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>TCT</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>N</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>S</italic>
</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold><italic>253</italic>
</bold>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">25166</td>
<td align="left" rowspan="1" colspan="1">1 (tv)</td>
<td align="left" rowspan="1" colspan="1">ACT</td>
<td align="left" rowspan="1" colspan="1">AAT</td>
<td align="left" rowspan="1" colspan="1">T</td>
<td align="left" rowspan="1" colspan="1">N</td>
<td align="left" rowspan="1" colspan="1">509</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">25267</td>
<td align="left" rowspan="1" colspan="1">0 (tv)</td>
<td align="left" rowspan="1" colspan="1">GCT</td>
<td align="left" rowspan="1" colspan="1">TCT</td>
<td align="left" rowspan="1" colspan="1">A</td>
<td align="left" rowspan="1" colspan="1">S</td>
<td align="left" rowspan="1" colspan="1">543</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">25351</td>
<td align="left" rowspan="1" colspan="1">0 (ts)</td>
<td align="left" rowspan="1" colspan="1">CAT</td>
<td align="left" rowspan="1" colspan="1">TAT</td>
<td align="left" rowspan="1" colspan="1">H</td>
<td align="left" rowspan="1" colspan="1">Y</td>
<td align="left" rowspan="1" colspan="1">571</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">25458</td>
<td align="left" rowspan="1" colspan="1">2 (tv)</td>
<td align="left" rowspan="1" colspan="1">TTG</td>
<td align="left" rowspan="1" colspan="1">TTT</td>
<td align="left" rowspan="1" colspan="1">L</td>
<td align="left" rowspan="1" colspan="1">F</td>
<td align="left" rowspan="1" colspan="1">606</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">25559</td>
<td align="left" rowspan="1" colspan="1">1 (ts)</td>
<td align="left" rowspan="1" colspan="1">ACT</td>
<td align="left" rowspan="1" colspan="1">ATA</td>
<td align="left" rowspan="1" colspan="1">T</td>
<td align="left" rowspan="1" colspan="1">I</td>
<td align="left" rowspan="1" colspan="1">640</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">25945</td>
<td align="left" rowspan="1" colspan="1">0 (ts)</td>
<td align="left" rowspan="1" colspan="1">TCG</td>
<td align="left" rowspan="1" colspan="1">GCG</td>
<td align="left" rowspan="1" colspan="1">S</td>
<td align="left" rowspan="1" colspan="1">A</td>
<td align="left" rowspan="1" colspan="1">769</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">26357</td>
<td align="left" rowspan="1" colspan="1">1 (ts)</td>
<td align="left" rowspan="1" colspan="1">GAT</td>
<td align="left" rowspan="1" colspan="1">GGT</td>
<td align="left" rowspan="1" colspan="1">D</td>
<td align="left" rowspan="1" colspan="1">G</td>
<td align="left" rowspan="1" colspan="1">906</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot><fn id="nt102"><label></label>
<p>UP = unpassaged, P = passaged, AA = amino acid, (ts) = transition, (tv) = transversion. Mutations located in N-terminal receptor binding domain in the spike protein (see <xref ref-type="fig" rid="pone-0052752-g003">Figure 3</xref>
) are in italics.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>All 42 passage markers were found as variants in unpassaged samples #1 and #27, but only 15 were identified in unpassaged #59. Since all of sample #1′s passage derivatives rapidly changed with passage to resemble the NEB strain, it was not surprising that all of the “acquired” NEB strain SNPs were present in the initial unpassaged sample. The genetic changes that occurred in the #27 and #59 derived samples showed a less clear pattern since only a subset of samples acquired NEB strain SNPs as the consensus. We checked if the lack of NEB variants in unpassaged sample #59 was an anomaly by comparing the variant pools in its passaged descendants to the NEB SNPs. There indeed was a large increase in NEB strain SNPs existing as variants in the passaged descendants of #59 compared to the parent sample, suggesting these SNPs were either independently acquired after a single cell passage or were undetected by sequencing of the original unpassaged sample. It seems likely that increased numbers of serial passages would result in samples #27 and #59 acquiring more NEB strain SNPs.</p>
</sec>
<sec id="s2f"><title>Many of the Passage Associated Mutations can be Mapped to the Surface of the Spike Protein</title>
<p>Of the 42 passage associated mutations, 17 of the mutations were non-synonymous and 16 of these were located in the spike gene (mutation sites are listed in <xref ref-type="supplementary-material" rid="pone.0052752.s004">Table S3</xref>
). All 17 non-synonymous mutations (9 transitions (ts) and 8 transversions (tv)) are shown in <xref ref-type="table" rid="pone-0052752-t002">Table 2</xref>
. A structural model of the N-terminal part of the receptor binding domain of the spike protein from the passaged BCoV was constructed using the AS2TS system 5 <xref rid="pone.0052752-Zemla1" ref-type="bibr">[24]</xref>
based on the homology (54% of sequence identity) to the spike protein from a murine coronavirus (PDB template 3r4d_B) <xref rid="pone.0052752-Peng1" ref-type="bibr">[12]</xref>
. As shown in <xref ref-type="fig" rid="pone-0052752-g003">Figure 3</xref>
, at least 7 out of 8 of the non-synonymous mutations (see positions colored in green in <xref ref-type="table" rid="pone-0052752-t002">Table 2</xref>
) occurred at surface-accessible residues in the receptor binding domain. One of these surface exposed amino-acid mutations, N253S required 2 transversions (AAT – TCT). Only one non-synonymous mutation, M248L (1 transversion), was located in a buried position on the spike protein.</p>
<fig id="pone-0052752-g003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0052752.g003</object-id>
<label>Figure 3</label>
<caption><title>Structural model of the receptor binding domain from the BCoV passaged sample.</title>
<p>Ribbons representation (left) and hydrophobic surface (right). Colors of the hydrophobic surface range from blue for the most polar residues to white to orange red for the most hydrophobic residues. Positions of non-synonymous mutations are colored in green (see <xref ref-type="table" rid="pone-0052752-t002">Table 2</xref>
).</p>
</caption>
<graphic xlink:href="pone.0052752.g003"></graphic>
</fig>
</sec>
<sec id="s2g"><title>Haplotype Reconstruction</title>
<p>Given that the 12 nt insert and the majority of the passage associated point mutations were located on the spike protein, we sought to determine whether these passage markers are part of a ‘passaged genotype’, and whether they were linked with the insert on the same haplotype. The open source software, ShoRAH, was used to carry out genotype reconstruction. ShoRAH estimates genetic diversity of a heterogeneous sample from deep sequencing data and has previously been applied to estimating HIV and hepatitis C virus quasispecies <xref rid="pone.0052752-Zagordi1" ref-type="bibr">[25]</xref>
, <xref rid="pone.0052752-Zagordi2" ref-type="bibr">[26]</xref>
. Given any two positions in the genome, ShoRAH returns genotypes found between the two positions based on all reads overlapping the region between the positions.</p>
<p>The possibility that the passage marker mutations were linked to the 12-nt insert on the same genotype was examined. To check for linkage between the SNPs and the insert (at nt position 26374), reads overlapping the two closest passage markers (positions 26361 and 26394) were divided into the “insert-containing” group and the “insert-lacking” group (separately for each sample). For ten of the samples (27.UP, 27.BO.1, 27.THP.1, 27.BO.4, 59.BO.1, 59.THP.1, 59.HRT.1, 59.THP.5, 59.HRT.5, and 1.UP), most of the insert-lacking reads contained the UP-SNPs ‘C’ at 26361 and 26394, and majority of the insert-containing reads showed P-SNPs ‘T’ at these two positions. However, in the case of all of the NEB passaged samples and all passaged samples for nasal sample #1, a different pattern was observed. Here, although the insert-lacking reads occurred at low levels (∼2.5%), they continued to exclusively show the P-SNP variant. The majority genetic population in nasal sample #1 lacked the insert and had the genotype 26361-C, 26394-C while its passaged descendants retained a minority population of insert-lacking genomes that contained the genotype 26361-T, 26394-T (identical to the insert containing dominant genotype). Thus, while there appeared to be two distinct genotypes circulating, it remained unclear whether the insert-lacking passaged variants of nasal sample #1 and NEB originated from a P-variant that lost the insert or rather originated from an insert lacking UP-variant and acquired the P-variant mutations independently.</p>
</sec>
</sec>
<sec id="s3"><title>Discussion</title>
<p>Deep sequencing of three naturally infected BCoV samples show that each sample differed in the composition of mutant spectra, even within the same herd, and this underscores the importance of including multiple naturally-infected samples in virus evolution studies when possible. Furthermore, identification of a relatively common genetic insert in the BCoV genome that is likely to impact phenotype suggests that this type of in-depth analysis of viral populations may be necessary to fully understand evolutionary mechanisms and interspecies transmission events. Indeed, a similar finding was recently described for foot-and-mouth disease virus (FMDV) where deep Illumina sequencing revealed a genotype present at the subconsensus level that enabled the virus to bind heparin sulfate receptors for adaptation to cell culture <xref rid="pone.0052752-Wright1" ref-type="bibr">[27]</xref>
.</p>
<p>The RT-PCR amplification results suggest that the naturally-infected viruses did not reach sufficiently high titers to allow prolonged serial passage in Bomac cell lines but did replicate relatively efficiently in human macrophage and enteric cell lines. Sequence data indicate that the virus population undergoes rapid change when grown in Bomac cells and this may have made the viral genome difficult to detect with our PCR primers. However, the fact that amplification failed for most regions of the genome simultaneously, even regions that are relatively conserved, indicate that the virus was not replicating efficiently in this environment. Alternatively, BCoV is known to persistently infect cell lines, shedding few viral particles in the supernatant <xref rid="pone.0052752-Hofmann1" ref-type="bibr">[28]</xref>
and because RNA was extracted from the supernatant but not the cells, intracellular infections would not have been detected.</p>
<p>Our data indicate that BCoV is able to infect and replicate in human macrophage after exposure to relatively small numbers of virus, and it is possible that this may lead to a more systemic infection due to movement of infected macrophage to the lymph nodes. More mutations were observed in the consensus sequence of the viral genome after passage in macrophage cells, particularly bovine macrophage. This could be a result of increased selective pressure in this cell environment. For example, the production of nitric oxide by macrophage has been shown to increase viral mutation <xref rid="pone.0052752-Akaike1" ref-type="bibr">[29]</xref>
.</p>
<sec id="s3a"><title>Potential Role of Insert in Host Range, Viral Infectivity, and Host Immune Response</title>
<p>Coronavirus infectivity is mediated by the spike protein, which consists of two subunits, the N-terminal S1 subunit and the C-terminal transmembrane S2 subunit. Although the S1 subunit of coronaviruses mediates receptor binding, the S2 subunit mediates fusion of the viral envelope to the host cell membrane <xref rid="pone.0052752-Watanabe1" ref-type="bibr">[30]</xref>
. Potential effects of the BCoV insert in the S2 subunit region of the spike gene include increased host cell range via trypsin-independent fusion host cell entry due to the creation of a furin cleavage site <xref rid="pone.0052752-Watanabe1" ref-type="bibr">[30]</xref>
, and enhanced binding to heparan sulfate on the host cell surface due to the addition of a multibasic region <xref rid="pone.0052752-deHaan1" ref-type="bibr">[31]</xref>
. Although the location of the insert in the S2 subunit may decrease the influence of the insert in heparin sulfate binding, it is interesting to note the presence of other basic residues present in the passaged genotype that may impact binding to heparin sulfate (<xref ref-type="fig" rid="pone-0052752-g003">Figure 3</xref>
).</p>
<p>Studies of other coronaviruses show that fusogenic activity of the spike protein may allow for trypsin-independent entry of the virus into a variety of cell types. For example, the addition of a multibasic motif in the S2 subunit of the BCoV spike protein may create a trypsin-independent spike protein activation site as was recently proposed for SARS-CoV <xref rid="pone.0052752-Belouzard1" ref-type="bibr">[23]</xref>
. Watanabe et al. <xref rid="pone.0052752-Watanabe1" ref-type="bibr">[30]</xref>
generated a SARS-CoV construct with a furin site at the SARS-CoV S2 position (793-KPTKR-797 to 793-KRRKR-797) and have shown SARS-CoV S activation at the cell surface in a trypsin-independent manner. This work highlights the importance of residue R797 in the context of SARS-CoV infection. The position of SARS-CoV residue R797 corresponds to residue 918 in passaged BCoV, and multiple sequence alignments show that this arginine residue is extremely well-conserved across coronaviruses (<xref ref-type="fig" rid="pone-0052752-g004">Figure 4</xref>
). Interestingly, the two arginines (R794 and R795) from the SARS-CoV construct overlap with arginines from the insert 912-SRRR-915 identified in passaged BCoV.</p>
<fig id="pone-0052752-g004" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0052752.g004</object-id>
<label>Figure 4</label>
<caption><title>Multiple sequence alignment of BCoV insert region compared with that of other coronaviruses.</title>
<p>Sequence region (894–927) including insert SRRR from the passaged BCoV (B2.27.BO.P1) was aligned with other coronavirus spike protein amino acid sequences. The inserted amino acids are underlined; amino acids conserved among coronavirus species are in blue text; RR motifs observed in other sequences are in red. In green is highlighted 793-KPTKR-797 region in SARS-CoV S2 domain described as a furin cleavage site allowed trypsin-independent cell-cell fusion <xref rid="pone.0052752-Belouzard1" ref-type="bibr">[<bold>23</bold>
]</xref>
.</p>
</caption>
<graphic xlink:href="pone.0052752.g004"></graphic>
</fig>
<p>Although the multibasic insert may enhance binding of BCoV to a diversity of cell types via binding of heparan sulfate on the host cell surface, heparan sulfate binding has been associated with viral attenuation in vivo due to increased clearance rate <xref rid="pone.0052752-Bernard1" ref-type="bibr">[32]</xref>
–<xref rid="pone.0052752-Spillmann1" ref-type="bibr">[34]</xref>
and this may provide an explanation for the rarity of the insert genotype in naturally-infected samples. Additionally, coronavirus fixation to heparan sulfate on the surface of macrophages has been associated with increased host inflammatory cytokine response via binding of Toll-like receptor 2 <xref rid="pone.0052752-Jacques1" ref-type="bibr">[35]</xref>
. Other environmental factors such as the presence of trypsin-like enzymes in the bovine airway could also influence selective pressures on the frequency of genotypes lacking the insert <xref rid="pone.0052752-Bertram1" ref-type="bibr">[36]</xref>
.</p>
</sec>
<sec id="s3b"><title>Comparison of Amino Acid Substitutions to Previous Study Results</title>
<p>Chouljenko et al. (1998) identified amino acid changes in the spike protein gene that differentiate respiratory BCoV (RBCV) strains from enteric BCoV (EBCV) strains, as well as amino acid changes that were “virulence specific”, meaning that they were not detected in extensively lab-passaged strains that had an attenuated phenotype in bovine infections <xref rid="pone.0052752-Chouljenko2" ref-type="bibr">[37]</xref>
. Our data show 10 of the 13 RBCV-specific changes present in the UP group samples but not the P group samples (<xref ref-type="supplementary-material" rid="pone.0052752.s006">Table S5</xref>
). UP group amino acid changes also correlated with all 7 amino acid changes considered “virulence markers” by Chouljenko et al. (1998). Similarly, 8 of the 10 changes identified in a study by Zhang et al., (1991) as being associated with virulence were present in the UP group but not in the P group <xref rid="pone.0052752-Zhang2" ref-type="bibr">[38]</xref>
. Six of these changes are nonconservative and all of these were reproduced in our data. These results indicate that laboratory passage causes predictable changes in the amino acid sequence, although the rate of change observed in the consensus sequence differs according to cell type and individual sample.</p>
<p>A number of the amino acid changes in the spike protein have been hypothesized to be involved in virulence or tissue tropism of the virus <xref rid="pone.0052752-Chouljenko1" ref-type="bibr">[14]</xref>
. In particular, amino acid change from Ser in the UP group to Ala in the P group was observed at amino acid 769 which is directly adjacent to the spike protein proteolytic cleavage site. A serine at this site has been hypothesized to be associated with cell fusion <xref rid="pone.0052752-Chouljenko2" ref-type="bibr">[37]</xref>
, and our data indicate loss of the serine is associated with gain of the positively-charged insert in the S2 subunit of the spike protein. Another region immediately upstream of the insert (899–907) has been suggested to be a fusogenic domain by Parker et al. <xref rid="pone.0052752-Parker1" ref-type="bibr">[22]</xref>
. Interestingly there is a nonconservative amino acid change in this sequence from a nonpolar amino acid, Ala in Parker et al. (1990) and Gly in Chouljenko et al. (1998) and P group, to an acidic amino acid, Asp, in our UP group. This region is composed almost entirely of nonpolar amino acids and the addition of a charged residue may influence function. Another mutation occurred at amino acid 1030 which is part of the heptad repeat, a domain associated with class I fusion proteins and host range <xref rid="pone.0052752-McRoy1" ref-type="bibr">[39]</xref>
, in this case a Gly in the UP group was replaced with an Asp in the P group. Ultimately, the role of these mutations in viral phenotype will need to be confirmed using reverse genetics.</p>
<p>Using high throughput sequencing, we tracked the evolution of variant genotypes over several passages through cell lines. The identification of variants through deep sequencing of multiple BCoV samples points to the RNA virus’ relatively high mutation rate. However, the majority of the identified variants were detected in relatively low abundance, leaving open the question of their functional significance. Comparing consensus point mutations with their variant counterparts was used as a way to find variants with greater potential for functional significance. Our results showed that rather than <italic>de novo</italic>
mutation being the driver of genetic change at the consensus population level, selection of pre-existing mutations was evident in the majority of cases. Presumably with a sufficiently large population size and extended evolution, novel mutation becomes a bigger factor. The results here illustrate that introduction of novel mutations is a non-trivial process. Moreover, the results indicate the importance of identifying the variants in the population as a potential key piece of information in assessing how a viral isolate would most likely adapt genetically to a new host environment.</p>
</sec>
</sec>
<sec sec-type="materials|methods" id="s4"><title>Materials and Methods</title>
<sec id="s4a"><title>Viruses and Cell Lines</title>
<p>Nebraska (NEB) strain of BCoV was obtained from NVSL (National Veterinary Services Laboratories, Ames, IA) and used as a reference (positive control) for cell culture infections and molecular assays. The stock of NEB strain received from NVSL had been passed >50 times in bovine turbinate and embryonic kidney cells. BCoV strain Mebus was obtained from BEI Resources (Manassas, VA). HRT-18 (human colorectal carcinoma cells) cells (ATCC CRL-11663) were grown on Advanced Dulbecco’s Minimum Essential Medium (D-MEM) (Gibco) with 2 mM L-Glutamine (Gibco), 100 IU/ml Penicillin/100 µg/ml Streptomycin (Gibco) and supplemented with 5% heat-inactivated fetal bovine serum (FBS) (Gibco). EBTr (bovine embryonic tracheal) cells (ATCC CCL-44) and Calu-3 (human lung adenocarcinoma) cells (ATCC HTB-55) were grown on Eagle’s Minimum Essential Medium (EMEM) (Gibco) with 2 mM L-Glutamine, Penicillin/Streptomycin, and supplemented with 10% heat-inactivated FBS. THP-1 (human acute monocytic leukemia) cells (ATCC TIB-202) and Bomac (bovine macrophage <xref rid="pone.0052752-Stabel1" ref-type="bibr">[40]</xref>
; a gift from Dr. Mark Estes, Univ. Texas Medical Branch), were grown on RPMI-1640 Medium (Gibco) with 2 mM L-Glutamine, Penicillin/Streptomycin, and supplemented with 10% heat-inactivated FBS. A549 (human alveolar basal epithelial) cells (obtained from NVSL) and EBL (embryonic bovine lung) cells (ATCC, CCL-185) were grown on Minimum Essential Medium (MEM) (Gibco) with 2 mM L-Glutamine, Penicillin/Streptomycin, and supplemented with 10% heat-inactivated FBS. Above cells were propagated with each optimum culture medium until confluent and maintained in 5% CO2 incubator at 37°C. Cells were detached for passaging using 0.25% Trypsin-EDTA (Gibco), except for THP-1 cells. The THP-1 cells were maintained by either addition of freshly made growth medium or by centrifugation and resuspension of the cells into fresh medium.</p>
</sec>
<sec id="s4b"><title>Collection and Identification of Naturally-infected BCoV Samples</title>
<p>The frequency of BCoV infection is not known for healthy cattle therefore a large number of calves were sampled to ensure that at least several positive samples would be obtained. This study was approved by the Institutional Animal Care and Use Committee at Lawrence Livermore National Laboratory (Protocol Number 2009-207). Approximately 200 nasal cultures swabs were collected from calves and analyzed. Sterile polyester swabs were used to collect nasal samples, placed in 2–3 mL of Eagle’s Minimum Essential Medium (Gibco) and transported on ice back to the laboratory. Collected nasal swabs were homogenized in Eagle’s Minimum Essential Medium supplemented with 1% antibiotic-antimycotic solution (Gibco). The sample suspensions were clarified by centrifugation at 2000×g for 30 minutes, filtered through a 0.22 µm filter and aliquots of about 500–1000 µl were stored at −80°C.</p>
<p>The first group of cattle tested was housed individually in outdoor pens on a dairy farm in Davis, CA. It is likely these calves were relatively stressed due to a several day period of very cold and rainy winter weather. Nasal swabs were obtained from 109 animals and 42 of these animals were positive for BCoV infection as determined by RT-PCR assay <xref rid="pone.0052752-Cho1" ref-type="bibr">[41]</xref>
. Another group of 67 calves were tested during a period of relatively mild winter weather. These calves were part of a herd that included cow-calf pairs that were grazing on pasture. Nasal swabs were obtained from the calves during branding and vaccination procedures. None of these samples were positive for BCoV RNA by RT-PCR assay.</p>
</sec>
<sec id="s4c"><title><italic>In vitro</italic>
Passage of Virus</title>
<p>Confluent cell lines in 12-well plates were infected with 200 µl nasal samples diluted in MEM. Virus samples were incubated for 6 hours, inoculums were removed and cells were rinsed twice in sterile PBS. In the case of THP-1 cells, cells were pelleted in between each step. The appropriate media was added to each cell line. Supernatant from infected cells were transferred to fresh cell monolayers every 72 h. NEB strain stock was used as a positive control for infection of each cell line.</p>
</sec>
<sec id="s4d"><title>RNA Extraction</title>
<p>A 500 µl sample was homogenized with 500 µl TRIzol reagent (Invitrogen) following manufacturer’s protocol except that incubations were performed at 30°C instead of room temperature.</p>
</sec>
<sec id="s4e"><title>Primer Design</title>
<p>Primers were designed to amplify eight regions of 1.5 to 2.5 kb and amplified regions overlapped to allow the genotypes to be reconstructed. Primers were designed to be as sensitive to target strain variants as possible, while still being specific enough to not cross-react with non-targets. In all, 24 primer sets (3 sets for each genome region) were tested to select the two best performing primer sets for amplification of each region (<xref ref-type="supplementary-material" rid="pone.0052752.s002">Table S1</xref>
).</p>
</sec>
<sec id="s4f"><title>Reverse Transcriptase (RT) Polymerase Chain Reaction (PCR)</title>
<p>RNA was converted to cDNA using random hexamers and the SuperScript III First-Strand Synthesis System for RT-PCR (Invitrogen) according to manufacturer’s instructions. In cases where viral titer appeared to be low, 10 µl RNA was added to the RT mixture. The PCR primers and conditions used to detect BCoV RNA present in nasal samples are described in Cho et al., 2001 <xref rid="pone.0052752-Cho1" ref-type="bibr">[41]</xref>
. High fidelity Phusion polymerase (New England BioLabs) was used to amplify regions of the BCoV genome for sequencing following manufacturer’s instructions, using 50 µl reactions with 5 µl of cDNA template. The primer sets used for genome amplification are listed in <xref ref-type="supplementary-material" rid="pone.0052752.s002">Table S1</xref>
. PCR conditions consisted of 94°C for 2 min, followed by 40 cycles of 98°C for 10 s, 60°C for 20 s, and 72°C for 1 min 20 s.</p>
</sec>
<sec id="s4g"><title>Sequencing of a 12 kb Region of the Genome</title>
<p>Each sample was amplified with each primer set and products were pooled, purified using the QIAquick PCR Purification kit (Qiagen), and quantified. Amplicons were submitted for 454 sequencing at DNA Sequencing Center at Brigham Young University and a subset of amplicons were sequenced on an Illumina GA IIx sequencer at Eureka Genomics Hercules, CA. The raw sequencer reads generated from the experiments and the reads aligned to the reference sequence is available at <ext-link ext-link-type="ftp" xlink:href="ftp://gdo-bioinformatics.ucllnl.org/pub/bcv1">ftp://gdo-bioinformatics.ucllnl.org/pub/bcv1</ext-link>
.</p>
</sec>
<sec id="s4h"><title>Taqman PCR</title>
<p>Taqman assays (qRT-PCR) were designed to detect a region of the virus polymerase gene as well as a region of the spike gene that was found to be the site of 12 nucleotide (nt) insert for some samples. RNA was extracted and reverse transcribed as described above. Taqman assay was performed on a 5 µl aliquot of cDNA in 50 µl reaction mix using the 7900 HT Sequence Detection System (Applied Biosystems) and Taqman Universal PCR Master Mix (Applied Biosystems) for amplification of the viral cDNA. Plasmids cloned with inserts from amplification of a larger region of the polymerase gene or spike insert region were used to generate the standard curve for quantitation. Taqman primers and probe for amplification of the polymerase gene were (5′ to 3′): Fwd- <named-content content-type="gene">CCATGTGTCATGCATTGGATT</named-content>
, Rvs- <named-content content-type="gene">CACCGATCATCCTGACAATCA</named-content>
, Probe- <named-content content-type="gene">CCGTGTTAGGATGGTATGGCATACTCCAGTG</named-content>
. Taqman primers and probe for amplification of the spike insert region were (5′ to 3′): Fwd- <named-content content-type="gene">GGTTGTTTAGGAAGCGGTTG</named-content>
, Rvs- <named-content content-type="gene">AGCCTCAACGAAACCGACAT</named-content>
, Probe- <named-content content-type="gene">CCGGCGTAGAAGTAGATCTGCTATAGAGGATTT</named-content>
. Cycling parameters were 50°C for 2 min, 95°C for 10 min, 45 cycles of 95°C for 15 sec and 60°C for 1 min.</p>
</sec>
<sec id="s4i"><title>Control for Analysis of PCR and Sequencing Error Rates</title>
<p>To determine the error rate of the PCR amplification and sequencing processes, a sequencing control plasmid was constructed by inserting a 1 kb fragment of the BCoV polymerase gene into a sequencing plasmid using an Invitrogen TOPO® Cloning Kit for Sequencing. The plasmid region was then Sanger sequenced to determine the sequence of the insert prior to PCR amplification and 454 or Illumina sequencing. The plasmid insert was amplified using the same protocol as the viral samples and submitted for 454 and Illumina sequencing along with the viral PCR products.</p>
</sec>
<sec id="s4j"><title>Viral Sequence Assembly</title>
<p>The sequencing control was used to model error by measuring observed miscall rates of non-consensus bases. Sample consensus sequences were iteratively constructed by first mapping to a closely related reference sequence (Genbank Identifier 15081544), constructing a new consensus sequence, mapping the reads to the new consensus sequence as the reference (in order to increase the number of mapped reads and reduce the potential for alignment error), and repeating the process until the consensus sequence converged to a single value.</p>
<p>The following error filtering rules were used to preclude a read that has been aligned to the consensus sequence from contributing to a base call at a given genome position due to likelihood of error: an indel (in the aligned read) cannot be present within five bases of the query base, the position within the read must be at least five bases away from the ends of the read (to further avoid misalignment); any variant call must be observed in reads that were sequenced from both strands. To further avoid miscalls from misalignment, the adjacent bases in the read (+/−1) must agree with the consensus sequence thus filtering out adjacent variant mutations. These rules eliminated a large percentage of the erroneous variant calls in the control sequences. Thus, rather than relying on a uniform observed error rate covering all positions in the control, we used the average per position, per base call error rate restricted to the cases where variant base calls were identified. The reasoning was that an alternative approach would lead to an artificial lowering of the error rate that includes cases where no variant base calls were observed. This led to an observed per genome position, per base call conditional error rate of 0.001, which reflects the average percentage of reads reporting an erroneous sub-consensus base call at each position in the genome.</p>
<p>The error rate was conditioned on the case when a non-consensus base call is made, rather the error rate averaged over all base calls including those where no-subconsensus mutation is attempted due to the read mapping and quality control exclusion rules and subsequent lack of observed variation. Since this is at least an order of magnitude greater than the mutation frequency expected of coronaviruses it is necessary to observe the sub-consensus mutation in multiple sequencer reads to ensure that a mutation occurs with a high enough frequency in the sample to be unlikely attributable to sequencing error.</p>
<p>To determine the minimum number of reads needed to support the presence of a sub-consensus mutation, the Binomial significance test described by Eriksson et al. (2008) was applied. This approach assumes a binomial error distribution so that for a given sequencer coverage amount, the number of reads containing sub-consensus mutations that would be expected from the assumed sequencing error rate can be calculated. When a mutation occurs in a number of reads that exceeds this amount for a P-value of 0.01 (and applying a Bonferonni correction), the mutation is inferred to be unlikely attributable to sequencing error. The detection sensitivity is therefore coverage and sequencer error dependent. For coverage 1,000× the mutation would be required to occur in 5 reads, which constitute 0.5 percent of the total with at least two reads occurring on opposite strands.</p>
<p>Illumina sequencing was carried out on four samples (<xref ref-type="supplementary-material" rid="pone.0052752.s004">Table S3</xref>
) using Illumina paired-end technology where read pairs of length 112 bp were created with an overlapping region of approximately 80 bp. Although the paired-end technology was originally intended to create read pairs with gaps of known size when the templates are cut to lengths much longer than the read length (e.g. 500 bp templates and 100 bp reads). The known gap size facilitates post-sequencing alignment that is particularly useful for more complex non-viral genomes. When the templates used are closer to the read length the paired-end technology creates read pairs with significant regions of overlap.</p>
<p>These overlapping read pairs essentially provide another mechanism of error checking, as each read pair came from the same template and therefore should be perfectly complementary. Any disagreement between the read pairs would be due to sequencing error. This disagreement between overlapping read pairs was used to correlate mismatch rates and the sequencer generated quality scores. A cutoff of Q30 was selected as the minimum base call quality score, which limits the mismatch error rate, while maintaining high amounts of sequencer coverage. Only base calls supported by both read pairs with the minimum quality score were used to infer an affective erorr rate in the control sequence, which was found to be 0.0005. This rate is the maximal percentage of reads that were observed to disagree at any position in the control. As with the 454 data, this error rate was used as input to the Binomial significance test to determine the minimum number of reads needed to separate sequencing error from true genetic variation.</p>
<p>Open source software, ShoRAH, was used to carry out genotype reconstruction <xref rid="pone.0052752-Zagordi1" ref-type="bibr">[25]</xref>
, <xref rid="pone.0052752-Zagordi2" ref-type="bibr">[26]</xref>
.</p>
<p>454 Pyrosequencing data were obtained for 30 BCoV samples: 7, 7, 8 and 8 samples of nasal samples #27, #59, #1, and the NEB strain, respectively. Across all 30 samples, 12,159 unique bases were sequenced, covering nsp1, nsp2, nsp3 and spike proteins between positions 2647 and 27809 on a reference genome, GenBank accession number NC_003045. Average coverage at the four protein regions are 531×, 2535×, 588× and 623×. The standard deviations of coverage at these four protein regions are 210×, 587×, 108× and 270×. Hence, coverage at the nsp2 genome was both highest and most variable among the four proteins. Four of the 30 samples were also sequenced using Illumina to obtain deeper coverage and clarify differences with the 454 samples. <xref ref-type="supplementary-material" rid="pone.0052752.s004">Table S3</xref>
summarizes the counts of variants detected and the average per base read coverage in each sample, and includes the number of genome positions with 50×coverage or greater, which was used as the minimum coverage cutoff for consensus base calls. Primer regions were excluded to avoid false positive variant calls from non-specific primer binding. An average of 10.7 kb consensus sequence was recovered.</p>
</sec>
</sec>
<sec sec-type="supplementary-material" id="s5"><title>Supporting Information</title>
<supplementary-material content-type="local-data" id="pone.0052752.s001"><label>Figure S1</label>
<caption><p><bold>Average coverage per sample (y-axis) versus count of rare variants detected in each sample, with a linear fit.</bold>
</p>
<p>(DOCX)</p>
</caption>
<media xlink:href="pone.0052752.s001.docx"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0052752.s002"><label>Table S1</label>
<caption><p><bold>Primers used to amplify regions of the BCoV genome.</bold>
</p>
<p>(DOCX)</p>
</caption>
<media xlink:href="pone.0052752.s002.docx"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0052752.s003"><label>Table S2</label>
<caption><p><bold>Amplification of BCoV RNA from NVSL by Taqman using Polymerase and Insert Primers.</bold>
</p>
<p>(DOCX)</p>
</caption>
<media xlink:href="pone.0052752.s003.docx"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0052752.s004"><label>Table S3</label>
<caption><p><bold>Summary deep sequencing data.</bold>
</p>
<p>(DOCX)</p>
</caption>
<media xlink:href="pone.0052752.s004.docx"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0052752.s005"><label>Table S4</label>
<caption><p><bold>Putative genetic markers of cell culture adaptation.</bold>
</p>
<p>(DOCX)</p>
</caption>
<media xlink:href="pone.0052752.s005.docx"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0052752.s006"><label>Table S5</label>
<caption><p><bold>Comparison of amino acid changes induced during passage.</bold>
Laboratory-adapted strains (Mebus and L9) and low passage strains or unpassaged strains</p>
<p>(DOCX)</p>
</caption>
<media xlink:href="pone.0052752.s006.docx"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
</sec>
</body>
<back><ack><p>This work was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344.</p>
</ack>
<ref-list><title>References</title>
<ref id="pone.0052752-Woolhouse1"><label>1</label>
<mixed-citation publication-type="journal"><name><surname>Woolhouse</surname>
<given-names>MEJ</given-names>
</name>
(<year>2006</year>
) <article-title>Where Do Emerging Pathogens Come from?</article-title>
<source>Microbe</source>
<volume>1</volume>
: <fpage>511</fpage>
–<lpage>515</lpage>
.</mixed-citation>
</ref>
<ref id="pone.0052752-Elena1"><label>2</label>
<mixed-citation publication-type="journal"><name><surname>Elena</surname>
<given-names>SF</given-names>
</name>
, <name><surname>Carrasco</surname>
<given-names>P</given-names>
</name>
, <name><surname>Daros</surname>
<given-names>J-A</given-names>
</name>
, <name><surname>Sanjuan</surname>
<given-names>R</given-names>
</name>
(<year>2006</year>
) <article-title>Mechanisms of genetic robustness in RNA viruses</article-title>
. <source>EMBO Rep</source>
<volume>7</volume>
: <fpage>168</fpage>
–<lpage>173</lpage>
.<pub-id pub-id-type="pmid">16452927</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Domingo1"><label>3</label>
<mixed-citation publication-type="journal"><name><surname>Domingo</surname>
<given-names>E</given-names>
</name>
, <name><surname>Sheldon</surname>
<given-names>J</given-names>
</name>
, <name><surname>Perales</surname>
<given-names>C</given-names>
</name>
(<year>2012</year>
) <article-title>Viral Quasispecies Evolution</article-title>
. <source>Microbiology and Molecular Biology Reviews</source>
<volume>76</volume>
: <fpage>159</fpage>
–<lpage>216</lpage>
.<pub-id pub-id-type="pmid">22688811</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Borderia1"><label>4</label>
<mixed-citation publication-type="journal"><name><surname>Borderia</surname>
<given-names>AV</given-names>
</name>
, <name><surname>Stapleford</surname>
<given-names>KA</given-names>
</name>
, <name><surname>Vignuzzi</surname>
<given-names>M</given-names>
</name>
(<year>2011</year>
) <article-title>RNA virus population diversity: implications for inter-species transmission</article-title>
. <source>Current Opinion in Virology</source>
<volume>1</volume>
: <fpage>643</fpage>
–<lpage>648</lpage>
.<pub-id pub-id-type="pmid">22440922</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Graham1"><label>5</label>
<mixed-citation publication-type="journal"><name><surname>Graham</surname>
<given-names>RL</given-names>
</name>
, <name><surname>Baric</surname>
<given-names>RS</given-names>
</name>
(<year>2010</year>
) <article-title>Recombination, Reservoirs, and the Modular Spike: Mechanisms of Coronavirus Cross-Species Transmission</article-title>
. <source>J Virol</source>
<volume>84</volume>
: <fpage>3134</fpage>
–<lpage>3146</lpage>
.<pub-id pub-id-type="pmid">19906932</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Vijgen1"><label>6</label>
<mixed-citation publication-type="journal"><name><surname>Vijgen</surname>
<given-names>L</given-names>
</name>
, <name><surname>Keyaerts</surname>
<given-names>E</given-names>
</name>
, <name><surname>Lemey</surname>
<given-names>P</given-names>
</name>
, <name><surname>Maes</surname>
<given-names>P</given-names>
</name>
, <name><surname>Van Reeth</surname>
<given-names>K</given-names>
</name>
, <etal>et al</etal>
(<year>2006</year>
) <article-title>Evolutionary history of the closely related group 2 coronaviruses: porcine hemagglutinating encephalomyelitis virus, bovine coronavirus, and human coronavirus OC43</article-title>
. <source>J Virol</source>
<volume>80</volume>
: <fpage>7270</fpage>
–<lpage>7274</lpage>
.<pub-id pub-id-type="pmid">16809333</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Vijgen2"><label>7</label>
<mixed-citation publication-type="journal"><name><surname>Vijgen</surname>
<given-names>L</given-names>
</name>
, <name><surname>Keyaerts</surname>
<given-names>E</given-names>
</name>
, <name><surname>Moes</surname>
<given-names>E</given-names>
</name>
, <name><surname>Thoelen</surname>
<given-names>I</given-names>
</name>
, <name><surname>Wollants</surname>
<given-names>E</given-names>
</name>
, <etal>et al</etal>
(<year>2005</year>
) <article-title>Complete genomic sequence of human coronavirus OC43: molecular clock analysis suggests a relatively recent zoonotic coronavirus transmission event</article-title>
. <source>J Virol</source>
<volume>79</volume>
: <fpage>1595</fpage>
–<lpage>1604</lpage>
.<pub-id pub-id-type="pmid">15650185</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Gonzlez1"><label>8</label>
<mixed-citation publication-type="journal"><name><surname>González</surname>
<given-names>JM</given-names>
</name>
, <name><surname>Gomez-Puertas</surname>
<given-names>P</given-names>
</name>
, <name><surname>Cavanagh</surname>
<given-names>D</given-names>
</name>
, <name><surname>Gorbalenya</surname>
<given-names>AE</given-names>
</name>
, <name><surname>Enjuanes</surname>
<given-names>L</given-names>
</name>
(<year>2003</year>
) <article-title>A comparative sequence analysis to revise the current taxonomy of the family <italic>Coronaviridae</italic>
</article-title>
. <source>Archives of Virology</source>
<volume>148</volume>
: <fpage>2207</fpage>
–<lpage>2235</lpage>
.<pub-id pub-id-type="pmid">14579179</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Han1"><label>9</label>
<mixed-citation publication-type="journal"><name><surname>Han</surname>
<given-names>MG</given-names>
</name>
, <name><surname>Cheon</surname>
<given-names>D-S</given-names>
</name>
, <name><surname>Zhang</surname>
<given-names>X</given-names>
</name>
, <name><surname>Saif</surname>
<given-names>LJ</given-names>
</name>
(<year>2006</year>
) <article-title>Cross-Protection against a Human Enteric Coronavirus and a Virulent Bovine Enteric Coronavirus in Gnotobiotic Calves</article-title>
. <source>J Virol</source>
<volume>80</volume>
: <fpage>12350</fpage>
–<lpage>12356</lpage>
.<pub-id pub-id-type="pmid">16971444</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Schultze1"><label>10</label>
<mixed-citation publication-type="journal"><name><surname>Schultze</surname>
<given-names>B</given-names>
</name>
, <name><surname>Herrler</surname>
<given-names>G</given-names>
</name>
(<year>1992</year>
) <article-title>Bovine coronavirus uses N-acetyl-9-O-acetylneuraminic acid as a receptor determinant to initiate the infection of cultured cells</article-title>
. <source>J Gen Virol 73 (Pt</source>
<volume>4)</volume>
: <fpage>901</fpage>
–<lpage>906</lpage>
.</mixed-citation>
</ref>
<ref id="pone.0052752-Li1"><label>11</label>
<mixed-citation publication-type="journal"><name><surname>Li</surname>
<given-names>W</given-names>
</name>
, <name><surname>Wong</surname>
<given-names>S-K</given-names>
</name>
, <name><surname>Li</surname>
<given-names>F</given-names>
</name>
, <name><surname>Kuhn</surname>
<given-names>JH</given-names>
</name>
, <name><surname>Huang</surname>
<given-names>IC</given-names>
</name>
, <etal>et al</etal>
(<year>2006</year>
) <article-title>Animal Origins of the Severe Acute Respiratory Syndrome Coronavirus: Insight from ACE2-S-Protein Interactions</article-title>
. <source>Journal of Virology</source>
<volume>80</volume>
: <fpage>4211</fpage>
–<lpage>4219</lpage>
.<pub-id pub-id-type="pmid">16611880</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Peng1"><label>12</label>
<mixed-citation publication-type="journal"><name><surname>Peng</surname>
<given-names>G</given-names>
</name>
, <name><surname>Sun</surname>
<given-names>D</given-names>
</name>
, <name><surname>Rajashankar</surname>
<given-names>KR</given-names>
</name>
, <name><surname>Qian</surname>
<given-names>Z</given-names>
</name>
, <name><surname>Holmes</surname>
<given-names>KV</given-names>
</name>
, <etal>et al</etal>
(<year>2011</year>
) <article-title>Crystal structure of mouse coronavirus receptor-binding domain complexed with its murine receptor</article-title>
. <source>Proceedings of the National Academy of Sciences</source>
<volume>108</volume>
: <fpage>10696</fpage>
–<lpage>10701</lpage>
.</mixed-citation>
</ref>
<ref id="pone.0052752-Zhang1"><label>13</label>
<mixed-citation publication-type="journal"><name><surname>Zhang</surname>
<given-names>X</given-names>
</name>
, <name><surname>Hasoksuz</surname>
<given-names>M</given-names>
</name>
, <name><surname>Spiro</surname>
<given-names>D</given-names>
</name>
, <name><surname>Halpin</surname>
<given-names>R</given-names>
</name>
, <name><surname>Wang</surname>
<given-names>S</given-names>
</name>
, <etal>et al</etal>
(<year>2007</year>
) <article-title>Quasispecies of bovine enteric and respiratory coronaviruses based on complete genome sequences and genetic changes after tissue culture adaptation</article-title>
. <source>Virology</source>
<volume>363</volume>
: <fpage>1</fpage>
–<lpage>10</lpage>
.<pub-id pub-id-type="pmid">17434558</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Chouljenko1"><label>14</label>
<mixed-citation publication-type="journal"><name><surname>Chouljenko</surname>
<given-names>VN</given-names>
</name>
, <name><surname>Lin</surname>
<given-names>XQ</given-names>
</name>
, <name><surname>Storz</surname>
<given-names>J</given-names>
</name>
, <name><surname>Kousoulas</surname>
<given-names>KG</given-names>
</name>
, <name><surname>Gorbalenya</surname>
<given-names>AE</given-names>
</name>
(<year>2001</year>
) <article-title>Comparison of genomic and predicted amino acid sequences of respiratory and enteric bovine coronaviruses isolated from the same animal with fatal shipping pneumonia</article-title>
. <source>J Gen Virol</source>
<volume>82</volume>
: <fpage>2927</fpage>
–<lpage>2933</lpage>
.<pub-id pub-id-type="pmid">11714968</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Perlman1"><label>15</label>
<mixed-citation publication-type="journal"><name><surname>Perlman</surname>
<given-names>S</given-names>
</name>
, <name><surname>Netland</surname>
<given-names>J</given-names>
</name>
(<year>2009</year>
) <article-title>Coronaviruses post-SARS: update on replication and pathogenesis</article-title>
. <source>Nat Rev Microbiol</source>
<volume>7</volume>
: <fpage>439</fpage>
–<lpage>450</lpage>
.<pub-id pub-id-type="pmid">19430490</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Devaraj1"><label>16</label>
<mixed-citation publication-type="journal"><name><surname>Devaraj</surname>
<given-names>SG</given-names>
</name>
, <name><surname>Wang</surname>
<given-names>N</given-names>
</name>
, <name><surname>Chen</surname>
<given-names>Z</given-names>
</name>
, <name><surname>Chen</surname>
<given-names>Z</given-names>
</name>
, <name><surname>Tseng</surname>
<given-names>M</given-names>
</name>
, <etal>et al</etal>
(<year>2007</year>
) <article-title>Regulation of IRF-3-dependent Innate Immunity by the Papain-like Protease Domain of the Severe Acute Respiratory Syndrome Coronavirus</article-title>
. <source>Journal of Biological Chemistry</source>
<volume>282</volume>
: <fpage>32208</fpage>
–<lpage>32221</lpage>
.<pub-id pub-id-type="pmid">17761676</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Sun1"><label>17</label>
<mixed-citation publication-type="journal"><name><surname>Sun</surname>
<given-names>L</given-names>
</name>
, <name><surname>Xing</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Chen</surname>
<given-names>X</given-names>
</name>
, <name><surname>Zheng</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Yang</surname>
<given-names>Y</given-names>
</name>
, <etal>et al</etal>
(<year>2012</year>
) <article-title>Coronavirus Papain-like Proteases Negatively Regulate Antiviral Innate Immune Response through Disruption of STING-Mediated Signaling</article-title>
. <source>PLoS ONE</source>
<volume>7</volume>
: <fpage>e30802</fpage>
.<pub-id pub-id-type="pmid">22312431</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Almeida1"><label>18</label>
<mixed-citation publication-type="journal"><name><surname>Almeida</surname>
<given-names>MS</given-names>
</name>
, <name><surname>Johnson</surname>
<given-names>MA</given-names>
</name>
, <name><surname>Herrmann</surname>
<given-names>T</given-names>
</name>
, <name><surname>Geralt</surname>
<given-names>M</given-names>
</name>
, <name><surname>Wüthrich</surname>
<given-names>K</given-names>
</name>
(<year>2007</year>
) <article-title>Novel β-Barrel Fold in the Nuclear Magnetic Resonance Structure of the Replicase Nonstructural Protein 1 from the Severe Acute Respiratory Syndrome Coronavirus</article-title>
. <source>Journal of Virology</source>
<volume>81</volume>
: <fpage>3151</fpage>
–<lpage>3161</lpage>
.<pub-id pub-id-type="pmid">17202208</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Brockway1"><label>19</label>
<mixed-citation publication-type="journal"><name><surname>Brockway</surname>
<given-names>SM</given-names>
</name>
, <name><surname>Lu</surname>
<given-names>XT</given-names>
</name>
, <name><surname>Peters</surname>
<given-names>TR</given-names>
</name>
, <name><surname>Dermody</surname>
<given-names>TS</given-names>
</name>
, <name><surname>Denison</surname>
<given-names>MR</given-names>
</name>
(<year>2004</year>
) <article-title>Intracellular Localization and Protein Interactions of the Gene 1 Protein p28 during Mouse Hepatitis Virus Replication</article-title>
. <source>Journal of Virology</source>
<volume>78</volume>
: <fpage>11551</fpage>
–<lpage>11562</lpage>
.<pub-id pub-id-type="pmid">15479796</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Saif1"><label>20</label>
<mixed-citation publication-type="journal"><name><surname>Saif</surname>
<given-names>LJ</given-names>
</name>
, <name><surname>Redman</surname>
<given-names>DR</given-names>
</name>
, <name><surname>Brock</surname>
<given-names>KV</given-names>
</name>
, <name><surname>Kohler</surname>
<given-names>EM</given-names>
</name>
, <name><surname>Heckert</surname>
<given-names>RA</given-names>
</name>
(<year>1988</year>
) <article-title>Winter dysentery in adult dairy cattle: detection of coronavirus in the faeces</article-title>
. <source>Vet Rec</source>
<volume>123</volume>
: <fpage>300</fpage>
–<lpage>301</lpage>
.<pub-id pub-id-type="pmid">2848350</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-StCyrCoats1"><label>21</label>
<mixed-citation publication-type="journal"><name><surname>St Cyr-Coats</surname>
<given-names>KS</given-names>
</name>
, <name><surname>Storz</surname>
<given-names>J</given-names>
</name>
, <name><surname>Hussain</surname>
<given-names>KA</given-names>
</name>
, <name><surname>Schnorr</surname>
<given-names>KL</given-names>
</name>
(<year>1988</year>
) <article-title>Structural proteins of bovine coronavirus strain L 9: effects of the host cell and trypsin treatment</article-title>
. <source>Arch Virol</source>
<volume>103</volume>
: <fpage>35</fpage>
–<lpage>45</lpage>
.<pub-id pub-id-type="pmid">2463821</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Parker1"><label>22</label>
<mixed-citation publication-type="journal"><name><surname>Parker</surname>
<given-names>MD</given-names>
</name>
, <name><surname>Yoo</surname>
<given-names>D</given-names>
</name>
, <name><surname>Cox</surname>
<given-names>GJ</given-names>
</name>
, <name><surname>Babiuk</surname>
<given-names>LA</given-names>
</name>
(<year>1990</year>
) <article-title>Primary structure of the S peplomer gene of bovine coronavirus and surface expression in insect cells</article-title>
. <source>Journal of General Virology</source>
<volume>71</volume>
: <fpage>263</fpage>
–<lpage>270</lpage>
.<pub-id pub-id-type="pmid">2155283</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Belouzard1"><label>23</label>
<mixed-citation publication-type="journal"><name><surname>Belouzard</surname>
<given-names>S</given-names>
</name>
, <name><surname>Chu</surname>
<given-names>VC</given-names>
</name>
, <name><surname>Whittaker</surname>
<given-names>GR</given-names>
</name>
(<year>2009</year>
) <article-title>Activation of the SARS coronavirus spike protein via sequential proteolytic cleavage at two distinct sites</article-title>
. <source>Proceedings of the National Academy of Sciences</source>
<volume>106</volume>
: <fpage>5871</fpage>
–<lpage>5876</lpage>
.</mixed-citation>
</ref>
<ref id="pone.0052752-Zemla1"><label>24</label>
<mixed-citation publication-type="journal"><name><surname>Zemla</surname>
<given-names>A</given-names>
</name>
, <name><surname>Zhou</surname>
<given-names>CE</given-names>
</name>
, <name><surname>Slezak</surname>
<given-names>T</given-names>
</name>
, <name><surname>Kuczmarski</surname>
<given-names>T</given-names>
</name>
, <name><surname>Rama</surname>
<given-names>D</given-names>
</name>
, <etal>et al</etal>
(<year>2005</year>
) <article-title>AS2TS system for protein structure modeling and analysis</article-title>
. <source>Nucleic Acids Research</source>
<volume>33</volume>
: <fpage>W111</fpage>
–<lpage>W115</lpage>
.<pub-id pub-id-type="pmid">15980437</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Zagordi1"><label>25</label>
<mixed-citation publication-type="journal"><name><surname>Zagordi</surname>
<given-names>O</given-names>
</name>
, <name><surname>Bhattacharya</surname>
<given-names>A</given-names>
</name>
, <name><surname>Eriksson</surname>
<given-names>N</given-names>
</name>
, <name><surname>Beerenwinkel</surname>
<given-names>N</given-names>
</name>
(<year>2011</year>
) <article-title>ShoRAH: estimating the genetic diversity of a mixed sample from next-generation sequencing data</article-title>
. <source>BMC Bioinformatics</source>
<volume>12</volume>
: <fpage>119</fpage>
.<pub-id pub-id-type="pmid">21521499</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Zagordi2"><label>26</label>
<mixed-citation publication-type="journal"><name><surname>Zagordi</surname>
<given-names>O</given-names>
</name>
, <name><surname>Klein</surname>
<given-names>R</given-names>
</name>
, <name><surname>Daumer</surname>
<given-names>M</given-names>
</name>
, <name><surname>Beerenwinkel</surname>
<given-names>N</given-names>
</name>
(<year>2010</year>
) <article-title>Error correction of next-generation sequencing data and reliable estimation of HIV quasispecies</article-title>
. <source>Nucleic Acids Res</source>
<volume>38</volume>
: <fpage>7400</fpage>
–<lpage>7409</lpage>
.<pub-id pub-id-type="pmid">20671025</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Wright1"><label>27</label>
<mixed-citation publication-type="journal"><name><surname>Wright</surname>
<given-names>CF</given-names>
</name>
, <name><surname>Morelli</surname>
<given-names>MJ</given-names>
</name>
, <name><surname>Thébaud</surname>
<given-names>G</given-names>
</name>
, <name><surname>Knowles</surname>
<given-names>NJ</given-names>
</name>
, <name><surname>Herzyk</surname>
<given-names>P</given-names>
</name>
, <etal>et al</etal>
(<year>2011</year>
) <article-title>Beyond the Consensus: Dissecting Within-Host Viral Population Diversity of Foot-and-Mouth Disease Virus by Using Next-Generation Genome Sequencing</article-title>
. <source>Journal of Virology</source>
<volume>85</volume>
: <fpage>2266</fpage>
–<lpage>2275</lpage>
.<pub-id pub-id-type="pmid">21159860</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Hofmann1"><label>28</label>
<mixed-citation publication-type="journal"><name><surname>Hofmann</surname>
<given-names>MA</given-names>
</name>
, <name><surname>Sethna</surname>
<given-names>PB</given-names>
</name>
, <name><surname>Brian</surname>
<given-names>DA</given-names>
</name>
(<year>1990</year>
) <article-title>Bovine coronavirus mRNA replication continues throughout persistent infection in cell culture</article-title>
. <source>J Virol</source>
<volume>64</volume>
: <fpage>4108</fpage>
–<lpage>4114</lpage>
.<pub-id pub-id-type="pmid">2384915</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Akaike1"><label>29</label>
<mixed-citation publication-type="journal"><name><surname>Akaike</surname>
<given-names>T</given-names>
</name>
, <name><surname>Fujii</surname>
<given-names>S</given-names>
</name>
, <name><surname>Kato</surname>
<given-names>A</given-names>
</name>
, <name><surname>Yoshitake</surname>
<given-names>J</given-names>
</name>
, <name><surname>Miyamoto</surname>
<given-names>Y</given-names>
</name>
, <etal>et al</etal>
(<year>2000</year>
) <article-title>Viral mutation accelerated by nitric oxide production during infection in vivo</article-title>
. <source>Faseb J</source>
<volume>14</volume>
: <fpage>1447</fpage>
–<lpage>1454</lpage>
.<pub-id pub-id-type="pmid">10877838</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Watanabe1"><label>30</label>
<mixed-citation publication-type="journal"><name><surname>Watanabe</surname>
<given-names>R</given-names>
</name>
, <name><surname>Matsuyama</surname>
<given-names>S</given-names>
</name>
, <name><surname>Shirato</surname>
<given-names>K</given-names>
</name>
, <name><surname>Maejima</surname>
<given-names>M</given-names>
</name>
, <name><surname>Fukushi</surname>
<given-names>S</given-names>
</name>
, <etal>et al</etal>
(<year>2008</year>
) <article-title>Entry from the Cell Surface of Severe Acute Respiratory Syndrome Coronavirus with Cleaved S Protein as Revealed by Pseudotype Virus Bearing Cleaved S Protein</article-title>
. <source>J Virol</source>
<volume>82</volume>
: <fpage>11985</fpage>
–<lpage>11991</lpage>
.<pub-id pub-id-type="pmid">18786990</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-deHaan1"><label>31</label>
<mixed-citation publication-type="journal"><name><surname>de Haan</surname>
<given-names>CAM</given-names>
</name>
, <name><surname>Haijema</surname>
<given-names>BJ</given-names>
</name>
, <name><surname>Schellen</surname>
<given-names>P</given-names>
</name>
, <name><surname>Schreur</surname>
<given-names>PW</given-names>
</name>
, <name><surname>te Lintelo</surname>
<given-names>E</given-names>
</name>
, <etal>et al</etal>
(<year>2008</year>
) <article-title>Cleavage of Group 1 Coronavirus Spike Proteins: How Furin Cleavage Is Traded Off against Heparan Sulfate Binding upon Cell Culture Adaptation</article-title>
. <source>J Virol</source>
<volume>82</volume>
: <fpage>6078</fpage>
–<lpage>6083</lpage>
.<pub-id pub-id-type="pmid">18400867</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Bernard1"><label>32</label>
<mixed-citation publication-type="journal"><name><surname>Bernard</surname>
<given-names>KA</given-names>
</name>
, <name><surname>Klimstra</surname>
<given-names>WB</given-names>
</name>
, <name><surname>Johnston</surname>
<given-names>RE</given-names>
</name>
(<year>2000</year>
) <article-title>Mutations in the E2 Glycoprotein of Venezuelan Equine Encephalitis Virus Confer Heparan Sulfate Interaction, Low Morbidity, and Rapid Clearance from Blood of Mice</article-title>
. <source>Virology</source>
<volume>276</volume>
: <fpage>93</fpage>
–<lpage>103</lpage>
.<pub-id pub-id-type="pmid">11021998</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-deHaan2"><label>33</label>
<mixed-citation publication-type="journal"><name><surname>de Haan</surname>
<given-names>CAM</given-names>
</name>
, <name><surname>Li</surname>
<given-names>Z</given-names>
</name>
, <name><surname>te Lintelo</surname>
<given-names>E</given-names>
</name>
, <name><surname>Bosch</surname>
<given-names>BJ</given-names>
</name>
, <name><surname>Haijema</surname>
<given-names>BJ</given-names>
</name>
, <etal>et al</etal>
(<year>2005</year>
) <article-title>Murine Coronavirus with an Extended Host Range Uses Heparan Sulfate as an Entry Receptor</article-title>
. <source>J Virol</source>
<volume>79</volume>
: <fpage>14451</fpage>
–<lpage>14456</lpage>
.<pub-id pub-id-type="pmid">16254381</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Spillmann1"><label>34</label>
<mixed-citation publication-type="journal"><name><surname>Spillmann</surname>
<given-names>D</given-names>
</name>
(<year>2001</year>
) <article-title>Heparan sulfate: Anchor for viral intruders?</article-title>
<source>Biochimie</source>
<volume>83</volume>
: <fpage>811</fpage>
–<lpage>817</lpage>
.<pub-id pub-id-type="pmid">11530214</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Jacques1"><label>35</label>
<mixed-citation publication-type="journal"><name><surname>Jacques</surname>
<given-names>A</given-names>
</name>
, <name><surname>Bleau</surname>
<given-names>C</given-names>
</name>
, <name><surname>Turbide</surname>
<given-names>C</given-names>
</name>
, <name><surname>Beauchemin</surname>
<given-names>N</given-names>
</name>
, <name><surname>Lamontagne</surname>
<given-names>L</given-names>
</name>
(<year>2009</year>
) <article-title>Macrophage interleukin-6 and tumour necrosis factor-alpha are induced by coronavirus fixation to Toll-like receptor 2/heparan sulphate receptors but not carcinoembryonic cell adhesion antigen 1a</article-title>
. <source>Immunology</source>
<volume>128</volume>
: <fpage>e181</fpage>
–<lpage>192</lpage>
.<pub-id pub-id-type="pmid">19740307</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Bertram1"><label>36</label>
<mixed-citation publication-type="journal"><name><surname>Bertram</surname>
<given-names>S</given-names>
</name>
, <name><surname>Glowacka</surname>
<given-names>I</given-names>
</name>
, <name><surname>Müller</surname>
<given-names>MA</given-names>
</name>
, <name><surname>Lavender</surname>
<given-names>H</given-names>
</name>
, <name><surname>Gnirss</surname>
<given-names>K</given-names>
</name>
, <etal>et al</etal>
(<year>2011</year>
) <article-title>Cleavage and Activation of the Severe Acute Respiratory Syndrome Coronavirus Spike Protein by Human Airway Trypsin-Like Protease</article-title>
. <source>Journal of Virology</source>
<volume>85</volume>
: <fpage>13363</fpage>
–<lpage>13372</lpage>
.<pub-id pub-id-type="pmid">21994442</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Chouljenko2"><label>37</label>
<mixed-citation publication-type="journal"><name><surname>Chouljenko</surname>
<given-names>VN</given-names>
</name>
, <name><surname>Kousoulas</surname>
<given-names>KG</given-names>
</name>
, <name><surname>Lin</surname>
<given-names>X</given-names>
</name>
, <name><surname>Storz</surname>
<given-names>J</given-names>
</name>
(<year>1998</year>
) <article-title>Nucleotide and predicted amino acid sequences of all genes encoded by the 3′ genomic portion (9.5 kb) of respiratory bovine coronaviruses and comparisons among respiratory and enteric coronaviruses</article-title>
. <source>Virus Genes</source>
<volume>17</volume>
: <fpage>33</fpage>
–<lpage>42</lpage>
.<pub-id pub-id-type="pmid">9778786</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Zhang2"><label>38</label>
<mixed-citation publication-type="journal"><name><surname>Zhang</surname>
<given-names>XM</given-names>
</name>
, <name><surname>Kousoulas</surname>
<given-names>KG</given-names>
</name>
, <name><surname>Storz</surname>
<given-names>J</given-names>
</name>
(<year>1991</year>
) <article-title>Comparison of the nucleotide and deduced amino acid sequences of the S genes specified by virulent and avirulent strains of bovine coronaviruses</article-title>
. <source>Virology</source>
<volume>183</volume>
: <fpage>397</fpage>
–<lpage>404</lpage>
.<pub-id pub-id-type="pmid">2053289</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-McRoy1"><label>39</label>
<mixed-citation publication-type="journal"><name><surname>McRoy</surname>
<given-names>WC</given-names>
</name>
, <name><surname>Baric</surname>
<given-names>RS</given-names>
</name>
(<year>2008</year>
) <article-title>Amino acid substitutions in the S2 subunit of mouse hepatitis virus variant V51 encode determinants of host range expansion</article-title>
. <source>J Virol</source>
<volume>82</volume>
: <fpage>1414</fpage>
–<lpage>1424</lpage>
.<pub-id pub-id-type="pmid">18032498</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Stabel1"><label>40</label>
<mixed-citation publication-type="journal"><name><surname>Stabel</surname>
<given-names>JR</given-names>
</name>
, <name><surname>Stabel</surname>
<given-names>TJ</given-names>
</name>
(<year>1995</year>
) <article-title>Immortalization and characterization of bovine peritoneal macrophages transfected with SV40 plasmid DNA</article-title>
. <source>Veterinary Immunology and Immunopathology</source>
<volume>45</volume>
: <fpage>211</fpage>
–<lpage>220</lpage>
.<pub-id pub-id-type="pmid">7676607</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Cho1"><label>41</label>
<mixed-citation publication-type="journal"><name><surname>Cho</surname>
<given-names>KO</given-names>
</name>
, <name><surname>Hasoksuz</surname>
<given-names>M</given-names>
</name>
, <name><surname>Nielsen</surname>
<given-names>PR</given-names>
</name>
, <name><surname>Chang</surname>
<given-names>KO</given-names>
</name>
, <name><surname>Lathrop</surname>
<given-names>S</given-names>
</name>
, <etal>et al</etal>
(<year>2001</year>
) <article-title>Cross-protection studies between respiratory and calf diarrhea and winter dysentery coronavirus strains in calves and RT-PCR and nested PCR for their detection</article-title>
. <source>Arch Virol</source>
<volume>146</volume>
: <fpage>2401</fpage>
–<lpage>2419</lpage>
.<pub-id pub-id-type="pmid">11811688</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Tamura1"><label>42</label>
<mixed-citation publication-type="journal"><name><surname>Tamura</surname>
<given-names>K</given-names>
</name>
, <name><surname>Nei</surname>
<given-names>M</given-names>
</name>
(<year>1993</year>
) <article-title>Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees</article-title>
. <source>Molecular Biology and Evolution</source>
<volume>10</volume>
: <fpage>512</fpage>
–<lpage>526</lpage>
.<pub-id pub-id-type="pmid">8336541</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0052752-Tamura2"><label>43</label>
<mixed-citation publication-type="journal"><name><surname>Tamura</surname>
<given-names>K</given-names>
</name>
, <name><surname>Peterson</surname>
<given-names>D</given-names>
</name>
, <name><surname>Peterson</surname>
<given-names>N</given-names>
</name>
, <name><surname>Stecher</surname>
<given-names>G</given-names>
</name>
, <name><surname>Nei</surname>
<given-names>M</given-names>
</name>
, <etal>et al</etal>
(<year>2011</year>
) <article-title>MEGA5: Molecular Evolutionary Genetics Analysis Using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods</article-title>
. <source>Molecular Biology and Evolution</source>
<volume>28</volume>
: <fpage>2731</fpage>
–<lpage>2739</lpage>
.<pub-id pub-id-type="pmid">21546353</pub-id>
</mixed-citation>
</ref>
</ref-list>
</back>
</pmc>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Sante/explor/SrasV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 001382 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 001382 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Sante |area= SrasV1 |flux= Pmc |étape= Corpus |type= RBID |clé= PMC:3538757 |texte= The Role of Viral Population Diversity in Adaptation of Bovine Coronavirus to New Host Environments }}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:23308119" \ | HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \ | NlmPubMed2Wicri -a SrasV1
![]() | This area was generated with Dilib version V0.6.33. | ![]() |