Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">The Roles of Competition and Mutation in Shaping Antigenic and Genetic Diversity in Influenza</title>
<author><name sortKey="Zinder, Daniel" sort="Zinder, Daniel" uniqKey="Zinder D" first="Daniel" last="Zinder">Daniel Zinder</name>
<affiliation><nlm:aff id="aff1"><addr-line>University of Michigan, Department of Computational Medicine and Bioinformatics, Department of Ecology and Evolutionary Biology, Ann Arbor, Michigan, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Bedford, Trevor" sort="Bedford, Trevor" uniqKey="Bedford T" first="Trevor" last="Bedford">Trevor Bedford</name>
<affiliation><nlm:aff id="aff2"><addr-line>University of Edinburgh, Institute of Evolutionary Biology, Edinburgh, United Kingdom</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Gupta, Sunetra" sort="Gupta, Sunetra" uniqKey="Gupta S" first="Sunetra" last="Gupta">Sunetra Gupta</name>
<affiliation><nlm:aff id="aff3"><addr-line>University of Oxford, Department of Zoology, Oxford, United Kingdom</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Pascual, Mercedes" sort="Pascual, Mercedes" uniqKey="Pascual M" first="Mercedes" last="Pascual">Mercedes Pascual</name>
<affiliation><nlm:aff id="aff1"><addr-line>University of Michigan, Department of Computational Medicine and Bioinformatics, Department of Ecology and Evolutionary Biology, Ann Arbor, Michigan, United States of America</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff4"><addr-line>Howard Hughes Medical Institute, University of Michigan, Ann Arbor, Michigan, United States of America</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff5"><addr-line>Santa Fe Institute, Santa Fe, New Mexico, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">23300455</idno>
<idno type="pmc">3536651</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3536651</idno>
<idno type="RBID">PMC:3536651</idno>
<idno type="doi">10.1371/journal.ppat.1003104</idno>
<date when="2013">2013</date>
<idno type="wicri:Area/Pmc/Corpus">000C11</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000C11</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">The Roles of Competition and Mutation in Shaping Antigenic and Genetic Diversity in Influenza</title>
<author><name sortKey="Zinder, Daniel" sort="Zinder, Daniel" uniqKey="Zinder D" first="Daniel" last="Zinder">Daniel Zinder</name>
<affiliation><nlm:aff id="aff1"><addr-line>University of Michigan, Department of Computational Medicine and Bioinformatics, Department of Ecology and Evolutionary Biology, Ann Arbor, Michigan, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Bedford, Trevor" sort="Bedford, Trevor" uniqKey="Bedford T" first="Trevor" last="Bedford">Trevor Bedford</name>
<affiliation><nlm:aff id="aff2"><addr-line>University of Edinburgh, Institute of Evolutionary Biology, Edinburgh, United Kingdom</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Gupta, Sunetra" sort="Gupta, Sunetra" uniqKey="Gupta S" first="Sunetra" last="Gupta">Sunetra Gupta</name>
<affiliation><nlm:aff id="aff3"><addr-line>University of Oxford, Department of Zoology, Oxford, United Kingdom</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Pascual, Mercedes" sort="Pascual, Mercedes" uniqKey="Pascual M" first="Mercedes" last="Pascual">Mercedes Pascual</name>
<affiliation><nlm:aff id="aff1"><addr-line>University of Michigan, Department of Computational Medicine and Bioinformatics, Department of Ecology and Evolutionary Biology, Ann Arbor, Michigan, United States of America</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff4"><addr-line>Howard Hughes Medical Institute, University of Michigan, Ann Arbor, Michigan, United States of America</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff5"><addr-line>Santa Fe Institute, Santa Fe, New Mexico, United States of America</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>Influenza A (H3N2) offers a well-studied, yet not fully understood, disease in terms of the interactions between pathogen population dynamics, epidemiology and genetics. A major open question is why the virus population is globally dominated by a single and very recently diverged (2–8 years) lineage. Classically, this has been modeled by limiting the generation of new successful antigenic variants, such that only a small subset of progeny acquire the necessary mutations to evade host immunity. An alternative approach was recently suggested by Recker et al. in which a limited number of antigenic variants are continuously generated, but most of these are suppressed by pre-existing host population immunity. Here we develop a framework spanning the regimes described above to explore the impact of rates of mutation and levels of competition on phylodynamic patterns. We find that the evolutionary dynamics of the subtype H3N2 influenza is most easily generated within this framework when it is mutation limited as well as being under strong immune selection at a number of epitope regions of limited diversity.</p>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Walsh, Ja" uniqKey="Walsh J">JA Walsh</name>
</author>
<author><name sortKey="Maher, C" uniqKey="Maher C">C Maher</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Fitch, Wm" uniqKey="Fitch W">WM Fitch</name>
</author>
<author><name sortKey="Bush, Rm" uniqKey="Bush R">RM Bush</name>
</author>
<author><name sortKey="Bender, Ca" uniqKey="Bender C">CA Bender</name>
</author>
<author><name sortKey="Cox, Nj" uniqKey="Cox N">NJ Cox</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bush, Rm" uniqKey="Bush R">RM Bush</name>
</author>
<author><name sortKey="Bender, Ca" uniqKey="Bender C">CA Bender</name>
</author>
<author><name sortKey="Subbarao, K" uniqKey="Subbarao K">K Subbarao</name>
</author>
<author><name sortKey="Cox, Nj" uniqKey="Cox N">NJ Cox</name>
</author>
<author><name sortKey="Fitch, Wm" uniqKey="Fitch W">WM Fitch</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Rambaut, A" uniqKey="Rambaut A">A Rambaut</name>
</author>
<author><name sortKey="Pybus, Og" uniqKey="Pybus O">OG Pybus</name>
</author>
<author><name sortKey="Nelson, Mi" uniqKey="Nelson M">MI Nelson</name>
</author>
<author><name sortKey="Viboud, C" uniqKey="Viboud C">C Viboud</name>
</author>
<author><name sortKey="Taubenberger, Jk" uniqKey="Taubenberger J">JK Taubenberger</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bedford, T" uniqKey="Bedford T">T Bedford</name>
</author>
<author><name sortKey="Cobey, S" uniqKey="Cobey S">S Cobey</name>
</author>
<author><name sortKey="Pascual, M" uniqKey="Pascual M">M Pascual</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gerhard, W" uniqKey="Gerhard W">W Gerhard</name>
</author>
<author><name sortKey="Webster, Rg" uniqKey="Webster R">RG Webster</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hay, Aj" uniqKey="Hay A">AJ Hay</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nakajima, K" uniqKey="Nakajima K">K Nakajima</name>
</author>
<author><name sortKey="Nobusawa, E" uniqKey="Nobusawa E">E Nobusawa</name>
</author>
<author><name sortKey="Nagy, A" uniqKey="Nagy A">A Nagy</name>
</author>
<author><name sortKey="Nakajima, S" uniqKey="Nakajima S">S Nakajima</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Fitch, Wm" uniqKey="Fitch W">WM Fitch</name>
</author>
<author><name sortKey="Leiter, Jm" uniqKey="Leiter J">JM Leiter</name>
</author>
<author><name sortKey="Li, Xq" uniqKey="Li X">XQ Li</name>
</author>
<author><name sortKey="Palese, P" uniqKey="Palese P">P Palese</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Suzuki, Y" uniqKey="Suzuki Y">Y Suzuki</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wolf, Yi" uniqKey="Wolf Y">YI Wolf</name>
</author>
<author><name sortKey="Viboud, C" uniqKey="Viboud C">C Viboud</name>
</author>
<author><name sortKey="Holmes, Ec" uniqKey="Holmes E">EC Holmes</name>
</author>
<author><name sortKey="Koonin, Ev" uniqKey="Koonin E">EV Koonin</name>
</author>
<author><name sortKey="Lipman, Dj" uniqKey="Lipman D">DJ Lipman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lindstrom, Se" uniqKey="Lindstrom S">SE Lindstrom</name>
</author>
<author><name sortKey="Hiromoto, Y" uniqKey="Hiromoto Y">Y Hiromoto</name>
</author>
<author><name sortKey="Nishimura, H" uniqKey="Nishimura H">H Nishimura</name>
</author>
<author><name sortKey="Saito, T" uniqKey="Saito T">T Saito</name>
</author>
<author><name sortKey="Nerome, R" uniqKey="Nerome R">R Nerome</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tria, F" uniqKey="Tria F">F Tria</name>
</author>
<author><name sortKey="Laessig, M" uniqKey="Laessig M">M Laessig</name>
</author>
<author><name sortKey="Peliti, L" uniqKey="Peliti L">L Peliti</name>
</author>
<author><name sortKey="Franz, S" uniqKey="Franz S">S Franz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ferguson, Nm" uniqKey="Ferguson N">NM Ferguson</name>
</author>
<author><name sortKey="Galvani, Ap" uniqKey="Galvani A">AP Galvani</name>
</author>
<author><name sortKey="Bush, Rm" uniqKey="Bush R">RM Bush</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Koelle, K" uniqKey="Koelle K">K Koelle</name>
</author>
<author><name sortKey="Cobey, S" uniqKey="Cobey S">S Cobey</name>
</author>
<author><name sortKey="Grenfell, B" uniqKey="Grenfell B">B Grenfell</name>
</author>
<author><name sortKey="Pascual, M" uniqKey="Pascual M">M Pascual</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Smith, Dj" uniqKey="Smith D">DJ Smith</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Andreasen, V" uniqKey="Andreasen V">V Andreasen</name>
</author>
<author><name sortKey="Sasaki, A" uniqKey="Sasaki A">A Sasaki</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kryazhimskiy, S" uniqKey="Kryazhimskiy S">S Kryazhimskiy</name>
</author>
<author><name sortKey="Dieckmann, U" uniqKey="Dieckmann U">U Dieckmann</name>
</author>
<author><name sortKey="Levin, Sa" uniqKey="Levin S">SA Levin</name>
</author>
<author><name sortKey="Dushoff, J" uniqKey="Dushoff J">J Dushoff</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bedford, T" uniqKey="Bedford T">T Bedford</name>
</author>
<author><name sortKey="Rambaut, A" uniqKey="Rambaut A">A Rambaut</name>
</author>
<author><name sortKey="Pascual, M" uniqKey="Pascual M">M Pascual</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Recker, M" uniqKey="Recker M">M Recker</name>
</author>
<author><name sortKey="Pybus, Og" uniqKey="Pybus O">OG Pybus</name>
</author>
<author><name sortKey="Nee, S" uniqKey="Nee S">S Nee</name>
</author>
<author><name sortKey="Gupta, S" uniqKey="Gupta S">S Gupta</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Russell, Ca" uniqKey="Russell C">CA Russell</name>
</author>
<author><name sortKey="Jones, Tc" uniqKey="Jones T">TC Jones</name>
</author>
<author><name sortKey="Barr, Ig" uniqKey="Barr I">IG Barr</name>
</author>
<author><name sortKey="Cox, Nj" uniqKey="Cox N">NJ Cox</name>
</author>
<author><name sortKey="Garten, Rj" uniqKey="Garten R">RJ Garten</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kuiken, T" uniqKey="Kuiken T">T Kuiken</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gupta, S" uniqKey="Gupta S">S Gupta</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Buckee, Co" uniqKey="Buckee C">CO Buckee</name>
</author>
<author><name sortKey="Jolley, Ka" uniqKey="Jolley K">KA Jolley</name>
</author>
<author><name sortKey="Recker, M" uniqKey="Recker M">M Recker</name>
</author>
<author><name sortKey="Penman, B" uniqKey="Penman B">B Penman</name>
</author>
<author><name sortKey="Kriz, P" uniqKey="Kriz P">P Kriz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="May, Rm" uniqKey="May R">RM May</name>
</author>
<author><name sortKey="Nowak, Ma" uniqKey="Nowak M">MA Nowak</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nowak, Ma" uniqKey="Nowak M">MA Nowak</name>
</author>
<author><name sortKey="May, Rm" uniqKey="May R">RM May</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Volz, Em" uniqKey="Volz E">EM Volz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nowak, Ma" uniqKey="Nowak M">MA Nowak</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mcdonald, Jh" uniqKey="Mcdonald J">JH McDonald</name>
</author>
<author><name sortKey="Kreitman, M" uniqKey="Kreitman M">M Kreitman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bush, Rm" uniqKey="Bush R">RM Bush</name>
</author>
<author><name sortKey="Fitch, Wm" uniqKey="Fitch W">WM Fitch</name>
</author>
<author><name sortKey="Bender, Ca" uniqKey="Bender C">CA Bender</name>
</author>
<author><name sortKey="Cox, Nj" uniqKey="Cox N">NJ Cox</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bhatt, S" uniqKey="Bhatt S">S Bhatt</name>
</author>
<author><name sortKey="Katzourakis, A" uniqKey="Katzourakis A">A Katzourakis</name>
</author>
<author><name sortKey="Pybus, Og" uniqKey="Pybus O">OG Pybus</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Chew, Ft" uniqKey="Chew F">FT Chew</name>
</author>
<author><name sortKey="Doraisingham, S" uniqKey="Doraisingham S">S Doraisingham</name>
</author>
<author><name sortKey="Ling, Ae" uniqKey="Ling A">AE Ling</name>
</author>
<author><name sortKey="Kumarasinghe, G" uniqKey="Kumarasinghe G">G Kumarasinghe</name>
</author>
<author><name sortKey="Lee, Bw" uniqKey="Lee B">BW Lee</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Huang, Jw" uniqKey="Huang J">JW Huang</name>
</author>
<author><name sortKey="Yang, Jm" uniqKey="Yang J">JM Yang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Fleming, Dm" uniqKey="Fleming D">DM Fleming</name>
</author>
<author><name sortKey="Ayres, Jg" uniqKey="Ayres J">JG Ayres</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Monto, As" uniqKey="Monto A">AS Monto</name>
</author>
<author><name sortKey="Sullivan, Km" uniqKey="Sullivan K">KM Sullivan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Karpova, Ls" uniqKey="Karpova L">LS Karpova</name>
</author>
<author><name sortKey="Marinich, Ig" uniqKey="Marinich I">IG Marinich</name>
</author>
<author><name sortKey="Krainova, Ti" uniqKey="Krainova T">TI Krainova</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Saltyte Benth, J" uniqKey="Saltyte Benth J">J Saltyte Benth</name>
</author>
<author><name sortKey="Hofoss, D" uniqKey="Hofoss D">D Hofoss</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sumi, A" uniqKey="Sumi A">A Sumi</name>
</author>
<author><name sortKey="Kamo, K" uniqKey="Kamo K">K Kamo</name>
</author>
<author><name sortKey="Ohtomo, N" uniqKey="Ohtomo N">N Ohtomo</name>
</author>
<author><name sortKey="Mise, K" uniqKey="Mise K">K Mise</name>
</author>
<author><name sortKey="Kobayashi, N" uniqKey="Kobayashi N">N Kobayashi</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Du, X" uniqKey="Du X">X Du</name>
</author>
<author><name sortKey="Dong, L" uniqKey="Dong L">L Dong</name>
</author>
<author><name sortKey="Lan, Y" uniqKey="Lan Y">Y Lan</name>
</author>
<author><name sortKey="Peng, Y" uniqKey="Peng Y">Y Peng</name>
</author>
<author><name sortKey="Wu, A" uniqKey="Wu A">A Wu</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bedford, T" uniqKey="Bedford T">T Bedford</name>
</author>
<author><name sortKey="Cobey, S" uniqKey="Cobey S">S Cobey</name>
</author>
<author><name sortKey="Beerli, P" uniqKey="Beerli P">P Beerli</name>
</author>
<author><name sortKey="Pascual, M" uniqKey="Pascual M">M Pascual</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Creanza, N" uniqKey="Creanza N">N Creanza</name>
</author>
<author><name sortKey="Schwarz, Js" uniqKey="Schwarz J">JS Schwarz</name>
</author>
<author><name sortKey="Cohen, Je" uniqKey="Cohen J">JE Cohen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Minayev, P" uniqKey="Minayev P">P Minayev</name>
</author>
<author><name sortKey="Ferguson, N" uniqKey="Ferguson N">N Ferguson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Yamashita, M" uniqKey="Yamashita M">M Yamashita</name>
</author>
<author><name sortKey="Krystal, M" uniqKey="Krystal M">M Krystal</name>
</author>
<author><name sortKey="Fitch, Wm" uniqKey="Fitch W">WM Fitch</name>
</author>
<author><name sortKey="Palese, P" uniqKey="Palese P">P Palese</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hensley, Se" uniqKey="Hensley S">SE Hensley</name>
</author>
<author><name sortKey="Das, Sr" uniqKey="Das S">SR Das</name>
</author>
<author><name sortKey="Bailey, Al" uniqKey="Bailey A">AL Bailey</name>
</author>
<author><name sortKey="Schmidt, Lm" uniqKey="Schmidt L">LM Schmidt</name>
</author>
<author><name sortKey="Hickman, Hd" uniqKey="Hickman H">HD Hickman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hancock, K" uniqKey="Hancock K">K Hancock</name>
</author>
<author><name sortKey="Veguilla, V" uniqKey="Veguilla V">V Veguilla</name>
</author>
<author><name sortKey="Lu, X" uniqKey="Lu X">X Lu</name>
</author>
<author><name sortKey="Zhong, W" uniqKey="Zhong W">W Zhong</name>
</author>
<author><name sortKey="Butler, En" uniqKey="Butler E">EN Butler</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Itoh, Y" uniqKey="Itoh Y">Y Itoh</name>
</author>
<author><name sortKey="Shinya, K" uniqKey="Shinya K">K Shinya</name>
</author>
<author><name sortKey="Kiso, M" uniqKey="Kiso M">M Kiso</name>
</author>
<author><name sortKey="Watanabe, T" uniqKey="Watanabe T">T Watanabe</name>
</author>
<author><name sortKey="Sakoda, Y" uniqKey="Sakoda Y">Y Sakoda</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bandaranayake, D" uniqKey="Bandaranayake D">D Bandaranayake</name>
</author>
<author><name sortKey="Huang, Qs" uniqKey="Huang Q">QS Huang</name>
</author>
<author><name sortKey="Bissielo, A" uniqKey="Bissielo A">A Bissielo</name>
</author>
<author><name sortKey="Wood, T" uniqKey="Wood T">T Wood</name>
</author>
<author><name sortKey="Mackereth, G" uniqKey="Mackereth G">G Mackereth</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chi, Cy" uniqKey="Chi C">CY Chi</name>
</author>
<author><name sortKey="Liu, Cc" uniqKey="Liu C">CC Liu</name>
</author>
<author><name sortKey="Lin, Cc" uniqKey="Lin C">CC Lin</name>
</author>
<author><name sortKey="Wang, Hc" uniqKey="Wang H">HC Wang</name>
</author>
<author><name sortKey="Cheng, Yt" uniqKey="Cheng Y">YT Cheng</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wikramaratna, Ps" uniqKey="Wikramaratna P">PS Wikramaratna</name>
</author>
<author><name sortKey="Sandeman, M" uniqKey="Sandeman M">M Sandeman</name>
</author>
<author><name sortKey="Recker, M" uniqKey="Recker M">M Recker</name>
</author>
<author><name sortKey="Gupta, S" uniqKey="Gupta S">S Gupta</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tsuchiya, E" uniqKey="Tsuchiya E">E Tsuchiya</name>
</author>
<author><name sortKey="Sugawara, K" uniqKey="Sugawara K">K Sugawara</name>
</author>
<author><name sortKey="Hongo, S" uniqKey="Hongo S">S Hongo</name>
</author>
<author><name sortKey="Matsuzaki, Y" uniqKey="Matsuzaki Y">Y Matsuzaki</name>
</author>
<author><name sortKey="Muraki, Y" uniqKey="Muraki Y">Y Muraki</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Reichert, T" uniqKey="Reichert T">T Reichert</name>
</author>
<author><name sortKey="Chowell, G" uniqKey="Chowell G">G Chowell</name>
</author>
<author><name sortKey="Nishiura, H" uniqKey="Nishiura H">H Nishiura</name>
</author>
<author><name sortKey="Christensen, Ra" uniqKey="Christensen R">RA Christensen</name>
</author>
<author><name sortKey="Mccullers, Ja" uniqKey="Mccullers J">JA McCullers</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bui, Hh" uniqKey="Bui H">HH Bui</name>
</author>
<author><name sortKey="Peters, B" uniqKey="Peters B">B Peters</name>
</author>
<author><name sortKey="Assarsson, E" uniqKey="Assarsson E">E Assarsson</name>
</author>
<author><name sortKey="Mbawuike, I" uniqKey="Mbawuike I">I Mbawuike</name>
</author>
<author><name sortKey="Sette, A" uniqKey="Sette A">A Sette</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wang, Tt" uniqKey="Wang T">TT Wang</name>
</author>
<author><name sortKey="Tan, Gs" uniqKey="Tan G">GS Tan</name>
</author>
<author><name sortKey="Hai, R" uniqKey="Hai R">R Hai</name>
</author>
<author><name sortKey="Pica, N" uniqKey="Pica N">N Pica</name>
</author>
<author><name sortKey="Petersen, E" uniqKey="Petersen E">E Petersen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Viboud, C" uniqKey="Viboud C">C Viboud</name>
</author>
<author><name sortKey="Alonso, Wj" uniqKey="Alonso W">WJ Alonso</name>
</author>
<author><name sortKey="Simonsen, L" uniqKey="Simonsen L">L Simonsen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Drummond, Aj" uniqKey="Drummond A">AJ Drummond</name>
</author>
<author><name sortKey="Suchard, Ma" uniqKey="Suchard M">MA Suchard</name>
</author>
<author><name sortKey="Xie, D" uniqKey="Xie D">D Xie</name>
</author>
<author><name sortKey="Rambaut, A" uniqKey="Rambaut A">A Rambaut</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lemey, P" uniqKey="Lemey P">P Lemey</name>
</author>
<author><name sortKey="Rambaut, A" uniqKey="Rambaut A">A Rambaut</name>
</author>
<author><name sortKey="Drummond, Aj" uniqKey="Drummond A">AJ Drummond</name>
</author>
<author><name sortKey="Suchard, Ma" uniqKey="Suchard M">MA Suchard</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">23300455</article-id>
<article-id pub-id-type="pmc">3536651</article-id>
<article-id pub-id-type="publisher-id">PPATHOGENS-D-12-01327</article-id>
<article-id pub-id-type="doi">10.1371/journal.ppat.1003104</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>Computational Biology</subject>
</subj-group>
<subj-group><subject>Ecology</subject>
<subj-group><subject>Evolutionary Ecology</subject>
</subj-group>
</subj-group>
<subj-group><subject>Evolutionary Biology</subject>
</subj-group>
<subj-group><subject>Microbiology</subject>
<subj-group><subject>Virology</subject>
<subj-group><subject>Viral Evolution</subject>
</subj-group>
</subj-group>
</subj-group>
<subj-group><subject>Population Biology</subject>
<subj-group><subject>Epidemiology</subject>
<subj-group><subject>Infectious Disease Epidemiology</subject>
</subj-group>
</subj-group>
</subj-group>
</subj-group>
</article-categories>
<title-group><article-title>The Roles of Competition and Mutation in Shaping Antigenic and Genetic Diversity in Influenza</article-title>
<alt-title alt-title-type="running-head">Competition and Mutation in H3N2 Phylodynamics</alt-title>
</title-group>
<contrib-group><contrib contrib-type="author"><name><surname>Zinder</surname>
<given-names>Daniel</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Bedford</surname>
<given-names>Trevor</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Gupta</surname>
<given-names>Sunetra</given-names>
</name>
<xref ref-type="aff" rid="aff3"><sup>3</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Pascual</surname>
<given-names>Mercedes</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff4"><sup>4</sup>
</xref>
<xref ref-type="aff" rid="aff5"><sup>5</sup>
</xref>
<xref ref-type="corresp" rid="cor1"><sup>*</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1"><label>1</label>
<addr-line>University of Michigan, Department of Computational Medicine and Bioinformatics, Department of Ecology and Evolutionary Biology, Ann Arbor, Michigan, United States of America</addr-line>
</aff>
<aff id="aff2"><label>2</label>
<addr-line>University of Edinburgh, Institute of Evolutionary Biology, Edinburgh, United Kingdom</addr-line>
</aff>
<aff id="aff3"><label>3</label>
<addr-line>University of Oxford, Department of Zoology, Oxford, United Kingdom</addr-line>
</aff>
<aff id="aff4"><label>4</label>
<addr-line>Howard Hughes Medical Institute, University of Michigan, Ann Arbor, Michigan, United States of America</addr-line>
</aff>
<aff id="aff5"><label>5</label>
<addr-line>Santa Fe Institute, Santa Fe, New Mexico, United States of America</addr-line>
</aff>
<contrib-group><contrib contrib-type="editor"><name><surname>Ferguson</surname>
<given-names>Neil</given-names>
</name>
<role>Editor</role>
<xref ref-type="aff" rid="edit1"></xref>
</contrib>
</contrib-group>
<aff id="edit1"><addr-line>Imperial College London, United Kingdom</addr-line>
</aff>
<author-notes><corresp id="cor1">* E-mail: <email>pascual@umich.edu</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: DZ TB SG MP. Performed the experiments: DZ TB. Analyzed the data: DZ TB MP. Contributed reagents/materials/analysis tools: DZ TB. Wrote the paper: DZ TB SG MP.</p>
</fn>
</author-notes>
<pub-date pub-type="collection"><month>1</month>
<year>2013</year>
</pub-date>
<pub-date pub-type="epub"><day>3</day>
<month>1</month>
<year>2013</year>
</pub-date>
<volume>9</volume>
<issue>1</issue>
<elocation-id>e1003104</elocation-id>
<history><date date-type="received"><day>4</day>
<month>6</month>
<year>2012</year>
</date>
<date date-type="accepted"><day>10</day>
<month>11</month>
<year>2012</year>
</date>
</history>
<permissions><copyright-statement>© 2013 Zinder et al</copyright-statement>
<copyright-year>2013</copyright-year>
<copyright-holder>Zinder 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>Influenza A (H3N2) offers a well-studied, yet not fully understood, disease in terms of the interactions between pathogen population dynamics, epidemiology and genetics. A major open question is why the virus population is globally dominated by a single and very recently diverged (2–8 years) lineage. Classically, this has been modeled by limiting the generation of new successful antigenic variants, such that only a small subset of progeny acquire the necessary mutations to evade host immunity. An alternative approach was recently suggested by Recker et al. in which a limited number of antigenic variants are continuously generated, but most of these are suppressed by pre-existing host population immunity. Here we develop a framework spanning the regimes described above to explore the impact of rates of mutation and levels of competition on phylodynamic patterns. We find that the evolutionary dynamics of the subtype H3N2 influenza is most easily generated within this framework when it is mutation limited as well as being under strong immune selection at a number of epitope regions of limited diversity.</p>
</abstract>
<abstract abstract-type="summary"><title>Author Summary</title>
<p>Influenza A (H3N2) has circulated in the human population since 1968 causing considerable annual morbidity and mortality worldwide. Despite the rapid evolution of the hemagglutinin (HA) protein and strong diversifying selection, the global virus population is characterized by a low standing diversity, evident in the serial replacement of antigenic types and in the ‘cactus-like’ structure of its genealogical tree. Elucidating the mechanisms behind these puzzling patterns is key to understanding the evolution of seasonal (H3N2) influenza. One recent epidemiological model proposes a restricted set of antigenic types whose waves of dominance result from frequency-dependent immune selection. Here we develop a model of limited antigenic diversity that explicitly incorporates mutational processes, and use it to address, first, whether this type of antigenic space is capable of generating the characteristic phylogeny of HA sequences, and second, whether the dynamics of (H3N2) influenza are primarily limited by the arrival of mutations or by the opening of antigenic niches. We conclude that a limited antigenic space can explain the observed phylogenetic patterns and that a limited mutation rate is a key property underlying the dynamics of (H3N2) influenza. Our study provides a general framework for assessing the relative roles of selection and mutation in a variety of infectious disease systems.</p>
</abstract>
<funding-group><funding-statement>MP was supported by a Centennial Fellowship of the James S. McDonnell Foundation (<ext-link ext-link-type="uri" xlink:href="http://www.jsmf.org/">http://www.jsmf.org/</ext-link>
), SG, by the European Research Council (ERC Advanced Grant - Diversity, <ext-link ext-link-type="uri" xlink:href="http://erc.europa.eu/advanced-grants">http://erc.europa.eu/advanced-grants</ext-link>
), and TB, by the Howard Hughes Medical Institute, a Long-Term Fellowship from the European Molecular Biology Organization and a Newton International Fellowship from the Royal Society. MP is an investigator of the Howard Hughes Medical Institute and SG is a Royal Society Wolfson Research Fellow (<ext-link ext-link-type="uri" xlink:href="http://royalsociety.org/grants/schemes/wolfson-research-merit/">http://royalsociety.org/grants/schemes/wolfson-research-merit/</ext-link>
). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.</funding-statement>
</funding-group>
<counts><page-count count="13"></page-count>
</counts>
</article-meta>
</front>
<body><sec id="s1"><title>Introduction</title>
<p>Influenza viruses are classified into types A–C, among which influenza A is the most pathogenic. These viruses cause between a quarter to half a million deaths worldwide <xref rid="ppat.1003104-1" ref-type="bibr">[1]</xref>
and tens of thousands of deaths in the US during annual epidemics <xref rid="ppat.1003104-CDC1" ref-type="bibr">[2]</xref>
. The economic burden of seasonal influenza in the US is estimated at more than ten billion dollars in healthcare costs alone <xref rid="ppat.1003104-Walsh1" ref-type="bibr">[3]</xref>
.</p>
<p>The major targets of humoral immunity against influenza A are its envelope glycoproteins, hemagglutinin (HA) and neuraminidase (NA); these form the basis of its crude classification into subtypes H1N1, H2N2 and H3N2 etc. Since its emergence in 1968, influenza A (H3N2) has continually circulated in the human population. The phylogeny of its HA protein (<xref ref-type="fig" rid="ppat-1003104-g001">Figure 1</xref>
) shows a distinctive ‘cactus-like’ shape with a narrow, usually single-trunked, tree <xref rid="ppat.1003104-Fitch1" ref-type="bibr">[4]</xref>
, . The ‘narrowness’ of the tree is derived from the fact that contemporaneous H3 proteins share a single common ancestor 2–8 years in the past <xref rid="ppat.1003104-Rambaut1" ref-type="bibr">[6]</xref>
, <xref rid="ppat.1003104-Bedford1" ref-type="bibr">[7]</xref>
. This short time is unique to H3N2 given its global spread and its high prevalence and incidence <xref rid="ppat.1003104-Bedford1" ref-type="bibr">[7]</xref>
.</p>
<fig id="ppat-1003104-g001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1003104.g001</object-id>
<label>Figure 1</label>
<caption><title>Phylogenetic tree reconstruction of H3 depicting major antigenic clusters and including the associated pairwise diversity.</title>
<p>Phylogenetic tree with highest posterior likelihood was reconstructed using 377 representative sequences sampled between 1968–2009. Colors represent estimated antigenic clusters (Hong-Kong 1968 – Perth 2009). Approximately half of the samples include an established cluster annotation <xref rid="ppat.1003104-Smith1" ref-type="bibr">[18]</xref>
and three additional clusters relating to: California 2004, Brisbane 2007, and Perth 2009. Additional sequences were sampled uniformly overtime on a bi-annual scale. Phylogenetic tree was reconstructed using Bayesian MCMC analysis <xref rid="ppat.1003104-Drummond1" ref-type="bibr">[58]</xref>
, <xref rid="ppat.1003104-Lemey1" ref-type="bibr">[59]</xref>
and includes state reconstruction for unannotated sequences and ancestral sequences. Diversity skyline was calculated for the same representative tree. Branches with colors differing from their main neighboring cluster represent uncertainty in the reconstruction, rather than actual cluster changes.</p>
</caption>
<graphic xlink:href="ppat.1003104.g001"></graphic>
</fig>
<p>The classical view of influenza evolution is one of antigenic drift <xref rid="ppat.1003104-Gerhard1" ref-type="bibr">[8]</xref>
, <xref rid="ppat.1003104-Hay1" ref-type="bibr">[9]</xref>
, <xref rid="ppat.1003104-Nakajima1" ref-type="bibr">[10]</xref>
in which antigenic change continually and gradually accumulates in the virus through the influence of selection by way of changes to the HA and NA proteins. By itself, the ‘cactus-like’ structure of the A/H3N2 phylogenetic tree suggests the presence of adaptive evolution <xref rid="ppat.1003104-Bedford1" ref-type="bibr">[7]</xref>
and several studies have provided evidence for positive selection <xref rid="ppat.1003104-Bush1" ref-type="bibr">[5]</xref>
, <xref rid="ppat.1003104-Fitch2" ref-type="bibr">[11]</xref>
, <xref rid="ppat.1003104-Suzuki1" ref-type="bibr">[12]</xref>
, <xref rid="ppat.1003104-Wolf1" ref-type="bibr">[13]</xref>
, <xref rid="ppat.1003104-Lindstrom1" ref-type="bibr">[14]</xref>
. However, it is difficult to explain the limited standing diversity of influenza <xref rid="ppat.1003104-Tria1" ref-type="bibr">[15]</xref>
, <xref rid="ppat.1003104-Ferguson1" ref-type="bibr">[16]</xref>
, and the empirical evidence for discontinuous antigenic change <xref rid="ppat.1003104-Koelle1" ref-type="bibr">[17]</xref>
, <xref rid="ppat.1003104-Smith1" ref-type="bibr">[18]</xref>
, under a general antigenic drift framework. Multiple epidemiological hypotheses have been advanced to reconcile the these observations with a process of continual antigenic divergence including short-lived strain-transcending immunity <xref rid="ppat.1003104-Tria1" ref-type="bibr">[15]</xref>
, <xref rid="ppat.1003104-Ferguson1" ref-type="bibr">[16]</xref>
, <xref rid="ppat.1003104-Andreasen1" ref-type="bibr">[19]</xref>
, epochal or punctuated evolution <xref rid="ppat.1003104-Wolf1" ref-type="bibr">[13]</xref>
, <xref rid="ppat.1003104-Koelle1" ref-type="bibr">[17]</xref>
, trait-space reduction <xref rid="ppat.1003104-Kryazhimskiy1" ref-type="bibr">[20]</xref>
and canalized evolution <xref rid="ppat.1003104-Bedford2" ref-type="bibr">[21]</xref>
.</p>
<p>A competing hypothesis advanced by Recker et al. <xref rid="ppat.1003104-Recker1" ref-type="bibr">[22]</xref>
eschews the paradigm of antigenic drift, instead considering that, owing to functional constraints on the defining epitopes, the virus population is limited phenotypically to a restricted set of antigenic types. Antigenic types replace each other with waves of dominance resulting from frequency-dependent immune mediated selection as “niches” in antigenic space are dynamically generated and are exploited by the existing virus population. In its original implementation, this model assumes that all antigenic types remain present in the population in low frequencies, as an approximation to the idea that they can be generated by mutation from preexisting strains at a sufficient rate as not to limit the emergence of a type favored by selection. Thus, the model describes in practice a case where influenza outbreaks are caused by host immune selection in a manner that is not limited by the rate of antigenic mutations. Although patterns of turnover are consistent with those observed for H3N2, it is not clear whether the characteristic phylogenetic trees can be generated by this model.</p>
<p>Here, we have attempted to resolve this question using a large-scale individual-based simulation of epidemiological and evolutionary dynamics that allows the complete phylogenetic tracking of a virus population characterized by defined repertoires of polymorphic epitopes. Our model is based on the multi-locus structure employed by <xref rid="ppat.1003104-Recker1" ref-type="bibr">[22]</xref>
with host immunity operating at an epitope-specific level. When contacted by a virus, a host's risk of infection is determined by the number of alleles/epitopes recognized by its immune system. We also introduce the possibility of a long-lived strain-transcending component to the model. Thus, competition between strains is determined both by the number of shared epitopes and a variable level of generalized immunity. Our model differs in this regard from that of Recker et al. <xref rid="ppat.1003104-Recker1" ref-type="bibr">[22]</xref>
which does not permit full cross-protection except in the case of having experienced the exact same combination of epitopes, a feature that implicitly accounts for the effect of highly variable epitopes unique to each strain.</p>
<p>This model structure allows us to make inferences about the roles of mutation and competition in a more general context. Models of antigenic dynamics tend to polarize between those in which the availability of antigenic types dictate the dynamics <xref rid="ppat.1003104-Wolf1" ref-type="bibr">[13]</xref>
, <xref rid="ppat.1003104-Koelle1" ref-type="bibr">[17]</xref>
, <xref rid="ppat.1003104-Russell1" ref-type="bibr">[23]</xref>
, <xref rid="ppat.1003104-Kuiken1" ref-type="bibr">[24]</xref>
, and those where host immune-mediated selection is the only driver <xref rid="ppat.1003104-Gupta1" ref-type="bibr">[25]</xref>
, <xref rid="ppat.1003104-Buckee1" ref-type="bibr">[26]</xref>
. We refer to the latter regime, where antigenic change is constrained by host population immunity, as selection limited, whereas the former, in which the availability of antigenic mutations pose the rate limiting step, is described as mutation limited. The approach we take in this paper offers a tool for locating influenza on this continuum and would easily generalize to other antigenically diverse infectious diseases.</p>
</sec>
<sec id="s2"><title>Results</title>
<p>We use an individual-based SIR model, explicitly tracking the chains of infection of viral lineages as well as the antigenic phenotype of every virus in the population (<xref ref-type="fig" rid="ppat-1003104-g002">Figure 2A</xref>
). Thus, our model explicitly tracks viral genealogy rather than conducting phylogenetic inference, and therefore does not include any genotype to phenotype map.</p>
<fig id="ppat-1003104-g002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1003104.g002</object-id>
<label>Figure 2</label>
<caption><title>Methods.</title>
<p>(<bold>A</bold>
) Virus genealogy is tracked at the inter-host level. The genealogy is periodically sampled and the resulting tree is used for analysis. <bold>(B)</bold>
Hosts acquire immunity to viral epitopes following infection. Fully naïve hosts are always infected following contact at rate β. The risk of reinfection is based on the similarity to previously encountered strains, as measured through the number of previously encountered epitopes: <inline-formula><inline-graphic xlink:href="ppat.1003104.e001.jpg"></inline-graphic>
</inline-formula>
and is at most 100%. Where <italic>f</italic>
is the fraction of previously encountered epitopes and σ is the strength of crossimmunity. Lower <italic>σ</italic>
values correspond to weaker competition between strains. A form of generalized immunity is attained for σ>0.8 in the five epitope case, relating to a reduced risk of reinfection following previous exposure to any strain. <bold>(C)</bold>
Mean pairwise genealogical diversity π is measured by averaging the pairwise distance in years between random contemporaneous samples on the genealogical tree. <bold>(D)</bold>
The MK related index is calculated as the ratio of the antigenic mutation rate on the trunk of the genealogy (red) versus the antigenic mutation rate on the sidebranches (black). The trunk of the genealogy was determined by tracing back viral lineages that survived until the end of the simulation and excluding the last 5 years. Antigenic changes are represented by color changes on tree branches (top-tree). The rate of antigenic change on the sidebranches is calculated as the number of antigenic changes on the sidebranches divided by the total length of the side branches in years. The rate of antigenic change on the trunk is calculated as the number of antigenic changes on the trunk divided by the total length of the trunk in years.</p>
</caption>
<graphic xlink:href="ppat.1003104.g002"></graphic>
</fig>
<p>A host population of constant size N is simulated with an equal birth and death rate of μ. Infected hosts randomly contact other hosts in the population at rate β. The risk of infection is based on crossimmunity interactions based on their previous unordered infection history <xref rid="ppat.1003104-May1" ref-type="bibr">[27]</xref>
, <xref rid="ppat.1003104-Nowak1" ref-type="bibr">[28]</xref>
. Following infection with a virus displaying a specific antigenic phenotype, hosts acquire partial protection against reinfection with viruses sharing epitopes, and full protection against the exact same strain. There is no super-infection in the model; while infected, hosts are protected from co-infection with other strains. Hosts recover from infection at a constant rate ν. (see <xref ref-type="sec" rid="s4">Methods</xref>
for details of the model).</p>
<sec id="s2a"><title>Antigenic-Phenotype Replacement and Single Phenotype Dominance in an Evolution-Free Framework</title>
<p>To explore the epidemiological dynamics of our model in isolation, we implemented a parameterization lacking mutation, in which extinction was preempted by maintaining at least one carrier for each antigenic-phenotype. Different colors (<xref ref-type="fig" rid="ppat-1003104-g003">Figure 3</xref>
) represent the prevalence of different antigenic-phenotypes. Here, the antigenic repertoire is derived from combinations of variants at 5 distinct epitopes (see <xref ref-type="sec" rid="s4">Methods</xref>
for full description of epidemiological parameters).</p>
<fig id="ppat-1003104-g003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1003104.g003</object-id>
<label>Figure 3</label>
<caption><title>Changes in the proportions of hosts that are infectious with different strains within a 2 variants per epitope, 5 epitope system in an “evolutionary free” framework: for all of the possible 32 strains the existence of at least one carrier was assured and no antigenic mutations were introduced.</title>
<p>The superimposed time series were smoothed and ordered back to front by peak prevalence, maintaining the least prevalent strain in the front. The 3th highest peaking strain was outlined as an example. Single strain dominance was calculated based on the quantity ε from <xref rid="ppat.1003104-Recker1" ref-type="bibr">[22]</xref>
. Major peaks of incidence are generally associated with one or two dominant antigenic-phenotypes with ε = 0.36±0.07 (mean ± standard-deviation across 5 simulations) and a myriad of lower prevalence ones. Antigenic-phenotypes reemerge with alternating frequency. This simulation includes a single homogeneously mixed host population of 40M hosts, contact rate β = 0.6 and a 4 day recovery rate. Each epitope unencountered by the host contributes to a 17.5% increase in the risk of infection with a different strain (see <xref ref-type="sec" rid="s4">methods</xref>
for full description of epidemiological parameters).</p>
</caption>
<graphic xlink:href="ppat.1003104.g003"></graphic>
</fig>
<p>A mutation-free model can result in different alternative dynamics based on model parameters ranging from stable coexistence of completely discordant antigenic-phenotypes to the successive replacement of strains through chaotic or cyclic behavior <xref rid="ppat.1003104-Gupta1" ref-type="bibr">[25]</xref>
. Not surprisingly, our model implementation with no explicit evolution also generates these waves of replacement (<xref ref-type="fig" rid="ppat-1003104-g003">Figure 3</xref>
), suggestive of H3N2 influenza as proposed by Recker et al. <xref rid="ppat.1003104-Recker1" ref-type="bibr">[22]</xref>
.</p>
</sec>
<sec id="s2b"><title>Phylodynamic Patterns across Different Mutation Rates</title>
<p>We examined the effects of mutation at different rates on the resulting phylodynamic patterns of the virus by seeding the population with a single strain and tracking antigenic and evolutionary changes. We measure diversity π as the average time separating two randomly selected contemporaneous viruses since their divergence from a common ancestor. Because branch lengths in our genealogies are measured in years, the resulting diversity is also measured in years.</p>
<p>In the absence of antigenic mutation, only a single strain persists, experiencing transient oscillatory dynamics between near extinction, and endemic equilibrium conditions (<xref ref-type="fig" rid="ppat-1003104-g004">Figure 4A</xref>
). As all viral traits are equal, there are no selective forces and the observed phylogeny and coalescence rates can be directly related to prevalence and incidence <xref rid="ppat.1003104-Volz1" ref-type="bibr">[29]</xref>
. This yields random coalescence within contemporaneous viral lineages and a large associated pairwise genetic diversity (π = 30±12 years; mean±std across 5 simulations).</p>
<fig id="ppat-1003104-g004" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1003104.g004</object-id>
<label>Figure 4</label>
<caption><title>Changes in the proportions of hosts that are infectious with different strains and the related phylogenetic behavior with increasing mutation rate.</title>
<p>Phylogenetic trees are based on samples of directly measured virus genealogy in the simulation, and only the last 40 years are visualized in the figure (for the complete genealogy over the whole time period see ). Diversity π is calculated as the mean distance, measured in years, for the coalescence of random pairs of contemporaneous samples in a tree. The MK related index (MKR) is calculated as the ratio of the antigenic mutation rate on the trunk (fixed) versus the antigenic mutation rate on the sidebranches. <bold>(A)</bold>
Model with no mutation, a single antigenic type persists under neutral evolution. When there are no antigenic mutations a genealogical tree which follows neutral viral evolution exists. Genetic diversity for this tree relates to population dynamics only – to incidence and the prevalence. <bold>(B)</bold>
Model with low mutation rate of ξ = 7.5×10<sup>−6</sup>
antigenic-mutations per day. Successive strain replacement with higher epidemic peaks is observed. Rare antigenic mutations are advantageous and are more likely to fix and have viable offsprings, consequently lowering genetic diversity. <bold>(C)</bold>
The introduction of a higher mutation rate ξ = 7.5×10<sup>−5</sup>
leads to antigenic and genetic divergence. Dynamics are ruled by endemic or cyclic behavior of discordant antigenic strains. Mutations are more likely to be deleterious, facing competition from the two prevalent strains. Phylogenetic patterns include two deep branches representing each strain and a low rate of coalescence between strains. <bold>(D)</bold>
For a mutation rate ξ = 7.5×10<sup>−4</sup>
epidemiological behavior resembles the evolution free framework (<xref ref-type="fig" rid="ppat-1003104-g003">Figure 3</xref>
). Phylogenetic patterns exhibit high genetic diversity and weak negative selection pressure. <bold>(E)</bold>
Loss of strain structure due to high mutation rate ξ = 7.5×10<sup>−3</sup>
. At this high mutation rate the antigenic traits are no longer heritable and each linage displays a constantly varying antigenic phenotype. No selection forces are measured and genetic diversity is expected to be determined by random coalescence. <bold>(F)</bold>
Summary statistics and typical trees for varying mutation rates and fixed crossimmunity (filled area within rectangles indicates 1σ confidence interval for 5 repeated runs). Simulation parameters are the same as those described in <xref ref-type="fig" rid="ppat-1003104-g003">Figure 3</xref>
, but include the possible extinction of strains, and mutations to individual epitopes at a specified rate ξ (see <xref ref-type="sec" rid="s4">methods</xref>
for full description of epidemiological parameters).</p>
</caption>
<graphic xlink:href="ppat.1003104.g004"></graphic>
</fig>
<p>For low mutation rates (<xref ref-type="fig" rid="ppat-1003104-g004">Figure 4B</xref>
) the introduction of new mutations is the critical determinant of strain dynamics. Each new variant outcompetes the one that came before, resulting in a spindly phylogenetic tree and therefore low diversity (π = 5.7±0.8 years). Temporally adjacent strains are antigenically similar, rather than discordant, forcing strong competitive exclusion and single strain dominance (ε = 0.93±0.02), where ε is the proportion of infections caused by the most common strain.</p>
<p>An increase in mutation rate (<xref ref-type="fig" rid="ppat-1003104-g004">Figure 4C</xref>
) leads to deeper branches with a corresponding increase in phylogenetic diversity (π = 220±100 years) and more pronounced antigenic divergence. Here, the population dynamics are ruled by the endemic or cyclic behavior of discordant antigenic sets. The emergence of new intermediate antigenic types is suppressed by competition from the two prevalent strains <xref rid="ppat.1003104-Gupta1" ref-type="bibr">[25]</xref>
.</p>
<p>At a relatively high mutation rate (<xref ref-type="fig" rid="ppat-1003104-g004">Figure 4D</xref>
), we approach population dynamics similar in appearance to those of the mutation-free model (<xref ref-type="fig" rid="ppat-1003104-g003">Figure 3</xref>
). Diversity is high (π = 120±30 years), with deep yet occasionally coalescing branches.</p>
<p>In general, a threshold exists at which mutation overwhelms selection resulting in a population drifting away from the fittest genotype <xref rid="ppat.1003104-Nowak2" ref-type="bibr">[30]</xref>
. For sufficiently high rates of antigenic mutation, all antigenic types reach near equal frequency in the population (<xref ref-type="fig" rid="ppat-1003104-g004">Figure 4E</xref>
). On a population scale, the high mutation rate weakens frequency-dependent selection and results in the breakage of antigenic strain structure; antigenic types do not cluster across the genealogy. The loss of selection forces breaks down phylogenetic structure and leads to a reduction in the depth of the branches (π = 120±70 years) compared to the one observed for discordant antigenic types.</p>
</sec>
<sec id="s2c"><title>Measuring Selection Strength and Direction</title>
<p>By comparing fixation versus extinction of antigenic mutations, using a quantity related to the McDonald-Kreitman (MK) index <xref rid="ppat.1003104-McDonald1" ref-type="bibr">[31]</xref>
, <xref rid="ppat.1003104-Bush2" ref-type="bibr">[32]</xref>
, <xref rid="ppat.1003104-Bhatt1" ref-type="bibr">[33]</xref>
we estimated the strength and direction of selection on antigenic mutations in our model (see <xref ref-type="sec" rid="s4">Methods</xref>
). Here we calculate an MK related (MKR) index as the ratio of the per-year rate of antigenic mutation on the trunk to the per-year rate of antigenic mutation on the side branches.</p>
<p>If antigenic mutations are advantageous for long-term virus persistence, an MKR ratio above 1 is expected. In this case, individuals exhibiting these antigenic mutations will be more likely to fix in the population and contribute to substitutions on the trunk of the phylogeny. Similarly, if antigenic mutations are deleterious to the long-term success of the virus, an MKR index of less than 1 is expected. This is because mutant individuals will tend to be lost from the population and side branches will show an excess rate of substitution.</p>
<p>We find that, when rare, antigenic mutations show highly increased rates of fixation (MKR = 19±11), and therefore evidence of strong positive selection (<xref ref-type="fig" rid="ppat-1003104-g004">Figure 4B</xref>
). Hence, we find that strong positive selection results in both a spindly tree and an overabundance of antigenic mutations of the trunk of the phylogeny. An increase in the mutation rate leads to the emergence of antigenically discordant types, and the suppression of other antigenic mutants (<xref ref-type="fig" rid="ppat-1003104-g004">Figure 4C</xref>
); here, we find strong negative selection mediated by host immunity with an MKR index of 0.47±0.22. At a still higher mutation rate (<xref ref-type="fig" rid="ppat-1003104-g004">Figure 4D</xref>
), we observe a balance of positive and negative selection resulting in MKR = 1.1±0.3. With saturating mutation rates (<xref ref-type="fig" rid="ppat-1003104-g004">Figure 4E</xref>
) we further lose the signature of selection (MKR = 1.0±0.1) on phylodynamic patterns.</p>
</sec>
<sec id="s2d"><title>Relationship between Cross-immunity, Mutation and Selection</title>
<p>Our model contains a cross-immunity parameter σ which allows us to explore a range of immune selection regimes: when σ = 1, we have full cross-protection (as might arise if each epitope elicited a very strong immune response) and when σ = 0, cross-protection between strains is only high if they share their entire variable repertoire.</p>
<p>In general, stronger cross-immunity results in lower prevalence as hosts fail to be re-infected (<xref ref-type="supplementary-material" rid="ppat.1003104.s002">Figure S1</xref>
). We find that, for most of the parameter space, genetic (genealogical) diversity π, increases with weaker cross-immunity and with more rapid mutation (<xref ref-type="fig" rid="ppat-1003104-g005">Figure 5A</xref>
). The (mostly) monotonic relationship between competition and diversity is broken at the threshold of limiting similarity <xref rid="ppat.1003104-Abrams1" ref-type="bibr">[34]</xref>
where, regardless of epitope differences, two strains suffer full cross-protection. This scenario, shown as a band on the right-hand side of <xref ref-type="fig" rid="ppat-1003104-g005">Figure 5A</xref>
where σ = 1, results in the disappearance of selective effects and greater levels of genetic diversity. Here, diversity rebounds to its neutral expectation due to random coalescence. Exceptions to the monotonic pattern of diversity with competition can also be found for intermediate mutation rates.</p>
<fig id="ppat-1003104-g005" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1003104.g005</object-id>
<label>Figure 5</label>
<caption><title>Changes in genetic diversity and the McDonald-Kreitman related index (MKR) for varying strengths of strain competition and antigenic mutation rates.</title>
<p><bold>(A)</bold>
Mean pairwise genetic diversity π is measured as the mean distance, measured in years, for the coalescence of random pairs of contemporaneous samples in a tree. Diversity measurement is capped by twice the simulation run length which amounts to 240 years. <bold>(B)</bold>
The MKR is measured as the ratio between the trunk antigenic mutation rate (fixed) and the sidebranches antigenic mutation rate. Evidence of positive selection is observed when the MKR index is significantly above one, and negative selection when it is significantly below one. Areas of strong positive selection are associated with lower genetic diversity as a small subset of the population contributes to long term viral evolution. Strong negative selection is associated with disruptive selection maintained by existing strains. <bold>(C)</bold>
Diversity for varying strengths and directions of selection as measured by the MKR index. Diversity decreases with stronger positive selection <italic>ρ</italic>
= −0.85 (Pearson's correlation right of dotted line), and increases for stronger negative selection <italic>ρ</italic>
= −0.28. The harmonic mean of the prevalence is also strongly correlated with genetic diversity <italic>ρ</italic>
= 0.76 (heat map). <bold>(D)</bold>
Typical trees for varying strengths of strain competition and antigenic mutation rate. Effective competition combined with a limited availability of antigenic mutations results in narrower trees with lower pairwise diversity. This figure was parameterized to use R0 = 3 and a population size of 50M to limit stochastic extinctions for a large parameter range (see <xref ref-type="sec" rid="s4">methods</xref>
for full description of epidemiological parameters).</p>
</caption>
<graphic xlink:href="ppat.1003104.g005"></graphic>
</fig>
<p>The relationship between mutation, cross-immunity, and the MKR index is less straight-forward (<xref ref-type="fig" rid="ppat-1003104-g005">Figure 5B</xref>
). Here, the highest levels of positive selection are present when cross-immunity is strong (σ = 0.8–0.9), and mutation is weak (ξ = 10<sup>−5</sup>
). When mutation rate is limiting (ξ<5×10<sup>−5</sup>
) then antigenic mutations are favored by natural selection (MKR>1). However, when mutation rates are higher (5×10<sup>−5</sup>
<ξ<10<sup>−3</sup>
), negative selection by-and-large predominates. The strongest negative selection occurs in a region of moderate cross-immunity (σ = 0.6) corresponding to previously observed discordant dynamics (<xref ref-type="fig" rid="ppat-1003104-g004">Figure 4C</xref>
).</p>
<p>There is also a clear relationship between diversity and selection as measured by the MKR index (<xref ref-type="fig" rid="ppat-1003104-g005">Figure 5C</xref>
). We observe a strong negative correlation between MKR and levels of diversity (<italic>ρ</italic>
= −0.86; Pearson's correlation). If we separate results into a regime of positive selection (MKR>1) and a regime of negative selection (MKR<1), we observe similar results within each regime. Stronger positive selection coincides with a decrease in genetic diversity (<italic>ρ</italic>
= −0.85 when MKR>1), and stronger negative selection tends further increase diversity through the persistence of discordant strains and associated deep branches (<italic>ρ</italic>
= −0.28 when MKR<1). As expected from population genetic theory <xref rid="ppat.1003104-Volz1" ref-type="bibr">[29]</xref>
, increases in viral prevalence also coincide with increases in viral diversity, however, the correlation is weaker under positive selection (<italic>ρ</italic>
= 0.74 when MKR>1) than under negative selection (<italic>ρ</italic>
= 0.87 when MKR<1) and cannot be trivially dissociated from the effects of selection.</p>
<p>Two additional strain diversity measurements based on the ecological dynamics are shown in <xref ref-type="supplementary-material" rid="ppat.1003104.s003">Figure S2</xref>
, the Shannon diversity index and the level of single strain dominance (<xref ref-type="sec" rid="s4">Methods</xref>
). Similar to genetic diversity, positive selection is correlated with an increase in single strain dominance (<italic>ρ</italic>
= 0.88 when MKR>1) and a decrease in Shannon diversity (<italic>ρ</italic>
= −0.88 when MKR>1). Negative selection decreases Shannon diversity (<italic>ρ</italic>
= 0.39 when MKR<1) and increases single strain dominance (<italic>ρ</italic>
= −0.23 when MKR<1). While negative selection lowers the number of circulating strains, it increases genetic diversity π through the existence of deep non coalescing branches.</p>
</sec>
<sec id="s2e"><title>H3N2-like Characteristics</title>
<p>The patterns described so far suggest that the dynamics of H3N2 influenza within this framework correspond to a regime in which host immune mediated selection is strong and the antigenic mutation rate is low. We now extend the model in order to examine other characteristics relevant to H3N2 in a more detailed epidemiological setting that includes seasonality and a basic global population structure.</p>
<p>In this analysis we include three demes representing the northern hemisphere, the southern hemisphere and the tropics. Northern and southern hemisphere demes experience an opposing seasonal modulation (with a 14% amplitude and six months phase difference) while tropical regions experience two weaker seasons annually <xref rid="ppat.1003104-Chew1" ref-type="bibr">[35]</xref>
(see <xref ref-type="supplementary-material" rid="ppat.1003104.s004">Figure S3</xref>
and <xref ref-type="sec" rid="s4">Methods</xref>
). In addition the southern hemisphere population is reduced in comparison to northern hemisphere and tropical populations (<xref ref-type="sec" rid="s4">Methods</xref>
).</p>
<p>In this model we use an antigenic repertoire with 4 epitopes differing in the number of alternative variants per epitope. A typical tree for this configuration together with the corresponding diversity skyline is depicted in <xref ref-type="fig" rid="ppat-1003104-g006">Figure 6A</xref>
. We observe 13±6 antigenic clusters that come to dominate the virus population over the course of the 40 year simulation (<xref ref-type="supplementary-material" rid="ppat.1003104.s005">Figure S4</xref>
-A) with an average duration of 4±2 years. Clusters are defined based on cumulative changes in two or more epitopes based on <xref rid="ppat.1003104-Huang1" ref-type="bibr">[36]</xref>
(<xref ref-type="sec" rid="s4">Methods</xref>
). The turnover of virus strains results in a characteristic spindly phylogenetic tree and low standing genetic diversity (π = 5.7±0.1 years). Over the course of the 40-year timespan, genetic diversity experiences a boom and bust pattern (<xref ref-type="fig" rid="ppat-1003104-g006">Figure 6A</xref>
) with a 10%–90% range of 3–9.5 years measured by combining diversity skylines of five repeated simulations. The repeat of exact antigenic types is uncommon in the model (<xref ref-type="supplementary-material" rid="ppat.1003104.s006">Figure S5</xref>
) while epitopes with more restricted variability (2–3 variants) frequently reemerge (<xref ref-type="supplementary-material" rid="ppat.1003104.s007">Figure S6</xref>
).</p>
<fig id="ppat-1003104-g006" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1003104.g006</object-id>
<label>Figure 6</label>
<caption><title>H3N2 like characteristics.</title>
<p>Using a model with an alternative epitope configuration and three global demes representing the northern hemisphere N, southern hemisphere S, and the tropics T. <bold>(A)</bold>
Top - typical tree colored by antigenic clusters (see <xref ref-type="sec" rid="s4">Methods</xref>
). On average 13±6 clusters lasting 4±2 years come to dominate the virus population over the course of the 40 year simulation. Bottom - Genealogical diversity (π) displaying “boom & bust” patterns (10%–90% range of 3–9.5 years) associated with H3 diversity with an average of 5.7±0.1 years. <bold>(B)</bold>
Top - typical tree colored by deme. The tropics metapopulation has a higher proportion in establishing the trunk (68%±9) of the phylogeny followed by the northern (20%±8) and the southern (12%±3) population. Bottom – proportion of hosts infected in the northern hemisphere, tropics and southern hemisphere. Epitope configuration for this figure was 5×4×3×2 variants per epitope, 4 epitope system. Average yearly incidence in the northern and southern hemisphere demes is (5.7%±0.1, 5.8%±0.1) respectively, while incidence in the tropics is slightly lower 5.5%±0.1. Annual epidemics are generated almost regularly yet display a high level of variability in peak size in the northern hemisphere (CoV = 1.1±0.1, coefficient of variation) and lower variability in the tropics (CoV = 0.7±0.1). In this simulation 42M hosts were divided to three demes with the tropics and the south having 16M hosts and the southern hemisphere having a lower 10M population. Annual seasonal patterns were established for the temperate demes, and biannual weaker seasonality in the tropics (see <xref ref-type="sec" rid="s4">methods</xref>
for full description of epidemiological parameters).</p>
</caption>
<graphic xlink:href="ppat.1003104.g006"></graphic>
</fig>
<p>Average yearly incidence in the northern and southern hemisphere demes is (5.7%±0.1, 5.8%±0.1) respectively (<xref ref-type="fig" rid="ppat-1003104-g006">Figure 6B</xref>
), while incidence in the tropics is slightly lower 5.5%±0.1. Annual epidemics are generated almost regularly yet display a high level of variability in peak size in the northern hemisphere (CoV = 1.1±0.1, coefficient of variation) and lower variability in the tropics (CoV = 0.7±0.1). The interquartile range in peak weekly cases ranges (IQR≈200–800 cases per 100000) in the northern hemisphere and (IQR≈400–1100) in the tropics.</p>
<p>In this model the tropics or lower and mixed seasonality populations exhibit a greater role (68%±9) in establishing the trunk of the influenza tree (<xref ref-type="fig" rid="ppat-1003104-g006">Figure 6B</xref>
). The southern hemisphere experiences a smaller (12%±3) part in establishing the trunk of the tree in comparison to the northern hemisphere (20%±8). In addition we find that antigenic variants are more likely to reach significant prevalence in the tropics earlier making the tropics “antigenically ahead” (<xref ref-type="supplementary-material" rid="ppat.1003104.s005">Figure S4</xref>
-B). Antigenic variants reach 5% of their total deme prevalence 2±1.5 months earlier in the tropics compared to the northern hemisphere and 3±2 months earlier in the tropics compared to the southern hemisphere (p<0.001 for the combined results). Similarly antigenic variants decline (reach 95% of their total prevalence) earlier (1.7±0.3 month, p<0.0005) in the tropics compared to the northern hemisphere, yet not significantly earlier or later than the southern hemisphere.</p>
<p>We find that strong competition and high R0 values generates more regular annual epidemic peaks while maintaining low genetic diversity. In addition we find that within antigenic cluster evolution also contributes to maintain low genetic diversity. An increase (0.005 compared to 0.001 of contacts) in the strength of the metapopulation coupling slightly improves the epidemiology by decreasing the likelihood of long periods without annual epidemics.</p>
<p>To establish whether low genetic diversity can be maintained when the number of epitopes or variants per epitope is increased we repeated the same parameterization with double the number of epitopes and with twice the number of variants per epitope and by keeping either the per-epitope or overall mutation rate. The results are summarized in <xref ref-type="supplementary-material" rid="ppat.1003104.s009">Table S1</xref>
. When doubling the number of epitopes, the model can attain similar results with respect to genetic diversity π and overall incidence, when the total mutation rate across epitopes is maintained and the cross-immunity decrease per–epitope change is halved. In contrast, it is not clear whether a model that includes an increase in the number of variants per epitope can maintain low genetic diversity levels and maintain similar or higher incidence levels.</p>
</sec>
</sec>
<sec id="s3"><title>Discussion</title>
<p>Herein, we implemented an individual-based model that allowed us to track both the ecological and evolutionary dynamics of a pathogen population, in which cross-immunity is orchestrated by a finite set of antigenic loci of limited variability <xref rid="ppat.1003104-Recker1" ref-type="bibr">[22]</xref>
. We used this model to compare phylodynamic patterns under a regime governed primarily by limitation on the introduction of antigenic mutations (mutation limited), to a regime determined by the availability of antigenic niches (selection limited), and under varying strengths of competition between strains. We use this framework to determine the conditions under which a limited antigenic repertoire could explain the observed phylodynamic patterns of H3N2 influenza.</p>
<p>Explicit modeling of evolution, through the introduction of antigenic mutation at different rates, allows us to consider phylogenetic trees in addition to epidemiological dynamics. Resulting phylodynamic patterns range from successive strain turnover, to discordant antigenic sets, to dynamics resembling those of a model lacking explicit evolution and finally to the collapse of antigenic structure. Each of these can be explained by the interplay of selection and mutation, as measured here through the MKR index, and by considering different strengths of immunity generating competition between strains.</p>
<p>The dynamics of our individual-based model are generally in good agreement with the epidemic behavior of influenza A (<xref ref-type="fig" rid="ppat-1003104-g006">Figure 6</xref>
). Like observed epidemiological patterns <xref rid="ppat.1003104-Koelle1" ref-type="bibr">[17]</xref>
, <xref rid="ppat.1003104-Russell1" ref-type="bibr">[23]</xref>
, <xref rid="ppat.1003104-Fleming1" ref-type="bibr">[37]</xref>
, <xref rid="ppat.1003104-Monto1" ref-type="bibr">[38]</xref>
, <xref rid="ppat.1003104-Karpova1" ref-type="bibr">[39]</xref>
, <xref rid="ppat.1003104-SaltyteBenth1" ref-type="bibr">[40]</xref>
, <xref rid="ppat.1003104-Sumi1" ref-type="bibr">[41]</xref>
, annual temperate climate epidemics occur almost regularly with substantial year-to-year variation in incidence (CoV = 1.1±0.1 compared to (CoV = 1.0±0.2) in literature survey. Observed temperate climate annual attack rates of influenza A (H3N2) are slightly higher, approximately 8% from 1976 to 1981 <xref rid="ppat.1003104-Monto1" ref-type="bibr">[38]</xref>
compared to 5.8%±0.1 in simulation, while peak epidemic weakly cases are higher in the simulation (IQR≈200–800 cases per 100000) in comparison to (IQR≈130–380, IQR≈80–240) in <xref rid="ppat.1003104-Fleming1" ref-type="bibr">[37]</xref>
and <xref rid="ppat.1003104-SaltyteBenth1" ref-type="bibr">[40]</xref>
respectively. The tropics exhibit lower and weaker seasonality (<xref ref-type="fig" rid="ppat-1003104-g006">Figure 6</xref>
, <xref ref-type="supplementary-material" rid="ppat.1003104.s004">Figure S3</xref>
) with slightly lower yearly attack rates (5.5%±0.1) and substantially lower prevalence (<xref ref-type="fig" rid="ppat-1003104-g006">Figure 6</xref>
). In agreement with antigenic cartography <xref rid="ppat.1003104-Smith1" ref-type="bibr">[18]</xref>
, <xref rid="ppat.1003104-Du1" ref-type="bibr">[42]</xref>
13±6 clusters dominate the global world population (<xref ref-type="supplementary-material" rid="ppat.1003104.s005">Figure S4</xref>
-A) with an average duration of 4±2 years, exhibiting mostly the dominance of 1–2 clusters globally. With respect to individual epitope changes we find the model reproduces the observation of the tropics being “antigenically ahead” <xref rid="ppat.1003104-Russell1" ref-type="bibr">[23]</xref>
, giving rise to antigenic changes 2.5±1.5 month ahead of the northern and southern hemisphere (<xref ref-type="supplementary-material" rid="ppat.1003104.s005">Figure S4</xref>
-B) and showing decline in antigenic variants 1.7±0.2 month earlier than the northern hemisphere. In agreement with observed phylodynamic patterns <xref rid="ppat.1003104-Bedford3" ref-type="bibr">[43]</xref>
the tropics metapopulation has a higher proportion in establishing the trunk (68%±9) of the phylogeny followed by the northern (20%±8) and the southern (12%±3) population. The higher contribution of East and South-East Asia as the origin of H3N2 globally circulating lineages is hypothesized to originate from lower and mixed seasonality in these regions and is consistent with our model <xref rid="ppat.1003104-Russell1" ref-type="bibr">[23]</xref>
. The key difference between the hemispheres in the model being, lower population size in the southern hemisphere with proportionally lower contact rate between the meta-populations.</p>
<p>Refinement of the epidemiological model, such as the inclusion of an exposed period, can further improve the comparison to empirical data. In particular, the above properties were obtained with a basic reproduction number of R0≈3.24, on the upper bounds of current estimates for seasonal influenza. This value can possibly be decreased by considering such an extension.</p>
<p>We find that a model with 4 epitopes and a low but variable number of variants per site, an antigenic mutation rate of ≈10<sup>−5</sup>
per day and a reduction of cross-immunity of 13% per epitope results in phylodynamic patterns broadly consistent with those seen in H3N2 influenza (<xref ref-type="fig" rid="ppat-1003104-g006">Figure 6</xref>
). When doubling the number of epitopes, the model maintained similar results with respect to genetic diversity π and overall incidence, when the total mutation rate across epitopes was maintained and the cross-immunity was modified to a 6.5% per epitope change. These parameters are quite comparable to parameters used in other models of influenza evolution. For example, Koelle et al. <xref rid="ppat.1003104-Smith1" ref-type="bibr">[18]</xref>
use 5 epitopes with mutations of either large or small antigenic effect. Small mutations reduce cross-immunity by 7% and occur at a rate of ∼5×10<sup>−4</sup>
per day, while large mutations reduce cross-immunity by 20% and occur at a rate of ∼10<sup>−5</sup>
per day. In the model of Bedford et al. <xref rid="ppat.1003104-Bedford2" ref-type="bibr">[21]</xref>
mutations reduce cross-immunity by between 1% and 11% (95% bounds), but occur at a faster rate of 10<sup>−4</sup>
per day. Ferguson et al. <xref rid="ppat.1003104-Koelle1" ref-type="bibr">[17]</xref>
find that a model with 12 codons, each with 20 amino acid variants, in which mutations occur at a rate of 3×10<sup>−5</sup>
per day and reduce cross-immunity by ∼7% gives restricted diversity without short-term strain-transcending immunity, and 1.2×10<sup>−4</sup>
per day, when transient immunity is included. From this, it seems clear that models involving a slow influx of antigenic mutants of around 10<sup>−5</sup>
per day are generally capable of producing influenza-like patterns of restricted diversity.</p>
<p>Increasing host population size in the model results in an increase in viral genetic diversity, as more opportunities for antigenic mutation arise within the larger host population. Thus, scaling competitive interactions between strains, and/or antigenic mutation rate, is required to maintain limitations on the effective exploration of antigenic space. In addition, other epidemiological phenomena, besides low antigenic mutation rates, may also contribute to limit the rate at which novel antigenic phenotypes emerge in the influenza population. These may be provided by population structure and the seasonality of transmission <xref rid="ppat.1003104-Russell1" ref-type="bibr">[23]</xref>
, <xref rid="ppat.1003104-Creanza1" ref-type="bibr">[44]</xref>
, as well as by short-term strain-transcending immunity, which was found capable of limiting genetic and antigenic diversity in a similar model with a much larger antigenic space <xref rid="ppat.1003104-Ferguson1" ref-type="bibr">[16]</xref>
and in a limited diversity antigenic model <xref rid="ppat.1003104-Minayev1" ref-type="bibr">[45]</xref>
. However, a global metapopulation structure is not expected to be the dominant cause behind the low standing genetic diversity of influenza. Influenza B exhibits similar epidemiological dynamics, and lower prevalence, yet it exhibits much higher genetic diversity through the co-circulation of multiple lineages <xref rid="ppat.1003104-Ferguson1" ref-type="bibr">[16]</xref>
, <xref rid="ppat.1003104-Yamashita1" ref-type="bibr">[46]</xref>
. Also, a more complex metapopulation structure with multiple patches can either increase genetic diversity by facilitating the coexistence of viruses at different weakly coupled patches, or decrease genetic diversity through the generation of population bottlenecks. The role of variation in viral fitness is an important consideration in future studies, particularly in light of recent observations linking binding properties of HA with antigenic escape <xref rid="ppat.1003104-Hensley1" ref-type="bibr">[47]</xref>
. The empirical finding of a non-trivial relationship between virus fitness in susceptible individuals and immune evasion was suggested as a possible alternative mechanism for generating positive selection pressure on antigenic sites and for limiting antigenic diversity <xref rid="ppat.1003104-Hensley1" ref-type="bibr">[47]</xref>
.</p>
<p>Future work should investigate quantitative patterns and statistical approaches for discriminating among the different models and associated hypotheses that currently exist in the literature and for inferring the relative importance of the mechanisms they represent, keeping in mind that the models are not necessarily mutually exclusive. At the same time, empirical advances on the molecular basis of immune evasion and recognition, on the genotype-to-phenotype map, and on epitope identification and population serology, will allow a better evaluation of the models' assumptions, including the representation of serological space.</p>
<p>In common with the Recker et al. <xref rid="ppat.1003104-Recker1" ref-type="bibr">[22]</xref>
model and in contrast with other phylodynamic models <xref rid="ppat.1003104-Ferguson1" ref-type="bibr">[16]</xref>
, <xref rid="ppat.1003104-Koelle1" ref-type="bibr">[17]</xref>
, <xref rid="ppat.1003104-Bedford2" ref-type="bibr">[21]</xref>
, we find here that antigenic epitopes are frequently recycled (<xref ref-type="supplementary-material" rid="ppat.1003104.s007">Figure S6</xref>
). Importantly, this does not mean that such recycling is observed for the antigenic types (epitope repertoires) themselves, since the same antigenic type only re-emerges at long intervals (<xref ref-type="supplementary-material" rid="ppat.1003104.s006">Figure S5</xref>
) and rarely in the course of 40 simulated years. It's possible that such reemergence could explain the antigenic cross-reactivity between sera from around the 1918 H1N1 pandemic and viruses emerging in the 2009 H1N1 pandemic <xref rid="ppat.1003104-Hancock1" ref-type="bibr">[48]</xref>
, <xref rid="ppat.1003104-Itoh1" ref-type="bibr">[49]</xref>
, <xref rid="ppat.1003104-Bandaranayake1" ref-type="bibr">[50]</xref>
, <xref rid="ppat.1003104-Chi1" ref-type="bibr">[51]</xref>
. However, antigenic stasis of the swine lineage leading to the 2009 pandemic could also explain these observations. Much further work on epitope identification and population-wide serological surveys is necessary to establish the validity of this model's prediction on the re-cycling of constituent low diversity epitope variants (<xref ref-type="supplementary-material" rid="ppat.1003104.s007">Figure S6</xref>
). Nevertheless, several empirical observations are becoming available that are consistent with such recycling and the subject is discussed in detail in the companion paper <xref rid="ppat.1003104-Wikramaratna1" ref-type="bibr">[52]</xref>
. For example, an antigenic analyses performed on H2N2 influenza, a number of monoclonal antibodies raised against a 1957 strain were shown to cross-react strongly with a strain isolated in 1964, yet not with the 1963 strain <xref rid="ppat.1003104-Tsuchiya1" ref-type="bibr">[53]</xref>
. In Reichert et al. <xref rid="ppat.1003104-Reichert1" ref-type="bibr">[54]</xref>
the hemagglutinin of both novel pandemic H1N1 and pre-1940 H1N1 lack specific glycosylation sites on the globular head of HA1. These reverse glycosylation patterns were suggested to possibly shield antigenic sites for a timescale of decades and in so doing, to effectively contribute to their recycling and to the age distribution of cases. In Bui et al. <xref rid="ppat.1003104-Bui1" ref-type="bibr">[55]</xref>
, several protective antigenic and T cell H3 epitopes show temporal variability across drift variants, with two of these specifically exhibiting a decrease and increase in conservancy consistent with epitope “recycling”. Post translational and conformational changes may hinder the validity of this analysis especially for epitope 1 which acquired two surrounding glycosylation sites. In Wang et al. <xref rid="ppat.1003104-Wang1" ref-type="bibr">[56]</xref>
, mice Anti-H3 mAbs were shown to neutralize H3 viruses that span 40 years, as measured by immunofluorescence against MDCK cells (see <xref ref-type="table" rid="ppat-1003104-t002">Table 2</xref>
in <xref rid="ppat.1003104-Wang1" ref-type="bibr">[56]</xref>
). All three mAbs (see <xref ref-type="fig" rid="ppat-1003104-g004">Figure 4</xref>
in <xref rid="ppat.1003104-Wang1" ref-type="bibr">[56]</xref>
) displayed variability in their ability to neutralize H3 viruses for lower concentrations (<15 µg/ml of 7A7 and <25 µg/ml for the other two) in plaque reduction assays. For example for mAb 7A7 neutralization was better for HK68, than diminished for BJ92 and then increases for PAN99 and BRIS07. This pattern could also be due to secondary effects of amino acid differences outside the actual epitope as well through structural effects, but effectively behaves as epitope recycling over substantial durations of many years.</p>
<p>In conclusion, within our framework, the rate of antigenic mutation was found to strongly influence whether selection was positive or negative, and hence, the topology of the tree and associated diversity of the virus. Strong positive selection is generated by effective competition under low mutation rates, and results in spindly trees with low genetic diversity. In this regime, antigenic mutations often fix in the virus population, lowering genetic diversity, as consistent with H3N2. An increase in mutation rate across a broad spectrum in competition strength, leads to negative selection and generates antigenic divergence. This can potentially result in the coexistence of discordant antigenic types repressing the emergence of antigenic hybrids, through strong negative selection on antigenic change, with each discordant antigenic type maintaining a deep phylogenetic branch. Although not strictly mirroring the assumptions about development of the Recker et al model, our framework strongly implies that limitations on antigenic architecture alone are unlikely to reliably reproduce “skinny' trees and some restrictions on mutation rate and/or other considerations such as fitness differences are likely to play a role. It is important to note that this exercise does not also privilege other hypotheses concerning diversity restriction in influenza as these also are strongly sensitive to mutation rate. Overall, it emphasizes that phylogenetic patterns do not serve as a discriminatory tool between these by no means mutually exclusive hypotheses. However, they can provide a basis to exclude specific hypotheses and offer a means by which the contributions of mutation and selection can be assessed. Needless to say, the latter has important implications for the updating of vaccines against influenza. Under a mutation-limited regime, a hypothetical vaccine should be effective until a new antigenic variant is introduced to the population through migration or mutation. Alternatively, when a limited number of alternative but conserved epitopes are continuously circulating with their abundance determined through competitive interactions and immune mediated selection, a vaccine against one of them may lead to rapid strain replacement, while a vaccine against all of them may result in effective intervention.</p>
</sec>
<sec sec-type="materials|methods" id="s4"><title>Materials and Methods</title>
<sec id="s4a"><title>Cross Immunity</title>
<p>We assume that a strain's antigenic attributes are determined by a set of separate epitopes and that each epitope contains a discrete number of alternative variants. Thus, a strain's antigenic properties are defined by an <italic>n</italic>
-tuple with <italic>k<sub>i</sub>
</italic>
variants per epitope giving a possible number of <inline-formula><inline-graphic xlink:href="ppat.1003104.e002.jpg"></inline-graphic>
</inline-formula>
antigenically distinct strains. Hosts acquire immunity to viral epitopes following infection. Fully naïve hosts are always infected following contact at rate β. The risk of reinfection with the exact same strain is always zero. The chance of reinfection with a different strain is based on the similarity to previously encountered strains, measured through the fraction of previously encountered epitopes: <inline-formula><inline-graphic xlink:href="ppat.1003104.e003.jpg"></inline-graphic>
</inline-formula>
and is at most 100%. Where <italic>f</italic>
is the fraction of previously encountered epitopes and σ is the strength of crossimmunity. Lower σ values correspond to weaker competition between strains. A form of generalized immunity is attained for σ>0.8 in the five epitope case, relating to a reduced risk of reinfection following previous exposure to any strain (<xref ref-type="fig" rid="ppat-1003104-g002">Figure 2B</xref>
).</p>
</sec>
<sec id="s4b"><title>Mutations</title>
<p>Phenotype changes are driven by mutation events. Mutations change the antigenic properties of a strain but do not influence the shape of the genealogical tree directly; the tree shape will be determined by selection, epidemiological dynamics and the stochastic processes favoring a specific isolate and its offsprings, implicitly rather than explicitly. Mutations involve changes in a single epitope; this makes some phenotypic changes more attainable in comparison to others, even for high mutation rates. The antigenic mutation rate ξ gives the per-day probability for a virus to mutate in a single epitope site.</p>
</sec>
<sec id="s4c"><title>Genealogy Tracking and Related Diversity Quantities</title>
<p>The genealogy of the virus population was tracked directly throughout the simulation (<xref ref-type="fig" rid="ppat-1003104-g002">Figure 2A</xref>
). Constant random sampling of viruses was performed periodically. Genealogical pairwise diversity (π) was measured by averaging the time unit distance on the tree between random contemporaneous sample pairs (<xref ref-type="fig" rid="ppat-1003104-g002">Figure 2C</xref>
). This quantity relates to pairwise genetic diversity, as measured on an accurately reconstructed phylogenetic tree. Diversity measures are limited in our simulation to a maximum total of twice the total running time of the simulation which amounts to 240 years. Although it's clear that some parameter ranges would show diversity greater than 240 years, they will saturate at this threshold (<xref ref-type="fig" rid="ppat-1003104-g005">Figure 5C</xref>
).</p>
<p>The MKR index was calculated by dividing the observed rate of occurrence of phenotypic mutations on the trunk of the tree, by the per-year mutation rate on the side branches (<xref ref-type="fig" rid="ppat-1003104-g002">Figure 2D</xref>
). This allows us to estimate the importance of antigenic mutations on the likelihood of fixation of a given viral linage.</p>
<p>Single strain dominance was calculated based on the quantity ε from <xref rid="ppat.1003104-Recker1" ref-type="bibr">[22]</xref>
and is calculated using the following formula:<disp-formula id="ppat.1003104.e004"><graphic xlink:href="ppat.1003104.e004.jpg" position="anchor" orientation="portrait"></graphic>
<label>(1)</label>
</disp-formula>
where y<sub>max</sub>
and y<sub>sub</sub>
are the prevalence of the most and 2<sup>nd</sup>
most prevalent strains. The normalized difference between the two is averaged across P epidemic peaks.</p>
</sec>
<sec id="s4d"><title>Increasing Mutation Rate Parameters</title>
<p>For <xref ref-type="fig" rid="ppat-1003104-g003">figures 3</xref>
, <xref ref-type="fig" rid="ppat-1003104-g004">4</xref>
and <xref ref-type="supplementary-material" rid="ppat.1003104.s008">S7</xref>
, R0 was set to ≈2.4. Antigenic diversity was limited to 2 variants per epitope to attain tractable results for a wide range of mutation and cross-immunity parameters, and for the same reason no metapopulation structure was established (<xref ref-type="table" rid="ppat-1003104-t001">Table 1</xref>
).</p>
<table-wrap id="ppat-1003104-t001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1003104.t001</object-id>
<label>Table 1</label>
<caption><title>Increasing mutation rate model parameters.</title>
</caption>
<alternatives><graphic id="ppat-1003104-t001-1" xlink:href="ppat.1003104.t001"></graphic>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead><tr><td align="left" rowspan="1" colspan="1">Parameter</td>
<td align="left" rowspan="1" colspan="1">Value</td>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1">population size N</td>
<td align="left" rowspan="1" colspan="1">40,000,000</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">contact rate β</td>
<td align="left" rowspan="1" colspan="1">0.6[1/day]</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">recovery rate ν</td>
<td align="left" rowspan="1" colspan="1">0.25[1/day]</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">birth/death rate μ</td>
<td align="left" rowspan="1" colspan="1">1/25[1/year]</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Epitopes</td>
<td align="left" rowspan="1" colspan="1">5</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">variants per epitope</td>
<td align="left" rowspan="1" colspan="1">2×2×2×2×2</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">cross-immunity σ</td>
<td align="left" rowspan="1" colspan="1">0.825</td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
</sec>
<sec id="s4e"><title>Mutation-Competition Simulation Parameters</title>
<p>For this set of simulations (<xref ref-type="fig" rid="ppat-1003104-g005">Figure 5</xref>
) R0 was set to ≈3 and no metapopulation structure was assumed. Higher R0 and longer duration of infection reduce the effect of critical community size over the range of parameters analyzed given the population size determined by computer resource use (<xref ref-type="table" rid="ppat-1003104-t002">Table 2</xref>
).</p>
<table-wrap id="ppat-1003104-t002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1003104.t002</object-id>
<label>Table 2</label>
<caption><title>Mutation-competition simulation parameters.</title>
</caption>
<alternatives><graphic id="ppat-1003104-t002-2" xlink:href="ppat.1003104.t002"></graphic>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead><tr><td align="left" rowspan="1" colspan="1">Parameter</td>
<td align="left" rowspan="1" colspan="1">Value</td>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1">population size N</td>
<td align="left" rowspan="1" colspan="1">50,000,000</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">contact rate β</td>
<td align="left" rowspan="1" colspan="1">0.6[1/day]</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">recovery rate ν</td>
<td align="left" rowspan="1" colspan="1">0.2[1/day]</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">birth/death rate μ</td>
<td align="left" rowspan="1" colspan="1">1/30[1/year]</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">epitopes</td>
<td align="left" rowspan="1" colspan="1">5</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">variants per epitope</td>
<td align="left" rowspan="1" colspan="1">2×2×2×2×2</td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
</sec>
<sec id="s4f"><title>H3N2 Simulation Parameters</title>
<p>This single simulation parameterization is intended to test whether H3N2-like phylogenetic trees can result from a model with a restricted set of antigenic phenotypes using an alternative epitope configuration and a basic metapopulation structure which includes seasonality. The basic reproduction number was set to R0≈3.24. In this case three demes were assumed, representing the northern hemisphere, the southern hemisphere and the tropics. Four epitopes with a variable number of variants per epitope were used. Contact rate was attenuated sinusoidally for southern and northern hemispheres to establish seasonal patterns. Tropical climate seasonality was sinusoidally modulated include two seasons <xref rid="ppat.1003104-Chew1" ref-type="bibr">[35]</xref>
of weaker amplitude <xref rid="ppat.1003104-Viboud1" ref-type="bibr">[57]</xref>
(<xref ref-type="table" rid="ppat-1003104-t003">Table 3</xref>
).</p>
<table-wrap id="ppat-1003104-t003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.ppat.1003104.t003</object-id>
<label>Table 3</label>
<caption><title>H3N2 simulation parameters.</title>
</caption>
<alternatives><graphic id="ppat-1003104-t003-3" xlink:href="ppat.1003104.t003"></graphic>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead><tr><td align="left" rowspan="1" colspan="1">Parameter</td>
<td align="left" rowspan="1" colspan="1">Value</td>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1">population size N North/Tropics/South</td>
<td align="left" rowspan="1" colspan="1">16M/16M/10M</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">contact rate β<sub>0</sub>
</td>
<td align="left" rowspan="1" colspan="1">0.6[1/day]</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Between deme contact</td>
<td align="left" rowspan="1" colspan="1">0.005[1/day]</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">recovery rate ν</td>
<td align="left" rowspan="1" colspan="1">0.185[1/day]</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">birth/death rate μ</td>
<td align="left" rowspan="1" colspan="1">1/30[1/year]</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">epitopes</td>
<td align="left" rowspan="1" colspan="1">4</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">variants per epitope</td>
<td align="left" rowspan="1" colspan="1">5×4×3×2</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Temperate Climate Seasonality North/South</td>
<td align="left" rowspan="1" colspan="1">14% Jan/July Six Month Phase</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Tropics Seasonality</td>
<td align="left" rowspan="1" colspan="1">7% Dec/June</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Mutation rate ξ</td>
<td align="left" rowspan="1" colspan="1">0.000008[1/day]</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">cross-immunity σ</td>
<td align="left" rowspan="1" colspan="1">0.87</td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
</sec>
<sec id="s4g"><title>Antigenic Cluster Transitions</title>
<p>It is not clear how changes in individual epitopes relate to the antigenic clusters as proposed by Smith et al. <xref rid="ppat.1003104-Smith1" ref-type="bibr">[18]</xref>
or Du et al. <xref rid="ppat.1003104-Du1" ref-type="bibr">[42]</xref>
. A new antigenic variant differing by 2 or more epitopes from any previous cluster strain was grouped in a new cluster in agreement with Huang et al. <xref rid="ppat.1003104-Huang1" ref-type="bibr">[36]</xref>
. This does not affect model dynamics, and relates only to the coloring of clusters in <xref ref-type="fig" rid="ppat-1003104-g006">Figure 6</xref>
and <xref ref-type="supplementary-material" rid="ppat.1003104.s005">Figure S4</xref>
-A.</p>
</sec>
</sec>
<sec sec-type="supplementary-material" id="s5"><title>Supporting Information</title>
<supplementary-material content-type="local-data" id="ppat.1003104.s001"><label>Accessions S1</label>
<caption><p><bold>GenBank accession numbers used in reconstruction of </bold>
<xref ref-type="fig" rid="ppat-1003104-g001"><bold>figure 1</bold>
</xref>
<bold> phylogenetic tree.</bold>
</p>
<p>(CSV)</p>
</caption>
<media xlink:href="ppat.1003104.s001.csv"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="ppat.1003104.s002"><label>Figure S1</label>
<caption><p><bold>Changes in prevalence for carrying strengths of strain competition and antigenic mutation rates.</bold>
Prevalence increases as crossimmunity between strains decreases, enabling multiple infections. When crossimmunity is 1, all strains are antigenically equal, and one lifetime infection is possible. When no crossimmunity is present, each antigenic-type can independently infect a host once.</p>
<p>(PDF)</p>
</caption>
<media xlink:href="ppat.1003104.s002.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="ppat.1003104.s003"><label>Figure S2</label>
<caption><p><bold>Changes in antigenic diversity and the McDonald-Kreitman related index (MKR) for varying strengths of strain competition and antigenic mutation rates.</bold>
<bold>(A)</bold>
Mean Shannon diversity was measured for 40 years of simulation. Shannon diversity ranges from zero when a single circulating antigenic variant is present at each time point, to approximately 3.5nats when all the possible antigenic variants are continuously present. For a large range of the parameter space, stronger competition and lower mutation rates decrease Shanon diversity as fewer circulating antigenic types co-exist. <bold>(B)</bold>
Single strain dominance based on the quantity ε from <xref rid="ppat.1003104-Recker1" ref-type="bibr">[22]</xref>
(see <xref ref-type="sec" rid="s4">Methods</xref>
). With stronger competition and lower mutation rates epidemics are contain a larger fraction of a single antigenic type <bold>(C)</bold>
Shannon diversity decreases with stronger positive selection (<italic>ρ</italic>
= −0.88 when MKR>1) and with stronger negative selection (<italic>ρ</italic>
= 0.39 when MKR<1). Positive selection roughly corresponds to lower mutation rates (ξ<10<sup>−3</sup>
), while negative selection corresponds to higher mutation rates (ξ>10<sup>−4</sup>
) <bold>(D)</bold>
Single strain dominance increases with stronger positive selection (<italic>ρ</italic>
= 0.88 when MKR>1) and an increases for stronger negative selection (<italic>ρ</italic>
= −0.23 when MKR<1). (see <xref ref-type="sec" rid="s4">methods</xref>
for full description of epidemiological parameters).</p>
<p>(PDF)</p>
</caption>
<media xlink:href="ppat.1003104.s003.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="ppat.1003104.s004"><label>Figure S3</label>
<caption><p><bold>Seasonal patterns.</bold>
For simulations including metapopulation sinusoidal seasonal forcing was used (see <xref ref-type="sec" rid="s4">Methods</xref>
). Contact rate was modulated sinusoidally with 14% amplitude in temperate demes, and lower biannual seasonal cycles of weaker (7%) amplitude in the tropics. The observed seasonal patterns in the simulation include annual peaks centered around Jan–Feb in the northern hemisphere, July in the southern hemisphere and weaker peaks centered around late July and mid January in the tropics.</p>
<p>(PDF)</p>
</caption>
<media xlink:href="ppat.1003104.s004.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="ppat.1003104.s005"><label>Figure S4</label>
<caption><p><bold>Dynamic changes in the percentage of antigenic clusters and the dominance of antigenic variants in different metapopulation demes.</bold>
<bold>(A)</bold>
Changes in the percentage of the population infected with a specific antigenic cluster variant for the northern hemisphere, the southern hemisphere and the tropics. On average we observe 13±6 antigenic clusters that come to dominate the virus population over the course of the 40 year simulation with an average duration of 4±2 years. One or two clusters usually dominate the deme population. Clusters are defined based on a threshold set when the cumulative change of two or more epitopes between any previous cluster antigenic-type is reached, based on <xref rid="ppat.1003104-Huang1" ref-type="bibr">[36]</xref>
(see <xref ref-type="sec" rid="s4">Methods</xref>
). Clusters are only used for coloring of strains and figures and do not affect the model dynamics. A time window of 25 years was selected for comparison with <xref rid="ppat.1003104-Du1" ref-type="bibr">[42]</xref>
<bold>(B)</bold>
Measurement of the onset time for all antigenic types (prior to cluster subdivision). Onset time was measured as the point where prevalence was estimated to reach 5% of its overall deme prevalence. Antigenic variants are more likely to reach significant prevalence in the tropics: 2±1.5 months earlier in the tropics compared to the northern hemisphere and 3±2 months earlier in the tropics compared to the southern hemisphere (p<0.001 for the combined results).</p>
<p>(PDF)</p>
</caption>
<media xlink:href="ppat.1003104.s005.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="ppat.1003104.s006"><label>Figure S5</label>
<caption><p><bold>Antigenic types across 40 years of simulated years.</bold>
A sample of antigenic types that emerge in the simulation are sampled and numbered sequentially. Red – Antigenic variants sampled from the trunk of the tree (fixed). Black – Antigenic variants sampled from sidebranches of the phylogenetic tree.</p>
<p>(PDF)</p>
</caption>
<media xlink:href="ppat.1003104.s006.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="ppat.1003104.s007"><label>Figure S6</label>
<caption><p><bold>Crossimmunity patterns for individual epitopes.</bold>
The thirteen most prevalent antigenic types from a span of 40 years of simulation were sampled (<xref ref-type="fig" rid="ppat-1003104-g006">Figure 6</xref>
) and ordered by year of introduction. Individual epitopes were compared between the strains. Epitopes with lower variability (2–3 variants per epitope) show a larger degree of reemergence while epitopes with higher variability (4–5 variants per epitope) show a lower degree of epitope reemergence.</p>
<p>(PDF)</p>
</caption>
<media xlink:href="ppat.1003104.s007.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="ppat.1003104.s008"><label>Figure S7</label>
<caption><p><bold>Full phylogenetic trees following the initial introduction of a virus with an increasing mutation rate.</bold>
Phylogenetic trees are based on samples of directly measured virus genealogy in the simulation (434 years). <xref ref-type="fig" rid="ppat-1003104-g004">Figure 4</xref>
in the main body of the paper shows the last 40 years of a simulation with the same parameters (see caption of that figure for details). For these longer sampling windows, extinction was prevented by maintaining at least 50 infected individuals. <bold>(A)</bold>
Model with no mutation. <bold>(B)</bold>
Model with low mutation rate of ξ = 7.5×10<sup>−6</sup>
. <bold>(C)</bold>
Mutation rate of ξ = 7.5×10<sup>−5</sup>
. <bold>(D)</bold>
Mutation rate of ξ = 7.5×10<sup>−4</sup>
. <bold>(E)</bold>
Mutation rate of ξ = 7.5×10<sup>−3</sup>
.</p>
<p>(PDF)</p>
</caption>
<media xlink:href="ppat.1003104.s008.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="ppat.1003104.s009"><label>Table S1</label>
<caption><p><bold>Model results for doubling the number of epitopes or variants per epitope.</bold>
</p>
<p>(PDF)</p>
</caption>
<media xlink:href="ppat.1003104.s009.pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
</sec>
</body>
<back><ack><p>We thank Pamela Martinez, Oliver Pybus and Adam Lauring for comments and discussion.</p>
</ack>
<ref-list><title>References</title>
<ref id="ppat.1003104-1"><label>1</label>
<mixed-citation publication-type="other">(2009) WHO | Influenza (Seasonal).</mixed-citation>
</ref>
<ref id="ppat.1003104-CDC1"><label>2</label>
<mixed-citation publication-type="journal"><collab>C.D.C</collab>
(<year>2010</year>
) <article-title>Estimates of Deaths Associated with Seasonal Influenza - United States, 1976–2007</article-title>
. <source>Morbidity and Mortality Weekly Report (MMWR)</source>
<volume>59</volume>
: <fpage>1057</fpage>
–<lpage>1062</lpage>
.<pub-id pub-id-type="pmid">20798667</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Walsh1"><label>3</label>
<mixed-citation publication-type="journal"><name><surname>Walsh</surname>
<given-names>JA</given-names>
</name>
, <name><surname>Maher</surname>
<given-names>C</given-names>
</name>
(<year>2011</year>
) <article-title>Economic implications of influenza and influenza vaccine</article-title>
. <source>Influenza Vaccines for the Future</source>
<fpage>425</fpage>
–<lpage>440</lpage>
.</mixed-citation>
</ref>
<ref id="ppat.1003104-Fitch1"><label>4</label>
<mixed-citation publication-type="journal"><name><surname>Fitch</surname>
<given-names>WM</given-names>
</name>
, <name><surname>Bush</surname>
<given-names>RM</given-names>
</name>
, <name><surname>Bender</surname>
<given-names>CA</given-names>
</name>
, <name><surname>Cox</surname>
<given-names>NJ</given-names>
</name>
(<year>1997</year>
) <article-title>Long term trends in the evolution of H(3) HA1 human influenza type A</article-title>
. <source>Proc Natl Acad Sci USA</source>
<volume>94</volume>
: <fpage>7712</fpage>
–<lpage>7718</lpage>
.<pub-id pub-id-type="pmid">9223253</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Bush1"><label>5</label>
<mixed-citation publication-type="journal"><name><surname>Bush</surname>
<given-names>RM</given-names>
</name>
, <name><surname>Bender</surname>
<given-names>CA</given-names>
</name>
, <name><surname>Subbarao</surname>
<given-names>K</given-names>
</name>
, <name><surname>Cox</surname>
<given-names>NJ</given-names>
</name>
, <name><surname>Fitch</surname>
<given-names>WM</given-names>
</name>
(<year>1999</year>
) <article-title>Predicting the evolution of human influenza A</article-title>
. <source>Science</source>
<volume>286</volume>
: <fpage>1921</fpage>
–<lpage>1925</lpage>
.<pub-id pub-id-type="pmid">10583948</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Rambaut1"><label>6</label>
<mixed-citation publication-type="journal"><name><surname>Rambaut</surname>
<given-names>A</given-names>
</name>
, <name><surname>Pybus</surname>
<given-names>OG</given-names>
</name>
, <name><surname>Nelson</surname>
<given-names>MI</given-names>
</name>
, <name><surname>Viboud</surname>
<given-names>C</given-names>
</name>
, <name><surname>Taubenberger</surname>
<given-names>JK</given-names>
</name>
, <etal>et al</etal>
(<year>2008</year>
) <article-title>The genomic and epidemiological dynamics of human influenza A virus</article-title>
. <source>Nature</source>
<volume>453</volume>
: <fpage>615</fpage>
–<lpage>619</lpage>
.<pub-id pub-id-type="pmid">18418375</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Bedford1"><label>7</label>
<mixed-citation publication-type="journal"><name><surname>Bedford</surname>
<given-names>T</given-names>
</name>
, <name><surname>Cobey</surname>
<given-names>S</given-names>
</name>
, <name><surname>Pascual</surname>
<given-names>M</given-names>
</name>
(<year>2011</year>
) <article-title>Strength and tempo of selection revealed in viral gene genealogies</article-title>
. <source>BMC Evolutionary Biology</source>
<volume>11</volume>
: <fpage>220</fpage>
–<lpage>220</lpage>
.<pub-id pub-id-type="pmid">21787390</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Gerhard1"><label>8</label>
<mixed-citation publication-type="journal"><name><surname>Gerhard</surname>
<given-names>W</given-names>
</name>
, <name><surname>Webster</surname>
<given-names>RG</given-names>
</name>
(<year>1978</year>
) <article-title>Antigenic drift in influenza A viruses. I. Selection and characterization of antigenic variants of A/PR/8/34 [HON1] influenza virus with monoclonal antibodies</article-title>
. <source>The Journal of experimental medicine</source>
<volume>148</volume>
: <fpage>383</fpage>
–<lpage>383</lpage>
.<pub-id pub-id-type="pmid">359746</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Hay1"><label>9</label>
<mixed-citation publication-type="journal"><name><surname>Hay</surname>
<given-names>AJ</given-names>
</name>
(<year>2001</year>
) <article-title>The evolution of human influenza viruses</article-title>
. <source>Philosophical Transactions of the Royal Society B: Biological Sciences</source>
<volume>356</volume>
: <fpage>1861</fpage>
–<lpage>1870</lpage>
.</mixed-citation>
</ref>
<ref id="ppat.1003104-Nakajima1"><label>10</label>
<mixed-citation publication-type="journal"><name><surname>Nakajima</surname>
<given-names>K</given-names>
</name>
, <name><surname>Nobusawa</surname>
<given-names>E</given-names>
</name>
, <name><surname>Nagy</surname>
<given-names>A</given-names>
</name>
, <name><surname>Nakajima</surname>
<given-names>S</given-names>
</name>
(<year>2005</year>
) <article-title>Accumulation of amino acid substitutions promotes irreversible structural changes in the hemagglutinin of human influenza AH3 virus during evolution</article-title>
. <source>J Virol</source>
<volume>79</volume>
: <fpage>6472</fpage>
–<lpage>6477</lpage>
.<pub-id pub-id-type="pmid">15858030</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Fitch2"><label>11</label>
<mixed-citation publication-type="journal"><name><surname>Fitch</surname>
<given-names>WM</given-names>
</name>
, <name><surname>Leiter</surname>
<given-names>JM</given-names>
</name>
, <name><surname>Li</surname>
<given-names>XQ</given-names>
</name>
, <name><surname>Palese</surname>
<given-names>P</given-names>
</name>
(<year>1991</year>
) <article-title>Positive Darwinian evolution in human influenza A viruses</article-title>
. <source>Proc Natl Acad Sci U S A</source>
<volume>88</volume>
: <fpage>4270</fpage>
–<lpage>4274</lpage>
.<pub-id pub-id-type="pmid">1840695</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Suzuki1"><label>12</label>
<mixed-citation publication-type="journal"><name><surname>Suzuki</surname>
<given-names>Y</given-names>
</name>
(<year>2008</year>
) <article-title>Positive selection operates continuously on hemagglutinin during evolution of H3N2 human influenza A virus</article-title>
. <source>Gene</source>
<volume>427</volume>
: <fpage>111</fpage>
–<lpage>116</lpage>
.<pub-id pub-id-type="pmid">18848975</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Wolf1"><label>13</label>
<mixed-citation publication-type="journal"><name><surname>Wolf</surname>
<given-names>YI</given-names>
</name>
, <name><surname>Viboud</surname>
<given-names>C</given-names>
</name>
, <name><surname>Holmes</surname>
<given-names>EC</given-names>
</name>
, <name><surname>Koonin</surname>
<given-names>EV</given-names>
</name>
, <name><surname>Lipman</surname>
<given-names>DJ</given-names>
</name>
(<year>2006</year>
) <article-title>Long intervals of stasis punctuated by bursts of positive selection in the seasonal evolution of influenza A virus</article-title>
. <source>Biology Direct</source>
<volume>1</volume>
: <fpage>34</fpage>
–<lpage>34</lpage>
.<pub-id pub-id-type="pmid">17067369</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Lindstrom1"><label>14</label>
<mixed-citation publication-type="journal"><name><surname>Lindstrom</surname>
<given-names>SE</given-names>
</name>
, <name><surname>Hiromoto</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Nishimura</surname>
<given-names>H</given-names>
</name>
, <name><surname>Saito</surname>
<given-names>T</given-names>
</name>
, <name><surname>Nerome</surname>
<given-names>R</given-names>
</name>
, <etal>et al</etal>
(<year>1999</year>
) <article-title>Comparative analysis of evolutionary mechanisms of the hemagglutinin and three internal protein genes of influenza B virus: multiple cocirculating lineages and frequent reassortment of the NP, M, and NS genes</article-title>
. <source>Journal of virology</source>
<volume>73</volume>
: <fpage>4413</fpage>
–<lpage>4413</lpage>
.<pub-id pub-id-type="pmid">10196339</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Tria1"><label>15</label>
<mixed-citation publication-type="journal"><name><surname>Tria</surname>
<given-names>F</given-names>
</name>
, <name><surname>Laessig</surname>
<given-names>M</given-names>
</name>
, <name><surname>Peliti</surname>
<given-names>L</given-names>
</name>
, <name><surname>Franz</surname>
<given-names>S</given-names>
</name>
(<year>2005</year>
) <article-title>A minimal stochastic model for influenza evolution</article-title>
. <source>Journal of Statistical Mechanics</source>
Theory and Experiment: P07008/doi:10.1088/1742-5468/2005/07/P07008.</mixed-citation>
</ref>
<ref id="ppat.1003104-Ferguson1"><label>16</label>
<mixed-citation publication-type="journal"><name><surname>Ferguson</surname>
<given-names>NM</given-names>
</name>
, <name><surname>Galvani</surname>
<given-names>AP</given-names>
</name>
, <name><surname>Bush</surname>
<given-names>RM</given-names>
</name>
(<year>2003</year>
) <article-title>Ecological and immunological determinants of influenza evolution</article-title>
. <source>Nature</source>
<volume>422</volume>
: <fpage>428</fpage>
–<lpage>433</lpage>
.<pub-id pub-id-type="pmid">12660783</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Koelle1"><label>17</label>
<mixed-citation publication-type="journal"><name><surname>Koelle</surname>
<given-names>K</given-names>
</name>
, <name><surname>Cobey</surname>
<given-names>S</given-names>
</name>
, <name><surname>Grenfell</surname>
<given-names>B</given-names>
</name>
, <name><surname>Pascual</surname>
<given-names>M</given-names>
</name>
(<year>2006</year>
) <article-title>Epochal Evolution Shapes the Phylodynamics of Interpandemic Influenza A (H3N2) in Humans</article-title>
. <source>Science</source>
<volume>314</volume>
: <fpage>1898</fpage>
–<lpage>1903</lpage>
.<pub-id pub-id-type="pmid">17185596</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Smith1"><label>18</label>
<mixed-citation publication-type="journal"><name><surname>Smith</surname>
<given-names>DJ</given-names>
</name>
(<year>2004</year>
) <article-title>Mapping the Antigenic and Genetic Evolution of Influenza Virus</article-title>
. <source>Science</source>
<volume>305</volume>
: <fpage>371</fpage>
–<lpage>376</lpage>
.<pub-id pub-id-type="pmid">15218094</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Andreasen1"><label>19</label>
<mixed-citation publication-type="journal"><name><surname>Andreasen</surname>
<given-names>V</given-names>
</name>
, <name><surname>Sasaki</surname>
<given-names>A</given-names>
</name>
(<year>2006</year>
) <article-title>Shaping the phylogenetic tree of influenza by cross-immunity</article-title>
. <source>Theoretical population biology</source>
<volume>70</volume>
: <fpage>164</fpage>
–<lpage>173</lpage>
.<pub-id pub-id-type="pmid">16723145</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Kryazhimskiy1"><label>20</label>
<mixed-citation publication-type="journal"><name><surname>Kryazhimskiy</surname>
<given-names>S</given-names>
</name>
, <name><surname>Dieckmann</surname>
<given-names>U</given-names>
</name>
, <name><surname>Levin</surname>
<given-names>SA</given-names>
</name>
, <name><surname>Dushoff</surname>
<given-names>J</given-names>
</name>
(<year>2007</year>
) <article-title>On State-Space Reduction in Multi-Strain Pathogen Models, with an Application to Antigenic Drift in Influenza A</article-title>
. <source>PLoS Computational Biology</source>
<volume>3</volume>
: <fpage>e159</fpage>
–<lpage>e159</lpage>
.<pub-id pub-id-type="pmid">17708677</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Bedford2"><label>21</label>
<mixed-citation publication-type="journal"><name><surname>Bedford</surname>
<given-names>T</given-names>
</name>
, <name><surname>Rambaut</surname>
<given-names>A</given-names>
</name>
, <name><surname>Pascual</surname>
<given-names>M</given-names>
</name>
(<year>2011</year>
) <article-title>Canalization of the evolutionary trajectory of the human influenza virus</article-title>
. <source>Arxiv preprint arXiv</source>
<fpage>11114579</fpage>
.</mixed-citation>
</ref>
<ref id="ppat.1003104-Recker1"><label>22</label>
<mixed-citation publication-type="journal"><name><surname>Recker</surname>
<given-names>M</given-names>
</name>
, <name><surname>Pybus</surname>
<given-names>OG</given-names>
</name>
, <name><surname>Nee</surname>
<given-names>S</given-names>
</name>
, <name><surname>Gupta</surname>
<given-names>S</given-names>
</name>
(<year>2007</year>
) <article-title>The generation of influenza outbreaks by a network of host immune responses against a limited set of antigenic types</article-title>
. <source>Proceedings of the National Academy of Sciences</source>
<volume>104</volume>
: <fpage>7711</fpage>
–<lpage>7711</lpage>
.</mixed-citation>
</ref>
<ref id="ppat.1003104-Russell1"><label>23</label>
<mixed-citation publication-type="journal"><name><surname>Russell</surname>
<given-names>CA</given-names>
</name>
, <name><surname>Jones</surname>
<given-names>TC</given-names>
</name>
, <name><surname>Barr</surname>
<given-names>IG</given-names>
</name>
, <name><surname>Cox</surname>
<given-names>NJ</given-names>
</name>
, <name><surname>Garten</surname>
<given-names>RJ</given-names>
</name>
, <etal>et al</etal>
(<year>2008</year>
) <article-title>The Global Circulation of Seasonal Influenza A (H3N2) Viruses</article-title>
. <source>Science</source>
<volume>320</volume>
: <fpage>340</fpage>
–<lpage>346</lpage>
.<pub-id pub-id-type="pmid">18420927</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Kuiken1"><label>24</label>
<mixed-citation publication-type="journal"><name><surname>Kuiken</surname>
<given-names>T</given-names>
</name>
(<year>2006</year>
) <article-title>Host Species Barriers to Influenza Virus Infections</article-title>
. <source>Science</source>
<volume>312</volume>
: <fpage>394</fpage>
–<lpage>397</lpage>
.<pub-id pub-id-type="pmid">16627737</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Gupta1"><label>25</label>
<mixed-citation publication-type="journal"><name><surname>Gupta</surname>
<given-names>S</given-names>
</name>
(<year>1998</year>
) <article-title>Chaos, Persistence, and Evolution of Strain Structure in Antigenically Diverse Infectious Agents</article-title>
. <source>Science</source>
<volume>280</volume>
: <fpage>912</fpage>
–<lpage>915</lpage>
.<pub-id pub-id-type="pmid">9572737</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Buckee1"><label>26</label>
<mixed-citation publication-type="journal"><name><surname>Buckee</surname>
<given-names>CO</given-names>
</name>
, <name><surname>Jolley</surname>
<given-names>KA</given-names>
</name>
, <name><surname>Recker</surname>
<given-names>M</given-names>
</name>
, <name><surname>Penman</surname>
<given-names>B</given-names>
</name>
, <name><surname>Kriz</surname>
<given-names>P</given-names>
</name>
, <etal>et al</etal>
(<year>2008</year>
) <article-title>Role of selection in the emergence of lineages and the evolution of virulence in Neisseria meningitidis</article-title>
. <source>Proceedings of the National Academy of Sciences</source>
<volume>105</volume>
: <fpage>15082</fpage>
–<lpage>15082</lpage>
.</mixed-citation>
</ref>
<ref id="ppat.1003104-May1"><label>27</label>
<mixed-citation publication-type="journal"><name><surname>May</surname>
<given-names>RM</given-names>
</name>
, <name><surname>Nowak</surname>
<given-names>MA</given-names>
</name>
(<year>1994</year>
) <article-title>Superinfection, Metapopulation Dynamics, and the Evolution of Diversity</article-title>
. <source>Journal of Theoretical Biology</source>
<volume>170</volume>
: <fpage>95</fpage>
–<lpage>114</lpage>
.<pub-id pub-id-type="pmid">7967636</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Nowak1"><label>28</label>
<mixed-citation publication-type="journal"><name><surname>Nowak</surname>
<given-names>MA</given-names>
</name>
, <name><surname>May</surname>
<given-names>RM</given-names>
</name>
(<year>1994</year>
) <article-title>Superinfection and the Evolution of Parasite Virulence</article-title>
. <source>Proceedings of the Royal Society of London Series B: Biological Sciences</source>
<volume>255</volume>
: <fpage>81</fpage>
–<lpage>89</lpage>
.<pub-id pub-id-type="pmid">8153140</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Volz1"><label>29</label>
<mixed-citation publication-type="journal"><name><surname>Volz</surname>
<given-names>EM</given-names>
</name>
(<year>2012</year>
) <article-title>Complex Population Dynamics and the Coalescent under Neutrality</article-title>
. <source>Genetics</source>
<volume>190</volume>
: <fpage>187</fpage>
–<lpage>201</lpage>
.<pub-id pub-id-type="pmid">22042576</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Nowak2"><label>30</label>
<mixed-citation publication-type="journal"><name><surname>Nowak</surname>
<given-names>MA</given-names>
</name>
(<year>1992</year>
) <article-title>What is a quasispecies?</article-title>
<source>Trends in Ecology & Evolution</source>
<volume>7</volume>
: <fpage>118</fpage>
–<lpage>121</lpage>
.<pub-id pub-id-type="pmid">21235976</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-McDonald1"><label>31</label>
<mixed-citation publication-type="journal"><name><surname>McDonald</surname>
<given-names>JH</given-names>
</name>
, <name><surname>Kreitman</surname>
<given-names>M</given-names>
</name>
(<year>1991</year>
) <article-title>Adaptive protein evolution at the Adh locus in Drosophila</article-title>
. <source>Nature</source>
<volume>351</volume>
: <fpage>652</fpage>
–<lpage>654</lpage>
.<pub-id pub-id-type="pmid">1904993</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Bush2"><label>32</label>
<mixed-citation publication-type="journal"><name><surname>Bush</surname>
<given-names>RM</given-names>
</name>
, <name><surname>Fitch</surname>
<given-names>WM</given-names>
</name>
, <name><surname>Bender</surname>
<given-names>CA</given-names>
</name>
, <name><surname>Cox</surname>
<given-names>NJ</given-names>
</name>
(<year>1999</year>
) <article-title>Positive selection on the H3 hemagglutinin gene of human influenza virus A</article-title>
. <source>Mol Biol Evol</source>
<volume>16</volume>
: <fpage>1457</fpage>
–<lpage>1465</lpage>
.<pub-id pub-id-type="pmid">10555276</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Bhatt1"><label>33</label>
<mixed-citation publication-type="journal"><name><surname>Bhatt</surname>
<given-names>S</given-names>
</name>
, <name><surname>Katzourakis</surname>
<given-names>A</given-names>
</name>
, <name><surname>Pybus</surname>
<given-names>OG</given-names>
</name>
(<year>2010</year>
) <article-title>Detecting natural selection in RNA virus populations using sequence summary statistics</article-title>
. <source>Infection, Genetics and Evolution</source>
<volume>10</volume>
: <fpage>421</fpage>
–<lpage>430</lpage>
.</mixed-citation>
</ref>
<ref id="ppat.1003104-Abrams1"><label>34</label>
<mixed-citation publication-type="other">Abrams P (1983) The theory of limiting similarity. Palo Alto, CA, ETATS-UNIS: Annual Reviews.</mixed-citation>
</ref>
<ref id="ppat.1003104-Chew1"><label>35</label>
<mixed-citation publication-type="journal"><name><surname>Chew</surname>
<given-names>FT</given-names>
</name>
, <name><surname>Doraisingham</surname>
<given-names>S</given-names>
</name>
, <name><surname>Ling</surname>
<given-names>AE</given-names>
</name>
, <name><surname>Kumarasinghe</surname>
<given-names>G</given-names>
</name>
, <name><surname>Lee</surname>
<given-names>BW</given-names>
</name>
(<year>1998</year>
) <article-title>Seasonal trends of viral respiratory tract infections in the tropics</article-title>
. <source>Epidemiol Infect</source>
<volume>121</volume>
: <fpage>121</fpage>
–<lpage>128</lpage>
.<pub-id pub-id-type="pmid">9747763</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Huang1"><label>36</label>
<mixed-citation publication-type="journal"><name><surname>Huang</surname>
<given-names>JW</given-names>
</name>
, <name><surname>Yang</surname>
<given-names>JM</given-names>
</name>
(<year>2011</year>
) <article-title>Changed epitopes drive the antigenic drift for influenza A (H3N2) viruses</article-title>
. <source>BMC Bioinformatics</source>
<volume>12 Suppl 1</volume>
: <fpage>S31</fpage>
.<pub-id pub-id-type="pmid">21342562</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Fleming1"><label>37</label>
<mixed-citation publication-type="journal"><name><surname>Fleming</surname>
<given-names>DM</given-names>
</name>
, <name><surname>Ayres</surname>
<given-names>JG</given-names>
</name>
(<year>1988</year>
) <article-title>Diagnosis and patterns of incidence of influenza, influenza-like illness and the common cold in general practice</article-title>
. <source>J R Coll Gen Pract</source>
<volume>38</volume>
: <fpage>159</fpage>
–<lpage>162</lpage>
.<pub-id pub-id-type="pmid">3265157</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Monto1"><label>38</label>
<mixed-citation publication-type="journal"><name><surname>Monto</surname>
<given-names>AS</given-names>
</name>
, <name><surname>Sullivan</surname>
<given-names>KM</given-names>
</name>
(<year>1993</year>
) <article-title>Acute respiratory illness in the community. Frequency of illness and the agents involved</article-title>
. <source>Epidemiology and Infection</source>
<volume>110</volume>
: <fpage>145</fpage>
–<lpage>160</lpage>
.<pub-id pub-id-type="pmid">8432318</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Karpova1"><label>39</label>
<mixed-citation publication-type="journal"><name><surname>Karpova</surname>
<given-names>LS</given-names>
</name>
, <name><surname>Marinich</surname>
<given-names>IG</given-names>
</name>
, <name><surname>Krainova</surname>
<given-names>TI</given-names>
</name>
(<year>2006</year>
) <article-title>[Etiology of influenza according to the data of serological diagnostics in hospital patients in St. Petersburg at the period of 1968–2000]</article-title>
. <source>Zh Mikrobiol Epidemiol Immunobiol</source>
<fpage>41</fpage>
–<lpage>45</lpage>
.<pub-id pub-id-type="pmid">16830588</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-SaltyteBenth1"><label>40</label>
<mixed-citation publication-type="journal"><name><surname>Saltyte Benth</surname>
<given-names>J</given-names>
</name>
, <name><surname>Hofoss</surname>
<given-names>D</given-names>
</name>
(<year>2008</year>
) <article-title>Modelling and prediction of weekly incidence of influenza A specimens in England and Wales</article-title>
. <source>Epidemiol Infect</source>
<volume>136</volume>
: <fpage>1658</fpage>
–<lpage>1666</lpage>
.<pub-id pub-id-type="pmid">18394206</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Sumi1"><label>41</label>
<mixed-citation publication-type="journal"><name><surname>Sumi</surname>
<given-names>A</given-names>
</name>
, <name><surname>Kamo</surname>
<given-names>K</given-names>
</name>
, <name><surname>Ohtomo</surname>
<given-names>N</given-names>
</name>
, <name><surname>Mise</surname>
<given-names>K</given-names>
</name>
, <name><surname>Kobayashi</surname>
<given-names>N</given-names>
</name>
(<year>2011</year>
) <article-title>Time series analysis of incidence data of influenza in Japan</article-title>
. <source>J Epidemiol</source>
<volume>21</volume>
: <fpage>21</fpage>
–<lpage>29</lpage>
.<pub-id pub-id-type="pmid">21088372</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Du1"><label>42</label>
<mixed-citation publication-type="journal"><name><surname>Du</surname>
<given-names>X</given-names>
</name>
, <name><surname>Dong</surname>
<given-names>L</given-names>
</name>
, <name><surname>Lan</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Peng</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Wu</surname>
<given-names>A</given-names>
</name>
, <etal>et al</etal>
(<year>2012</year>
) <article-title>Mapping of H3N2 influenza antigenic evolution in China reveals a strategy for vaccine strain recommendation</article-title>
. <source>Nat Commun</source>
<volume>3</volume>
: <fpage>709</fpage>
.<pub-id pub-id-type="pmid">22426230</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Bedford3"><label>43</label>
<mixed-citation publication-type="journal"><name><surname>Bedford</surname>
<given-names>T</given-names>
</name>
, <name><surname>Cobey</surname>
<given-names>S</given-names>
</name>
, <name><surname>Beerli</surname>
<given-names>P</given-names>
</name>
, <name><surname>Pascual</surname>
<given-names>M</given-names>
</name>
(<year>2010</year>
) <article-title>Global migration dynamics underlie evolution and persistence of human influenza A (H3N2)</article-title>
. <source>PLoS Pathog</source>
<volume>6</volume>
: <fpage>e1000918</fpage>
.<pub-id pub-id-type="pmid">20523898</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Creanza1"><label>44</label>
<mixed-citation publication-type="journal"><name><surname>Creanza</surname>
<given-names>N</given-names>
</name>
, <name><surname>Schwarz</surname>
<given-names>JS</given-names>
</name>
, <name><surname>Cohen</surname>
<given-names>JE</given-names>
</name>
(<year>2010</year>
) <article-title>Intraseasonal dynamics and dominant sequences in H3N2 influenza</article-title>
. <source>PloS one</source>
<volume>5</volume>
: <fpage>e8544</fpage>
–<lpage>e8544</lpage>
.<pub-id pub-id-type="pmid">20049091</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Minayev1"><label>45</label>
<mixed-citation publication-type="journal"><name><surname>Minayev</surname>
<given-names>P</given-names>
</name>
, <name><surname>Ferguson</surname>
<given-names>N</given-names>
</name>
(<year>2009</year>
) <article-title>Incorporating demographic stochasticity into multi-strain epidemic models: application to influenza A</article-title>
. <source>J R Soc Interface</source>
<volume>6</volume>
: <fpage>989</fpage>
–<lpage>996</lpage>
.<pub-id pub-id-type="pmid">19158010</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Yamashita1"><label>46</label>
<mixed-citation publication-type="journal"><name><surname>Yamashita</surname>
<given-names>M</given-names>
</name>
, <name><surname>Krystal</surname>
<given-names>M</given-names>
</name>
, <name><surname>Fitch</surname>
<given-names>WM</given-names>
</name>
, <name><surname>Palese</surname>
<given-names>P</given-names>
</name>
(<year>1988</year>
) <article-title>Influenza B virus evolution: co-circulating lineages and comparison of evolutionary pattern with those of influenza A and C viruses</article-title>
. <source>Virology</source>
<volume>163</volume>
: <fpage>112</fpage>
–<lpage>122</lpage>
.<pub-id pub-id-type="pmid">3267218</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Hensley1"><label>47</label>
<mixed-citation publication-type="journal"><name><surname>Hensley</surname>
<given-names>SE</given-names>
</name>
, <name><surname>Das</surname>
<given-names>SR</given-names>
</name>
, <name><surname>Bailey</surname>
<given-names>AL</given-names>
</name>
, <name><surname>Schmidt</surname>
<given-names>LM</given-names>
</name>
, <name><surname>Hickman</surname>
<given-names>HD</given-names>
</name>
, <etal>et al</etal>
(<year>2009</year>
) <article-title>Hemagglutinin Receptor Binding Avidity Drives Influenza A Virus Antigenic Drift</article-title>
. <source>Science</source>
<volume>326</volume>
: <fpage>734</fpage>
–<lpage>736</lpage>
.<pub-id pub-id-type="pmid">19900932</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Hancock1"><label>48</label>
<mixed-citation publication-type="journal"><name><surname>Hancock</surname>
<given-names>K</given-names>
</name>
, <name><surname>Veguilla</surname>
<given-names>V</given-names>
</name>
, <name><surname>Lu</surname>
<given-names>X</given-names>
</name>
, <name><surname>Zhong</surname>
<given-names>W</given-names>
</name>
, <name><surname>Butler</surname>
<given-names>EN</given-names>
</name>
, <etal>et al</etal>
(<year>2009</year>
) <article-title>Cross-reactive antibody responses to the 2009 pandemic H1N1 influenza virus</article-title>
. <source>N Engl J Med</source>
<volume>361</volume>
: <fpage>1945</fpage>
–<lpage>1952</lpage>
.<pub-id pub-id-type="pmid">19745214</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Itoh1"><label>49</label>
<mixed-citation publication-type="journal"><name><surname>Itoh</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Shinya</surname>
<given-names>K</given-names>
</name>
, <name><surname>Kiso</surname>
<given-names>M</given-names>
</name>
, <name><surname>Watanabe</surname>
<given-names>T</given-names>
</name>
, <name><surname>Sakoda</surname>
<given-names>Y</given-names>
</name>
, <etal>et al</etal>
(<year>2009</year>
) <article-title>In vitro and in vivo characterization of new swine-origin H1N1 influenza viruses</article-title>
. <source>Nature</source>
<volume>460</volume>
: <fpage>1021</fpage>
–<lpage>1025</lpage>
.<pub-id pub-id-type="pmid">19672242</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Bandaranayake1"><label>50</label>
<mixed-citation publication-type="journal"><name><surname>Bandaranayake</surname>
<given-names>D</given-names>
</name>
, <name><surname>Huang</surname>
<given-names>QS</given-names>
</name>
, <name><surname>Bissielo</surname>
<given-names>A</given-names>
</name>
, <name><surname>Wood</surname>
<given-names>T</given-names>
</name>
, <name><surname>Mackereth</surname>
<given-names>G</given-names>
</name>
, <etal>et al</etal>
(<year>2010</year>
) <article-title>Risk factors and immunity in a nationally representative population following the 2009 influenza A(H1N1) pandemic</article-title>
. <source>PloS one</source>
<volume>5</volume>
: <fpage>e13211</fpage>
.<pub-id pub-id-type="pmid">20976224</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Chi1"><label>51</label>
<mixed-citation publication-type="journal"><name><surname>Chi</surname>
<given-names>CY</given-names>
</name>
, <name><surname>Liu</surname>
<given-names>CC</given-names>
</name>
, <name><surname>Lin</surname>
<given-names>CC</given-names>
</name>
, <name><surname>Wang</surname>
<given-names>HC</given-names>
</name>
, <name><surname>Cheng</surname>
<given-names>YT</given-names>
</name>
, <etal>et al</etal>
(<year>2010</year>
) <article-title>Preexisting antibody response against 2009 pandemic influenza H1N1 viruses in the Taiwanese population</article-title>
. <source>Clin Vaccine Immunol</source>
<volume>17</volume>
: <fpage>1958</fpage>
–<lpage>1962</lpage>
.<pub-id pub-id-type="pmid">20876823</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Wikramaratna1"><label>52</label>
<mixed-citation publication-type="journal"><name><surname>Wikramaratna</surname>
<given-names>PS</given-names>
</name>
, <name><surname>Sandeman</surname>
<given-names>M</given-names>
</name>
, <name><surname>Recker</surname>
<given-names>M</given-names>
</name>
, <name><surname>Gupta</surname>
<given-names>S</given-names>
</name>
(<year>2012</year>
) <article-title>The Antigenic Evolution of Influenza: Drift or Thrift?</article-title>
<source>Phil Trans Royal Society B</source>
in press.</mixed-citation>
</ref>
<ref id="ppat.1003104-Tsuchiya1"><label>53</label>
<mixed-citation publication-type="journal"><name><surname>Tsuchiya</surname>
<given-names>E</given-names>
</name>
, <name><surname>Sugawara</surname>
<given-names>K</given-names>
</name>
, <name><surname>Hongo</surname>
<given-names>S</given-names>
</name>
, <name><surname>Matsuzaki</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Muraki</surname>
<given-names>Y</given-names>
</name>
, <etal>et al</etal>
(<year>2001</year>
) <article-title>Antigenic structure of the haemagglutinin of human influenza A/H2N2 virus</article-title>
. <source>J Gen Virol</source>
<volume>82</volume>
: <fpage>2475</fpage>
–<lpage>2484</lpage>
.<pub-id pub-id-type="pmid">11562540</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Reichert1"><label>54</label>
<mixed-citation publication-type="journal"><name><surname>Reichert</surname>
<given-names>T</given-names>
</name>
, <name><surname>Chowell</surname>
<given-names>G</given-names>
</name>
, <name><surname>Nishiura</surname>
<given-names>H</given-names>
</name>
, <name><surname>Christensen</surname>
<given-names>RA</given-names>
</name>
, <name><surname>McCullers</surname>
<given-names>JA</given-names>
</name>
(<year>2010</year>
) <article-title>Does Glycosylation as a modifier of Original Antigenic Sin explain the case age distribution and unusual toxicity in pandemic novel H1N1 influenza?</article-title>
<source>BMC Infect Dis</source>
<volume>10</volume>
: <fpage>5</fpage>
.<pub-id pub-id-type="pmid">20059763</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Bui1"><label>55</label>
<mixed-citation publication-type="journal"><name><surname>Bui</surname>
<given-names>HH</given-names>
</name>
, <name><surname>Peters</surname>
<given-names>B</given-names>
</name>
, <name><surname>Assarsson</surname>
<given-names>E</given-names>
</name>
, <name><surname>Mbawuike</surname>
<given-names>I</given-names>
</name>
, <name><surname>Sette</surname>
<given-names>A</given-names>
</name>
(<year>2007</year>
) <article-title>Ab and T cell epitopes of influenza A virus, knowledge and opportunities</article-title>
. <source>Proc Natl Acad Sci U S A</source>
<volume>104</volume>
: <fpage>246</fpage>
–<lpage>251</lpage>
.<pub-id pub-id-type="pmid">17200302</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Wang1"><label>56</label>
<mixed-citation publication-type="journal"><name><surname>Wang</surname>
<given-names>TT</given-names>
</name>
, <name><surname>Tan</surname>
<given-names>GS</given-names>
</name>
, <name><surname>Hai</surname>
<given-names>R</given-names>
</name>
, <name><surname>Pica</surname>
<given-names>N</given-names>
</name>
, <name><surname>Petersen</surname>
<given-names>E</given-names>
</name>
, <etal>et al</etal>
(<year>2010</year>
) <article-title>Broadly protective monoclonal antibodies against H3 influenza viruses following sequential immunization with different hemagglutinins</article-title>
. <source>PLoS Pathog</source>
<volume>6</volume>
: <fpage>e1000796</fpage>
.<pub-id pub-id-type="pmid">20195520</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Viboud1"><label>57</label>
<mixed-citation publication-type="journal"><name><surname>Viboud</surname>
<given-names>C</given-names>
</name>
, <name><surname>Alonso</surname>
<given-names>WJ</given-names>
</name>
, <name><surname>Simonsen</surname>
<given-names>L</given-names>
</name>
(<year>2006</year>
) <article-title>Influenza in tropical regions</article-title>
. <source>PLoS Med</source>
<volume>3</volume>
: <fpage>e89</fpage>
.<pub-id pub-id-type="pmid">16509764</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Drummond1"><label>58</label>
<mixed-citation publication-type="journal"><name><surname>Drummond</surname>
<given-names>AJ</given-names>
</name>
, <name><surname>Suchard</surname>
<given-names>MA</given-names>
</name>
, <name><surname>Xie</surname>
<given-names>D</given-names>
</name>
, <name><surname>Rambaut</surname>
<given-names>A</given-names>
</name>
(<year>2012</year>
) <article-title>Bayesian phylogenetics with BEAUti and the BEAST 1.7</article-title>
. <source>Molecular Biology and Evolution</source>
<volume>29</volume>
: <fpage>1969</fpage>
–<lpage>1973</lpage>
.<pub-id pub-id-type="pmid">22367748</pub-id>
</mixed-citation>
</ref>
<ref id="ppat.1003104-Lemey1"><label>59</label>
<mixed-citation publication-type="journal"><name><surname>Lemey</surname>
<given-names>P</given-names>
</name>
, <name><surname>Rambaut</surname>
<given-names>A</given-names>
</name>
, <name><surname>Drummond</surname>
<given-names>AJ</given-names>
</name>
, <name><surname>Suchard</surname>
<given-names>MA</given-names>
</name>
(<year>2009</year>
) <article-title>Bayesian phylogeography finds its roots</article-title>
. <source>PLoS computational biology</source>
<volume>5</volume>
: <fpage>e1000520</fpage>
–<lpage>e1000520</lpage>
.<pub-id pub-id-type="pmid">19779555</pub-id>
</mixed-citation>
</ref>
</ref-list>
</back>
</pmc>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Sante/explor/H2N2V1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000C11 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000C11 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Sante |area= H2N2V1 |flux= Pmc |étape= Corpus |type= RBID |clé= |texte= }}
![]() | This area was generated with Dilib version V0.6.33. | ![]() |