Canalization of the evolutionary trajectory of the human influenza virus
Identifieur interne : 000904 ( Pmc/Corpus ); précédent : 000903; suivant : 000905Canalization of the evolutionary trajectory of the human influenza virus
Auteurs : Trevor Bedford ; Andrew Rambaut ; Mercedes PascualSource :
- BMC Biology [ 1741-7007 ] ; 2012.
Abstract
Since its emergence in 1968, influenza A (H3N2) has evolved extensively in genotype and antigenic phenotype. However, despite strong pressure to evolve away from human immunity and to diversify in antigenic phenotype, H3N2 influenza shows paradoxically limited genetic and antigenic diversity present at any one time. Here, we propose a simple model of antigenic evolution in the influenza virus that accounts for this apparent discrepancy.
In this model, antigenic phenotype is represented by a
Selective dynamics induce a canalized evolutionary trajectory, in which the evolutionary fate of the influenza population is surprisingly repeatable. In the model, the influenza population shows a 1- to 2-year timescale of repeatability, suggesting a window in which evolutionary dynamics could be, in theory, predictable.
Url:
DOI: 10.1186/1741-7007-10-38
PubMed: 22546494
PubMed Central: 3373370
Links to Exploration step
PMC:3373370Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Canalization of the evolutionary trajectory of the human influenza virus</title><author><name sortKey="Bedford, Trevor" sort="Bedford, Trevor" uniqKey="Bedford T" first="Trevor" last="Bedford">Trevor Bedford</name><affiliation><nlm:aff id="I1">Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI, 48109, USA</nlm:aff></affiliation><affiliation><nlm:aff id="I2">Howard Hughes Medical Institute, University of Michigan, Ann Arbor, MI, 48109, USA</nlm:aff></affiliation><affiliation><nlm:aff id="I3">Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, EH9 3JT, UK</nlm:aff></affiliation></author><author><name sortKey="Rambaut, Andrew" sort="Rambaut, Andrew" uniqKey="Rambaut A" first="Andrew" last="Rambaut">Andrew Rambaut</name><affiliation><nlm:aff id="I3">Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, EH9 3JT, UK</nlm:aff></affiliation><affiliation><nlm:aff id="I4">Fogarty International Center, National Institutes of Health, Bethesda, MD, 20892, USA</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="I1">Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI, 48109, USA</nlm:aff></affiliation><affiliation><nlm:aff id="I2">Howard Hughes Medical Institute, University of Michigan, Ann Arbor, MI, 48109, USA</nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">22546494</idno><idno type="pmc">3373370</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3373370</idno><idno type="RBID">PMC:3373370</idno><idno type="doi">10.1186/1741-7007-10-38</idno><date when="2012">2012</date><idno type="wicri:Area/Pmc/Corpus">000904</idno><idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000904</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Canalization of the evolutionary trajectory of the human influenza virus</title><author><name sortKey="Bedford, Trevor" sort="Bedford, Trevor" uniqKey="Bedford T" first="Trevor" last="Bedford">Trevor Bedford</name><affiliation><nlm:aff id="I1">Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI, 48109, USA</nlm:aff></affiliation><affiliation><nlm:aff id="I2">Howard Hughes Medical Institute, University of Michigan, Ann Arbor, MI, 48109, USA</nlm:aff></affiliation><affiliation><nlm:aff id="I3">Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, EH9 3JT, UK</nlm:aff></affiliation></author><author><name sortKey="Rambaut, Andrew" sort="Rambaut, Andrew" uniqKey="Rambaut A" first="Andrew" last="Rambaut">Andrew Rambaut</name><affiliation><nlm:aff id="I3">Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, EH9 3JT, UK</nlm:aff></affiliation><affiliation><nlm:aff id="I4">Fogarty International Center, National Institutes of Health, Bethesda, MD, 20892, USA</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="I1">Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI, 48109, USA</nlm:aff></affiliation><affiliation><nlm:aff id="I2">Howard Hughes Medical Institute, University of Michigan, Ann Arbor, MI, 48109, USA</nlm:aff></affiliation></author></analytic><series><title level="j">BMC Biology</title><idno type="eISSN">1741-7007</idno><imprint><date when="2012">2012</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><sec><title>Background</title><p>Since its emergence in 1968, influenza A (H3N2) has evolved extensively in genotype and antigenic phenotype. However, despite strong pressure to evolve away from human immunity and to diversify in antigenic phenotype, H3N2 influenza shows paradoxically limited genetic and antigenic diversity present at any one time. Here, we propose a simple model of antigenic evolution in the influenza virus that accounts for this apparent discrepancy.</p></sec><sec><title>Results</title><p>In this model, antigenic phenotype is represented by a <italic>N</italic>-dimensional vector, and virus mutations perturb phenotype within this continuous Euclidean space. We implement this model in a large-scale individual-based simulation, and in doing so, we find a remarkable correspondence between model behavior and observed influenza dynamics. This model displays rapid evolution but low standing diversity and simultaneously accounts for the epidemiological, genetic, antigenic, and geographical patterns displayed by the virus. We find that evolution away from existing human immunity results in rapid population turnover in the influenza virus and that this population turnover occurs primarily along a single antigenic axis.</p></sec><sec><title>Conclusions</title><p>Selective dynamics induce a canalized evolutionary trajectory, in which the evolutionary fate of the influenza population is surprisingly repeatable. In the model, the influenza population shows a 1- to 2-year timescale of repeatability, suggesting a window in which evolutionary dynamics could be, in theory, predictable.</p></sec></div></front><back><div1 type="bibliography"><listBibl><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Nelson, Mi" uniqKey="Nelson M">MI Nelson</name></author><author><name sortKey="Holmes, Ec" uniqKey="Holmes E">EC Holmes</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="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="Tria, F" uniqKey="Tria F">F Tria</name></author><author><name sortKey="L Ssig, M" uniqKey="L Ssig M">M Lässig</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="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="Recker, M" uniqKey="Recker M">M Recker</name></author><author><name sortKey="Pybus, O" uniqKey="Pybus O">O 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="Hirst, G" uniqKey="Hirst G">G Hirst</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Smith, Dj" uniqKey="Smith D">DJ Smith</name></author><author><name sortKey="Lapedes, As" uniqKey="Lapedes A">AS Lapedes</name></author><author><name sortKey="De Jong, Jc" uniqKey="De Jong J">JC de Jong</name></author><author><name sortKey="Bestebroer, Tm" uniqKey="Bestebroer T">TM Bestebroer</name></author><author><name sortKey="Rimmelzwaan, Gf" uniqKey="Rimmelzwaan G">GF Rimmelzwaan</name></author><author><name sortKey="Osterhaus, Adme" uniqKey="Osterhaus A">ADME Osterhaus</name></author><author><name sortKey="Fouch Ier, Ram" uniqKey="Fouch Ier R">RAM Fouch-ier</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Monto, A" uniqKey="Monto A">A Monto</name></author><author><name sortKey="Sullivan, K" uniqKey="Sullivan K">K Sullivan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Koelle, K" uniqKey="Koelle K">K Koelle</name></author><author><name sortKey="Kamradt, M" uniqKey="Kamradt M">M Kamradt</name></author><author><name sortKey="Pascual, M" uniqKey="Pascual M">M Pascual</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="Bush, R" uniqKey="Bush R">R Bush</name></author><author><name sortKey="Fitch, W" uniqKey="Fitch W">W Fitch</name></author><author><name sortKey="Bender, C" uniqKey="Bender C">C Bender</name></author><author><name sortKey="Cox, N" uniqKey="Cox N">N Cox</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="Park, S" uniqKey="Park S">S Park</name></author><author><name sortKey="Krug, J" uniqKey="Krug J">J Krug</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><author><name sortKey="Gregory, V" uniqKey="Gregory V">V Gregory</name></author><author><name sortKey="Gust, Id" uniqKey="Gust I">ID Gust</name></author><author><name sortKey="Hampson, Aw" uniqKey="Hampson A">AW Hampson</name></author><author><name sortKey="Hay, Aj" uniqKey="Hay A">AJ Hay</name></author><author><name sortKey="Hurt, Ac" uniqKey="Hurt A">AC Hurt</name></author><author><name sortKey="De Jong, Jc" uniqKey="De Jong J">JC de Jong</name></author><author><name sortKey="Kelso, A" uniqKey="Kelso A">A Kelso</name></author><author><name sortKey="Klimov, Ai" uniqKey="Klimov A">AI Klimov</name></author><author><name sortKey="Kageyama, T" uniqKey="Kageyama T">T Kageyama</name></author><author><name sortKey="Komadina, N" uniqKey="Komadina N">N Komadina</name></author><author><name sortKey="Lapedes, As" uniqKey="Lapedes A">AS Lapedes</name></author><author><name sortKey="Lin, Yp" uniqKey="Lin Y">YP Lin</name></author><author><name sortKey="Mosterin, A" uniqKey="Mosterin A">A Mosterin</name></author><author><name sortKey="Obuchi, M" uniqKey="Obuchi M">M Obuchi</name></author><author><name sortKey="Odagiri, T" uniqKey="Odagiri T">T Odagiri</name></author><author><name sortKey="Osterhaus, Adme" uniqKey="Osterhaus A">ADME Osterhaus</name></author><author><name sortKey="Rimmelzwaan, Gf" uniqKey="Rimmelzwaan G">GF Rimmelzwaan</name></author><author><name sortKey="Shaw, Mw" uniqKey="Shaw M">MW Shaw</name></author><author><name sortKey="Skepner, E" uniqKey="Skepner E">E Skepner</name></author><author><name sortKey="Stohr, K" uniqKey="Stohr K">K Stohr</name></author><author><name sortKey="Tashiro, M" uniqKey="Tashiro M">M Tashiro</name></author><author><name sortKey="Fouchier, Ram" uniqKey="Fouchier R">RAM Fouchier</name></author><author><name sortKey="Smith, Dj" uniqKey="Smith D">DJ Smith</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="Adams, B" uniqKey="Adams B">B Adams</name></author><author><name sortKey="Mchardy, A" uniqKey="Mchardy A">A McHardy</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Koelle, K" uniqKey="Koelle K">K Koelle</name></author><author><name sortKey="Ratmann, O" uniqKey="Ratmann O">O Ratmann</name></author><author><name sortKey="Rasmussen, D" uniqKey="Rasmussen D">D Rasmussen</name></author><author><name sortKey="Pasour, V" uniqKey="Pasour V">V Pasour</name></author><author><name sortKey="Mattingly, J" uniqKey="Mattingly J">J Mattingly</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Finkelman, B" uniqKey="Finkelman B">B Finkelman</name></author><author><name sortKey="Viboud, C" uniqKey="Viboud C">C Viboud</name></author><author><name sortKey="Koelle, K" uniqKey="Koelle K">K Koelle</name></author><author><name sortKey="Ferrari, M" uniqKey="Ferrari M">M Ferrari</name></author><author><name sortKey="Bharti, N" uniqKey="Bharti N">N Bharti</name></author><author><name sortKey="Grenfell, B" uniqKey="Grenfell B">B Grenfell</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bhatt, S" uniqKey="Bhatt S">S Bhatt</name></author><author><name sortKey="Holmes, E" uniqKey="Holmes E">E Holmes</name></author><author><name sortKey="Pybus, O" uniqKey="Pybus O">O Pybus</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Goldstein, E" uniqKey="Goldstein E">E Goldstein</name></author><author><name sortKey="Cobey, S" uniqKey="Cobey S">S Cobey</name></author><author><name sortKey="Takahashi, S" uniqKey="Takahashi S">S Takahashi</name></author><author><name sortKey="Miller, J" uniqKey="Miller J">J Miller</name></author><author><name sortKey="Lipsitch, M" uniqKey="Lipsitch M">M Lipsitch</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Voeten, J" uniqKey="Voeten J">J Voeten</name></author><author><name sortKey="Bestebroer, T" uniqKey="Bestebroer T">T Bestebroer</name></author><author><name sortKey="Nieuwkoop, N" uniqKey="Nieuwkoop N">N Nieuwkoop</name></author><author><name sortKey="Fouchier, R" uniqKey="Fouchier R">R Fouchier</name></author><author><name sortKey="Osterhaus, A" uniqKey="Osterhaus A">A Osterhaus</name></author><author><name sortKey="Rimmelzwaan, G" uniqKey="Rimmelzwaan G">G Rimmelzwaan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="O Brien, J" uniqKey="O Brien J">J O'Brien</name></author><author><name sortKey="Minin, V" uniqKey="Minin V">V Minin</name></author><author><name sortKey="Suchard, M" uniqKey="Suchard M">M Suchard</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yang, Z" uniqKey="Yang Z">Z Yang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gould, S" uniqKey="Gould S">S Gould</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Barr, Ig" uniqKey="Barr I">IG Barr</name></author><author><name sortKey="Mccauley, J" uniqKey="Mccauley J">J McCauley</name></author><author><name sortKey="Cox, N" uniqKey="Cox N">N Cox</name></author><author><name sortKey="Daniels, R" uniqKey="Daniels R">R Daniels</name></author><author><name sortKey="Engelhardt, Og" uniqKey="Engelhardt O">OG Engelhardt</name></author><author><name sortKey="Fukuda, K" uniqKey="Fukuda K">K Fukuda</name></author><author><name sortKey="Grohmann, G" uniqKey="Grohmann G">G Grohmann</name></author><author><name sortKey="Hay, A" uniqKey="Hay A">A Hay</name></author><author><name sortKey="Kelso, A" uniqKey="Kelso A">A Kelso</name></author><author><name sortKey="Klimov, A" uniqKey="Klimov A">A Klimov</name></author><author><name sortKey="Odagiri, T" uniqKey="Odagiri T">T Odagiri</name></author><author><name sortKey="Smith, D" uniqKey="Smith D">D Smith</name></author><author><name sortKey="Russell, C" uniqKey="Russell C">C Russell</name></author><author><name sortKey="Tashiro, M" uniqKey="Tashiro M">M Tashiro</name></author><author><name sortKey="Webby, R" uniqKey="Webby R">R Webby</name></author><author><name sortKey="Wood, J" uniqKey="Wood J">J Wood</name></author><author><name sortKey="Ye, Z" uniqKey="Ye Z">Z Ye</name></author><author><name sortKey="Zhang, W" uniqKey="Zhang W">W Zhang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Weinreich, D" uniqKey="Weinreich D">D Weinreich</name></author><author><name sortKey="Delaney, N" uniqKey="Delaney N">N Delaney</name></author><author><name sortKey="Depristo, M" uniqKey="Depristo M">M DePristo</name></author><author><name sortKey="Hartl, D" uniqKey="Hartl D">D Hartl</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lin, J" uniqKey="Lin J">J Lin</name></author><author><name sortKey="Andreasen, V" uniqKey="Andreasen V">V Andreasen</name></author><author><name sortKey="Levin, S" uniqKey="Levin S">S Levin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gog, J" uniqKey="Gog J">J Gog</name></author><author><name sortKey="Grenfell, B" uniqKey="Grenfell B">B Grenfell</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, S" uniqKey="Levin S">S Levin</name></author><author><name sortKey="Dushoff, J" uniqKey="Dushoff J">J Dushoff</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Koelle, K" uniqKey="Koelle K">K Koelle</name></author><author><name sortKey="Khatri, P" uniqKey="Khatri P">P Khatri</name></author><author><name sortKey="Kamradt, M" uniqKey="Kamradt M">M Kamradt</name></author><author><name sortKey="Kepler, T" uniqKey="Kepler T">T Kepler</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Volz, E" uniqKey="Volz E">E Volz</name></author><author><name sortKey="Kosakovsky Pond, S" uniqKey="Kosakovsky Pond S">S Kosakovsky Pond</name></author><author><name sortKey="Ward, M" uniqKey="Ward M">M Ward</name></author><author><name sortKey="Leigh Brown, A" uniqKey="Leigh Brown A">A Leigh Brown</name></author><author><name sortKey="Frost, S" uniqKey="Frost S">S Frost</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="O Dea, E" uniqKey="O Dea E">E O'Dea</name></author><author><name sortKey="Wilke, C" uniqKey="Wilke C">C Wilke</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chowell, G" uniqKey="Chowell G">G Chowell</name></author><author><name sortKey="Miller, M" uniqKey="Miller M">M Miller</name></author><author><name sortKey="Viboud, C" uniqKey="Viboud C">C Viboud</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Viboud, C" uniqKey="Viboud C">C Viboud</name></author><author><name sortKey="Bj Rnstad, O" uniqKey="Bj Rnstad O">O Bjørnstad</name></author><author><name sortKey="Smith, D" uniqKey="Smith D">D Smith</name></author><author><name sortKey="Simonsen, L" uniqKey="Simonsen L">L Simonsen</name></author><author><name sortKey="Miller, M" uniqKey="Miller M">M Miller</name></author><author><name sortKey="Grenfell, B" uniqKey="Grenfell B">B Grenfell</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Carrat, F" uniqKey="Carrat F">F Carrat</name></author><author><name sortKey="Vergu, E" uniqKey="Vergu E">E Vergu</name></author><author><name sortKey="Ferguson, Nm" uniqKey="Ferguson N">NM Ferguson</name></author><author><name sortKey="Lemaitre, M" uniqKey="Lemaitre M">M Lemaitre</name></author><author><name sortKey="Cauchemez, S" uniqKey="Cauchemez S">S Cauchemez</name></author><author><name sortKey="Leach, S" uniqKey="Leach S">S Leach</name></author><author><name sortKey="Valleron, Aj" uniqKey="Valleron A">AJ Valleron</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Park, A" uniqKey="Park A">A Park</name></author><author><name sortKey="Daly, J" uniqKey="Daly J">J Daly</name></author><author><name sortKey="Lewis, N" uniqKey="Lewis N">N Lewis</name></author><author><name sortKey="Smith, D" uniqKey="Smith D">D Smith</name></author><author><name sortKey="Wood, J" uniqKey="Wood J">J Wood</name></author><author><name sortKey="Grenfell, B" uniqKey="Grenfell B">B Grenfell</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gupta, V" uniqKey="Gupta V">V Gupta</name></author><author><name sortKey="Earl, D" uniqKey="Earl D">D Earl</name></author><author><name sortKey="Deem, M" uniqKey="Deem M">M Deem</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Truscott, J" uniqKey="Truscott J">J Truscott</name></author><author><name sortKey="Fraser, C" uniqKey="Fraser C">C Fraser</name></author><author><name sortKey="Cauchemez, S" uniqKey="Cauchemez S">S Cauchemez</name></author><author><name sortKey="Meeyai, A" uniqKey="Meeyai A">A Meeyai</name></author><author><name sortKey="Hinsley, W" uniqKey="Hinsley W">W Hinsley</name></author><author><name sortKey="Donnelly, C" uniqKey="Donnelly C">C Donnelly</name></author><author><name sortKey="Ghani, A" uniqKey="Ghani A">A Ghani</name></author><author><name sortKey="Ferguson, N" uniqKey="Ferguson N">N Ferguson</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><author><name sortKey="Holmes, Ec" uniqKey="Holmes E">EC Holmes</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kucharski, A" uniqKey="Kucharski A">A Kucharski</name></author><author><name sortKey="Gog, Jr" uniqKey="Gog J">JR Gog</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">BMC Biol</journal-id><journal-id journal-id-type="iso-abbrev">BMC Biol</journal-id><journal-title-group><journal-title>BMC Biology</journal-title></journal-title-group><issn pub-type="epub">1741-7007</issn><publisher><publisher-name>BioMed Central</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">22546494</article-id><article-id pub-id-type="pmc">3373370</article-id><article-id pub-id-type="publisher-id">1741-7007-10-38</article-id><article-id pub-id-type="doi">10.1186/1741-7007-10-38</article-id><article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject></subj-group></article-categories><title-group><article-title>Canalization of the evolutionary trajectory of the human influenza virus</article-title></title-group><contrib-group><contrib contrib-type="author" corresp="yes" id="A1"><name><surname>Bedford</surname><given-names>Trevor</given-names></name><xref ref-type="aff" rid="I1">1</xref><xref ref-type="aff" rid="I2">2</xref><xref ref-type="aff" rid="I3">3</xref><email>t.bedford@ed.ac.uk</email></contrib><contrib contrib-type="author" id="A2"><name><surname>Rambaut</surname><given-names>Andrew</given-names></name><xref ref-type="aff" rid="I3">3</xref><xref ref-type="aff" rid="I4">4</xref><email>a.rambaut@ed.ac.uk</email></contrib><contrib contrib-type="author" id="A3"><name><surname>Pascual</surname><given-names>Mercedes</given-names></name><xref ref-type="aff" rid="I1">1</xref><xref ref-type="aff" rid="I2">2</xref><email>pascual@umich.edu</email></contrib></contrib-group><aff id="I1"><label>1</label>Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI, 48109, USA</aff><aff id="I2"><label>2</label>Howard Hughes Medical Institute, University of Michigan, Ann Arbor, MI, 48109, USA</aff><aff id="I3"><label>3</label>Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, EH9 3JT, UK</aff><aff id="I4"><label>4</label>Fogarty International Center, National Institutes of Health, Bethesda, MD, 20892, USA</aff><pub-date pub-type="collection"><year>2012</year></pub-date><pub-date pub-type="epub"><day>30</day><month>4</month><year>2012</year></pub-date><volume>10</volume><fpage>38</fpage><lpage>38</lpage><history><date date-type="received"><day>13</day><month>1</month><year>2012</year></date><date date-type="accepted"><day>30</day><month>4</month><year>2012</year></date></history><permissions><copyright-statement>Copyright ©2012 Bedford et al; licensee BioMed Central Ltd.</copyright-statement><copyright-year>2012</copyright-year><copyright-holder>Bedford et al; licensee BioMed Central Ltd.</copyright-holder><license license-type="open-access" xlink:href="http://creativecommons.org/licenses/by/2.0"><license-p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0">http://creativecommons.org/licenses/by/2.0</ext-link>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</license-p></license></permissions><self-uri xlink:href="http://www.biomedcentral.com/1741-7007/10/38"></self-uri><abstract><sec><title>Background</title><p>Since its emergence in 1968, influenza A (H3N2) has evolved extensively in genotype and antigenic phenotype. However, despite strong pressure to evolve away from human immunity and to diversify in antigenic phenotype, H3N2 influenza shows paradoxically limited genetic and antigenic diversity present at any one time. Here, we propose a simple model of antigenic evolution in the influenza virus that accounts for this apparent discrepancy.</p></sec><sec><title>Results</title><p>In this model, antigenic phenotype is represented by a <italic>N</italic>-dimensional vector, and virus mutations perturb phenotype within this continuous Euclidean space. We implement this model in a large-scale individual-based simulation, and in doing so, we find a remarkable correspondence between model behavior and observed influenza dynamics. This model displays rapid evolution but low standing diversity and simultaneously accounts for the epidemiological, genetic, antigenic, and geographical patterns displayed by the virus. We find that evolution away from existing human immunity results in rapid population turnover in the influenza virus and that this population turnover occurs primarily along a single antigenic axis.</p></sec><sec><title>Conclusions</title><p>Selective dynamics induce a canalized evolutionary trajectory, in which the evolutionary fate of the influenza population is surprisingly repeatable. In the model, the influenza population shows a 1- to 2-year timescale of repeatability, suggesting a window in which evolutionary dynamics could be, in theory, predictable.</p></sec></abstract></article-meta></front><body><sec><title>Background</title><p>Epidemic influenza is responsible for between 250,000 and 500,000 global deaths annually, with influenza A (and in particular, A/H3N2) having caused the bulk of human mortality and morbidity [<xref ref-type="bibr" rid="B1">1</xref>]. Influenza A (H3N2) has continually circulated within the human population since its introduction in 1968, exhibiting recurrent seasonal epidemics in temperate regions and less periodic transmission in the tropics. During this time, H3N2 influenza has continually evolved both genetically and antigenically. Most antigenic drift is thought to be driven by changes to epitopes in the hemagglutinin (HA) protein [<xref ref-type="bibr" rid="B2">2</xref>]. Phylogenetic analysis of the relationships among HA sequences has revealed a distinctive genealogical tree showing a single predominant trunk lineage and side branches that persist for only 1 to 5 years before going extinct [<xref ref-type="bibr" rid="B3">3</xref>]. This tree shape is indicative of serial replacement of strains over time; H3N2 influenza shows rapid evolution, but low standing genetic diversity.</p><p>This observation has remained puzzling from an epidemiological standpoint. Antigenic evolution occurs rapidly, and strong diversifying selection exists to escape from human immunity; why then do we see serial replacement of strains rather than continual genetic and antigenic diversification? Indeed, simple epidemiological models show explosive diversity of genotype and phenotype over time [<xref ref-type="bibr" rid="B4">4</xref>,<xref ref-type="bibr" rid="B5">5</xref>]. Previous work has sought model-based explanations of the limited diversity of influenza, relying on short-lived strain-transcending immunity [<xref ref-type="bibr" rid="B4">4</xref>,<xref ref-type="bibr" rid="B5">5</xref>], complex genotype-to-phenotype maps [<xref ref-type="bibr" rid="B6">6</xref>], or a limited repertoire of antigenic phenotypes [<xref ref-type="bibr" rid="B7">7</xref>].</p><p>Experimental characterization of antigenic phenotype is possible through the hemagglutination inhibition (HI) assay, which measures the cross-reactivity of HA from one virus strain to serum raised against another strain [<xref ref-type="bibr" rid="B8">8</xref>]. The results of many HI assays can be combined to yield a two-dimensional map, representing antigenic similarity and distance between strains as an easily visualized and quantified measure [<xref ref-type="bibr" rid="B9">9</xref>]. The path traced across this map by influenza A (H3N2) from 1968 until present is largely linear, showing serial replacement of one strain by another; there are no major bifurcations of antigenic phenotype [<xref ref-type="bibr" rid="B9">9</xref>].</p><p>Herein, we seek to simultaneously model the genetic and antigenic evolution of the influenza virus. We represent antigenic phenotypes as points in a <italic>N</italic>-dimensional Euclidean space. Based on the finding that a two-dimensional map adequately explains observed antigenic distance between strains [<xref ref-type="bibr" rid="B9">9</xref>], we begin with antigenic phenotypes as points on a plane, but relax this assumption later on in the analysis. After exposure to a virus, a host's risk of infection is proportional to the Euclidean distance between the infecting phenotype and the closest phenotype in the host's immune history. Mutations perturb antigenic phenotype, moving phenotype in a random radial direction and for a randomly distributed distance. We implemented this geometrical model in a large-scale individual-based simulation intended to directly model the antigenic map and genealogical tree of the global influenza population. The simulation includes multiple host populations with different seasonal forcing, hosts with complete immune histories of infection, and viruses with antigenic phenotypes. As the simulation proceeds, infections are tracked and a complete genealogy connecting virus samples is constructed. Results shown here are for simulations of 40 years of virus evolution in a population of 90 million hosts.</p></sec><sec><title>Results and discussion</title><sec><title>Antigenic evolution and genealogical patterns</title><p>The virus persists over the course of the 40-year simulation, and at the end of most simulations, there remain only a few closely related viral lineages, indicating that genealogical diversity is restricted by evolution in the two-dimensional antigenic landscape. Reduced diversity is substantially more common in models with less mutation or models with less variable mutation effects (Figure <xref ref-type="fig" rid="F1">1</xref>). At higher mutation rates, viruses may move apart in antigenic phenotype too rapidly for competition to always eliminate the weaker of two diverging lineages. Similarly, with high variance in mutational effect, there can sometimes emerge new antigenic types, too distant from the existing population to suffer limiting competitive pressure. Both these scenarios lead to coexistence of multiple antigenic phenotypes. We thus restrict the model to parameter regimes with lower mutation rates and lower mutation effect variances. We primarily focus on the model with 10<sup>-4 </sup>mutations per infection per day and mutation effects with a standard deviation of 0.4 antigenic units. In this model, 80 out of the 100 replicate simulations show reduced genealogical diversity (defined as less than 9 years of evolution separating contemporaneous viruses). We conditioned the following analysis on these 80 simulations, compiling summary statistics across this pool and presenting a detailed analysis of a single representative simulation.</p><fig id="F1" position="float"><label>Figure 1</label><caption><p><bold>Genealogical diversity at the end of 40 years across 100 simulations for varying mutational parameters</bold>. Genealogical diversity varies with (<bold>A</bold>) mutation rate and with (<bold>B</bold>) standard deviation of mutation effect. Points represent individual simulation outcomes, and gray bars represent medians and interquartile ranges across replicate simulations. Outcomes with diversity greater than 9 years are shown in blue, and outcomes with diversity less than 9 years are shown in black. Counts of these two classes are shown in blue and black, respectively. Genealogical diversity is measured in years, mutation rate is measured in mutations per infection per day, and standard deviation of mutation effect is measured in antigenic units. Diversity less than 9 years is chosen as a cutoff based on observed patterns of branching in the H3N2 influenza genealogy.</p></caption><graphic xlink:href="1741-7007-10-38-1"></graphic></fig><p>The model exhibits annual winter epidemics in temperate regions and less periodic epidemics in the tropics (Figure <xref ref-type="fig" rid="F2">2A</xref>). Across replicate simulations, we observe average yearly attack rates of 6.8% in temperate regions and rates of 7.1% in the tropics, comparable with estimated attack rates of influenza A (H3N2) of 3% to 8% per year [<xref ref-type="bibr" rid="B10">10</xref>,<xref ref-type="bibr" rid="B11">11</xref>]. Over the course of the simulation, the virus population evolves in antigenic phenotype exhibiting, at any point, a handful of highly abundant phenotypes sampled repeatedly and a large number of phenotypes appearing at low abundance (Figure <xref ref-type="fig" rid="F2">2B</xref>). The observed antigenic map of H3N2 influenza includes substantial experimental noise; replicate strains appear in diverse positions on the observed map. By including measurement noise on antigenic locations (see Methods), we approximate an experimental antigenic map of H3N2 influenza (Figure <xref ref-type="fig" rid="F2">2D</xref>). Over the 40-year simulation, antigenic drift moves the virus population at an average rate across replicate simulations of 1.05 antigenic units per year, corresponding closely to the empirical rate of 1.2 units per year [<xref ref-type="bibr" rid="B9">9</xref>]. The appearance of clusters in the antigenic map comes from the regular spacing of high abundance phenotypes combined with measurement noise. Over time, clusters of antigenically similar strains are replaced by novel clusters of more advanced strains (Figure <xref ref-type="fig" rid="F3">3A</xref>). Across replicate simulations, clusters persist for an average of 5.0 years, measured as the time it takes for a new cluster to reach 10% frequency, peak, and decline to 10% frequency. The transition between clusters occurs quickly, taking an average of 1.8 years.</p><fig id="F2" position="float"><label>Figure 2</label><caption><p><bold>Simulation results showing epidemiological, antigenic, and genealogical dynamics</bold>. (<bold>A</bold>) Weekly time series of incidence of viral infection in north and tropics regions. (<bold>B</bold>) Two-dimensional antigenic phenotypes of 5,943 viruses sampled over the course of the simulation. Each discrete virus phenotype is shown as a bubble, with the bubble area proportional to the number of times this phenotype was sampled. (<bold>C</bold>) Genealogical tree depicting the infection history of 376 samples from the virus population. Parent/offspring relationships were tracked over the course of the simulation, giving a direct observation of the genealogy rather than a phylogenetic inference. (<bold>D</bold>) Antigenic map depicting phenotypes of 5,943 viruses sampled over the course of the simulation. To approximate experimental noise present in the observed antigenic map of H3N2 influenza, noise was added to each sample, and the resulting observations were grouped into 11 clusters and colored accordingly. Grid lines show single units of antigenic distance.</p></caption><graphic xlink:href="1741-7007-10-38-2"></graphic></fig><fig id="F3" position="float"><label>Figure 3</label><caption><p><bold>Antigenic evolution over the course of the 40-year simulation</bold>. (<bold>A</bold>) Proportion of virus population composed of each antigenic cluster through time. (<bold>B</bold>) Antigenic distance from the initial phenotype (<italic>x </italic>= 0, <italic>y </italic>= 0) for each of 5,943 virus samples relative to time of virus sampling. Viruses were sampled at a constant rate proportional to prevalence, and coloring was determined from the antigenic map in Figure 2D.</p></caption><graphic xlink:href="1741-7007-10-38-3"></graphic></fig><p>Remarkably, although antigenic phenotype is free to mutate in any direction in the two-dimensional space, selection pressures force the virus population to move in nearly a straight line in antigenic space (Figure <xref ref-type="fig" rid="F2">2B</xref>). Across replicate simulations, 94% of the variance of antigenic phenotype can be explained by a single dimension of variation. This mirrors the empirical results showing a largely linear antigenic map for H3N2 influenza isolates from 1968 to 2003 [<xref ref-type="bibr" rid="B9">9</xref>]. Because of the primarily one-dimensional movement, antigenic distance from the original phenotype increases nearly linearly with time (Figure <xref ref-type="fig" rid="F3">3B</xref>). Antigenic evolution occurs in a punctuated fashion; periods of relative stasis are interspersed with more rapid antigenic change (Figure <xref ref-type="fig" rid="F3">3B</xref>). Antigenic and epidemiological dynamics show a fundamental linkage so that large jumps of antigenic phenotype result in increased rates of infection (Figure <xref ref-type="fig" rid="F2">2</xref>).</p><p>The genealogical tree connecting the evolving virus population appears characteristically sparse with pronounced trunk and short-lived side branches (Figure <xref ref-type="fig" rid="F2">2C</xref>). This tree shape is reflected in low levels of standing diversity; across replicates, an average of 5.68 years of evolution separates two randomly sampled viruses in the population. This result matches well with the average diversity observed in influenza A (H3N2) of 5.65 years, separating randomly sampled contemporaneous viruses [<xref ref-type="bibr" rid="B12">12</xref>]. A spindly genealogical tree is indicative of population turnover, wherein novel antigenic phenotypes continually replace more primitive 'spent' phenotypes, purging their genealogical diversity. In general, natural selection reduces effective population size and genealogical diversity [<xref ref-type="bibr" rid="B12">12</xref>].</p><p>Selective pressures can be examined by comparing which mutations fix, that is, incorporated into the progenitor trunk lineage, and which mutations are lost, that is, incorporated into side branches bound for extinction. This approach has shown that, in influenza A (H3N2), natural selection promotes mutations to epitope sites in the HA1 region [<xref ref-type="bibr" rid="B13">13</xref>,<xref ref-type="bibr" rid="B14">14</xref>]. By examining antigenic mutations, we find a corresponding effect in simulated evolutionary trajectories (Table <xref ref-type="table" rid="T1">1</xref>). Additionally, we find that trunk mutations occur at strikingly regular intervals, with less variation of waiting times than expected under a simple random process (Figure <xref ref-type="fig" rid="F4">4</xref>). There is a relative scarcity of mutation events occurring in intervals under 1 year and a relative excess of mutation events occurring in 2- to 3-year intervals (Figure <xref ref-type="fig" rid="F4">4</xref>). This pattern arises from clonal interference between competing mutations which reduces variability in the fixation process of adaptive substitutions [<xref ref-type="bibr" rid="B15">15</xref>].</p><table-wrap id="T1" position="float"><label>Table 1</label><caption><p>Rates of mutation and phenotypic change on trunk and side branches and mutational expectation</p></caption><table frame="hsides" rules="groups"><thead><tr><th></th><th align="center">Baseline</th><th align="center">Side branch</th><th align="center">Trunk</th><th align="center">Trunk/side branch</th></tr></thead><tbody><tr><td align="left">Mutation size (antigenic units)</td><td align="center">0.60</td><td align="center">0.79</td><td align="center">1.58</td><td align="center">1.99×</td></tr><tr><td align="left">Mutation rate (mutations per year)</td><td align="center">0.04</td><td align="center">0.06</td><td align="center">0.81</td><td align="center">13.23×</td></tr><tr><td align="left">Antigenic flux (antigenic units per year)</td><td align="center">0.02</td><td align="center">0.05</td><td align="center">1.27</td><td align="center">26.25×</td></tr></tbody></table></table-wrap><fig id="F4" position="float"><label>Figure 4</label><caption><p><bold>Observed vs. expected distributions of waiting times between phenotypic mutations along genealogy trunk</bold>. (<bold>A</bold>) Histogram bins show the observed distribution of waiting times in years across 80 replicate simulations representing 1,584 mutations. The mean of this distribution is 1.76 years. The dashed line shows the Poisson process expectation of exponentially distributed waiting times. (<bold>B</bold>) The density distribution of waiting times is transformed into a hazard function, representing the rate of trunk mutation after a specific waiting time. The dashed line shows the memoryless hazard function of the Poisson process expectation.</p></caption><graphic xlink:href="1741-7007-10-38-4"></graphic></fig></sec><sec><title>Spatial dynamics</title><p>The genealogical tree also contains detailed information on the history of migration between regions. We find that, consistent with empirical estimates [<xref ref-type="bibr" rid="B16">16</xref>,<xref ref-type="bibr" rid="B17">17</xref>], the trunk resides primarily within the tropics, where seasonal dynamics are less prevalent (Figure <xref ref-type="fig" rid="F5">5A</xref>). Across replicate simulations, we observe 72% of the trunk's history within the tropics and 28% within temperate regions. With symmetric host contact rates and equal host population sizes, and without seasonal forcing, we would expect trunk proportions of one third for each region. We calculated rates of migration based on observed event counts across replicate simulations, separating region-specific rates on side branches from region-specific rates on trunk branches. We find that migration patterns on side branches are close to symmetric, with similar rates between all regions, while migration patterns on trunk branches are highly asymmetric, with high rates of movement between temperate regions and from temperate regions into the tropics (Figure <xref ref-type="fig" rid="F5">5B</xref>). Extrapolating from these rates, we arrive at an expected stationary distribution of trunk location of 76% tropics and 24% temperate regions, in line with the observed residency patterns of the trunk. It may at first seem counterintuitive to see higher rates of movement from the temperate regions into the tropics along trunk branches, but it makes sense when thought of in terms of conditional probability. Only those lineages that remain in the tropics, migrate into the tropics, or which rapidly migrate between the north and south have a chance at becoming the trunk lineage, while lineages that remain within the temperate regions are doomed to extinction. Along similar lines, Adams and McHardy [<xref ref-type="bibr" rid="B18">18</xref>] use a modeling approach to show the importance of nonseasonal transmission to the evolution of the virus.</p><fig id="F5" position="float"><label>Figure 5</label><caption><p><bold>Patterns of geographic movement of virus lineages</bold>. (<bold>A</bold>) Evolutionary relationships among 376 viruses sampled evenly through time colored by geographic location. Lineages residing in the north (N), south (S), and tropics (T) are colored yellow, red, and blue, respectively. (<bold>B</bold>) Observed migration rates between regions on side branch lineages (left) and on trunk lineages (right). Arrows denote movement of lineages, and arrow width is proportional to migration rate. Circle area is proportional to the expected stationary frequency of a region given the observed migration rates. In both cases, migration rates are calculated across 80 replicate simulations.</p></caption><graphic xlink:href="1741-7007-10-38-5"></graphic></fig><p>These findings suggest that persistence and migration are fundamentally connected and have important implications for future phylogeographic analyses. Russell <italic>et al. </italic>[<xref ref-type="bibr" rid="B16">16</xref>] emphasize a source-sink model of movement of the HA protein of influenza A (H3N2) based on their finding of a trunk lineage residing within China and the Southeast Asian tropics, whereas Bedford <italic>et al. </italic>[<xref ref-type="bibr" rid="B17">17</xref>] emphasize a global metapopulation model based on phylogenetic inference of migration rates across the entire tree. Our results suggest that both scenarios are simultaneously possible; side branches may be highly volatile, moving rapidly and symmetrically between regions, while the trunk lineage may be more stable remaining within a region (or within a highly connected network of regions) that has more persistent transmission. In light of these results, we suggest that future work on the phylogeography of influenza take into account trunk vs. side branch differences in migration patterns.</p></sec><sec><title>Correspondence between model and data</title><p>In our model, antigenic evolution is driven by the appearance of novel antigenic variants that best escape existing human immunity. Although multiple epidemiological/evolutionary mechanisms have been proposed to explain the restricted genetic diversity and rapid population turnover of influenza A (H3N2) [4-7], our results show that a simple model coupling antigenic and genealogical evolution exhibits broad explanatory power. We find a strong correspondence between the antigenic and genealogical patterns generated by our model (Figure <xref ref-type="fig" rid="F2">2</xref>) and patterns of genetic and antigenic evolution exhibited by influenza A (H3N2) [<xref ref-type="bibr" rid="B3">3</xref>,<xref ref-type="bibr" rid="B9">9</xref>]. Our model simultaneously captures seasonal attack rates, the rate and pattern of antigenic drift, genealogical diversity, and geographic migration patterns.</p><p>Our model predicts that detailed classification of influenza strains will support a relatively small number of predominant phenotypes. Rather than each influenza strain possessing a unique antigenic location, many strains group together with shared antigenic phenotypes (Figure <xref ref-type="fig" rid="F2">2B</xref>). We suggest that a large proportion of intra-cluster variation in the observed antigenic map is due to experimental noise rather than each strain possessing a unique antigenic location. The relationship between Figure <xref ref-type="fig" rid="F2">2B,D</xref> illustrates this effect, where a large number of antigenic locations emerge from a comparatively small number of unique antigenic phenotypes. Additionally, our model accurately predicts the contrasting dynamics of other types/subtypes of influenza. We find that lowering mutation size/effect or lowering intrinsic <italic>R</italic><sub>0 </sub>results in decreased incidence, slower antigenic movement, and greater genealogical diversity, all distinguishing characteristics of H1N1 influenza and influenza B [see Figure S1 in Additional file <xref ref-type="supplementary-material" rid="S1">1</xref>].</p><p>The historical record of influenza evolution suggests that bifurcation of viral lineages is rare, but possible. We have observed no bifurcations in H1N1 influenza from 1918 to 1957 and again from 1977 to 2010, no bifurcations in H2N2 influenza from 1957 to 1968, and no bifurcations of H3N2 influenza from 1968 to 2010. We have observed one bifurcation in influenza B from 1940 to 2010. Thus, ignoring differences between influenza types and subtypes, we have very roughly observed a rate of one major bifurcation in 195 years of evolution. In our model, in 20 out of 100 replicate simulations, we observe a deep bifurcation in the viral genealogy, which translates to observing one deep bifurcation in 200 years of evolution. Thus, we suggest that the 20 of 100 simulations where deep branching occurs are not necessarily evidence of poor model t. Similar to Koelle <italic>et al. </italic>[<xref ref-type="bibr" rid="B19">19</xref>], we assume that although the historical evolution of H3N2 influenza followed the path of a single lineage, it could have included a major bifurcation.</p><p>In our model, when antigenic phenotype remains static, there may be multiple consecutive seasons without appreciable incidence (Figure <xref ref-type="fig" rid="F2">2A</xref>), a pattern apparently absent from H3N2 influenza [<xref ref-type="bibr" rid="B20">20</xref>]. Additionally, we observe antigenic trajectories that are more linear and deterministic than the highly clustered trajectory observed by Smith <italic>et al. </italic>[<xref ref-type="bibr" rid="B9">9</xref>]. We suggest that any model exhibiting punctuated evolution broadly consistent with the punctuated change seen in the experimental antigenic map will show similar patterns of incidence. We can 'fix' the incidence patterns but at the cost of too smooth an antigenic map [see Figure S2 in Additional file <xref ref-type="supplementary-material" rid="S1">1</xref>]. Evolutionary patterns of the neuraminidase (NA) protein may provide an explanation. Epitopes in the HA and NA proteins are jointly responsible for determining antigenicity [<xref ref-type="bibr" rid="B2">2</xref>], and it is now clear that levels of adaptive evolution are similar between HA and NA [<xref ref-type="bibr" rid="B21">21</xref>]. Thus, changes in NA may be driving incidence patterns as well, resulting in an observed time series of incidence partially divorced from the antigenic map of HA. Incorporating antigenic evolution of NA could thus yield a rougher antigenic map for HA, more closely matching experimental results, while simultaneously yielding smoother year-to-year incidence patterns.</p><p>It remains a central question as to the extent that short-lived strain-transcending immunity is responsible for influenza's limited diversity and spindly genealogical tree [<xref ref-type="bibr" rid="B4">4</xref>,<xref ref-type="bibr" rid="B5">5</xref>]. Our findings suggest a possible resolution. Although lacking short-lived immunity, our model shows a detailed correspondence to both the antigenic map and genealogical tree of H3N2 influenza. If an antigenic map were to show a deep bifurcation, where two viral lineages move in different antigenic directions, then we would expect the same bifurcation to be evident in the genealogical tree. Short-lived strain-transcending immunity provides a mechanism by which lineages may diverge in antigenic phenotype but still show epidemiological interference. This mechanism would explain a situation where bifurcations emerge in the antigenic map, but competition results in the extinction of divergent antigenic lineages. In our model, one cluster leads to another cluster in orderly succession, and there is never competition between antigenically distant clusters. Thus, short-lived strain-transcending immunity is not required to limit diversity in the model. This is not to say that short-lived strain-transcending immunity is not present; observed interference between subtypes [<xref ref-type="bibr" rid="B4">4</xref>,<xref ref-type="bibr" rid="B22">22</xref>], evolution at CTL epitopes [<xref ref-type="bibr" rid="B23">23</xref>], and the exclusion of the Beijing/89 cluster by the antigenically distant Beijing/92 cluster [<xref ref-type="bibr" rid="B9">9</xref>] all suggest some form of more general interaction between influenza viruses.</p></sec><sec><title>Linear antigenic movement</title><p>It might seem reasonable for one viral lineage to move in one antigenic direction, while another lineage moves tangentially, eventually resulting in two non-interacting viral lineages. Instead, we find that movement in a single antigenic direction is favored, resulting in most replicate simulations showing low standing diversity (Figure <xref ref-type="fig" rid="F1">1</xref>). The origins of this pattern can be seen in the interaction between virus evolution and host immunity (Figure <xref ref-type="fig" rid="F6">6</xref>). As the virus population evolves forward, it leaves a wake of immunity in the host population, and evolution away from this immunity results in the canalization of the antigenic phenotype; mutations that continue along the line of primary antigenic variation will show a transmission advantage compared to more tangential mutations.</p><fig id="F6" position="float"><label>Figure 6</label><caption><p><bold>Host immunity and antigenic history of the virus population</bold>. Contour lines represent the state of host immunity at the end of the 40-year simulation. They show the mean risk of infection (as a percentage) after a random host in the population encounters a virus bearing a particular antigenic phenotype. Contour lines are spaced in intervals of 2.5%. Bubbles represent a sample of antigenic phenotypes present at the end of the 40-year simulation. The area of each bubble is proportional to the number of samples with this phenotype. Lines leading into these bubbles show past antigenic history. The current phenotypes rapidly coalesce to a trunk phenotype. The movement of the virus population from the left to the center of the figure can be seen from the antigenic history. At the end of the simulation, several virus phenotypes exist with similar antigenic locations; all of these phenotypes lie significantly ahead of the peak of host immunity.</p></caption><graphic xlink:href="1741-7007-10-38-6"></graphic></fig><p>Following the work of Smith <italic>et al. </italic>[<xref ref-type="bibr" rid="B9">9</xref>], it remained an open question of why a two-dimensional map should explain the antigenic variation of H3N2 influenza, although the authors astutely speculated that 'there is a selective advantage for clusters that move away linearly from previous clusters as they most effectively escape existing population-level immunity, and this is a plausible explanation for the somewhat linear antigenic evolution in regions of the antigenic map.' This hypothesis remained to be tested. Here, we show from a simple model of epidemiology and evolution that a linear trajectory of antigenic evolution dynamically emerges due to basic selective pressures. This result simultaneously explains the linear pattern of antigenic drift [<xref ref-type="bibr" rid="B9">9</xref>] and the characteristically spindly genealogical tree [<xref ref-type="bibr" rid="B3">3</xref>] exhibited by influenza A (H3N2).</p><p>For this process to take hold, the virus population needs to be somewhat mutationally limited; if functional antigenic variants of novel phenotype emerge too quickly, then antigenic change will occur too rapidly for competition to winnow down the virus population to a single lineage (Figure <xref ref-type="fig" rid="F1">1</xref>). Assuming that antigenic mutations have an average effect of 0.6 antigenic units and a standard deviation of 0.4 units, then the rate of new antigenic mutations cannot be greater than approximately 10<sup>-4 </sup>mutations per day (Figure <xref ref-type="fig" rid="F1">1</xref>). Thus, it is important that the rate of 10<sup>-4 </sup>mutations per day be biological plausible. Here, we take the rate of synonymous substitution as a proxy for the neutral rate of mutation. The rate of synonymous change has been estimated at 2.5 × 10<sup>-6 </sup>per site per day [<xref ref-type="bibr" rid="B24">24</xref>]. As there are approximately two nonsynonymous sites per codon in influenza [<xref ref-type="bibr" rid="B25">25</xref>], this gives a neutral rate of amino acid change of approximately 5 × 10<sup>-6 </sup>per site per day. Other work has shown that there appear to be approximately 18 amino acid sites implicated in the majority of adaptive change [<xref ref-type="bibr" rid="B13">13</xref>]. These sites evolve along the trunk of the phylogeny at rate of 0.053 substitutions per site per year or at a combined rate of 0.95 substitutions per year [<xref ref-type="bibr" rid="B4">4</xref>]. This result agrees well with our finding of 0.81 antigenic mutations per year on the phylogeny trunk (Table <xref ref-type="table" rid="T1">1</xref>). If we assume 18 sites involved in antigenic change, this gives an overall rate of antigenic mutation of 9 × 10<sup>-5 </sup>per day. Thus, we believe that 10<sup>-4 </sup>mutations per day represents a biologically reasonable estimate.</p><p>To consider to what extent these results were contingent on the dimensionality of the underlying antigenic model, we further implemented our model in a 10-dimensional antigenic space. Here, mutations occur as 10-spheres, but the distance moved by a mutation is the same as in the previous two-dimensional formulation. We arrive at nearly the same results with this model; principal components analysis shows that the first and second dimensions of variation account for 87% and 7%, respectively, of the total variance [see Figure S3 in Additional file <xref ref-type="supplementary-material" rid="S1">1</xref>]. Thus, our model predicts that future work probing mutational effects will support an underlying high-dimensional antigenic space, even though a two-dimensional map is sufficient to explain observed antigenic relationships among evolving strains.</p></sec><sec><title>Winding back the tape</title><p>The 40-year simulation of influenza dynamics shows broad correspondence with observed patterns. However, year-to-year details are not captured, for example, years that undergo antigenic transitions in the 40-year simulation do not match up with observed years of antigenic transitions. Over long time spans, year-to-year correspondence seems impossible to achieve in this sort of stochastic system, where evolution is often driven by chance mutations of large antigenic effect. However, correspondence in the shorter term may be possible. To test this, we examined repeatability in replicate simulations, showing what happens when we 'wind back the tape' [<xref ref-type="bibr" rid="B26">26</xref>] on the evolution of the virus. We ran 100 replicate simulations, each starting from the endpoint of the representative 40-year simulation shown in Figure <xref ref-type="fig" rid="F2">2</xref>. The starting point for these replicate simulations was the exact end state of the 40-year simulation, including the frequencies of every virus strain and the entire host immune profile. These replicate simulations were run for an additional 6 years, and all evolutionary and epidemiological parameters were identical to the initial 40-year simulation.</p><p>Initially, we find a great detail of repeatability; during the first year of evolution, every replicate virus population undergoes a similar antigenic transition (Figure <xref ref-type="fig" rid="F7">7</xref>), resulting in a repeatable peak in northern hemisphere incidence (Figure <xref ref-type="fig" rid="F8">8</xref>). After 3 years, repeatability has mostly disappeared, with antigenic phenotype and incidence appearing highly variable across replicates (Figures <xref ref-type="fig" rid="F7">7</xref> and <xref ref-type="fig" rid="F8">8</xref>). The 1- to 2-year timescale of repeatability can be explained by the presence of standing antigenic variation. In the initial virus population, there are several novel antigenic variants present at low frequency (Figure <xref ref-type="fig" rid="F6">6</xref>), one of which, without fail, comes to predominate the virus population.</p><fig id="F7" position="float"><label>Figure 7</label><caption><p><bold>Antigenic phenotypes over the course of 4 years of evolution across 100 replicate simulations starting from identical initial conditions</bold>. Replicate simulations were initialized with the end state of the 40-year simulation shown in Figure 2. Each panel shows an additional year of evolution, with black points representing the mean antigenic phenotypes of the 100 replicate simulations and gray lines representing the history of each mean antigenic phenotype.</p></caption><graphic xlink:href="1741-7007-10-38-7"></graphic></fig><fig id="F8" position="float"><label>Figure 8</label><caption><p><bold>Time series of incidence across 100 replicate simulations with identical initial conditions</bold>. Panels show incidence in the north, tropics, and south regions over the course of 6 years. Solid black lines represent the median weekly incidence across the 100 replicate simulations, while gray intervals represent the interquartile range across simulations. There is little variability for the first year of replicate simulations. Replicate simulations were initialized with the end state of the 40-year simulation shown in Figure 2.</p></caption><graphic xlink:href="1741-7007-10-38-8"></graphic></fig><p>We see that the initial evolutionary trajectory, during which time standing variation plays out, is highly repeatable and thus predictable given enough information and the right methods of analysis. However, prediction of longer-term evolutionary scenarios will necessarily be difficult or impossible except in a vague sense. Through careful surveillance e orts and genetic and antigenic characterization of influenza strains, the World Health Organization makes twice-yearly vaccine strain recommendations [<xref ref-type="bibr" rid="B27">27</xref>]. It should be possible to combine these sorts of modeling approaches with surveillance data to gauge the likelihood that a sampled variant will spread through the population.</p></sec></sec><sec sec-type="conclusions"><title>Conclusions</title><p>Recent work on empirical fitness landscapes has shown that natural selection follows few mutational paths [<xref ref-type="bibr" rid="B28">28</xref>]. The spindly genealogical tree and the almost linear serial replacement of influenza strains have remained a puzzling phenomenon. We suggest that the evolutionary and epidemiological dynamics displayed by the influenza virus may simply be explained as an outgrowth of selection to avoid host immunity. Natural selection can only 'see' one step ahead and so sacrifices long-term gains for short-term advantages. The result is a canalized evolutionary trajectory lacking antigenic diversification.</p></sec><sec sec-type="methods"><title>Methods</title><sec><title>Transmission model</title><p>To characterize the joint epidemiological, genealogical, antigenic, and spatial patterns of influenza, we implemented a large-scale individual-based model (full simulation source code available at <ext-link ext-link-type="uri" xlink:href="http://www.trevorbedford.com/antigen/">http://www.trevorbedford.com/antigen/</ext-link>). This model consists of daily time steps, in which the states of hosts and viruses are updated. Hosts may be born, may die, may contact other hosts allowing viral transmission, or may recover from infection. Viruses may mutate in antigenic phenotype. Each simulation ran for 40 years of model time.</p><p>Hosts in this model are divided between three regions: north, south, and tropics. There are 30 million hosts within each of the three regions, giving N = 9 × 10<sup>7 </sup>hosts. Host population size remains fixed at this number, but vital dynamics cause births and deaths of hosts at a rate of 1/30 years = 9.1 × 10<sup>-5 </sup>per host per day. Within each region, transmission proceeds through mass action with contacts between hosts occurring at a rate of <italic>β </italic>= 0.36 per host per day. Regional transmission rates in temperate regions vary according to sinusoidal seasonal forcing with amplitude <italic>∊ </italic>= 0.15 and opposite phase in the north and in the south. Transmission rate does not vary over time in the tropics. Transmission between region <italic>i </italic>and region <italic>j </italic>occurs at rate <italic>m</italic><italic>β<sub>i</sub></italic>, where <italic>m </italic>= 0.001 and is the same between each pair of regions and <italic>β<sub>i </sub></italic>is the within-region contact rate. Hosts recover from infection at rate <italic>v </italic>= 0.2 per host per day so that <italic>R</italic><sub>0 </sub>in a naive host population is 1.8. An individual cannot be simultaneously infected with multiple virus lineages; there is no superinfection in the model.</p><p>Each virus possesses an antigenic phenotype, represented as a location in Euclidean space. Here, we primarily use a two-dimensional antigenic location. After recovery, a host 'remembers' the phenotype of its infecting virus as part of its immune history. When a contact event occurs and a virus attempts to infect a host, the Euclidean distance from the infecting phenotype <italic>ϕ<sub>v </sub></italic>is calculated to each of the phenotypes in the host immune history <inline-formula><mml:math id="M1" name="1741-7007-10-38-i1" overflow="scroll"><mml:msub><mml:mrow><mml:mi>ϕ</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>h</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:msub><mml:mo class="MathClass-punc">,</mml:mo><mml:mo class="MathClass-op">…</mml:mo><mml:mo class="MathClass-punc">,</mml:mo><mml:msub><mml:mrow><mml:mi>ϕ</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>h</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub></mml:math></inline-formula>. Here, 1 unit of antigenic distance is designed to correspond to a twofold dilution of antiserum in a HI assay [<xref ref-type="bibr" rid="B9">9</xref>]. The probability that infection occurs after exposure is proportional to the distance <italic>d </italic>to the closest phenotype in the host immune history. Risk of infection follows the form <italic>ρ </italic>= min{<italic>d s</italic>, 1}, where <italic>s </italic>= 0.07. Cross-immunity <italic>σ </italic>equals 1 - <italic>ρ</italic>. With no initial immunity in the host population, the virus undergoes a severe trough in prevalence after the initial pandemic increase. With this number of host individuals, the virus population usually stochastically dies out during this trough. To prevent this, we gave the initial host population enough immunity to slow down the initial virus upswing and place the dynamics closer to their equilibrium state; initial <italic>R </italic>was 1.28. Future work should attempt to more accurately model initial evolutionary dynamics.</p><p>As in the current analysis, previous studies [<xref ref-type="bibr" rid="B4">4</xref>,<xref ref-type="bibr" rid="B6">6</xref>,<xref ref-type="bibr" rid="B7">7</xref>] have assumed that host immune response is dictated by the closest phenotype in the immune history. Due to original antigenic sin, other phenotypes in the host immune history may have a disproportionate effect on host immunity. This may be an important aspect to modeling influenza and should be addressed in future studies. Our model follows that of Lin <italic>et al. </italic>[<xref ref-type="bibr" rid="B29">29</xref>] and Gog and Grenfell [<xref ref-type="bibr" rid="B30">30</xref>] in representing antigenic distance as the distance between points in a geometric space. By forcing one dimension to directly modulate <italic>β</italic>, Gog and Grenfell find an orderly replacement of strains. Kryazhimskiy <italic>et al. </italic>[<xref ref-type="bibr" rid="B31">31</xref>] use a two-dimensional strain space, in which cross-immunity between two strains is proportional to their distance in one dimension or the other, whichever is closer. This cross-immunity kernel directly favors moving along a diagonal line away from previous strains. Our model more closely approximates how HI distance is incorporated into the antigenic map of Smith <italic>et al. </italic>[<xref ref-type="bibr" rid="B9">9</xref>]; HI between strains is projected as a Euclidean distance, rather than as the closest distance between strains in either dimension one or two.</p><p>The initial virus population consisted of 10 infections each with the identical antigenic phenotype of {0, 0}. Over time, viruses evolve in antigenic phenotype. Each day, there is a chance <italic>μ </italic>= 10<sup>-4 </sup>that an infection mutates to a new phenotype. This mutation rate represents a phenotypic rate, rather than genetic mutation rate, and can be thought of as arising from multiple genetic sources. When a mutation occurs, the virus's phenotype is moved in a random radial direction. Mutation size is sampled from a gamma distribution, with distribution parameters chosen to give a mean mutation size of δ<sub>avg </sub>= 0.6 units and a standard deviation of δ<sub>sd </sub>= 0.4 units. This distribution is parameterized so that mutation usually has little effect on antigenic phenotype but occasionally has a larger effect. This is similar to the neutral networks implemented by Koelle <italic>et al. </italic>[<xref ref-type="bibr" rid="B6">6</xref>], wherein most amino acid changes result in little decrease to cross-immunity between strains, but some changes result in large jumps in cross-immunity. This mutation model implicitly assumes independence between mutations; each mutation's effect is independent of genetic background.</p></sec><sec><title>Model output</title><p>Daily incidence and prevalence are recorded for each region. During the course of the simulation, samples of current infections are taken from the evolving virus population at a rate proportional to prevalence. Each viral infection is assigned a unique identification (ID), and in addition, infections have their phenotypes, locations, and dates of infection recorded. In this model, viruses lack sequences and so standard phylogenetic inference of the evolutionary relationships among strains is impossible. Instead, the viral genealogy is directly recorded. This is made possible by tracking transmission events connecting infections during the simulation; infections record the ID of their 'parent' infection. Proceeding from a sample of infections, their genealogical history can be reconstructed by following consecutive links to parental infections. Following this procedure, lineages coalesce to the ancestral lineages shared by the sampled infections, eventually arriving at the initial infection introduced at the beginning of the simulation. Commonly, phylodynamic simulations generate sequences that are subsequently analyzed with a phylogenetic software to produce an estimated genealogy [<xref ref-type="bibr" rid="B4">4</xref>,<xref ref-type="bibr" rid="B6">6</xref>,<xref ref-type="bibr" rid="B32">32</xref>]. This step of phylogenetic inference is imperfect and computationally intensive, and by side-stepping phylogenetic reconstruction, we arrive at genealogies quickly and accurately. Other authors have implemented similar tracking of infection trees [<xref ref-type="bibr" rid="B33">33</xref>,<xref ref-type="bibr" rid="B34">34</xref>]. This genealogy-centric approach makes many otherwise difficult calculations transparent, such as calculating lineage-specific region-specific migration rates (Figure <xref ref-type="fig" rid="F5">5</xref>) and lineage-specific mutation effects (Table <xref ref-type="table" rid="T1">1</xref>).</p><p>Infections are sampled at a rate designed to give approximately 6,000 samples over the course of the simulation, with genealogies constructed from a subsample of approximately 300 samples. The results presented in Figure <xref ref-type="fig" rid="F2">2</xref> represent a single representative model output; 100 replicate simulations were conducted to arrive at statistical estimates.</p></sec><sec><title>Parameter selection and sensitivity analysis</title><p>Estimating what the basic reproductive number <italic>R</italic><sub>0 </sub>for seasonal influenza would be in a naive population is notoriously difficult. Season-to-season estimates of effective reproductive number <italic>R </italic>for the USA and France gathered from mortality time series display an interquartile range of 0.9 to 1.8 [<xref ref-type="bibr" rid="B35">35</xref>]. Geographic spread within the USA suggests an average seasonal <italic>R </italic>of 1.35 [<xref ref-type="bibr" rid="B36">36</xref>]. These estimates of <italic>R </italic>will be lower than the <italic>R</italic><sub>0 </sub>of influenza due to the effects of human immunity. We assumed a <italic>R</italic><sub>0 </sub>of 1.8, consistent with the upper range of seasonal estimates. Duration of infection was chosen based on patterns of viral shedding shown during challenge studies [<xref ref-type="bibr" rid="B37">37</xref>]. The linear form of the risk of infection and its increase as a function of antigenic distance <italic>s </italic>= 0.07 were based on experimental work on equine influenza [<xref ref-type="bibr" rid="B38">38</xref>] and from studies of vaccine effectiveness [<xref ref-type="bibr" rid="B39">39</xref>]. Between-region contact rate <italic>m </italic>was chosen to yield a trunk lineage that resides predominantly in the tropics. With much higher rates of mixing, the trunk lineage ceases to show a preference in the tropics, and with much lower rates of mixing, particular seasons in the north and the south will often be skipped. The amplitude of seasonal forcing <italic>ε </italic>was chosen to be just large enough to get consistent fade-outs in the summer months and is consistent with empirical estimates [<xref ref-type="bibr" rid="B40">40</xref>].</p><p>Mutational parameters were based, in part, on model behavior. We assumed 10 amino acid sites involved in antigenicity, each mutating at a rate of 10<sup>-5 </sup>[<xref ref-type="bibr" rid="B41">41</xref>] to give a phenotypic mutation rate <italic>μ </italic>= 10<sup>-4</sup>per infection per day. We chose mutational effect parameters (<italic>δ</italic><sub>avg </sub>= 0.6, <italic>δ</italic><sub>sd </sub>= 0.4) that would give suitably fast rates of antigenic evolution corresponding to approximately 1.2 units of antigenic change per year while simultaneously giving clustered patterns of antigenic evolution [<xref ref-type="bibr" rid="B9">9</xref>]. Similar outcomes are possible under a variety of parameterizations. If mutations are more common (<italic>μ </italic>= 3 × 10<sup>-4</sup>) and show less variation in effect size (<italic>δ</italic><sub>avg </sub>= 0.6, <italic>δ</italic><sub>sd </sub>= 0.2), then antigenic drift occurs in a more continuous fashion, resulting in less variation in seasonal incidence and a smoother distribution of antigenic phenotypes [see Figure S2A,C in Additional file <xref ref-type="supplementary-material" rid="S1">1</xref>]. If mutations are less common (<italic>μ </italic>= 5 × 10<sup>-5</sup>) and show more variance in effect (<italic>δ</italic><sub>avg </sub>= 0.7, <italic>δ</italic><sub>sd </sub>= 0.5), then antigenic change occurs in a more punctuated fashion [see Figure S2B,D in Additional file <xref ref-type="supplementary-material" rid="S1">1</xref>]. Basic reproductive number <italic>R</italic><sub>0 </sub>can be traded off with mutational parameters to some extent. Less mutational input and higher <italic>R</italic><sub>0 </sub>will give similar patterns of antigenic drift and seasonal incidence. Similarly, Kucharski and Gog [<xref ref-type="bibr" rid="B42">42</xref>] find that increasing <italic>R</italic><sub>0 </sub>results in increased rates of emergence of antigenically novel strains.</p></sec><sec><title>Antigenic map</title><p>Antigenic phenotypes are modeled as discrete entities on the Euclidean plane; multiple samples have the same antigenic location. However, in the empirical antigenic map of influenza A (H3N2), each strain appears in a unique location [<xref ref-type="bibr" rid="B9">9</xref>]. We would argue that some of this pattern comes from experimental noise. Indeed, Smith <italic>et al. </italic>[<xref ref-type="bibr" rid="B9">9</xref>] find that observed measurements and measurements predicted from the map differ by an average of 0.83 antigenic units with a standard deviation of 0.67 antigenic units. We take this as a proxy for experimental noise and add jitter to each sampled antigenic phenotype by moving it in a random direction for an exponentially distributed distance with mean of 0.53 antigenic units. If two samples with the same underlying antigenic phenotype are jittered in this fashion, the distance between them averages 0.83 antigenic units with a standard deviation of 0.64 units.</p><p>We added noise to each of the 5,943 sampled viruses in this fashion, resulting in an approximated antigenic map (Figure <xref ref-type="fig" rid="F2">2D</xref>). Virus samples were then clustered following standard clustering algorithms. We tried clustering by the <italic>k</italic>-means algorithm and also agglomerative hierarchical clustering with a variety of linkage criterion. We found that clustering by Ward's criterion consistently outperformed other methods when measured in terms of within-cluster and between-cluster variances. However, the exact clustering algorithm had little effect on our overall results.</p></sec></sec><sec><title>Competing interests</title><p>The authors declare that they have no competing interests.</p></sec><sec><title>Authors' contributions</title><p>TB and MP conceived the project. TB conducted the simulation experiments and performed statistical analyses. TB, AR, and MP analyzed the results and wrote the paper. All authors have read and approved the final manuscript.</p></sec><sec><title>Authors' information</title><p>TB holds a Long-Term Fellowship from the European Molecular Biology Organization. AR works as a part of the Interdisciplinary Centre for Human and Avian Influenza Research (ICHAIR) and the University of Edinburgh's Centre for Immunity, Infection and Evolution (CIIE). MP is an investigator of the Howard Hughes Medical Institute.</p></sec><sec sec-type="supplementary-material"><title>Supplementary Material</title><supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p><bold>Supplementary figures</bold>. Includes supplementary figures 1 to 3 detailing results from simulations with alternative parameter values.</p></caption><media xlink:href="1741-7007-10-38-S1.PDF" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media></supplementary-material></sec></body><back><sec><title>Acknowledgements</title><p>We would like to thank Sarah Cobey, Aaron King, Pejman Rohani, and the attendees of the 2011 RAPIDD Workshop on Phylodynamics for the helpful discussion. We would also like to thank Ed Baskerville and Daniel Zinder for the programming advice. The term 'canalization' was originally suggested by Micaela Martinez-Bakker.</p></sec><ref-list><ref id="B1"><mixed-citation publication-type="other"><collab>WHO</collab><article-title>Influenza Fact Sheet</article-title><ext-link ext-link-type="uri" xlink:href="http://www.who.int/mediacentre/factsheets/fs211/en/">http://www.who.int/mediacentre/factsheets/fs211/en/</ext-link></mixed-citation></ref><ref id="B2"><mixed-citation publication-type="journal"><name><surname>Nelson</surname><given-names>MI</given-names></name><name><surname>Holmes</surname><given-names>EC</given-names></name><article-title>The evolution of epidemic influenza</article-title><source>Nat Rev Genet</source><year>2007</year><volume>8</volume><issue>3</issue><fpage>196</fpage><lpage>205</lpage><pub-id pub-id-type="doi">10.1038/nrg2053</pub-id><pub-id pub-id-type="pmid">17262054</pub-id></mixed-citation></ref><ref id="B3"><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><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><year>1997</year><volume>94</volume><issue>15</issue><fpage>7712</fpage><lpage>7718</lpage><pub-id pub-id-type="doi">10.1073/pnas.94.15.7712</pub-id><pub-id pub-id-type="pmid">9223253</pub-id></mixed-citation></ref><ref id="B4"><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><article-title>Ecological and immunological determinants of influenza evolution</article-title><source>Nature</source><year>2003</year><volume>422</volume><issue>6930</issue><fpage>428</fpage><lpage>433</lpage><pub-id pub-id-type="doi">10.1038/nature01509</pub-id><pub-id pub-id-type="pmid">12660783</pub-id></mixed-citation></ref><ref id="B5"><mixed-citation publication-type="other"><name><surname>Tria</surname><given-names>F</given-names></name><name><surname>Lässig</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><article-title>A minimal stochastic model for influenza evolution</article-title><source>J Stat Mech</source><year>2005</year><comment>doi:10.1088/1742-5468/2005/07/P07008</comment></mixed-citation></ref><ref id="B6"><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><article-title>Epochal evolution shapes the phylodynamics of interpandemic influenza A (H3N2) in humans</article-title><source>Science</source><year>2006</year><volume>314</volume><issue>5807</issue><fpage>1898</fpage><lpage>1903</lpage><pub-id pub-id-type="doi">10.1126/science.1132745</pub-id><pub-id pub-id-type="pmid">17185596</pub-id></mixed-citation></ref><ref id="B7"><mixed-citation publication-type="journal"><name><surname>Recker</surname><given-names>M</given-names></name><name><surname>Pybus</surname><given-names>O</given-names></name><name><surname>Nee</surname><given-names>S</given-names></name><name><surname>Gupta</surname><given-names>S</given-names></name><article-title>The generation of influenza outbreaks by a network of host immune responses against a limited set of antigenic types</article-title><source>Proc Natl Acad Sci USA</source><year>2007</year><volume>104</volume><issue>18</issue><fpage>7711</fpage><lpage>7716</lpage><pub-id pub-id-type="doi">10.1073/pnas.0702154104</pub-id><pub-id pub-id-type="pmid">17460037</pub-id></mixed-citation></ref><ref id="B8"><mixed-citation publication-type="journal"><name><surname>Hirst</surname><given-names>G</given-names></name><article-title>Studies of antigenic differences among strains of influenza A by means of red cell agglutination</article-title><source>J Exp Med</source><year>1943</year><volume>78</volume><issue>5</issue><fpage>407</fpage><lpage>423</lpage><pub-id pub-id-type="doi">10.1084/jem.78.5.407</pub-id><pub-id pub-id-type="pmid">19871338</pub-id></mixed-citation></ref><ref id="B9"><mixed-citation publication-type="journal"><name><surname>Smith</surname><given-names>DJ</given-names></name><name><surname>Lapedes</surname><given-names>AS</given-names></name><name><surname>de Jong</surname><given-names>JC</given-names></name><name><surname>Bestebroer</surname><given-names>TM</given-names></name><name><surname>Rimmelzwaan</surname><given-names>GF</given-names></name><name><surname>Osterhaus</surname><given-names>ADME</given-names></name><name><surname>Fouch-ier</surname><given-names>RAM</given-names></name><article-title>Mapping the antigenic and genetic evolution of influenza virus</article-title><source>Science</source><year>2004</year><volume>305</volume><issue>5682</issue><fpage>371</fpage><lpage>376</lpage><pub-id pub-id-type="doi">10.1126/science.1097211</pub-id><pub-id pub-id-type="pmid">15218094</pub-id></mixed-citation></ref><ref id="B10"><mixed-citation publication-type="journal"><name><surname>Monto</surname><given-names>A</given-names></name><name><surname>Sullivan</surname><given-names>K</given-names></name><article-title>Acute respiratory illness in the community. Frequency of illness and the agents involved</article-title><source>Epidemiol Infect</source><year>1993</year><volume>110</volume><issue>01</issue><fpage>145</fpage><lpage>160</lpage><pub-id pub-id-type="doi">10.1017/S0950268800050779</pub-id><pub-id pub-id-type="pmid">8432318</pub-id></mixed-citation></ref><ref id="B11"><mixed-citation publication-type="journal"><name><surname>Koelle</surname><given-names>K</given-names></name><name><surname>Kamradt</surname><given-names>M</given-names></name><name><surname>Pascual</surname><given-names>M</given-names></name><article-title>Understanding the dynamics of rapidly evolving pathogens through modeling the tempo of antigenic change: influenza as a case study</article-title><source>Epidemics</source><year>2009</year><volume>1</volume><issue>2</issue><fpage>129</fpage><lpage>137</lpage><pub-id pub-id-type="doi">10.1016/j.epidem.2009.05.003</pub-id><pub-id pub-id-type="pmid">21352760</pub-id></mixed-citation></ref><ref id="B12"><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><article-title>Strength and temp of selection revealed in viral gene genealogies</article-title><source>BMC Evol Biol</source><year>2011</year><volume>11</volume><fpage>220</fpage><pub-id pub-id-type="doi">10.1186/1471-2148-11-220</pub-id><pub-id pub-id-type="pmid">21787390</pub-id></mixed-citation></ref><ref id="B13"><mixed-citation publication-type="journal"><name><surname>Bush</surname><given-names>R</given-names></name><name><surname>Fitch</surname><given-names>W</given-names></name><name><surname>Bender</surname><given-names>C</given-names></name><name><surname>Cox</surname><given-names>N</given-names></name><article-title>Positive selection on the H3 hemagglutinin gene of human influenza virus A</article-title><source>Mol Biol Evol</source><year>1999</year><volume>16</volume><issue>11</issue><fpage>1457</fpage><lpage>1465</lpage><pub-id pub-id-type="doi">10.1093/oxfordjournals.molbev.a026057</pub-id><pub-id pub-id-type="pmid">10555276</pub-id></mixed-citation></ref><ref id="B14"><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><article-title>Long intervals of stasis punctuated by bursts of positive selection in the seasonal evolution of influenza A virus</article-title><source>Biol Direct</source><year>2006</year><volume>1</volume><fpage>34</fpage><pub-id pub-id-type="doi">10.1186/1745-6150-1-34</pub-id><pub-id pub-id-type="pmid">17067369</pub-id></mixed-citation></ref><ref id="B15"><mixed-citation publication-type="journal"><name><surname>Park</surname><given-names>S</given-names></name><name><surname>Krug</surname><given-names>J</given-names></name><article-title>Clonal interference in large populations</article-title><source>Proc Natl Acad Sci USA</source><year>2007</year><volume>104</volume><issue>46</issue><fpage>18135</fpage><lpage>18140</lpage><pub-id pub-id-type="doi">10.1073/pnas.0705778104</pub-id><pub-id pub-id-type="pmid">17984061</pub-id></mixed-citation></ref><ref id="B16"><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><name><surname>Gregory</surname><given-names>V</given-names></name><name><surname>Gust</surname><given-names>ID</given-names></name><name><surname>Hampson</surname><given-names>AW</given-names></name><name><surname>Hay</surname><given-names>AJ</given-names></name><name><surname>Hurt</surname><given-names>AC</given-names></name><name><surname>de Jong</surname><given-names>JC</given-names></name><name><surname>Kelso</surname><given-names>A</given-names></name><name><surname>Klimov</surname><given-names>AI</given-names></name><name><surname>Kageyama</surname><given-names>T</given-names></name><name><surname>Komadina</surname><given-names>N</given-names></name><name><surname>Lapedes</surname><given-names>AS</given-names></name><name><surname>Lin</surname><given-names>YP</given-names></name><name><surname>Mosterin</surname><given-names>A</given-names></name><name><surname>Obuchi</surname><given-names>M</given-names></name><name><surname>Odagiri</surname><given-names>T</given-names></name><name><surname>Osterhaus</surname><given-names>ADME</given-names></name><name><surname>Rimmelzwaan</surname><given-names>GF</given-names></name><name><surname>Shaw</surname><given-names>MW</given-names></name><name><surname>Skepner</surname><given-names>E</given-names></name><name><surname>Stohr</surname><given-names>K</given-names></name><name><surname>Tashiro</surname><given-names>M</given-names></name><name><surname>Fouchier</surname><given-names>RAM</given-names></name><name><surname>Smith</surname><given-names>DJ</given-names></name><article-title>The global circulation of seasonal influenza A (H3N2) viruses</article-title><source>Science</source><year>2008</year><volume>320</volume><issue>5874</issue><fpage>340</fpage><lpage>346</lpage><pub-id pub-id-type="doi">10.1126/science.1154137</pub-id><pub-id pub-id-type="pmid">18420927</pub-id></mixed-citation></ref><ref id="B17"><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><article-title>Global migration dynamics underlie evolution and persistence in human influenza A (H3N2)</article-title><source>PLoS Pathog</source><year>2010</year><volume>6</volume><fpage>e1000918</fpage><pub-id pub-id-type="doi">10.1371/journal.ppat.1000918</pub-id><pub-id pub-id-type="pmid">20523898</pub-id></mixed-citation></ref><ref id="B18"><mixed-citation publication-type="journal"><name><surname>Adams</surname><given-names>B</given-names></name><name><surname>McHardy</surname><given-names>A</given-names></name><article-title>The impact of seasonal and year-round transmission regimes on the evolution of influenza A virus</article-title><source>Proc R Soc B</source><year>2011</year><volume>278</volume><issue>1716</issue><fpage>2249</fpage><lpage>2256</lpage><pub-id pub-id-type="doi">10.1098/rspb.2010.2191</pub-id></mixed-citation></ref><ref id="B19"><mixed-citation publication-type="other"><name><surname>Koelle</surname><given-names>K</given-names></name><name><surname>Ratmann</surname><given-names>O</given-names></name><name><surname>Rasmussen</surname><given-names>D</given-names></name><name><surname>Pasour</surname><given-names>V</given-names></name><name><surname>Mattingly</surname><given-names>J</given-names></name><article-title>A dimensionless number for under-standing the evolutionary dynamics of antigenically variable RNA viruses</article-title><source>Proc R Soc Lond B Biol Sci</source><year>2011</year><comment>doi: 10.1098/rspb.2011.0435</comment></mixed-citation></ref><ref id="B20"><mixed-citation publication-type="journal"><name><surname>Finkelman</surname><given-names>B</given-names></name><name><surname>Viboud</surname><given-names>C</given-names></name><name><surname>Koelle</surname><given-names>K</given-names></name><name><surname>Ferrari</surname><given-names>M</given-names></name><name><surname>Bharti</surname><given-names>N</given-names></name><name><surname>Grenfell</surname><given-names>B</given-names></name><article-title>Global patterns in seasonal activity of influenza A/H3N2, A/H1N1, and B from 1997 to 2005: viral coexistence and latitudinal gradients</article-title><source>PLoS One</source><year>2007</year><volume>2</volume><issue>12</issue><fpage>e1296</fpage><pub-id pub-id-type="doi">10.1371/journal.pone.0001296</pub-id><pub-id pub-id-type="pmid">18074020</pub-id></mixed-citation></ref><ref id="B21"><mixed-citation publication-type="other"><name><surname>Bhatt</surname><given-names>S</given-names></name><name><surname>Holmes</surname><given-names>E</given-names></name><name><surname>Pybus</surname><given-names>O</given-names></name><article-title>The genomic rate of molecular adaptation of the human influenza A virus</article-title><source>Mol Biol Evol</source><year>2011</year><comment>doi: 10.1093/molbev/msr044</comment></mixed-citation></ref><ref id="B22"><mixed-citation publication-type="journal"><name><surname>Goldstein</surname><given-names>E</given-names></name><name><surname>Cobey</surname><given-names>S</given-names></name><name><surname>Takahashi</surname><given-names>S</given-names></name><name><surname>Miller</surname><given-names>J</given-names></name><name><surname>Lipsitch</surname><given-names>M</given-names></name><article-title>Predicting the epidemic sizes of influenza A/H1N1, A/H3N2, and B: a statistical method</article-title><source>PLoS Med</source><year>2011</year><volume>8</volume><issue>7</issue><fpage>e1001051</fpage><pub-id pub-id-type="doi">10.1371/journal.pmed.1001051</pub-id><pub-id pub-id-type="pmid">21750666</pub-id></mixed-citation></ref><ref id="B23"><mixed-citation publication-type="journal"><name><surname>Voeten</surname><given-names>J</given-names></name><name><surname>Bestebroer</surname><given-names>T</given-names></name><name><surname>Nieuwkoop</surname><given-names>N</given-names></name><name><surname>Fouchier</surname><given-names>R</given-names></name><name><surname>Osterhaus</surname><given-names>A</given-names></name><name><surname>Rimmelzwaan</surname><given-names>G</given-names></name><article-title>Antigenic drift in the influenza A virus (H3N2) nucleoprotein and escape from recognition by cytotoxic T lymphocytes</article-title><source>J Virol</source><year>2000</year><volume>74</volume><issue>15</issue><fpage>6800</fpage><lpage>6807</lpage><pub-id pub-id-type="doi">10.1128/JVI.74.15.6800-6807.2000</pub-id><pub-id pub-id-type="pmid">10888619</pub-id></mixed-citation></ref><ref id="B24"><mixed-citation publication-type="journal"><name><surname>O'Brien</surname><given-names>J</given-names></name><name><surname>Minin</surname><given-names>V</given-names></name><name><surname>Suchard</surname><given-names>M</given-names></name><article-title>Learning to count: robust estimates for labeled distances between molecular sequences</article-title><source>Mol Biol Evol</source><year>2009</year><volume>26</volume><issue>4</issue><fpage>801</fpage><lpage>814</lpage><pub-id pub-id-type="doi">10.1093/molbev/msp003</pub-id><pub-id pub-id-type="pmid">19131426</pub-id></mixed-citation></ref><ref id="B25"><mixed-citation publication-type="journal"><name><surname>Yang</surname><given-names>Z</given-names></name><article-title>Maximum likelihood estimation on large phylogenies and analysis of adaptive evolution in human influenza virus A</article-title><source>Journal of Molecular Evolution</source><year>2000</year><volume>51</volume><issue>5</issue><fpage>423</fpage><lpage>432</lpage><pub-id pub-id-type="pmid">11080365</pub-id></mixed-citation></ref><ref id="B26"><mixed-citation publication-type="book"><name><surname>Gould</surname><given-names>S</given-names></name><source>Wonderful Life</source><year>1989</year><publisher-name>New York: Norton</publisher-name></mixed-citation></ref><ref id="B27"><mixed-citation publication-type="journal"><name><surname>Barr</surname><given-names>IG</given-names></name><name><surname>McCauley</surname><given-names>J</given-names></name><name><surname>Cox</surname><given-names>N</given-names></name><name><surname>Daniels</surname><given-names>R</given-names></name><name><surname>Engelhardt</surname><given-names>OG</given-names></name><name><surname>Fukuda</surname><given-names>K</given-names></name><name><surname>Grohmann</surname><given-names>G</given-names></name><name><surname>Hay</surname><given-names>A</given-names></name><name><surname>Kelso</surname><given-names>A</given-names></name><name><surname>Klimov</surname><given-names>A</given-names></name><name><surname>Odagiri</surname><given-names>T</given-names></name><name><surname>Smith</surname><given-names>D</given-names></name><name><surname>Russell</surname><given-names>C</given-names></name><name><surname>Tashiro</surname><given-names>M</given-names></name><name><surname>Webby</surname><given-names>R</given-names></name><name><surname>Wood</surname><given-names>J</given-names></name><name><surname>Ye</surname><given-names>Z</given-names></name><name><surname>Zhang</surname><given-names>W</given-names></name><collab>Writing Committee of the World Health Organization Consultation on Northern Hemisphere Influenza Vaccine Composition for 2009-2010</collab><article-title>Epidemiological, antigenic and genetic characteristics of seasonal influenza A (H1N1), A (H3N2) and B influenza viruses: basis for the WHO recommendation on the composition of influenza vaccines for use in the 2009-2010 Northern Hemisphere season</article-title><source>Vaccine</source><year>2010</year><volume>28</volume><issue>5</issue><fpage>1156</fpage><lpage>1167</lpage><pub-id pub-id-type="doi">10.1016/j.vaccine.2009.11.043</pub-id><pub-id pub-id-type="pmid">20004635</pub-id></mixed-citation></ref><ref id="B28"><mixed-citation publication-type="journal"><name><surname>Weinreich</surname><given-names>D</given-names></name><name><surname>Delaney</surname><given-names>N</given-names></name><name><surname>DePristo</surname><given-names>M</given-names></name><name><surname>Hartl</surname><given-names>D</given-names></name><article-title>Darwinian evolution can follow only very few mutational paths to fitter proteins</article-title><source>Science</source><year>2006</year><volume>312</volume><issue>5770</issue><fpage>111</fpage><lpage>114</lpage><pub-id pub-id-type="doi">10.1126/science.1123539</pub-id><pub-id pub-id-type="pmid">16601193</pub-id></mixed-citation></ref><ref id="B29"><mixed-citation publication-type="journal"><name><surname>Lin</surname><given-names>J</given-names></name><name><surname>Andreasen</surname><given-names>V</given-names></name><name><surname>Levin</surname><given-names>S</given-names></name><article-title>Dynamics of influenza A drift: the linear three-strain model</article-title><source>Math Biosci</source><year>1999</year><volume>162</volume><issue>1-2</issue><fpage>33</fpage><lpage>51</lpage><pub-id pub-id-type="doi">10.1016/S0025-5564(99)00042-5</pub-id><pub-id pub-id-type="pmid">10616279</pub-id></mixed-citation></ref><ref id="B30"><mixed-citation publication-type="journal"><name><surname>Gog</surname><given-names>J</given-names></name><name><surname>Grenfell</surname><given-names>B</given-names></name><article-title>Dynamics and selection of many-strain pathogens</article-title><source>Proc Natl Acad Sci USA</source><year>2002</year><volume>99</volume><issue>26</issue><fpage>17209</fpage><lpage>17214</lpage><pub-id pub-id-type="doi">10.1073/pnas.252512799</pub-id><pub-id pub-id-type="pmid">12481034</pub-id></mixed-citation></ref><ref id="B31"><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>S</given-names></name><name><surname>Dushoff</surname><given-names>J</given-names></name><article-title>On state-space reduction in multi-strain pathogen models, with an application to antigenic drift in influenza A</article-title><source>PLoS Comp Biol</source><year>2007</year><volume>3</volume><issue>8</issue><fpage>e159</fpage><pub-id pub-id-type="doi">10.1371/journal.pcbi.0030159</pub-id></mixed-citation></ref><ref id="B32"><mixed-citation publication-type="journal"><name><surname>Koelle</surname><given-names>K</given-names></name><name><surname>Khatri</surname><given-names>P</given-names></name><name><surname>Kamradt</surname><given-names>M</given-names></name><name><surname>Kepler</surname><given-names>T</given-names></name><article-title>A two-tiered model for simulating the ecological and evolutionary dynamics of rapidly evolving viruses, with an application to influenza</article-title><source>J R Soc Interface</source><year>2010</year><volume>7</volume><issue>50</issue><fpage>1257</fpage><lpage>1274</lpage><pub-id pub-id-type="doi">10.1098/rsif.2010.0007</pub-id><pub-id pub-id-type="pmid">20335193</pub-id></mixed-citation></ref><ref id="B33"><mixed-citation publication-type="journal"><name><surname>Volz</surname><given-names>E</given-names></name><name><surname>Kosakovsky Pond</surname><given-names>S</given-names></name><name><surname>Ward</surname><given-names>M</given-names></name><name><surname>Leigh Brown</surname><given-names>A</given-names></name><name><surname>Frost</surname><given-names>S</given-names></name><article-title>Phylodynamics of infectious disease epidemics</article-title><source>Genetics</source><year>2009</year><volume>183</volume><issue>4</issue><fpage>1421</fpage><lpage>1430</lpage><pub-id pub-id-type="doi">10.1534/genetics.109.106021</pub-id><pub-id pub-id-type="pmid">19797047</pub-id></mixed-citation></ref><ref id="B34"><mixed-citation publication-type="journal"><name><surname>O'Dea</surname><given-names>E</given-names></name><name><surname>Wilke</surname><given-names>C</given-names></name><article-title>Contact heterogeneity and phylodynamics: how contact networks shape parasite evolutionary trees</article-title><source>Interdiscip Perspect Infect Dis</source><year>2011</year><volume>2011</volume><fpage>238743</fpage><pub-id pub-id-type="pmid">21151699</pub-id></mixed-citation></ref><ref id="B35"><mixed-citation publication-type="journal"><name><surname>Chowell</surname><given-names>G</given-names></name><name><surname>Miller</surname><given-names>M</given-names></name><name><surname>Viboud</surname><given-names>C</given-names></name><article-title>Seasonal influenza in the United States, France, and Australia: transmission and prospects for control</article-title><source>Epidemiol Infect</source><year>2008</year><volume>136</volume><issue>06</issue><fpage>852</fpage><lpage>864</lpage><pub-id pub-id-type="pmid">17634159</pub-id></mixed-citation></ref><ref id="B36"><mixed-citation publication-type="journal"><name><surname>Viboud</surname><given-names>C</given-names></name><name><surname>Bjørnstad</surname><given-names>O</given-names></name><name><surname>Smith</surname><given-names>D</given-names></name><name><surname>Simonsen</surname><given-names>L</given-names></name><name><surname>Miller</surname><given-names>M</given-names></name><name><surname>Grenfell</surname><given-names>B</given-names></name><article-title>Synchrony, waves, and spatial hierarchies in the spread of influenza</article-title><source>Science</source><year>2006</year><volume>312</volume><issue>5772</issue><fpage>447</fpage><lpage>451</lpage><pub-id pub-id-type="doi">10.1126/science.1125237</pub-id><pub-id pub-id-type="pmid">16574822</pub-id></mixed-citation></ref><ref id="B37"><mixed-citation publication-type="journal"><name><surname>Carrat</surname><given-names>F</given-names></name><name><surname>Vergu</surname><given-names>E</given-names></name><name><surname>Ferguson</surname><given-names>NM</given-names></name><name><surname>Lemaitre</surname><given-names>M</given-names></name><name><surname>Cauchemez</surname><given-names>S</given-names></name><name><surname>Leach</surname><given-names>S</given-names></name><name><surname>Valleron</surname><given-names>AJ</given-names></name><article-title>Time lines of infection and disease in human influenza: a review of volunteer challenge studies</article-title><source>Am J Epidemiol</source><year>2008</year><volume>167</volume><issue>7</issue><fpage>775</fpage><lpage>785</lpage><pub-id pub-id-type="doi">10.1093/aje/kwm375</pub-id><pub-id pub-id-type="pmid">18230677</pub-id></mixed-citation></ref><ref id="B38"><mixed-citation publication-type="journal"><name><surname>Park</surname><given-names>A</given-names></name><name><surname>Daly</surname><given-names>J</given-names></name><name><surname>Lewis</surname><given-names>N</given-names></name><name><surname>Smith</surname><given-names>D</given-names></name><name><surname>Wood</surname><given-names>J</given-names></name><name><surname>Grenfell</surname><given-names>B</given-names></name><article-title>Quantifying the impact of immune escape on transmission dynamics of influenza</article-title><source>Science</source><year>2009</year><volume>326</volume><issue>5953</issue><fpage>726</fpage><lpage>728</lpage><pub-id pub-id-type="doi">10.1126/science.1175980</pub-id><pub-id pub-id-type="pmid">19900931</pub-id></mixed-citation></ref><ref id="B39"><mixed-citation publication-type="journal"><name><surname>Gupta</surname><given-names>V</given-names></name><name><surname>Earl</surname><given-names>D</given-names></name><name><surname>Deem</surname><given-names>M</given-names></name><article-title>Quantifying influenza vaccine efficacy and antigenic distance</article-title><source>Vaccine</source><year>2006</year><volume>24</volume><issue>18</issue><fpage>3881</fpage><lpage>3888</lpage><pub-id pub-id-type="doi">10.1016/j.vaccine.2006.01.010</pub-id><pub-id pub-id-type="pmid">16460844</pub-id></mixed-citation></ref><ref id="B40"><mixed-citation publication-type="other"><name><surname>Truscott</surname><given-names>J</given-names></name><name><surname>Fraser</surname><given-names>C</given-names></name><name><surname>Cauchemez</surname><given-names>S</given-names></name><name><surname>Meeyai</surname><given-names>A</given-names></name><name><surname>Hinsley</surname><given-names>W</given-names></name><name><surname>Donnelly</surname><given-names>C</given-names></name><name><surname>Ghani</surname><given-names>A</given-names></name><name><surname>Ferguson</surname><given-names>N</given-names></name><article-title>Essential epidemiological mechanisms underpinning the transmission dynamics of seasonal influenza</article-title><source>J R Soc Interface</source><year>2011</year><comment>doi: 10.1098/rsif.2011.0309</comment></mixed-citation></ref><ref id="B41"><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><name><surname>Holmes</surname><given-names>EC</given-names></name><article-title>The genomic and epidemiological dynamics of human influenza A virus</article-title><source>Nature</source><year>2008</year><volume>453</volume><issue>7195</issue><fpage>615</fpage><lpage>619</lpage><pub-id pub-id-type="doi">10.1038/nature06945</pub-id><pub-id pub-id-type="pmid">18418375</pub-id></mixed-citation></ref><ref id="B42"><mixed-citation publication-type="journal"><name><surname>Kucharski</surname><given-names>A</given-names></name><name><surname>Gog</surname><given-names>JR</given-names></name><article-title>Influenza emergence in the face of evolutionary constraints</article-title><source>Proc R Soc B</source><year>2012</year><volume>279</volume><issue>1729</issue><fpage>645</fpage><lpage>652</lpage><pub-id pub-id-type="doi">10.1098/rspb.2011.1168</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 000904 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000904 | 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é= PMC:3373370
|texte= Canalization of the evolutionary trajectory of the human influenza virus
}}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:22546494" \
| HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \
| NlmPubMed2Wicri -a H2N2V1
| This area was generated with Dilib version V0.6.33. |  |