Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">A Perspective on Multiple Waves of Influenza Pandemics</title>
<author><name sortKey="Mummert, Anna" sort="Mummert, Anna" uniqKey="Mummert A" first="Anna" last="Mummert">Anna Mummert</name>
<affiliation><nlm:aff id="aff1"><addr-line>Department of Mathematics, Marshall University, Huntington, West Virginia, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Weiss, Howard" sort="Weiss, Howard" uniqKey="Weiss H" first="Howard" last="Weiss">Howard Weiss</name>
<affiliation><nlm:aff id="aff2"><addr-line>School of Mathematics, Georgia Institute of Technology, Atlanta, Georgia, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Long, Li Ping" sort="Long, Li Ping" uniqKey="Long L" first="Li-Ping" last="Long">Li-Ping Long</name>
<affiliation><nlm:aff id="aff3"><addr-line>Department of Basic Sciences, Mississippi State University, Mississippi State, Mississippi, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Amig, Jose M" sort="Amig, Jose M" uniqKey="Amig J" first="José M." last="Amig">José M. Amig</name>
<affiliation><nlm:aff id="aff4"><addr-line>Centro de Investigación Operativa, Universidad Miguel Hernández, Elche, Spain</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Wan, Xiu Feng" sort="Wan, Xiu Feng" uniqKey="Wan X" first="Xiu-Feng" last="Wan">Xiu-Feng Wan</name>
<affiliation><nlm:aff id="aff3"><addr-line>Department of Basic Sciences, Mississippi State University, Mississippi State, Mississippi, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">23637746</idno>
<idno type="pmc">3634039</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3634039</idno>
<idno type="RBID">PMC:3634039</idno>
<idno type="doi">10.1371/journal.pone.0060343</idno>
<date when="2013">2013</date>
<idno type="wicri:Area/Pmc/Corpus">000B60</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000B60</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">A Perspective on Multiple Waves of Influenza Pandemics</title>
<author><name sortKey="Mummert, Anna" sort="Mummert, Anna" uniqKey="Mummert A" first="Anna" last="Mummert">Anna Mummert</name>
<affiliation><nlm:aff id="aff1"><addr-line>Department of Mathematics, Marshall University, Huntington, West Virginia, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Weiss, Howard" sort="Weiss, Howard" uniqKey="Weiss H" first="Howard" last="Weiss">Howard Weiss</name>
<affiliation><nlm:aff id="aff2"><addr-line>School of Mathematics, Georgia Institute of Technology, Atlanta, Georgia, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Long, Li Ping" sort="Long, Li Ping" uniqKey="Long L" first="Li-Ping" last="Long">Li-Ping Long</name>
<affiliation><nlm:aff id="aff3"><addr-line>Department of Basic Sciences, Mississippi State University, Mississippi State, Mississippi, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Amig, Jose M" sort="Amig, Jose M" uniqKey="Amig J" first="José M." last="Amig">José M. Amig</name>
<affiliation><nlm:aff id="aff4"><addr-line>Centro de Investigación Operativa, Universidad Miguel Hernández, Elche, Spain</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Wan, Xiu Feng" sort="Wan, Xiu Feng" uniqKey="Wan X" first="Xiu-Feng" last="Wan">Xiu-Feng Wan</name>
<affiliation><nlm:aff id="aff3"><addr-line>Department of Basic Sciences, Mississippi State University, Mississippi State, Mississippi, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">PLoS ONE</title>
<idno type="eISSN">1932-6203</idno>
<imprint><date when="2013">2013</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><sec><title>Background</title>
<p>A striking characteristic of the past four influenza pandemic outbreaks in the United States has been the multiple waves of infections. However, the mechanisms responsible for the multiple waves of influenza or other acute infectious diseases are uncertain. Understanding these mechanisms could provide knowledge for health authorities to develop and implement prevention and control strategies.</p>
</sec>
<sec><title>Materials and Methods</title>
<p>We exhibit five distinct mechanisms, each of which can generate two waves of infections for an acute infectious disease. The first two mechanisms capture changes in virus transmissibility and behavioral changes. The third mechanism involves population heterogeneity (e.g., demography, geography), where each wave spreads through one sub-population. The fourth mechanism is virus mutation which causes delayed susceptibility of individuals. The fifth mechanism is waning immunity. Each mechanism is incorporated into separate mathematical models, and outbreaks are then simulated. We use the models to examine the effects of the initial number of infected individuals (e.g., border control at the beginning of the outbreak) and the timing of and amount of available vaccinations.</p>
</sec>
<sec><title>Results</title>
<p>Four models, individually or in any combination, reproduce the two waves of the 2009 H1N1 pandemic in the United States, both qualitatively and quantitatively. One model reproduces the two waves only qualitatively. All models indicate that significantly reducing or delaying the initial numbers of infected individuals would have little impact on the attack rate. Instead, this reduction or delay results in a single wave as opposed to two waves. Furthermore, four of these models also indicate that a vaccination program started earlier than October 2009 (when the H1N1 vaccine was initially distributed) could have eliminated the second wave of infection, while more vaccine available starting in October would not have eliminated the second wave.</p>
</sec>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Miller, Ma" uniqKey="Miller M">MA Miller</name>
</author>
<author><name sortKey="Viboud, C" uniqKey="Viboud C">C Viboud</name>
</author>
<author><name sortKey="Balinska, M" uniqKey="Balinska M">M Balinska</name>
</author>
<author><name sortKey="Simonsen, L" uniqKey="Simonsen L">L Simonsen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chao, Dl" uniqKey="Chao D">DL Chao</name>
</author>
<author><name sortKey="Halloran, Me" uniqKey="Halloran M">ME Halloran</name>
</author>
<author><name sortKey="Longini, Im" uniqKey="Longini I">IM Longini</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chowell, G" uniqKey="Chowell G">G Chowell</name>
</author>
<author><name sortKey="Echevarria Zuno, S" uniqKey="Echevarria Zuno S">S Echevarria-Zuno</name>
</author>
<author><name sortKey="Viboud, C" uniqKey="Viboud C">C Viboud</name>
</author>
<author><name sortKey="Simonsen, L" uniqKey="Simonsen L">L Simonsen</name>
</author>
<author><name sortKey="Tamerius, J" uniqKey="Tamerius J">J Tamerius</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Monto, As" uniqKey="Monto A">AS Monto</name>
</author>
<author><name sortKey="Koopman, Js" uniqKey="Koopman J">JS Koopman</name>
</author>
<author><name sortKey="Longini, Im" uniqKey="Longini I">IM Longini</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cauchemez, S" uniqKey="Cauchemez S">S Cauchemez</name>
</author>
<author><name sortKey="Ferguson, Nm" uniqKey="Ferguson N">NM Ferguson</name>
</author>
<author><name sortKey="Wachtel, C" uniqKey="Wachtel C">C Wachtel</name>
</author>
<author><name sortKey="Tegnell, A" uniqKey="Tegnell A">A Tegnell</name>
</author>
<author><name sortKey="Saour, G" uniqKey="Saour G">G Saour</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cauchemez, S" uniqKey="Cauchemez S">S Cauchemez</name>
</author>
<author><name sortKey="Valleron, Aj" uniqKey="Valleron A">AJ Valleron</name>
</author>
<author><name sortKey="Boelle, Py" uniqKey="Boelle P">PY Boelle</name>
</author>
<author><name sortKey="Flahault, A" uniqKey="Flahault A">A Flahault</name>
</author>
<author><name sortKey="Ferguson, Nm" uniqKey="Ferguson N">NM Ferguson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Melidou, A" uniqKey="Melidou A">A Melidou</name>
</author>
<author><name sortKey="Gioula, G" uniqKey="Gioula G">G Gioula</name>
</author>
<author><name sortKey="Exindari, M" uniqKey="Exindari M">M Exindari</name>
</author>
<author><name sortKey="Chatzidimitriou, D" uniqKey="Chatzidimitriou D">D Chatzidimitriou</name>
</author>
<author><name sortKey="Diza, E" uniqKey="Diza E">E Diza</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chan, Pk" uniqKey="Chan P">PK Chan</name>
</author>
<author><name sortKey="Lee, N" uniqKey="Lee N">N Lee</name>
</author>
<author><name sortKey="Joynt, Gm" uniqKey="Joynt G">GM Joynt</name>
</author>
<author><name sortKey="Choi, Kw" uniqKey="Choi K">KW Choi</name>
</author>
<author><name sortKey="Cheung, Jl" uniqKey="Cheung J">JL Cheung</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pan, C" uniqKey="Pan C">C Pan</name>
</author>
<author><name sortKey="Cheung, B" uniqKey="Cheung B">B Cheung</name>
</author>
<author><name sortKey="Tan, S" uniqKey="Tan S">S Tan</name>
</author>
<author><name sortKey="Li, C" uniqKey="Li C">C Li</name>
</author>
<author><name sortKey="Li, L" uniqKey="Li L">L Li</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Morlighem, Je" uniqKey="Morlighem J">JE Morlighem</name>
</author>
<author><name sortKey="Aoki, S" uniqKey="Aoki S">S Aoki</name>
</author>
<author><name sortKey="Kishima, M" uniqKey="Kishima M">M Kishima</name>
</author>
<author><name sortKey="Hanami, M" uniqKey="Hanami M">M Hanami</name>
</author>
<author><name sortKey="Ogawa, C" uniqKey="Ogawa C">C Ogawa</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wang, B" uniqKey="Wang B">B Wang</name>
</author>
<author><name sortKey="Dwyer, De" uniqKey="Dwyer D">DE Dwyer</name>
</author>
<author><name sortKey="Soedjono, M" uniqKey="Soedjono M">M Soedjono</name>
</author>
<author><name sortKey="Shi, H" uniqKey="Shi H">H Shi</name>
</author>
<author><name sortKey="Matlho, K" uniqKey="Matlho K">K Matlho</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Drews, Sj" uniqKey="Drews S">SJ Drews</name>
</author>
<author><name sortKey="Pabbaraju, K" uniqKey="Pabbaraju K">K Pabbaraju</name>
</author>
<author><name sortKey="Wong, S" uniqKey="Wong S">S Wong</name>
</author>
<author><name sortKey="Tokaryk, Kl" uniqKey="Tokaryk K">KL Tokaryk</name>
</author>
<author><name sortKey="May Hadford, J" uniqKey="May Hadford J">J May-Hadford</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Miller, Rr" uniqKey="Miller R">RR Miller</name>
</author>
<author><name sortKey="Maclean, Ar" uniqKey="Maclean A">AR MacLean</name>
</author>
<author><name sortKey="Gunson, Rn" uniqKey="Gunson R">RN Gunson</name>
</author>
<author><name sortKey="Carman, Wf" uniqKey="Carman W">WF Carman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Valli, Mb" uniqKey="Valli M">MB Valli</name>
</author>
<author><name sortKey="Selleri, M" uniqKey="Selleri M">M Selleri</name>
</author>
<author><name sortKey="Meschi, S" uniqKey="Meschi S">S Meschi</name>
</author>
<author><name sortKey="Zaccaro, P" uniqKey="Zaccaro P">P Zaccaro</name>
</author>
<author><name sortKey="Vincenti, D" uniqKey="Vincenti D">D Vincenti</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tse, H" uniqKey="Tse H">H Tse</name>
</author>
<author><name sortKey="To, Kk" uniqKey="To K">KK To</name>
</author>
<author><name sortKey="Wen, X" uniqKey="Wen X">X Wen</name>
</author>
<author><name sortKey="Chen, H" uniqKey="Chen H">H Chen</name>
</author>
<author><name sortKey="Chan, Kh" uniqKey="Chan K">KH Chan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Puzelli, S" uniqKey="Puzelli S">S Puzelli</name>
</author>
<author><name sortKey="Facchini, M" uniqKey="Facchini M">M Facchini</name>
</author>
<author><name sortKey="Spagnolo, D" uniqKey="Spagnolo D">D Spagnolo</name>
</author>
<author><name sortKey="De Marco, Ma" uniqKey="De Marco M">MA De Marco</name>
</author>
<author><name sortKey="Calzoletti, L" uniqKey="Calzoletti L">L Calzoletti</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mak, Gc" uniqKey="Mak G">GC Mak</name>
</author>
<author><name sortKey="Au, Kw" uniqKey="Au K">KW Au</name>
</author>
<author><name sortKey="Tai, Ls" uniqKey="Tai L">LS Tai</name>
</author>
<author><name sortKey="Chuang, Kc" uniqKey="Chuang K">KC Chuang</name>
</author>
<author><name sortKey="Cheng, Kc" uniqKey="Cheng K">KC Cheng</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kilander, A" uniqKey="Kilander A">A Kilander</name>
</author>
<author><name sortKey="Rykkvin, R" uniqKey="Rykkvin R">R Rykkvin</name>
</author>
<author><name sortKey="Dudman, Sg" uniqKey="Dudman S">SG Dudman</name>
</author>
<author><name sortKey="Hungnes, O" uniqKey="Hungnes O">O Hungnes</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Belser, Ja" uniqKey="Belser J">JA Belser</name>
</author>
<author><name sortKey="Jayaraman, A" uniqKey="Jayaraman A">A Jayaraman</name>
</author>
<author><name sortKey="Raman, R" uniqKey="Raman R">R Raman</name>
</author>
<author><name sortKey="Pappas, C" uniqKey="Pappas C">C Pappas</name>
</author>
<author><name sortKey="Zeng, H" uniqKey="Zeng H">H Zeng</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Abed, Y" uniqKey="Abed Y">Y Abed</name>
</author>
<author><name sortKey="Pizzorno, A" uniqKey="Pizzorno A">A Pizzorno</name>
</author>
<author><name sortKey="Hamelin, Me" uniqKey="Hamelin M">ME Hamelin</name>
</author>
<author><name sortKey="Leung, A" uniqKey="Leung A">A Leung</name>
</author>
<author><name sortKey="Joubert, P" uniqKey="Joubert P">P Joubert</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Liu, Y" uniqKey="Liu Y">Y Liu</name>
</author>
<author><name sortKey="Childs, Ra" uniqKey="Childs R">RA Childs</name>
</author>
<author><name sortKey="Matrosovich, T" uniqKey="Matrosovich T">T Matrosovich</name>
</author>
<author><name sortKey="Wharton, S" uniqKey="Wharton S">S Wharton</name>
</author>
<author><name sortKey="Palma, As" uniqKey="Palma A">AS Palma</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wang, W" uniqKey="Wang W">W Wang</name>
</author>
<author><name sortKey="Castelan Vega, Ja" uniqKey="Castelan Vega J">JA Castelan-Vega</name>
</author>
<author><name sortKey="Jimenez Alberto, A" uniqKey="Jimenez Alberto A">A Jimenez-Alberto</name>
</author>
<author><name sortKey="Vassell, R" uniqKey="Vassell R">R Vassell</name>
</author>
<author><name sortKey="Ye, Z" uniqKey="Ye Z">Z Ye</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zheng, B" uniqKey="Zheng B">B Zheng</name>
</author>
<author><name sortKey="Chan, Kh" uniqKey="Chan K">KH Chan</name>
</author>
<author><name sortKey="Zhang, Aj" uniqKey="Zhang A">AJ Zhang</name>
</author>
<author><name sortKey="Zhou, J" uniqKey="Zhou J">J Zhou</name>
</author>
<author><name sortKey="Chan, Cc" uniqKey="Chan C">CC Chan</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Hadeler, Kp" uniqKey="Hadeler K">KP Hadeler</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mummert, A" uniqKey="Mummert A">A Mummert</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Greenberg, Me" uniqKey="Greenberg M">ME Greenberg</name>
</author>
<author><name sortKey="Lai, Mh" uniqKey="Lai M">MH Lai</name>
</author>
<author><name sortKey="Hartel, Gf" uniqKey="Hartel G">GF Hartel</name>
</author>
<author><name sortKey="Wichems, Ch" uniqKey="Wichems C">CH Wichems</name>
</author>
<author><name sortKey="Gittleson, C" uniqKey="Gittleson C">C Gittleson</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Jackson, C" uniqKey="Jackson C">C Jackson</name>
</author>
<author><name sortKey="Vynnycky, E" uniqKey="Vynnycky E">E Vynnycky</name>
</author>
<author><name sortKey="Mangtani, P" uniqKey="Mangtani P">P Mangtani</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wood, Jg" uniqKey="Wood J">JG Wood</name>
</author>
<author><name sortKey="Zamani, N" uniqKey="Zamani N">N Zamani</name>
</author>
<author><name sortKey="Macintyre, Cr" uniqKey="Macintyre C">CR MacIntyre</name>
</author>
<author><name sortKey="Beckert, Ng" uniqKey="Beckert N">NG Beckert</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bishop, Jf" uniqKey="Bishop J">JF Bishop</name>
</author>
<author><name sortKey="Murnane, Mp" uniqKey="Murnane M">MP Murnane</name>
</author>
<author><name sortKey="Owen, R" uniqKey="Owen R">R Owen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cowling, Bj" uniqKey="Cowling B">BJ Cowling</name>
</author>
<author><name sortKey="Lau, Ll" uniqKey="Lau L">LL Lau</name>
</author>
<author><name sortKey="Wu, P" uniqKey="Wu P">P Wu</name>
</author>
<author><name sortKey="Wong, Hw" uniqKey="Wong H">HW Wong</name>
</author>
<author><name sortKey="Fang, Vj" uniqKey="Fang V">VJ Fang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Miller, E" uniqKey="Miller E">E Miller</name>
</author>
<author><name sortKey="Hoschler, K" uniqKey="Hoschler K">K Hoschler</name>
</author>
<author><name sortKey="Hardelid, P" uniqKey="Hardelid P">P Hardelid</name>
</author>
<author><name sortKey="Stanford, E" uniqKey="Stanford E">E Stanford</name>
</author>
<author><name sortKey="Andrews, N" uniqKey="Andrews N">N Andrews</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Bao, Y" uniqKey="Bao Y">Y Bao</name>
</author>
<author><name sortKey="Bolotov, P" uniqKey="Bolotov P">P Bolotov</name>
</author>
<author><name sortKey="Dernovoy, D" uniqKey="Dernovoy D">D Dernovoy</name>
</author>
<author><name sortKey="Kiryutin, B" uniqKey="Kiryutin B">B Kiryutin</name>
</author>
<author><name sortKey="Zaslavsky, L" uniqKey="Zaslavsky L">L Zaslavsky</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Long, Lp" uniqKey="Long L">LP Long</name>
</author>
<author><name sortKey="Yuan, C" uniqKey="Yuan C">C Yuan</name>
</author>
<author><name sortKey="Cai, Z" uniqKey="Cai Z">Z Cai</name>
</author>
<author><name sortKey="Xu, H" uniqKey="Xu H">H Xu</name>
</author>
<author><name sortKey="Wan, X F" uniqKey="Wan X">X-F Wan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wan, Xf" uniqKey="Wan X">XF Wan</name>
</author>
<author><name sortKey="Wu, X" uniqKey="Wu X">X Wu</name>
</author>
<author><name sortKey="Lin, G" uniqKey="Lin G">G Lin</name>
</author>
<author><name sortKey="Holton, Sb" uniqKey="Holton S">SB Holton</name>
</author>
<author><name sortKey="Desmone, Ra" uniqKey="Desmone R">RA Desmone</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Huson, Dh" uniqKey="Huson D">DH Huson</name>
</author>
<author><name sortKey="Auch, Af" uniqKey="Auch A">AF Auch</name>
</author>
<author><name sortKey="Qi, J" uniqKey="Qi J">J Qi</name>
</author>
<author><name sortKey="Schuster, Sc" uniqKey="Schuster S">SC Schuster</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Fourment, M" uniqKey="Fourment M">M Fourment</name>
</author>
<author><name sortKey="Gibbs, Mj" uniqKey="Gibbs M">MJ Gibbs</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Fraser, C" uniqKey="Fraser C">C Fraser</name>
</author>
<author><name sortKey="Donnelly, Ca" uniqKey="Donnelly C">CA Donnelly</name>
</author>
<author><name sortKey="Cauchemez, S" uniqKey="Cauchemez S">S Cauchemez</name>
</author>
<author><name sortKey="Hanage, Wp" uniqKey="Hanage W">WP Hanage</name>
</author>
<author><name sortKey="Van Kerkhove, Md" uniqKey="Van Kerkhove M">MD Van Kerkhove</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hancock, K" uniqKey="Hancock K">K Hancock</name>
</author>
<author><name sortKey="Veguilla, V" uniqKey="Veguilla V">V Veguilla</name>
</author>
<author><name sortKey="Lu, X" uniqKey="Lu X">X Lu</name>
</author>
<author><name sortKey="Zhong, W" uniqKey="Zhong W">W Zhong</name>
</author>
<author><name sortKey="Butler, En" uniqKey="Butler E">EN Butler</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Dietz, K" uniqKey="Dietz K">K Dietz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Towers, S" uniqKey="Towers S">S Towers</name>
</author>
<author><name sortKey="Feng, Z" uniqKey="Feng Z">Z Feng</name>
</author>
</analytic>
</biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">PLoS One</journal-id>
<journal-id journal-id-type="iso-abbrev">PLoS ONE</journal-id>
<journal-id journal-id-type="publisher-id">plos</journal-id>
<journal-id journal-id-type="pmc">plosone</journal-id>
<journal-title-group><journal-title>PLoS ONE</journal-title>
</journal-title-group>
<issn pub-type="epub">1932-6203</issn>
<publisher><publisher-name>Public Library of Science</publisher-name>
<publisher-loc>San Francisco, USA</publisher-loc>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">23637746</article-id>
<article-id pub-id-type="pmc">3634039</article-id>
<article-id pub-id-type="publisher-id">PONE-D-12-33244</article-id>
<article-id pub-id-type="doi">10.1371/journal.pone.0060343</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject>
</subj-group>
<subj-group subj-group-type="Discipline-v2"><subject>Mathematics</subject>
<subj-group><subject>Applied Mathematics</subject>
</subj-group>
<subj-group><subject>Mathematical Computing</subject>
</subj-group>
</subj-group>
<subj-group subj-group-type="Discipline-v2"><subject>Medicine</subject>
<subj-group><subject>Epidemiology</subject>
</subj-group>
<subj-group><subject>Infectious Diseases</subject>
<subj-group><subject>Viral Diseases</subject>
<subj-group><subject>Influenza</subject>
</subj-group>
</subj-group>
<subj-group><subject>Infectious Disease Control</subject>
<subject>Infectious Disease Modeling</subject>
</subj-group>
</subj-group>
<subj-group><subject>Public Health</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group><article-title>A Perspective on Multiple Waves of Influenza Pandemics</article-title>
<alt-title alt-title-type="running-head">Multiple Waves of Influenza Pandemics</alt-title>
</title-group>
<contrib-group><contrib contrib-type="author"><name><surname>Mummert</surname>
<given-names>Anna</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Weiss</surname>
<given-names>Howard</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Long</surname>
<given-names>Li-Ping</given-names>
</name>
<xref ref-type="aff" rid="aff3"><sup>3</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Amigó</surname>
<given-names>José M.</given-names>
</name>
<xref ref-type="aff" rid="aff4"><sup>4</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Wan</surname>
<given-names>Xiu-Feng</given-names>
</name>
<xref ref-type="aff" rid="aff3"><sup>3</sup>
</xref>
<xref ref-type="corresp" rid="cor1"><sup>*</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1"><label>1</label>
<addr-line>Department of Mathematics, Marshall University, Huntington, West Virginia, United States of America</addr-line>
</aff>
<aff id="aff2"><label>2</label>
<addr-line>School of Mathematics, Georgia Institute of Technology, Atlanta, Georgia, United States of America</addr-line>
</aff>
<aff id="aff3"><label>3</label>
<addr-line>Department of Basic Sciences, Mississippi State University, Mississippi State, Mississippi, United States of America</addr-line>
</aff>
<aff id="aff4"><label>4</label>
<addr-line>Centro de Investigación Operativa, Universidad Miguel Hernández, Elche, Spain</addr-line>
</aff>
<contrib-group><contrib contrib-type="editor"><name><surname>Roberts</surname>
<given-names>Michael George</given-names>
</name>
<role>Editor</role>
<xref ref-type="aff" rid="edit1"></xref>
</contrib>
</contrib-group>
<aff id="edit1"><addr-line>Massey University, New Zealand</addr-line>
</aff>
<author-notes><corresp id="cor1">* E-mail: <email>wan@cvm.msstate.edu</email>
</corresp>
<fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>
The authors have declared that no competing interests exist.</p>
</fn>
<fn fn-type="con"><p>Conceived and designed the experiments: XFW HW. Performed the experiments: AM HW XFW LPL. Analyzed the data: AM HW XFW JA. Contributed reagents/materials/analysis tools: XFW. Wrote the paper: AM HW XFW.</p>
</fn>
</author-notes>
<pub-date pub-type="collection"><year>2013</year>
</pub-date>
<pub-date pub-type="epub"><day>23</day>
<month>4</month>
<year>2013</year>
</pub-date>
<volume>8</volume>
<issue>4</issue>
<elocation-id>e60343</elocation-id>
<history><date date-type="received"><day>21</day>
<month>10</month>
<year>2012</year>
</date>
<date date-type="accepted"><day>25</day>
<month>2</month>
<year>2013</year>
</date>
</history>
<permissions><copyright-statement>© 2013 Mummert et al</copyright-statement>
<copyright-year>2013</copyright-year>
<copyright-holder>Mummert et al</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><license-p>This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are properly credited.</license-p>
</license>
</permissions>
<abstract><sec><title>Background</title>
<p>A striking characteristic of the past four influenza pandemic outbreaks in the United States has been the multiple waves of infections. However, the mechanisms responsible for the multiple waves of influenza or other acute infectious diseases are uncertain. Understanding these mechanisms could provide knowledge for health authorities to develop and implement prevention and control strategies.</p>
</sec>
<sec><title>Materials and Methods</title>
<p>We exhibit five distinct mechanisms, each of which can generate two waves of infections for an acute infectious disease. The first two mechanisms capture changes in virus transmissibility and behavioral changes. The third mechanism involves population heterogeneity (e.g., demography, geography), where each wave spreads through one sub-population. The fourth mechanism is virus mutation which causes delayed susceptibility of individuals. The fifth mechanism is waning immunity. Each mechanism is incorporated into separate mathematical models, and outbreaks are then simulated. We use the models to examine the effects of the initial number of infected individuals (e.g., border control at the beginning of the outbreak) and the timing of and amount of available vaccinations.</p>
</sec>
<sec><title>Results</title>
<p>Four models, individually or in any combination, reproduce the two waves of the 2009 H1N1 pandemic in the United States, both qualitatively and quantitatively. One model reproduces the two waves only qualitatively. All models indicate that significantly reducing or delaying the initial numbers of infected individuals would have little impact on the attack rate. Instead, this reduction or delay results in a single wave as opposed to two waves. Furthermore, four of these models also indicate that a vaccination program started earlier than October 2009 (when the H1N1 vaccine was initially distributed) could have eliminated the second wave of infection, while more vaccine available starting in October would not have eliminated the second wave.</p>
</sec>
</abstract>
<funding-group><funding-statement>L.P. Long and X.-F. Wan were supported by 2010-DD-BX-0596 from the United States Department of Justice to X.-F. Wan. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</funding-statement>
</funding-group>
<counts><page-count count="12"></page-count>
</counts>
</article-meta>
</front>
<body><sec id="s1"><title>Introduction</title>
<p>An influenza pandemic occurs when a new strain of the influenza virus emerges, usually through antigenic shift, for which there is little or no immunity in the human population. The mutation is such that it is able to cause illness in a single individual and then spreads person-to-person worldwide. In the 20th century, the world experienced three influenza pandemics: the 1918 H1N1 "Spanish flu'', the 1957 H2N2 "Asian flu'', and the 1968 H3N2 "Hong Kong flu''. The first influenza pandemic of the 21st century occurred in 2009 and was caused by a swine-origin H1N1 influenza A virus <xref rid="pone.0060343-Miller1" ref-type="bibr">[1]</xref>
.</p>
<p>During each of these four pandemics, the United States experienced multiple waves of infections, where the numbers of infections and deaths exhibited well-separated temporal peaks with a separation time-scale of months <xref rid="pone.0060343-Miller1" ref-type="bibr">[1]</xref>
. For example, the first wave of the 2009 pandemic in the United States began in March and peaked in late June and early July. There were markedly fewer cases throughout August, and the second larger wave peaked in late October, early November (<xref ref-type="fig" rid="pone-0060343-g001">Figure 1</xref>
, Table S1 in <xref ref-type="supplementary-material" rid="pone.0060343.s001">File S1</xref>
). While many countries such as the United States experienced at least two waves of infections during the 2009 pandemic, other countries such as China experienced only a single wave of infection.</p>
<fig id="pone-0060343-g001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0060343.g001</object-id>
<label>Figure 1</label>
<caption><title>The H1N1 positive tests reported to the CDC in the United States from April 24, 2009 (week 17) to November 27, 2009 (week 48).</title>
</caption>
<graphic xlink:href="pone.0060343.g001"></graphic>
</fig>
<p>The underlying mechanisms leading to single or multiple waves of acute infectious diseases are not well understood. In this paper, we describe and explore several mechanisms that can produce multiple waves.</p>
<p>For several infectious diseases, including influenza, the timing of school vacations is thought to be a mechanism for multiple waves because children have reduced transmission due to far fewer contacts during vacations than when school is in session. Parts of either the summer or winter school breaks occurred during the gaps between the two waves of the past four influenza pandemics in the United States. Evidence supporting this potential mechanism for the second wave of the 2009 US H1N1 pandemic include the surge in hospital outpatient visits for influenza-like illness approximately two weeks after schools re-opened in the fall of 2009 <xref rid="pone.0060343-Chao1" ref-type="bibr">[2]</xref>
. However, school vacations clearly cannot be the sole cause of multiple waves <xref rid="pone.0060343-Chao1" ref-type="bibr">[2]</xref>
, <xref rid="pone.0060343-Chowell1" ref-type="bibr">[3]</xref>
, <xref rid="pone.0060343-Monto1" ref-type="bibr">[4]</xref>
, <xref rid="pone.0060343-Cauchemez1" ref-type="bibr">[5]</xref>
, <xref rid="pone.0060343-Cauchemez2" ref-type="bibr">[6]</xref>
. For instance, during the 1968 pandemic in the United States, the second wave began in November, two months after school resumed. In addition, children have school vacations every year, while seasonal influenzas have only single waves of infections.</p>
<p>Another possible mechanism for multiple waves is due to rapid mutations of the H1N1 virus. Many RNA viruses, such as influenza, quickly acquire genetic variants through random mutations, some of which lead to non-synonymous mutations in protein sequences. The more genetic variants that there are, the higher the odds are that one of them carries a mutation upon which selection can act. Through random mutation and subsequent selection, RNA viruses evolve into a form better adapted for human-to-human transmission.</p>
<p>The protein sequence analyses for the 2009 H1N1 HA genes showed the average number of mutations increased slightly from April (2.75±0.71, based on the comparison with the A/California/07/2009 strain) to November of 2009 (4.26±1.53). The maximum number of mutations occurred in November 2009 (11 residues when compared with A/California/07/2009, versus two residues when compared with the major strain circulating in April 2009). Virus diversity, which we introduce in the next section, quantifies the extent of the mutations.</p>
<p>Some of these mutations occurred in the receptor binding sites of HA genes <xref rid="pone.0060343-Melidou1" ref-type="bibr">[7]</xref>
, <xref rid="pone.0060343-Chan1" ref-type="bibr">[8]</xref>
and in other segments of the virus <xref rid="pone.0060343-Pan1" ref-type="bibr">[9]</xref>
. Some mutations first appeared after the first wave ended (e.g., HA-S220T NA-N248D in Japan) <xref rid="pone.0060343-Morlighem1" ref-type="bibr">[10]</xref>
. As the pandemic progressed, the number of mutations at the receptor binding site position 222 increased around the world <xref rid="pone.0060343-Melidou1" ref-type="bibr">[7]</xref>
, <xref rid="pone.0060343-Wang1" ref-type="bibr">[11]</xref>
, <xref rid="pone.0060343-Drews1" ref-type="bibr">[12]</xref>
, <xref rid="pone.0060343-Miller2" ref-type="bibr">[13]</xref>
, <xref rid="pone.0060343-Valli1" ref-type="bibr">[14]</xref>
. The mutation D222G/N was observed in 90% of the blood samples of the A/H1N1/2009 viremia cases <xref rid="pone.0060343-Chan1" ref-type="bibr">[8]</xref>
, <xref rid="pone.0060343-Tse1" ref-type="bibr">[15]</xref>
, <xref rid="pone.0060343-Puzelli1" ref-type="bibr">[16]</xref>
, <xref rid="pone.0060343-Mak1" ref-type="bibr">[17]</xref>
, <xref rid="pone.0060343-Kilander1" ref-type="bibr">[18]</xref>
. The D222G mutation has been shown to increase virulence and to increase the virus transmissibility in animal models, as well in human airway epithelial cell lines <xref rid="pone.0060343-Chutinimitkul1" ref-type="bibr">[19]</xref>
, <xref rid="pone.0060343-Belser1" ref-type="bibr">[20]</xref>
, <xref rid="pone.0060343-Abed1" ref-type="bibr">[21]</xref>
, <xref rid="pone.0060343-Liu1" ref-type="bibr">[22]</xref>
. Other mutations in the receptor binding sites of the HA gene <xref rid="pone.0060343-Wang2" ref-type="bibr">[23]</xref>
, <xref rid="pone.0060343-Zheng1" ref-type="bibr">[24]</xref>
, and other mutations at other segments (e.g., PA, PB1-F2, PB2, NP, and NS1) have been found in clinical isolates and have been shown to increase the replication efficiency and pathogenesis <italic>in vitro</italic>
in animal models. Increased virulence and pathogenesis of such mutations is a mechanism that may lead to a second wave of infection.</p>
<p>In this study, five mathematical models are formulated to explore complementary mechanisms which can produce two waves of acute infections. Four reproduce the two waves of the 2009 H1N1 influenza pandemic, quantitatively. The other reproduces the waves qualitatively. The first two mechanisms capture changing contact rates and changing (or evolving) virus transmissibility. In these models, the changes manifest as a time varying transmission rate. The transmission rate is the per capita rate of infection given contact, and depends in a highly complex way on both the contact rate between susceptible and infected individuals and the transmissibility of the infection. The first mechanism uses a periodic transmission rate to capture the seasonal contact rate. The second mechanism incorporates all sources of variability in contacts and transmissibility into the time-dependent transmission rate.</p>
<p>The third mechanism incorporates a heterogeneous population with individuals split into two weakly interacting sub-populations. The split could be based on demographics, geography, or variations in immunity, among others. The two populations only very weakly interact, that is, the transmission between the two groups is low. Two waves appear with each sub-population experiencing only one wave of infection.</p>
<p>The fourth mechanism is virus mutation which causes delayed susceptibility of some individuals. As the 2009 H1N1 pandemic progressed in the United States, the virus mutated rapidly and created new viable quasispecies. We hypothesize that some of these mutations increased the transmissibility of the influenza virus and that this increased viral transmissibility caused some individuals who were not able to be infected during the first wave to become infected during the second wave.</p>
<p>The fifth mechanism is waning immunity, where recovered individuals lose immunity to the influenza virus at which point they again become susceptible to infection. The second wave appears because some individuals become infected for a second time.</p>
<p>The five mechanisms are incorporated into appropriate extensions of the standard (Susceptible-Exposed-Infected-Removed) transmission model <xref rid="pone.0060343-Towers1" ref-type="bibr">[47]</xref>
. Models 1–4 reproduce the two waves during the 2009 pandemic in the United States, qualitatively and quantitatively. Model 5 reproduces the two waves qualitatively.</p>
<p>Countries such as China implemented border control strategies and experienced only a single wave of infections. We implement border control in each model to determine whether the US could also have experienced one wave with border control. In addition, the US vaccination production and distribution data is incorporated into the models to test the effect of increasing vaccine availability and altering the timing of the vaccination effort.</p>
</sec>
<sec id="s2"><title>Results</title>
<sec id="s2a"><title>Models 1 and 2: SEIR models with variable transmission rate</title>
<p>Models 1 and 2 explore the ability of a variable transmission rate to generate the multiple waves of infection in a homogenous population. To do this we use a time-dependent, or variable, transmission rate β<sub>t</sub>
in a classical SEIR compartment transmission model (see Supplementary Information in <xref ref-type="supplementary-material" rid="pone.0060343.s001">File S1</xref>
, equations (S1)–(S4)). We provide two strikingly different transmission rates, each of which can reproduce the two waves. Model 1 exploits seasonality. We demonstrate that the school year hypothesis of a higher transmission rate during the school year and a lower rate during the summer and school vacations can generate the two waves. Clearly, this type of seasonality alone provides no insight into why two waves occur only during pandemic years and not every year; there must be additional factors driving the pandemic waves. Therefore, we provide Model 2, also with a variable transmission rate that captures the intricate interplay between the contact rate and transmissibility of the infection.</p>
</sec>
<sec id="s2b"><title>Model 1: Periodic transmission rate</title>
<p>Many experts believe that contacts between school age children play a significant role in influenza transmission and in producing pandemic waves. The contact rate is thought to be significantly higher for school age children during school terms than during school vacations <xref rid="pone.0060343-Chao1" ref-type="bibr">[2]</xref>
, <xref rid="pone.0060343-Chowell1" ref-type="bibr">[3]</xref>
, <xref rid="pone.0060343-Monto1" ref-type="bibr">[4]</xref>
, <xref rid="pone.0060343-Cauchemez1" ref-type="bibr">[5]</xref>
, <xref rid="pone.0060343-Cauchemez2" ref-type="bibr">[6]</xref>
. The presumption is that the first wave begins in the winter or early spring when school is in session and wanes when the children have significantly reduced contact during their summer breaks. The second wave begins soon after school resumes in late summer or early fall, when the contact rate increases.</p>
<p>To capture the seasonality of school contacts, the transmission rate β<sub>t</sub>
is set to be the periodic function (with period one year)<disp-formula id="pone.0060343.e001"><graphic xlink:href="pone.0060343.e001.jpg" position="anchor" orientation="portrait"></graphic>
</disp-formula>
(<xref ref-type="fig" rid="pone-0060343-g002">Figure 2A</xref>
), where there is higher transmission when students are in school, and lower transmission over the summer months.</p>
<fig id="pone-0060343-g002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0060343.g002</object-id>
<label>Figure 2</label>
<caption><title>Reproducing multiple waves using Model 1 with periodic time-dependent transmission rate, β<sub>t</sub>
.</title>
<p><bold></bold>
(A) The periodic transmission rate, which is low during the summer break and higher when school is in session. Summer is approximately June and July (weeks 23–31 and days 161–217). (B) The model generated disease prevalence, compared with the CDC data scaled for underreporting.</p>
</caption>
<graphic xlink:href="pone.0060343.g002"></graphic>
</fig>
<p>The resulting model was simulated using parameters from Table S9 in <xref ref-type="supplementary-material" rid="pone.0060343.s001">File S1</xref>
. <xref ref-type="fig" rid="pone-0060343-g002">Figure 2B</xref>
shows the output of Model 1 superimposed on the derived CDC incidence data for the United States in 2009. The model clearly captures the two waves of infection.</p>
</sec>
<sec id="s2c"><title>Model 2: Derived time-dependent transmission rate</title>
<p>The transmission rate in this model captures all sources of variability of the contact rate between susceptible and infected individuals and the virus transmissibility. This variability could occur because of public health interventions, seasonality, or the evolution of the virus, among other factors.</p>
<p>The transmission rate β<sub>t</sub>
is determined using a new algorithm <xref rid="pone.0060343-Pollicott1" ref-type="bibr">[25]</xref>
, <xref rid="pone.0060343-Hadeler1" ref-type="bibr">[26]</xref>
, <xref rid="pone.0060343-Mummert1" ref-type="bibr">[27]</xref>
that ensures the output <italic>I</italic>
(<italic>t</italic>
) of the model <italic>perfectly</italic>
fits a smooth interpolation of the reported cases by the Centers for Disease Prevention and Control (CDC)<xref rid="pone.0060343-CDC1" ref-type="bibr">[28]</xref>
. There is no error between the model output of <italic>I</italic>
(<italic>t</italic>
) and the data, and only the initial value β<sub>0</sub>
needs to be specified.</p>
<p>The resulting model was simulated using the parameters in Table S10 in <xref ref-type="supplementary-material" rid="pone.0060343.s001">File S1</xref>
. <xref ref-type="fig" rid="pone-0060343-g003">Figure 3A</xref>
shows the output of Model 2 superimposed on the derived CDC incidence data for the United States in 2009. Recall that the near perfect agreement was ensured by our choice of β<sub>t</sub>
. Errors in the simulation are due to the accumulation of small numerical errors.</p>
<fig id="pone-0060343-g003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0060343.g003</object-id>
<label>Figure 3</label>
<caption><title>Reproducing multiple waves using Model 2 with extracted transmission rate, βt.</title>
<p>(A) The model generated disease prevalence generated for the entire range 0.54 ≤β<sub>0</sub>
≤1.03, compared with the CDC data scaled for underreporting. (B) The extracted transmission rates recovered for three β<sub>0</sub>
values. Top green curve β<sub>0</sub>
= 0.73; middle blue curve β<sub>0</sub>
= 0.63; bottom red curve β<sub>0</sub>
= 0.54.</p>
</caption>
<graphic xlink:href="pone.0060343.g003"></graphic>
</fig>
<p><xref ref-type="fig" rid="pone-0060343-g003">Figure 3A</xref>
was generated using β<sub>0</sub>
= 0.63, however any β<sub>0</sub>
in the range 0.54≤β<sub>0</sub>
≤1.03 will produce the same model output. Any β<sub>0</sub>
outside of this range produces unrealistic outputs. Since β<sub>t</sub>
depends on β<sub>0</sub>
in a complicated way, in <xref ref-type="fig" rid="pone-0060343-g003">Figure 3B</xref>
we show representative graphs of β<sub>t</sub>
across a range of values for β<sub>0</sub>
.</p>
<p>In conclusion, Model 2 reproduces the data, and thus the two waves of infection. The mechanism of Model 2 is the construction of a time dependent transmission rate β<sub>t</sub>
, capturing all of the fluxuations in transmission due to changing contact rates, virus transmissibility, environment, and so on. The resulting transmission function indicates that the second wave cannot be explained as simply a drop in contacts due to school closures. A drop in contact rate due to school closures would appear as a corresponding drop in β<sub>t</sub>
. However, no such drop is seen in the transmission function during the summer months.</p>
</sec>
<sec id="s2d"><title>Model 3: Two weakly interacting sub-populations</title>
<p>For this model, we hypothesize that the population is split into two groups with minimal interaction. Each group experiences one wave of infection resulting in two waves for the entire population.</p>
<p>In Model 3, the total population is split into two sub-populations; population 1 is 22% of the total and the rest are in population 2. The SEIR model (Supplementary Information equations (S1)–(S4) in <xref ref-type="supplementary-material" rid="pone.0060343.s001">File S1</xref>
) is extended to eight equations, one S, E, I, and R equation for each sub-population. The sub-populations weakly interact, meaning that the transmission rate between populations is small. This model has three constant transmission rates: β<sub>1</sub>
for the transmission among sub-population 1, β<sub>2</sub>
for the transmission among sub-population 2, and β<sub>3</sub>
for the transmission between the two sub-populations. This model assumes that only sub-population 1 is initially infected and sub-population 2 does not acquire the infection until some time during the summer at which point some infected individuals appear (perhaps arrive home from traveling) who are infected.</p>
<p>Model 3 was simulated using the parameters in Table S11 in <xref ref-type="supplementary-material" rid="pone.0060343.s001">File S1</xref>
. Model 3 reproduced the two waves qualitatively (<xref ref-type="fig" rid="pone-0060343-g004">Figure 4B</xref>
), with each sub-population experiencing only one wave (<xref ref-type="fig" rid="pone-0060343-g004">Figure 4A</xref>
). Its attack rate is just below that of the 2009 pandemic H1N1 (20.08% compared with the CDC estimate of 25%). These simulations verify that two weakly interacting populations can experience two waves.</p>
<fig id="pone-0060343-g004" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0060343.g004</object-id>
<label>Figure 4</label>
<caption><title>Reproducing multiple waves using Model 3 with two weakly interacting sub-populations.</title>
<p>(A) Infection wave of sub-population 1 (blue) and infection wave of sub-population 2 (black). (B) The model generated disease prevalence compared with the CDC data scaled for underreporting.</p>
</caption>
<graphic xlink:href="pone.0060343.g004"></graphic>
</fig>
</sec>
<sec id="s2e"><title>Model 4: Virus mutation and delayed susceptibility</title>
<p>During the course of the 2009 influenza pandemic, the genetic diversity of the virus in both the United States and around the world increased significantly. For this model, we hypothesize that the increase in genetic diversity results in strains with higher transmissibility in the human population. This hypothesis supports the notion that some individuals who were not susceptible to the virus during the first wave become susceptible when new strains with greater transmissibility emerge as the pandemic progresses.</p>
<p>We incorporate this idea into the standard SEIR model (with constant transmission rate <italic>β</italic>
) by including a reserved class <italic>N</italic>
consisting of individuals who are initially non-susceptible (see equations (S5)–(S9) in the Supplementary Information, <xref ref-type="supplementary-material" rid="pone.0060343.s001">File S1</xref>
). These individuals become susceptible as the viral genetic diversity increases at a rate that is proportional to the time-dependent genetic diversity <italic>d</italic>
(t).</p>
<p>Genetic diversity of virus quasispecies is defined as the mean of pairwise genetic distances among the genomes of the viruses isolated in one month. In bioinformatics there are various measures of genetic diversity. Here, genetic diversity was measured using three distances: the <italic>p-distance</italic>
, the <italic>patristic</italic>
, and the <italic>CCV distance</italic>
between nucleotide sequences. These distances were normalized to attain values between 0 and 1, and the genetic diversity function was generated assuming that the diversity was constant throughout the entire month. In the United States data both of these genetic diversities increased slowly during the first wave and dramatically during the second (<xref ref-type="fig" rid="pone-0060343-g005">Figure 5</xref>
, Tables S2, S3, S4, S5, S6, S7, and S8 in <xref ref-type="supplementary-material" rid="pone.0060343.s001">File S1</xref>
).</p>
<fig id="pone-0060343-g005" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0060343.g005</object-id>
<label>Figure 5</label>
<caption><title>Dynamics of genetic diversity of the HA gene of 2009 H1N1 pandemic virus in the United States.</title>
<p>Genetic diversity is defined as the mean of pairwise genetic distances among the genomes of the viruses isolated in one month. The distances were normalized to attain values between 0 and 1. Two genetic distances were measured: (A) p-distance, (B) patristic distance, (C) CCV distance. The genetic diversity of NA gene and other internal genes are shown in Tables S1, S2, S3, S4, S5, S6, and S7 in File S1.</p>
</caption>
<graphic xlink:href="pone.0060343.g005"></graphic>
</fig>
<p>Model 4 reproduced the two waves using all three distances (<xref ref-type="fig" rid="pone-0060343-g006">Figure 6</xref>
), using the parameters in Table S12 in <xref ref-type="supplementary-material" rid="pone.0060343.s001">File S1</xref>
. Here, we construct the genetic diversity function as a step function of order 0. However, Model 4 will also reproduce the second wave using linear, quadratic, etc., interpolations of the monthly distance data. Thus, the model is robust in terms of the distance measure used and the way that these distance measurements are combined to determine the diversity function <italic>d</italic>
(<italic>t</italic>
).</p>
<fig id="pone-0060343-g006" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0060343.g006</object-id>
<label>Figure 6</label>
<caption><title>Reproducing multiple waves using Model 4 with virus mutation causing delayed susceptibility of some individuals.</title>
<p>The model generated disease prevalence using the (A) p-distance, with <italic>γ</italic>
= 0.11, (B) patristic distance, with <italic>γ</italic>
= 0.09, (C) CCV distance, with <italic>γ</italic>
= 0.08, compared with the CDC data scaled for underreporting.</p>
</caption>
<graphic xlink:href="pone.0060343.g006"></graphic>
</fig>
<p>In conclusion, Model 4 qualitatively reproduces the US incidence data and quantitatively reproduces the attack rate for the 2009 pandemic H1N1. This model lends support to the supposition that virus mutations can drive the second wave through the addition of susceptible individuals who were not susceptible to previous virus quasispecies. The unknown factor of Model 4 is the initial fraction of reserve non-susceptible individuals.</p>
<p>Here the value <italic>N</italic>
(<italic>0</italic>
) = 0.35 was chosen so that the simulation recreated the data.</p>
</sec>
<sec id="s2f"><title>Model 5: Waning immunity</title>
<p>A fifth mechanism is known to produce two waves of infection – waning immunity <xref rid="pone.0060343-Chan1" ref-type="bibr">[8]</xref>
,<xref rid="pone.0060343-CDC3" ref-type="bibr">[36]</xref>
. We incorporate waning immunity to the H1N1 virus with an SEIR model where individuals in the removed class lose immunity at a constant rate and rejoin the susceptible class (see Supplementary Information equations (S10)–(S13), in <xref ref-type="supplementary-material" rid="pone.0060343.s001">File S1</xref>
).</p>
<p>This model was simulated using the parameters in Table S13 in <xref ref-type="supplementary-material" rid="pone.0060343.s001">File S1</xref>
. It qualitatively reproduces the two waves, but always (necessarily) with a smaller amplitude second wave (<xref ref-type="fig" rid="pone-0060343-g007">Figure 7</xref>
).</p>
<fig id="pone-0060343-g007" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0060343.g007</object-id>
<label>Figure 7</label>
<caption><title>Reproducing multiple waves using Model 5 with waning immunity.</title>
<p>The model generated disease prevalence, compared with the CDC data scaled for underreporting.</p>
</caption>
<graphic xlink:href="pone.0060343.g007"></graphic>
</fig>
</sec>
<sec id="s2g"><title>The effect of initial infections on the second wave</title>
<p>The models were used to investigate the consequences of a smaller number of initially infected individuals in the United States and to identify factors that may affect the on-set of waves.</p>
<p>We assumed that an early intervention, e.g., border control, would have reduced the initial number of imported cases from Mexico. We included this in our models by reducing the initial fraction of infected individuals, <italic>I</italic>
(0). For Model 3, we reduce the number of individuals of both sub-population 1 and 2. The simulation outputs (<xref ref-type="fig" rid="pone-0060343-g008">Figure 8</xref>
) illustrate the result of two different intervention levels. The first plot for each model shows the resulting number of infected individuals, assuming that interventions reduced the initial number of infected individuals by 25% of the reported value adjusted for underreporting, while the second plot corresponds to a simulation output that reduced the number of individuals initially infected by 90%.</p>
<fig id="pone-0060343-g008" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0060343.g008</object-id>
<label>Figure 8</label>
<caption><title>Effects of reduced and delayed initial infected individuals on influenza prevalence using (A) Model 1 with periodic transmission rate (B) Model 2 with extracted transmission rate, β<sub>t</sub>
(C) Model 3 with two weakly interacting sub-populations (D) Model 4 with virus mutation causing delayed susceptibility of some individuals (E) Model 5 with waning immunity.</title>
<p>The model generated disease prevalence with the control strength indicated, compared with the CDC data scaled for underreporting, where 25% strength of control means that the initial fraction of infected individuals is 75% of the adjusted reported fraction.</p>
</caption>
<graphic xlink:href="pone.0060343.g008"></graphic>
</fig>
<p>A value, called the attack rate (<italic>I<sub>total</sub>
</italic>
), was computed for all models. The attack rate is the total fraction of individuals who become infected during the disease outbreak. With no intervention, the attack rate is <italic>I<sub>total</sub>
</italic>
= 0.2457, 0.2258, 0.2008, and 0.2307 for Models 1, 2, 3, and 4, respectively. The attack rate barely changes with the addition of interventions, even strong interventions. For example, with a 90% effective intervention, the attack rate is <italic>I<sub>total</sub>
</italic>
= 0.2379, 0.2324, 0.1978, and 0.2132 for Models 1, 2, 3, and 4, respectively. In addition, in the simulations of Models 1, 2, and 4, interventions of over 80% result in a loss of the second wave. For Model 5, individuals may lose immunity to the disease and become infected a second time. Thus the notion of attack rate is not as well defined. We compute instead the total number of individuals who become infected once throughout the outbreak as 0.7395, and 0.7361 with 90% effective intervention.</p>
</sec>
<sec id="s2h"><title>The effect of vaccinations on the second wave</title>
<p>The models were used to investigate the consequences of a vaccination program started earlier or with an increased availability of vaccine in the United States on the outbreak progression. We assumed that vaccination would have reduced the number of susceptible (and reserve non-susceptible) individuals by the number of available vaccines, with a vaccine efficacy of 90% <xref rid="pone.0060343-Greenberg1" ref-type="bibr">[29]</xref>
. We have accounted for individuals receiving the vaccine according to the dates and amounts of distributed vaccine, along with a two week delay <xref rid="pone.0060343-CDC2" ref-type="bibr">[30]</xref>
. During the delay the immune system is generating antibodies and is still vulnerable to infection. We investigated both the scenarios of earlier available vaccine (0–5 months earlier) and the scenarios of more available vaccine (1–5 times as much available). <xref ref-type="fig" rid="pone-0060343-g009">Figure 9</xref>
illustrates the result of different vaccination efforts for the models.</p>
<fig id="pone-0060343-g009" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0060343.g009</object-id>
<label>Figure 9</label>
<caption><title>Effects of timing and amount of available vaccine on influenza prevalence using (A) Model 1 with periodic transmission rate (B) Model 2 with extracted transmission rate, β<sub>t</sub>
(C) Model 3 with two weakly interacting sub-populations (D) Model 4 with virus mutation causing delayed susceptibility of some individuals (E) Model 5 with waning immunity.</title>
<p>The model generated disease prevalence with the vaccination strategy indicated, compared with the CDC data scaled for underreporting. Here original timing begins vaccination on October 14, 2009, and one month earlier begins vaccination thirty days prior; original amount available is the daily amount of shipped vaccine with 90% efficacy, and twice as much available is twice the shipped amount with 90% efficacy.</p>
</caption>
<graphic xlink:href="pone.0060343.g009"></graphic>
</fig>
<p>Earlier vaccinations can both eliminate the second wave of infections and significantly reduce the attack rate. Simulations indicate that vaccinations must be distributed approximately 2 months prior to the actual distribution dates to eliminate the second wave in Models 1, 2, 4, and 5. A larger availability of vaccine can only reduce the attack rate when combined with earlier vaccination. For example, with the original amount of vaccine distributed in October 2009, the attack rate is <italic>I<sub>total</sub>
</italic>
= 0.2142, 0.2072, 0.1931, and 0.2143 for Models 1, 2, 3, and 4, respectively, while with five times the original amount distributed at the same time, the attack rate is <italic>I<sub>total</sub>
</italic>
= 0.1788, 0.1845, 0.1689, and 0.1931 for Models 1, 2, 3, and 4, respectively. For Model 5, the total number who are infected once with the original vaccine distribution is 0.6951, while it is 0.6670 with five times the original amount distributed with the original timing.</p>
<p><xref ref-type="table" rid="pone-0060343-t001">Table 1</xref>
provides a summary of the five models, including the number of peaks of infection resulting from border control and vaccination, the attack rate with border control, and a reference for the model, if available.</p>
<table-wrap id="pone-0060343-t001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0060343.t001</object-id>
<label>Table 1</label>
<caption><title>Summary of border control and vaccination on multiple waves of 2009 H1N1 influenza pandemics.</title>
</caption>
<alternatives><graphic id="pone-0060343-t001-1" xlink:href="pone.0060343.t001"></graphic>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead><tr><td align="left" rowspan="1" colspan="1">Model</td>
<td align="left" rowspan="1" colspan="1">Description</td>
<td align="left" rowspan="1" colspan="1">Result of Border Control (90%)</td>
<td align="left" rowspan="1" colspan="1">Attack Rate with Border Control (CDC est 25% with no Border Control)</td>
<td align="left" rowspan="1" colspan="1">Result of Vaccinating Two Months Earlier</td>
<td align="left" rowspan="1" colspan="1">Reference</td>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1">1</td>
<td align="left" rowspan="1" colspan="1">Periodic transmission rate</td>
<td align="left" rowspan="1" colspan="1">One wave</td>
<td align="left" rowspan="1" colspan="1">24%</td>
<td align="left" rowspan="1" colspan="1">One wave</td>
<td align="left" rowspan="1" colspan="1"><xref rid="pone.0060343-Dietz1" ref-type="bibr">[46]</xref>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">2</td>
<td align="left" rowspan="1" colspan="1">Extracted transmission rate</td>
<td align="left" rowspan="1" colspan="1">One wave</td>
<td align="left" rowspan="1" colspan="1">23%</td>
<td align="left" rowspan="1" colspan="1">One wave</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">3</td>
<td align="left" rowspan="1" colspan="1">Two weakly interacting populations</td>
<td align="left" rowspan="1" colspan="1">One wave (BC 99.9%)</td>
<td align="left" rowspan="1" colspan="1">20%</td>
<td align="left" rowspan="1" colspan="1">Two waves</td>
<td align="left" rowspan="1" colspan="1"><xref rid="pone.0060343-Towers1" ref-type="bibr">[47]</xref>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">4</td>
<td align="left" rowspan="1" colspan="1">Genetic diversity of flu quasispecies</td>
<td align="left" rowspan="1" colspan="1">One wave</td>
<td align="left" rowspan="1" colspan="1">23%</td>
<td align="left" rowspan="1" colspan="1">One wave</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"><italic>5</italic>
</td>
<td align="left" rowspan="1" colspan="1">Waning immunity</td>
<td align="left" rowspan="1" colspan="1">Two waves</td>
<td align="left" rowspan="1" colspan="1">73%</td>
<td align="left" rowspan="1" colspan="1">One wave</td>
<td align="left" rowspan="1" colspan="1"><xref rid="pone.0060343-Chan1" ref-type="bibr">[8]</xref>
</td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
</sec>
</sec>
<sec id="s3"><title>Discussion</title>
<p>We study five mechanisms that can produce two waves of infection, and may explain the two waves experienced in the United States during the 2009 pandemic H1N1. These mechanisms are implemented and explored in mathematical models. Four models perform well, recreating the two waves (Figures 2, 3, 5 and 7) qualitatively and quantitatively, while one recreates the two waves only qualitatively. The two waves can be reproduced using the models individually or in any combination.</p>
<p>Two novel mechanisms (Models 2 and 4) are proposed. Models 1 and 2 incorporate a variable transmission rate, which was assumed constant in Models 3, 4 and 5. Model 4 explicitly incorporates the genetic diversity of the influenza quasispecies. We deem these models (2 and 4) complementary since changes in the transmission rate and transmission patterns among sub-populations could have been due, in part, to the emergence of more transmissible H1N1 mutants, which increases the genetic diversity of the quasispecies.</p>
<p>Model 1 is an SEIR component model of transmission with a time-dependent transmission rate reflecting the contact rate of school age children. It produces two waves with realistic parameter values, but provides no insight into why two waves of infection appear only with pandemics and not every year.</p>
<p>Model 2 is an SEIR model of transmission with a time-dependent transmission rate. Such a time-dependent transmission rate can account for temporal changes in virulence, seasonality, and contact rates, among others. The transmission rate is obtained from an algorithm that guarantees that the model output perfectly agrees with the infection data, even if the assumptions of the model do not apply. Thus one must be careful applying this algorithm and not over fit the data. We do not use our models to predict numbers of infections; we use them to exhibit possible mechanisms for two waves. In addition, we believe that the assumptions of a time-dependent transmission rate along with the structure of an SEIR model are reasonable for influenza, thus the potential weakness of over fitting is avoided.</p>
<p>The extracted transmission rate for Model 2 seems to be at odds with the common belief that summer affords fewer interactions amongst children <xref rid="pone.0060343-Chao1" ref-type="bibr">[2]</xref>
, <xref rid="pone.0060343-Chowell1" ref-type="bibr">[3]</xref>
, <xref rid="pone.0060343-Monto1" ref-type="bibr">[4]</xref>
, <xref rid="pone.0060343-Cauchemez1" ref-type="bibr">[5]</xref>
, <xref rid="pone.0060343-Cauchemez2" ref-type="bibr">[6]</xref>
and therefore causes a drop in the transmission rate. For all graphs in <xref ref-type="fig" rid="pone-0060343-g003">Figure 3B</xref>
, there was no significant drop in β<sub>t</sub>
after the school year ended; the transmission rate fluctuated around a nearly constant level in the first wave from mid-May through the end of August. It increased when the fall school term and the second wave began. Although at first surprising, there is no contradiction with a decrease in the number of infected individuals and a (basically) constant transmission rate, as is seen over the summer months. A constant transmission rate implies that the ratio of new infections to the number of susceptible individuals remains constant, which happens over the summer when both the numbers of new infections and of susceptible individuals are decreasing.</p>
<p>From days 150 to 250, the transmission rate β<sub>t</sub>
oscillated with mean period 19.9. These oscillations appear more rapid than the contact rate and environmental factors can account for. We do not see how this oscillation can be a manifestation of reporting conventions. However, these oscillations seem essential to reproduce the reported data. Figure S1A in <xref ref-type="supplementary-material" rid="pone.0060343.s001">File S1</xref>
shows a 21-day Gaussian filtered (smoothed) β<sub>t</sub>
function, and the corresponding model simulation (Figure S1B in <xref ref-type="supplementary-material" rid="pone.0060343.s001">File S1</xref>
), which is a substantially lower quality fit to the data. All attempts at smoothing the oscillations using different types of filters resulted in a worse fit to the data. Presumably, smoothing the data would eliminate these oscillations, but this goes in a very different direction.</p>
<p>Model 3 consists of two weakly coupled SEIR models of transmission, with a constant transmission rate. We have found that this model is the least robust, since it is the most sensitive to changes in the model parameter values. By construction, the model will always exhibit two waves regardless of the initial number of infected individuals (e.g., with a strong border control).</p>
<p>Model 4 involves virus mutation causing delayed susceptibility of some individuals. This new mechanism posits that as the genetic structure of the influenza virus changed during the pandemic, more of the general population became susceptible to infection. Such a dynamic of transmissibility of influenza A virus has been hypothesized to generate the second wave <xref rid="pone.0060343-Jackson1" ref-type="bibr">[31]</xref>
. Model 4 assumes that the initially non-susceptible individuals become susceptible at a rate proportional to the genetic diversity of the flu quasispecies. The effects of vaccination and border control are considered for Model 4 (<xref ref-type="fig" rid="pone-0060343-g008">Figure 8D</xref>
, <xref ref-type="fig" rid="pone-0060343-g009">Figure 9D</xref>
) using the CCV distance, and similar results hold for both the p-distance and the patristic distance.</p>
<p>Model 5 incorporates waning immunity to the influenza virus causing some individuals to become infected again. This model qualitatively reproduces the two waves of the 2009 H1N1 pandemic in the US; however, the first wave will always be larger than the second. In addition, waning immunity requires an unrealistically short period of immunity (six months) to accurately reproduce the timing of the two waves of the 2009 H1N1 pandemic. For these reasons, we believe that this mechanism played at most a minor role in the second wave.</p>
<p>While all combinations of all five models will generate the two waves, it is impossible to determine the actual mechanism(s) from any model. This is a fundamental limitation of mathematical modeling. Models can suggest plausible mechanisms, but they cannot prove them. There is so little data from other countries and other pandemics that any type of statistical investigation trying to narrow down the true explanation is impossible.</p>
<p>Although age structure probably plays a role in flu transmission, our models do not explicitly include such structure and they do not need to. We are able to illustrate that all five mechanisms can reproduce two waves without explicit age structure in the model. We point out that the two populations in Model 3 could be two age groups, though other groups are also possible.</p>
<sec id="s3a"><title>The transmission rate during the summer school vacation</title>
<p>For all graphs in <xref ref-type="fig" rid="pone-0060343-g003">Figure 3B</xref>
, there was no significant drop in β<sub>t</sub>
after the school year ended; the transmission rate fluxuated around a nearly constant level in the first wave from mid-May through the end of August. It increased when the fall school term and the second wave began. Although at first surprising, there is no contradiction with a decrease in the number of infected individuals and a (basically) constant transmission rate, as is seen over the summer months. A constant transmission rate implies that the ratio of new infections to the number of susceptible individuals remains constant, which happens over the summer when both the numbers of new infections and of susceptible individuals are decreasing.</p>
</sec>
<sec id="s3b"><title>The effects of border control</title>
<p>During the 2009 H1N1 pandemic, countries including China, Japan, Hong Kong, Singapore, and Australia implemented various border control strategies to prevent arriving airline passengers from infecting the countries' citizens <xref rid="pone.0060343-Wood1" ref-type="bibr">[32]</xref>
, <xref rid="pone.0060343-Bishop1" ref-type="bibr">[33]</xref>
, <xref rid="pone.0060343-Cowling1" ref-type="bibr">[34]</xref>
, <xref rid="pone.0060343-Miller3" ref-type="bibr">[35]</xref>
. Health authorities screened arriving passengers for flu-like symptoms using thermal scanners. In China, passengers suspected of being infected, along with passengers seated nearby, were placed into quarantine for seven days. China's border control began when the first cases of H1N1 were confirmed in California; China continued screening passengers for more than two months.</p>
<p>There has been debate as to whether the United States should have implemented some form of border control with Mexico after the first novel influenza cases appeared. The results from simulations indicate that strong intervention at the beginning of the outbreak would have had negligible impact on the attack rate in the United States, but it could have resulted in a single wave of infections. Our simulations imply that by significantly decreasing the initial number of infected individuals, the two waves would collapse into a single wave of infections, and the peak number of infections would occur slightly earlier with one wave. This suggests that China's strong border control could be a potential mechanism explaining their one wave of infection.</p>
<p>In practice, the efficacy of such border control measures to reduce transmission is uncertain. These screenings will miss asymptomatic individuals <xref rid="pone.0060343-Miller3" ref-type="bibr">[35]</xref>
, and in particular infected individuals with no fever. However border control has been shown to delay the local community transmission <xref rid="pone.0060343-Wood1" ref-type="bibr">[32]</xref>
, <xref rid="pone.0060343-Bishop1" ref-type="bibr">[33]</xref>
, <xref rid="pone.0060343-Cowling1" ref-type="bibr">[34]</xref>
, which may allow sufficient time for health interventions, such as vaccinations. On the other hand, border control may only delay transmission and could still result in a second wave, as seen in the simulations.</p>
</sec>
<sec id="s3c"><title>The effects of vaccination</title>
<p>The first doses of 2009 H1N1 vaccine were available in the United States in early October 2009, which was just before the peak of the 2009 H1N1 pandemic. Conventional wisdom is that this vaccination had little effect on the course of the disease spread. Though we were unable to find any empirical evidence in the literature, our simulations support this idea. Simulations of all models indicate that the actual vaccine schedule did not significantly reduce the total number of infections and even a significantly larger amount of vaccine available in October would have little effect on the attack rate. However, an earlier vaccination program would significantly decrease the attack rate and can eliminate the second wave of infections. Simulations show that the attack rate would be significantly reduced had the vaccine been distributed two months earlier than it was (with the same availability).</p>
<p>Recall that by design, Model 3 will produce two waves. Extremely strict border control (99.9% control) can prevent the second wave; however starting the vaccination program five months earlier with the original amount of vaccine cannot prevent the second wave.</p>
</sec>
</sec>
<sec sec-type="materials|methods" id="s4"><title>Materials and Methods</title>
<sec id="s4a"><title>Data sets</title>
<p>The epidemiologic dataset of 2009 weekly H1N1 positive tests reported to the CDC was downloaded from the CDC online weekly reports (<ext-link ext-link-type="uri" xlink:href="http://www.cdc.gov/flu/">www.cdc.gov/flu/</ext-link>
). Clearly, these weekly numbers are highly affected by the testing efforts. An estimated 25% of the U.S. population was infected with the pandemic H1N1 throughout the two waves <xref rid="pone.0060343-CDC3" ref-type="bibr">[36]</xref>
. Using this along with the CDC compiled data, the influenza confirmed case data was adjusted to reflect underreporting by a scale factor of 700∶1. The authors are aware that a team of investigators estimated an underreporting scale factor of 79∶1 at the beginning (April through July) of the pandemic <xref rid="pone.0060343-Wood1" ref-type="bibr">[32]</xref>
. This is not a contradiction, since, very quickly, many physicians and hospitals stopped testing and the underreporting factor skyrocketed. Since we could find so little information about the underreporting rate during the entire pandemic in the U.S., we chose a uniform scaling factor.</p>
<p>The genomic sequence dataset was downloaded from the Influenza Virus Resource Database <xref rid="pone.0060343-Bao1" ref-type="bibr">[37]</xref>
, which was updated on November 16, 2011. This dataset includes 1,502 HA genes of 2009 pH1N1 pandemic influenza viruses (229 were isolated in April of 2009; 273 in May; 247 in June; 74 in July; 46 in August; 119 in September; 195 in October; 214 in November; and 105 in December). All HA genes were fully or nearly fully sequenced, and only a single HA sequence for the same strain was selected.</p>
</sec>
<sec id="s4b"><title>Measurement of genetic distance and genetic diversity</title>
<p>The genetic diversity of the influenza virus quasispecies is defined as the mean of the pairwise genetic distances among the genes of the virus. We employ three notions of genetic distance: the p-distance, the patristic distance, and the CCV distance <xref rid="pone.0060343-Long1" ref-type="bibr">[38]</xref>
, <xref rid="pone.0060343-Wan1" ref-type="bibr">[39]</xref>
. P-distance measures the proportion of nucleotide sites between two sequences. This is the simplest measure and does not take into consideration multiple substitutions at the same site, substitution rate biases, or differences in the evolutionary rates among sites <xref rid="pone.0060343-Nei1" ref-type="bibr">[40]</xref>
. The p-distances were measured using the software package MEGAN <xref rid="pone.0060343-Huson1" ref-type="bibr">[41]</xref>
. The patristic distance encodes the total amount of genetic change that exists between genetic sequences by summing the lengths of branches in the phylogenetic tree between the two viruses. The patristic distances were calculated using PATRISTIC <xref rid="pone.0060343-Fourment1" ref-type="bibr">[42]</xref>
using the maximum likelihood phylogenetic trees constructed by GARLI <xref rid="pone.0060343-Zwickl1" ref-type="bibr">[43]</xref>
. CCV is a method for measuring genetic distance based on the probability of strings [Wan et al. 2007], and this method has been shown to be effective in measuring the genetic distance between influenza genes [Wan et al. 2007b; Wan et al. 2008a; Wan et al. 2008b; Long et al. 2012].</p>
</sec>
<sec id="s4c"><title>Model parameter selection</title>
<p>Every model parameter value is consistent with the range commonly found in the literature. Within the common ranges, parameters for Models 1–4 were found that quantitatively and qualitatively reproduce the two waves of the 2009 H1N1 pandemic in the US. For Model 5, no parameters were found from the common ranges that reproduce the two waves quantitatively.</p>
<p>The incubation period for influenza can range between one and four days. We follow modeling convention and use a one day incubation period in all models. Adults shed influenza virus for up to 11 days; however, the amount shed is significantly lower after the first 3 to 5 days. Several previous modeling studies of pandemic H1N1 have assumed the infectivity period (1/<italic>v</italic>
) is the lower bound of 3 days, as is done in all models here.</p>
<p>The transmission rate <italic>β</italic>
(or the initial transmission rate for Models 1 and 2) is chosen to ensure that the basic reproduction number R<sub>0</sub>
is in the allowable range for influenza. The basic reproduction number is the average number of secondary cases caused by one infected individual into an entirely susceptible population. However, it is well known that approximately 20% of the US population had prior immunity to the pandemic influenza. Therefore, it is unreasonable to compute R<sub>0</sub>
assuming that the entire population was susceptible. We use instead the standard definition where it is not assumed that all individuals are initially susceptible. For the 2009 H1N1 influenza, R<sub>0</sub>
is estimated to be in the range 1.2–2.25 <xref rid="pone.0060343-Fraser1" ref-type="bibr">[44]</xref>
. For Model 1, the transmission rate at the start of the outbreak (t = 140) is 0.5267, resulting in an R<sub>0</sub>
of 1.58. For Model 2, simulations restrict β<sub>0</sub>
to the range 0.54 to 1.03. <xref ref-type="fig" rid="pone-0060343-g003">Figure 3</xref>
shows β<sub>0</sub>
in the range 0.54 to 0.74, which gives a basic reproduction number in the range 1.39 to 1.88. For Model 3, at the beginning of the infection β<sub>1</sub>
= 2.1, which gives a basic reproduction number of R<sub>0</sub>
= 1.386, and for the second wave, β<sub>2</sub>
= 0.565, which gives a basic reproduction number of R<sub>0</sub>
= 1.322. For Model 4, β = 0.9, and so the basic reproduction number is R<sub>0</sub>
= 1.20. For Model 5, β = 0.63, and so the basic reproduction number is R<sub>0</sub>
= 1.89.</p>
<p>For Models 1, 3, 4, and 5, the initial fraction of infected individuals, <italic>I</italic>
(<italic>0</italic>
) (I<sub>1</sub>
(0) for Model 3) is set by the US incidence data and the initial fraction of exposed individuals <italic>E</italic>
(<italic>0</italic>
) is set to be 0. For Models 1, 4 and 5, the initial fraction of removed individuals <italic>R</italic>
(<italic>0</italic>
) is set to 0.2, corresponding to elderly individuals with acquired immunity <xref rid="pone.0060343-Hancock1" ref-type="bibr">[45]</xref>
; in Model 3, this value is set to 0. For Model 3, the infection in sub-population 2 begins at day t = 210 and I<sub>2</sub>
(210) is the US incidence data scaled by 0.06.</p>
<p>For Model 2, the initial values <italic>S</italic>
(<italic>0</italic>
), <italic>E</italic>
(<italic>0</italic>
), <italic>I</italic>
(<italic>0</italic>
), <italic>and R</italic>
(<italic>0</italic>
) are determined by the algorithm which computes the time-dependent transmission rate β<sub>t</sub>
. Therefore the choice of β<sub>0</sub>
determines all four of the initial values. For β<sub>0</sub>
= 0.63 the initial conditions agree reasonably well with the initial conditions of the other models.</p>
<p>For Model 4, N(0) was chosen to ensure that the basic reproduction rate R<sub>0</sub>
was within the known bounds. With β<sub>0</sub>
= 0.9, the value of N(0) can range from 0 to 0.35. In addition, model simulations show that to generate two waves, the value of N(0) must be larger than 0.25. With the chosen value of N(0) = 0.35 the value of S(0) is forced to be 0.44.</p>
</sec>
</sec>
<sec sec-type="supplementary-material" id="s5"><title>Supporting Information</title>
<supplementary-material content-type="local-data" id="pone.0060343.s001"><label>File S1</label>
<caption><p>Supplementary Equations, Tables, and Figures.</p>
<p>(DOC)</p>
</caption>
<media xlink:href="pone.0060343.s001.doc"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
</sec>
</body>
<back><ack><p>We are grateful to Marc Lipsitch and Rustom Antia for critical comments on a preliminary draft of this manuscript, and to Jialiang Yang for his assistance in the genetic distance measurements.</p>
</ack>
<ref-list><title>References</title>
<ref id="pone.0060343-Miller1"><label>1</label>
<mixed-citation publication-type="journal"><name><surname>Miller</surname>
<given-names>MA</given-names>
</name>
, <name><surname>Viboud</surname>
<given-names>C</given-names>
</name>
, <name><surname>Balinska</surname>
<given-names>M</given-names>
</name>
, <name><surname>Simonsen</surname>
<given-names>L</given-names>
</name>
(<year>2009</year>
) <article-title>The signature features of influenza pandemics-implications for policy</article-title>
. <source>N Engl J Med</source>
<volume>360</volume>
: <fpage>2595</fpage>
–<lpage>2598</lpage>
.<pub-id pub-id-type="pmid">19423872</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Chao1"><label>2</label>
<mixed-citation publication-type="journal"><name><surname>Chao</surname>
<given-names>DL</given-names>
</name>
, <name><surname>Halloran</surname>
<given-names>ME</given-names>
</name>
, <name><surname>Longini</surname>
<given-names>IM</given-names>
<suffix>Jr</suffix>
</name>
(<year>2010</year>
) <article-title>School opening dates predict pandemic influenza A(H1N1) outbreaks in the United States</article-title>
. <source>J Infect Dis</source>
<volume>202</volume>
: <fpage>877</fpage>
–<lpage>880</lpage>
.<pub-id pub-id-type="pmid">20704486</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Chowell1"><label>3</label>
<mixed-citation publication-type="journal"><name><surname>Chowell</surname>
<given-names>G</given-names>
</name>
, <name><surname>Echevarria-Zuno</surname>
<given-names>S</given-names>
</name>
, <name><surname>Viboud</surname>
<given-names>C</given-names>
</name>
, <name><surname>Simonsen</surname>
<given-names>L</given-names>
</name>
, <name><surname>Tamerius</surname>
<given-names>J</given-names>
</name>
, <etal>et al</etal>
(<year>2011</year>
) <article-title>Characterizing the epidemiology of the 2009 influenza A/H1N1 pandemic in Mexico</article-title>
. <source>PLoS Med</source>
<volume>8</volume>
: <fpage>e1000436</fpage>
.<pub-id pub-id-type="pmid">21629683</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Monto1"><label>4</label>
<mixed-citation publication-type="journal"><name><surname>Monto</surname>
<given-names>AS</given-names>
</name>
, <name><surname>Koopman</surname>
<given-names>JS</given-names>
</name>
, <name><surname>Longini</surname>
<given-names>IM</given-names>
<suffix>Jr</suffix>
</name>
(<year>1985</year>
) <article-title>Tecumseh study of illness. XIII. Influenza infection and disease, 1976–1981</article-title>
. <source>Am J Epidemiol</source>
<volume>121</volume>
: <fpage>811</fpage>
–<lpage>822</lpage>
.<pub-id pub-id-type="pmid">4014174</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Cauchemez1"><label>5</label>
<mixed-citation publication-type="journal"><name><surname>Cauchemez</surname>
<given-names>S</given-names>
</name>
, <name><surname>Ferguson</surname>
<given-names>NM</given-names>
</name>
, <name><surname>Wachtel</surname>
<given-names>C</given-names>
</name>
, <name><surname>Tegnell</surname>
<given-names>A</given-names>
</name>
, <name><surname>Saour</surname>
<given-names>G</given-names>
</name>
, <etal>et al</etal>
(<year>2009</year>
) <article-title>Closure of schools during an influenza pandemic</article-title>
. <source>Lancet Infect Dis</source>
<volume>9</volume>
: <fpage>473</fpage>
–<lpage>481</lpage>
.<pub-id pub-id-type="pmid">19628172</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Cauchemez2"><label>6</label>
<mixed-citation publication-type="journal"><name><surname>Cauchemez</surname>
<given-names>S</given-names>
</name>
, <name><surname>Valleron</surname>
<given-names>AJ</given-names>
</name>
, <name><surname>Boelle</surname>
<given-names>PY</given-names>
</name>
, <name><surname>Flahault</surname>
<given-names>A</given-names>
</name>
, <name><surname>Ferguson</surname>
<given-names>NM</given-names>
</name>
(<year>2008</year>
) <article-title>Estimating the impact of school closure on influenza transmission from Sentinel data</article-title>
. <source>Nature</source>
<volume>452</volume>
: <fpage>750</fpage>
–<lpage>754</lpage>
.<pub-id pub-id-type="pmid">18401408</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Melidou1"><label>7</label>
<mixed-citation publication-type="journal"><name><surname>Melidou</surname>
<given-names>A</given-names>
</name>
, <name><surname>Gioula</surname>
<given-names>G</given-names>
</name>
, <name><surname>Exindari</surname>
<given-names>M</given-names>
</name>
, <name><surname>Chatzidimitriou</surname>
<given-names>D</given-names>
</name>
, <name><surname>Diza</surname>
<given-names>E</given-names>
</name>
, <etal>et al</etal>
(<year>2010</year>
) <article-title>Molecular and phylogenetic analysis of the haemagglutinin gene of pandemic influenza H1N1 2009 viruses associated with severe and fatal infections</article-title>
. <source>Virus Res</source>
<volume>151</volume>
: <fpage>192</fpage>
–<lpage>199</lpage>
.<pub-id pub-id-type="pmid">20493216</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Chan1"><label>8</label>
<mixed-citation publication-type="journal"><name><surname>Chan</surname>
<given-names>PK</given-names>
</name>
, <name><surname>Lee</surname>
<given-names>N</given-names>
</name>
, <name><surname>Joynt</surname>
<given-names>GM</given-names>
</name>
, <name><surname>Choi</surname>
<given-names>KW</given-names>
</name>
, <name><surname>Cheung</surname>
<given-names>JL</given-names>
</name>
, <etal>et al</etal>
(<year>2011</year>
) <article-title>Clinical and virological course of infection with haemagglutinin D222G mutant strain of 2009 pandemic influenza A (H1N1) virus</article-title>
. <source>J Clin Virol</source>
<volume>50</volume>
: <fpage>320</fpage>
–<lpage>324</lpage>
.<pub-id pub-id-type="pmid">21330192</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Pan1"><label>9</label>
<mixed-citation publication-type="journal"><name><surname>Pan</surname>
<given-names>C</given-names>
</name>
, <name><surname>Cheung</surname>
<given-names>B</given-names>
</name>
, <name><surname>Tan</surname>
<given-names>S</given-names>
</name>
, <name><surname>Li</surname>
<given-names>C</given-names>
</name>
, <name><surname>Li</surname>
<given-names>L</given-names>
</name>
, <etal>et al</etal>
(<year>2010</year>
) <article-title>Genomic signature and mutation trend analysis of pandemic (H1N1) 2009 influenza A virus</article-title>
. <source>PLoS One</source>
<volume>5</volume>
: <fpage>e9549</fpage>
.<pub-id pub-id-type="pmid">20221396</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Morlighem1"><label>10</label>
<mixed-citation publication-type="journal"><name><surname>Morlighem</surname>
<given-names>JE</given-names>
</name>
, <name><surname>Aoki</surname>
<given-names>S</given-names>
</name>
, <name><surname>Kishima</surname>
<given-names>M</given-names>
</name>
, <name><surname>Hanami</surname>
<given-names>M</given-names>
</name>
, <name><surname>Ogawa</surname>
<given-names>C</given-names>
</name>
, <etal>et al</etal>
(<year>2011</year>
) <article-title>Mutation analysis of 2009 pandemic influenza A(H1N1) viruses collected in Japan during the peak phase of the pandemic</article-title>
. <source>PLoS One</source>
<volume>6</volume>
: <fpage>e18956</fpage>
.<pub-id pub-id-type="pmid">21572517</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Wang1"><label>11</label>
<mixed-citation publication-type="journal"><name><surname>Wang</surname>
<given-names>B</given-names>
</name>
, <name><surname>Dwyer</surname>
<given-names>DE</given-names>
</name>
, <name><surname>Soedjono</surname>
<given-names>M</given-names>
</name>
, <name><surname>Shi</surname>
<given-names>H</given-names>
</name>
, <name><surname>Matlho</surname>
<given-names>K</given-names>
</name>
, <etal>et al</etal>
(<year>2011</year>
) <article-title>Evidence of the circulation of pandemic influenza (H1N1) 2009 with D222D/G/N/S hemagglutinin polymorphisms during the first wave of the 2009 influenza pandemic</article-title>
. <source>J Clin Virol</source>
<volume>52</volume>
: <fpage>304</fpage>
–<lpage>306</lpage>
.<pub-id pub-id-type="pmid">21925936</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Drews1"><label>12</label>
<mixed-citation publication-type="journal"><name><surname>Drews</surname>
<given-names>SJ</given-names>
</name>
, <name><surname>Pabbaraju</surname>
<given-names>K</given-names>
</name>
, <name><surname>Wong</surname>
<given-names>S</given-names>
</name>
, <name><surname>Tokaryk</surname>
<given-names>KL</given-names>
</name>
, <name><surname>May-Hadford</surname>
<given-names>J</given-names>
</name>
, <etal>et al</etal>
(<year>2011</year>
) <article-title>Surveillance of autopsy cases for D222G substitutions in haemagglutinin of the pandemic (H1N1) 2009 virus in Alberta, Canada</article-title>
. <source>Clin Microbiol Infect</source>
<volume>17</volume>
: <fpage>582</fpage>
–<lpage>584</lpage>
.<pub-id pub-id-type="pmid">20718801</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Miller2"><label>13</label>
<mixed-citation publication-type="journal"><name><surname>Miller</surname>
<given-names>RR</given-names>
</name>
, <name><surname>MacLean</surname>
<given-names>AR</given-names>
</name>
, <name><surname>Gunson</surname>
<given-names>RN</given-names>
</name>
, <name><surname>Carman</surname>
<given-names>WF</given-names>
</name>
(<year>2010</year>
) <article-title>Occurrence of haemagglutinin mutation D222G in pandemic influenza A(H1N1) infected patients in the West of Scotland, United Kingdom, 2009–10</article-title>
. <source>Euro Surveill</source>
<volume>15</volume>
.</mixed-citation>
</ref>
<ref id="pone.0060343-Valli1"><label>14</label>
<mixed-citation publication-type="journal"><name><surname>Valli</surname>
<given-names>MB</given-names>
</name>
, <name><surname>Selleri</surname>
<given-names>M</given-names>
</name>
, <name><surname>Meschi</surname>
<given-names>S</given-names>
</name>
, <name><surname>Zaccaro</surname>
<given-names>P</given-names>
</name>
, <name><surname>Vincenti</surname>
<given-names>D</given-names>
</name>
, <etal>et al</etal>
(<year>2011</year>
) <article-title>Hemagglutinin 222 variants in pandemic (H1N1) 2009 virus</article-title>
. <source>Emerg Infect Dis</source>
<volume>17</volume>
: <fpage>749</fpage>
–<lpage>751</lpage>
.<pub-id pub-id-type="pmid">21470481</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Tse1"><label>15</label>
<mixed-citation publication-type="journal"><name><surname>Tse</surname>
<given-names>H</given-names>
</name>
, <name><surname>To</surname>
<given-names>KK</given-names>
</name>
, <name><surname>Wen</surname>
<given-names>X</given-names>
</name>
, <name><surname>Chen</surname>
<given-names>H</given-names>
</name>
, <name><surname>Chan</surname>
<given-names>KH</given-names>
</name>
, <etal>et al</etal>
(<year>2011</year>
) <article-title>Clinical and virological factors associated with viremia in pandemic influenza A/H1N1/2009 virus infection</article-title>
. <source>PLoS One</source>
<volume>6</volume>
: <fpage>e22534</fpage>
.<pub-id pub-id-type="pmid">21980333</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Puzelli1"><label>16</label>
<mixed-citation publication-type="journal"><name><surname>Puzelli</surname>
<given-names>S</given-names>
</name>
, <name><surname>Facchini</surname>
<given-names>M</given-names>
</name>
, <name><surname>Spagnolo</surname>
<given-names>D</given-names>
</name>
, <name><surname>De Marco</surname>
<given-names>MA</given-names>
</name>
, <name><surname>Calzoletti</surname>
<given-names>L</given-names>
</name>
, <etal>et al</etal>
(<year>2010</year>
) <article-title>Transmission of hemagglutinin D222G mutant strain of pandemic (H1N1) 2009 virus</article-title>
. <source>Emerg Infect Dis</source>
<volume>16</volume>
: <fpage>863</fpage>
–<lpage>865</lpage>
.<pub-id pub-id-type="pmid">20409386</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Mak1"><label>17</label>
<mixed-citation publication-type="journal"><name><surname>Mak</surname>
<given-names>GC</given-names>
</name>
, <name><surname>Au</surname>
<given-names>KW</given-names>
</name>
, <name><surname>Tai</surname>
<given-names>LS</given-names>
</name>
, <name><surname>Chuang</surname>
<given-names>KC</given-names>
</name>
, <name><surname>Cheng</surname>
<given-names>KC</given-names>
</name>
, <etal>et al</etal>
(<year>2010</year>
) <article-title>Association of D222G substitution in haemagglutinin of 2009 pandemic influenza A (H1N1) with severe disease</article-title>
. <source>Euro Surveill</source>
<volume>15</volume>
.</mixed-citation>
</ref>
<ref id="pone.0060343-Kilander1"><label>18</label>
<mixed-citation publication-type="journal"><name><surname>Kilander</surname>
<given-names>A</given-names>
</name>
, <name><surname>Rykkvin</surname>
<given-names>R</given-names>
</name>
, <name><surname>Dudman</surname>
<given-names>SG</given-names>
</name>
, <name><surname>Hungnes</surname>
<given-names>O</given-names>
</name>
(<year>2010</year>
) <article-title>Observed association between the HA1 mutation D222G in the 2009 pandemic influenza A(H1N1) virus and severe clinical outcome, Norway 2009–2010</article-title>
. <source>Euro Surveill</source>
<volume>15</volume>
.</mixed-citation>
</ref>
<ref id="pone.0060343-Chutinimitkul1"><label>19</label>
<mixed-citation publication-type="other">Chutinimitkul S, Herfst S, Steel J, Lowen AC, Ye J, <etal>et al</etal>
.. (2010) Virulence-associated substitution D222G in hemagglutinin of 2009 pandemic influenza A(H1N1) virus affects receptor binding. J Virol.</mixed-citation>
</ref>
<ref id="pone.0060343-Belser1"><label>20</label>
<mixed-citation publication-type="journal"><name><surname>Belser</surname>
<given-names>JA</given-names>
</name>
, <name><surname>Jayaraman</surname>
<given-names>A</given-names>
</name>
, <name><surname>Raman</surname>
<given-names>R</given-names>
</name>
, <name><surname>Pappas</surname>
<given-names>C</given-names>
</name>
, <name><surname>Zeng</surname>
<given-names>H</given-names>
</name>
, <etal>et al</etal>
(<year>2011</year>
) <article-title>Effect of D222G mutation in the hemagglutinin protein on receptor binding, pathogenesis and transmissibility of the 2009 pandemic H1N1 influenza virus</article-title>
. <source>PLoS One</source>
<volume>6</volume>
: <fpage>e25091</fpage>
.<pub-id pub-id-type="pmid">21966421</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Abed1"><label>21</label>
<mixed-citation publication-type="journal"><name><surname>Abed</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Pizzorno</surname>
<given-names>A</given-names>
</name>
, <name><surname>Hamelin</surname>
<given-names>ME</given-names>
</name>
, <name><surname>Leung</surname>
<given-names>A</given-names>
</name>
, <name><surname>Joubert</surname>
<given-names>P</given-names>
</name>
, <etal>et al</etal>
(<year>2011</year>
) <article-title>The 2009 pandemic H1N1 D222G hemagglutinin mutation alters receptor specificity and increases virulence in mice but not in ferrets</article-title>
. <source>J Infect Dis</source>
<volume>204</volume>
: <fpage>1008</fpage>
–<lpage>1016</lpage>
.<pub-id pub-id-type="pmid">21881115</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Liu1"><label>22</label>
<mixed-citation publication-type="journal"><name><surname>Liu</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Childs</surname>
<given-names>RA</given-names>
</name>
, <name><surname>Matrosovich</surname>
<given-names>T</given-names>
</name>
, <name><surname>Wharton</surname>
<given-names>S</given-names>
</name>
, <name><surname>Palma</surname>
<given-names>AS</given-names>
</name>
, <etal>et al</etal>
(<year>2010</year>
) <article-title>Altered receptor specificity and cell tropism of D222G hemagglutinin mutants isolated from fatal cases of pandemic A(H1N1) 2009 influenza virus</article-title>
. <source>J Virol</source>
<volume>84</volume>
: <fpage>12069</fpage>
–<lpage>12074</lpage>
.<pub-id pub-id-type="pmid">20826688</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Wang2"><label>23</label>
<mixed-citation publication-type="journal"><name><surname>Wang</surname>
<given-names>W</given-names>
</name>
, <name><surname>Castelan-Vega</surname>
<given-names>JA</given-names>
</name>
, <name><surname>Jimenez-Alberto</surname>
<given-names>A</given-names>
</name>
, <name><surname>Vassell</surname>
<given-names>R</given-names>
</name>
, <name><surname>Ye</surname>
<given-names>Z</given-names>
</name>
, <etal>et al</etal>
(<year>2010</year>
) <article-title>A mutation in the receptor binding site enhances infectivity of 2009 H1N1 influenza hemagglutinin pseudotypes without changing antigenicity</article-title>
. <source>Virology</source>
<volume>407</volume>
: <fpage>374</fpage>
–<lpage>380</lpage>
.<pub-id pub-id-type="pmid">20869738</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Zheng1"><label>24</label>
<mixed-citation publication-type="journal"><name><surname>Zheng</surname>
<given-names>B</given-names>
</name>
, <name><surname>Chan</surname>
<given-names>KH</given-names>
</name>
, <name><surname>Zhang</surname>
<given-names>AJ</given-names>
</name>
, <name><surname>Zhou</surname>
<given-names>J</given-names>
</name>
, <name><surname>Chan</surname>
<given-names>CC</given-names>
</name>
, <etal>et al</etal>
(<year>2010</year>
) <article-title>D225G mutation in hemagglutinin of pandemic influenza H1N1 (2009) virus enhances virulence in mice</article-title>
. <source>Exp Biol Med (Maywood)</source>
<volume>235</volume>
: <fpage>981</fpage>
–<lpage>988</lpage>
.<pub-id pub-id-type="pmid">20660098</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Pollicott1"><label>25</label>
<mixed-citation publication-type="other">Pollicott M, Wang H, Weiss H (2009) Recovering the Time-dependent Transmission Rate from Infection Data. Available from: arXiv:09073529.</mixed-citation>
</ref>
<ref id="pone.0060343-Hadeler1"><label>26</label>
<mixed-citation publication-type="journal"><name><surname>Hadeler</surname>
<given-names>KP</given-names>
</name>
(<year>2011</year>
) <article-title>Parameter identification in epidemic models</article-title>
. <source>Math Biosci</source>
<volume>229</volume>
: <fpage>185</fpage>
–<lpage>189</lpage>
.<pub-id pub-id-type="pmid">21192953</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Mummert1"><label>27</label>
<mixed-citation publication-type="journal"><name><surname>Mummert</surname>
<given-names>A</given-names>
</name>
(<year>2012</year>
) <article-title>Studying the Recovery Procedure for the Time-dependent Transmission Rate(s) in Epidemic Models</article-title>
. <source>Journal of Mathmatical Biology</source>
<source>Available from DOI 10.1007/s00285-012-0558-1</source>
</mixed-citation>
</ref>
<ref id="pone.0060343-CDC1"><label>28</label>
<mixed-citation publication-type="other">CDC (2010) 2009 H1N1 Flu. CDC.</mixed-citation>
</ref>
<ref id="pone.0060343-Greenberg1"><label>29</label>
<mixed-citation publication-type="journal"><name><surname>Greenberg</surname>
<given-names>ME</given-names>
</name>
, <name><surname>Lai</surname>
<given-names>MH</given-names>
</name>
, <name><surname>Hartel</surname>
<given-names>GF</given-names>
</name>
, <name><surname>Wichems</surname>
<given-names>CH</given-names>
</name>
, <name><surname>Gittleson</surname>
<given-names>C</given-names>
</name>
, <etal>et al</etal>
(<year>2009</year>
) <article-title>Response to a monovalent 2009 influenza A (H1N1) vaccine</article-title>
. <source>N Engl J Med</source>
<volume>361</volume>
: <fpage>2405</fpage>
–<lpage>2413</lpage>
.<pub-id pub-id-type="pmid">19745216</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-CDC2"><label>30</label>
<mixed-citation publication-type="other">CDC (2010) 2009 H1N1 Vaccine Doses Allocated, Ordered, and Shipped by Project Area CDC.</mixed-citation>
</ref>
<ref id="pone.0060343-Jackson1"><label>31</label>
<mixed-citation publication-type="journal"><name><surname>Jackson</surname>
<given-names>C</given-names>
</name>
, <name><surname>Vynnycky</surname>
<given-names>E</given-names>
</name>
, <name><surname>Mangtani</surname>
<given-names>P</given-names>
</name>
(<year>2010</year>
) <article-title>Estimates of the transmissibility of the 1968 (Hong Kong) influenza pandemic: evidence of increased transmissibility between successive waves</article-title>
. <source>Am J Epidemiol</source>
<volume>171</volume>
: <fpage>465</fpage>
–<lpage>478</lpage>
.<pub-id pub-id-type="pmid">20007674</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Wood1"><label>32</label>
<mixed-citation publication-type="journal"><name><surname>Wood</surname>
<given-names>JG</given-names>
</name>
, <name><surname>Zamani</surname>
<given-names>N</given-names>
</name>
, <name><surname>MacIntyre</surname>
<given-names>CR</given-names>
</name>
, <name><surname>Beckert</surname>
<given-names>NG</given-names>
</name>
(<year>2007</year>
) <article-title>Effects of internal border control on spread of pandemic influenza</article-title>
. <source>Emerg Infect Dis</source>
<volume>13</volume>
: <fpage>1038</fpage>
–<lpage>1045</lpage>
.<pub-id pub-id-type="pmid">18214176</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Bishop1"><label>33</label>
<mixed-citation publication-type="journal"><name><surname>Bishop</surname>
<given-names>JF</given-names>
</name>
, <name><surname>Murnane</surname>
<given-names>MP</given-names>
</name>
, <name><surname>Owen</surname>
<given-names>R</given-names>
</name>
(<year>2009</year>
) <article-title>Australia's winter with the 2009 pandemic influenza A (H1N1) virus</article-title>
. <source>N Engl J Med</source>
<volume>361</volume>
: <fpage>2591</fpage>
–<lpage>2594</lpage>
.<pub-id pub-id-type="pmid">19940287</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Cowling1"><label>34</label>
<mixed-citation publication-type="journal"><name><surname>Cowling</surname>
<given-names>BJ</given-names>
</name>
, <name><surname>Lau</surname>
<given-names>LL</given-names>
</name>
, <name><surname>Wu</surname>
<given-names>P</given-names>
</name>
, <name><surname>Wong</surname>
<given-names>HW</given-names>
</name>
, <name><surname>Fang</surname>
<given-names>VJ</given-names>
</name>
, <etal>et al</etal>
(<year>2010</year>
) <article-title>Entry screening to delay local transmission of 2009 pandemic influenza A (H1N1)</article-title>
. <source>BMC Infect Dis</source>
<volume>10</volume>
: <fpage>82</fpage>
.<pub-id pub-id-type="pmid">20353566</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Miller3"><label>35</label>
<mixed-citation publication-type="journal"><name><surname>Miller</surname>
<given-names>E</given-names>
</name>
, <name><surname>Hoschler</surname>
<given-names>K</given-names>
</name>
, <name><surname>Hardelid</surname>
<given-names>P</given-names>
</name>
, <name><surname>Stanford</surname>
<given-names>E</given-names>
</name>
, <name><surname>Andrews</surname>
<given-names>N</given-names>
</name>
, <etal>et al</etal>
(<year>2010</year>
) <article-title>Incidence of 2009 pandemic influenza A H1N1 infection in England: a cross-sectional serological study</article-title>
. <source>Lancet</source>
<volume>375</volume>
: <fpage>1100</fpage>
–<lpage>1108</lpage>
.<pub-id pub-id-type="pmid">20096450</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-CDC3"><label>36</label>
<mixed-citation publication-type="other">CDC (2010) Updated CDC Estimates of 2009 H1N1 Influenza Cases, Hospitalizations and Deaths in the United States, April 2009 – April 10, 2010. CDC.</mixed-citation>
</ref>
<ref id="pone.0060343-Bao1"><label>37</label>
<mixed-citation publication-type="journal"><name><surname>Bao</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Bolotov</surname>
<given-names>P</given-names>
</name>
, <name><surname>Dernovoy</surname>
<given-names>D</given-names>
</name>
, <name><surname>Kiryutin</surname>
<given-names>B</given-names>
</name>
, <name><surname>Zaslavsky</surname>
<given-names>L</given-names>
</name>
, <etal>et al</etal>
(<year>2008</year>
) <article-title>The influenza virus resource at the National Center for Biotechnology Information</article-title>
. <source>J Virol</source>
<volume>82</volume>
: <fpage>596</fpage>
–<lpage>601</lpage>
.<pub-id pub-id-type="pmid">17942553</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Long1"><label>38</label>
<mixed-citation publication-type="journal"><name><surname>Long</surname>
<given-names>LP</given-names>
</name>
, <name><surname>Yuan</surname>
<given-names>C</given-names>
</name>
, <name><surname>Cai</surname>
<given-names>Z</given-names>
</name>
, <name><surname>Xu</surname>
<given-names>H</given-names>
</name>
, <name><surname>Wan</surname>
<given-names>X-F</given-names>
</name>
(<year>2012</year>
) <article-title>Mixture model analysis reflecting dynamics of the population diversity of 2009 pandemic H1N1 influenza virus</article-title>
. <source>In Silico Biology</source>
<volume>11</volume>
: <fpage>225</fpage>
–<lpage>236</lpage>
.</mixed-citation>
</ref>
<ref id="pone.0060343-Wan1"><label>39</label>
<mixed-citation publication-type="journal"><name><surname>Wan</surname>
<given-names>XF</given-names>
</name>
, <name><surname>Wu</surname>
<given-names>X</given-names>
</name>
, <name><surname>Lin</surname>
<given-names>G</given-names>
</name>
, <name><surname>Holton</surname>
<given-names>SB</given-names>
</name>
, <name><surname>Desmone</surname>
<given-names>RA</given-names>
</name>
, <etal>et al</etal>
(<year>2007</year>
) <article-title>Computational identification of reassortments in avian influenza viruses</article-title>
. <source>Avian Dis</source>
<volume>51</volume>
: <fpage>434</fpage>
–<lpage>439</lpage>
.<pub-id pub-id-type="pmid">17494602</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Nei1"><label>40</label>
<mixed-citation publication-type="book">Nei M, Kumar S (2000) Molecular Evolution and Phylogenetics. New York: Oxford University Press.</mixed-citation>
</ref>
<ref id="pone.0060343-Huson1"><label>41</label>
<mixed-citation publication-type="journal"><name><surname>Huson</surname>
<given-names>DH</given-names>
</name>
, <name><surname>Auch</surname>
<given-names>AF</given-names>
</name>
, <name><surname>Qi</surname>
<given-names>J</given-names>
</name>
, <name><surname>Schuster</surname>
<given-names>SC</given-names>
</name>
(<year>2007</year>
) <article-title>MEGAN analysis of metagenomic data</article-title>
. <source>Genome Res</source>
<volume>17</volume>
: <fpage>377</fpage>
–<lpage>386</lpage>
.<pub-id pub-id-type="pmid">17255551</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Fourment1"><label>42</label>
<mixed-citation publication-type="journal"><name><surname>Fourment</surname>
<given-names>M</given-names>
</name>
, <name><surname>Gibbs</surname>
<given-names>MJ</given-names>
</name>
(<year>2006</year>
) <article-title>PATRISTIC: a program for calculating patristic distances and graphically comparing the components of genetic change</article-title>
. <source>BMC Evol Biol</source>
<volume>6</volume>
: <fpage>1</fpage>
.<pub-id pub-id-type="pmid">16388682</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Zwickl1"><label>43</label>
<mixed-citation publication-type="other">Zwickl DJ (2006) Genetic algorithm approaches for the phylogenetic analysis of large biological sequence datasets under the maximum likelihood criterion: The University of Texas at Austin.</mixed-citation>
</ref>
<ref id="pone.0060343-Fraser1"><label>44</label>
<mixed-citation publication-type="journal"><name><surname>Fraser</surname>
<given-names>C</given-names>
</name>
, <name><surname>Donnelly</surname>
<given-names>CA</given-names>
</name>
, <name><surname>Cauchemez</surname>
<given-names>S</given-names>
</name>
, <name><surname>Hanage</surname>
<given-names>WP</given-names>
</name>
, <name><surname>Van Kerkhove</surname>
<given-names>MD</given-names>
</name>
, <etal>et al</etal>
(<year>2009</year>
) <article-title>Pandemic potential of a strain of influenza A (H1N1): early findings</article-title>
. <source>Science</source>
<volume>324</volume>
: <fpage>1557</fpage>
–<lpage>1561</lpage>
.<pub-id pub-id-type="pmid">19433588</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Hancock1"><label>45</label>
<mixed-citation publication-type="journal"><name><surname>Hancock</surname>
<given-names>K</given-names>
</name>
, <name><surname>Veguilla</surname>
<given-names>V</given-names>
</name>
, <name><surname>Lu</surname>
<given-names>X</given-names>
</name>
, <name><surname>Zhong</surname>
<given-names>W</given-names>
</name>
, <name><surname>Butler</surname>
<given-names>EN</given-names>
</name>
, <etal>et al</etal>
(<year>2009</year>
) <article-title>Cross-reactive antibody responses to the 2009 pandemic H1N1 influenza virus</article-title>
. <source>N Engl J Med</source>
<volume>361</volume>
: <fpage>1945</fpage>
–<lpage>1952</lpage>
.<pub-id pub-id-type="pmid">19745214</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0060343-Dietz1"><label>46</label>
<mixed-citation publication-type="journal"><name><surname>Dietz</surname>
<given-names>K</given-names>
</name>
(<year>1976</year>
) <article-title>The incidence of infectious diseases under the influence of seasonal fluctuation</article-title>
. <source>Lect Notes Biomath</source>
<volume>11</volume>
: <fpage>1</fpage>
.</mixed-citation>
</ref>
<ref id="pone.0060343-Towers1"><label>47</label>
<mixed-citation publication-type="journal"><name><surname>Towers</surname>
<given-names>S</given-names>
</name>
, <name><surname>Feng</surname>
<given-names>Z</given-names>
</name>
(<year>2009</year>
) <article-title>Pandemic H1N1 influenza: predicting the course of a pandemic and assessing the efficacy of the planned vaccination programme in the United States</article-title>
. <source>Euro Surveill</source>
<volume>14</volume>
: <fpage>19358</fpage>
.<pub-id pub-id-type="pmid">19883540</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 000B609 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000B609 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Sante |area= H2N2V1 |flux= Pmc |étape= Corpus |type= RBID |clé= |texte= }}
![]() | This area was generated with Dilib version V0.6.33. | ![]() |