Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Random Codon Re-encoding Induces Stable Reduction of Replicative Fitness of Chikungunya Virus in Primate and Mosquito Cells</title>
<author><name sortKey="Nougairede, Antoine" sort="Nougairede, Antoine" uniqKey="Nougairede A" first="Antoine" last="Nougairede">Antoine Nougairede</name>
<affiliation><nlm:aff id="aff1"><addr-line>Aix Marseille Univ, IRD French Institute of Research for Development, EHESP French School of Public Health, UMR_D 190 “Emergence des Pathologies Virales,” Marseille, France</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="De Fabritus, Lauriane" sort="De Fabritus, Lauriane" uniqKey="De Fabritus L" first="Lauriane" last="De Fabritus">Lauriane De Fabritus</name>
<affiliation><nlm:aff id="aff1"><addr-line>Aix Marseille Univ, IRD French Institute of Research for Development, EHESP French School of Public Health, UMR_D 190 “Emergence des Pathologies Virales,” Marseille, France</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Aubry, Fabien" sort="Aubry, Fabien" uniqKey="Aubry F" first="Fabien" last="Aubry">Fabien Aubry</name>
<affiliation><nlm:aff id="aff1"><addr-line>Aix Marseille Univ, IRD French Institute of Research for Development, EHESP French School of Public Health, UMR_D 190 “Emergence des Pathologies Virales,” Marseille, France</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Gould, Ernest A" sort="Gould, Ernest A" uniqKey="Gould E" first="Ernest A." last="Gould">Ernest A. Gould</name>
<affiliation><nlm:aff id="aff1"><addr-line>Aix Marseille Univ, IRD French Institute of Research for Development, EHESP French School of Public Health, UMR_D 190 “Emergence des Pathologies Virales,” Marseille, France</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Holmes, Edward C" sort="Holmes, Edward C" uniqKey="Holmes E" first="Edward C." last="Holmes">Edward C. Holmes</name>
<affiliation><nlm:aff id="aff2"><addr-line>Sydney Emerging Infections and Biosecurity Institute, School of Biological Sciences and Sydney Medical School, The University of Sydney, Sydney, New South Wales, Australia</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff3"><addr-line>Fogarty International Center, National Institutes of Health, Bethesda, Maryland, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="De Lamballerie, Xavier" sort="De Lamballerie, Xavier" uniqKey="De Lamballerie X" first="Xavier" last="De Lamballerie">Xavier De Lamballerie</name>
<affiliation><nlm:aff id="aff1"><addr-line>Aix Marseille Univ, IRD French Institute of Research for Development, EHESP French School of Public Health, UMR_D 190 “Emergence des Pathologies Virales,” Marseille, France</addr-line>
</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">23436995</idno>
<idno type="pmc">3578757</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3578757</idno>
<idno type="RBID">PMC:3578757</idno>
<idno type="doi">10.1371/journal.ppat.1003172</idno>
<date when="2013">2013</date>
<idno type="wicri:Area/Pmc/Corpus">002C59</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">002C59</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Random Codon Re-encoding Induces Stable Reduction of Replicative Fitness of Chikungunya Virus in Primate and Mosquito Cells</title>
<author><name sortKey="Nougairede, Antoine" sort="Nougairede, Antoine" uniqKey="Nougairede A" first="Antoine" last="Nougairede">Antoine Nougairede</name>
<affiliation><nlm:aff id="aff1"><addr-line>Aix Marseille Univ, IRD French Institute of Research for Development, EHESP French School of Public Health, UMR_D 190 “Emergence des Pathologies Virales,” Marseille, France</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="De Fabritus, Lauriane" sort="De Fabritus, Lauriane" uniqKey="De Fabritus L" first="Lauriane" last="De Fabritus">Lauriane De Fabritus</name>
<affiliation><nlm:aff id="aff1"><addr-line>Aix Marseille Univ, IRD French Institute of Research for Development, EHESP French School of Public Health, UMR_D 190 “Emergence des Pathologies Virales,” Marseille, France</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Aubry, Fabien" sort="Aubry, Fabien" uniqKey="Aubry F" first="Fabien" last="Aubry">Fabien Aubry</name>
<affiliation><nlm:aff id="aff1"><addr-line>Aix Marseille Univ, IRD French Institute of Research for Development, EHESP French School of Public Health, UMR_D 190 “Emergence des Pathologies Virales,” Marseille, France</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Gould, Ernest A" sort="Gould, Ernest A" uniqKey="Gould E" first="Ernest A." last="Gould">Ernest A. Gould</name>
<affiliation><nlm:aff id="aff1"><addr-line>Aix Marseille Univ, IRD French Institute of Research for Development, EHESP French School of Public Health, UMR_D 190 “Emergence des Pathologies Virales,” Marseille, France</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Holmes, Edward C" sort="Holmes, Edward C" uniqKey="Holmes E" first="Edward C." last="Holmes">Edward C. Holmes</name>
<affiliation><nlm:aff id="aff2"><addr-line>Sydney Emerging Infections and Biosecurity Institute, School of Biological Sciences and Sydney Medical School, The University of Sydney, Sydney, New South Wales, Australia</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff3"><addr-line>Fogarty International Center, National Institutes of Health, Bethesda, Maryland, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="De Lamballerie, Xavier" sort="De Lamballerie, Xavier" uniqKey="De Lamballerie X" first="Xavier" last="De Lamballerie">Xavier De Lamballerie</name>
<affiliation><nlm:aff id="aff1"><addr-line>Aix Marseille Univ, IRD French Institute of Research for Development, EHESP French School of Public Health, UMR_D 190 “Emergence des Pathologies Virales,” Marseille, France</addr-line>
</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">PLoS Pathogens</title>
<idno type="ISSN">1553-7366</idno>
<idno type="eISSN">1553-7374</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>Large-scale codon re-encoding represents a powerful method of attenuating viruses to generate safe and cost-effective vaccines. In contrast to specific approaches of codon re-encoding which modify genome-scale properties, we evaluated the effects of random codon re-encoding on the re-emerging human pathogen Chikungunya virus (CHIKV), and assessed the stability of the resultant viruses during serial <italic>in cellulo</italic>
passage. Using different combinations of three 1.4 kb randomly re-encoded regions located throughout the CHIKV genome six codon re-encoded viruses were obtained. Introducing a large number of slightly deleterious synonymous mutations reduced the replicative fitness of CHIKV in both primate and arthropod cells, demonstrating the impact of synonymous mutations on fitness. Decrease of replicative fitness correlated with the extent of re-encoding, an observation that may assist in the modulation of viral attenuation. The wild-type and two re-encoded viruses were passaged 50 times either in primate or insect cells, or in each cell line alternately. These viruses were analyzed using detailed fitness assays, complete genome sequences and the analysis of intra-population genetic diversity. The response to codon re-encoding and adaptation to culture conditions occurred simultaneously, resulting in significant replicative fitness increases for both re-encoded and wild type viruses. Importantly, however, the most re-encoded virus failed to recover its replicative fitness. Evolution of these viruses in response to codon re-encoding was largely characterized by the emergence of both synonymous and non-synonymous mutations, sometimes located in genomic regions other than those involving re-encoding, and multiple convergent and compensatory mutations. However, there was a striking absence of codon reversion (<0.4%). Finally, multiple mutations were rapidly fixed in primate cells, whereas mosquito cells acted as a brake on evolution. In conclusion, random codon re-encoding provides important information on the evolution and genetic stability of CHIKV viruses and could be exploited to develop a safe, live attenuated CHIKV vaccine.</p>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Guzman, Mg" uniqKey="Guzman M">MG Guzman</name>
</author>
<author><name sortKey="Halstead, Sb" uniqKey="Halstead S">SB Halstead</name>
</author>
<author><name sortKey="Artsob, H" uniqKey="Artsob H">H Artsob</name>
</author>
<author><name sortKey="Buchy, P" uniqKey="Buchy P">P Buchy</name>
</author>
<author><name sortKey="Farrar, J" uniqKey="Farrar J">J Farrar</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lanciotti, Rs" uniqKey="Lanciotti R">RS Lanciotti</name>
</author>
<author><name sortKey="Roehrig, Jt" uniqKey="Roehrig J">JT Roehrig</name>
</author>
<author><name sortKey="Deubel, V" uniqKey="Deubel V">V Deubel</name>
</author>
<author><name sortKey="Smith, J" uniqKey="Smith J">J Smith</name>
</author>
<author><name sortKey="Parker, M" uniqKey="Parker M">M Parker</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Staples, Je" uniqKey="Staples J">JE Staples</name>
</author>
<author><name sortKey="Breiman, Rf" uniqKey="Breiman R">RF Breiman</name>
</author>
<author><name sortKey="Powers, Am" uniqKey="Powers A">AM Powers</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Boshra, H" uniqKey="Boshra H">H Boshra</name>
</author>
<author><name sortKey="Lorenzo, G" uniqKey="Lorenzo G">G Lorenzo</name>
</author>
<author><name sortKey="Busquets, N" uniqKey="Busquets N">N Busquets</name>
</author>
<author><name sortKey="Brun, A" uniqKey="Brun A">A Brun</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Randolph, Se" uniqKey="Randolph S">SE Randolph</name>
</author>
<author><name sortKey="Rogers, Dj" uniqKey="Rogers D">DJ Rogers</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Weaver, Sc" uniqKey="Weaver S">SC Weaver</name>
</author>
<author><name sortKey="Reisen, Wk" uniqKey="Reisen W">WK Reisen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hershberg, R" uniqKey="Hershberg R">R Hershberg</name>
</author>
<author><name sortKey="Petrov, Da" uniqKey="Petrov D">DA Petrov</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jenkins, Gm" uniqKey="Jenkins G">GM Jenkins</name>
</author>
<author><name sortKey="Holmes, Ec" uniqKey="Holmes E">EC Holmes</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sanchez, G" uniqKey="Sanchez G">G Sanchez</name>
</author>
<author><name sortKey="Bosch, A" uniqKey="Bosch A">A Bosch</name>
</author>
<author><name sortKey="Pinto, Rm" uniqKey="Pinto R">RM Pinto</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mueller, S" uniqKey="Mueller S">S Mueller</name>
</author>
<author><name sortKey="Coleman, Jr" uniqKey="Coleman J">JR Coleman</name>
</author>
<author><name sortKey="Papamichail, D" uniqKey="Papamichail D">D Papamichail</name>
</author>
<author><name sortKey="Ward, Cb" uniqKey="Ward C">CB Ward</name>
</author>
<author><name sortKey="Nimnual, A" uniqKey="Nimnual A">A Nimnual</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Burns, Cc" uniqKey="Burns C">CC Burns</name>
</author>
<author><name sortKey="Campagnoli, R" uniqKey="Campagnoli R">R Campagnoli</name>
</author>
<author><name sortKey="Shaw, J" uniqKey="Shaw J">J Shaw</name>
</author>
<author><name sortKey="Vincent, A" uniqKey="Vincent A">A Vincent</name>
</author>
<author><name sortKey="Jorba, J" uniqKey="Jorba J">J Jorba</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Burns, Cc" uniqKey="Burns C">CC Burns</name>
</author>
<author><name sortKey="Shaw, J" uniqKey="Shaw J">J Shaw</name>
</author>
<author><name sortKey="Campagnoli, R" uniqKey="Campagnoli R">R Campagnoli</name>
</author>
<author><name sortKey="Jorba, J" uniqKey="Jorba J">J Jorba</name>
</author>
<author><name sortKey="Vincent, A" uniqKey="Vincent A">A Vincent</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Coleman, Jr" uniqKey="Coleman J">JR Coleman</name>
</author>
<author><name sortKey="Papamichail, D" uniqKey="Papamichail D">D Papamichail</name>
</author>
<author><name sortKey="Skiena, S" uniqKey="Skiena S">S Skiena</name>
</author>
<author><name sortKey="Futcher, B" uniqKey="Futcher B">B Futcher</name>
</author>
<author><name sortKey="Wimmer, E" uniqKey="Wimmer E">E Wimmer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mueller, S" uniqKey="Mueller S">S Mueller</name>
</author>
<author><name sortKey="Papamichail, D" uniqKey="Papamichail D">D Papamichail</name>
</author>
<author><name sortKey="Coleman, Jr" uniqKey="Coleman J">JR Coleman</name>
</author>
<author><name sortKey="Skiena, S" uniqKey="Skiena S">S Skiena</name>
</author>
<author><name sortKey="Wimmer, E" uniqKey="Wimmer E">E Wimmer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bull, Jj" uniqKey="Bull J">JJ Bull</name>
</author>
<author><name sortKey="Molineux, Ij" uniqKey="Molineux I">IJ Molineux</name>
</author>
<author><name sortKey="Wilke, Co" uniqKey="Wilke C">CO Wilke</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Coffey, Ll" uniqKey="Coffey L">LL Coffey</name>
</author>
<author><name sortKey="Vasilakis, N" uniqKey="Vasilakis N">N Vasilakis</name>
</author>
<author><name sortKey="Brault, Ac" uniqKey="Brault A">AC Brault</name>
</author>
<author><name sortKey="Powers, Am" uniqKey="Powers A">AM Powers</name>
</author>
<author><name sortKey="Tripet, F" uniqKey="Tripet F">F Tripet</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Strauss, Jh" uniqKey="Strauss J">JH Strauss</name>
</author>
<author><name sortKey="Strauss, Eg" uniqKey="Strauss E">EG Strauss</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Robinson, Mc" uniqKey="Robinson M">MC Robinson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Schuffenecker, I" uniqKey="Schuffenecker I">I Schuffenecker</name>
</author>
<author><name sortKey="Iteman, I" uniqKey="Iteman I">I Iteman</name>
</author>
<author><name sortKey="Michault, A" uniqKey="Michault A">A Michault</name>
</author>
<author><name sortKey="Murri, S" uniqKey="Murri S">S Murri</name>
</author>
<author><name sortKey="Frangeul, L" uniqKey="Frangeul L">L Frangeul</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pialoux, G" uniqKey="Pialoux G">G Pialoux</name>
</author>
<author><name sortKey="Gauzere, Ba" uniqKey="Gauzere B">BA Gauzere</name>
</author>
<author><name sortKey="Jaureguiberry, S" uniqKey="Jaureguiberry S">S Jaureguiberry</name>
</author>
<author><name sortKey="Strobel, M" uniqKey="Strobel M">M Strobel</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Volk, Sm" uniqKey="Volk S">SM Volk</name>
</author>
<author><name sortKey="Chen, R" uniqKey="Chen R">R Chen</name>
</author>
<author><name sortKey="Tsetsarkin, Ka" uniqKey="Tsetsarkin K">KA Tsetsarkin</name>
</author>
<author><name sortKey="Adams, Ap" uniqKey="Adams A">AP Adams</name>
</author>
<author><name sortKey="Garcia, Ti" uniqKey="Garcia T">TI Garcia</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tsetsarkin, Ka" uniqKey="Tsetsarkin K">KA Tsetsarkin</name>
</author>
<author><name sortKey="Vanlandingham, Dl" uniqKey="Vanlandingham D">DL Vanlandingham</name>
</author>
<author><name sortKey="Mcgee, Ce" uniqKey="Mcgee C">CE McGee</name>
</author>
<author><name sortKey="Higgs, S" uniqKey="Higgs S">S Higgs</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="De Lamballerie, X" uniqKey="De Lamballerie X">X de Lamballerie</name>
</author>
<author><name sortKey="Leroy, E" uniqKey="Leroy E">E Leroy</name>
</author>
<author><name sortKey="Charrel, Rn" uniqKey="Charrel R">RN Charrel</name>
</author>
<author><name sortKey="Ttsetsarkin, K" uniqKey="Ttsetsarkin K">K Ttsetsarkin</name>
</author>
<author><name sortKey="Higgs, S" uniqKey="Higgs S">S Higgs</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gerardin, P" uniqKey="Gerardin P">P Gerardin</name>
</author>
<author><name sortKey="Guernier, V" uniqKey="Guernier V">V Guernier</name>
</author>
<author><name sortKey="Perrau, J" uniqKey="Perrau J">J Perrau</name>
</author>
<author><name sortKey="Fianu, A" uniqKey="Fianu A">A Fianu</name>
</author>
<author><name sortKey="Le Roux, K" uniqKey="Le Roux K">K Le Roux</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Economopoulou, A" uniqKey="Economopoulou A">A Economopoulou</name>
</author>
<author><name sortKey="Dominguez, M" uniqKey="Dominguez M">M Dominguez</name>
</author>
<author><name sortKey="Helynck, B" uniqKey="Helynck B">B Helynck</name>
</author>
<author><name sortKey="Sissoko, D" uniqKey="Sissoko D">D Sissoko</name>
</author>
<author><name sortKey="Wichmann, O" uniqKey="Wichmann O">O Wichmann</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gerardin, P" uniqKey="Gerardin P">P Gerardin</name>
</author>
<author><name sortKey="Barau, G" uniqKey="Barau G">G Barau</name>
</author>
<author><name sortKey="Michault, A" uniqKey="Michault A">A Michault</name>
</author>
<author><name sortKey="Bintner, M" uniqKey="Bintner M">M Bintner</name>
</author>
<author><name sortKey="Randrianaivo, H" uniqKey="Randrianaivo H">H Randrianaivo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gould, Ea" uniqKey="Gould E">EA Gould</name>
</author>
<author><name sortKey="Gallian, P" uniqKey="Gallian P">P Gallian</name>
</author>
<author><name sortKey="De Lamballerie, X" uniqKey="De Lamballerie X">X De Lamballerie</name>
</author>
<author><name sortKey="Charrel, Rn" uniqKey="Charrel R">RN Charrel</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Rezza, G" uniqKey="Rezza G">G Rezza</name>
</author>
<author><name sortKey="Nicoletti, L" uniqKey="Nicoletti L">L Nicoletti</name>
</author>
<author><name sortKey="Angelini, R" uniqKey="Angelini R">R Angelini</name>
</author>
<author><name sortKey="Romi, R" uniqKey="Romi R">R Romi</name>
</author>
<author><name sortKey="Finarelli, Ac" uniqKey="Finarelli A">AC Finarelli</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Casals, J" uniqKey="Casals J">J Casals</name>
</author>
<author><name sortKey="Brown, Lv" uniqKey="Brown L">LV Brown</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Clarke, Dh" uniqKey="Clarke D">DH Clarke</name>
</author>
<author><name sortKey="Casals, J" uniqKey="Casals J">J Casals</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Roux, L" uniqKey="Roux L">L Roux</name>
</author>
<author><name sortKey="Simon, Ae" uniqKey="Simon A">AE Simon</name>
</author>
<author><name sortKey="Holland, Jj" uniqKey="Holland J">JJ Holland</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kam, Yw" uniqKey="Kam Y">YW Kam</name>
</author>
<author><name sortKey="Lum, Fm" uniqKey="Lum F">FM Lum</name>
</author>
<author><name sortKey="Teo, Th" uniqKey="Teo T">TH Teo</name>
</author>
<author><name sortKey="Lee, Ww" uniqKey="Lee W">WW Lee</name>
</author>
<author><name sortKey="Simarmata, D" uniqKey="Simarmata D">D Simarmata</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cho, B" uniqKey="Cho B">B Cho</name>
</author>
<author><name sortKey="Jeon, By" uniqKey="Jeon B">BY Jeon</name>
</author>
<author><name sortKey="Kim, J" uniqKey="Kim J">J Kim</name>
</author>
<author><name sortKey="Noh, J" uniqKey="Noh J">J Noh</name>
</author>
<author><name sortKey="Park, M" uniqKey="Park M">M Park</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cho, B" uniqKey="Cho B">B Cho</name>
</author>
<author><name sortKey="Kim, J" uniqKey="Kim J">J Kim</name>
</author>
<author><name sortKey="Cho, Je" uniqKey="Cho J">JE Cho</name>
</author>
<author><name sortKey="Jeon, By" uniqKey="Jeon B">BY Jeon</name>
</author>
<author><name sortKey="Park, S" uniqKey="Park S">S Park</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Yap, G" uniqKey="Yap G">G Yap</name>
</author>
<author><name sortKey="Pok, Ky" uniqKey="Pok K">KY Pok</name>
</author>
<author><name sortKey="Lai, Yl" uniqKey="Lai Y">YL Lai</name>
</author>
<author><name sortKey="Hapuarachchi, Hc" uniqKey="Hapuarachchi H">HC Hapuarachchi</name>
</author>
<author><name sortKey="Chow, A" uniqKey="Chow A">A Chow</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Greene, Ip" uniqKey="Greene I">IP Greene</name>
</author>
<author><name sortKey="Wang, E" uniqKey="Wang E">E Wang</name>
</author>
<author><name sortKey="Deardorff, Er" uniqKey="Deardorff E">ER Deardorff</name>
</author>
<author><name sortKey="Milleron, R" uniqKey="Milleron R">R Milleron</name>
</author>
<author><name sortKey="Domingo, E" uniqKey="Domingo E">E Domingo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Taylor, Wp" uniqKey="Taylor W">WP Taylor</name>
</author>
<author><name sortKey="Marshall, Id" uniqKey="Marshall I">ID Marshall</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Vasilakis, N" uniqKey="Vasilakis N">N Vasilakis</name>
</author>
<author><name sortKey="Deardorff, Er" uniqKey="Deardorff E">ER Deardorff</name>
</author>
<author><name sortKey="Kenney, Jl" uniqKey="Kenney J">JL Kenney</name>
</author>
<author><name sortKey="Rossi, Sl" uniqKey="Rossi S">SL Rossi</name>
</author>
<author><name sortKey="Hanley, Ka" uniqKey="Hanley K">KA Hanley</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Weaver, Sc" uniqKey="Weaver S">SC Weaver</name>
</author>
<author><name sortKey="Brault, Ac" uniqKey="Brault A">AC Brault</name>
</author>
<author><name sortKey="Kang, W" uniqKey="Kang W">W Kang</name>
</author>
<author><name sortKey="Holland, Jj" uniqKey="Holland J">JJ Holland</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zarate, S" uniqKey="Zarate S">S Zarate</name>
</author>
<author><name sortKey="Novella, Is" uniqKey="Novella I">IS Novella</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Malet, H" uniqKey="Malet H">H Malet</name>
</author>
<author><name sortKey="Coutard, B" uniqKey="Coutard B">B Coutard</name>
</author>
<author><name sortKey="Jamal, S" uniqKey="Jamal S">S Jamal</name>
</author>
<author><name sortKey="Dutartre, H" uniqKey="Dutartre H">H Dutartre</name>
</author>
<author><name sortKey="Papageorgiou, N" uniqKey="Papageorgiou N">N Papageorgiou</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Coffin, Jm" uniqKey="Coffin J">JM Coffin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Li, Yg" uniqKey="Li Y">YG Li</name>
</author>
<author><name sortKey="Siripanyaphinyo, U" uniqKey="Siripanyaphinyo U">U Siripanyaphinyo</name>
</author>
<author><name sortKey="Tumkosit, U" uniqKey="Tumkosit U">U Tumkosit</name>
</author>
<author><name sortKey="Noranate, N" uniqKey="Noranate N">N Noranate</name>
</author>
<author><name sortKey="A Nuegoonpipat, A" uniqKey="A Nuegoonpipat A">A A-nuegoonpipat</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Brackney, De" uniqKey="Brackney D">DE Brackney</name>
</author>
<author><name sortKey="Scott, Jc" uniqKey="Scott J">JC Scott</name>
</author>
<author><name sortKey="Sagawa, F" uniqKey="Sagawa F">F Sagawa</name>
</author>
<author><name sortKey="Woodward, Je" uniqKey="Woodward J">JE Woodward</name>
</author>
<author><name sortKey="Miller, Na" uniqKey="Miller N">NA Miller</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wikan, N" uniqKey="Wikan N">N Wikan</name>
</author>
<author><name sortKey="Sakoonwatanyoo, P" uniqKey="Sakoonwatanyoo P">P Sakoonwatanyoo</name>
</author>
<author><name sortKey="Ubol, S" uniqKey="Ubol S">S Ubol</name>
</author>
<author><name sortKey="Yoksan, S" uniqKey="Yoksan S">S Yoksan</name>
</author>
<author><name sortKey="Smith, Dr" uniqKey="Smith D">DR Smith</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cuevas, Jm" uniqKey="Cuevas J">JM Cuevas</name>
</author>
<author><name sortKey="Domingo Calap, P" uniqKey="Domingo Calap P">P Domingo-Calap</name>
</author>
<author><name sortKey="Sanjuan, R" uniqKey="Sanjuan R">R Sanjuan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chen, Sl" uniqKey="Chen S">SL Chen</name>
</author>
<author><name sortKey="Lee, W" uniqKey="Lee W">W Lee</name>
</author>
<author><name sortKey="Hottes, Ak" uniqKey="Hottes A">AK Hottes</name>
</author>
<author><name sortKey="Shapiro, L" uniqKey="Shapiro L">L Shapiro</name>
</author>
<author><name sortKey="Mcadams, Hh" uniqKey="Mcadams H">HH McAdams</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Rima, Bk" uniqKey="Rima B">BK Rima</name>
</author>
<author><name sortKey="Mcferran, Nv" uniqKey="Mcferran N">NV McFerran</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zhou, J" uniqKey="Zhou J">J Zhou</name>
</author>
<author><name sortKey="Liu, Wj" uniqKey="Liu W">WJ Liu</name>
</author>
<author><name sortKey="Peng, Sw" uniqKey="Peng S">SW Peng</name>
</author>
<author><name sortKey="Sun, Xy" uniqKey="Sun X">XY Sun</name>
</author>
<author><name sortKey="Frazer, I" uniqKey="Frazer I">I Frazer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Goodfellow, I" uniqKey="Goodfellow I">I Goodfellow</name>
</author>
<author><name sortKey="Chaudhry, Y" uniqKey="Chaudhry Y">Y Chaudhry</name>
</author>
<author><name sortKey="Richardson, A" uniqKey="Richardson A">A Richardson</name>
</author>
<author><name sortKey="Meredith, J" uniqKey="Meredith J">J Meredith</name>
</author>
<author><name sortKey="Almond, Jw" uniqKey="Almond J">JW Almond</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tuplin, A" uniqKey="Tuplin A">A Tuplin</name>
</author>
<author><name sortKey="Evans, Dj" uniqKey="Evans D">DJ Evans</name>
</author>
<author><name sortKey="Buckley, A" uniqKey="Buckley A">A Buckley</name>
</author>
<author><name sortKey="Jones, Im" uniqKey="Jones I">IM Jones</name>
</author>
<author><name sortKey="Gould, Ea" uniqKey="Gould E">EA Gould</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ou, Jh" uniqKey="Ou J">JH Ou</name>
</author>
<author><name sortKey="Strauss, Eg" uniqKey="Strauss E">EG Strauss</name>
</author>
<author><name sortKey="Strauss, Jh" uniqKey="Strauss J">JH Strauss</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Watts, Jm" uniqKey="Watts J">JM Watts</name>
</author>
<author><name sortKey="Dang, Kk" uniqKey="Dang K">KK Dang</name>
</author>
<author><name sortKey="Gorelick, Rj" uniqKey="Gorelick R">RJ Gorelick</name>
</author>
<author><name sortKey="Leonard, Cw" uniqKey="Leonard C">CW Leonard</name>
</author>
<author><name sortKey="Bess, Jw" uniqKey="Bess J">JW Bess</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Song, Y" uniqKey="Song Y">Y Song</name>
</author>
<author><name sortKey="Liu, Y" uniqKey="Liu Y">Y Liu</name>
</author>
<author><name sortKey="Ward, Cb" uniqKey="Ward C">CB Ward</name>
</author>
<author><name sortKey="Mueller, S" uniqKey="Mueller S">S Mueller</name>
</author>
<author><name sortKey="Futcher, B" uniqKey="Futcher B">B Futcher</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wichman, Ha" uniqKey="Wichman H">HA Wichman</name>
</author>
<author><name sortKey="Badgett, Mr" uniqKey="Badgett M">MR Badgett</name>
</author>
<author><name sortKey="Scott, La" uniqKey="Scott L">LA Scott</name>
</author>
<author><name sortKey="Boulianne, Cm" uniqKey="Boulianne C">CM Boulianne</name>
</author>
<author><name sortKey="Bull, Jj" uniqKey="Bull J">JJ Bull</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Igarashi, A" uniqKey="Igarashi A">A Igarashi</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zhang, Sp" uniqKey="Zhang S">SP Zhang</name>
</author>
<author><name sortKey="Zubay, G" uniqKey="Zubay G">G Zubay</name>
</author>
<author><name sortKey="Goldman, E" uniqKey="Goldman E">E Goldman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tsetsarkin, K" uniqKey="Tsetsarkin K">K Tsetsarkin</name>
</author>
<author><name sortKey="Higgs, S" uniqKey="Higgs S">S Higgs</name>
</author>
<author><name sortKey="Mcgee, Ce" uniqKey="Mcgee C">CE McGee</name>
</author>
<author><name sortKey="De Lamballerie, X" uniqKey="De Lamballerie X">X De Lamballerie</name>
</author>
<author><name sortKey="Charrel, Rn" uniqKey="Charrel R">RN Charrel</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Baronti, C" uniqKey="Baronti C">C Baronti</name>
</author>
<author><name sortKey="Sire, J" uniqKey="Sire J">J Sire</name>
</author>
<author><name sortKey="De Lamballerie, X" uniqKey="De Lamballerie X">X de Lamballerie</name>
</author>
<author><name sortKey="Querat, G" uniqKey="Querat G">G Querat</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hall, Ta" uniqKey="Hall T">TA Hall</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Thompson, Jd" uniqKey="Thompson J">JD Thompson</name>
</author>
<author><name sortKey="Higgins, Dg" uniqKey="Higgins D">DG Higgins</name>
</author>
<author><name sortKey="Gibson, Tj" uniqKey="Gibson T">TJ Gibson</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>
<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="Clement, M" uniqKey="Clement M">M Clement</name>
</author>
<author><name sortKey="Posada, D" uniqKey="Posada D">D Posada</name>
</author>
<author><name sortKey="Crandall, Ka" uniqKey="Crandall K">KA Crandall</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Peris, Jb" uniqKey="Peris J">JB Peris</name>
</author>
<author><name sortKey="Davis, P" uniqKey="Davis P">P Davis</name>
</author>
<author><name sortKey="Cuevas, Jm" uniqKey="Cuevas J">JM Cuevas</name>
</author>
<author><name sortKey="Nebot, Mr" uniqKey="Nebot M">MR Nebot</name>
</author>
<author><name sortKey="Sanjuan, R" uniqKey="Sanjuan R">R Sanjuan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Richardson, Ba" uniqKey="Richardson B">BA Richardson</name>
</author>
<author><name sortKey="Overbaugh, J" uniqKey="Overbaugh J">J Overbaugh</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 Pathog</journal-id>
<journal-id journal-id-type="iso-abbrev">PLoS Pathog</journal-id>
<journal-id journal-id-type="publisher-id">plos</journal-id>
<journal-id journal-id-type="pmc">plospath</journal-id>
<journal-title-group><journal-title>PLoS Pathogens</journal-title>
</journal-title-group>
<issn pub-type="ppub">1553-7366</issn>
<issn pub-type="epub">1553-7374</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">23436995</article-id>
<article-id pub-id-type="pmc">3578757</article-id>
<article-id pub-id-type="publisher-id">PPATHOGENS-D-12-02336</article-id>
<article-id pub-id-type="doi">10.1371/journal.ppat.1003172</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>Microbiology</subject>
<subj-group><subject>Virology</subject>
<subj-group><subject>Viral Evolution</subject>
<subject>Viral Vaccines</subject>
</subj-group>
</subj-group>
<subj-group><subject>Emerging Infectious Diseases</subject>
</subj-group>
</subj-group>
</subj-group>
<subj-group subj-group-type="Discipline-v2"><subject>Medicine</subject>
<subj-group><subject>Infectious Diseases</subject>
<subj-group><subject>Neglected Tropical Diseases</subject>
<subj-group><subject>Arboviral Infections</subject>
</subj-group>
</subj-group>
<subj-group><subject>Infectious Disease Control</subject>
</subj-group>
</subj-group>
</subj-group>
</article-categories>
<title-group><article-title>Random Codon Re-encoding Induces Stable Reduction of Replicative Fitness of Chikungunya Virus in Primate and Mosquito Cells</article-title>
<alt-title alt-title-type="running-head">Random Codon Re-encoding on Chikungunya Virus</alt-title>
</title-group>
<contrib-group><contrib contrib-type="author"><name><surname>Nougairede</surname>
<given-names>Antoine</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
<xref ref-type="corresp" rid="cor1"><sup>*</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>De Fabritus</surname>
<given-names>Lauriane</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Aubry</surname>
<given-names>Fabien</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Gould</surname>
<given-names>Ernest A.</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Holmes</surname>
<given-names>Edward C.</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff3"><sup>3</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>de Lamballerie</surname>
<given-names>Xavier</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1"><label>1</label>
<addr-line>Aix Marseille Univ, IRD French Institute of Research for Development, EHESP French School of Public Health, UMR_D 190 “Emergence des Pathologies Virales,” Marseille, France</addr-line>
</aff>
<aff id="aff2"><label>2</label>
<addr-line>Sydney Emerging Infections and Biosecurity Institute, School of Biological Sciences and Sydney Medical School, The University of Sydney, Sydney, New South Wales, Australia</addr-line>
</aff>
<aff id="aff3"><label>3</label>
<addr-line>Fogarty International Center, National Institutes of Health, Bethesda, Maryland, United States of America</addr-line>
</aff>
<contrib-group><contrib contrib-type="editor"><name><surname>Vignuzzi</surname>
<given-names>Marco</given-names>
</name>
<role>Editor</role>
<xref ref-type="aff" rid="edit1"></xref>
</contrib>
</contrib-group>
<aff id="edit1"><addr-line>Institut Pasteur, France</addr-line>
</aff>
<author-notes><corresp id="cor1">* E-mail: <email>antoine.nougairede@univ-amu.fr</email>
</corresp>
<fn fn-type="COI-statement"><p>The authors have declared that no competing interests exist.</p>
</fn>
<fn fn-type="con"><p>Conceived and designed the experiments: XDL AN. Performed the experiments: AN LDF FA. Analyzed the data: AN LDF FA EAG ECH. Contributed reagents/materials/analysis tools: AN LDF FA ECH. Conducted the experiments: AN. Wrote the paper: AN. Supervised the manuscript writing: EAG ECH XDL.</p>
</fn>
</author-notes>
<pub-date pub-type="collection"><month>2</month>
<year>2013</year>
</pub-date>
<pub-date pub-type="epub"><day>21</day>
<month>2</month>
<year>2013</year>
</pub-date>
<volume>9</volume>
<issue>2</issue>
<elocation-id>e1003172</elocation-id>
<history><date date-type="received"><day>18</day>
<month>9</month>
<year>2012</year>
</date>
<date date-type="accepted"><day>18</day>
<month>12</month>
<year>2012</year>
</date>
</history>
<permissions><copyright-statement>© 2013 Nougairede et al</copyright-statement>
<copyright-year>2013</copyright-year>
<copyright-holder>Nougairede et al</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><license-p>This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are properly credited.</license-p>
</license>
</permissions>
<abstract><p>Large-scale codon re-encoding represents a powerful method of attenuating viruses to generate safe and cost-effective vaccines. In contrast to specific approaches of codon re-encoding which modify genome-scale properties, we evaluated the effects of random codon re-encoding on the re-emerging human pathogen Chikungunya virus (CHIKV), and assessed the stability of the resultant viruses during serial <italic>in cellulo</italic>
passage. Using different combinations of three 1.4 kb randomly re-encoded regions located throughout the CHIKV genome six codon re-encoded viruses were obtained. Introducing a large number of slightly deleterious synonymous mutations reduced the replicative fitness of CHIKV in both primate and arthropod cells, demonstrating the impact of synonymous mutations on fitness. Decrease of replicative fitness correlated with the extent of re-encoding, an observation that may assist in the modulation of viral attenuation. The wild-type and two re-encoded viruses were passaged 50 times either in primate or insect cells, or in each cell line alternately. These viruses were analyzed using detailed fitness assays, complete genome sequences and the analysis of intra-population genetic diversity. The response to codon re-encoding and adaptation to culture conditions occurred simultaneously, resulting in significant replicative fitness increases for both re-encoded and wild type viruses. Importantly, however, the most re-encoded virus failed to recover its replicative fitness. Evolution of these viruses in response to codon re-encoding was largely characterized by the emergence of both synonymous and non-synonymous mutations, sometimes located in genomic regions other than those involving re-encoding, and multiple convergent and compensatory mutations. However, there was a striking absence of codon reversion (<0.4%). Finally, multiple mutations were rapidly fixed in primate cells, whereas mosquito cells acted as a brake on evolution. In conclusion, random codon re-encoding provides important information on the evolution and genetic stability of CHIKV viruses and could be exploited to develop a safe, live attenuated CHIKV vaccine.</p>
</abstract>
<abstract abstract-type="summary"><title>Author Summary</title>
<p>Emerging arthropod-borne viruses (arboviruses) are a major cause of human and animal morbidity and mortality. Climatic and anthropological activities are responsible for the dispersal of arbovirus transmission vectors into new territories. Chikungunya virus (CHIKV) is an important example of a re-emerging pathogen for which no licensed vaccine exists. One of the vectors of CHIKV, the mosquito <italic>Aedes albopictus</italic>
, has dispersed into new temperate regions resulting in outbreaks where they had not been previously observed. Here, we demonstrate that random codon re-encoding, a method that modifies the nucleic acid composition of large coding regions without modifying the encoded proteins, can significantly decrease the replicative fitness of CHIKV. This powerful method of attenuating viruses has several potential advantages for vaccine development, including the possibility to modulate precisely the degree of replicative fitness loss and to generate safe, live-attenuated vaccines that confer long-term protection, in a cost effective manner. Our studies also demonstrate that these re-encoded viruses exhibit a stable phenotype, and that the response to codon re-encoding was largely compensatory in nature, with little reversion of mutations. Finally, we provide further evidence that many synonymous sites in RNA viruses are not neutral and clearly impact viral fitness.</p>
</abstract>
<funding-group><funding-statement>This work was supported by Aix-Marseille University. 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="18"></page-count>
</counts>
</article-meta>
</front>
<body><sec id="s1"><title>Introduction</title>
<p>Many emerging infectious diseases are caused by arthropod-borne viruses (arboviruses), almost all of which are single strand RNA viruses. The major outbreaks of dengue fever <xref rid="ppat.1003172-Guzman1" ref-type="bibr">[1]</xref>
, West Nile encephalitis <xref rid="ppat.1003172-Lanciotti1" ref-type="bibr">[2]</xref>
, Chikungunya fever <xref rid="ppat.1003172-Staples1" ref-type="bibr">[3]</xref>
, and Rift Valley fever <xref rid="ppat.1003172-Boshra1" ref-type="bibr">[4]</xref>
that have occurred in recent decades, each with a significant impact on human health, highlight the urgent need to understand the factors that allow these viruses to invade new territories or adapt to new host or vector species <xref rid="ppat.1003172-Randolph1" ref-type="bibr">[5]</xref>
–<xref rid="ppat.1003172-Weaver1" ref-type="bibr">[6]</xref>
. Understanding the factors that shape the adaptability of these rapidly evolving infectious agents may provide new opportunities for their eventual control.</p>
<p>Codon usage bias is an important indicator of the evolutionary forces shaping genomes and could arise either through neutral mutational pressure or because specific synonymous codons are selectively advantageous; for example, by increasing the efficiency and/or accuracy of protein expression by maximizing the match to cellular tRNA abundance. Thus, determining the underlying causes of codon bias has become a key topic in evolutionary genetics <xref rid="ppat.1003172-Hershberg1" ref-type="bibr">[7]</xref>
. RNA viruses often exhibit codon biases that match the nucleotide biases across viral genomes as a whole, suggesting that background mutational pressure is the dominant factor shaping codon choice. However, natural selection may still act at the scale of overall nucleotide composition <xref rid="ppat.1003172-Jenkins1" ref-type="bibr">[8]</xref>
. Direct selection for specific codon biases has also been documented. For example, in hepatitis A virus rare codons that utilize non-abundant tRNAs are preferred, slowing down the translation process to ensure proper protein folding <xref rid="ppat.1003172-Sanchez1" ref-type="bibr">[9]</xref>
. The large-scale re-encoding of codon usage in poliovirus, influenza A virus and bacterial virus T7 by reverse genetic methods resulted in virus attenuation <xref rid="ppat.1003172-Mueller1" ref-type="bibr">[10]</xref>
–<xref rid="ppat.1003172-Bull1" ref-type="bibr">[15]</xref>
, demonstrating that mutations at synonymous sites can indeed have a major impact on viral fitness. To date, all studies of codon re-coding have employed a specific approach such as codon de-optimisation <xref rid="ppat.1003172-Burns2" ref-type="bibr">[12]</xref>
, <xref rid="ppat.1003172-Mueller2" ref-type="bibr">[14]</xref>
, codon pair de-optimisation <xref rid="ppat.1003172-Mueller1" ref-type="bibr">[10]</xref>
, <xref rid="ppat.1003172-Coleman1" ref-type="bibr">[13]</xref>
and increase of CpG/UpA dinucleotide frequency <xref rid="ppat.1003172-Burns1" ref-type="bibr">[11]</xref>
; but all result in a reduction of viral fitness. Currently, there are no studies of codon re-coding in arboviruses. However, studies of these viruses are of special interest because their evolution is strongly constrained by host alternation, such that mutations that may be advantageous in one host type (e.g. mosquitoes or mammals) are deleterious in another <xref rid="ppat.1003172-Coffey1" ref-type="bibr">[16]</xref>
.</p>
<p>Chikungunya virus (CHIKV; <italic>Togaviridae</italic>
; <italic>Alphavirus</italic>
) is a small (60–70 nm), enveloped, single-strand positive-sense RNA virus. Its genome of approximately 12 kb contains two open reading frames (ORFs) encoding non-structural and structural proteins, respectively <xref rid="ppat.1003172-Strauss1" ref-type="bibr">[17]</xref>
. In Swahili, “Chikungunya” means “bent walker”, reflecting the severe arthralgia associated with CHIKV infections. First isolated in Tanzania in 1952 <xref rid="ppat.1003172-Robinson1" ref-type="bibr">[18]</xref>
, CHIKV is transmitted by mosquito vectors of the <italic>Aedes (Stegomyia)</italic>
subgenus and has caused a number of outbreaks in Africa and Asia during the last 50 years <xref rid="ppat.1003172-Schuffenecker1" ref-type="bibr">[19]</xref>
. It is believed that the original natural history of viral transmission relies on virus maintenance in a yellow fever-like zoonotic sylvatic cycle involving non-peridomestic mosquitoes and nonhuman primates, as previously described in Africa. However, explosive urban outbreaks were associated with a dengue-like direct “human-mosquito-human” transmission cycle implicating <italic>A. aegypti</italic>
or more recently <italic>A. albopictus</italic>
mosquitoes <xref rid="ppat.1003172-Pialoux1" ref-type="bibr">[20]</xref>
. Particularly large CHIKV outbreaks have occurred in Indian Ocean islands, in India, and in Southeast Asia since 2005 <xref rid="ppat.1003172-Weaver1" ref-type="bibr">[6]</xref>
. All these epidemics originated from east Africa and were associated with the East-Central-South African genotype <xref rid="ppat.1003172-Volk1" ref-type="bibr">[21]</xref>
. Significantly, these recent epidemics are also associated with viral transmission by <italic>A. albopictus</italic>
, although this was previously considered a secondary vector, and convergent adaptation to this mosquito has been observed in different geographical regions <xref rid="ppat.1003172-Tsetsarkin1" ref-type="bibr">[22]</xref>
–<xref rid="ppat.1003172-deLamballerie1" ref-type="bibr">[23]</xref>
. Approximately 40% of the population of Reunion Island was infected (around 300,000 people) <xref rid="ppat.1003172-Gerardin1" ref-type="bibr">[24]</xref>
, including a relatively small proportion of severe disease in adults and newborns, as well as some vertical transmission <xref rid="ppat.1003172-Economopoulou1" ref-type="bibr">[25]</xref>
–<xref rid="ppat.1003172-Gerardin2" ref-type="bibr">[26]</xref>
. Worryingly, the increasingly widespread distribution of <italic>A. albopictus</italic>
may result in CHIKV epidemics in more temperate regions. Indeed, in 2007, a CHIKV outbreak occurred in Italy, and the detection of two autochthonous cases in France in 2010 confirmed fears of the possible expansion of this important viral disease <xref rid="ppat.1003172-Gould1" ref-type="bibr">[27]</xref>
–<xref rid="ppat.1003172-Rezza1" ref-type="bibr">[28]</xref>
, particularly as there is currently no commercialized antiviral or vaccine for this virus.</p>
<p>Codon re-encoding may represent an important tool for the design of effective vaccines for CHIKV and other arboviruses. However, for such a strategy to succeed it is essential to determine how the virus might respond to this profound change in selection pressure, particularly given its reliance on both mosquitoes and mammals for transmission. Key questions include: to what extent will codon re-encoding reduce viral fitness, how rapidly will CHIKV recover fitness following re-encoding, and will this fitness recovery involve direct reversion at synonymous sites? To address these questions, we studied the <italic>in cellulo</italic>
replicative fitness of codon re-encoded CHIKV and the <italic>in cellulo</italic>
evolution of re-encoded viruses by combinations of either alternate or continuous passage of each virus in primate or insect cells.</p>
</sec>
<sec id="s2"><title>Results</title>
<p>Using a large scale random re-encoding method which randomly attributed nucleotide codons based on their corresponding amino acid sequence (e.g. the amino acid valine was randomly replaced by GTT, GTC, GTA or GTG; see <xref ref-type="sec" rid="s4">materials and methods</xref>
for more details), we designed, synthesized, and incorporated into a WT infectious clone (IC; derived from the LR2006 CHIKV strain) three re-encoded cassettes of around 1.4 Kb, located in the non-structural proteins nsP1, nsP4, and in a region overlapping the structural proteins E2 and E1. A total of 264, 298 and 320 synonymous mutations were present in these re-encoded regions, respectively. Six re-encoded viruses were studied using combinations of these re-encoded regions: Φnsp1, Φnsp4 and Φenv with one re-encoded region; Φnsp1 Φnsp4, and Φnsp4 Φenv with two re-encoded regions and Φnsp1 Φnsp4 Φenv with three re-encoded regions (<bold>Figure S1 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
). Their genetic characteristics, such as G+C%, codon usage bias, frequencies of rare codons and of the dinucleotides CpG/UpA, remained comparable with WT virus and 132 other CHIKVs from GenBank (<bold>Table S1 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
).</p>
<sec id="s2a"><title>Replicative fitness of the re-encoded viruses</title>
<p>The WT virus and the six re-encoded viruses were derived following transfection of the corresponding infectious DNA clones into Vero cells. Viruses were then passaged once in Vero cells and their replicative fitness was studied. All the viruses produced a cytopathic effect (CPE) which was delayed proportionally with the degree of re-encoding (from 2 days for the WT virus to 5–6 days for the Φnsp1 Φnsp4 Φenv virus).</p>
<sec id="s2a1"><title>(i) Single cycle replication kinetics and plaque morphology</title>
<p>Using primate (Vero and HEK293) and mosquito cell lines (C6/36), replication kinetics were investigated using a high estimated multiplicity of infection (MOI = 5) (<xref ref-type="fig" rid="ppat-1003172-g001"><bold>Figure 1</bold>
</xref>
) during the first 28 hours post-infection (pi). Quantitative analyses were performed at 8 and 14 hours pi (<bold>Table S2 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
), i.e. a period shorter than two complete replication cycles on the basis of 8 hours per replication cycle <xref rid="ppat.1003172-Strauss1" ref-type="bibr">[17]</xref>
.</p>
<fig id="ppat-1003172-g001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1003172.g001</object-id>
<label>Figure 1</label>
<caption><title>Single cycle replication kinetics and plaque morphologies using WT and re-encoded viruses.</title>
<p>To study replicative fitness of the re-encoded viruses, all viruses (including the WT) were used to perform single cycle replication kinetics in Vero (<bold>A–B</bold>
) HEK293 (<bold>C–D</bold>
) and C6/36 cells (<bold>E–F</bold>
). An estimated MOI of 5 was used to infect confluent cells. Cells were washed 30 minutes after the infection and 1 ml of cell supernatant was sampled at several time pi. Clarified (centrifugation) cell supernatants were used to measure (i) the infectious titres using a TCID50 assay (<bold>A, C, E</bold>
) and (ii) the amount of viral RNA using a real time RT-PCR assay (<bold>B, D, F</bold>
). The plaque morphology of these viruses was also analysed comparatively in Vero cells (<bold>G</bold>
).</p>
</caption>
<graphic xlink:href="ppat.1003172.g001"></graphic>
</fig>
<p>When analysis of virus replication in primate or mosquito cells was based on infectious titres in cell supernatants (as measured using TCID50 method), an inverse linear relationship was observed between infectious titres at 14 hours pi and the number of synonymous mutations in re-encoded regions (<bold>Figure S2 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
). Detailed results are reported in <bold>Table S2 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
and summarised below: For the Φnsp1 Φnsp4 Φenv virus, infectious titres were significantly lower than those of the wild type virus. The difference reached 1 to 1.5 and 2 to 2.5 log<sub>10</sub>
TCID50/ml at 8 and 14 hours pi, respectively. For viruses with two regions re-encoded (Φnsp1 Φnsp4 and Φnsp1 Φenv), significant infectious titre reductions compared with WT virus (0.5 to1 log<sub>10</sub>
TCID50/ml) were observed for all constructs propagated in C6/36 cells and for the Φnsp1 Φnsp4 virus in Vero cells at 8 hours pi. All titres were significantly lower at 14 hours pi (1–1.5 log<sub>10</sub>
TCID50/ml). For viruses with one region re-encoded (Φnsp1, Φnsp4 and Φenv), significant infectious titre reductions, compared with WT, were observed for the Φenv virus in HEK293 cells (∼0.5–1 log<sub>10</sub>
TCID50/ml at 8 and 14 hours pi) and the Φnsp1 and Φenv viruses in C6/36 cells (∼0.75–1.25 log<sub>10</sub>
TCID50/ml at 14 hours pi). However, compared with WT virus, there were no significant reductions of infectious titres in Vero cells at any time.</p>
<p>When replication was assessed on viral RNA yields (as measured using quantitative real time RT-PCR and cell supernatants), strikingly different results were observed in mosquito and primate cells (<xref ref-type="fig" rid="ppat-1003172-g001"><bold>Figure 1</bold>
</xref>
and <bold>Figure S2 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
). In C6/36 mosquito cells, viral RNA yields remained identical regardless of the degree of re-encoding. In primate cells, an inverse linear relationship between RNA yield at 14 hours pi and the number of synonymous mutations in re-encoded region(s) was observed.</p>
<p>To assess the relationship between viral RNA yields and the number of viral particles (i.e. infectious and non infectious particles) in the clarified supernatants, results of the quantitative RT-PCR assay, the TCID50 assay and a haemagglutination assay were compared. The latter provided an estimate of the relative quantity of CHIKV haemagglutinin. This, in turn, reflects the number of viral particles <xref rid="ppat.1003172-Casals1" ref-type="bibr">[29]</xref>
–<xref rid="ppat.1003172-Clarke1" ref-type="bibr">[30]</xref>
, and provides a direct comparison between WT and re-encoded viruses which encode identical proteins. At 16 hours pi, using an estimated MOI of 5, viral RNA yields and haemagglutinin titres from both Vero and C6/36 cell supernatants showed parallel production kinetics. On the other hand, TCID50 titres progressively decreased with increasing levels of re-encoding in both primate and mosquito cells (<xref ref-type="fig" rid="ppat-1003172-g002"><bold>Figure 2</bold>
</xref>
). This result implies that the real-time RT-PCR assay reflects the number of viral particles and hence that the re-encoding does not appear to impact significantly on the incorporation of viral RNA into viral particles.</p>
<fig id="ppat-1003172-g002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1003172.g002</object-id>
<label>Figure 2</label>
<caption><title>Comparison of the results from haemagglutination, real time RT-PCR and TCID50 assays.</title>
<p>To assess the relationship between viral RNA yields and the number of viral particles in the cell supernatants, the results of quantitative RT-PCR, TCID50 and haemagglutination assay were compared. An estimated MOI of 5 was used to infect confluent Vero (<bold>A</bold>
) or C6/36 (<bold>B</bold>
) cells. Cells were washed after the infection and cell supernatants were sampled at 16 hours pi. Clarified (centrifugation) cell supernatants were used to perform haemagglutination, real time RT-PCR and TCID50 assays. Results are expressed in the graphs as percentage of WT and are the mean/standard deviation of three independent experiments.</p>
</caption>
<graphic xlink:href="ppat.1003172.g002"></graphic>
</fig>
<p>In C6/36 cells, as the degree of codon re-encoding increased, the number of viral particles remained relatively stable. However, infectious titre values gradually decreased up to 2.5 log<sub>10</sub>
TCID50/ml (<xref ref-type="fig" rid="ppat-1003172-g001"><bold>Figure 1</bold>
</xref>
<bold> and Table S2 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
). This indicates that re-encoding was associated with a significant modification of the infectivity of viral particles (i.e. the ratio of the number of infectious particles [calculated here using a TCID50 assay] over the number of viral particles [calculated here using a quantitative PCR assay]). When compared with the WT virus, we observed a 500 fold decrease of infectivity for the most re-encoded (Φnsp1 Φnsp4 Φenv) virus, 30–50 fold for viruses with two regions re-encoded and 5–25 times for viruses with one region re-encoded (<bold>Table S2 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
).</p>
<p>Primate cells also showed a decrease of viral infectivity as the degree of re-encoding increased. However, this level of decrease was more limited, reaching a 10–20 fold reduction for the most re-encoded virus (<bold>Table S2 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
). By contrast with the C6/36 cells, the production of viral particles was significantly decreased at 14 hours pi for viruses with 3 regions (∼1 log<sub>10</sub>
RNA copies/ml) and 2 regions re-encoded (0.5–1 log<sub>10</sub>
RNA copies/ml) (<xref ref-type="fig" rid="ppat-1003172-g001"><bold>Figure 1</bold>
</xref>
<bold>and Table S2 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
). In addition, each virus with one region re-encoded showed a significant reduction at 14 hours pi in one primate cell type at least (approximately 0.5 log<sub>10</sub>
RNA copies/ml). Therefore, whilst the major modification observed in C6/36 cells was a reduction of virus infectivity yields, in primate cells there was also a reduction of the number of viral particles.</p>
<p>The plaque morphology in Vero cells confirmed the results of the single cycle replication kinetics: globally, as the degree of re-encoding increased, there was a corresponding decrease in plaque size (<xref ref-type="fig" rid="ppat-1003172-g001"><bold>Figure 1</bold>
</xref>
). Interestingly, the Φenv virus plaque size most closely resembled that of viruses with two re-encoded regions rather than Φnsp1 or Φnsp4 virus.</p>
</sec>
<sec id="s2a2"><title>(ii) Replication kinetics with low MOI viral infection or infectious DNA clones</title>
<p>Defective interfering (DI) particles are known to impact on viral replication cycles <xref rid="ppat.1003172-Roux1" ref-type="bibr">[31]</xref>
. To assess their potential role in the differences observed in our single cycle replication studies, we analysed the kinetics of virus production using a low estimated MOI (0.01) in Vero and C6/36 cells (<xref ref-type="fig" rid="ppat-1003172-g003"><bold>Figure 3</bold>
</xref>
). The infectious titre values decreased with increasing degrees of re-encoding, regardless of the observation time pi, indicating that DI particles did not play a major role in the effect observed using high estimated MOI. Globally, the degree of reduction of virus infectivity in comparison with the WT virus was higher than that observed in the single cycle replication kinetics (approximately 3–4 log<sub>10</sub>
TCID50/ml for the most re-encoded virus 24 hours pi). This presumably reflects the fact that these analyses were performed following several viral replication cycles, thus amplifying the difference observed in the single-cycle experiments reported earlier (<xref ref-type="fig" rid="ppat-1003172-g003"><bold>Figure 3</bold>
</xref>
).</p>
<fig id="ppat-1003172-g003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1003172.g003</object-id>
<label>Figure 3</label>
<caption><title>Replication kinetics using low MOI with WT and re-encoded viruses.</title>
<p>An estimated MOI of 0.01 was used to infect confluent Vero (<bold>A</bold>
) or C6/36 (<bold>B</bold>
) cells. Cells were washed and cell supernatant was sampled at several times pi. Clarified (centrifugation) cell supernatants were used to measure the infectious titres using a TCID50 assay.</p>
</caption>
<graphic xlink:href="ppat.1003172.g003"></graphic>
</fig>
<p>Using a different method we confirmed that the observed effect could not be due to higher numbers of DI particles in cell supernatants derived from re-encoded viruses: ICs were directly used to transfect HEK293 cells. Infectious titres and number of viral particles were then measured in cell supernatants (<xref ref-type="fig" rid="ppat-1003172-g004"><bold>Figure 4</bold>
</xref>
). The infectious titre values decreased with the degree of re-encoding, regardless of the time pi. At 16 hours pi, the infectious titre of the most re-encoded virus was undetectable whilst values of the WT virus were around 4 log<sub>10</sub>
TCID50/ml. As observed in single cycle replication kinetics, the number of viral particles using RNA yields as the indicator of particle numbers decreased less than the infectious titres (<xref ref-type="fig" rid="ppat-1003172-g004"><bold>Figure 4</bold>
</xref>
).</p>
<fig id="ppat-1003172-g004" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1003172.g004</object-id>
<label>Figure 4</label>
<caption><title>Replication kinetics using infectious DNA clones with WT and re-encoded viruses.</title>
<p>Equal quantities of ICs were used to transfect culture flasks of HEK293 cells. Cells were washed and cell supernatant was sampled at several time pi. Clarified (centrifugation) cell supernatants were used to measure (i) the infectious titres using a TCID50 assay (<bold>A</bold>
) and (ii) the amount of viral RNA using a real time RT-PCR assay (<bold>B</bold>
).</p>
</caption>
<graphic xlink:href="ppat.1003172.g004"></graphic>
</fig>
<p>In comparison with the single cycle replication experiments, when compared with WT virus, both methods identified a significant reduction of infectious titres for each of the three viruses with one region re-encoded (Student's <italic>t</italic>
test, all <italic>p</italic>
value<0.01 with low estimated MOI 48 hours pi and DNA clones 24 hours pi). As with plaque morphology in Vero cells, we observed that the replicative fitness of the Φenv virus more closely resembled viruses with two regions re-encoded than that of Φnsp1 or Φnsp4 virus (<xref ref-type="fig" rid="ppat-1003172-g003"><bold>Figure 3</bold>
</xref>
<bold>–</bold>
<xref ref-type="fig" rid="ppat-1003172-g004"><bold>4</bold>
</xref>
).</p>
</sec>
<sec id="s2a3"><title>(iii) Competition experiments</title>
<p>Using two real time RT-PCR assays, specific for the WT virus or the viruses with the nsP4 region re-encoded (Φnsp4, Φnsp1 Φnsp4, Φnsp4 Φenv and Φnsp1 Φnsp4 Φenv), competition experiments were performed in Vero and C6/36 cells. Five initial PFU ratios (WT/re-encoded virus: 1/99, 20/80, 50/50, 80/20, 99/1) were used. As described in the Methods section, this procedure provided quantitative estimates of the proportion of each viral genome in the viral population (expressed as log<sub>10</sub>
WT/re-encoded ratio in <xref ref-type="fig" rid="ppat-1003172-g005"><bold>Figure 5</bold>
</xref>
). In addition, cell supernatants were passaged 10 times to enable long term follow-up of each viral population (viruses completed an estimated 60 replication cycles).</p>
<fig id="ppat-1003172-g005" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1003172.g005</object-id>
<label>Figure 5</label>
<caption><title>Competition experiments.</title>
<p>Using two real time RT-PCR assays, specific for the WT virus or the viruses with the nsP4 region re-encoded (Φnsp4, Φnsp1 Φnsp4, Φnsp4 Φenv and Φnsp1 Φnsp4 Φenv viruses), competition experiments were performed in Vero (<bold>A–D</bold>
) and C6/36 (<bold>E–H</bold>
) cells. Five initial PFU ratios (WT/re-encoded virus: 1/99, 20/80, 50/50, 80/20, 99/1) were used. This enabled us to quantify the proportion of each virus in the viral population (expressed as log<sub>10</sub>
WT/re-encoded ratio). In addition, cell supernatants were passaged 10 times to enable long term follow-up of each viral population.</p>
</caption>
<graphic xlink:href="ppat.1003172.g005"></graphic>
</fig>
<p>In Vero cells (<xref ref-type="fig" rid="ppat-1003172-g005"><bold>Figures 5A–D</bold>
</xref>
), WT virus was proportionally more fit than re-encoded viruses depending on the degree of re-encoding. Re-encoded viruses disappeared completely after 8–10, 5–7 and only 3–5 passages for the Φnsp4 virus, the viruses with two re-encoded regions and the most re-encoded virus respectively. Two extreme initial PFU ratios were tested in Vero cells with the Φnsp4 virus (0.1/99.9 and 0.01/99.99) and viral populations followed a similar pattern, with WT virus rapidly supplanting the re-encoded virus (Data not shown). In C6/36 cells (<xref ref-type="fig" rid="ppat-1003172-g005"><bold>Figures 5E–H</bold>
</xref>
), WT virus was also proportionally more fit than re-encoded viruses depending on the degree of re-encoding but this process occurred more slowly than in Vero cells: WT/Φnsp4 ratios remained stable during 10 passages, whereas ratios for the viruses with two re-encoded regions increased by ∼2–2.5 log<sub>10</sub>
and the Φnsp1 Φnsp4 Φenv virus was rapidly supplanted by the WT virus.</p>
</sec>
<sec id="s2a4"><title>(iv) Synthesis of viral proteins and intracellular viral RNA</title>
<p>All the experiments described above incorporated clarified infectious cell supernatant medium to study the replicative fitness of re-encoded viruses. However, we also studied these parameters at the intracellular level by assessing how the codon re-encoding method modified intracellular viral RNA and protein yields. Because the presence of stable levels of viral RNA and HA in the supernatant of the C6/36 cells indicated that RNA and protein synthesis was unaffected by re-encoding, only human cells were studied in this analysis. HEK293 cells were infected with a high estimated MOI (5.0), and intracellular RNA and proteins were extracted at 8 hours pi. Viral protein levels were measured by western blot analysis using anti-nsP1/nsP2 rabbit pAb and ELISA using a CHIKV-specific immune human serum (recognizing especially the structural proteins <xref rid="ppat.1003172-Kam1" ref-type="bibr">[32]</xref>
–<xref rid="ppat.1003172-Yap1" ref-type="bibr">[35]</xref>
; <bold>Protocol S1</bold>
). Viral RNA yields were measured using a real time RT-PCR assay (<xref ref-type="fig" rid="ppat-1003172-g006"><bold>Figure 6</bold>
</xref>
).</p>
<fig id="ppat-1003172-g006" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1003172.g006</object-id>
<label>Figure 6</label>
<caption><title>Synthesis of viral proteins and intracellular viral RNA.</title>
<p>To study, at the intracellular level, how our codon re-encoding method modified viral RNA and protein synthesis, 12-well-plates of HEK293 cells were infected with virus at an estimated MOI of 5. Cells were washed after the infection and intracellular RNA and proteins were extracted at 8 hours pi simultaneously or individually (cf. <xref ref-type="sec" rid="s4">materials and methods</xref>
). Viral protein levels were measured using two methods: a western blot analysis using anti-nsP1/nsP2 rabbit pAb (<bold>A</bold>
) and an ELISA (<bold>B</bold>
) using a CHIKV-specific immune human serum (recognizing especially the structural proteins <xref rid="ppat.1003172-Kam1" ref-type="bibr">[32]</xref>
–<xref rid="ppat.1003172-Yap1" ref-type="bibr">[35]</xref>
). In panel <bold>B</bold>
, results are expressed as percentage of WT and are the mean/standard deviation of three independent experiments. Viral RNA yields were quantified by real time comparative RT-PCR using actin mRNA as a normalizer and cells infected by the WT virus as a calibrator (<bold>C</bold>
); Results are the mean relative quantity (RQ) from three independent experiments and errors bars are RQ max and RQ min values generated by the Stratagene Mx3005P software at 99% confidence level; Results are the mean and standard deviation from three independent experiments.</p>
</caption>
<graphic xlink:href="ppat.1003172.g006"></graphic>
</fig>
<p>Both methods revealed a quantitative decrease in the level of structural and non-structural viral protein production with increasing degrees of codon re-encoding (around five times for the most re-encoded virus using ELISA). Viral RNA yields also decreased with increasing degrees of codon re-encoding (around 25 times for the most re-encoded virus) (<xref ref-type="fig" rid="ppat-1003172-g006"><bold>Figure 6</bold>
</xref>
). These results confirmed the supernatant medium analysis (<xref ref-type="fig" rid="ppat-1003172-g001"><bold>Figure 1</bold>
</xref>
) and implied that in human cells, re-encoding compromised the replication complex resulting in decreased levels of viral RNA and viral protein in cells and supernatant medium.</p>
</sec>
</sec>
<sec id="s2b"><title>Experimental passages of the WT and two re-encoded viruses</title>
<p>After their recovery by transfection, two re-encoded viruses (Φnsp4 and Φnsp1 Φnsp4 Φenv) and the WT virus were passaged using three different protocols: 50 serial passages in non-human primate (Vero) or mosquito (C6/36) cells, and 50 alternate passages in Vero and C6/36 cells (<italic>i.e.</italic>
, 25 double passages (Vero/C6/36 cells)). At the time of each passage, the estimated MOI was bottlenecked at approximately 0.1 to minimize the generation of defective interfering particles without generating a major population bottleneck. Each passage was terminated after 48 hours. It is therefore estimated that each virus completed ∼300 replication cycles after fifty passages on the basis of 8 hours per replication cycle <xref rid="ppat.1003172-Strauss1" ref-type="bibr">[17]</xref>
.</p>
<sec id="s2b1"><title>Fitness of passaged CHIKV in non-human primate and insect cells</title>
<p>To study the evolution of replicative fitness during passages in response to codon re-encoding, we examined the replicative kinetics of each passaged virus at the 1<sup>st</sup>
, 12<sup>th</sup>
, 25<sup>th</sup>
, 37<sup>th</sup>
and 50<sup>th</sup>
passage in Vero and C6/36 cells (<bold>Figure S3, S4, S5 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
). We first measured the viral growth rate at 24 hours pi, based on TCID50 values. The corresponding relative fitness effect values are detailed in <xref ref-type="fig" rid="ppat-1003172-g007"><bold>Figure 7</bold>
</xref>
. We then performed a global analysis of TCID50 values at 24, 48 and 72 hours pi by performing two-way repeated-measures ANOVA and tukey's HDS post-hoc comparisons (<bold>Table S3 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
), which provided very similar results.</p>
<fig id="ppat-1003172-g007" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1003172.g007</object-id>
<label>Figure 7</label>
<caption><title>Fitness modifications observed during the passage experiments.</title>
<p>The replicative fitness of each virus effect is represented here using relative fitness effect values. These values represent: [(the viral growth rate of a passaged virus at 24 hours pi, based on TCID50 values)/(the viral growth rate of the corresponding virus at the first passage)]−1 (See <xref ref-type="sec" rid="s4"><bold>Materials and Methods</bold>
</xref>
). The relative fitness effect was calculated using replication kinetics experiments in Vero cells (<bold>A–C</bold>
) and C6/36 cells (<bold>D–F</bold>
). Panels <bold>A</bold>
and <bold>D</bold>
represent viruses serially passaged in Vero cells, <bold>B</bold>
and <bold>E</bold>
viruses serially passaged in C6/36, and <bold>C</bold>
and <bold>F</bold>
viruses alternately passaged. Significant values (comparison made with the first passage of the same virus in the same cells, Student's <italic>t</italic>
test) are indicated by * (0.01<<italic>p</italic>
<0.05), ** (0.001<<italic>p</italic>
<0.01 and *** (<italic>p</italic>
<0.001).</p>
</caption>
<graphic xlink:href="ppat.1003172.g007"></graphic>
</fig>
<p>First, we showed that the fitness trajectories of WT virus changed with the passaging methods. This phenomenon has been widely reported in the arbovirus literature <xref rid="ppat.1003172-Greene1" ref-type="bibr">[36]</xref>
–<xref rid="ppat.1003172-Zarate1" ref-type="bibr">[40]</xref>
. Replicative fitness in Vero cells significantly increased following serial passage in Vero cells, whereas it decreased following serial passage in C6/36 cells and slightly decreased in the case of alternate passages (<xref ref-type="fig" rid="ppat-1003172-g007"><bold>Figure 7</bold>
</xref>
). Replicative fitness in C6/36 cells remained globally unchanged following serial passage in C6/36 cells, in Vero cells or both. Φnsp4 behaviour appeared to be similar to that of the WT virus, with the exception of replicative fitness in C6/36 cells which significantly increased after serial passage in these cells.</p>
<p>In the case of Φnsp1 Φnsp4 Φenv, the evolution of the replicative fitness was strikingly different: it increased in all cases except in Vero cells after serial passage in C6/36 where a significant decrease was observed (<xref ref-type="fig" rid="ppat-1003172-g007"><bold>Figure 7</bold>
</xref>
). Thus, in most cases, the replicative fitness of the re-encoded viruses increased significantly during serial passage whereas that of the WT virus remained unchanged. In addition, after 25 serial passages in Vero cells, both re-encoded viruses showed an increase of their replicative fitness in Vero cells which is significantly higher than that observed for the WT virus (Student's <italic>t</italic>
test, all <italic>p</italic>
value <0.01).</p>
<p>Despite these specific replicative fitness gains, the infectious titres of the Φnsp1 Φnsp4 Φenv virus in each cell type did not approach those obtained with the first passage of the WT virus: the smallest difference observed at 24, 48 and 72 h pi was 2.6, 1.9 and 0.7 log<sub>10</sub>
TCID50 in Vero cells, and 1.7, 0.9 and 0.5 log<sub>10</sub>
TCID50 in C6/36 cell, respectively (<bold>Figure S3, S4, S5 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
).</p>
</sec>
<sec id="s2b2"><title>Genome sequence analysis of passaged CHIKV</title>
<p>To understand which genomic modifications accompanied these changes in replicative fitness, the complete genome consensus nucleotide sequences of Φnsp4 and Φnsp1 Φnsp4 Φenv were established at passage 50. Intra-population viral genetic diversity was assessed in three coding regions (nsP2/nsP3, E2 and E1): approximately 20 clones were sequenced for each PCR product studied (mean: 19.2+/−4.2; range from 14 to 44). Because each virus passage was not repeated, and in view of the globally low number of mutations detected, no statistical analysis of the data was performed in this section.</p>
<p>In our experimental model of evolution, selection pressure applied to the re-encoded viruses arose from both codon re-encoding and adaptation to culture conditions. All the viruses, re-encoded or not, followed similar evolutionary pathways according to passaging method. First, we noted an association between the number of observed mutations in complete genome consensus sequences and the passage method: highest after serial passage in Vero cells (range 8 to 16), lowest after serial passage in C6/36 cells (range 3 to 6), and intermediate after alternate passage (range 6 to 10) (<xref ref-type="table" rid="ppat-1003172-t001"><bold>Table 1</bold>
</xref>
). Very few mutations were detected during serial passage in C6/36 cells, all of these were non-synonymous and most (85%) were located in the structural protein coding region (<xref ref-type="table" rid="ppat-1003172-t001"><bold>Table 1</bold>
</xref>
and <bold>Table S4 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
). Moreover, among the few mutations observed after the serial passages in C6/36 cells, only 1/13 became fixed in the population (based on complete genome consensus sequences), while among those observed after serial passage in Vero cells and alternate passage,19/34 and 9/23 mutations, respectively, were fixed (<xref ref-type="table" rid="ppat-1003172-t001"><bold>Table 1</bold>
</xref>
). This observation was confirmed by analysis of intra-population genetic diversity: (i) fluctuations in the size of variant sub-populations were observed with viruses serially passaged in C6/36, and (ii) mutations observed after serial passage in Vero cells and alternate passage tended to be fixed on a more regular and rapid basis (<xref ref-type="fig" rid="ppat-1003172-g008"><bold>Figure 8</bold>
</xref>
and <bold>Figure S6 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
). Furthermore, from the 6<sup>th</sup>
passage onwards, viruses serially passaged in Vero cells exhibited a continuous appearance of mutations, whereas those passaged in C6/36 cells exhibited a very low mutation frequency, and a delayed emergence of mutations was observed for alternately passaged viruses; in the latter case mutations were first detected at the 25<sup>th</sup>
passage for both the WT and Φnsp4 viruses (<xref ref-type="fig" rid="ppat-1003172-g009"><bold>Figure 9</bold>
</xref>
). Our experiments also revealed frequent convergent evolution; a total of 15 and 24 convergent mutations were observed using complete genome consensus and clonal sequences, respectively (<xref ref-type="table" rid="ppat-1003172-t002"><bold>Table 2</bold>
</xref>
). Approximately 50% of these convergent mutations were shared between re-encoded viruses and the WT virus. Notably, all those convergent mutations found in complete genome consensus sequences occurred in viruses passaged using the same method or between serially and alternatively passaged viruses (<xref ref-type="table" rid="ppat-1003172-t002"><bold>Table 2</bold>
</xref>
), denoting adaptation to cell culture conditions.</p>
<fig id="ppat-1003172-g008" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1003172.g008</object-id>
<label>Figure 8</label>
<caption><title>Intra-population genetic diversity of CHIKV.</title>
<p>Minimum spanning trees were used to depict the intra-population genetic diversity of CHIKV. Each panel represents all the clones sequenced in one region of one virus passage from one method: (<bold>A</bold>
) E3/E2 region of the Φnsp4 virus passaged in C6/36 cells, (<bold>B</bold>
) E3/E2 region of the Φnsp1 Φnsp4 Φenv virus alternatively passaged, and (<bold>C</bold>
) E1 region of the Φnsp1 Φnsp4 Φenv virus serially passaged in Vero cells. Each circle represents one variant and its size corresponds to the number of clones with the same nt sequence. The original sequence is represented by the biggest circle except in panel (c) where it is the circle at the top. Mutation positions are indicated in each branch. For the point mutations, the nt present in each viral population is shown. For the 6 nt deletion found in the viruses serially passaged in C6/36 cells at nt positions 8556/61 (panel <bold>A</bold>
) which was considered as a unique event, (i) the word ‘Del’ indicates that it is present in the nearest viral population, (ii) no modification compared to the original sequence are represented by U, (iii) C and A indicate the presence of the 8566u>c and 8566u>a mutations, respectively.</p>
</caption>
<graphic xlink:href="ppat.1003172.g008"></graphic>
</fig>
<fig id="ppat-1003172-g009" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1003172.g009</object-id>
<label>Figure 9</label>
<caption><title>Mutation emergence kinetic curves.</title>
<p>For (<bold>A</bold>
) the virus serially passaged in Vero cells, (<bold>B</bold>
) serially passaged in C6/36 cells, and (<bold>C</bold>
) alternatively passaged, curves represent the number of mutations detected in complete genome consensus sequences during the passages. These graphics take into account the fixation rate of each mutation corresponding to the presence of mixed viral populations detected when double peaks were observed on the sequencing chromatograms (<bold>Table S4 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
). Numbers are as follows: 1 for a fixed mutation, 0.75 for <, 0.5 for ≈, 0.25 for > and zero when the mutation was not yet detected.</p>
</caption>
<graphic xlink:href="ppat.1003172.g009"></graphic>
</fig>
<table-wrap id="ppat-1003172-t001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1003172.t001</object-id>
<label>Table 1</label>
<caption><title>Characteristics of the mutations detected in the complete genome consensus sequences of the passaged viruses.</title>
</caption>
<alternatives><graphic id="ppat-1003172-t001-1" xlink:href="ppat.1003172.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>
<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>
<col align="center" span="1"></col>
</colgroup>
<thead><tr><td align="left" rowspan="1" colspan="1">Virus</td>
<td align="left" rowspan="1" colspan="1">Passage method</td>
<td align="left" rowspan="1" colspan="1">No. of mutations</td>
<td colspan="4" align="left" rowspan="1">Point mutations</td>
<td align="left" rowspan="1" colspan="1">No. of deletions</td>
<td colspan="3" align="left" rowspan="1">Reversions</td>
<td align="left" rowspan="1" colspan="1">No. of fixed mutations</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">total</td>
<td align="left" rowspan="1" colspan="1">5′UTR</td>
<td align="left" rowspan="1" colspan="1">S</td>
<td align="left" rowspan="1" colspan="1">NS</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">No.</td>
<td align="left" rowspan="1" colspan="1">Freq.</td>
<td align="left" rowspan="1" colspan="1">Rate</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1">WT</td>
<td align="left" rowspan="1" colspan="1">Vero</td>
<td align="left" rowspan="1" colspan="1">10</td>
<td align="left" rowspan="1" colspan="1">10</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">2</td>
<td align="left" rowspan="1" colspan="1">8</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">5</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">C6/36</td>
<td align="left" rowspan="1" colspan="1">3</td>
<td align="left" rowspan="1" colspan="1">2</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">2</td>
<td align="left" rowspan="1" colspan="1">1</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">0</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">Alternate</td>
<td align="left" rowspan="1" colspan="1">6</td>
<td align="left" rowspan="1" colspan="1">6</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">1</td>
<td align="left" rowspan="1" colspan="1">5</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">3</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Φnsp4</td>
<td align="left" rowspan="1" colspan="1">Vero</td>
<td align="left" rowspan="1" colspan="1">8</td>
<td align="left" rowspan="1" colspan="1">8</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">4</td>
<td align="left" rowspan="1" colspan="1">4</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">1</td>
<td align="left" rowspan="1" colspan="1">12.5%</td>
<td align="left" rowspan="1" colspan="1">0.34%</td>
<td align="left" rowspan="1" colspan="1">6</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">C6/36</td>
<td align="left" rowspan="1" colspan="1">6</td>
<td align="left" rowspan="1" colspan="1">4</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">4</td>
<td align="left" rowspan="1" colspan="1">2</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">0%</td>
<td align="left" rowspan="1" colspan="1">0%</td>
<td align="left" rowspan="1" colspan="1">0</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">Alternate</td>
<td align="left" rowspan="1" colspan="1">7</td>
<td align="left" rowspan="1" colspan="1">7</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">1</td>
<td align="left" rowspan="1" colspan="1">6</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">0%</td>
<td align="left" rowspan="1" colspan="1">0%</td>
<td align="left" rowspan="1" colspan="1">0</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Reenc<sup>3</sup>
</td>
<td align="left" rowspan="1" colspan="1">Vero</td>
<td align="left" rowspan="1" colspan="1">16</td>
<td align="left" rowspan="1" colspan="1">16</td>
<td align="left" rowspan="1" colspan="1">2</td>
<td align="left" rowspan="1" colspan="1">7</td>
<td align="left" rowspan="1" colspan="1">7</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">2</td>
<td align="left" rowspan="1" colspan="1">12.5%</td>
<td align="left" rowspan="1" colspan="1">0.23%</td>
<td align="left" rowspan="1" colspan="1">8</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">C6/36</td>
<td align="left" rowspan="1" colspan="1">4</td>
<td align="left" rowspan="1" colspan="1">4</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">4</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">0%</td>
<td align="left" rowspan="1" colspan="1">0%</td>
<td align="left" rowspan="1" colspan="1">1</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">Alternate</td>
<td align="left" rowspan="1" colspan="1">10</td>
<td align="left" rowspan="1" colspan="1">10</td>
<td align="left" rowspan="1" colspan="1">2</td>
<td align="left" rowspan="1" colspan="1">2</td>
<td align="left" rowspan="1" colspan="1">6</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">0%</td>
<td align="left" rowspan="1" colspan="1">0%</td>
<td align="left" rowspan="1" colspan="1">6</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot><fn id="nt101"><label></label>
<p>WT, Φnsp4 and Reenc<sup>3</sup>
represent the WT, Φnsp4 and Φnsp1 Φnsp4 Φenv viruses passaged 50 times using the method described in the column entitled ‘Passage method’.</p>
</fn>
<fn id="nt102"><label></label>
<p>5′UTR designates mutations in the 5′ untranslated region, S synonymous mutations and NS non-synonymous mutations.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<table-wrap id="ppat-1003172-t002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1003172.t002</object-id>
<label>Table 2</label>
<caption><title>Shared mutations detected in complete genome consensus sequences and in viral sub-populations via cloning.</title>
</caption>
<alternatives><graphic id="ppat-1003172-t002-2" xlink:href="ppat.1003172.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>
</colgroup>
<thead><tr><td align="left" rowspan="1" colspan="1">Nt position</td>
<td align="left" rowspan="1" colspan="1">Region</td>
<td align="left" rowspan="1" colspan="1">Nt change</td>
<td align="left" rowspan="1" colspan="1">AA change</td>
<td colspan="2" align="left" rowspan="1">Viruses sharing the mutation when based on:</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">Consensus sequences</td>
<td align="left" rowspan="1" colspan="1">Sub-populations</td>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1">22</td>
<td align="left" rowspan="1" colspan="1">5′UTR</td>
<td align="left" rowspan="1" colspan="1">A→G</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">Reenc<sup>3</sup>
-VERO, Reenc<sup>3</sup>
-ALT</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">66</td>
<td align="left" rowspan="1" colspan="1">5′UTR</td>
<td align="left" rowspan="1" colspan="1">A→G</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">Reenc<sup>3</sup>
-VERO, Reenc<sup>3</sup>
-ALT</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">822</td>
<td align="left" rowspan="1" colspan="1">nsP1</td>
<td align="left" rowspan="1" colspan="1">C→U</td>
<td align="left" rowspan="1" colspan="1">P→L</td>
<td align="left" rowspan="1" colspan="1">Reenc<sup>3</sup>
-VERO, Reenc<sup>3</sup>
-ALT</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">3218</td>
<td align="left" rowspan="1" colspan="1">nsP2</td>
<td align="left" rowspan="1" colspan="1">U→G</td>
<td align="left" rowspan="1" colspan="1">S→A</td>
<td align="left" rowspan="1" colspan="1">Φnsp4-ALT, Reenc<sup>3</sup>
-ALT</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">4167</td>
<td align="left" rowspan="1" colspan="1">nsP3</td>
<td align="left" rowspan="1" colspan="1">G→U/A</td>
<td align="left" rowspan="1" colspan="1">G→V/D</td>
<td align="left" rowspan="1" colspan="1">Φnsp4-ALT</td>
<td align="left" rowspan="1" colspan="1">WT-ALT, Reenc<sup>3</sup>
-C6/36</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">6330</td>
<td align="left" rowspan="1" colspan="1">nsP4</td>
<td align="left" rowspan="1" colspan="1">U→C</td>
<td align="left" rowspan="1" colspan="1">V→A</td>
<td align="left" rowspan="1" colspan="1">Φnsp4-VERO, Reenc<sup>3</sup>
-ALT</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">6426</td>
<td align="left" rowspan="1" colspan="1">nsP4</td>
<td align="left" rowspan="1" colspan="1">C→U</td>
<td align="left" rowspan="1" colspan="1">A→V</td>
<td align="left" rowspan="1" colspan="1">WT-VERO, WT-ALT, Reenc<sup>3</sup>
-ALT</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">8563-8</td>
<td align="left" rowspan="1" colspan="1">E2</td>
<td align="left" rowspan="1" colspan="1">Del-GUCUAU</td>
<td align="left" rowspan="1" colspan="1">Del-VY</td>
<td align="left" rowspan="1" colspan="1">WT-C6/36, Φnsp4-C6/36</td>
<td align="left" rowspan="1" colspan="1">Reenc<sup>3</sup>
-C6/36</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">8566</td>
<td align="left" rowspan="1" colspan="1">E2</td>
<td align="left" rowspan="1" colspan="1">U→C/A</td>
<td align="left" rowspan="1" colspan="1">Y→H/N</td>
<td align="left" rowspan="1" colspan="1">WT-C6/36, Φnsp4-C6/36</td>
<td align="left" rowspan="1" colspan="1">Reenc<sup>3</sup>
-C6/36</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">8731</td>
<td align="left" rowspan="1" colspan="1">E2</td>
<td align="left" rowspan="1" colspan="1">U→C/A</td>
<td align="left" rowspan="1" colspan="1">W→R</td>
<td align="left" rowspan="1" colspan="1">WT-VERO, WT-ALT</td>
<td align="left" rowspan="1" colspan="1">WT-C6/36, Φnsp4-ALT</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">8831</td>
<td align="left" rowspan="1" colspan="1">E2</td>
<td align="left" rowspan="1" colspan="1">U→C</td>
<td align="left" rowspan="1" colspan="1">M→T</td>
<td align="left" rowspan="1" colspan="1">Reenc<sup>3</sup>
-C6/36</td>
<td align="left" rowspan="1" colspan="1">WT-C6/36</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">8836</td>
<td align="left" rowspan="1" colspan="1">E2</td>
<td align="left" rowspan="1" colspan="1">C→U</td>
<td align="left" rowspan="1" colspan="1">H→>Y</td>
<td align="left" rowspan="1" colspan="1">WT-VERO</td>
<td align="left" rowspan="1" colspan="1">WT-C6/36</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">10,509</td>
<td align="left" rowspan="1" colspan="1">E1</td>
<td align="left" rowspan="1" colspan="1">C→U</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">Reenc<sup>3</sup>
-VERO</td>
<td align="left" rowspan="1" colspan="1">Reenc<sup>3</sup>
-C6/36</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">10,542</td>
<td align="left" rowspan="1" colspan="1">E1</td>
<td align="left" rowspan="1" colspan="1">C→U</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">Reenc<sup>3</sup>
-VERO</td>
<td align="left" rowspan="1" colspan="1">Reenc<sup>3</sup>
-ALT</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">10,670</td>
<td align="left" rowspan="1" colspan="1">E1</td>
<td align="left" rowspan="1" colspan="1">U→C</td>
<td align="left" rowspan="1" colspan="1">V→A</td>
<td align="left" rowspan="1" colspan="1">WT-C6/36, WT-ALT, Φnsp4-C6/36, Φnsp4-ALT, Reenc<sup>3</sup>
-C6/36, Reenc<sup>3</sup>
-ALT</td>
<td align="left" rowspan="1" colspan="1">Reenc<sup>3</sup>
-VERO</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">3641</td>
<td align="left" rowspan="1" colspan="1">nsP2</td>
<td align="left" rowspan="1" colspan="1">G→A/C</td>
<td align="left" rowspan="1" colspan="1">V→L</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">Φnsp4-C6/36, Reenc<sup>3</sup>
-ALT</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">3897</td>
<td align="left" rowspan="1" colspan="1">nsP2</td>
<td align="left" rowspan="1" colspan="1">G→A</td>
<td align="left" rowspan="1" colspan="1">C→>Y</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">WT-C6/36, WT-ALT</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">8394</td>
<td align="left" rowspan="1" colspan="1">E3</td>
<td align="left" rowspan="1" colspan="1">C→U</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">WT-C6/36, Reenc<sup>3</sup>
-ALT</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">8558</td>
<td align="left" rowspan="1" colspan="1">E2</td>
<td align="left" rowspan="1" colspan="1">U→C</td>
<td align="left" rowspan="1" colspan="1">F→S</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">WT-C6/36, Φnsp4-C6/36</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">8573</td>
<td align="left" rowspan="1" colspan="1">E2</td>
<td align="left" rowspan="1" colspan="1">C→U</td>
<td align="left" rowspan="1" colspan="1">A→V</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">Φnsp4-C6/36, Φnsp4-ALT, Reenc<sup>3</sup>
-C6/36</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">8750</td>
<td align="left" rowspan="1" colspan="1">E2</td>
<td align="left" rowspan="1" colspan="1">U→C/A</td>
<td align="left" rowspan="1" colspan="1">M→K/T</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">WT-C6/36, Φnsp4-C6/36</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">8830</td>
<td align="left" rowspan="1" colspan="1">E2</td>
<td align="left" rowspan="1" colspan="1">A→G</td>
<td align="left" rowspan="1" colspan="1">M→V</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">ΦNSP4-C6/36, Reenc<sup>3</sup>
-ALT</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">8898</td>
<td align="left" rowspan="1" colspan="1">E2</td>
<td align="left" rowspan="1" colspan="1">G→U</td>
<td align="left" rowspan="1" colspan="1">R→S</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">ΦNSP4-VERO, ΦNSP4-C6/36</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">8783</td>
<td align="left" rowspan="1" colspan="1">E2</td>
<td align="left" rowspan="1" colspan="1">C→U</td>
<td align="left" rowspan="1" colspan="1">A→V</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">ΦNSP4-ALT, Reenc<sup>3</sup>
-ALT</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot><fn id="nt103"><label></label>
<p>WT, Φnsp4 and Reenc<sup>3</sup>
represent the passaged WT, Φnsp4 and Φnsp1 Φnsp4 Φenv viruses. Viruses are represented by their name followed by ‘VERO’ for serial passages in Vero, ‘C6/36’ for serial passages in C6/36, and ‘ALT’ for alternate passages. All the mutations detected in complete genome consensus sequences during the passages were used (<bold>Table S4 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
) as well as those used for the analysis of intra-population genetic diversity (<xref ref-type="fig" rid="ppat-1003172-g008"><bold>Figure 8</bold>
</xref>
<bold> and Figure S6 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
).</p>
</fn>
<fn id="nt104"><label></label>
<p>5′UTR designates the 5′ untranslated region.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>Of note, a six nucleotide deletion at positions 8562/72 which was found in the three viruses serially passaged in C6/36 (twice in the complete genome consensus sequences and once using cloning methods), and there was a rapid increase in frequency of a 9 nt deletion at positions 4139/47 in Φnsp4 (from 0% to 60% of the overall population in only 13 passages). This indicates that acquisition of deletion mutations constituted an important mechanism for adaptation to C6/36 cells (<bold>Table S4 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
). In addition, complex evolutionary pathways were observed during passages in C6/36 cells: some sub-populations acquired mutations – such as 8573c>u, 8831u>c and 8563_8568del – which then disappeared partially or completely during subsequent passages (<xref ref-type="fig" rid="ppat-1003172-g008"><bold>Figure 8</bold>
</xref>
and <bold>Figure S6 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
).</p>
<p>The re-encoded viruses, particularly the Φnsp1 Φnsp4 Φenv virus, were also associated with specific evolutionary patterns. Serial passages in Vero cells resulted in a higher number of mutations in the Φnsp1 Φnsp4 Φenv virus (16 versus 10 and 8 for the WT and Φnsp4 viruses, respectively) and an increased proportion of synonymous mutations in both re-encoded viruses (50% versus 20% for the WT virus) (<xref ref-type="table" rid="ppat-1003172-t001"><bold>Table 1</bold>
</xref>
). During alternate passages for the Φnsp1 Φnsp4 Φenv virus, there was also an early appearance and fixation of mutations (<xref ref-type="fig" rid="ppat-1003172-g009"><bold>Figure 9</bold>
</xref>
), whilst at the 50<sup>th</sup>
passage, all alternately passaged viruses displayed a similar total number of mutations (<xref ref-type="table" rid="ppat-1003172-t001"><bold>Table 1</bold>
</xref>
). These specific evolutionary patterns indicate that despite the adaptation to culture conditions shared by the WT and re-encoded viruses, a proportion of the mutations that emerged during serial passage in Vero cells appeared in response to codon re-encoding. The fact that the Φnsp1 Φnsp4 Φenv virus serially passaged in Vero cells showed significantly increased replicative fitness in both Vero cells and C6/36 cells, reinforces this conclusion.</p>
<p>To detect the mutations specific to codon re-encoding, we first hypothesized that adaptation to culture conditions and response to codon re-re-encoding may have occurred sequentially. However, the fact that some mutations associated with codon re-encoding appeared very early (<italic>e.g.</italic>
, 22a>g, 66au and 10542 c>u, see below) suggested that both phenomena occurred concurrently.</p>
<p>When those convergent mutations shared only between re-encoded viruses were analysed (<xref ref-type="table" rid="ppat-1003172-t002"><bold>Table 2</bold>
</xref>
), the early emergence of two identical mutations located in the 5′untranslated region (5′UTR) (22a>g and 66 a>g) was identified for the Φnsp1 Φnsp4 Φenv viruses serially passaged either in Vero cells or alternately in both cells. This suggested that these mutations induced major fitness changes in response to codon re-encoding (both mutations were already detected at the third passage for both viruses). In addition, two convergent non-synonymous mutations were identified within the re-encoded nsP1 (822c>u) and nsP4 (6330u>c) regions, suggesting that non-synonymous mutations could arise in response to codon re-encoding (<xref ref-type="table" rid="ppat-1003172-t003"><bold>Table 3</bold>
</xref>
) in the case of viruses passaged in Vero cells or alternately in both cells. No convergent mutation was identified in C6/36 passaged viruses.</p>
<table-wrap id="ppat-1003172-t003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1003172.t003</object-id>
<label>Table 3</label>
<caption><title>Mutations that emerged in response to codon re-encoding and detected using complete genome consensus sequences.</title>
</caption>
<alternatives><graphic id="ppat-1003172-t003-3" xlink:href="ppat.1003172.t003"></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">Nucleotide position</td>
<td align="left" rowspan="1" colspan="1">Region</td>
<td align="left" rowspan="1" colspan="1">Nucleotide change</td>
<td align="left" rowspan="1" colspan="1">AA change</td>
<td align="left" rowspan="1" colspan="1">Reversion</td>
<td align="left" rowspan="1" colspan="1">In a re-encoding region</td>
<td align="left" rowspan="1" colspan="1">Virus(es) harboring the mutation</td>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1">22</td>
<td align="left" rowspan="1" colspan="1">5′UTR</td>
<td align="left" rowspan="1" colspan="1">A→G</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">Reenc<sup>3</sup>
-VERO, Reenc<sup>3</sup>
-ALT</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">66</td>
<td align="left" rowspan="1" colspan="1">5′UTR</td>
<td align="left" rowspan="1" colspan="1">A→G</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">Reenc<sup>3</sup>
-VERO, Reenc<sup>3</sup>
-ALT</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">742</td>
<td align="left" rowspan="1" colspan="1">Nsp1</td>
<td align="left" rowspan="1" colspan="1">A→U</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">+</td>
<td align="left" rowspan="1" colspan="1">Reenc<sup>3</sup>
-ALT</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">822</td>
<td align="left" rowspan="1" colspan="1">Nsp1</td>
<td align="left" rowspan="1" colspan="1">C→U</td>
<td align="left" rowspan="1" colspan="1">P→L</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">+</td>
<td align="left" rowspan="1" colspan="1">Reenc<sup>3</sup>
-VERO, Reenc<sup>3</sup>
-ALT</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">6330</td>
<td align="left" rowspan="1" colspan="1">Nsp4</td>
<td align="left" rowspan="1" colspan="1">U→C</td>
<td align="left" rowspan="1" colspan="1">V→A</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">+</td>
<td align="left" rowspan="1" colspan="1">Φnsp4-VERO, Reenc<sup>3</sup>
-ALT</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">6670</td>
<td align="left" rowspan="1" colspan="1">Nsp4</td>
<td align="left" rowspan="1" colspan="1">C→A</td>
<td align="left" rowspan="1" colspan="1">N→K</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">+</td>
<td align="left" rowspan="1" colspan="1">Φnsp4-C6/36</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">6717</td>
<td align="left" rowspan="1" colspan="1">Nsp4</td>
<td align="left" rowspan="1" colspan="1">A→C</td>
<td align="left" rowspan="1" colspan="1">E→A</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">+</td>
<td align="left" rowspan="1" colspan="1">Reenc<sup>3</sup>
-VERO</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">6761</td>
<td align="left" rowspan="1" colspan="1">Nsp4</td>
<td align="left" rowspan="1" colspan="1">A→G</td>
<td align="left" rowspan="1" colspan="1">T→A</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">+</td>
<td align="left" rowspan="1" colspan="1">Reenc<sup>3</sup>
-VERO</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">6971</td>
<td align="left" rowspan="1" colspan="1">Nsp4</td>
<td align="left" rowspan="1" colspan="1">U→C</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">+</td>
<td align="left" rowspan="1" colspan="1">+</td>
<td align="left" rowspan="1" colspan="1">Φnsp4-VERO</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">9591</td>
<td align="left" rowspan="1" colspan="1">E2</td>
<td align="left" rowspan="1" colspan="1">U→G</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">+</td>
<td align="left" rowspan="1" colspan="1">Reenc<sup>3</sup>
-ALT</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">9855</td>
<td align="left" rowspan="1" colspan="1">6K</td>
<td align="left" rowspan="1" colspan="1">A→G</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">+</td>
<td align="left" rowspan="1" colspan="1">Reenc<sup>3</sup>
-VERO</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">10419</td>
<td align="left" rowspan="1" colspan="1">E1</td>
<td align="left" rowspan="1" colspan="1">U→C</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">+</td>
<td align="left" rowspan="1" colspan="1">+</td>
<td align="left" rowspan="1" colspan="1">Reenc<sup>3</sup>
-VERO</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">10509</td>
<td align="left" rowspan="1" colspan="1">E1</td>
<td align="left" rowspan="1" colspan="1">C→U</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1">+</td>
<td align="left" rowspan="1" colspan="1">+</td>
<td align="left" rowspan="1" colspan="1">Reenc<sup>3</sup>
-VERO</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">10542</td>
<td align="left" rowspan="1" colspan="1">E1</td>
<td align="left" rowspan="1" colspan="1">C→U</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">+</td>
<td align="left" rowspan="1" colspan="1">Reenc<sup>3</sup>
-VERO</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">10896</td>
<td align="left" rowspan="1" colspan="1">E1</td>
<td align="left" rowspan="1" colspan="1">C→U</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">+</td>
<td align="left" rowspan="1" colspan="1">Reenc<sup>3</sup>
-VERO</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot><fn id="nt105"><label></label>
<p>Φnsp4 and Reenc<sup>3</sup>
represent the passaged Φnsp4 and Φnsp1 Φnsp4 Φenv viruses. Viruses are represented by their name followed by ‘VERO’ for serial passages in Vero, ‘C6/36’ for serial passages in C6/36, and ‘ALT’ for alternate passages.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>Eleven mutations were identified specifically in the re-encoded regions of genetically modified viruses (<xref ref-type="table" rid="ppat-1003172-t003"><bold>Table 3</bold>
</xref>
). These included (i) 3 reversions: 2 for Φnsp1 Φnsp4 Φenv virus and 1 for Φnsp4 virus passaged in Vero cells; (ii) 5 additional synonymous mutations observed for Φnsp1 Φnsp4 Φenv virus passaged in Vero cells or alternately in both cells; (iii) 3 non-synonymous mutations, of which 2 appeared in Φnsp1 Φnsp4 Φenv virus passaged in Vero cells and 1 in Φnsp4 virus passaged in C6/36 cells.</p>
<p>In summary, despite the substantial specific replicative fitness gains observed for the re-encoded viruses, there was a marked absence of reversion mutations for viruses passaged in C6/36 cells and alternately in vertebrate/invertebrate cells, and only a low reversion rate (<0.4%) in Vero cells (<xref ref-type="table" rid="ppat-1003172-t001"><bold>Table 1</bold>
</xref>
). Rather than re-introducing their original codons through synonymous back mutations, the re-encoded viruses largely acquired non-synonymous mutations (61%) or novel synonymous mutations (22%).</p>
<p>In C6/36 cells, only one mutation candidate in response to codon re-encoding (6670c>a) was identified for the Φnsp4 virus despite the significant replicative fitness changes observed (<xref ref-type="table" rid="ppat-1003172-t003"><bold>Table 3</bold>
</xref>
). By contrast, in Vero cells 12 mutation candidates could be identified, including a cluster of six synonymous mutations in the Φenv re-encoded region. This specific response to codon re-encoding probably reflects the greater decrease in replicative fitness observed for the Φenv virus in comparison to other viruses with one region re-encoded (see above).</p>
<p>The distribution of mutations along the genome observed for re-encoded viruses partly matched that of the WT virus (a probable consequence of the adaptation to the same culture conditions), but also that of genetic variability observed in CHIKV genomes retrieved from GenBank (<xref ref-type="fig" rid="ppat-1003172-g010"><bold>Figure 10</bold>
</xref>
). Passaged viruses exhibited the same highly variable regions and some of these mutational hot spots (<italic>e.g.</italic>
, the hypervariable C-terminal region of the nsP3 <xref rid="ppat.1003172-Strauss1" ref-type="bibr">[17]</xref>
and the E2 region) were also highly variable amongst GenBank CHIKV sequences. These regions were associated with high levels of amino acid diversity in passaged viruses (<bold>Table S4 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
) and other CHIKVs (<xref ref-type="fig" rid="ppat-1003172-g010"><bold>Figure 10</bold>
</xref>
). Conversely, regions of low variability amongst the GenBank CHIKV genomes were associated with a reduction in synonymous variability, particularly within nsP1, nsP2 and E1, and were also invariable in all passaged viruses (except the nsP1 region). Nevertheless, some mutational hot spots found in all passaged viruses were relatively invariable amongst GenBank CHIKV genomes, notably the X-domain of the nsP3 <xref rid="ppat.1003172-Malet1" ref-type="bibr">[41]</xref>
, and specific regions in nsP1, E2/E1 and the capsid binding RNA region <xref rid="ppat.1003172-Strauss1" ref-type="bibr">[17]</xref>
(<xref ref-type="fig" rid="ppat-1003172-g010"><bold>Figure 10</bold>
</xref>
and <bold>Table S4 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
).</p>
<fig id="ppat-1003172-g010" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1003172.g010</object-id>
<label>Figure 10</label>
<caption><title>Mutation distributions.</title>
<p>Patterns of mutation distribution in 132 CHIKV extracted from GenBank (<bold>A–C</bold>
) and of the passaged viruses (<bold>D–G</bold>
), all of them aligned with the CHIKV complete genome sequence (<bold>H</bold>
). Using the alignment of the CHIKVs from GenBank, the variability at (<bold>A</bold>
) 1<sup>st</sup>
+2<sup>nd</sup>
+3<sup>rd</sup>
(<bold>B</bold>
) 1<sup>st</sup>
+2<sup>nd</sup>
codon and (<bold>C</bold>
) 3<sup>rd</sup>
codon positions were determined and represented using 150, 100 and 50 nt sliding windows, respectively. (<bold>A</bold>
): 1 represents the mean value from all sites; (<bold>B–C</bold>
): values represent the number of variable nt positions in the sliding window interval; mean +/− standard deviations were also represented). Using complete genome consensus sequences from (<bold>D</bold>
) all passaged viruses, (<bold>E</bold>
) viruses serially passaged in Vero cells, (<bold>F</bold>
) viruses passaged in C6/36 cells or (<bold>G</bold>
) viruses alternatively passaged, the sum of mutations detected at each nt position (deletions were considered as one unique mutation) was determined and represented using 100 nt sliding window. All the mutations are represented in light blue and mutations from re-encoded viruses were superimposed in dark blue. Shaded regions indicate the three re-encoded regions. From left to right, brackets indicate the X-domain and the hypervariable C-terminal region of the nsP3 and the capsid binding RNA region. The grey arrow indicates the E1-V226A mutation.</p>
</caption>
<graphic xlink:href="ppat.1003172.g010"></graphic>
</fig>
<p>Finally, the mutation at position 10670 (E1-V226A) was observed in all viruses passaged in C6/36 cells and alternately in Vero and C6/36 cells but was fixed only once (during serial passages in C6/36 of Φnsp1 Φnsp4 Φenv virus). This confirms that the A226V mutation, a strong inducer of adaptation to <italic>A. albopictus</italic>
(<xref ref-type="fig" rid="ppat-1003172-g010"><bold>Figure 10</bold>
</xref>
), does not provide a significant fitness advantage for replication in the C6/36 <italic>A. albopictus</italic>
cell culture system <xref rid="ppat.1003172-Schuffenecker1" ref-type="bibr">[19]</xref>
, <xref rid="ppat.1003172-Tsetsarkin1" ref-type="bibr">[22]</xref>
.</p>
</sec>
</sec>
</sec>
<sec id="s3"><title>Discussion</title>
<p>We have evaluated the effect on replicative fitness and cytopathogenicity of large-scale re-encoding of CHIKV, a re-emerging Old World pathogenic arbovirus. The generation of attenuated viruses by large-scale re-encoding represents an exciting and potentially important route to vaccine development, and also to understanding the basis of the evolution of viral pathogenicity. Site-directed re-encoding, associated with no modification of amino acid sequences, alleviates the likelihood of novel phenotypic properties, allows us to modulate fitness by altering the length of the codon replacement interval, but additionally provides benefits to the generic development of live attenuated vaccines, including reduced costs and single dose induction of long-term immunity <xref rid="ppat.1003172-Coffin1" ref-type="bibr">[42]</xref>
.</p>
<p>A key result was the observation that our random re-encoding method decreased the replicative fitness of CHIKV in both primate and arthropod cells. The diminution of CHIKV replicative fitness correlated directly with the degree of re-encoding. As reported in previous studies with unrelated viruses, we found that during one replicative cycle in mosquito cells, codon re-encoding profoundly reduced the infectious titre of released virus whilst the number of viral particles remained stable <xref rid="ppat.1003172-Burns2" ref-type="bibr">[12]</xref>
, <xref rid="ppat.1003172-Mueller2" ref-type="bibr">[14]</xref>
. This implies that the maturation process (i.e. the formation of ribonucleoproteins and their insertion into plasma membranes that contain HA) could be at fault when viruses are re-encoded. In contrast, in primate cells, this decline in infectivity of the viral particles was associated with the reduced generation of viral RNA and proteins probably due to a compromised replication complex. Because they can be identified at different stages of the CHIKV replication cycle, these results imply that the observed decrease in replicative fitness is probably the consequence of several independent re-encoding induced events. However, it is important to note that (i) alphaviruses produce very different kinds of infection in mosquito and primate cells (i.e. persistent infection in mosquito cells and cytolytic infection in vertebrate cells), and (ii) there are differences in host cell response and innate immunity between mosquito and primate cells, and which could also explain these observed differences in the different cell lines <xref rid="ppat.1003172-Strauss1" ref-type="bibr">[17]</xref>
, <xref rid="ppat.1003172-Li1" ref-type="bibr">[43]</xref>
–<xref rid="ppat.1003172-Wikan1" ref-type="bibr">[45]</xref>
.</p>
<p>There is mounting evidence that synonymous mutations in viral genomes may have major fitness effects and not only in the small number of cis-acting elements described previously <xref rid="ppat.1003172-Cuevas1" ref-type="bibr">[46]</xref>
. In the current study, six re-encoded viruses were produced of which the most re-encoded virus modified in three regions that encode different proteins (together, 882 synonymous mutations were introduced spanning 4,212 nt). In support of previous studies which demonstrated that re-encoded poliovirus and influenza A viruses are attenuated <xref rid="ppat.1003172-Mueller1" ref-type="bibr">[10]</xref>
–<xref rid="ppat.1003172-Mueller2" ref-type="bibr">[14]</xref>
, our observations of a reduction in replicative fitness strongly suggest that a proportion of synonymous mutations are not neutral in RNA viruses. Indeed, it is likely that some synonymous mutations were positively selected during the passaging process, reinforcing the idea that synonymous sites are central to viral fitness. In conclusion, it is likely that synonymous mutations can be either neutral, beneficial or deleterious as is the case for non-synonymous mutations.</p>
<p>Evolutionary patterns at synonymous sites could be shaped by genome-wide mutational processes, such as G+C%, codon usage bias and dinucleotide frequency <xref rid="ppat.1003172-Jenkins1" ref-type="bibr">[8]</xref>
, <xref rid="ppat.1003172-Chen1" ref-type="bibr">[47]</xref>
–<xref rid="ppat.1003172-Zhou1" ref-type="bibr">[49]</xref>
. These global constraints, which theoretically produce a subset of viable genomes, were assessed by previous studies of codon re-encoding in poliovirus, influenza A virus and bacterial virus T7 which applied specific modification of codon usage bias, codon pair bias or CpG/UpA frequencies <xref rid="ppat.1003172-Mueller1" ref-type="bibr">[10]</xref>
–<xref rid="ppat.1003172-Bull1" ref-type="bibr">[15]</xref>
. Using a large-scale random re-encoding method, which only slightly modified these global properties, we still observed replicative fitness reductions in both primate and arthropod cells. Our results suggest that local constraints may also provide significant selection pressure on synonymous sites in RNA viruses, for example by disrupting RNA secondary structures. Since numerous functional secondary structures are present in coding regions of RNA viruses, and hence include synonymous sites (with notable examples in poliovirus <xref rid="ppat.1003172-Goodfellow1" ref-type="bibr">[50]</xref>
, tick-borne encephalitis virus <xref rid="ppat.1003172-Tuplin1" ref-type="bibr">[51]</xref>
, alphaviruses <xref rid="ppat.1003172-Strauss1" ref-type="bibr">[17]</xref>
, <xref rid="ppat.1003172-Ou1" ref-type="bibr">[52]</xref>
and HIV-1 <xref rid="ppat.1003172-Watts1" ref-type="bibr">[53]</xref>
), it is likely that similar structures are common in CHIKV. Recently, it was demonstrated that a similar re-encoding strategy applied to the noncapsid regions of the poliovirus resulted in the identification of two novel functional RNA elements <xref rid="ppat.1003172-Song1" ref-type="bibr">[54]</xref>
. The concept of large-scale random re-encoding, as described here, is also supported by the report of the negative impact of random single synonymous mutations (which did not modify the genetic characteristics of the genome) on viral replicative fitness <xref rid="ppat.1003172-Cuevas1" ref-type="bibr">[46]</xref>
.</p>
<p>Finally, our results and those of previous studies of re-encoded viruses <xref rid="ppat.1003172-Mueller1" ref-type="bibr">[10]</xref>
–<xref rid="ppat.1003172-Mueller2" ref-type="bibr">[14]</xref>
suggest that the reduction of viral replicative fitness is driven by a variety of factors. First, the nature of the virus studied is an important parameter: we found that introducing up to 882 random synonymous mutations clearly affected the replicative fitness of the CHIKV, whilst two previous studies demonstrated that comparable random re-encoding methods applied to the capsid precursor (P1) region of the poliovirus did not significantly affect replicative fitness (934 <xref rid="ppat.1003172-Mueller2" ref-type="bibr">[14]</xref>
and 153 <xref rid="ppat.1003172-Burns1" ref-type="bibr">[11]</xref>
synonymous substitutions were introduced, respectively). The location of the re-encoded region constitutes the second factor of importance: re-encoding in the E2/E1 region resulted in a greater loss of fitness than in other genomic regions. The analysis of complete wild type CHIKV genomes revealed naturally low levels of synonymous diversity in this re-encoded region (<xref ref-type="fig" rid="ppat-1003172-g010"><bold>Figure 10</bold>
</xref>
) indicating that this region is subject to specific local evolutionary constraints which in part explain the significant impact of re-encoding in this region. The re-encoding method applied is obviously an additional important parameter: previous studies with polioviruses showed that reduction of replicative fitness was strongly dependent on the method used to re-encode the genome <xref rid="ppat.1003172-Burns1" ref-type="bibr">[11]</xref>
–<xref rid="ppat.1003172-Mueller2" ref-type="bibr">[14]</xref>
.</p>
<p>The average impact of one mutation is clearly likely to be less important in random re-encoding than in specific approaches <xref rid="ppat.1003172-Burns1" ref-type="bibr">[11]</xref>
–<xref rid="ppat.1003172-Mueller2" ref-type="bibr">[14]</xref>
. This suggests that random large-scale re-encoding could be advantageous in several aspects when designing future vaccine candidates, namely: <italic>(i)</italic>
reversion to wild-type should be intrinsically more difficult, given the high number of mutations produced; <italic>(ii)</italic>
since in our experiments the reduction of replicative fitness decreased with the degree of re-encoding, the method opens the door to finely tuning fitness reduction through modulation of the length of re-encoded regions and the number of synonymous mutations introduced; <italic>(iii)</italic>
the use of a combination of several re-encoded regions located throughout the viral genome may prevent complete phenotypic reversion due to recombination between WT and re-encoded viruses: large scale sequence modification may render recombination intrinsically more difficult, and in the case of recombination, the part of the genome representing the re-encoded strain would likely still carry some mutations associated with fitness reduction.</p>
<p>Taken together, these observations suggest that, following large-scale random re-encoding, recovery of the original replicative fitness should require a large number of reversion mutations. Consequently these re-encoded viruses should be very stable <xref rid="ppat.1003172-Burns2" ref-type="bibr">[12]</xref>
, <xref rid="ppat.1003172-Mueller2" ref-type="bibr">[14]</xref>
. To test this hypothesis and to study the constraints that shape CHIKV codon usage, we passaged the wild type and two re-encoded CHIKVs <italic>in cellulo</italic>
. It is commonly stated that arboviruses are subject to strong evolutionary trade-offs, such that mutations that are favoured in one host are deleterious in another, and that this imposes constraints on viral evolution <xref rid="ppat.1003172-Coffey1" ref-type="bibr">[16]</xref>
, <xref rid="ppat.1003172-Taylor1" ref-type="bibr">[37]</xref>
, <xref rid="ppat.1003172-Weaver2" ref-type="bibr">[39]</xref>
. This implies that a strain that has been adapted to a mammalian cell line should have a reduced replicative fitness in mosquito cells, and <italic>vice versa</italic>
. Consequently, we could have initiated our experiments <italic>(i)</italic>
with a strain previously adapted to a given cell line, with the disadvantage of introducing an obvious bias in other cell lines (<italic>i.e.</italic>
, follow a strategy inspired by Bull <italic>et al.</italic>
who used a pre-adapted bacterial virus T7 <xref rid="ppat.1003172-Burns2" ref-type="bibr">[12]</xref>
, <xref rid="ppat.1003172-Bull1" ref-type="bibr">[15]</xref>
), or <italic>(ii)</italic>
with a clinical strain, <italic>i.e.</italic>
isolated from previous alternate passages in mosquitoes and humans, with the disadvantage of possibly observing a simultaneous response to codon re-encoding and adaptation to culture conditions (as performed by Burns <italic>et al.</italic>
who used a non-adapted poliovirus <xref rid="ppat.1003172-Burns2" ref-type="bibr">[12]</xref>
, <xref rid="ppat.1003172-Bull1" ref-type="bibr">[15]</xref>
). We chose the latter, based on the observation that no criteria existed for defining the adaptation period in the case of arboviruses. This choice was retrospectively justified by the observation that the follow-up of replicative fitness and molecular evolution could not distinguish the criteria necessary for differentiating both phenomena.</p>
<p>A key observation of our study was that few reversion mutations occurred, despite specific replicative fitness enhancements in response to codon re-encoding. This suggests that the effectiveness of large-scale re-encoding methods results from the accumulation of slightly deleterious mutations that push the virus into a fitness valley, and that there are multiple opportunities through diverse mutational pathways, sometimes in genomic regions other than those we re-encoded, in which these viruses can partially restore their fitness. For example, amongst the mutation candidates emerging in response to codon re-encoding, two were in the 5′ UTR. That these mutations were fixed so rapidly is strongly suggestive of their selectively beneficial effect <xref rid="ppat.1003172-Wichman1" ref-type="bibr">[55]</xref>
. Moreover, it was notable that although we only made modifications to synonymous sites, some of the mutations observed were non-synonymous. Therefore, the evolution of these viruses in response to codon re-encoding was largely compensatory in nature and very few mutations were the result of reversion. Moreover, even with the specific fitness enhancements in response to codon re-encoding observed, the most re-encoded virus failed to reach fitness levels (infectious titres) equivalent to those observed at the first passage of the WT virus. Interestingly, these fitness improvements were not always accompanied by greater numbers of mutations or specific molecular modifications. For example, during alternate passage, the most re-encoded virus succeeded in increasing its fitness without accumulating more mutations than other viruses (<italic>i.e.</italic>
, WT and the less re-encoded viruses).</p>
<p>During serial passage of the re-encoded viruses, we observed that the response to codon re-encoding and adaptation to culture conditions occurred simultaneously. However, the high levels of observed convergent evolution between the WT virus and the re-encoded viruses indicates that selection arising from codon re-encoding was likely weaker than that for adaptation to culture conditions, and/or that the beneficial mutations to restore the cost of re-encoding were less likely to arise. Therefore, this indirect insight into the difficulty of reversing the effects of re-encoding further highlights the stability of these re-encoded viruses.</p>
<p>Our experiments also confirm that mutations acquired in one host can be deleterious in a different host type (serial passages in primate cells increased viral replicative fitness in primate cells, whilst serial passages in mosquito cells decreased viral fitness in primate cells) and, with the exception of the most re-encoded virus, that alternate passages seriously (i) limit replicative fitness enhancement, and (ii) delay the appearance of the mutations. It is noteworthy that replicative fitness in mosquito cells remained globally unchanged following serial passage in mosquito cells, in Vero cells or in both. Moreover, as described previously with dengue virus <xref rid="ppat.1003172-Vasilakis1" ref-type="bibr">[38]</xref>
, replication in mosquito cells appears to act as a brake on viral evolution; in our case, very few mutations were detected after 50 passages and only one was fixed, suggesting that the majority of the emerging mutants have a deleterious effect on viral fitness. In addition, we observed weaker selection pressure in these cells during competition experiments. The fact that the C6/36 cell line was selected more than 30 years ago for its capacity to replicate CHIKVs and dengue viruses <xref rid="ppat.1003172-Igarashi1" ref-type="bibr">[56]</xref>
could explain this weaker selection pressure observed. Conversely, the rapid adaptation of CHIKV to <italic>A. albopictus</italic>
was accompanied by multiple appearances of the E1-A226V mutation <xref rid="ppat.1003172-Schuffenecker1" ref-type="bibr">[19]</xref>
, <xref rid="ppat.1003172-Tsetsarkin1" ref-type="bibr">[22]</xref>
–<xref rid="ppat.1003172-deLamballerie1" ref-type="bibr">[23]</xref>
and the appearance of amino-acid deletions in the nsP3 and E2 genes.</p>
<p>In conclusion, this study demonstrates that random codon re-encoding significantly decreases the replicative fitness of CHIKV. Although all these results are important and encouraging, they cannot be easily extended to RNA viruses producing chronic infections. Thus, studies in animal models are obviously needed to evaluate the potential of these new generation attenuation methods for producing vaccine candidates. However, this approach could assist in the development of future RNA virus vaccines, including those for arboviruses. Introducing a large number of slightly deleterious synonymous mutations reduced the replicative fitness of CHIKV by orders of magnitude in both primate and arthropod cells. This strategy resulted in limited reversion and recovery of fitness after intensive serial subculture of the viruses, and is likely to reduce the risk of complete phenotypic reversion if recombination with wild type virus occurs. Our results encourage us that such modified viruses would find it difficult to return to their natural arboviral cycle in the real world. Furthermore, the decrease of the replicative fitness correlated with the extent of re-encoding, an observation that may be advantageous in the development of future strategies to modulate viral attenuation.</p>
</sec>
<sec sec-type="materials|methods" id="s4"><title>Materials and Methods</title>
<p>Cell culture protocols for virus stock production and virus titration (TCID50, plaque and haemagglutination assays) are detailed in <bold>Protocol S1</bold>
.</p>
<sec id="s4a"><title>In s<italic>ilico</italic>
re-encoding method</title>
<p>Three regions of the CHIKV genome were re-encoded using a computer program that randomly attributed nucleotide codons based on their corresponding amino acid sequence: for example, the amino acid valine was randomly replaced by GTT, GTC, GTA or GTG. To minimize the influence of rare codons in primate cell lines, the number and the position of such rare codons in primate genomes <xref rid="ppat.1003172-Zhang1" ref-type="bibr">[57]</xref>
(i.e. CGU, CGC, CGA, CGG, UCG, CCG, GCG, ACG) were not modified. In addition, unique restriction sites were conserved by correcting synonymous mutations at some sites. The location of the re-encoded cassettes, first based on the availability of unique restriction sites was adjusted to avoid overlap with known RNA secondary structures <xref rid="ppat.1003172-Strauss1" ref-type="bibr">[17]</xref>
, <xref rid="ppat.1003172-Ou1" ref-type="bibr">[52]</xref>
. Finally, three cassettes of 1302, 1410 and 1500 bases and located in the nsP1, nsP4 and E2/E1 regions, respectively, were designed using this method (<bold><xref ref-type="supplementary-material" rid="ppat.1003172.s003">Text S2</xref>
</bold>
).</p>
</sec>
<sec id="s4b"><title>Construction of CHIKV infectious clones (ICs)</title>
<p>We modified a previously described IC of the LR2006 strain of CHIKV <xref rid="ppat.1003172-Tsetsarkin2" ref-type="bibr">[58]</xref>
(GenBank accession EU224268) and all the re-encoded regions were synthesized (GenScript) and then inserted into ICs as described in <bold>Protocol S1</bold>
. Using a combination of re-encoded regions, six re-encoded ICs were generated: Φnsp1, Φnsp4 and Φenv with one re-encoded region; Φnsp1 Φnsp4, and Φnsp4 Φenv with two re-encoded regions and Φnsp1 Φnsp4 Φenv with three re-encoded regions (<bold>Figure S1 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
).</p>
</sec>
<sec id="s4c"><title>Real time RT-PCR assays</title>
<p>A fragment of 179 nt located in the nsP2 region (nucleotide position 2631 to 2809) was used to detect the genomic RNA (plus strand) of all the CHIKVs (universal assay), re-encoded or not. Another fragment of 168 nt located in the nsP4 region (nucleotide position 6804 to 6971) was used to analyze cell supernatants from competition experiments: two sets of primers and probes allowed us to specifically detect either the viruses re-encoded in the nsP4 region or the viruses without modification in the same region. Primer and probe sequences and the real time PCR protocol are detailed in <bold>Table S6 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
and <bold>Protocol S1</bold>
, respectively.</p>
</sec>
<sec id="s4d"><title>Replication kinetics</title>
<p>The replicative fitness of each virus was determined using the results of replication kinetics studies, performed in triplicate in Vero, HEK293 or C6/36 cells. For comparison of the seven viruses from the seven ICs (the WT virus and the 6 re-encoded viruses), one experiment was performed with all the viruses. Virus stock (see <bold>Protocol S1</bold>
) or ICs were used to infect or transfect cells respectively. For the evaluation of replicative fitness of the passaged viruses, we performed one experiment for each virus (WT, Φnsp4 and Φnsp1 Φnsp4 Φenv viruses) with the first passage in Vero and the 12<sup>th</sup>
, 25<sup>th</sup>
, 37<sup>th</sup>
and 50<sup>th</sup>
passages for each passage regimen (13 supernatants tested in triplicate). For the single cycle replication kinetics (<xref ref-type="fig" rid="ppat-1003172-g001"><bold>Figure 1</bold>
</xref>
), an estimated MOI of 5 was used to infect a 75 cm<sup>2</sup>
culture flask of confluent Vero, C6/36 or HEK293 cells. Cells were washed twice (HBSS) 30 minutes after the infection and 20 ml of medium was added. 1 ml of cell supernatant was sampled just before the washes and at 2, 8, 14, 20 and 28 hours pi. For the replication kinetics with low estimated MOI (<xref ref-type="fig" rid="ppat-1003172-g003"><bold>Figure 3</bold>
</xref>
) and the evaluation of the replicative fitness of the passaged viruses (<bold>Figure S3, S3, S5 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
), an estimated MOI of 0.01 was used to infect a 25 cm<sup>2</sup>
culture flask of confluent Vero or C6/36 cells. Cells were washed twice (HBSS) 2 hours after infection and 8 ml of medium was added. 1 ml of cell supernatant was sampled after the washes (T0) and at 24, 48 and 72 hours pi. For the replication kinetics using infectious DNA clones (<xref ref-type="fig" rid="ppat-1003172-g004"><bold>Figure 4</bold>
</xref>
), a 75 cm<sup>2</sup>
culture flask of subconfluent HEK293 cells was transfected with the ICs using Lipofectamine 2000 (Invitrogen) according to the manufacturer's instructions. Cells were washed twice (HBSS) 4 hours after the transfection and 20 ml of medium was added. 1 ml of cell supernatant was sampled after the washes (T0) and at 16, 24 and 48 hours pi.</p>
<p>All the sampled cell supernatants were clarified by centrifugation, aliquoted and stored at −80°C. They were then analysed using a TCID50 assay and a real time RT-PCR assay (not performed systematically, see figure legends). Nucleic acids were extracted from clarified cell supernatants using the EZ1 Virus Mini Kit v2 on the EZ1 Biorobot (both from Qiagen).</p>
</sec>
<sec id="s4e"><title>Virus competition experiments</title>
<p>WT virus was grown in competition with one of four re-encoded viruses (Φnsp4, Φnsp1 Φnsp4, Φnsp1 Φenv or Φnsp1 Φnsp4 Φenv) using five different PFU ratios (WT/re-encoded virus 1/99, 20/80, 50/50, 80/20, 99/1). A global estimated MOI of 0.5 was used for the first inoculation. For each experiment, a 25 cm<sup>2</sup>
flask culture of confluent cells was infected for 2 hours, washed (HBSS) and then incubated for 48 h after the addition of 7 ml of medium. Viruses from each experiment were then passaged nine times as follows: a 25 cm<sup>2</sup>
flask culture of confluent cells was infected for 2 hours with the purified culture supernatant (centrifugation), washed (HBSS) and then incubated for 48 h after the addition of 7 ml of medium. At each passage, the estimated MOI was bottlenecked at approximately 0.5. After each infection, nucleic acids were extracted from the clarified culture supernatant using the EZ1 Virus Mini Kit v2 on the EZ1 Biorobot (both from Qiagen). Using two specific real time RT-PCR assays targeting the ΦnsP4 region (see above), the amount of each virus was assessed and the ratio of the two values (WT/re-encoded) was calculated.</p>
</sec>
<sec id="s4f"><title>Quantification of intracellular RNA and viral proteins</title>
<p>A global estimated MOI of 5 was used to infect confluent 12 well-plates of HEK293 cells with virus stock (see <bold>Protocol S1</bold>
). Cells were washed once (HBSS) 30 minutes after the infection and 2 ml of media was added. At 8 hours pi, the absence of cytopathic effect was checked, culture supernatants were discarded, and cells were washed once (HBSS). All experiments were performed in triplicate. For Western blot analysis and intracellular viral RNA quantification, total RNA and protein isolation was performed using the same well with the Nucleospin RNA/protein kit according to the manufacturer's instructions (Macherey-Nagel). Protein extracts were resolved on 10% polyacrylamide gels containing SDS and transferred to PVDF membrane. Anti-Nsp1/2 rabbit pAb (see <bold>Protocol S1</bold>
), anti-actin C-2 mAb (Santa Cruz Biotechnology) and the corresponding HRP-conjugated secondary antibody were used. Protein bands were revealed using Immobilon (Millipore) followed by exposure of blot to radiographic film. Real time RT-PCR assay (see above) was performed to assess viral intracellular RNA (mRNA actin was used as a normalizer to account for differences in cells number and/or quality of extracted RNA as described previously <xref rid="ppat.1003172-Baronti1" ref-type="bibr">[59]</xref>
). For the quantification of viral proteins by ELISA, cells were mechanically harvested using a cell scraper, resuspended in 800 µL of PBS, vortexed and disrupted by sonication (30 seconds at 20 KHz, Misonix Sonicator XL). Pre-treated CHIKV-specific immune human serum was used to detect viral proteins. ELISA protocol is detailed in <bold>Protocol S1</bold>
.</p>
</sec>
<sec id="s4g"><title>Experimental passage of viruses <italic>in cellulo</italic>
</title>
<p>The WT and two re-encoded viruses (Φnsp4 and Φnsp1 Φnsp4 Φenv) were passaged 50 times following three regimens: serial passages in Vero or C6/36 cells and alternate passages between Vero and C6/36. For each passage, a 25 cm<sup>2</sup>
culture flask of confluent cells was infected for 2 hours with the diluted clarified cell supernatant, washed (HBSS) and incubated for 48 hours after the addition of 7 ml of medium. Cell supernatant was then harvested, clarified by centrifugation, aliquoted and stored at −80°C. For each passage, the estimated MOI was bottlenecked at approximately 0.1. To avoid contamination, virus passages were performed in three phases: serial passages of WT and Φnsp4 viruses, alternate passages of the same viruses and passages of the Φnsp1 Φnsp4 Φenv virus. All the viruses passaged at the same time were manipulated sequentially and in different laminar flow cabinets.</p>
</sec>
<sec id="s4h"><title>Consensus sequencing</title>
<p>Whole genome nucleotide sequences (excluding the first 18 nucleotides of the 5′UTR and the 22 nucleotides upstream of the polyA tail) were determined for all the 50<sup>th</sup>
passage viruses (nine viruses in total). The timing of emergence of each mutation found in the 50<sup>th</sup>
passage was then determined by sequencing with appropriate primer pairs for the 6<sup>th</sup>
, 12<sup>th</sup>
, 18<sup>th</sup>
, 25<sup>th</sup>
, 31<sup>st</sup>
, 37<sup>th</sup>
and 43<sup>rd</sup>
passages. To avoid contamination by PCR products and plasmids, we utilized a molecular biology laboratory that is specifically designed for clinical diagnosis using molecular techniques, and which includes separate laboratories dedicated to perform nucleic acid extraction, PCR/RT mix, RNA/cDNA manipulations and PCR products/plasmids manipulations. In addition, each step from extraction to sequencing or cloning of PCR products (see below) was performed in separate experiments for each passaged virus. Nucleic acids were extracted from the purified culture supernatant using the EZ1 Virus Mini Kit v2 on the EZ1 Biorobot (both from Qiagen). A set of specific primer pairs (<bold>Table S5 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
) was used to generate amplicons with the Access RT-PCR System (Promega) according to the manufacturer's instructions. PCR products were then purified and sequenced using forward and reverse primers with the BigDye Terminator v3.1 Cycle Sequencing Kit on an ABI Prism 31310X Genetic Analyser sequencer (both from Applied Biosystems). Analysis of sequencing chromatogram and combination of sequences was performed using the Sequencher 4.9 software (Gene Codes Corporation).</p>
</sec>
<sec id="s4i"><title>Intra-population genetic diversity</title>
<p>Three genomic regions were chosen for analysis: one of 619 nt overlapping the nsP2 and nsP3 regions (positions 3601 to 4220), one of 619 nt overlapping the E3 and E2 regions (positions 8351 to 8970), and one of 740 nt located in the E1 region (positions 10140 to 10880). The 1<sup>st</sup>
, 12<sup>th</sup>
, 25<sup>th</sup>
, 37<sup>th</sup>
, 43<sup>rd</sup>
and 50<sup>th</sup>
passages of each passaged virus were chosen to assess the intra-population genetic diversity in the E3/E2 region, while the nsP2/nsP3 and E1 regions were analyzed in the passaged Φnsp1 Φnsp4 Φenv viruses (using the same passage numbers as above). In addition, the 12<sup>th</sup>
, 25<sup>th</sup>
, 37<sup>th</sup>
, 43<sup>rd</sup>
and 50<sup>th</sup>
passages of the Φnsp4 virus serially passaged in C6/36 were also analyzed in the nsP2/nsP3 region to study the emergence of the 4139_4147del mutation. The AccuScript <italic>PfuUltra</italic>
II RT-PCR Kit (Agilent) was used according to the manufacturer's instructions to generate amplicons from extracted nucleic acids (see above). PCR products were cloned after purification into the StrataClone PCR Cloning Vector and transformed into competent cells (StrataClone PCR Cloning Kit; Agilent). A plasmid extraction was performed from bacterial colonies with the correct insert which had been previously cultured and plasmid DNA was automatically sequenced with a T7 primer (GATC Biotech). Finally, approximately 20 clones were sequenced for each PCR product analyzed (mean: 19.2+/−4.2; range from 14 to 44), resulting in a total of 315, 933 and 227 sequenced clones for the nsP2/nsP3, E3/E2 and E1 regions, respectively.</p>
</sec>
<sec id="s4j"><title>Evolutionary analysis</title>
<p>Complete viral genome consensus sequences were manually constructed and aligned. Base ambiguity symbols were used to represent all the mixed viral populations when double peaks were observed on the sequencing chromatograms (as represented ≈, > and < in <bold>Table S4 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
). The 132 CHIKV complete genome nucleotide sequences already available in GenBank (<bold><xref ref-type="supplementary-material" rid="ppat.1003172.s003">Text S2</xref>
</bold>
), along with the outgroup O'Nyong-Nyong virus (ONNV) strain SG650 (GenBank accession AF079456) sequence, were extracted from GenBank/NCBI. The two ORFs of these viruses were manually extracted and concatenated using the Bioedit v7.0.9 program <xref rid="ppat.1003172-Hall1" ref-type="bibr">[60]</xref>
and aligned with ClustalW <xref rid="ppat.1003172-Thompson1" ref-type="bibr">[61]</xref>
according to the amino acid sequence. Ambiguously aligned regions were removed manually. This alignment was used to estimate the variability at each nt position using the Mega5 software <xref rid="ppat.1003172-Tamura1" ref-type="bibr">[62]</xref>
–<xref rid="ppat.1003172-Tamura2" ref-type="bibr">[63]</xref>
.</p>
<p>All the sequences obtained to assess the intra-population genetic diversity were edited using the Sequencher 4.9 software (Gene Codes Corporation). Using the Mega5 software <xref rid="ppat.1003172-Tamura1" ref-type="bibr">[62]</xref>
, all the sequences from the same region were aligned and each mutation which was found only once in the comparative alignment (i.e. singletons) was removed to ensure that no mutations introduced during the RT-PCR were included in the analysis: 22, 49 and 15 singletons were removed from the nsP2/nsP3, E3/E2 and E1 alignments, respectively. Minimum spanning trees of clonal CHIKV data were constructed using the TCS1.21 software <xref rid="ppat.1003172-Clement1" ref-type="bibr">[64]</xref>
with a 99% connection limit. For each virus and primer pair, we used the alignment of all the sequences (see above) after removing singletons. Each deletion (6 or 9 nt) was considered to be a unique event irrespective of length. When a mutation occurred in a deleted region (detected in another clone), it was assigned the position of the corresponding deletion.</p>
</sec>
<sec id="s4k"><title>Estimation of the relative fitness effect and statistical analysis</title>
<p>To study the evolution of replicative fitness during the serial subculture of viruses we performed replicative kinetics studies at the 1<sup>st</sup>
, 12<sup>th</sup>
, 25<sup>th</sup>
, 37<sup>th</sup>
and 50<sup>th</sup>
passage of each virus in Vero and C6/36 cells (<bold>Figure S3, S3, S5 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
). Data from these kinetic studies were analyzed using two different methods.</p>
<p>We first estimated the relative fitness effect values (<xref ref-type="fig" rid="ppat-1003172-g007"><bold>Figure 7</bold>
</xref>
): we measured the viral growth rate at 24 hours pi, based on TCID50 values, to calculate the relative fitness effect as described previously <xref rid="ppat.1003172-Peris1" ref-type="bibr">[65]</xref>
. This measure has to be performed during the viral exponential growth phase. Because most of the infectious titre values obtained with the Φnsp1 Φnsp4 Φenv virus at 24 hours pi were close to the detection threshold of our TCID50 assay, we chose the 24 hours pi values for all the viruses, even though some of the values of the WT and Φnsp4 viruses reached a plateau before 24 hours (<bold>Figure S3, S4, S5 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
). Titres at time <italic>t<sub>0</sub>
</italic>
(all the values were under the detection threshold of our TCID50 assay and then considered as zero) and at <italic>t<sub>1</sub>
</italic>
(24 hours pi) were used to calculate the growth rate (<italic>r</italic>
) as the increase in log-titre per 24 hours. Relative fitness (<italic>W</italic>
) was defined as the growth rate ratio and the relative fitness effect as <italic>s = W−1</italic>
. The relative fitness for each experiment <italic>i</italic>
was calculated as <inline-formula><inline-graphic xlink:href="ppat.1003172.e001.jpg"></inline-graphic>
</inline-formula>
where <inline-formula><inline-graphic xlink:href="ppat.1003172.e002.jpg"></inline-graphic>
</inline-formula>
is the average of three determinations for the first passage of the virus in the same cells. We then performed a global analysis of TCID50 values at 24, 48 and 72 hours pi by performing two-way repeated-measures ANOVA <xref rid="ppat.1003172-Richardson1" ref-type="bibr">[66]</xref>
and tukey's HDS post-hoc comparisons (<bold>Table S3 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
). This method was used before to analyze similar results <xref rid="ppat.1003172-Vasilakis1" ref-type="bibr">[38]</xref>
. A tukey's HDS post-hoc test was used to compare, once in Vero cells and once in C6/36 cells, the replicative fitness of the first passage for each virus (WT, Φnsp4 or Φnsp1 Φnsp4 Φenv virus) with that of the corresponding passaged virus. Comparison of the significance of the results from both methods showed that they gave very similar results (67/72 were concordant; <bold>Table S3 in <xref ref-type="supplementary-material" rid="ppat.1003172.s002">Text S1</xref>
</bold>
). Infectious titres, viral RNA yields and relative fitness effect values were compared using a Student's <italic>t</italic>
test. For all statistical tests used, all <italic>p</italic>
values below 0.05 were considered significant.</p>
</sec>
</sec>
<sec sec-type="supplementary-material" id="s5"><title>Supporting Information</title>
<supplementary-material content-type="local-data" id="ppat.1003172.s001"><label>Protocol S1</label>
<caption><p>Cells and antibodies, construction of CHIKV infectious clones, plasmid transfection/virus stock production, plaque assay, tissue culture infectious dose 50 (TCID50) assay, real time RT-PCR assay, haemagglutination assay and quantification of intracellular viral proteins (ELISA method).</p>
<p>(PDF)</p>
</caption>
<media xlink:href="ppat.1003172.s001.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="ppat.1003172.s002"><label>Text S1</label>
<caption><p>Figure S1: Schematic representation of the CHIKV re-encoded viruses. Figure S2: Relationship between either infectious titres or viral RNA yields and the number of synonymous mutations in re-encoded region(s). Figure S3: Replication curves with WT passaged viruses. Figure S4: Replication curves with Φnsp4 passaged viruses. Figure S5: Replication curves with Φnsp1 Φnsp4 Φenv passaged viruses. Figure S6: Intra-population genetic diversity of CHIKV revealed using minimum spanning trees. Figure S7: Schematic representation of the CHIKV infectious clones (IC). Table S1: Genetic characteristics of the coding regions (concatenated ORFs) of the re-encoded viruses, the WT virus and 132 other CHIKVs extracted from GenBank. Table S2: Summary of Single cycle replication kinetics values at 8 and 14 hours pi. Table S3: Replicative fitness modifications observed during the passage experiments: comparison of the results from both analysis methods. Table S4: Mutations detected in the CHIKV consensus sequences during experimental passage. Table S5: Primers used for the sequencing of CHIKVs. Table S6: Primers and probes used for the real time RT-PCR assays.</p>
<p>(PDF)</p>
</caption>
<media xlink:href="ppat.1003172.s002.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="ppat.1003172.s003"><label>Text S2</label>
<caption><p>The re-encoded sequences, the list of the 132 CHIKV sequences extracted from GenBank, the nucleotide sequence of the synthetic RNA transcript used as standard for the universal real time RT-PCR assay and the amino-acid sequence of the recombinant protein used to immunized rabbits.</p>
<p>(PDF)</p>
</caption>
<media xlink:href="ppat.1003172.s003.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
</sec>
</body>
<back><ack><p>We thank Karine Barthelemy and Reine De Chesse from the UMR190-Emergence des Pathologies Virales (Marseille, France) for their technical contribution; and Andrew Kitchen, Etienne Simon-Loriere and Israel Pagán from the Center for Infectious Disease (The Pennsylvania State University, USA) for their advice on evolutionary analysis; and Bruno Coutard from the AFMB laboratory (Architecture et Fonction des Macromolécules Biologiques, UMR 6098, Marseille, France) for providing the recombinant protein (fusion between the C-terminal region of the nsP2 and the N-terminal region of the nsP3).</p>
</ack>
<ref-list><title>References</title>
<ref id="ppat.1003172-Guzman1"><label>1</label>
<mixed-citation publication-type="journal"><name><surname>Guzman</surname>
<given-names>MG</given-names>
</name>
, <name><surname>Halstead</surname>
<given-names>SB</given-names>
</name>
, <name><surname>Artsob</surname>
<given-names>H</given-names>
</name>
, <name><surname>Buchy</surname>
<given-names>P</given-names>
</name>
, <name><surname>Farrar</surname>
<given-names>J</given-names>
</name>
, <etal>et al</etal>
(<year>2010</year>
) <article-title>Dengue: a continuing global threat</article-title>
. <source>Nat Rev Microbiol</source>
<volume>8</volume>
: <fpage>S7</fpage>
–<lpage>16</lpage>
.<pub-id pub-id-type="pmid">21079655</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Lanciotti1"><label>2</label>
<mixed-citation publication-type="journal"><name><surname>Lanciotti</surname>
<given-names>RS</given-names>
</name>
, <name><surname>Roehrig</surname>
<given-names>JT</given-names>
</name>
, <name><surname>Deubel</surname>
<given-names>V</given-names>
</name>
, <name><surname>Smith</surname>
<given-names>J</given-names>
</name>
, <name><surname>Parker</surname>
<given-names>M</given-names>
</name>
, <etal>et al</etal>
(<year>1999</year>
) <article-title>Origin of the West Nile virus responsible for an outbreak of encephalitis in the northeastern United States</article-title>
. <source>Science</source>
<volume>286</volume>
: <fpage>2333</fpage>
–<lpage>2337</lpage>
.<pub-id pub-id-type="pmid">10600742</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Staples1"><label>3</label>
<mixed-citation publication-type="journal"><name><surname>Staples</surname>
<given-names>JE</given-names>
</name>
, <name><surname>Breiman</surname>
<given-names>RF</given-names>
</name>
, <name><surname>Powers</surname>
<given-names>AM</given-names>
</name>
(<year>2009</year>
) <article-title>Chikungunya fever: an epidemiological review of a re-emerging infectious disease</article-title>
. <source>Clin Infect Dis</source>
<volume>49</volume>
: <fpage>942</fpage>
–<lpage>948</lpage>
.<pub-id pub-id-type="pmid">19663604</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Boshra1"><label>4</label>
<mixed-citation publication-type="journal"><name><surname>Boshra</surname>
<given-names>H</given-names>
</name>
, <name><surname>Lorenzo</surname>
<given-names>G</given-names>
</name>
, <name><surname>Busquets</surname>
<given-names>N</given-names>
</name>
, <name><surname>Brun</surname>
<given-names>A</given-names>
</name>
(<year>2011</year>
) <article-title>Rift valley fever: recent insights into pathogenesis and prevention</article-title>
. <source>J Virol</source>
<volume>85</volume>
: <fpage>6098</fpage>
–<lpage>6105</lpage>
.<pub-id pub-id-type="pmid">21450816</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Randolph1"><label>5</label>
<mixed-citation publication-type="journal"><name><surname>Randolph</surname>
<given-names>SE</given-names>
</name>
, <name><surname>Rogers</surname>
<given-names>DJ</given-names>
</name>
(<year>2010</year>
) <article-title>The arrival, establishment and spread of exotic diseases: patterns and predictions</article-title>
. <source>Nat Rev Microbiol</source>
<volume>8</volume>
: <fpage>361</fpage>
–<lpage>371</lpage>
.<pub-id pub-id-type="pmid">20372156</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Weaver1"><label>6</label>
<mixed-citation publication-type="journal"><name><surname>Weaver</surname>
<given-names>SC</given-names>
</name>
, <name><surname>Reisen</surname>
<given-names>WK</given-names>
</name>
(<year>2010</year>
) <article-title>Present and future arboviral threats</article-title>
. <source>Antiviral Res</source>
<volume>85</volume>
: <fpage>328</fpage>
–<lpage>345</lpage>
.<pub-id pub-id-type="pmid">19857523</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Hershberg1"><label>7</label>
<mixed-citation publication-type="journal"><name><surname>Hershberg</surname>
<given-names>R</given-names>
</name>
, <name><surname>Petrov</surname>
<given-names>DA</given-names>
</name>
(<year>2008</year>
) <article-title>Selection on codon bias</article-title>
. <source>Annu Rev Genet</source>
<volume>42</volume>
: <fpage>287</fpage>
–<lpage>299</lpage>
.<pub-id pub-id-type="pmid">18983258</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Jenkins1"><label>8</label>
<mixed-citation publication-type="journal"><name><surname>Jenkins</surname>
<given-names>GM</given-names>
</name>
, <name><surname>Holmes</surname>
<given-names>EC</given-names>
</name>
(<year>2003</year>
) <article-title>The extent of codon usage bias in human RNA viruses and its evolutionary origin</article-title>
. <source>Virus Res</source>
<volume>92</volume>
: <fpage>1</fpage>
–<lpage>7</lpage>
.<pub-id pub-id-type="pmid">12606071</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Sanchez1"><label>9</label>
<mixed-citation publication-type="journal"><name><surname>Sanchez</surname>
<given-names>G</given-names>
</name>
, <name><surname>Bosch</surname>
<given-names>A</given-names>
</name>
, <name><surname>Pinto</surname>
<given-names>RM</given-names>
</name>
(<year>2003</year>
) <article-title>Genome variability and capsid structural constraints of hepatitis a virus</article-title>
. <source>J Virol</source>
<volume>77</volume>
: <fpage>452</fpage>
–<lpage>459</lpage>
.<pub-id pub-id-type="pmid">12477850</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Mueller1"><label>10</label>
<mixed-citation publication-type="journal"><name><surname>Mueller</surname>
<given-names>S</given-names>
</name>
, <name><surname>Coleman</surname>
<given-names>JR</given-names>
</name>
, <name><surname>Papamichail</surname>
<given-names>D</given-names>
</name>
, <name><surname>Ward</surname>
<given-names>CB</given-names>
</name>
, <name><surname>Nimnual</surname>
<given-names>A</given-names>
</name>
, <etal>et al</etal>
(<year>2010</year>
) <article-title>Live attenuated influenza virus vaccines by computer-aided rational design</article-title>
. <source>Nat Biotechnol</source>
<volume>28</volume>
: <fpage>723</fpage>
–<lpage>726</lpage>
.<pub-id pub-id-type="pmid">20543832</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Burns1"><label>11</label>
<mixed-citation publication-type="journal"><name><surname>Burns</surname>
<given-names>CC</given-names>
</name>
, <name><surname>Campagnoli</surname>
<given-names>R</given-names>
</name>
, <name><surname>Shaw</surname>
<given-names>J</given-names>
</name>
, <name><surname>Vincent</surname>
<given-names>A</given-names>
</name>
, <name><surname>Jorba</surname>
<given-names>J</given-names>
</name>
, <etal>et al</etal>
(<year>2009</year>
) <article-title>Genetic inactivation of poliovirus infectivity by increasing the frequencies of CpG and UpA dinucleotides within and across synonymous capsid region codons</article-title>
. <source>J Virol</source>
<volume>83</volume>
: <fpage>9957</fpage>
–<lpage>9969</lpage>
.<pub-id pub-id-type="pmid">19605476</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Burns2"><label>12</label>
<mixed-citation publication-type="journal"><name><surname>Burns</surname>
<given-names>CC</given-names>
</name>
, <name><surname>Shaw</surname>
<given-names>J</given-names>
</name>
, <name><surname>Campagnoli</surname>
<given-names>R</given-names>
</name>
, <name><surname>Jorba</surname>
<given-names>J</given-names>
</name>
, <name><surname>Vincent</surname>
<given-names>A</given-names>
</name>
, <etal>et al</etal>
(<year>2006</year>
) <article-title>Modulation of poliovirus replicative fitness in HeLa cells by deoptimization of synonymous codon usage in the capsid region</article-title>
. <source>J Virol</source>
<volume>80</volume>
: <fpage>3259</fpage>
–<lpage>3272</lpage>
.<pub-id pub-id-type="pmid">16537593</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Coleman1"><label>13</label>
<mixed-citation publication-type="journal"><name><surname>Coleman</surname>
<given-names>JR</given-names>
</name>
, <name><surname>Papamichail</surname>
<given-names>D</given-names>
</name>
, <name><surname>Skiena</surname>
<given-names>S</given-names>
</name>
, <name><surname>Futcher</surname>
<given-names>B</given-names>
</name>
, <name><surname>Wimmer</surname>
<given-names>E</given-names>
</name>
, <etal>et al</etal>
(<year>2008</year>
) <article-title>Virus attenuation by genome-scale changes in codon pair bias</article-title>
. <source>Science</source>
<volume>320</volume>
: <fpage>1784</fpage>
–<lpage>1787</lpage>
.<pub-id pub-id-type="pmid">18583614</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Mueller2"><label>14</label>
<mixed-citation publication-type="journal"><name><surname>Mueller</surname>
<given-names>S</given-names>
</name>
, <name><surname>Papamichail</surname>
<given-names>D</given-names>
</name>
, <name><surname>Coleman</surname>
<given-names>JR</given-names>
</name>
, <name><surname>Skiena</surname>
<given-names>S</given-names>
</name>
, <name><surname>Wimmer</surname>
<given-names>E</given-names>
</name>
(<year>2006</year>
) <article-title>Reduction of the rate of poliovirus protein synthesis through large-scale codon deoptimization causes attenuation of viral virulence by lowering specific infectivity</article-title>
. <source>J Virol</source>
<volume>80</volume>
: <fpage>9687</fpage>
–<lpage>9696</lpage>
.<pub-id pub-id-type="pmid">16973573</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Bull1"><label>15</label>
<mixed-citation publication-type="journal"><name><surname>Bull</surname>
<given-names>JJ</given-names>
</name>
, <name><surname>Molineux</surname>
<given-names>IJ</given-names>
</name>
, <name><surname>Wilke</surname>
<given-names>CO</given-names>
</name>
(<year>2012</year>
) <article-title>Slow fitness recovery in a codon-modified viral genome</article-title>
. <source>Mol Biol Evol</source>
<volume>29</volume>
: <fpage>2997</fpage>
–<lpage>3004</lpage>
.<pub-id pub-id-type="pmid">22532576</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Coffey1"><label>16</label>
<mixed-citation publication-type="journal"><name><surname>Coffey</surname>
<given-names>LL</given-names>
</name>
, <name><surname>Vasilakis</surname>
<given-names>N</given-names>
</name>
, <name><surname>Brault</surname>
<given-names>AC</given-names>
</name>
, <name><surname>Powers</surname>
<given-names>AM</given-names>
</name>
, <name><surname>Tripet</surname>
<given-names>F</given-names>
</name>
, <etal>et al</etal>
(<year>2008</year>
) <article-title>Arbovirus evolution in vivo is constrained by host alternation</article-title>
. <source>Proc Natl Acad Sci U S A</source>
<volume>105</volume>
: <fpage>6970</fpage>
–<lpage>6975</lpage>
.<pub-id pub-id-type="pmid">18458341</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Strauss1"><label>17</label>
<mixed-citation publication-type="journal"><name><surname>Strauss</surname>
<given-names>JH</given-names>
</name>
, <name><surname>Strauss</surname>
<given-names>EG</given-names>
</name>
(<year>1994</year>
) <article-title>The alphaviruses: gene expression, replication, and evolution</article-title>
. <source>Microbiol Rev</source>
<volume>58</volume>
: <fpage>491</fpage>
–<lpage>562</lpage>
.<pub-id pub-id-type="pmid">7968923</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Robinson1"><label>18</label>
<mixed-citation publication-type="journal"><name><surname>Robinson</surname>
<given-names>MC</given-names>
</name>
(<year>1955</year>
) <article-title>An epidemic of virus disease in Southern Province, Tanganyika Territory, in 1952–53. I. Clinical features</article-title>
. <source>Trans R Soc Trop Med Hyg</source>
<volume>49</volume>
: <fpage>28</fpage>
–<lpage>32</lpage>
.<pub-id pub-id-type="pmid">14373834</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Schuffenecker1"><label>19</label>
<mixed-citation publication-type="journal"><name><surname>Schuffenecker</surname>
<given-names>I</given-names>
</name>
, <name><surname>Iteman</surname>
<given-names>I</given-names>
</name>
, <name><surname>Michault</surname>
<given-names>A</given-names>
</name>
, <name><surname>Murri</surname>
<given-names>S</given-names>
</name>
, <name><surname>Frangeul</surname>
<given-names>L</given-names>
</name>
, <etal>et al</etal>
(<year>2006</year>
) <article-title>Genome microevolution of chikungunya viruses causing the Indian Ocean outbreak</article-title>
. <source>PLoS Med</source>
<volume>3</volume>
: <fpage>e263</fpage>
.<pub-id pub-id-type="pmid">16700631</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Pialoux1"><label>20</label>
<mixed-citation publication-type="journal"><name><surname>Pialoux</surname>
<given-names>G</given-names>
</name>
, <name><surname>Gauzere</surname>
<given-names>BA</given-names>
</name>
, <name><surname>Jaureguiberry</surname>
<given-names>S</given-names>
</name>
, <name><surname>Strobel</surname>
<given-names>M</given-names>
</name>
(<year>2007</year>
) <article-title>Chikungunya, an epidemic arbovirosis</article-title>
. <source>Lancet Infect Dis</source>
<volume>7</volume>
: <fpage>319</fpage>
–<lpage>327</lpage>
.<pub-id pub-id-type="pmid">17448935</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Volk1"><label>21</label>
<mixed-citation publication-type="journal"><name><surname>Volk</surname>
<given-names>SM</given-names>
</name>
, <name><surname>Chen</surname>
<given-names>R</given-names>
</name>
, <name><surname>Tsetsarkin</surname>
<given-names>KA</given-names>
</name>
, <name><surname>Adams</surname>
<given-names>AP</given-names>
</name>
, <name><surname>Garcia</surname>
<given-names>TI</given-names>
</name>
, <etal>et al</etal>
(<year>2010</year>
) <article-title>Genome-scale phylogenetic analyses of chikungunya virus reveal independent emergences of recent epidemics and various evolutionary rates</article-title>
. <source>J Virol</source>
<volume>84</volume>
: <fpage>6497</fpage>
–<lpage>6504</lpage>
.<pub-id pub-id-type="pmid">20410280</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Tsetsarkin1"><label>22</label>
<mixed-citation publication-type="journal"><name><surname>Tsetsarkin</surname>
<given-names>KA</given-names>
</name>
, <name><surname>Vanlandingham</surname>
<given-names>DL</given-names>
</name>
, <name><surname>McGee</surname>
<given-names>CE</given-names>
</name>
, <name><surname>Higgs</surname>
<given-names>S</given-names>
</name>
(<year>2007</year>
) <article-title>A single mutation in chikungunya virus affects vector specificity and epidemic potential</article-title>
. <source>PLoS Pathog</source>
<volume>3</volume>
: <fpage>e201</fpage>
.<pub-id pub-id-type="pmid">18069894</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-deLamballerie1"><label>23</label>
<mixed-citation publication-type="journal"><name><surname>de Lamballerie</surname>
<given-names>X</given-names>
</name>
, <name><surname>Leroy</surname>
<given-names>E</given-names>
</name>
, <name><surname>Charrel</surname>
<given-names>RN</given-names>
</name>
, <name><surname>Ttsetsarkin</surname>
<given-names>K</given-names>
</name>
, <name><surname>Higgs</surname>
<given-names>S</given-names>
</name>
, <etal>et al</etal>
(<year>2008</year>
) <article-title>Chikungunya virus adapts to tiger mosquito via evolutionary convergence: a sign of things to come?</article-title>
<source>Virol J</source>
<volume>5</volume>
: <fpage>33</fpage>
.<pub-id pub-id-type="pmid">18304328</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Gerardin1"><label>24</label>
<mixed-citation publication-type="journal"><name><surname>Gerardin</surname>
<given-names>P</given-names>
</name>
, <name><surname>Guernier</surname>
<given-names>V</given-names>
</name>
, <name><surname>Perrau</surname>
<given-names>J</given-names>
</name>
, <name><surname>Fianu</surname>
<given-names>A</given-names>
</name>
, <name><surname>Le Roux</surname>
<given-names>K</given-names>
</name>
, <etal>et al</etal>
(<year>2008</year>
) <article-title>Estimating Chikungunya prevalence in La Reunion Island outbreak by serosurveys: two methods for two critical times of the epidemic</article-title>
. <source>BMC Infect Dis</source>
<volume>8</volume>
: <fpage>99</fpage>
.<pub-id pub-id-type="pmid">18662384</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Economopoulou1"><label>25</label>
<mixed-citation publication-type="journal"><name><surname>Economopoulou</surname>
<given-names>A</given-names>
</name>
, <name><surname>Dominguez</surname>
<given-names>M</given-names>
</name>
, <name><surname>Helynck</surname>
<given-names>B</given-names>
</name>
, <name><surname>Sissoko</surname>
<given-names>D</given-names>
</name>
, <name><surname>Wichmann</surname>
<given-names>O</given-names>
</name>
, <etal>et al</etal>
(<year>2009</year>
) <article-title>Atypical Chikungunya virus infections: clinical manifestations, mortality and risk factors for severe disease during the 2005–2006 outbreak on Reunion</article-title>
. <source>Epidemiol Infect</source>
<volume>137</volume>
: <fpage>534</fpage>
–<lpage>541</lpage>
.<pub-id pub-id-type="pmid">18694529</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Gerardin2"><label>26</label>
<mixed-citation publication-type="journal"><name><surname>Gerardin</surname>
<given-names>P</given-names>
</name>
, <name><surname>Barau</surname>
<given-names>G</given-names>
</name>
, <name><surname>Michault</surname>
<given-names>A</given-names>
</name>
, <name><surname>Bintner</surname>
<given-names>M</given-names>
</name>
, <name><surname>Randrianaivo</surname>
<given-names>H</given-names>
</name>
, <etal>et al</etal>
(<year>2008</year>
) <article-title>Multidisciplinary prospective study of mother-to-child chikungunya virus infections on the island of La Reunion</article-title>
. <source>PLoS Med</source>
<volume>5</volume>
: <fpage>e60</fpage>
.<pub-id pub-id-type="pmid">18351797</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Gould1"><label>27</label>
<mixed-citation publication-type="journal"><name><surname>Gould</surname>
<given-names>EA</given-names>
</name>
, <name><surname>Gallian</surname>
<given-names>P</given-names>
</name>
, <name><surname>De Lamballerie</surname>
<given-names>X</given-names>
</name>
, <name><surname>Charrel</surname>
<given-names>RN</given-names>
</name>
(<year>2010</year>
) <article-title>First cases of autochthonous dengue fever and chikungunya fever in France: from bad dream to reality!</article-title>
. <source>Clin Microbiol Infect</source>
<volume>16</volume>
: <fpage>1702</fpage>
–<lpage>1704</lpage>
.<pub-id pub-id-type="pmid">21040155</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Rezza1"><label>28</label>
<mixed-citation publication-type="journal"><name><surname>Rezza</surname>
<given-names>G</given-names>
</name>
, <name><surname>Nicoletti</surname>
<given-names>L</given-names>
</name>
, <name><surname>Angelini</surname>
<given-names>R</given-names>
</name>
, <name><surname>Romi</surname>
<given-names>R</given-names>
</name>
, <name><surname>Finarelli</surname>
<given-names>AC</given-names>
</name>
, <etal>et al</etal>
(<year>2007</year>
) <article-title>Infection with chikungunya virus in Italy: an outbreak in a temperate region</article-title>
. <source>Lancet</source>
<volume>370</volume>
: <fpage>1840</fpage>
–<lpage>1846</lpage>
.<pub-id pub-id-type="pmid">18061059</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Casals1"><label>29</label>
<mixed-citation publication-type="journal"><name><surname>Casals</surname>
<given-names>J</given-names>
</name>
, <name><surname>Brown</surname>
<given-names>LV</given-names>
</name>
(<year>1954</year>
) <article-title>Hemagglutination with arthropod-borne viruses</article-title>
. <source>J Exp Med</source>
<volume>99</volume>
: <fpage>429</fpage>
–<lpage>449</lpage>
.<pub-id pub-id-type="pmid">13163320</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Clarke1"><label>30</label>
<mixed-citation publication-type="journal"><name><surname>Clarke</surname>
<given-names>DH</given-names>
</name>
, <name><surname>Casals</surname>
<given-names>J</given-names>
</name>
(<year>1958</year>
) <article-title>Techniques for hemagglutination and hemagglutination-inhibition with arthropod-borne viruses</article-title>
. <source>Am J Trop Med Hyg</source>
<volume>7</volume>
: <fpage>561</fpage>
–<lpage>573</lpage>
.<pub-id pub-id-type="pmid">13571577</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Roux1"><label>31</label>
<mixed-citation publication-type="journal"><name><surname>Roux</surname>
<given-names>L</given-names>
</name>
, <name><surname>Simon</surname>
<given-names>AE</given-names>
</name>
, <name><surname>Holland</surname>
<given-names>JJ</given-names>
</name>
(<year>1991</year>
) <article-title>Effects of defective interfering viruses on virus replication and pathogenesis in vitro and in vivo</article-title>
. <source>Adv Virus Res</source>
<volume>40</volume>
: <fpage>181</fpage>
–<lpage>211</lpage>
.<pub-id pub-id-type="pmid">1957718</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Kam1"><label>32</label>
<mixed-citation publication-type="journal"><name><surname>Kam</surname>
<given-names>YW</given-names>
</name>
, <name><surname>Lum</surname>
<given-names>FM</given-names>
</name>
, <name><surname>Teo</surname>
<given-names>TH</given-names>
</name>
, <name><surname>Lee</surname>
<given-names>WW</given-names>
</name>
, <name><surname>Simarmata</surname>
<given-names>D</given-names>
</name>
, <etal>et al</etal>
(<year>2012</year>
) <article-title>Early neutralizing IgG response to Chikungunya virus in infected patients targets a dominant linear epitope on the E2 glycoprotein</article-title>
. <source>EMBO Mol Med</source>
<volume>4</volume>
: <fpage>330</fpage>
–<lpage>343</lpage>
.<pub-id pub-id-type="pmid">22389221</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Cho1"><label>33</label>
<mixed-citation publication-type="journal"><name><surname>Cho</surname>
<given-names>B</given-names>
</name>
, <name><surname>Jeon</surname>
<given-names>BY</given-names>
</name>
, <name><surname>Kim</surname>
<given-names>J</given-names>
</name>
, <name><surname>Noh</surname>
<given-names>J</given-names>
</name>
, <name><surname>Park</surname>
<given-names>M</given-names>
</name>
, <etal>et al</etal>
(<year>2008</year>
) <article-title>Expression and evaluation of Chikungunya virus E1 and E2 envelope proteins for serodiagnosis of Chikungunya virus infection</article-title>
. <source>Yonsei Med J</source>
<volume>49</volume>
: <fpage>828</fpage>
–<lpage>835</lpage>
.<pub-id pub-id-type="pmid">18972604</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Cho2"><label>34</label>
<mixed-citation publication-type="journal"><name><surname>Cho</surname>
<given-names>B</given-names>
</name>
, <name><surname>Kim</surname>
<given-names>J</given-names>
</name>
, <name><surname>Cho</surname>
<given-names>JE</given-names>
</name>
, <name><surname>Jeon</surname>
<given-names>BY</given-names>
</name>
, <name><surname>Park</surname>
<given-names>S</given-names>
</name>
(<year>2008</year>
) <article-title>Expression of the capsid protein of Chikungunya virus in a baculovirus for serodiagnosis of Chikungunya disease</article-title>
. <source>J Virol Methods</source>
<volume>154</volume>
: <fpage>154</fpage>
–<lpage>159</lpage>
.<pub-id pub-id-type="pmid">18761376</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Yap1"><label>35</label>
<mixed-citation publication-type="journal"><name><surname>Yap</surname>
<given-names>G</given-names>
</name>
, <name><surname>Pok</surname>
<given-names>KY</given-names>
</name>
, <name><surname>Lai</surname>
<given-names>YL</given-names>
</name>
, <name><surname>Hapuarachchi</surname>
<given-names>HC</given-names>
</name>
, <name><surname>Chow</surname>
<given-names>A</given-names>
</name>
, <etal>et al</etal>
(<year>2010</year>
) <article-title>Evaluation of Chikungunya diagnostic assays: differences in sensitivity of serology assays in two independent outbreaks</article-title>
. <source>PLoS Negl Trop Dis</source>
<volume>4</volume>
: <fpage>e753</fpage>
.<pub-id pub-id-type="pmid">20651930</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Greene1"><label>36</label>
<mixed-citation publication-type="journal"><name><surname>Greene</surname>
<given-names>IP</given-names>
</name>
, <name><surname>Wang</surname>
<given-names>E</given-names>
</name>
, <name><surname>Deardorff</surname>
<given-names>ER</given-names>
</name>
, <name><surname>Milleron</surname>
<given-names>R</given-names>
</name>
, <name><surname>Domingo</surname>
<given-names>E</given-names>
</name>
, <etal>et al</etal>
(<year>2005</year>
) <article-title>Effect of alternating passage on adaptation of sindbis virus to vertebrate and invertebrate cells</article-title>
. <source>J Virol</source>
<volume>79</volume>
: <fpage>14253</fpage>
–<lpage>14260</lpage>
.<pub-id pub-id-type="pmid">16254360</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Taylor1"><label>37</label>
<mixed-citation publication-type="journal"><name><surname>Taylor</surname>
<given-names>WP</given-names>
</name>
, <name><surname>Marshall</surname>
<given-names>ID</given-names>
</name>
(<year>1975</year>
) <article-title>Adaptation studies with Ross River virus: retention of field level virulence</article-title>
. <source>J Gen Virol</source>
<volume>28</volume>
: <fpage>73</fpage>
–<lpage>83</lpage>
.<pub-id pub-id-type="pmid">808589</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Vasilakis1"><label>38</label>
<mixed-citation publication-type="journal"><name><surname>Vasilakis</surname>
<given-names>N</given-names>
</name>
, <name><surname>Deardorff</surname>
<given-names>ER</given-names>
</name>
, <name><surname>Kenney</surname>
<given-names>JL</given-names>
</name>
, <name><surname>Rossi</surname>
<given-names>SL</given-names>
</name>
, <name><surname>Hanley</surname>
<given-names>KA</given-names>
</name>
, <etal>et al</etal>
(<year>2009</year>
) <article-title>Mosquitoes put the brake on arbovirus evolution: experimental evolution reveals slower mutation accumulation in mosquito than vertebrate cells</article-title>
. <source>PLoS Pathog</source>
<volume>5</volume>
: <fpage>e1000467</fpage>
.<pub-id pub-id-type="pmid">19503824</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Weaver2"><label>39</label>
<mixed-citation publication-type="journal"><name><surname>Weaver</surname>
<given-names>SC</given-names>
</name>
, <name><surname>Brault</surname>
<given-names>AC</given-names>
</name>
, <name><surname>Kang</surname>
<given-names>W</given-names>
</name>
, <name><surname>Holland</surname>
<given-names>JJ</given-names>
</name>
(<year>1999</year>
) <article-title>Genetic and fitness changes accompanying adaptation of an arbovirus to vertebrate and invertebrate cells</article-title>
. <source>J Virol</source>
<volume>73</volume>
: <fpage>4316</fpage>
–<lpage>4326</lpage>
.<pub-id pub-id-type="pmid">10196330</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Zarate1"><label>40</label>
<mixed-citation publication-type="journal"><name><surname>Zarate</surname>
<given-names>S</given-names>
</name>
, <name><surname>Novella</surname>
<given-names>IS</given-names>
</name>
(<year>2004</year>
) <article-title>Vesicular stomatitis virus evolution during alternation between persistent infection in insect cells and acute infection in mammalian cells is dominated by the persistence phase</article-title>
. <source>J Virol</source>
<volume>78</volume>
: <fpage>12236</fpage>
–<lpage>12242</lpage>
.<pub-id pub-id-type="pmid">15507610</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Malet1"><label>41</label>
<mixed-citation publication-type="journal"><name><surname>Malet</surname>
<given-names>H</given-names>
</name>
, <name><surname>Coutard</surname>
<given-names>B</given-names>
</name>
, <name><surname>Jamal</surname>
<given-names>S</given-names>
</name>
, <name><surname>Dutartre</surname>
<given-names>H</given-names>
</name>
, <name><surname>Papageorgiou</surname>
<given-names>N</given-names>
</name>
, <etal>et al</etal>
(<year>2009</year>
) <article-title>The crystal structures of Chikungunya and Venezuelan equine encephalitis virus nsP3 macro domains define a conserved adenosine binding pocket</article-title>
. <source>J Virol</source>
<volume>83</volume>
: <fpage>6534</fpage>
–<lpage>6545</lpage>
.<pub-id pub-id-type="pmid">19386706</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Coffin1"><label>42</label>
<mixed-citation publication-type="journal"><name><surname>Coffin</surname>
<given-names>JM</given-names>
</name>
(<year>2008</year>
) <article-title>Attenuation by a thousand cuts</article-title>
. <source>N Engl J Med</source>
<volume>359</volume>
: <fpage>2283</fpage>
–<lpage>2285</lpage>
.<pub-id pub-id-type="pmid">19020330</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Li1"><label>43</label>
<mixed-citation publication-type="journal"><name><surname>Li</surname>
<given-names>YG</given-names>
</name>
, <name><surname>Siripanyaphinyo</surname>
<given-names>U</given-names>
</name>
, <name><surname>Tumkosit</surname>
<given-names>U</given-names>
</name>
, <name><surname>Noranate</surname>
<given-names>N</given-names>
</name>
, <name><surname>A-nuegoonpipat</surname>
<given-names>A</given-names>
</name>
, <etal>et al</etal>
(<year>2012</year>
) <article-title>Chikungunya Virus Induces a More Moderate Cytopathic Effect in Mosquito Cells than in Mammalian Cells</article-title>
. <source>Intervirology</source>
<volume>56</volume>
: <fpage>6</fpage>
–<lpage>12</lpage>
.<pub-id pub-id-type="pmid">22907160</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Brackney1"><label>44</label>
<mixed-citation publication-type="journal"><name><surname>Brackney</surname>
<given-names>DE</given-names>
</name>
, <name><surname>Scott</surname>
<given-names>JC</given-names>
</name>
, <name><surname>Sagawa</surname>
<given-names>F</given-names>
</name>
, <name><surname>Woodward</surname>
<given-names>JE</given-names>
</name>
, <name><surname>Miller</surname>
<given-names>NA</given-names>
</name>
, <etal>et al</etal>
(<year>2010</year>
) <article-title>C6/36 Aedes albopictus cells have a dysfunctional antiviral RNA interference response</article-title>
. <source>PLoS Negl Trop Dis</source>
<volume>4</volume>
: <fpage>e856</fpage>
.<pub-id pub-id-type="pmid">21049065</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Wikan1"><label>45</label>
<mixed-citation publication-type="journal"><name><surname>Wikan</surname>
<given-names>N</given-names>
</name>
, <name><surname>Sakoonwatanyoo</surname>
<given-names>P</given-names>
</name>
, <name><surname>Ubol</surname>
<given-names>S</given-names>
</name>
, <name><surname>Yoksan</surname>
<given-names>S</given-names>
</name>
, <name><surname>Smith</surname>
<given-names>DR</given-names>
</name>
(<year>2012</year>
) <article-title>Chikungunya virus infection of cell lines: analysis of the East, Central and South African lineage</article-title>
. <source>PLoS One</source>
<volume>7</volume>
: <fpage>e31102</fpage>
.<pub-id pub-id-type="pmid">22299053</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Cuevas1"><label>46</label>
<mixed-citation publication-type="journal"><name><surname>Cuevas</surname>
<given-names>JM</given-names>
</name>
, <name><surname>Domingo-Calap</surname>
<given-names>P</given-names>
</name>
, <name><surname>Sanjuan</surname>
<given-names>R</given-names>
</name>
(<year>2011</year>
) <article-title>The fitness effects of synonymous mutations in DNA and RNA viruses</article-title>
. <source>Mol Biol Evol</source>
<volume>29</volume>
: <fpage>17</fpage>
–<lpage>20</lpage>
.<pub-id pub-id-type="pmid">21771719</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Chen1"><label>47</label>
<mixed-citation publication-type="journal"><name><surname>Chen</surname>
<given-names>SL</given-names>
</name>
, <name><surname>Lee</surname>
<given-names>W</given-names>
</name>
, <name><surname>Hottes</surname>
<given-names>AK</given-names>
</name>
, <name><surname>Shapiro</surname>
<given-names>L</given-names>
</name>
, <name><surname>McAdams</surname>
<given-names>HH</given-names>
</name>
(<year>2004</year>
) <article-title>Codon usage between genomes is constrained by genome-wide mutational processes</article-title>
. <source>Proc Natl Acad Sci U S A</source>
<volume>101</volume>
: <fpage>3480</fpage>
–<lpage>3485</lpage>
.<pub-id pub-id-type="pmid">14990797</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Rima1"><label>48</label>
<mixed-citation publication-type="journal"><name><surname>Rima</surname>
<given-names>BK</given-names>
</name>
, <name><surname>McFerran</surname>
<given-names>NV</given-names>
</name>
(<year>1997</year>
) <article-title>Dinucleotide and stop codon frequencies in single-stranded RNA viruses</article-title>
. <source>J Gen Virol</source>
<volume>78</volume>
(<issue>Pt 11</issue>
)
<fpage>2859</fpage>
–<lpage>2870</lpage>
.<pub-id pub-id-type="pmid">9367373</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Zhou1"><label>49</label>
<mixed-citation publication-type="journal"><name><surname>Zhou</surname>
<given-names>J</given-names>
</name>
, <name><surname>Liu</surname>
<given-names>WJ</given-names>
</name>
, <name><surname>Peng</surname>
<given-names>SW</given-names>
</name>
, <name><surname>Sun</surname>
<given-names>XY</given-names>
</name>
, <name><surname>Frazer</surname>
<given-names>I</given-names>
</name>
(<year>1999</year>
) <article-title>Papillomavirus capsid protein expression level depends on the match between codon usage and tRNA availability</article-title>
. <source>J Virol</source>
<volume>73</volume>
: <fpage>4972</fpage>
–<lpage>4982</lpage>
.<pub-id pub-id-type="pmid">10233959</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Goodfellow1"><label>50</label>
<mixed-citation publication-type="journal"><name><surname>Goodfellow</surname>
<given-names>I</given-names>
</name>
, <name><surname>Chaudhry</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Richardson</surname>
<given-names>A</given-names>
</name>
, <name><surname>Meredith</surname>
<given-names>J</given-names>
</name>
, <name><surname>Almond</surname>
<given-names>JW</given-names>
</name>
, <etal>et al</etal>
(<year>2000</year>
) <article-title>Identification of a cis-acting replication element within the poliovirus coding region</article-title>
. <source>J Virol</source>
<volume>74</volume>
: <fpage>4590</fpage>
–<lpage>4600</lpage>
.<pub-id pub-id-type="pmid">10775595</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Tuplin1"><label>51</label>
<mixed-citation publication-type="journal"><name><surname>Tuplin</surname>
<given-names>A</given-names>
</name>
, <name><surname>Evans</surname>
<given-names>DJ</given-names>
</name>
, <name><surname>Buckley</surname>
<given-names>A</given-names>
</name>
, <name><surname>Jones</surname>
<given-names>IM</given-names>
</name>
, <name><surname>Gould</surname>
<given-names>EA</given-names>
</name>
, <etal>et al</etal>
(<year>2011</year>
) <article-title>Replication enhancer elements within the open reading frame of tick-borne encephalitis virus and their evolution within the Flavivirus genus</article-title>
. <source>Nucleic Acids Res</source>
<volume>39</volume>
: <fpage>7034</fpage>
–<lpage>48</lpage>
.<pub-id pub-id-type="pmid">21622960</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Ou1"><label>52</label>
<mixed-citation publication-type="journal"><name><surname>Ou</surname>
<given-names>JH</given-names>
</name>
, <name><surname>Strauss</surname>
<given-names>EG</given-names>
</name>
, <name><surname>Strauss</surname>
<given-names>JH</given-names>
</name>
(<year>1983</year>
) <article-title>The 5′-terminal sequences of the genomic RNAs of several alphaviruses</article-title>
. <source>J Mol Biol</source>
<volume>168</volume>
: <fpage>1</fpage>
–<lpage>15</lpage>
.<pub-id pub-id-type="pmid">6308269</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Watts1"><label>53</label>
<mixed-citation publication-type="journal"><name><surname>Watts</surname>
<given-names>JM</given-names>
</name>
, <name><surname>Dang</surname>
<given-names>KK</given-names>
</name>
, <name><surname>Gorelick</surname>
<given-names>RJ</given-names>
</name>
, <name><surname>Leonard</surname>
<given-names>CW</given-names>
</name>
, <name><surname>Bess</surname>
<given-names>JW</given-names>
<suffix>Jr</suffix>
</name>
, <etal>et al</etal>
(<year>2009</year>
) <article-title>Architecture and secondary structure of an entire HIV-1 RNA genome</article-title>
. <source>Nature</source>
<volume>460</volume>
: <fpage>711</fpage>
–<lpage>716</lpage>
.<pub-id pub-id-type="pmid">19661910</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Song1"><label>54</label>
<mixed-citation publication-type="journal"><name><surname>Song</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Liu</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Ward</surname>
<given-names>CB</given-names>
</name>
, <name><surname>Mueller</surname>
<given-names>S</given-names>
</name>
, <name><surname>Futcher</surname>
<given-names>B</given-names>
</name>
, <etal>et al</etal>
(<year>2012</year>
) <article-title>Identification of two functionally redundant RNA elements in the coding sequence of poliovirus using computer-generated design</article-title>
. <source>Proc Natl Acad Sci U S A</source>
<volume>109</volume>
: <fpage>14301</fpage>
–<lpage>14307</lpage>
.<pub-id pub-id-type="pmid">22886087</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Wichman1"><label>55</label>
<mixed-citation publication-type="journal"><name><surname>Wichman</surname>
<given-names>HA</given-names>
</name>
, <name><surname>Badgett</surname>
<given-names>MR</given-names>
</name>
, <name><surname>Scott</surname>
<given-names>LA</given-names>
</name>
, <name><surname>Boulianne</surname>
<given-names>CM</given-names>
</name>
, <name><surname>Bull</surname>
<given-names>JJ</given-names>
</name>
(<year>1999</year>
) <article-title>Different trajectories of parallel evolution during viral adaptation</article-title>
. <source>Science</source>
<volume>285</volume>
: <fpage>422</fpage>
–<lpage>424</lpage>
.<pub-id pub-id-type="pmid">10411508</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Igarashi1"><label>56</label>
<mixed-citation publication-type="journal"><name><surname>Igarashi</surname>
<given-names>A</given-names>
</name>
(<year>1978</year>
) <article-title>Isolation of a Singh's Aedes albopictus cell clone sensitive to Dengue and Chikungunya viruses</article-title>
. <source>J Gen Virol</source>
<volume>40</volume>
: <fpage>531</fpage>
–<lpage>544</lpage>
.<pub-id pub-id-type="pmid">690610</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Zhang1"><label>57</label>
<mixed-citation publication-type="journal"><name><surname>Zhang</surname>
<given-names>SP</given-names>
</name>
, <name><surname>Zubay</surname>
<given-names>G</given-names>
</name>
, <name><surname>Goldman</surname>
<given-names>E</given-names>
</name>
(<year>1991</year>
) <article-title>Low-usage codons in Escherichia coli, yeast, fruit fly and primates</article-title>
. <source>Gene</source>
<volume>105</volume>
: <fpage>61</fpage>
–<lpage>72</lpage>
.<pub-id pub-id-type="pmid">1937008</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Tsetsarkin2"><label>58</label>
<mixed-citation publication-type="journal"><name><surname>Tsetsarkin</surname>
<given-names>K</given-names>
</name>
, <name><surname>Higgs</surname>
<given-names>S</given-names>
</name>
, <name><surname>McGee</surname>
<given-names>CE</given-names>
</name>
, <name><surname>De Lamballerie</surname>
<given-names>X</given-names>
</name>
, <name><surname>Charrel</surname>
<given-names>RN</given-names>
</name>
, <etal>et al</etal>
(<year>2006</year>
) <article-title>Infectious clones of Chikungunya virus (La Reunion isolate) for vector competence studies</article-title>
. <source>Vector Borne Zoonotic Dis</source>
<volume>6</volume>
: <fpage>325</fpage>
–<lpage>337</lpage>
.<pub-id pub-id-type="pmid">17187566</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Baronti1"><label>59</label>
<mixed-citation publication-type="journal"><name><surname>Baronti</surname>
<given-names>C</given-names>
</name>
, <name><surname>Sire</surname>
<given-names>J</given-names>
</name>
, <name><surname>de Lamballerie</surname>
<given-names>X</given-names>
</name>
, <name><surname>Querat</surname>
<given-names>G</given-names>
</name>
(<year>2010</year>
) <article-title>Nonstructural NS1 proteins of several mosquito-borne Flavivirus do not inhibit TLR3 signaling</article-title>
. <source>Virology</source>
<volume>404</volume>
: <fpage>319</fpage>
–<lpage>330</lpage>
.<pub-id pub-id-type="pmid">20554300</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Hall1"><label>60</label>
<mixed-citation publication-type="journal"><name><surname>Hall</surname>
<given-names>TA</given-names>
</name>
(<year>1999</year>
) <article-title>BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT</article-title>
. <source>Nucleic Acids Symp Ser</source>
<volume>41</volume>
: <fpage>95</fpage>
–<lpage>98</lpage>
.</mixed-citation>
</ref>
<ref id="ppat.1003172-Thompson1"><label>61</label>
<mixed-citation publication-type="journal"><name><surname>Thompson</surname>
<given-names>JD</given-names>
</name>
, <name><surname>Higgins</surname>
<given-names>DG</given-names>
</name>
, <name><surname>Gibson</surname>
<given-names>TJ</given-names>
</name>
(<year>1994</year>
) <article-title>CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice</article-title>
. <source>Nucleic Acids Res</source>
<volume>22</volume>
: <fpage>4673</fpage>
–<lpage>4680</lpage>
.<pub-id pub-id-type="pmid">7984417</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Tamura1"><label>62</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>Mol Biol Evol</source>
<volume>28</volume>
: <fpage>2731</fpage>
–<lpage>9</lpage>
.<pub-id pub-id-type="pmid">21546353</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Tamura2"><label>63</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>Mol Biol Evol</source>
<volume>10</volume>
: <fpage>512</fpage>
–<lpage>526</lpage>
.<pub-id pub-id-type="pmid">8336541</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Clement1"><label>64</label>
<mixed-citation publication-type="journal"><name><surname>Clement</surname>
<given-names>M</given-names>
</name>
, <name><surname>Posada</surname>
<given-names>D</given-names>
</name>
, <name><surname>Crandall</surname>
<given-names>KA</given-names>
</name>
(<year>2000</year>
) <article-title>TCS: a computer program to estimate gene genealogies</article-title>
. <source>Mol Ecol</source>
<volume>9</volume>
: <fpage>1657</fpage>
–<lpage>1659</lpage>
.<pub-id pub-id-type="pmid">11050560</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Peris1"><label>65</label>
<mixed-citation publication-type="journal"><name><surname>Peris</surname>
<given-names>JB</given-names>
</name>
, <name><surname>Davis</surname>
<given-names>P</given-names>
</name>
, <name><surname>Cuevas</surname>
<given-names>JM</given-names>
</name>
, <name><surname>Nebot</surname>
<given-names>MR</given-names>
</name>
, <name><surname>Sanjuan</surname>
<given-names>R</given-names>
</name>
(<year>2010</year>
) <article-title>Distribution of fitness effects caused by single-nucleotide substitutions in bacteriophage f1</article-title>
. <source>Genetics</source>
<volume>185</volume>
: <fpage>603</fpage>
–<lpage>609</lpage>
.<pub-id pub-id-type="pmid">20382832</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003172-Richardson1"><label>66</label>
<mixed-citation publication-type="journal"><name><surname>Richardson</surname>
<given-names>BA</given-names>
</name>
, <name><surname>Overbaugh</surname>
<given-names>J</given-names>
</name>
(<year>2005</year>
) <article-title>Basic statistical considerations in virological experiments</article-title>
. <source>J Virol</source>
<volume>79</volume>
: <fpage>669</fpage>
–<lpage>676</lpage>
.<pub-id pub-id-type="pmid">15613294</pub-id>
</mixed-citation>
</ref>
</ref-list>
</back>
</pmc>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Wicri/Asie/explor/AustralieFrV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 002C598 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 002C598 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Wicri/Asie |area= AustralieFrV1 |flux= Pmc |étape= Corpus |type= RBID |clé= |texte= }}
![]() | This area was generated with Dilib version V0.6.33. | ![]() |