Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Age-specific contacts and travel patterns in the spatial spread of 2009 H1N1 influenza pandemic</title>
<author><name sortKey="Apolloni, Andrea" sort="Apolloni, Andrea" uniqKey="Apolloni A" first="Andrea" last="Apolloni">Andrea Apolloni</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="GRID">grid.8991.9</institution-id>
<institution-id institution-id-type="ISNI">000000040425469X</institution-id>
<institution>Department of Infectious Disease Epidemiology,</institution>
<institution>London School of Hygiene and Tropical Medicine,</institution>
</institution-wrap>
London, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Poletto, Chiara" sort="Poletto, Chiara" uniqKey="Poletto C" first="Chiara" last="Poletto">Chiara Poletto</name>
<affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="GRID">grid.418750.f</institution-id>
<institution-id institution-id-type="ISNI">0000000417593658</institution-id>
<institution>Computational Epidemiology Laboratory,</institution>
<institution>Institute for Scientific Interchange (ISI),</institution>
</institution-wrap>
Torino, Italy</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff3"><institution-wrap><institution-id institution-id-type="GRID">grid.7429.8</institution-id>
<institution-id institution-id-type="ISNI">0000000121866389</institution-id>
<institution>INSERM, U707,</institution>
</institution-wrap>
Paris, France</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff4"><institution-wrap><institution-id institution-id-type="GRID">grid.5805.8</institution-id>
<institution-id institution-id-type="ISNI">0000000119553500</institution-id>
<institution>Faculté de Médecine Pierre et Marie Curie,</institution>
<institution>UPMC Université Paris 06, UMR S 707,</institution>
</institution-wrap>
Paris, France</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Colizza, Vittoria" sort="Colizza, Vittoria" uniqKey="Colizza V" first="Vittoria" last="Colizza">Vittoria Colizza</name>
<affiliation><nlm:aff id="Aff3"><institution-wrap><institution-id institution-id-type="GRID">grid.7429.8</institution-id>
<institution-id institution-id-type="ISNI">0000000121866389</institution-id>
<institution>INSERM, U707,</institution>
</institution-wrap>
Paris, France</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff4"><institution-wrap><institution-id institution-id-type="GRID">grid.5805.8</institution-id>
<institution-id institution-id-type="ISNI">0000000119553500</institution-id>
<institution>Faculté de Médecine Pierre et Marie Curie,</institution>
<institution>UPMC Université Paris 06, UMR S 707,</institution>
</institution-wrap>
Paris, France</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff5"><institution-wrap><institution-id institution-id-type="GRID">grid.418750.f</institution-id>
<institution-id institution-id-type="ISNI">0000000417593658</institution-id>
<institution>Institute for Scientific Interchange (ISI),</institution>
</institution-wrap>
Torino, Italy</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">23587010</idno>
<idno type="pmc">3644502</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3644502</idno>
<idno type="RBID">PMC:3644502</idno>
<idno type="doi">10.1186/1471-2334-13-176</idno>
<date when="2013">2013</date>
<idno type="wicri:Area/Pmc/Corpus">000182</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000182</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Age-specific contacts and travel patterns in the spatial spread of 2009 H1N1 influenza pandemic</title>
<author><name sortKey="Apolloni, Andrea" sort="Apolloni, Andrea" uniqKey="Apolloni A" first="Andrea" last="Apolloni">Andrea Apolloni</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="GRID">grid.8991.9</institution-id>
<institution-id institution-id-type="ISNI">000000040425469X</institution-id>
<institution>Department of Infectious Disease Epidemiology,</institution>
<institution>London School of Hygiene and Tropical Medicine,</institution>
</institution-wrap>
London, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Poletto, Chiara" sort="Poletto, Chiara" uniqKey="Poletto C" first="Chiara" last="Poletto">Chiara Poletto</name>
<affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="GRID">grid.418750.f</institution-id>
<institution-id institution-id-type="ISNI">0000000417593658</institution-id>
<institution>Computational Epidemiology Laboratory,</institution>
<institution>Institute for Scientific Interchange (ISI),</institution>
</institution-wrap>
Torino, Italy</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff3"><institution-wrap><institution-id institution-id-type="GRID">grid.7429.8</institution-id>
<institution-id institution-id-type="ISNI">0000000121866389</institution-id>
<institution>INSERM, U707,</institution>
</institution-wrap>
Paris, France</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff4"><institution-wrap><institution-id institution-id-type="GRID">grid.5805.8</institution-id>
<institution-id institution-id-type="ISNI">0000000119553500</institution-id>
<institution>Faculté de Médecine Pierre et Marie Curie,</institution>
<institution>UPMC Université Paris 06, UMR S 707,</institution>
</institution-wrap>
Paris, France</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Colizza, Vittoria" sort="Colizza, Vittoria" uniqKey="Colizza V" first="Vittoria" last="Colizza">Vittoria Colizza</name>
<affiliation><nlm:aff id="Aff3"><institution-wrap><institution-id institution-id-type="GRID">grid.7429.8</institution-id>
<institution-id institution-id-type="ISNI">0000000121866389</institution-id>
<institution>INSERM, U707,</institution>
</institution-wrap>
Paris, France</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff4"><institution-wrap><institution-id institution-id-type="GRID">grid.5805.8</institution-id>
<institution-id institution-id-type="ISNI">0000000119553500</institution-id>
<institution>Faculté de Médecine Pierre et Marie Curie,</institution>
<institution>UPMC Université Paris 06, UMR S 707,</institution>
</institution-wrap>
Paris, France</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff5"><institution-wrap><institution-id institution-id-type="GRID">grid.418750.f</institution-id>
<institution-id institution-id-type="ISNI">0000000417593658</institution-id>
<institution>Institute for Scientific Interchange (ISI),</institution>
</institution-wrap>
Torino, Italy</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">BMC Infectious Diseases</title>
<idno type="eISSN">1471-2334</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>Confirmed H1N1 cases during late spring and summer 2009 in various countries showed a substantial age shift between importations and local transmission cases, with adults mainly responsible for seeding unaffected regions and children most frequently driving community outbreaks.</p>
</sec>
<sec><title>Methods</title>
<p>We introduce a multi-host stochastic metapopulation model with two age classes to analytically investigate the role of a heterogeneously mixing population and its associated non-homogeneous travel behaviors on the risk of a major epidemic. We inform the model with demographic data, contact data and travel statistics of Europe and Mexico, and calibrate it to the 2009 H1N1 pandemic early outbreak. We allow for variations of the model parameters to explore the conditions of invasion under different scenarios.</p>
</sec>
<sec><title>Results</title>
<p>We derive the expression for the potential of global invasion of the epidemic that depends on the transmissibility of the pathogen, the transportation network and mobility features, the demographic profile and the mixing pattern. Higher assortativity in the contact pattern greatly increases the probability of spatial containment of the epidemic, this effect being contrasted by an increase in the social activity of adults vs. children. Heterogeneous features of the mobility network characterizing its topology and traffic flows strongly favor the invasion of the pathogen at the spatial level, as also a larger fraction of children traveling. Variations in the demographic profile and mixing habits across countries lead to heterogeneous outbreak situations. Model results are compatible with the H1N1 spatial transmission dynamics observed.</p>
</sec>
<sec><title>Conclusions</title>
<p>This work illustrates the importance of considering age-dependent mixing profiles and mobility features coupled together to study the conditions for the spatial invasion of an emerging influenza pandemic. Its results allow the immediate assessment of the risk of a major epidemic for a specific scenario upon availability of data, and the evaluation of the potential effectiveness of public health interventions targeting specific age groups, their interactions and mobility behaviors. The approach provides a general modeling framework that can be used for other types of partitions of the host population and applied to different settings.</p>
</sec>
<sec><title>Electronic supplementary material</title>
<p>The online version of this article (doi:10.1186/1471-2334-13-176) contains supplementary material, which is available to authorized users.</p>
</sec>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Khan, K" uniqKey="Khan K">K Khan</name>
</author>
</analytic>
</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>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Balcan, D" uniqKey="Balcan D">D Balcan</name>
</author>
<author><name sortKey="Hu, H" uniqKey="Hu H">H Hu</name>
</author>
<author><name sortKey="Goncalves, B" uniqKey="Goncalves B">B Goncalves</name>
</author>
<author><name sortKey="Bajardi, P" uniqKey="Bajardi P">P Bajardi</name>
</author>
<author><name sortKey="Poletto, C" uniqKey="Poletto C">C Poletto</name>
</author>
<author><name sortKey="Ramasco, Jj" uniqKey="Ramasco J">JJ Ramasco</name>
</author>
<author><name sortKey="Paolotti, D" uniqKey="Paolotti D">D Paolotti</name>
</author>
<author><name sortKey="Perra, N" uniqKey="Perra N">N Perra</name>
</author>
<author><name sortKey="Tizzoni, M" uniqKey="Tizzoni M">M Tizzoni</name>
</author>
<author><name sortKey="Van Den Broeck, W" uniqKey="Van Den Broeck W">W Van Den Broeck</name>
</author>
<author><name sortKey="Colizza, V" uniqKey="Colizza V">V Colizza</name>
</author>
<author><name sortKey="Vespignani, A" uniqKey="Vespignani A">A Vespignani</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nishiura, H" uniqKey="Nishiura H">H Nishiura</name>
</author>
<author><name sortKey="Castillo Chavez, C" uniqKey="Castillo Chavez C">C Castillo-Chavez</name>
</author>
<author><name sortKey="Safan, M" uniqKey="Safan M">M Safan</name>
</author>
<author><name sortKey="Chowell, G" uniqKey="Chowell G">G Chowell</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Cutler, J" uniqKey="Cutler J">J Cutler</name>
</author>
<author><name sortKey="Schleihauf, E" uniqKey="Schleihauf E">E Schleihauf</name>
</author>
<author><name sortKey="Hatchette, Tf" uniqKey="Hatchette T">TF Hatchette</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="Stanford, E" uniqKey="Stanford E">E Stanford</name>
</author>
<author><name sortKey="Andrews, N" uniqKey="Andrews N">N Andrews</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hahne, S" uniqKey="Hahne S">S Hahné</name>
</author>
<author><name sortKey="Donker, T" uniqKey="Donker T">T Donker</name>
</author>
<author><name sortKey="Meijer, A" uniqKey="Meijer A">A Meijer</name>
</author>
<author><name sortKey="Timen, A" uniqKey="Timen A">A Timen</name>
</author>
<author><name sortKey="Van Steenbergen, J" uniqKey="Van Steenbergen J">J Van Steenbergen</name>
</author>
<author><name sortKey="Osterhaus, A" uniqKey="Osterhaus A">A Osterhaus</name>
</author>
<author><name sortKey="Van Der Sande, M" uniqKey="Van Der Sande M">M Van Der Sande</name>
</author>
<author><name sortKey="Koopmans, M" uniqKey="Koopmans M">M Koopmans</name>
</author>
<author><name sortKey="Wallinga, J" uniqKey="Wallinga J">J Wallinga</name>
</author>
<author><name sortKey="Coutinho, R" uniqKey="Coutinho R">R Coutinho</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Rizzo, C" uniqKey="Rizzo C">C Rizzo</name>
</author>
<author><name sortKey="Declich, S" uniqKey="Declich S">S Declich</name>
</author>
<author><name sortKey="Bella, A" uniqKey="Bella A">A Bella</name>
</author>
<author><name sortKey="Caporali, Mg" uniqKey="Caporali M">MG Caporali</name>
</author>
<author><name sortKey="Lana, S" uniqKey="Lana S">S Lana</name>
</author>
<author><name sortKey="Pompa, Mg" uniqKey="Pompa M">MG Pompa</name>
</author>
<author><name sortKey="Vellucci, L" uniqKey="Vellucci L">L Vellucci</name>
</author>
<author><name sortKey="Salmaso, S" uniqKey="Salmaso S">S Salmaso</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nishiura, H" uniqKey="Nishiura H">H Nishiura</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nishiura, H" uniqKey="Nishiura H">H Nishiura</name>
</author>
<author><name sortKey="Cook, Ar" uniqKey="Cook A">AR Cook</name>
</author>
<author><name sortKey="Cowling, Bj" uniqKey="Cowling B">BJ Cowling</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Barrat, A" uniqKey="Barrat A">A Barrat</name>
</author>
<author><name sortKey="Barthelemy, M" uniqKey="Barthelemy M">M Barthelemy</name>
</author>
<author><name sortKey="Pastor Satorras, R" uniqKey="Pastor Satorras R">R Pastor-Satorras</name>
</author>
<author><name sortKey="Vespignani, A" uniqKey="Vespignani A">A Vespignani</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Guimera, R" uniqKey="Guimera R">R Guimera</name>
</author>
<author><name sortKey="Mossa, S" uniqKey="Mossa S">S Mossa</name>
</author>
<author><name sortKey="Turtschi, A" uniqKey="Turtschi A">A Turtschi</name>
</author>
<author><name sortKey="Amaral, Lan" uniqKey="Amaral L">LAN Amaral</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Demontis, A" uniqKey="Demontis A">A DeMontis</name>
</author>
<author><name sortKey="Barthelemy, M" uniqKey="Barthelemy M">M Barthelemy</name>
</author>
<author><name sortKey="Chessa, A" uniqKey="Chessa A">A Chessa</name>
</author>
<author><name sortKey="Vespignani, A" uniqKey="Vespignani A">A Vespignani</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chowell, G" uniqKey="Chowell G">G Chowell</name>
</author>
<author><name sortKey="Hyman, Jm" uniqKey="Hyman J">JM Hyman</name>
</author>
<author><name sortKey="Eubank, S" uniqKey="Eubank S">S Eubank</name>
</author>
<author><name sortKey="Castillo Chavez, C" uniqKey="Castillo Chavez C">C Castillo-Chavez</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Barrett, Cl" uniqKey="Barrett C">CL Barrett</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Mossong, J" uniqKey="Mossong J">J Mossong</name>
</author>
<author><name sortKey="Hens, N" uniqKey="Hens N">N Hens</name>
</author>
<author><name sortKey="Jit, M" uniqKey="Jit M">M Jit</name>
</author>
<author><name sortKey="Beutels, P" uniqKey="Beutels P">P Beutels</name>
</author>
<author><name sortKey="Auranen, K" uniqKey="Auranen K">K Auranen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Anderson, Rm" uniqKey="Anderson R">RM Anderson</name>
</author>
<author><name sortKey="May, Rm" uniqKey="May R">RM May</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wallinga, J" uniqKey="Wallinga J">J Wallinga</name>
</author>
<author><name sortKey="Teunis, P" uniqKey="Teunis P">P Teunis</name>
</author>
<author><name sortKey="Kretzschmar, M" uniqKey="Kretzschmar M">M Kretzschmar</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Del Valle, Sy" uniqKey="Del Valle S">SY Del Valle</name>
</author>
<author><name sortKey="Hyman, J" uniqKey="Hyman J">J Hyman</name>
</author>
<author><name sortKey="Hethcote, Hw" uniqKey="Hethcote H">HW Hethcote</name>
</author>
<author><name sortKey="Eubank, Sg" uniqKey="Eubank S">SG Eubank</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kiss, I" uniqKey="Kiss I">I Kiss</name>
</author>
<author><name sortKey="Simon, Pl" uniqKey="Simon P">PL Simon</name>
</author>
<author><name sortKey="Kao, R" uniqKey="Kao R">R Kao</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Goeyvaerts, N" uniqKey="Goeyvaerts N">N Goeyvaerts</name>
</author>
<author><name sortKey="Hens, N" uniqKey="Hens N">N Hens</name>
</author>
<author><name sortKey="Aerts, M" uniqKey="Aerts M">M Aerts</name>
</author>
<author><name sortKey="Beutels, P" uniqKey="Beutels P">P Beutels</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Colizza, V" uniqKey="Colizza V">V Colizza</name>
</author>
<author><name sortKey="Vespignani, A" uniqKey="Vespignani A">A Vespignani</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pastor Satorras, R" uniqKey="Pastor Satorras R">R Pastor-Satorras</name>
</author>
<author><name sortKey="Vespignani, A" uniqKey="Vespignani A">A Vespignani</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Colizza, V" uniqKey="Colizza V">V Colizza</name>
</author>
<author><name sortKey="Pastor Satorras, R" uniqKey="Pastor Satorras R">R Pastor-Satorras</name>
</author>
<author><name sortKey="Vespignani, A" uniqKey="Vespignani A">A Vespignani</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Goeyvaerts, N" uniqKey="Goeyvaerts N">N Goeyvaerts</name>
</author>
<author><name sortKey="Hens, N" uniqKey="Hens N">N Hens</name>
</author>
<author><name sortKey="Ogunjimi, B" uniqKey="Ogunjimi B">B Ogunjimi</name>
</author>
<author><name sortKey="Aerts, M" uniqKey="Aerts M">M Aerts</name>
</author>
<author><name sortKey="Shkedy, Z" uniqKey="Shkedy Z">Z Shkedy</name>
</author>
<author><name sortKey="Van Damme, P" uniqKey="Van Damme P">P Van Damme</name>
</author>
<author><name sortKey="Beutels, P" uniqKey="Beutels P">P Beutels</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Goldstein, E" uniqKey="Goldstein E">E Goldstein</name>
</author>
<author><name sortKey="Apolloni, A" uniqKey="Apolloni A">A Apolloni</name>
</author>
<author><name sortKey="Lewis, B" uniqKey="Lewis B">B Lewis</name>
</author>
<author><name sortKey="Miller, Jc" uniqKey="Miller J">JC Miller</name>
</author>
<author><name sortKey="Macauley, M" uniqKey="Macauley M">M Macauley</name>
</author>
<author><name sortKey="Eubank, S" uniqKey="Eubank S">S Eubank</name>
</author>
<author><name sortKey="Lipsitch, Mj" uniqKey="Lipsitch M">MJ Lipsitch</name>
</author>
<author><name sortKey="Wallinga, J" uniqKey="Wallinga J">J Wallinga</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Diekmann, O" uniqKey="Diekmann O">O Diekmann</name>
</author>
<author><name sortKey="Heesterbeek, Jap" uniqKey="Heesterbeek J">JAP Heesterbeek</name>
</author>
<author><name sortKey="Roberts, Mg" uniqKey="Roberts M">MG Roberts</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Leung, Gm" uniqKey="Leung G">GM Leung</name>
</author>
<author><name sortKey="Hedley, Aj" uniqKey="Hedley A">AJ Hedley</name>
</author>
<author><name sortKey="Ho, Lm" uniqKey="Ho L">LM Ho</name>
</author>
<author><name sortKey="Chau, P" uniqKey="Chau P">P Chau</name>
</author>
<author><name sortKey="Wong, Io" uniqKey="Wong I">IO Wong</name>
</author>
<author><name sortKey="Thach, Tq" uniqKey="Thach T">TQ Thach</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nishiura, H" uniqKey="Nishiura H">H Nishiura</name>
</author>
<author><name sortKey="Chowell, G" uniqKey="Chowell G">G Chowell</name>
</author>
<author><name sortKey="Safan, M" uniqKey="Safan M">M Safan</name>
</author>
<author><name sortKey="Castillo Chavez, C" uniqKey="Castillo Chavez C">C Castillo-Chavez</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Eames, Ktd" uniqKey="Eames K">KTD Eames</name>
</author>
<author><name sortKey="Tilston, Nl" uniqKey="Tilston N">NL Tilston</name>
</author>
<author><name sortKey="Brooks Pollock, E" uniqKey="Brooks Pollock E">E Brooks-Pollock</name>
</author>
<author><name sortKey="Edmunds, Wj" uniqKey="Edmunds W">WJ Edmunds</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hens, N" uniqKey="Hens N">N Hens</name>
</author>
<author><name sortKey="Ayele, Gm" uniqKey="Ayele G">GM Ayele</name>
</author>
<author><name sortKey="Goeyvaerts, N" uniqKey="Goeyvaerts N">N Goeyvaerts</name>
</author>
<author><name sortKey="Aerts, M" uniqKey="Aerts M">M Aerts</name>
</author>
<author><name sortKey="Mossong, J" uniqKey="Mossong J">J Mossong</name>
</author>
<author><name sortKey="Edmunds, Jw" uniqKey="Edmunds J">JW Edmunds</name>
</author>
<author><name sortKey="Beutels, P" uniqKey="Beutels P">P Beutels</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Diekmann, O" uniqKey="Diekmann O">O Diekmann</name>
</author>
<author><name sortKey="Heesterbeek, Jap" uniqKey="Heesterbeek J">JAP Heesterbeek</name>
</author>
<author><name sortKey="Metz, Jaj" uniqKey="Metz J">JAJ Metz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bajardi, P" uniqKey="Bajardi P">P Bajardi</name>
</author>
<author><name sortKey="Poletto, C" uniqKey="Poletto C">C Poletto</name>
</author>
<author><name sortKey="Ramasco, Jj" uniqKey="Ramasco J">JJ Ramasco</name>
</author>
<author><name sortKey="Tizzoni, M" uniqKey="Tizzoni M">M Tizzoni</name>
</author>
<author><name sortKey="Colizza, V" uniqKey="Colizza V">V Colizza</name>
</author>
<author><name sortKey="Vespignani, A" uniqKey="Vespignani A">A Vespignani</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ikonen, N" uniqKey="Ikonen N">N Ikonen</name>
</author>
<author><name sortKey="Strengell, M" uniqKey="Strengell M">M Strengell</name>
</author>
<author><name sortKey="Kinnunen, L" uniqKey="Kinnunen L">L Kinnunen</name>
</author>
<author><name sortKey="Osterlund, P" uniqKey="Osterlund P">P Osterlund</name>
</author>
<author><name sortKey="Pirhonen, J" uniqKey="Pirhonen J">J Pirhonen</name>
</author>
<author><name sortKey="Broman, M" uniqKey="Broman M">M Broman</name>
</author>
<author><name sortKey="Davidkin, I" uniqKey="Davidkin I">I Davidkin</name>
</author>
<author><name sortKey="Ziegler, T" uniqKey="Ziegler T">T Ziegler</name>
</author>
<author><name sortKey="Julkunen, I" uniqKey="Julkunen I">I Julkunen</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>
<author><name sortKey="Sun, H" uniqKey="Sun H">H Sun</name>
</author>
<author><name sortKey="Liu, F" uniqKey="Liu F">F Liu</name>
</author>
<author><name sortKey="Dong, L" uniqKey="Dong L">L Dong</name>
</author>
<author><name sortKey="Devos, Jr" uniqKey="Devos J">JR DeVos</name>
</author>
<author><name sortKey="Gargiullo, Pm" uniqKey="Gargiullo P">PM Gargiullo</name>
</author>
<author><name sortKey="Brammer, Tl" uniqKey="Brammer T">TL Brammer</name>
</author>
<author><name sortKey="Cox, Nj" uniqKey="Cox N">NJ Cox</name>
</author>
<author><name sortKey="Tumpey, Tm" uniqKey="Tumpey T">TM Tumpey</name>
</author>
<author><name sortKey="Katz, Jm" uniqKey="Katz J">JM Katz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Allwinn, R" uniqKey="Allwinn R">R Allwinn</name>
</author>
<author><name sortKey="Geiler, J" uniqKey="Geiler J">J Geiler</name>
</author>
<author><name sortKey="Berger, A" uniqKey="Berger A">A Berger</name>
</author>
<author><name sortKey="Cinatl, J" uniqKey="Cinatl J">J Cinatl</name>
</author>
<author><name sortKey="Doer, Hw" uniqKey="Doer H">HW Doer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ball, F" uniqKey="Ball F">F Ball</name>
</author>
<author><name sortKey="Mollison, D" uniqKey="Mollison D">D Mollison</name>
</author>
<author><name sortKey="Scalia Tomba, G" uniqKey="Scalia Tomba G">G Scalia-Tomba</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cross, P" uniqKey="Cross P">P Cross</name>
</author>
<author><name sortKey="Lloyd Smith, Jo" uniqKey="Lloyd Smith J">JO Lloyd-Smith</name>
</author>
<author><name sortKey="Johnson, Plf" uniqKey="Johnson P">PLF Johnson</name>
</author>
<author><name sortKey="Wayne, Mg" uniqKey="Wayne M">MG Wayne</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cross, P" uniqKey="Cross P">P Cross</name>
</author>
<author><name sortKey="Johnson, Plf" uniqKey="Johnson P">PLF Johnson</name>
</author>
<author><name sortKey="Lloyd Smith, Jo" uniqKey="Lloyd Smith J">JO Lloyd-Smith</name>
</author>
<author><name sortKey="Wayne, Mg" uniqKey="Wayne M">MG Wayne</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Colizza, V" uniqKey="Colizza V">V Colizza</name>
</author>
<author><name sortKey="Vespignani, A" uniqKey="Vespignani A">A Vespignani</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Balcan, D" uniqKey="Balcan D">D Balcan</name>
</author>
<author><name sortKey="Vespignani, A" uniqKey="Vespignani A">A Vespignani</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Belik, V" uniqKey="Belik V">V Belik</name>
</author>
<author><name sortKey="Geisel, T" uniqKey="Geisel T">T Geisel</name>
</author>
<author><name sortKey="Brockmann, D" uniqKey="Brockmann D">D Brockmann</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Meloni, S" uniqKey="Meloni S">S Meloni</name>
</author>
<author><name sortKey="Perra, N" uniqKey="Perra N">N Perra</name>
</author>
<author><name sortKey="Arenas, A" uniqKey="Arenas A">A Arenas</name>
</author>
<author><name sortKey="Gomes, S" uniqKey="Gomes S">S Gomes</name>
</author>
<author><name sortKey="Moreno, Y" uniqKey="Moreno Y">Y Moreno</name>
</author>
<author><name sortKey="Vespignani, A" uniqKey="Vespignani A">A Vespignani</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Poletto, C" uniqKey="Poletto C">C Poletto</name>
</author>
<author><name sortKey="Tizzoni, M" uniqKey="Tizzoni M">M Tizzoni</name>
</author>
<author><name sortKey="Colizza, V" uniqKey="Colizza V">V Colizza</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Harris, Te" uniqKey="Harris T">TE Harris</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ball, F" uniqKey="Ball F">F Ball</name>
</author>
<author><name sortKey="Clancy, D" uniqKey="Clancy D">D Clancy</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Flasche, S" uniqKey="Flasche S">S Flasche</name>
</author>
<author><name sortKey="Hens, N" uniqKey="Hens N">N Hens</name>
</author>
<author><name sortKey="Boelle, Py" uniqKey="Boelle P">PY Boëlle</name>
</author>
<author><name sortKey="Mossong, J" uniqKey="Mossong J">J Mossong</name>
</author>
<author><name sortKey="Van Ballegooijen, Wm" uniqKey="Van Ballegooijen W">WM van Ballegooijen</name>
</author>
<author><name sortKey="Nunes, B" uniqKey="Nunes B">B Nunes</name>
</author>
<author><name sortKey="Rizzo, C" uniqKey="Rizzo C">C Rizzo</name>
</author>
<author><name sortKey="Popovici, F" uniqKey="Popovici F">F Popovici</name>
</author>
<author><name sortKey="Santa Olalla, P" uniqKey="Santa Olalla P">P Santa-Olalla</name>
</author>
<author><name sortKey="Hruba, F" uniqKey="Hruba F">F Hrubá</name>
</author>
<author><name sortKey="Parmakova, K" uniqKey="Parmakova K">K Parmakova</name>
</author>
<author><name sortKey="Baguelin, M" uniqKey="Baguelin M">M Baguelin</name>
</author>
<author><name sortKey="Van Hoek, Aj" uniqKey="Van Hoek A">AJ van Hoek</name>
</author>
<author><name sortKey="Desenclos, Jc" uniqKey="Desenclos J">JC Desenclos</name>
</author>
<author><name sortKey="Bernillon, P" uniqKey="Bernillon P">P Bernillon</name>
</author>
<author><name sortKey="Camara, Al" uniqKey="Camara A">AL Cámara</name>
</author>
<author><name sortKey="Wallinga, J" uniqKey="Wallinga J">J Wallinga</name>
</author>
<author><name sortKey="Asikainen, T" uniqKey="Asikainen T">T Asikainen</name>
</author>
<author><name sortKey="White, Pj" uniqKey="White P">PJ White</name>
</author>
<author><name sortKey="Edmunds, Wj" uniqKey="Edmunds W">WJ Edmunds</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mills, Ce" uniqKey="Mills C">CE Mills</name>
</author>
<author><name sortKey="Robins, Jm" uniqKey="Robins J">JM Robins</name>
</author>
<author><name sortKey="Lipsitch, M" uniqKey="Lipsitch M">M Lipsitch</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="Obenchain, Vj" uniqKey="Obenchain V">VJ Obenchain</name>
</author>
<author><name sortKey="Longini, Im" uniqKey="Longini I">IM Longini</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Watts, Dj" uniqKey="Watts D">DJ Watts</name>
</author>
<author><name sortKey="Muhamad, R" uniqKey="Muhamad R">R Muhamad</name>
</author>
<author><name sortKey="Medina, Dc" uniqKey="Medina D">DC Medina</name>
</author>
<author><name sortKey="Dodds, Ps" uniqKey="Dodds P">PS Dodds</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Epstein, Jm" uniqKey="Epstein J">JM Epstein</name>
</author>
<author><name sortKey="Goedecke, Dm" uniqKey="Goedecke D">DM Goedecke</name>
</author>
<author><name sortKey="Yu, F" uniqKey="Yu F">F Yu</name>
</author>
<author><name sortKey="Morris, Rj" uniqKey="Morris R">RJ Morris</name>
</author>
<author><name sortKey="Wagener, Dk" uniqKey="Wagener D">DK Wagener</name>
</author>
<author><name sortKey="Bobashev, Gv" uniqKey="Bobashev G">GV Bobashev</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cooper, Bs" uniqKey="Cooper B">BS Cooper</name>
</author>
<author><name sortKey="Pitman, Rj" uniqKey="Pitman R">RJ Pitman</name>
</author>
<author><name sortKey="Edmunds, Wj" uniqKey="Edmunds W">WJ Edmunds</name>
</author>
<author><name sortKey="Gay, Nj" uniqKey="Gay N">NJ Gay</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hollingsworth, Td" uniqKey="Hollingsworth T">TD Hollingsworth</name>
</author>
<author><name sortKey="Ferguson, Nm" uniqKey="Ferguson N">NM Ferguson</name>
</author>
<author><name sortKey="Anderson, Rm" uniqKey="Anderson R">RM Anderson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lam, Ehy" uniqKey="Lam E">EHY Lam</name>
</author>
<author><name sortKey="Cowling, Bj" uniqKey="Cowling B">BJ Cowling</name>
</author>
<author><name sortKey="Cook, Ar" uniqKey="Cook A">AR Cook</name>
</author>
<author><name sortKey="Wong, Jyt" uniqKey="Wong J">JYT Wong</name>
</author>
<author><name sortKey="Lau, Msy" uniqKey="Lau M">MSY Lau</name>
</author>
<author><name sortKey="Nishiura, H" uniqKey="Nishiura H">H Nishiura</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Merler, S" uniqKey="Merler S">S Merler</name>
</author>
<author><name sortKey="Ajelli, M" uniqKey="Ajelli M">M Ajelli</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Cowling, Bj" uniqKey="Cowling B">BJ Cowling</name>
</author>
<author><name sortKey="Lau, Llh" uniqKey="Lau L">LLH Lau</name>
</author>
<author><name sortKey="Wu, P" uniqKey="Wu P">P Wu</name>
</author>
<author><name sortKey="Wong, Hwc" uniqKey="Wong H">HWC Wong</name>
</author>
<author><name sortKey="Fang, Vj" uniqKey="Fang V">VJ Fang</name>
</author>
<author><name sortKey="Riley, S" uniqKey="Riley S">S Riley</name>
</author>
<author><name sortKey="Nishiura, H" uniqKey="Nishiura H">H Nishiura</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Eubank, S" uniqKey="Eubank S">S Eubank</name>
</author>
<author><name sortKey="Guclu, H" uniqKey="Guclu H">H Guclu</name>
</author>
<author><name sortKey="Anil Kumar, Vs" uniqKey="Anil Kumar V">VS Anil Kumar</name>
</author>
<author><name sortKey="Marathe, Mv" uniqKey="Marathe M">MV Marathe</name>
</author>
<author><name sortKey="Srinivasan, A" uniqKey="Srinivasan A">A Srinivasan</name>
</author>
<author><name sortKey="Toroczkai, Z" uniqKey="Toroczkai Z">Z Toroczkai</name>
</author>
<author><name sortKey="Wang, N" uniqKey="Wang N">N Wang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Halloran, Me" uniqKey="Halloran M">ME Halloran</name>
</author>
<author><name sortKey="Ferguson, Nm" uniqKey="Ferguson N">NM Ferguson</name>
</author>
<author><name sortKey="Eubank, S" uniqKey="Eubank S">S Eubank</name>
</author>
<author><name sortKey="Longini, Im" uniqKey="Longini I">IM Longini</name>
</author>
<author><name sortKey="Cummings, Dat" uniqKey="Cummings D">DAT Cummings</name>
</author>
<author><name sortKey="Lewis, B" uniqKey="Lewis B">B Lewis</name>
</author>
<author><name sortKey="Xu, S" uniqKey="Xu S">S Xu</name>
</author>
<author><name sortKey="Fraser, C" uniqKey="Fraser C">C Fraser</name>
</author>
<author><name sortKey="Vullikanti, A" uniqKey="Vullikanti A">A Vullikanti</name>
</author>
<author><name sortKey="Germann, Tc" uniqKey="Germann T">TC Germann</name>
</author>
<author><name sortKey="Wagener, D" uniqKey="Wagener D">D Wagener</name>
</author>
<author><name sortKey="Beckman, R" uniqKey="Beckman R">R Beckman</name>
</author>
<author><name sortKey="Kadau, K" uniqKey="Kadau K">K Kadau</name>
</author>
<author><name sortKey="Macken, Bds" uniqKey="Macken B">BDS Macken</name>
</author>
<author><name sortKey="Cooley, P" uniqKey="Cooley P">P Cooley</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Keeling, Mj" uniqKey="Keeling M">MJ Keeling</name>
</author>
<author><name sortKey="Grenfell, Bt" uniqKey="Grenfell B">BT Grenfell</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Vergu, E" uniqKey="Vergu E">E Vergu</name>
</author>
<author><name sortKey="Busson, H" uniqKey="Busson H">H Busson</name>
</author>
<author><name sortKey="Ezanno, P" uniqKey="Ezanno P">P Ezanno</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tizzoni, M" uniqKey="Tizzoni M">M Tizzoni</name>
</author>
<author><name sortKey="Bajardi, P" uniqKey="Bajardi P">P Bajardi</name>
</author>
<author><name sortKey="Poletto, C" uniqKey="Poletto C">C Poletto</name>
</author>
<author><name sortKey="Ramasco, Jj" uniqKey="Ramasco J">JJ Ramasco</name>
</author>
<author><name sortKey="Balcan, D" uniqKey="Balcan D">D Balcan</name>
</author>
<author><name sortKey="Goncalves, B" uniqKey="Goncalves B">B Goncalves</name>
</author>
<author><name sortKey="Perra, N" uniqKey="Perra N">N Perra</name>
</author>
<author><name sortKey="Colizza, V" uniqKey="Colizza V">V Colizza</name>
</author>
<author><name sortKey="Vespignani, A" uniqKey="Vespignani A">A Vespignani</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Funk, S" uniqKey="Funk S">S Funk</name>
</author>
<author><name sortKey="Salathe, M" uniqKey="Salathe M">M Salathe</name>
</author>
<author><name sortKey="Jansen, Vaa" uniqKey="Jansen V">VAA Jansen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Salathe, M" uniqKey="Salathe M">M Salathé</name>
</author>
<author><name sortKey="Khandelwal, S" uniqKey="Khandelwal S">S Khandelwal</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Perra, N" uniqKey="Perra N">N Perra</name>
</author>
<author><name sortKey="Balcan, D" uniqKey="Balcan D">D Balcan</name>
</author>
<author><name sortKey="Goncalves, B" uniqKey="Goncalves B">B Goncalves</name>
</author>
<author><name sortKey="Vespignani, A" uniqKey="Vespignani A">A Vespignani</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bagnoli, F" uniqKey="Bagnoli F">F Bagnoli</name>
</author>
<author><name sortKey="Lio, P" uniqKey="Lio P">P Lio</name>
</author>
<author><name sortKey="Sguanci, L" uniqKey="Sguanci L">L Sguanci</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Poletti, P" uniqKey="Poletti P">P Poletti</name>
</author>
<author><name sortKey="Caprile, B" uniqKey="Caprile B">B Caprile</name>
</author>
<author><name sortKey="Ajelli, M" uniqKey="Ajelli M">M Ajelli</name>
</author>
<author><name sortKey="Pugliese, A" uniqKey="Pugliese A">A Pugliese</name>
</author>
<author><name sortKey="Merler, S" uniqKey="Merler S">S Merler</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Fumanelli, L" uniqKey="Fumanelli L">L Fumanelli</name>
</author>
<author><name sortKey="Ajelli, M" uniqKey="Ajelli M">M Ajelli</name>
</author>
<author><name sortKey="Manfredi, P" uniqKey="Manfredi P">P Manfredi</name>
</author>
<author><name sortKey="Vespignani, A" uniqKey="Vespignani A">A Vespignani</name>
</author>
<author><name sortKey="Merler, S" uniqKey="Merler S">S Merler</name>
</author>
</analytic>
</biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">BMC Infect Dis</journal-id>
<journal-id journal-id-type="iso-abbrev">BMC Infect. Dis</journal-id>
<journal-title-group><journal-title>BMC Infectious Diseases</journal-title>
</journal-title-group>
<issn pub-type="epub">1471-2334</issn>
<publisher><publisher-name>BioMed Central</publisher-name>
<publisher-loc>London</publisher-loc>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">23587010</article-id>
<article-id pub-id-type="pmc">3644502</article-id>
<article-id pub-id-type="publisher-id">2399</article-id>
<article-id pub-id-type="doi">10.1186/1471-2334-13-176</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject>
</subj-group>
</article-categories>
<title-group><article-title>Age-specific contacts and travel patterns in the spatial spread of 2009 H1N1 influenza pandemic</article-title>
</title-group>
<contrib-group><contrib contrib-type="author" equal-contrib="yes"><name><surname>Apolloni</surname>
<given-names>Andrea</given-names>
</name>
<address><email>andrea.apolloni@lshtm.ac.uk</email>
</address>
<xref ref-type="aff" rid="Aff1">1</xref>
</contrib>
<contrib contrib-type="author" equal-contrib="yes"><name><surname>Poletto</surname>
<given-names>Chiara</given-names>
</name>
<address><email>chiara.poletto@isi.it</email>
</address>
<xref ref-type="aff" rid="Aff2">2</xref>
<xref ref-type="aff" rid="Aff3">3</xref>
<xref ref-type="aff" rid="Aff4">4</xref>
</contrib>
<contrib contrib-type="author" corresp="yes"><name><surname>Colizza</surname>
<given-names>Vittoria</given-names>
</name>
<address><email>vittoria.colizza@inserm.fr</email>
</address>
<xref ref-type="aff" rid="Aff3">3</xref>
<xref ref-type="aff" rid="Aff4">4</xref>
<xref ref-type="aff" rid="Aff5">5</xref>
</contrib>
<aff id="Aff1"><label>1</label>
<institution-wrap><institution-id institution-id-type="GRID">grid.8991.9</institution-id>
<institution-id institution-id-type="ISNI">000000040425469X</institution-id>
<institution>Department of Infectious Disease Epidemiology,</institution>
<institution>London School of Hygiene and Tropical Medicine,</institution>
</institution-wrap>
London, UK</aff>
<aff id="Aff2"><label>2</label>
<institution-wrap><institution-id institution-id-type="GRID">grid.418750.f</institution-id>
<institution-id institution-id-type="ISNI">0000000417593658</institution-id>
<institution>Computational Epidemiology Laboratory,</institution>
<institution>Institute for Scientific Interchange (ISI),</institution>
</institution-wrap>
Torino, Italy</aff>
<aff id="Aff3"><label>3</label>
<institution-wrap><institution-id institution-id-type="GRID">grid.7429.8</institution-id>
<institution-id institution-id-type="ISNI">0000000121866389</institution-id>
<institution>INSERM, U707,</institution>
</institution-wrap>
Paris, France</aff>
<aff id="Aff4"><label>4</label>
<institution-wrap><institution-id institution-id-type="GRID">grid.5805.8</institution-id>
<institution-id institution-id-type="ISNI">0000000119553500</institution-id>
<institution>Faculté de Médecine Pierre et Marie Curie,</institution>
<institution>UPMC Université Paris 06, UMR S 707,</institution>
</institution-wrap>
Paris, France</aff>
<aff id="Aff5"><label>5</label>
<institution-wrap><institution-id institution-id-type="GRID">grid.418750.f</institution-id>
<institution-id institution-id-type="ISNI">0000000417593658</institution-id>
<institution>Institute for Scientific Interchange (ISI),</institution>
</institution-wrap>
Torino, Italy</aff>
</contrib-group>
<pub-date pub-type="epub"><day>15</day>
<month>4</month>
<year>2013</year>
</pub-date>
<pub-date pub-type="pmc-release"><day>15</day>
<month>4</month>
<year>2013</year>
</pub-date>
<pub-date pub-type="collection"><year>2013</year>
</pub-date>
<volume>13</volume>
<elocation-id>176</elocation-id>
<history><date date-type="received"><day>31</day>
<month>12</month>
<year>2012</year>
</date>
<date date-type="accepted"><day>27</day>
<month>3</month>
<year>2013</year>
</date>
</history>
<permissions><copyright-statement>© Apolloni et al.; licensee BioMed Central Ltd. 2013</copyright-statement>
<license><license-p>This article is made available via the PMC Open Access Subset for unrestricted research re-use and secondary analysis in any form or by any means with acknowledgement of the original source. These permissions are granted for the duration of the World Health Organization (WHO) declaration of COVID-19 as a global pandemic.</license-p>
</license>
</permissions>
<abstract id="Abs1"><sec><title>Background</title>
<p>Confirmed H1N1 cases during late spring and summer 2009 in various countries showed a substantial age shift between importations and local transmission cases, with adults mainly responsible for seeding unaffected regions and children most frequently driving community outbreaks.</p>
</sec>
<sec><title>Methods</title>
<p>We introduce a multi-host stochastic metapopulation model with two age classes to analytically investigate the role of a heterogeneously mixing population and its associated non-homogeneous travel behaviors on the risk of a major epidemic. We inform the model with demographic data, contact data and travel statistics of Europe and Mexico, and calibrate it to the 2009 H1N1 pandemic early outbreak. We allow for variations of the model parameters to explore the conditions of invasion under different scenarios.</p>
</sec>
<sec><title>Results</title>
<p>We derive the expression for the potential of global invasion of the epidemic that depends on the transmissibility of the pathogen, the transportation network and mobility features, the demographic profile and the mixing pattern. Higher assortativity in the contact pattern greatly increases the probability of spatial containment of the epidemic, this effect being contrasted by an increase in the social activity of adults vs. children. Heterogeneous features of the mobility network characterizing its topology and traffic flows strongly favor the invasion of the pathogen at the spatial level, as also a larger fraction of children traveling. Variations in the demographic profile and mixing habits across countries lead to heterogeneous outbreak situations. Model results are compatible with the H1N1 spatial transmission dynamics observed.</p>
</sec>
<sec><title>Conclusions</title>
<p>This work illustrates the importance of considering age-dependent mixing profiles and mobility features coupled together to study the conditions for the spatial invasion of an emerging influenza pandemic. Its results allow the immediate assessment of the risk of a major epidemic for a specific scenario upon availability of data, and the evaluation of the potential effectiveness of public health interventions targeting specific age groups, their interactions and mobility behaviors. The approach provides a general modeling framework that can be used for other types of partitions of the host population and applied to different settings.</p>
</sec>
<sec><title>Electronic supplementary material</title>
<p>The online version of this article (doi:10.1186/1471-2334-13-176) contains supplementary material, which is available to authorized users.</p>
</sec>
</abstract>
<kwd-group xml:lang="en"><title>Keywords</title>
<kwd>H1N1 pandemic</kwd>
<kwd>Host heterogeneities</kwd>
<kwd>Spatial spread</kwd>
<kwd>Age structure</kwd>
<kwd>Mobility</kwd>
<kwd>Mathematical modeling</kwd>
</kwd-group>
<custom-meta-group><custom-meta><meta-name>issue-copyright-statement</meta-name>
<meta-value>© The Author(s) 2013</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
</front>
<body><sec id="Sec1"><title>Background</title>
<p>The data collected during and after the 2009 H1N1 pandemic has contributed to achieve major insights regarding key factors of the transmission dynamics of the novel strain of influenza. Two aspects emerging from surveillance and serological data during the initial phase of the outbreak became strikingly clear: (i) international movements of passengers by air travel drove the spatial dissemination of the pathogen at the global level [<xref ref-type="bibr" rid="CR1">1</xref>
–<xref ref-type="bibr" rid="CR3">3</xref>
], and (ii) initial local epidemics mainly occurred in schools [<xref ref-type="bibr" rid="CR2">2</xref>
, <xref ref-type="bibr" rid="CR4">4</xref>
–<xref ref-type="bibr" rid="CR7">7</xref>
]. Country surveillance data from summer 2009 [<xref ref-type="bibr" rid="CR8">8</xref>
–<xref ref-type="bibr" rid="CR13">13</xref>
] show a dramatic difference in the age distributions of imported cases and indigenous cases, with imported cases on average older than the ones generated by local transmission (see Figure <xref rid="Fig1" ref-type="fig">1</xref>
A). This result strengthens the previous observations and highlights the presence of an age shift between the seeding events established by travelers from affected areas – mainly adults – and the local outbreaks mainly occurring among school-aged children. Travel statistics indeed confirm that the vast majority of passengers flying are aged ≥18 years (see Figure <xref rid="Fig1" ref-type="fig">1</xref>
B).<fig id="Fig1"><label>Figure 1</label>
<caption><p><bold>Imported vs. indigenous H1N1 cases and age-specific travel statistics for various countries.</bold>
(<bold>A</bold>
) Fraction of indigenous cases and of imported cases during the initial phase of the H1N1 pandemic outbreak in the [0–19] years age class, calculated from surveillance data for the following countries: The Netherlands [<xref ref-type="bibr" rid="CR8">8</xref>
], Belgium [<xref ref-type="bibr" rid="CR9">9</xref>
], UK [<xref ref-type="bibr" rid="CR10">10</xref>
], France [<xref ref-type="bibr" rid="CR11">11</xref>
], Japan [<xref ref-type="bibr" rid="CR13">13</xref>
], Italy [<xref ref-type="bibr" rid="CR12">12</xref>
]. (<bold>B</bold>
) Percentage of air-travel passengers in the younger age classes for a set of airports around the world. The age classification used by the demographic statistics vary across countries (Helsinki<sup>1</sup>
, Finland; Teheran<sup>3</sup>
, Iran; Los Angeles<sup>2</sup>
, USA; Amsterdam<sup>5</sup>
, The Netherlands; Heathrow<sup>4</sup>
, Gatwick<sup>4</sup>
, Stansted<sup>4</sup>
, Luton<sup>4</sup>
, UK; Venice<sup>6</sup>
, Italy; Hannover<sup>7</sup>
, Frankfurt<sup>8</sup>
, Hamburg<sup>8</sup>
, Munich<sup>8</sup>
, Germany) with the corresponding age brackets for the children class (expressed in years): 1=[0,15]; 2=[0,18]; 3=[0,19]; 4=[0,20]; 5=[0,21]; 6=[0,25]; 7=[0,26]; 8=[0,30]. Sources of the data are reported in the Additional file <xref rid="MOESM1" ref-type="media">1</xref>
. The statistics found for Italy and Germany correspond to larger age brackets. If we rescale the data as indicated in the Demographic and travel data subsection, we obtain the following estimates for the percentage of travelers in the [0–18] years old class: 1.05% (Venice), 0.49% (Hannover), 2.31% (Frankfurt), 2.31% (Hamburg), and 2.86% (Munich). These estimated values are consistent with the data.</p>
</caption>
<graphic xlink:href="12879_2012_Article_2399_Fig1_HTML" id="d29e405"></graphic>
</fig>
</p>
<p>The role of children or adults in driving a major epidemic can be assessed with a simplified modeling approach expressed into two age classes and quantifying the probability of temporary extinction of seed individuals in each class depending on their mixing patterns [<xref ref-type="bibr" rid="CR14">14</xref>
]. Here we aim at fully addressing the interplay between the two factors emerged as empirical evidence during the initial outbreak of H1N1 pandemic – namely, age-specific seeding events and age-specific epidemic dynamics – by considering a spatially explicit model that integrates the mobility patterns of flights connecting different populated areas, as well as age-dependent travel behavior in addition to age-specific mixing. We introduce a multi-host stochastic metapopulation model that includes heterogeneous features in: the spatial structure of the population; the travel behavior of individuals depending on their age; the corresponding mixing patterns. We assume a simple two age-groups classification as this allows the analytical treatment of the model to obtain the expression of the conditions for a major epidemic in terms of the age-specific contacts and travel features. By exploring theoretical contact matrices with varying assortativity levels (i.e. within-group mixing), we assess the role of assortativity on the risk of a major epidemic and compare it to scenarios informed with estimates from the H1N1 pandemic. We find that the assortativity observed in real data tends to drive the system to extinction, and is counterbalanced by the heterogeneity of the air mobility network structure, by the travel behavior of adults and by the relative proportion of contacts established by adults, all aspects that favor the virus spread. Through a systematic exploration of the role of the various ingredients considered in the system, the presented results allow the risk assessment analysis of a specific epidemic scenario and can be extended to other infectious diseases where the population partition may play a relevant role.</p>
</sec>
<sec id="Sec2"><title>Methods</title>
<sec id="Sec3"><title>Demographic and travel data</title>
<p>We consider a metapopulation framework to simulate the spread of an infectious disease across subpopulations of individuals through mobility connections. The approach is generically applicable to various real-world systems and here we focus on modeling an emerging influenza pandemic across urban areas through air travel. We consider the distinct cases of 8 countries in Europe, and Mexico, for which data needed to inform the model are available.</p>
<p>The network specifying the coupling between different populations in real systems is in many cases very heterogeneous, and examples range from transportation infrastructures to mobility patterns of various type [<xref ref-type="bibr" rid="CR15">15</xref>
–<xref ref-type="bibr" rid="CR19">19</xref>
]. In the case of air travel the coupling is given by the direct flights connecting different airports and the number of passengers flying on those connections. Analyses of air transportation data have shown the presence of large variability in the number of connections per airport, and of broad fluctuations in the traffic handled by each airport or flowing on a given connection between origin and destination [<xref ref-type="bibr" rid="CR15">15</xref>
, <xref ref-type="bibr" rid="CR16">16</xref>
]. In the Additional File we report the example of the European air transportation network [<xref ref-type="bibr" rid="CR20">20</xref>
], showing the probability distributions of the number of connections <italic>k</italic>
per airport (i.e. the subpopulation’s degree) and of the flows of passengers <italic>w</italic>
<sub><italic>ij</italic>
</sub>
travelling between any pair of linked airports <italic>i</italic>
and <italic>j</italic>
. These flows can also be expressed in terms of the number of connections of the origin/destination airports, where the average weight <italic>w</italic>
<sub><italic>ij</italic>
</sub>
along the link connecting airports <italic>i</italic>
and <italic>j</italic>
is a function of their degrees <italic>k</italic>
<sub><italic>i</italic>
</sub>
and <italic>k</italic>
<sub><italic>j</italic>
</sub>
<italic>,</italic>
〈 <italic>w</italic>
<sub><italic>ij</italic>
</sub>
〉 ∝ (<italic>k</italic>
<sub><italic>i</italic>
</sub>
<italic>k</italic>
<sub><italic>j</italic>
</sub>
)<sup><italic>θ</italic>
</sup>
, with <italic>θ</italic>
≅ 0.5 in the worldwide air-transportation network [<xref ref-type="bibr" rid="CR15">15</xref>
]. Such features, obtained from empirical evidence, will be used in the following to define the metapopulation model, by creating a realistic synthetic network of mobility connecting a number <italic>V</italic>
of urban areas.</p>
<p>Since we are interested in exploring how non-homogeneous travel habits, coupled with non-homogeneous mixing patterns, may drive the conditions for the spatial invasion of an epidemic, we collected age travel statistics across different countries. These are typically obtained from travel surveys at airports, and collect a variety of information about passengers and their travel behavior including demographic data. Figure <xref rid="Fig1" ref-type="fig">1</xref>
B shows the percentage of travelers in the younger age class for a variety of sources and for different airports in various countries. If we consider a classification of the younger age group that includes individuals up to 18–21 years old (where the upper value of the range depends on the specific limits imposed by the age classification adopted by each survey), the fraction of children traveling by air is on average equal to 7%, with a maximum variation between 0.7% (observed for Teheran airport, Iran) and 9.2% (Luton airport, UK). Other statistical sources have larger age brackets, with no breakdown below 25 or 30 years old, as shown in the Figure. In these cases, we can still estimate the fraction of traveling individuals in the age class [0–18] years old, by rescaling the statistics assuming a constant ratio across countries between the percentages of travelers in the [0–18] years old and those of travelers in the [0–24] or [0–29] years old classes, as obtained from sources with a finer age classification. The estimated values are reported in the caption of Figure <xref rid="Fig1" ref-type="fig">1</xref>
B and are consistent with the data.</p>
<p>Demographic data for the age distribution of the population by country was obtained from Eurostat [<xref ref-type="bibr" rid="CR20">20</xref>
] for European countries, and from the U.N. database [<xref ref-type="bibr" rid="CR21">21</xref>
] for Mexico. The data is provided by yearly age groups for European countries and 5 years age groups for Mexico, and it allows the calculation of the fraction of people under a given age, corresponding to the classification used in the model. If we consider the younger age class up to 18 years old, we obtain relative sizes of the population ranging in Europe from 17.1% for Italy to 22.2% for United Kingdom and Luxembourg, with an average value of 19.7%. The corresponding value in Mexico increases to 32.3% [<xref ref-type="bibr" rid="CR21">21</xref>
], with a classification up to 15 years old for data availability reasons.</p>
</sec>
<sec id="Sec4"><title>Theoretical and data-driven age mixing patterns</title>
<p>In addition to the travel behavior and demographic features of the population, we need to consider the mixing pattern among population classes. For the purpose of the study, we consider data-driven mixing patterns by country [<xref ref-type="bibr" rid="CR2">2</xref>
, <xref ref-type="bibr" rid="CR14">14</xref>
, <xref ref-type="bibr" rid="CR22">22</xref>
] and we also define a theoretical contact matrix, dependent on a set of parameters that we vary in a range of plausible values in order to systematically explore the behavior of the global invasion under different mixing conditions.</p>
<p>We consider the population divided into two classes, children and adults identified by subscripts <italic>c</italic>
and <italic>α</italic>
, respectively. Our aim is indeed to characterize the impact on the invasion dynamics of the two coupled phenomena, namely the early phase of the epidemic outbreak locally driven mainly by children and the spatial dissemination of the epidemic mainly driven by the adults traveling. While demographic and travel data allow for the consideration of a larger number of age classes, our choice of two classes is motivated by the simpler parameterization of the model in terms of the mixing patterns that allows us to formulate the invasion conditions analytically.</p>
<p>Children are assumed to represent a fraction <italic>α</italic>
of the population <italic>N</italic>
, so that the population size of children is expressed as <italic>N</italic>
<sub><italic>c</italic>
</sub>
= <italic>αN;</italic>
the population size of adults being thus expressed as <italic>N</italic>
<sub><italic>a</italic>
</sub>
= (1 - <italic>α</italic>
)<italic>N</italic>
. <italic>α</italic>
is a parameter defined between 0 and 1, which we assume to be homogeneous across the metapopulation system. This is a plausible assumption if we consider a metapopulation model applied to a single country or to regions that are rather homogeneous in terms of their demography, as these values are not expected to substantially vary. In the following sections we will explore different values of <italic>α</italic>
, and compare these results to the ones obtained in the case the model is informed with the data-driven values of <italic>α</italic>
from the demographic statistics of the European countries under study and of Mexico.</p>
<p>We define the contact matrix <bold><italic>C</italic>
</bold>
= {<italic>C</italic>
<sub><italic>ij</italic>
</sub>
} capturing the mixing between different age classes as<disp-formula id="Equ1"><label>1</label>
<alternatives><mml:math id="M1"><mml:mrow><mml:mfenced close=")" open="("><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:mtext mathvariant="italic">cc</mml:mtext>
<mml:mspace width="0.5em"></mml:mspace>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:mtext mathvariant="italic">ca</mml:mtext>
<mml:mspace width="0.5em"></mml:mspace>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:msub><mml:mi>C</mml:mi>
<mml:mtext mathvariant="italic">ac</mml:mtext>
</mml:msub>
<mml:mspace width="0.5em"></mml:mspace>
<mml:msub><mml:mi>C</mml:mi>
<mml:mtext mathvariant="italic">aa</mml:mtext>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mfenced>
<mml:mspace width="0.5em"></mml:mspace>
<mml:mo>=</mml:mo>
<mml:mspace width="0.5em"></mml:mspace>
<mml:mfenced close=")" open="("><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mspace width="0.5em"></mml:mspace>
<mml:msub><mml:mi>p</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:msub><mml:mi>q</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mfrac><mml:mi>N</mml:mi>
<mml:msub><mml:mi>N</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
</mml:mfrac>
<mml:mspace width="2.8em"></mml:mspace>
<mml:mfenced close=")" open="("><mml:mrow><mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:msub><mml:mi>p</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
<mml:msub><mml:mi>q</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
<mml:mfrac><mml:mi>N</mml:mi>
<mml:msub><mml:mi>N</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
</mml:mfrac>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mfenced close=")" open="("><mml:mrow><mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:msub><mml:mi>p</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
<mml:msub><mml:mi>q</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mfrac><mml:mi>N</mml:mi>
<mml:msub><mml:mi>N</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
</mml:mfrac>
<mml:mspace width="2.8em"></mml:mspace>
<mml:msub><mml:mi>p</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
<mml:msub><mml:mi>q</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
<mml:mfrac><mml:mi>N</mml:mi>
<mml:msub><mml:mi>N</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
</mml:mfrac>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mspace width="2.9em"></mml:mspace>
</mml:mfenced>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<graphic xlink:href="12879_2012_Article_2399_Equ1_HTML.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
</p>
<p>where <italic>q</italic>
<sub><italic>c</italic>
</sub>
and <italic>q</italic>
<sub><italic>c</italic>
</sub>
are the average number of contacts per unit time established by individuals in the children and adult classes, respectively, and <italic>p</italic>
<sub><italic>c</italic>
</sub>
and <italic>p</italic>
<sub><italic>a</italic>
</sub>
are the fractions of contacts that occur between individuals of the same age class. The variables <italic>q</italic>
<sub><italic>c</italic>
</sub>
and <italic>q</italic>
<sub><italic>a</italic>
</sub>
represent a measure of social activity of the individuals, whereas <italic>p</italic>
<sub><italic>c</italic>
</sub>
and <italic>p</italic>
<sub><italic>a</italic>
</sub>
describe how these contacts are established among classes. A variety of different assumptions can be done to describe the social mixing pattern, as the variables of the contact matrix of the above expression are not uniquely defined and need to be parameterized through available demographic and serologic data [<xref ref-type="bibr" rid="CR23">23</xref>
]. Here we focus our attention on assortative mixing, indicating the tendency of individuals in a given class to preferably interact with other individuals of the same class. This is indeed observed in real data where contacts patterns are found to be highly assortative with age [<xref ref-type="bibr" rid="CR22">22</xref>
, <xref ref-type="bibr" rid="CR24">24</xref>
, <xref ref-type="bibr" rid="CR25">25</xref>
], with a remarkable similarity across different European countries [<xref ref-type="bibr" rid="CR22">22</xref>
]. By indicating with <italic>ε</italic>
<sub><italic>c</italic>
</sub>
(<italic>ε</italic>
<sub><italic>a</italic>
</sub>
) the average fraction of contacts that a child (adult) establishes with an adult (child), i.e. across age groups, we can simply express <italic>p</italic>
<sub><italic>c</italic>
</sub>
(<italic>p</italic>
<sub><italic>a</italic>
</sub>
) in terms of <italic>ε</italic>
<sub><italic>c</italic>
</sub>
(<italic>ε</italic>
<sub><italic>a</italic>
</sub>
) as <italic>p</italic>
<sub><italic>c</italic>
</sub>
= 1 - <italic>ε</italic>
<sub><italic>c</italic>
</sub>
(<italic>p</italic>
<sub><italic>a</italic>
</sub>
<italic>=</italic>
1 - ε<sub>a</sub>
), thus the contact matrix can be rewritten as<disp-formula id="Equ2"><label>2</label>
<alternatives><mml:math id="M2"><mml:mrow><mml:mfenced close=")" open="("><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mi>C</mml:mi>
<mml:mtext mathvariant="italic">cc</mml:mtext>
</mml:msub>
<mml:mspace width="0.5em"></mml:mspace>
<mml:msub><mml:mi>C</mml:mi>
<mml:mtext mathvariant="italic">ca</mml:mtext>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:mtext mathvariant="italic">ac</mml:mtext>
<mml:mspace width="0.5em"></mml:mspace>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mi>C</mml:mi>
<mml:mtext mathvariant="italic">aa</mml:mtext>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mfenced>
<mml:mo>=</mml:mo>
<mml:msub><mml:mi>q</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mfenced close=")" open="("><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mfrac><mml:mrow><mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:msub><mml:mi>ε</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mi>α</mml:mi>
</mml:mfrac>
<mml:mspace width="1em"></mml:mspace>
<mml:mfrac><mml:mrow><mml:mi>η</mml:mi>
<mml:msub><mml:mi>ε</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mi>α</mml:mi>
</mml:mfrac>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mfrac><mml:msub><mml:mi>ε</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mrow><mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:mi>α</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mspace width="1em"></mml:mspace>
<mml:mfrac><mml:mrow><mml:mi>η</mml:mi>
<mml:mfenced close=")" open="("><mml:mrow><mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:msub><mml:mi>ε</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow><mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:mi>α</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mspace width="0.5em"></mml:mspace>
</mml:mfenced>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<graphic xlink:href="12879_2012_Article_2399_Equ2_HTML.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
</p>
<p>where we have used the relation <italic>N</italic>
<sub><italic>c</italic>
</sub>
= <italic>αN</italic>
and have indicated with <italic>η</italic>
the ratio between the average contact numbers per age class, <italic>η</italic>
= <italic>q</italic>
<sub><italic>a</italic>
</sub>
/<italic>q</italic>
<sub><italic>c</italic>
</sub>
. Interactions are reciprocal such that the number of contacts between children and adults is the same as the number of contacts between adults and children, requiring the matrix to be symmetric, i.e. <italic>C</italic>
<sub><italic>ca</italic>
</sub>
= <italic>C</italic>
<sub><italic>ac</italic>
</sub>
or <italic>ε</italic>
<sub><italic>c</italic>
</sub>
<italic>α</italic>
= <italic>ηε</italic>
<sub><italic>a</italic>
</sub>
(1 − <italic>α</italic>
). We indicate with <italic>ε</italic>
= <italic>ε</italic>
<sub><italic>c</italic>
</sub>
<italic>α</italic>
= <italic>ηε</italic>
<sub><italic>a</italic>
</sub>
(1 − <italic>α</italic>
) the total fraction of contacts across age classes, so that the contact matrix finally reads<disp-formula id="Equ3"><label>3</label>
<alternatives><mml:math id="M3"><mml:mrow><mml:mfenced close=")" open="("><mml:mtable columnalign="center"><mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:msub><mml:mi>C</mml:mi>
<mml:mtext mathvariant="italic">cc</mml:mtext>
</mml:msub>
</mml:mtd>
<mml:mtd columnalign="center"><mml:msub><mml:mi>C</mml:mi>
<mml:mtext mathvariant="italic">ca</mml:mtext>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:msub><mml:mi>C</mml:mi>
<mml:mtext mathvariant="italic">ac</mml:mtext>
</mml:msub>
</mml:mtd>
<mml:mtd columnalign="center"><mml:msub><mml:mi>C</mml:mi>
<mml:mtext mathvariant="italic">aa</mml:mtext>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mfenced>
<mml:mo>=</mml:mo>
<mml:msub><mml:mi>q</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mfenced close=")" open="("><mml:mtable columnalign="center"><mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mfrac><mml:mrow><mml:mi>α</mml:mi>
<mml:mo>−</mml:mo>
<mml:mi>ε</mml:mi>
</mml:mrow>
<mml:msup><mml:mi>α</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mfrac>
</mml:mtd>
<mml:mtd columnalign="center"><mml:mfrac><mml:mi>ε</mml:mi>
<mml:mrow><mml:mi>α</mml:mi>
<mml:mfenced close=")" open="("><mml:mrow><mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:mi>α</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mfrac>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mfrac><mml:mi>ε</mml:mi>
<mml:mrow><mml:mi>α</mml:mi>
<mml:mfenced close=")" open="("><mml:mrow><mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:mi>α</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mfrac>
</mml:mtd>
<mml:mtd columnalign="center"><mml:mfrac><mml:mrow><mml:mi>η</mml:mi>
<mml:mfenced close=")" open="("><mml:mrow><mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:mi>α</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:mo>−</mml:mo>
<mml:mi>ε</mml:mi>
</mml:mrow>
<mml:msup><mml:mfenced close=")" open="("><mml:mrow><mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:mi>α</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mfrac>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mfenced>
</mml:mrow>
</mml:math>
<graphic xlink:href="12879_2012_Article_2399_Equ3_HTML.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
</p>
<p>The matrix is fully expressed in terms of the average number of contacts in the children class, <italic>q</italic>
<sub><italic>c</italic>
</sub>
, the fraction of children <italic>α</italic>
, the ratio <italic>η</italic>
between the average number of contacts in the adult class and the one in the children class, and the parameter <italic>ε</italic>
. The latter is defined between 0 and min{<italic>α</italic>
, <italic>η</italic>
(1 − <italic>α</italic>
)}, and regulates the degree of assortativity between the age classes [<xref ref-type="bibr" rid="CR26">26</xref>
], with <italic>ε</italic>
→ 0 indicating high assortativity and <italic>ε</italic>
→ min{<italic>α</italic>
, <italic>η</italic>
(1 − <italic>α</italic>
)} indicating low assortativity, as schematically shown in Figure <xref rid="Fig2" ref-type="fig">2</xref>
. Table <xref rid="Tab1" ref-type="table">1</xref>
reports the list of variables and parameters used to introduce the age classes. A more general modeling framework that include also other mixing pattern types and that is applicable to different social classifications besides age is the object of future work.<fig id="Fig2"><label>Figure 2</label>
<caption><p><bold>Schematic example of different assortativity levels in mixing patterns.</bold>
Throughout the paper we use <italic>ε</italic>
= <italic>ε</italic>
<sub><italic>c</italic>
</sub>
<italic>α</italic>
as the parameter referring to the assortativity of the mixing pattern, since it represents the total fraction of across-groups contacts. In this scheme we show three examples of different assortativity levels. <bold>A</bold>
: maximum assortativity, corresponding to no mixing between the two classes (<italic>ε</italic>
<sub><italic>a</italic>
</sub>
= <italic>ε</italic>
<sub><italic>c</italic>
</sub>
= 0); <bold>B</bold>
: intermediate assortativity, i.e. a given fraction of the children contacts are directed to adults (like e.g. a random mixing scenario), the others being of the child-child type; <bold>C</bold>
: no assortativity in the children age class, as all contacts established by children are directed to the adults class (<italic>ε</italic>
<sub><italic>c</italic>
</sub>
= 1 and thus <italic>ε</italic>
= <italic>α</italic>
).</p>
</caption>
<graphic xlink:href="12879_2012_Article_2399_Fig2_HTML" id="d29e1348"></graphic>
</fig>
</p>
<table-wrap id="Tab1"><label>Table 1</label>
<caption><p><bold>Variables used to define the age classes in the epidemic model; c = children, a = adults</bold>
</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><th>Variable</th>
<th>Definition</th>
</tr>
</thead>
<tbody><tr><td><italic>α</italic>
</td>
<td>Children fraction of the population</td>
</tr>
<tr><td><italic>q</italic>
<sub><italic>a</italic>
</sub>
, <italic>q</italic>
<sub><italic>c</italic>
</sub>
</td>
<td>Average number of contacts per unit time established by individuals in the children and adult classes, respectively</td>
</tr>
<tr><td><italic>η</italic>
= <italic>q</italic>
<sub><italic>a</italic>
</sub>
/<italic>q</italic>
<sub><italic>c</italic>
</sub>
</td>
<td>Ratio of the average number of contacts</td>
</tr>
<tr><td><italic>ε</italic>
= <italic>ε</italic>
<sub><italic>c</italic>
</sub>
<italic>α</italic>
= <italic>ηε</italic>
<sub><italic>a</italic>
</sub>
(1 − <italic>α</italic>
)</td>
<td>Total fraction of contacts across age classes</td>
</tr>
<tr><td><italic>r</italic>
</td>
<td>Children fraction of the traveling population</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>To compare the theoretical results to realistic situations, we estimate the parameters <italic>ε</italic>
and <italic>η</italic>
from data-driven contact matrices <bold><italic>C</italic>
</bold>
= {<italic>C</italic>
<sub><italic>ij</italic>
</sub>
}, obtained from a population-based prospective survey of mixing patterns in eight European countries using a common paper-diary methodology [<xref ref-type="bibr" rid="CR22">22</xref>
]. We used smoothed daily contact rates based on a bivariate smoother [<xref ref-type="bibr" rid="CR27">27</xref>
] defined for age classes of 1 year interval, relative to all contacts, i.e. including both physical and non-physical contacts [<xref ref-type="bibr" rid="CR22">22</xref>
] (Additional file <xref rid="MOESM1" ref-type="media">1</xref>
). The European countries considered include: Belgium (BE), Germany (DE), Finland (FI), Great Britain (GB), Italy (IT), Luxembourg (LU), The Netherlands (NL), and Poland (PL).</p>
<p>In addition we also informed our model with the mixing patterns for Mexico, obtained from studies on the early outbreak of the 2009 H1N1 pandemic in the country [<xref ref-type="bibr" rid="CR2">2</xref>
]. The values for the parameters {<italic>α</italic>
, <italic>η</italic>
, <italic>ε</italic>
} obtained for all countries under study are reported in Table <xref rid="Tab2" ref-type="table">2</xref>
.<table-wrap id="Tab2"><label>Table 2</label>
<caption><p><bold>Values of the age classes parameters obtained from country-specific statistics and data (see main text for references)</bold>
</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><th>Country</th>
<th><italic>α</italic>
</th>
<th><italic>η</italic>
</th>
<th><bold><italic>ε</italic>
</bold>
</th>
</tr>
</thead>
<tbody><tr><td>Belgium</td>
<td>0.21</td>
<td>1.13</td>
<td>0.125</td>
</tr>
<tr><td>Germany</td>
<td>0.18</td>
<td>0.75</td>
<td>0.098</td>
</tr>
<tr><td>Finland</td>
<td>0.21</td>
<td>0.79</td>
<td>0.091</td>
</tr>
<tr><td>Great Britain</td>
<td>0.22</td>
<td>0.75</td>
<td>0.115</td>
</tr>
<tr><td>Italy</td>
<td>0.17</td>
<td>0.62</td>
<td>0.083</td>
</tr>
<tr><td>Luxembourg</td>
<td>0.22</td>
<td>0.93</td>
<td>0.107</td>
</tr>
<tr><td>The Netherlands</td>
<td>0.22</td>
<td>0.83</td>
<td>0.094</td>
</tr>
<tr><td>Poland</td>
<td>0.21</td>
<td>0.97</td>
<td>0.100</td>
</tr>
<tr><td><italic>Europe (average values)</italic>
</td>
<td><italic>0.20</italic>
</td>
<td><italic>0.79</italic>
</td>
<td><italic>0.097</italic>
</td>
</tr>
<tr><td>Mexico</td>
<td>0.32</td>
<td>0.32</td>
<td>0.063</td>
</tr>
</tbody>
</table>
</table-wrap>
</p>
</sec>
<sec id="Sec5"><title>Spatial metapopulation model with age structure</title>
<p>We consider a population of individuals that is spatially structured into <italic>V</italic>
subpopulations coupled by human mobility patterns, representing a metapopulation network where nodes correspond to subpopulations where the infection dynamics takes place, and links correspond to the mobility processes among them. We initially present this approach considering only one class of individuals, to highlight its main features and the assumptions we make, and in the following we will introduce the age structure of the population. The model is informed with the statistical laws empirically observed in real data on human population and mobility by air-travel, and discussed in the Demographic and travel data subsection. The metapopulation structure is characterized by a random connectivity pattern described by an arbitrary degree distribution <italic>P</italic>
(<italic>k</italic>
). In the following we will explore the role of realistic heterogeneous network structures, adopting power-law degree distributions <italic>P</italic>
(<italic>k</italic>
) ∝ <italic>k</italic>
<sup>−<italic>γ</italic>
</sup>
for analytical convenience, mimicking in this way the airline network as drawn from realistic data. Following the scaling properties observed in real-world mobility data, we define the number of individuals moving from the subpopulation of degree <italic>k</italic>
to the subpopulation of degree <italic>k</italic>
<sup>'</sup>
as <italic>w</italic>
<sub><italic>kk</italic>
'</sub>
= <italic>w</italic>
<sub>0</sub>
(<italic>kk</italic>
<sup>'</sup>
)<sup><italic>θ</italic>
</sup>
. We fix the exponent <italic>θ</italic>
to 0.5 and the scaling factor <italic>w</italic>
<sub>0</sub>
to 1.0, based on the empirical findings [<xref ref-type="bibr" rid="CR15">15</xref>
]. Smaller values of <italic>w</italic>
<sub>0</sub>
are also explored to simulate the implementation of travel-related intervention strategies such as reductions of the travel flows following the start of the outbreak.</p>
<p>We model the travel diffusion process to match the patterns <italic>w</italic>
<sub><italic>kk</italic>
'</sub>
, assuming that travelers are randomly chosen in the population with the per capita diffusion rate <italic>d</italic>
<sub><italic>kk</italic>
'</sub>
= <italic>w</italic>
<sub><italic>kk</italic>
'</sub>
/<italic>N</italic>
<sub><italic>k</italic>
</sub>
where <italic>N</italic>
<sub><italic>k</italic>
</sub>
is a variable indicating the population size [<xref ref-type="bibr" rid="CR28">28</xref>
].</p>
<p>The variables introduced to define the metapopulation model solely depend on the degree <italic>k</italic>
of each subpopulation, therefore we introduce a degree-block notation [<xref ref-type="bibr" rid="CR29">29</xref>
] that assumes statistical equivalence for subpopulations of equal degree [<xref ref-type="bibr" rid="CR30">30</xref>
]. It corresponds to assuming that all subpopulations having the same number of connections are considered statistically identical regarding the features of the metapopulation system relevant for the mobility and disease spreading processes (such as for instance the population size and the traffic of passengers). While disregarding more specific properties of each individual subpopulation – that may be related for instance to local, geographical or cultural aspects – this mean field approximation is able to account for the large degree fluctuations empirically observed, to capture the degree dependence of the system’s properties as found in the data, and also to allow for an analytical treatment of the system’s behavior [<xref ref-type="bibr" rid="CR30">30</xref>
]. The full list of variables used to define the metapopulation model with one class of individuals is provided in Table <xref rid="Tab3" ref-type="table">3</xref>
.<table-wrap id="Tab3"><label>Table 3</label>
<caption><p><bold>Variables defining the metapopulation model with one class of individuals in the degree-block approximation</bold>
</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><th>Variable</th>
<th>Definition</th>
</tr>
</thead>
<tbody><tr><td><italic>k</italic>
</td>
<td>Degree of a subpopulation, i.e. number of airline connections to other subpopulations</td>
</tr>
<tr><td><italic>V</italic>
, <italic>V</italic>
<sub><italic>k</italic>
</sub>
</td>
<td>Total number of subpopulations, number of subpopulations with degree <italic>k</italic>
</td>
</tr>
<tr><td><italic>N</italic>
<sub><italic>k</italic>
</sub>
</td>
<td>Population size of subpopulations with degree <italic>k</italic>
</td>
</tr>
<tr><td><italic>w</italic>
<sub><italic>kk</italic>
'</sub>
= <italic>w</italic>
<sub>0</sub>
(<italic>kk</italic>
<sup>'</sup>
)<sup><italic>θ</italic>
</sup>
</td>
<td>Number of passengers flying from a subpopulation with degree <italic>k</italic>
to a subpopulation with degree <italic>k’</italic>
</td>
</tr>
<tr><td><italic>w</italic>
<sub>0</sub>
</td>
<td>Mobility scale</td>
</tr>
<tr><td><italic>d</italic>
<sub><italic>kk</italic>
'</sub>
= <italic>w</italic>
<sub><italic>kk</italic>
'</sub>
/<italic>N</italic>
<sub><italic>k</italic>
</sub>
</td>
<td>Diffusion rate of passengers flying from a subpopulation with degree to a subpopulation with degree <italic>k’</italic>
</td>
</tr>
</tbody>
</table>
</table-wrap>
</p>
<p>The model defined so far considers a single class of individuals who homogeneously mix in the population. Here we introduce the age structure in order to consider different mixing patterns and travel probabilities depending on the age class, based on the results presented in the previous subsections. We consider a simple SIR compartmental model to describe the infection dynamics of the influenza epidemic, where individuals are assigned to mutually exclusive compartments – susceptible (S), infectious (I), and recovered (R) individuals [<xref ref-type="bibr" rid="CR23">23</xref>
]. Susceptible individuals may contract the infection from infectious individuals and enter the infectious compartment; all infectious individuals then recover permanently and enter the recovered compartment. The disease dynamics is encoded in the next generation matrix <bold><italic>R</italic>
</bold>
= {<italic>R</italic>
<sub><italic>ij</italic>
</sub>
} i.e. the average number of secondary infections in age group <italic>j</italic>
generated by a single primary case in age group <italic>i</italic>
in a fully susceptible population, with <italic>i</italic>
, <italic>j</italic>
∈ {<italic>c</italic>
, <italic>a</italic>
}. We assume that the infectious period is exponentially distributed (with average value <italic>μ</italic>
<sup>− 1</sup>
) and that both transmission rate and recovery rate are independent of the age group, thus neglecting age-specific susceptibility or infectiousness, as in [<xref ref-type="bibr" rid="CR14">14</xref>
]. We consider a fully susceptible population, and also the case of an age-specific prior immunity, as detailed in the next subsection. If we assume that disease transmission may only occur along the contacts captured by the matrix <bold><italic>C</italic>
</bold>
= {<italic>C</italic>
<sub><italic>ij</italic>
</sub>
} [<xref ref-type="bibr" rid="CR24">24</xref>
, <xref ref-type="bibr" rid="CR31">31</xref>
], we can express the generic entry of the next generation matrix as <inline-formula id="IEq1"><alternatives><mml:math id="M4"><mml:mrow><mml:msub><mml:mi>R</mml:mi>
<mml:mtext mathvariant="italic">ia</mml:mtext>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac><mml:mi>β</mml:mi>
<mml:mi>μ</mml:mi>
</mml:mfrac>
<mml:msub><mml:mi>C</mml:mi>
<mml:mtext mathvariant="italic">ia</mml:mtext>
</mml:msub>
<mml:mi>α</mml:mi>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="12879_2012_Article_2399_IEq1_HTML.gif"></inline-graphic>
</alternatives>
</inline-formula>
and <inline-formula id="IEq2"><alternatives><mml:math id="M5"><mml:mrow><mml:msub><mml:mi>R</mml:mi>
<mml:mtext mathvariant="italic">ic</mml:mtext>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac><mml:mi>β</mml:mi>
<mml:mi>μ</mml:mi>
</mml:mfrac>
<mml:msub><mml:mi>C</mml:mi>
<mml:mtext mathvariant="italic">ic</mml:mtext>
</mml:msub>
<mml:mfenced close=")" open="("><mml:mrow><mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:mi>α</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="12879_2012_Article_2399_IEq2_HTML.gif"></inline-graphic>
</alternatives>
</inline-formula>
[<xref ref-type="bibr" rid="CR32">32</xref>
, <xref ref-type="bibr" rid="CR33">33</xref>
] with <italic>i</italic>
∈ {<italic>a</italic>
, <italic>c</italic>
} and <italic>C</italic>
<sub><italic>ij</italic>
</sub>
the contact matrix defined in Eq. (<xref rid="Equ1" ref-type="">1</xref>
). In the application of the contact matrix of Eq. (<xref rid="Equ3" ref-type="">3</xref>
) to the metapopulation model, we assume that all parameters used to define the partition of the population into age classes are independent of spatial features, being therefore constant across subpopulations (e.g. the children population fraction of a subpopulation with degree <italic>k</italic>
is given by <italic>N</italic>
<sub><italic>k</italic>
,<italic>c</italic>
</sub>
= <italic>αN</italic>
<sub><italic>k</italic>
</sub>
). Similarly, we assume that the age-specific travel behavior does not change with the subpopulation of the system. The travel behavior is thus modeled by rescaling the per-capita diffusion rates <italic>d</italic>
<sub><italic>kk</italic>
'</sub>
with the corresponding age-specific probability of travel and the population sizes of the age classes. If we indicate with <italic>r</italic>
the fraction of children traveling (see Table <xref rid="Tab1" ref-type="table">1</xref>
), the per-capita diffusion rate for children traveling from a subpopulation with degree <italic>k</italic>
to a subpopulation with degree <italic>k</italic>
', <inline-formula id="IEq3"><alternatives><mml:math id="M6"><mml:msub><mml:mi>d</mml:mi>
<mml:mrow><mml:mi>k</mml:mi>
<mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mo>,</mml:mo>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
<inline-graphic xlink:href="12879_2012_Article_2399_IEq3_HTML.gif"></inline-graphic>
</alternatives>
</inline-formula>
, is expressed as<disp-formula id="Equ4"><label>4</label>
<alternatives><mml:math id="M7"><mml:mrow><mml:msub><mml:mi>d</mml:mi>
<mml:mrow><mml:mi>k</mml:mi>
<mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mo>,</mml:mo>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mi>r</mml:mi>
<mml:mfrac><mml:mrow><mml:msub><mml:mi>w</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:msup><mml:mfenced close=")" open="("><mml:mrow><mml:mi>k</mml:mi>
<mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
</mml:mrow>
</mml:mfenced>
<mml:mi>θ</mml:mi>
</mml:msup>
</mml:mrow>
<mml:msub><mml:mi>N</mml:mi>
<mml:mrow><mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mfrac>
<mml:mo>=</mml:mo>
<mml:mi>r</mml:mi>
<mml:msub><mml:mi>d</mml:mi>
<mml:mrow><mml:mi>k</mml:mi>
<mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
</mml:mrow>
</mml:msub>
<mml:mfrac><mml:msub><mml:mi>N</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:msub><mml:mi>N</mml:mi>
<mml:mrow><mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mfrac>
<mml:mo>=</mml:mo>
<mml:mfrac><mml:mrow><mml:mi>r</mml:mi>
<mml:msub><mml:mi>d</mml:mi>
<mml:mrow><mml:mi>k</mml:mi>
<mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mi>α</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:math>
<graphic xlink:href="12879_2012_Article_2399_Equ4_HTML.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
</p>
<p>Analogously, the per-capita diffusion rate for adults is given by <inline-formula id="IEq4"><alternatives><mml:math id="M8"><mml:mrow><mml:msub><mml:mi>d</mml:mi>
<mml:mrow><mml:mi>k</mml:mi>
<mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mo>,</mml:mo>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfenced close=")" open="("><mml:mrow><mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:msub><mml:mi>d</mml:mi>
<mml:mrow><mml:mi>k</mml:mi>
<mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
</mml:mrow>
</mml:msub>
<mml:mo stretchy="true">/</mml:mo>
<mml:mfenced close=")" open="("><mml:mrow><mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:mi>α</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="12879_2012_Article_2399_IEq4_HTML.gif"></inline-graphic>
</alternatives>
</inline-formula>
. This allows us to consider different traveling rates depending on the age classes, and to explore a range of values of <italic>r</italic>
, including <italic>r</italic>
= 0, i.e. only adults travel, in addition to its real values obtained from travel statistics.</p>
<p>In the Additional File we also consider a more refined compartmentalization that incorporates a latency period of duration <italic>τ</italic>
<sup>− 1</sup>
to account for the time elapsing from exposure to infectiousness. This corresponds to a more accurate approximation for the description of the disease etiology of influenza. Moreover, it also allows us to extend the applicability of our modeling framework to other infectious diseases where this period may last several days, therefore being non-negligible, as in the case of the severe acute respiratory syndrome (SARS) [<xref ref-type="bibr" rid="CR34">34</xref>
].</p>
</sec>
<sec id="Sec6"><title>2009 H1N1 pandemic case study</title>
<p>To clarify the impact of the study’s findings in a practical situation, we apply our framework to the case study of the 2009 H1N1 pandemic influenza. We parameterize the metapopulation model to the available epidemiological estimates of the outbreak, to the demographic and travel statistics, and to the contact pattern data. We assume a fully susceptible population, and consider values of the reproductive number <italic>R</italic>
<sub>0</sub>
consistent with the estimations obtained for the pandemic through several methods. In particular, we focus on the range from <italic>R</italic>
<sub>0</sub>
= 1.05, corresponding to the lower bound of the estimate obtained from global modeling approaches for the countries in the Northern hemisphere during summer, once seasonal rescaling is taken into account [<xref ref-type="bibr" rid="CR3">3</xref>
], to <italic>R</italic>
<sub>0</sub>
= 1.2 corresponding to the estimates available for Japan [<xref ref-type="bibr" rid="CR35">35</xref>
] and obtained from genetic studies [<xref ref-type="bibr" rid="CR2">2</xref>
], up to <italic>R</italic>
<sub>0</sub>
∈ [1.4, 1.6] as estimated from the early outbreak data of the H1N1 pandemic [<xref ref-type="bibr" rid="CR2">2</xref>
].</p>
<p>The consideration of <italic>R</italic>
<sub>0</sub>
= 1.05 is also important for two additional reasons. First, this value is used here to provide a comparison between the effects that school holidays may have had in the transmission scenario in Europe during Summer 2009 due to the altered contact pattern with respect to school term [<xref ref-type="bibr" rid="CR36">36</xref>
, <xref ref-type="bibr" rid="CR37">37</xref>
]. Indeed, contact data in the UK collected during school term and school holiday periods for 2009 showed that changing mixing patterns resulted in a decrease of approximately 25% to 35% in the reproductive number of influenza during the holidays, considering physical or conversational contacts, respectively [<xref ref-type="bibr" rid="CR36">36</xref>
]. With an estimated <italic>R</italic>
<sub>0</sub>
for the UK in the range [1.4, 1.6] when schools were open, such reductions would thus correspond to approximately <italic>R</italic>
<sub>0</sub>
= 1.05 during school holidays [<xref ref-type="bibr" rid="CR36">36</xref>
]. A second reason for considering the value <italic>R</italic>
<sub>0</sub>
= 1.05, is also for testing scenarios where a higher value of <italic>R</italic>
<sub>0</sub>
may have been reduced following the application of intervention strategies that do not alter the interaction or travel behaviors of individuals, such as e.g. through vaccination or antiviral treatment.</p>
<p>For each value of <italic>R</italic>
<sub>0</sub>
considered, the transmission rate <italic>β</italic>
is calculated from the dominant eigenvalue of the next generation matrix [<xref ref-type="bibr" rid="CR38">38</xref>
], where we set the infectious period <italic>μ</italic>
<sup>− 1</sup>
to 2.5 days [<xref ref-type="bibr" rid="CR3">3</xref>
].</p>
<p>In addition to the value <italic>w</italic>
<sub>0</sub>
= 1, corresponding to the estimate for the mobility scale regulating the travel fluxes obtained from the air mobility network, we also explore <italic>w</italic>
<sub>0</sub>
= 0.5 to simulate the travel-related controls imposed by some countries associated with the self-imposed travel limitations that contributed to a decline of about half the international air traffic to/from Mexico following the international alert in April 2009 (see [<xref ref-type="bibr" rid="CR39">39</xref>
] and references therein).</p>
<p>Results are obtained for eight European countries for which demographic, travel, and contact data are available, and for Mexico, this latter case informed with the age-dependent transmission matrix of Refs. [<xref ref-type="bibr" rid="CR2">2</xref>
, <xref ref-type="bibr" rid="CR14">14</xref>
]. The children age class is defined up to 18 years old for Europe and up to 15 years old for Mexico, to match available data.</p>
<p>We also consider the case of pre-existing immunity in older population and parameterize our model based on the serological evidence indicating that about 30 to 37% of the individuals aged ≥ 60 years had an initial degree of immunity prior to exposure [<xref ref-type="bibr" rid="CR40">40</xref>
–<xref ref-type="bibr" rid="CR42">42</xref>
]. We assume that 33% of individuals aged ≥ 60 years are immune and completely protected against H1N1 pandemic virus. We use the data from the different national age profiles [<xref ref-type="bibr" rid="CR20">20</xref>
] to estimate the corresponding fraction of the adult age class of each country with pre-exposure immunity.</p>
</sec>
</sec>
<sec id="Sec7"><title>Results and discussion</title>
<sec id="Sec8"><title>Calculation of the global invasion threshold</title>
<p>The reproductive number <italic>R</italic>
<sub>0</sub>
provides a threshold condition for a local outbreak in the community; if <italic>R</italic>
<sub>0</sub>
> 1 the epidemic will occur and will affect a finite fraction of the local population, otherwise the disease will die out [<xref ref-type="bibr" rid="CR23">23</xref>
]. The condition for global invasion is however made more complicated by the interplay between the local transmission dynamics and the mobility process that is responsible to seed non-infected subpopulations. Even in the occurrence of a local outbreak, given <italic>R</italic>
<sub>0</sub>
> 1, the epidemic may indeed fail to spread spatially if the mobility rate is not large enough to ensure the travel of infected individuals to other subpopulations before the end of the local outbreak, or if the amount of seeding cases is not large enough to ensure the start of an outbreak in the reached subpopulation due to local extinction events. All these processes have a clear stochastic nature and they are captured by the definition of an additional predictor of the disease dynamics, <italic>R</italic>
<sub>*</sub>
> 1, regulating the number of subpopulations that become infected from a single initially infected subpopulation [<xref ref-type="bibr" rid="CR28">28</xref>
, <xref ref-type="bibr" rid="CR43">43</xref>
–<xref ref-type="bibr" rid="CR46">46</xref>
], analogously to the reproductive number <italic>R</italic>
<sub>0</sub>
at the individual level. An expression for <italic>R</italic>
<sub>*</sub>
has been found in the case of metapopulation epidemic models with different types of mobility processes, including homogeneous, traffic-driven, and population-driven diffusion rates [<xref ref-type="bibr" rid="CR28">28</xref>
, <xref ref-type="bibr" rid="CR46">46</xref>
], commuting-like processes [<xref ref-type="bibr" rid="CR47">47</xref>
, <xref ref-type="bibr" rid="CR48">48</xref>
] and origin–destination processes with adaptive behavior [<xref ref-type="bibr" rid="CR49">49</xref>
] or heterogeneous length of stay at destination [<xref ref-type="bibr" rid="CR50">50</xref>
]. In all those cases the population is assumed to mix homogeneously and to travel according to rates that are uniform across individuals. Here we go beyond those assumptions and examine the relationship between the occurrence of a global outbreak and the age-dependent transmission dynamics coupled with the age-specific travel behavior, through the calculation of the global invasion threshold <italic>R</italic>
<sub>*</sub>
.</p>
<p>Let us consider the invasion process of the epidemic spread at the metapopulation level, by using the subpopulations as our elements of the description of the system. We assume that the outbreak starts in a single initially infected subpopulation of a given degree <italic>k</italic>
and describe the spread from one subpopulation to the neighbor subpopulations through a branching process approximation [<xref ref-type="bibr" rid="CR51">51</xref>
]. We denote by <inline-formula id="IEq5"><alternatives><mml:math id="M9"><mml:msubsup><mml:mi>D</mml:mi>
<mml:mi>k</mml:mi>
<mml:mi>n</mml:mi>
</mml:msubsup>
</mml:math>
<inline-graphic xlink:href="12879_2012_Article_2399_IEq5_HTML.gif"></inline-graphic>
</alternatives>
</inline-formula>
the number of infected subpopulations of degree <italic>k</italic>
at generation <italic>n</italic>
, with <inline-formula id="IEq6"><alternatives><mml:math id="M10"><mml:msubsup><mml:mi>D</mml:mi>
<mml:mi>k</mml:mi>
<mml:mn>0</mml:mn>
</mml:msubsup>
</mml:math>
<inline-graphic xlink:href="12879_2012_Article_2399_IEq6_HTML.gif"></inline-graphic>
</alternatives>
</inline-formula>
being the initially seeded subpopulation, <inline-formula id="IEq7"><alternatives><mml:math id="M11"><mml:msubsup><mml:mi>D</mml:mi>
<mml:mi>k</mml:mi>
<mml:mn>1</mml:mn>
</mml:msubsup>
</mml:math>
<inline-graphic xlink:href="12879_2012_Article_2399_IEq7_HTML.gif"></inline-graphic>
</alternatives>
</inline-formula>
the subpopulations of degree <italic>k</italic>
of generation 1 directly infected by <inline-formula id="IEq8"><alternatives><mml:math id="M12"><mml:msubsup><mml:mi>D</mml:mi>
<mml:mi>k</mml:mi>
<mml:mn>0</mml:mn>
</mml:msubsup>
</mml:math>
<inline-graphic xlink:href="12879_2012_Article_2399_IEq8_HTML.gif"></inline-graphic>
</alternatives>
</inline-formula>
through the mobility process, and so on. By iterating the seeding events, it is possible to describe the evolution of the number <inline-formula id="IEq9"><alternatives><mml:math id="M13"><mml:msubsup><mml:mi>D</mml:mi>
<mml:mi>k</mml:mi>
<mml:mi>n</mml:mi>
</mml:msubsup>
</mml:math>
<inline-graphic xlink:href="12879_2012_Article_2399_IEq9_HTML.gif"></inline-graphic>
</alternatives>
</inline-formula>
of infected subpopulations as follows:<disp-formula id="Equ5"><label>5</label>
<alternatives><mml:math id="M14"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msubsup><mml:mi>D</mml:mi>
<mml:mi>k</mml:mi>
<mml:mi>n</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:munder><mml:mstyle displaystyle="true"><mml:mi mathsize="big">∑</mml:mi>
</mml:mstyle>
<mml:mrow><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:mrow>
</mml:munder>
<mml:msubsup><mml:mi>D</mml:mi>
<mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mrow><mml:mi>n</mml:mi>
<mml:mo>−</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mfenced close=")" open="("><mml:mrow><mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mo>−</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfenced>
<mml:mi>P</mml:mi>
<mml:mfenced close=")" open="("><mml:mrow><mml:mi>k</mml:mi>
<mml:mo stretchy="true">|</mml:mo>
<mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
</mml:mrow>
</mml:mfenced>
<mml:mspace width="1pt"></mml:mspace>
<mml:mfenced close=")" open="("><mml:mrow><mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:munderover><mml:mstyle displaystyle="true"><mml:mi mathsize="big">∑</mml:mi>
</mml:mstyle>
<mml:mrow><mml:mi>m</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mrow><mml:mi>n</mml:mi>
<mml:mo>−</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:munderover>
<mml:mfrac><mml:msubsup><mml:mi>D</mml:mi>
<mml:mi>k</mml:mi>
<mml:mi>m</mml:mi>
</mml:msubsup>
<mml:msub><mml:mi>V</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
<mml:mo>·</mml:mo>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mspace width="8.3em"></mml:mspace>
<mml:msub><mml:mi mathvariant="normal">Ω</mml:mi>
<mml:mrow><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mfenced close=")" open="("><mml:mrow><mml:msub><mml:mi>λ</mml:mi>
<mml:mrow><mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub><mml:mi>λ</mml:mi>
<mml:mrow><mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
<graphic xlink:href="12879_2012_Article_2399_Equ5_HTML.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
</p>
<p>The r.h.s. of equation (<xref rid="Equ5" ref-type="">5</xref>
) describes the contribution of the subpopulations <inline-formula id="IEq10"><alternatives><mml:math id="M15"><mml:msubsup><mml:mi>D</mml:mi>
<mml:mrow><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:mrow>
<mml:mrow><mml:mi>n</mml:mi>
<mml:mo>−</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
</mml:math>
<inline-graphic xlink:href="12879_2012_Article_2399_IEq10_HTML.gif"></inline-graphic>
</alternatives>
</inline-formula>
of degree <italic>k</italic>
' at generation <italic>n</italic>
− 1 to the infection of subpopulations of degree <italic>k</italic>
at generation <italic>n</italic>
. Each of the <inline-formula id="IEq11"><alternatives><mml:math id="M16"><mml:msubsup><mml:mi>D</mml:mi>
<mml:mrow><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:mrow>
<mml:mrow><mml:mi>n</mml:mi>
<mml:mo>−</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
</mml:math>
<inline-graphic xlink:href="12879_2012_Article_2399_IEq11_HTML.gif"></inline-graphic>
</alternatives>
</inline-formula>
subpopulations has <italic>k</italic>
<sup>'</sup>
− 1 possible connections along which the infection can spread. The infection from <inline-formula id="IEq12"><alternatives><mml:math id="M17"><mml:msubsup><mml:mi>D</mml:mi>
<mml:mrow><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:mrow>
<mml:mrow><mml:mi>n</mml:mi>
<mml:mo>−</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
</mml:math>
<inline-graphic xlink:href="12879_2012_Article_2399_IEq12_HTML.gif"></inline-graphic>
</alternatives>
</inline-formula>
to <inline-formula id="IEq13"><alternatives><mml:math id="M18"><mml:msubsup><mml:mi>D</mml:mi>
<mml:mi>k</mml:mi>
<mml:mi>n</mml:mi>
</mml:msubsup>
</mml:math>
<inline-graphic xlink:href="12879_2012_Article_2399_IEq13_HTML.gif"></inline-graphic>
</alternatives>
</inline-formula>
occurs if: (i) the connections departing from nodes with degree <italic>k</italic>
' point to subpopulations with degree <italic>k</italic>
, as ensured by the conditional probability <italic>P</italic>
(<italic>k</italic>
|<italic>k</italic>
<sup>'</sup>
); (ii) the reached subpopulations are not yet infected, as indicated by the probability <inline-formula id="IEq14"><alternatives><mml:math id="M19"><mml:mfenced close=")" open="("><mml:mrow><mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:msubsup><mml:mstyle displaystyle="true"><mml:mi mathsize="big">∑</mml:mi>
</mml:mstyle>
<mml:mrow><mml:mi>m</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mrow><mml:mi>n</mml:mi>
<mml:mo>−</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mfrac><mml:msubsup><mml:mi>D</mml:mi>
<mml:mi>k</mml:mi>
<mml:mi>m</mml:mi>
</mml:msubsup>
<mml:msub><mml:mi>V</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:math>
<inline-graphic xlink:href="12879_2012_Article_2399_IEq14_HTML.gif"></inline-graphic>
</alternatives>
</inline-formula>
, where <italic>V</italic>
<sub><italic>k</italic>
</sub>
is the number of subpopulations with degree <italic>k</italic>
; (iii) the outbreak seeded by <inline-formula id="IEq15"><alternatives><mml:math id="M20"><mml:msub><mml:mi>λ</mml:mi>
<mml:mrow><mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
<inline-graphic xlink:href="12879_2012_Article_2399_IEq15_HTML.gif"></inline-graphic>
</alternatives>
</inline-formula>
and <inline-formula id="IEq16"><alternatives><mml:math id="M21"><mml:msub><mml:mi>λ</mml:mi>
<mml:mrow><mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
<inline-graphic xlink:href="12879_2012_Article_2399_IEq16_HTML.gif"></inline-graphic>
</alternatives>
</inline-formula>
infectious individuals, children and adults, respectively, traveling from subpopulation <italic>k’</italic>
to subpopulation <italic>k</italic>
takes place, and the probability of such event is given by <inline-formula id="IEq17"><alternatives><mml:math id="M22"><mml:mrow><mml:msub><mml:mi mathvariant="normal">Ω</mml:mi>
<mml:mrow><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mfenced close=")" open="("><mml:mrow><mml:msub><mml:mi>λ</mml:mi>
<mml:mrow><mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub><mml:mi>λ</mml:mi>
<mml:mrow><mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="12879_2012_Article_2399_IEq17_HTML.gif"></inline-graphic>
</alternatives>
</inline-formula>
. The latter term is the one that relates the dynamics of the local infection at the individual level to the coarse-grained view that describes the disease invasion at the metapopulation level. It also provides the contribution of children and adults age classes, thus including the effects of non-homogeneous travel behaviors and mixing patterns. The numbers of infectious individuals of each type flying from a subpopulation with degree <italic>k</italic>
' and arriving to a subpopulation with degree <italic>k</italic>
during the entire duration of the outbreak are given by:<disp-formula id="Equ6"><label>6</label>
<alternatives><mml:math id="M23"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mi>λ</mml:mi>
<mml:mrow><mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msub><mml:mi>d</mml:mi>
<mml:mrow><mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>·</mml:mo>
<mml:mfrac><mml:mrow><mml:msub><mml:mi mathvariant="italic">z</mml:mi>
<mml:mrow><mml:mi mathvariant="italic">c</mml:mi>
<mml:mspace width="0.25em"></mml:mspace>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mi mathvariant="italic">N</mml:mi>
<mml:mrow><mml:msup><mml:mi mathvariant="italic">k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="italic">c</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mi>μ</mml:mi>
</mml:mfrac>
<mml:mo>=</mml:mo>
<mml:mfrac><mml:mrow><mml:mi mathvariant="italic">r</mml:mi>
<mml:msub><mml:mi mathvariant="italic">d</mml:mi>
<mml:mrow><mml:mi mathvariant="italic">k</mml:mi>
<mml:msup><mml:mi mathvariant="italic">k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mi>α</mml:mi>
</mml:mfrac>
<mml:mo>·</mml:mo>
<mml:mfrac><mml:mrow><mml:msub><mml:mi mathvariant="italic">z</mml:mi>
<mml:mrow><mml:mi mathvariant="italic">c</mml:mi>
<mml:mspace width="0.25em"></mml:mspace>
</mml:mrow>
</mml:msub>
<mml:mi>α</mml:mi>
<mml:msub><mml:mi mathvariant="italic">N</mml:mi>
<mml:msup><mml:mi mathvariant="italic">k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
</mml:msub>
</mml:mrow>
<mml:mi>μ</mml:mi>
</mml:mfrac>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mspace width="3.4em"></mml:mspace>
<mml:mo>=</mml:mo>
<mml:mi mathvariant="italic">r</mml:mi>
<mml:msub><mml:mi mathvariant="italic">d</mml:mi>
<mml:mrow><mml:mi mathvariant="italic">k</mml:mi>
<mml:msup><mml:mi mathvariant="italic">k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
</mml:mrow>
</mml:msub>
<mml:mfrac><mml:mrow><mml:msub><mml:mi mathvariant="italic">z</mml:mi>
<mml:mrow><mml:mi mathvariant="italic">c</mml:mi>
<mml:mspace width="0.25em"></mml:mspace>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mi mathvariant="italic">N</mml:mi>
<mml:msup><mml:mi mathvariant="italic">k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
</mml:msub>
</mml:mrow>
<mml:mi>μ</mml:mi>
</mml:mfrac>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
<graphic xlink:href="12879_2012_Article_2399_Equ6_HTML.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
<disp-formula id="Equ7"><alternatives><mml:math id="M24"><mml:mrow><mml:msub><mml:mi>λ</mml:mi>
<mml:mrow><mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msub><mml:mi>d</mml:mi>
<mml:mrow><mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>·</mml:mo>
<mml:mfrac><mml:mrow><mml:msub><mml:mi mathvariant="italic">z</mml:mi>
<mml:mrow><mml:mi mathvariant="italic">a</mml:mi>
<mml:mspace width="0.25em"></mml:mspace>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mi mathvariant="italic">N</mml:mi>
<mml:mrow><mml:msup><mml:mi mathvariant="italic">k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="italic">a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mi>μ</mml:mi>
</mml:mfrac>
<mml:mo>=</mml:mo>
<mml:mfenced close=")" open="("><mml:mrow><mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:mi mathvariant="italic">r</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:msub><mml:mi mathvariant="italic">d</mml:mi>
<mml:mrow><mml:mi mathvariant="italic">k</mml:mi>
<mml:msup><mml:mi mathvariant="italic">k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
</mml:mrow>
</mml:msub>
<mml:mfrac><mml:mrow><mml:msub><mml:mi mathvariant="italic">z</mml:mi>
<mml:mrow><mml:mi mathvariant="italic">a</mml:mi>
<mml:mspace width="0.25em"></mml:mspace>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mi mathvariant="italic">N</mml:mi>
<mml:msup><mml:mi mathvariant="italic">k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
</mml:msub>
</mml:mrow>
<mml:mi>μ</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:math>
<graphic xlink:href="12879_2012_Article_2399_Equa_HTML.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
</p>
<p>i.e. the final proportion <italic>z</italic>
<sub><italic>i</italic>
</sub>
of the <inline-formula id="IEq18"><alternatives><mml:math id="M25"><mml:msub><mml:mi>N</mml:mi>
<mml:mrow><mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
<inline-graphic xlink:href="12879_2012_Article_2399_IEq18_HTML.gif"></inline-graphic>
</alternatives>
</inline-formula>
hosts who contract the infection and diffuse with rate <inline-formula id="IEq19"><alternatives><mml:math id="M26"><mml:msub><mml:mi>d</mml:mi>
<mml:mrow><mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
<inline-graphic xlink:href="12879_2012_Article_2399_IEq19_HTML.gif"></inline-graphic>
</alternatives>
</inline-formula>
during their infectious period, with <italic>i</italic>
= <italic>c</italic>
, <italic>a</italic>
. <italic>z</italic>
<sub><italic>c</italic>
</sub>
and <italic>z</italic>
<sub><italic>a</italic>
</sub>
indicate the attack rates in the children and adult age classes, respectively, and they are given by the solution to <inline-formula id="IEq20"><alternatives><mml:math id="M27"><mml:mrow><mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:msub><mml:mi>z</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mtext mathvariant="normal">exp</mml:mtext>
<mml:mspace width="1pt"></mml:mspace>
<mml:mfenced close=")" open="("><mml:mrow><mml:mo>−</mml:mo>
<mml:msub><mml:mstyle displaystyle="true"><mml:mi mathsize="big">∑</mml:mi>
</mml:mstyle>
<mml:mi>j</mml:mi>
</mml:msub>
<mml:msub><mml:mi>R</mml:mi>
<mml:mtext mathvariant="italic">ij</mml:mtext>
</mml:msub>
<mml:msub><mml:mi>z</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="12879_2012_Article_2399_IEq20_HTML.gif"></inline-graphic>
</alternatives>
</inline-formula>
[<xref ref-type="bibr" rid="CR52">52</xref>
]. Figure <xref rid="Fig3" ref-type="fig">3</xref>
A shows the behavior of the final attack rates <italic>z</italic>
<sub><italic>c</italic>
</sub>
and <italic>z</italic>
<sub><italic>a</italic>
</sub>
as a function of the reproductive number <italic>R</italic>
<sub>0</sub>
, considering the partition of the population in children and adults observed in the 8 European countries of the Polymod dataset [<xref ref-type="bibr" rid="CR22">22</xref>
] and their mixing properties. Variations between countries, depending on the age profile and the mixing patterns, are observed. Countries are generally predicted to have significantly larger epidemic sizes in children, as shown, for example, by the case of Italy. This would correspond, on average, to a larger number of individuals in the children class that could potentially seed other subpopulations and sustain the spatial invasion. However this effect is counterbalanced by the age-specific traveling probabilities that are much lower in the children class. In the case of Belgium the final epidemic sizes in the two classes are found to be almost equal, with the size of the epidemic in the adult population being slightly larger than the one in the children population. This is the only country in the dataset under consideration that has a ratio <italic>η</italic>
larger than 1, indicating a larger average number of contacts established by adults with respect to children, likely induced by the specific survey methodology adopted [<xref ref-type="bibr" rid="CR22">22</xref>
].<fig id="Fig3"><label>Figure 3</label>
<caption><p><bold>Final size, extinction probability, and global invasion threshold vs.</bold>
<bold><italic>R</italic>
</bold>
<sub><bold>0</bold>
</sub>
<bold>. A</bold>
-<bold>B</bold>
: Final sizes and extinction probabilities per age class as functions of the reproductive number <italic>R</italic>
<sub>0</sub>
. The various curves for the eight European countries under study are shown by means of a shaded area, with the exception of Belgium, see below. The maximum value for the epidemic size in children (and minimum for the epidemic size in adults) is obtained for Italy; the opposite is obtained for Poland. The situation is reversed for the extinction probabilities – the maximum value for the extinction probability in children (and minimum for the extinction probability in adults) is obtained for Poland; the opposite is obtained for Italy. In both plots, Belgium is a standalone example, with <italic>z</italic>
<sub><italic>c</italic>
</sub>
> <italic>z</italic>
<sub><italic>a</italic>
</sub>
and <italic>π</italic>
<sub><italic>c</italic>
</sub>
> <italic>π</italic>
<sub><italic>a</italic>
</sub>
, differently from all other countries and due to the fact that it is the only population in the dataset to have <italic>η</italic>
> 1, as discussed in the main text. The dashed line represents the case of homogeneous mixing when no partition of the population is considered (in panel B it corresponds to the function <italic>R</italic>
<sub>0</sub>
<sup>− 1</sup>
). <bold>C</bold>
: Global invasion threshold <italic>R</italic>
<sub>*</sub>
as a function of the reproductive number <italic>R</italic>
<sub>0</sub>
, for different values of the parameters describing the mobility process. Air mobility networks having degree distributions <italic>P</italic>
(<italic>k</italic>
) ∝ <italic>k</italic>
<sup>− <italic>γ</italic>
</sup>
with <italic>γ</italic>
= 2 and <italic>γ</italic>
= 3 are shown to consider different levels of heterogeneity. The results obtained in the two cases are compared to the scenarios with homogeneous diffusion rates <italic>d</italic>
<sub><italic>kk</italic>
'</sub>
= <italic>w</italic>
<sub>0</sub>
(<italic>kk</italic>
<sup>'</sup>
)<sup><italic>θ</italic>
</sup>
obtained for <italic>θ</italic>
= 0. All curves are obtained by setting the fraction of passengers in the children class equal to the observed data, i.e. <italic>r</italic>
= 7<italic>%</italic>
, and informing the model with the European average values for <italic>α</italic>
, <italic>η</italic>
, <italic>ε</italic>
.</p>
</caption>
<graphic xlink:href="12879_2012_Article_2399_Fig3_HTML" id="d29e3699"></graphic>
</fig>
</p>
<p>If we indicate with <italic>π</italic>
<sub><italic>c</italic>
</sub>
(<italic>π</italic>
<sub><italic>a</italic>
</sub>
) the probability of extinction given that a single infected individual of class <italic>C</italic>
(<italic>α</italic>
) is introduced in the population, the probability <inline-formula id="IEq21"><alternatives><mml:math id="M28"><mml:mrow><mml:msub><mml:mi mathvariant="normal">Ω</mml:mi>
<mml:mrow><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mfenced close=")" open="("><mml:mrow><mml:msub><mml:mi>λ</mml:mi>
<mml:mrow><mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub><mml:mi>λ</mml:mi>
<mml:mrow><mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
<mml:mspace width="0.25em"></mml:mspace>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="12879_2012_Article_2399_IEq21_HTML.gif"></inline-graphic>
</alternatives>
</inline-formula>
that <inline-formula id="IEq22"><alternatives><mml:math id="M29"><mml:msub><mml:mi>λ</mml:mi>
<mml:mrow><mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
<inline-graphic xlink:href="12879_2012_Article_2399_IEq22_HTML.gif"></inline-graphic>
</alternatives>
</inline-formula>
and <inline-formula id="IEq23"><alternatives><mml:math id="M30"><mml:msub><mml:mi>λ</mml:mi>
<mml:mrow><mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
<inline-graphic xlink:href="12879_2012_Article_2399_IEq23_HTML.gif"></inline-graphic>
</alternatives>
</inline-formula>
infectious individuals traveling from subpopulation <italic>k’</italic>
to subpopulation <italic>k</italic>
would start an outbreak can be expressed as <inline-formula id="IEq24"><alternatives><mml:math id="M31"><mml:mrow><mml:mspace width="0.25em"></mml:mspace>
<mml:msub><mml:mi mathvariant="normal">Ω</mml:mi>
<mml:mrow><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mfenced close=")" open="("><mml:mrow><mml:msub><mml:mi>λ</mml:mi>
<mml:mrow><mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub><mml:mi>λ</mml:mi>
<mml:mrow><mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:msub><mml:mi>π</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:msup><mml:mrow></mml:mrow>
<mml:msub><mml:mi>λ</mml:mi>
<mml:mrow><mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
</mml:msup>
<mml:msub><mml:mi>π</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
<mml:msup><mml:mrow></mml:mrow>
<mml:msub><mml:mi>λ</mml:mi>
<mml:mrow><mml:msup><mml:mi>k</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:msup>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="12879_2012_Article_2399_IEq24_HTML.gif"></inline-graphic>
</alternatives>
</inline-formula>
. Here the two processes are considered as independent since we assume a multi-type branching process approximation. The probabilities of epidemic extinction given the introduction of a single case are determined as the solutions of the quadratic system dependent on the elements <italic>R</italic>
<sub><italic>ij</italic>
</sub>
of the next generation matrix [<xref ref-type="bibr" rid="CR14">14</xref>
], <italic>π</italic>
<sub><italic>i</italic>
</sub>
= [1 + <italic>R</italic>
<sub><italic>ii</italic>
</sub>
(1 − <italic>π</italic>
<sub><italic>i</italic>
</sub>
) + <italic>R</italic>
<sub><italic>ji</italic>
</sub>
(1 − <italic>π</italic>
<sub><italic>j</italic>
</sub>
)]<sup>− 1</sup>
, with <italic>i</italic>
= <italic>c</italic>
, <italic>a</italic>
. If <italic>R</italic>
<sub>0</sub>
< 1, the only solution is <italic>π</italic>
<sub><italic>c</italic>
</sub>
= <italic>π</italic>
<sub><italic>a</italic>
</sub>
= 1, i.e. the epidemic dies off. Otherwise, the system has solutions (<italic>π</italic>
<sub><italic>c</italic>
</sub>
, <italic>π</italic>
<sub><italic>a</italic>
</sub>
) in the range [0,1], as shown in Figure <xref rid="Fig3" ref-type="fig">3</xref>
B. All countries (except Belgium) display a larger probability of extinction related to the introduction of a single adult case in the population, with respect to the introduction of a single infectious child. Given the mixing patterns, children are therefore more likely to start an outbreak than adults. <italic>π</italic>
<sub><italic>a</italic>
</sub>
ranges between 96.5% in the case of Italy and 99% in the case of Poland for <italic>R</italic>
<sub>0</sub>
= 1.05, and between 83.5% and 88% for <italic>R</italic>
<sub>0</sub>
= 1.2, showing that there were small chances for the H1N1 pandemic outbreak to start in the summer in those countries, in agreement with the observed unfolding [<xref ref-type="bibr" rid="CR53">53</xref>
]. In general, a large dependence of the probability of extinction on the reproductive number is observed. Analogously to the behavior discussed for the epidemic size, also for the probability of extinction Belgium represents a special case for which <italic>π</italic>
<sub><italic>c</italic>
</sub>
is slightly larger than <italic>π</italic>
<sub><italic>a</italic>
</sub>
, again induced by the larger average number of contacts in the adult class, differently from all other countries for which the data is available. Finally, in the case of homogeneous mixing in a non-partitioned population, we recover the probability of extinction to be equal to <italic>R</italic>
<sub>0</sub>
<sup>−1</sup>
(dashed line in the figure).</p>
<p>The quantities reported in panels A,B of Figure <xref rid="Fig3" ref-type="fig">3</xref>
represent the ingredients to assess the risk of a global epidemic as driven by the partition of the population into age classes and the non-homogeneous mixing pattern considered, as also discussed in Ref. [<xref ref-type="bibr" rid="CR14">14</xref>
]. The additional role of the non-homogeneous travel behavior is considered explicitly by modeling the invasion process through Eq. (<xref rid="Equ5" ref-type="">5</xref>
). To solve the system analytically, we simplify the recursive relation of Eq. (<xref rid="Equ5" ref-type="">5</xref>
) in the assumption of mild epidemics (i.e. in the limit of <italic>R</italic>
<sub>0</sub>
close to 1), introduced or emerged in the system through a localized seeding event (so that the number of infected subpopulations can be neglected at the early stage of the spatial invasion), and considering the case of a mobility network lacking topological correlations (in this approximation the conditional probability <italic>P</italic>
(<italic>k</italic>
|<italic>k</italic>
<sup>'</sup>
) can be simplified, see the Additional file <xref rid="MOESM1" ref-type="media">1</xref>
). By manipulating the Equation (see the Additional file <xref rid="MOESM1" ref-type="media">1</xref>
for the full details of the calculation), we obtain a condition allowing the increase of the number of infected subpopulations and a global epidemic in the metapopulation system only if<disp-formula id="Equ8"><label>7</label>
<alternatives><mml:math id="M32"><mml:mrow><mml:msub><mml:mi>R</mml:mi>
<mml:mo>*</mml:mo>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfenced close="]" open="["><mml:mrow><mml:mfenced close=")" open="("><mml:mrow><mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:msub><mml:mi>π</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
<mml:mi>r</mml:mi>
<mml:msub><mml:mi>z</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:mfenced close=")" open="("><mml:mrow><mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:msub><mml:mi>π</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
<mml:mfenced close=")" open="("><mml:mrow><mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:msub><mml:mi>z</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
<mml:mfrac><mml:msub><mml:mi>w</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mi>μ</mml:mi>
</mml:mfrac>
<mml:mfrac><mml:mrow><mml:mfenced close="〉" open="〈"><mml:msup><mml:mi>k</mml:mi>
<mml:mrow><mml:mn>2</mml:mn>
<mml:mo>+</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>θ</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mfenced>
<mml:mo>−</mml:mo>
<mml:mfenced close="〉" open="〈"><mml:msup><mml:mi>k</mml:mi>
<mml:mrow><mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>θ</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mfenced>
</mml:mrow>
<mml:mfenced close="〉" open="〈"><mml:mi>k</mml:mi>
</mml:mfenced>
</mml:mfrac>
<mml:mo>></mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
<graphic xlink:href="12879_2012_Article_2399_Equ7_HTML.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
</p>
<p>thus defining the global invasion threshold of the metapopulation system. Equation (<xref rid="Equ7" ref-type="">7</xref>
) defines the threshold condition for the global invasion: if <italic>R</italic>
<sub>*</sub>
assumes values larger than 1, the epidemic starting from a given subpopulation will reach global proportion affecting a finite fraction of the subpopulations of the system; if instead <italic>R</italic>
<sub>*</sub>
< 1, the epidemic will be contained at its source and will not spread further to other locations. The global invasion threshold is a complex function of the disease history parameters, and of the parameters describing the age-specific mixing patterns and travel behavior through <italic>π</italic>
<sub><italic>c</italic>
</sub>
, <italic>π</italic>
<sub><italic>a</italic>
</sub>
, <italic>r</italic>
, <italic>z</italic>
<sub><italic>c</italic>
</sub>
, <italic>z</italic>
<sub><italic>a</italic>
</sub>
. Its dependence on the population spatial structure is embodied by travel fluxes and the topology of the mobility network, through <italic>w</italic>
<sub>0</sub>
, <italic>θ</italic>
, and the degree moments 〈<italic>k</italic>
<sup><italic>a</italic>
</sup>
〉. However it is important to note that this indicator does not depend on the number of subpopulations <italic>V</italic>
of the system, therefore it may be applicable to a variety of countries for which data is available, independently of their size. In the following we explore the dependence of <italic>R</italic>
<sub>*</sub>
on these multiple factors, examine their role in driving the pandemic extinction or invasion, and provide possible applications examples for a set of countries considering the 2009 H1N1 pandemic.</p>
</sec>
<sec id="Sec9"><title>Impact of air mobility</title>
<p>The term <italic>w</italic>
<sub>0</sub>
(<italic>k</italic>
<sup>2 + 2 <italic>θ</italic>
</sup>
− <italic>k</italic>
<sup>1 + 2 <italic>θ</italic>
</sup>
)/<italic>k</italic>
represents the contribution of the air mobility network to the global invasion threshold, with <italic>w</italic>
<sub>0</sub>
and <italic>θ</italic>
regulating the travel fluxes, and the various moments of <italic>k</italic>
, <inline-formula id="IEq25"><alternatives><mml:math id="M33"><mml:mrow><mml:msup><mml:mi>k</mml:mi>
<mml:mi>m</mml:mi>
</mml:msup>
<mml:mo>=</mml:mo>
<mml:msub><mml:mstyle displaystyle="true"><mml:mi mathsize="big">∑</mml:mi>
</mml:mstyle>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:msup><mml:mi>k</mml:mi>
<mml:mi>m</mml:mi>
</mml:msup>
<mml:mi>P</mml:mi>
<mml:mfenced close=")" open="("><mml:mi>k</mml:mi>
</mml:mfenced>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="12879_2012_Article_2399_IEq25_HTML.gif"></inline-graphic>
</alternatives>
</inline-formula>
, expressing the dependence on the network structure as encoded in its degree distribution <italic>P</italic>
(<italic>k</italic>
) = <italic>k</italic>
<sup>− <italic>γ</italic>
</sup>
. The large degree fluctuations found in real transportation systems and mobility patterns [<xref ref-type="bibr" rid="CR15">15</xref>
–<xref ref-type="bibr" rid="CR19">19</xref>
] constitute one of the mechanisms responsible for driving <italic>R</italic>
<sub>*</sub>
to considerable high values, even when small values of the reproductive number are considered. Figure <xref rid="Fig3" ref-type="fig">3</xref>
C shows the dependence of <italic>R</italic>
<sub>*</sub>
on the reproductive number accounting for changes in the topological heterogeneity of the mobility network (<italic>γ</italic>
= 2 and <italic>γ</italic>
= 3) and in the traffic heterogeneity along the air connections (<italic>θ</italic>
= 0.5 and <italic>θ</italic>
= 0, the latter being the homogeneous traffic case). The larger degree fluctuations obtained in the case <italic>γ</italic>
= 2 strongly increase the ratio (<italic>k</italic>
<sup>2 + 2 <italic>θ</italic>
</sup>
− <italic>k</italic>
<sup>1 + 2 <italic>θ</italic>
</sup>
)/<italic>k</italic>
, leading to values of <italic>R</italic>
<sub>*</sub>
ranging from ~10 for <italic>R</italic>
<sub>0</sub>
= 1.05 up to approximately 10<sup>3</sup>
for <italic>R</italic>
<sub>0</sub>
= 2, a value of the reproductive number consistent with the estimate for the 1918–1919 pandemic [<xref ref-type="bibr" rid="CR54">54</xref>
]. It is important to note that, besides the threshold condition <italic>R</italic>
<sub>*</sub>
> 1, the absolute value of the estimator <italic>R</italic>
<sub>*</sub>
provides a quantitative indication of the effective reduction that needs to be reached through public health interventions in order to bring <italic>R</italic>
<sub>*</sub>
below its threshold value, i.e. the difference <italic>R</italic>
<sub>*</sub>
− 1. The <italic>R</italic>
<sub>*</sub>
values for <italic>γ</italic>
= 2 are roughly one order of magnitude greater than the values recovered in the case <italic>γ</italic>
= 3, and in addition the condition for the epidemic invasion is more sensitive to variations in <italic>R</italic>
<sub>0</sub>
in the case of larger heterogeneities in the air mobility patterns (<italic>γ</italic>
= 2 vs. <italic>γ</italic>
= 3).</p>
<p>For both network topologies considered, the epidemic is above the threshold value of 1, and the outbreak is predicted to spread globally in the system for all diseases considered (<italic>R</italic>
<sub>0</sub>
≥ 1.05), consistently with the H1N1 influenza virus invasion at the global level. The partition into classes, though lowering the epidemic sizes and increasing the probability of extinction [<xref ref-type="bibr" rid="CR35">35</xref>
, <xref ref-type="bibr" rid="CR55">55</xref>
, <xref ref-type="bibr" rid="CR56">56</xref>
], is not able to drive the system below the threshold of the global invasion for the range of values explored. In addition, the contribution of the topological heterogeneity of the mobility network (<italic>k</italic>
<sup>2 + 2 <italic>θ</italic>
</sup>
− <italic>k</italic>
<sup>1 + 2 <italic>θ</italic>
</sup>
)/<italic>k</italic>
in increasing the value of <italic>R</italic>
<sub>*</sub>
is so large that it cannot be easily counterbalanced by reductions of the mobility scale <italic>w</italic>
<sub>0</sub>
corresponding to interventions through air travel restrictions. This was already observed in numerical results obtained from data-driven modeling approaches and in analytical predictions based on simple homogeneous mixing among individuals within the subpopulation of the system [<xref ref-type="bibr" rid="CR28">28</xref>
, <xref ref-type="bibr" rid="CR46">46</xref>
, <xref ref-type="bibr" rid="CR50">50</xref>
, <xref ref-type="bibr" rid="CR57">57</xref>
–<xref ref-type="bibr" rid="CR60">60</xref>
]. We will see in the following subsections how differences across countries may impact the conditions for invasion, and will quantitatively assess within this framework the role of travel reductions consistent with the traffic drop observed in the traffic to/from Mexico.</p>
<p>Epidemic containment is instead reached for <italic>R</italic>
<sub>0</sub>
< 1.2 when exploring homogeneous traffic flows (i.e. <italic>θ</italic>
= 0 differently from the realistic scenario <italic>θ</italic>
= 0.5) in the case <italic>γ</italic>
= 3 (see Figure <xref rid="Fig3" ref-type="fig">3</xref>
C), showing how the large variations observed in the traffic flows along air travel connections represent an additional element favoring the spatial spread of the pathogen. On the other hand, extensive measures aimed at radically altering the amount of passenger traveling or the distribution of traffic on the air connections can hardly be achieved in reality.</p>
</sec>
<sec id="Sec10"><title>Impact of the contacts ratio <italic>η</italic>
</title>
<p>The global invasion threshold is a function of the extinction probabilities and of the epidemic sizes per age class for which an explicit solution cannot be obtained in the general case. An approximate solution can be recovered for small <italic>ε</italic>
and in the two limit cases of the contacts ratio <italic>η</italic>
= <italic>q</italic>
<sub><italic>a</italic>
</sub>
/<italic>q</italic>
<sub><italic>c</italic>
</sub>
: <italic>η</italic>
→ 0, i.e. a regime in which almost all contacts are established by children, and <italic>η</italic>
→ 1, i.e. a situation that is almost homogeneous in the distribution of the number of contacts per individual. The approximate solutions are reported in the Additional file <xref rid="MOESM1" ref-type="media">1</xref>
.</p>
<p>Besides these two limit cases, we investigate in Figure <xref rid="Fig4" ref-type="fig">4</xref>
the behavior of <italic>R</italic>
<sub>*</sub>
as a function of the contacts ratio <italic>η</italic>
, exploring the interval [0.25,1] to include the estimates from the Polymod data (range from 0.62 to 0.97, Belgium excluded as discussed before) and to satisfy the existence condition on <italic>ε</italic>
. The value of <italic>α</italic>
is fixed to the European average, <italic>α</italic>
= 0.197, whereas each country is represented with its corresponding assortativity parameter <italic>ε</italic>
, ranging from <italic>ε</italic>
= 0.083 in Italy to <italic>ε</italic>
= 0.125 in Belgium.<fig id="Fig4"><label>Figure 4</label>
<caption><p><bold>Impact of contacts ratio</bold>
<bold><italic>η.</italic>
</bold>
Global invasion threshold <italic>R</italic>
<sub>*</sub>
as a function of the contacts ratio <italic>η</italic>
, for the case of a mobility air network structure having <italic>P</italic>
(<italic>k</italic>
) ∝ <italic>k</italic>
<sup>− <italic>γ</italic>
</sup>
with <italic>γ</italic>
= 2 (panel <bold>A</bold>
) and <italic>γ</italic>
= 3 (panel <bold>B</bold>
). Each plot considers three values of the reproductive number – <italic>R</italic>
<sub>0</sub>
= 1.05, 1.2, 1.4. The various curves for the eight European countries under study are shown by means of a shaded area, for the sake of visualization. The distribution of the population into the children class is set to the average European value for all countries, whereas the assortativity level is left country-specific. Maximum assortativity (and thus minimum <italic>ε</italic>
) is reached for Italy, the opposite observed for Belgium. Here we assume that <italic>r</italic>
= 0. The dashed line indicates the threshold condition <italic>R</italic>
<sub>*</sub>
= 1.</p>
</caption>
<graphic xlink:href="12879_2012_Article_2399_Fig4_HTML" id="d29e4745"></graphic>
</fig>
</p>
<p>The global invasion threshold is predicted to increase with <italic>η</italic>
, showing how a larger number of contacts established by adults, regardless of the across-groups mixing, would favor the spatial propagation. The variations observed among countries are induced by the country-specific assortativity profiles and decrease with <italic>η</italic>
. Therefore, if <italic>R</italic>
<sub>*</sub>
reaches its critical level for relatively small values of <italic>η</italic>
, for which variations across countries are still large, we may reach a heterogeneous outbreak situation in which some countries would experience spatial transmission, while others would be able to contain the outbreak simply due to the role of country-specific age profiles and of the level of assortativity in each population. This may be for example the case with <italic>γ</italic>
= 3 and <italic>R</italic>
<sub>0</sub>
< 1.2 (Figure <xref rid="Fig4" ref-type="fig">4</xref>
B). Results are consistent with the numerical evidence obtained from a data-driven agent-based model in Europe [<xref ref-type="bibr" rid="CR61">61</xref>
].</p>
<p>Countries characterized by particularly low values of <italic>η</italic>
, for cultural, behavioral, and/or social reasons, would be at a lower risk of invasion. Therefore control measures aimed at reducing the contacts ratio <italic>η</italic>
of a specific country may represent an effective policy option to consider. This could be achieved through the application of workplace interventions, including for instance working at home, reducing or avoiding work meetings, and shifting the working timing to reduce overlap at workplaces and at break hours, as well as crowding on transports [<xref ref-type="bibr" rid="CR62">62</xref>
].</p>
<p>Finally, we report on the effects of the inclusion of the latency period in the compartmental model. The results reported in the Additional file <xref rid="MOESM1" ref-type="media">1</xref>
uncover the dependence of the global invasion threshold on the generation time of the disease, i.e. the sum of the latency and infectious periods in the compartmental approximation considered. A simple addition of the latency period to the description of the disease would therefore increase the generation time and, consequently, the global threshold parameter <italic>R</italic>
<sub>*</sub>
, while keeping the qualitative picture unchanged (Additional file <xref rid="MOESM1" ref-type="media">1</xref>
). Individuals would indeed have a longer time span available to travel and potentially spread the disease while carrying an infection.</p>
</sec>
<sec id="Sec11"><title>Impact of assortativity and age profile</title>
<p>We now examine the dependence of the critical condition for the global invasion on the assortativity level of the population partition, by plotting in Figure <xref rid="Fig5" ref-type="fig">5</xref>
<italic>R</italic>
<sub>*</sub>
as a function of the parameters <italic>ε</italic>
and <italic>η</italic>
, where we have set the children fraction <italic>α</italic>
equal to the European average value, <italic>α</italic>
= 0.197. <italic>R</italic>
<sub>*</sub>
is an increasing function of the across-groups mixing <italic>ε</italic>
, indicating that a decrease in the assortativity level of the mixing pattern (corresponding to an increase of <italic>ε</italic>
) favors the spatial invasion. If a larger fraction of contacts is indeed established between adults and children, the local transmission dynamics mainly driven by children is able to spread to a larger fraction of the adult population, thus increasing the chances for the spatial dissemination of the pathogen. Here we study a situation in which <italic>r</italic>
= 0, i.e. only adults travel, in order to isolate the effect of changes in the local transmission while neglecting the role of children in the mobility process.<fig id="Fig5"><label>Figure 5</label>
<caption><p><bold>Impact of assortativity and of school term vs. school holidays.</bold>
Global invasion threshold <italic>R</italic>
<sub>*</sub>
as a function of the across-groups mixing <italic>ε</italic>
and of contacts ratio <italic>η</italic>
, for the reproductive number estimated during school holidays (<italic>R</italic>
<sub>0</sub>
= 1.05, panel <bold>A</bold>
) and school term (<italic>R</italic>
<sub>0</sub>
= 1.4, panel <bold>B</bold>
), based on contact data in the UK [<xref ref-type="bibr" rid="CR62">62</xref>
]. Here we fix the children fraction of the population <italic>α</italic>
to its European average value. The grey area indicates the extinction phase where <italic>R</italic>
<sub>*</sub>
< 1, whereas the colored area refers to the region of the parameter phase space that is above the threshold condition. The rectangular box indicates the area corresponding to the European intervals for the parameters <italic>ε</italic>
and <italic>η</italic>
. The cases for Italy (<italic>ε</italic>
= 0.08, <italic>η</italic>
= 0.62), United Kingdom (<italic>ε</italic>
= 0.11, <italic>η</italic>
= 0.75), and Finland (<italic>ε</italic>
= 0.09, <italic>η</italic>
= 0.79) are highlighted.</p>
</caption>
<graphic xlink:href="12879_2012_Article_2399_Fig5_HTML" id="d29e4928"></graphic>
</fig>
</p>
<p>The rectangle shown in the panels indicates the ranges of the country values for <italic>ε</italic>
and <italic>η</italic>
in Europe. A variation of <italic>η</italic>
from the smallest to the largest of the country values produces a stronger effect on <italic>R</italic>
<sub>*</sub>
with respect to a variation of <italic>ε</italic>
within the European range. As such, <italic>η</italic>
represents an important source of country heterogeneity in the epidemic outcome, as already discussed.</p>
<p>The two panels differ for the value of the reproductive number considered that takes into account the effective reduction in the transmission potential observed during school holidays (panel A, corresponding to <italic>R</italic>
<sub>0</sub>
= 1.05) with respect to school term (panel B, <italic>R</italic>
<sub>0</sub>
= 1.4). These values are estimated for the UK on the basis of contact data for the country in the two periods [<xref ref-type="bibr" rid="CR36">36</xref>
], as illustrated in detail in the Methods section, and here we generalize this result to all countries under study to provide a comparison between school opening and school closure in terms of the predicted risk of a major outbreak. During school holiday period (panel A), European countries are found to be close to the critical level, with a portion of the parameter space for Europe in the extinction region, thus confirming the results presented before regarding variations in the country-specific risks of major epidemics. These findings are compatible with the heterogeneous transmission of H1N1 influenza virus in summer 2009 that was not sustained across continental Europe [<xref ref-type="bibr" rid="CR53">53</xref>
] and identify the main mechanisms responsible of these effects in the interplay between demographic/mixing features and the seeding during school holidays (due to similar school calendars, with the exception of the UK). If instead we consider that the influenza pandemic arrived in the UK during school term, our predictions indicate that the risk of a major epidemic with community and spatial transmission was expected to be very high even during summer period (Figure <xref rid="Fig5" ref-type="fig">5</xref>
B), as observed in the country. Other factors undoubtedly played a role in producing the different transmission across countries, and they include humidity conditions [<xref ref-type="bibr" rid="CR53">53</xref>
], and the timing and magnitude of the coupling of the European continent to Mexico and the United States through international travel [<xref ref-type="bibr" rid="CR3">3</xref>
], both factors being country-specific and not considered here.</p>
<p>Additional countries with similar age profile may be mapped onto the two-dimensional plots of Figure <xref rid="Fig5" ref-type="fig">5</xref>
to gather immediate analytical insight regarding the risk of a major epidemic for the pathogen under consideration, given data availability on the country-specific demographic profiles and mixing habits. This would provide valuable predictions on the conditions for spatio-temporal transmission and informed recommendations for effective control strategies at the start of an outbreak. As an example, we also explored the situation for the United States, considering synthetic information on individual contact networks built from activity surveys and simulations for the city of Portland [<xref ref-type="bibr" rid="CR25">25</xref>
]. By assuming the validity of this synthetic information for the whole country, we can compare the global invasion threshold in the US with the one studied in Europe from real country-specific data. Similar properties are found in the air transportation networks of the two regions, and the slight differences in the age partition (<italic>α</italic>
= 0.24 in the US vs. <italic>α</italic>
= 0.20 in Europe) do not lead to great discrepancies in the behavior of the global threshold as a function of the parameter <italic>η</italic>
(Additional file <xref rid="MOESM1" ref-type="media">1</xref>
).</p>
<p>If we assume that children travel according to the statistics obtained from travel data (i.e. <italic>r</italic>
is set to 7%), the increase in <italic>r</italic>
leads to an increase of <italic>R</italic>
<sub>*</sub>
, as expected, given that a fraction of the individuals driving the local outbreaks also represent potential seeds in new locations not yet affected by the epidemic (see Figure <xref rid="Fig6" ref-type="fig">6</xref>
A). Such a small change in <italic>r</italic>
may also be enough to drive the system from the extinction phase to the major epidemic phase, for certain values of the other parameters. For example, an epidemic starting in a subpopulation in Italy, where the across-groups mixing is equal to <italic>ε</italic>
= 0.08, would reach pandemic proportion if a small fraction of children would travel by air, with respect to the case in which such fraction is neglected (see the vertical line in Figure <xref rid="Fig6" ref-type="fig">6</xref>
A). Given the specific seeding role of children vs. adults, age-targeted entry screening of travelers at airports may be envisaged in an attempt to prevent spatial transmission. The mild and self-limiting nature of most influenza infections, in addition to the presence of asymptomatic infections, is however likely to prevent a perfect identification and notification of cases at entry ports. Given that even a small percentage of children traveling would considerably increase the risk of spatial transmission, leading only to small delays in the invasion process [<xref ref-type="bibr" rid="CR63">63</xref>
], such interventions should be evaluated with respect to their actual efficacy in specific epidemic emergencies and balanced against the resources required for their implementation.<fig id="Fig6"><label>Figure 6</label>
<caption><p><bold>Impact of age-specific travel behavior and age profile. A</bold>
: Global invasion threshold <italic>R</italic>
<sub>*</sub>
as a function of the across-groups mixing <italic>ε</italic>
for the cases of Italy, United Kingdom, and Finland, assuming <italic>R</italic>
<sub>0</sub>
= 1.05. The solid colored lines correspond to the cases when only adults travel (<italic>r</italic>
= 0) and the dashed colored lines to the cases when a percentage of 7% of passengers belongs to the children class. The continuous horizontal line indicates the threshold condition <italic>R</italic>
<sub>*</sub>
= 1. <bold>B</bold>
: Global invasion threshold <italic>R</italic>
<sub>*</sub>
as a function of the across-groups mixing <italic>ε</italic>
: comparison between Europe (<italic>α</italic>
= 0.20, <italic>η</italic>
= 0.79) and Mexico (<italic>α</italic>
= 0.32, <italic>η</italic>
= 0.32). Here we consider <italic>R</italic>
<sub>0</sub>
= 1.05 and <italic>R</italic>
<sub>0</sub>
= 1.4, assuming <italic>r</italic>
= 0.</p>
</caption>
<graphic xlink:href="12879_2012_Article_2399_Fig6_HTML" id="d29e5112"></graphic>
</fig>
</p>
<p>Besides mixing patterns and travel statistics, countries also differ for their age profiles (considered in the model through the parameter <italic>α</italic>
), an easy to access statistics for all countries in the world [<xref ref-type="bibr" rid="CR21">21</xref>
]. The country seed of the 2009 H1N1 pandemic, Mexico, has for instance a larger fraction of the population in the younger age class compared to the population of Europe, <italic>α</italic>
= 0.32 vs. <italic>α</italic>
= 0.20 (where the age classification used for Mexico is up to 15 years old instead of up to 18 years old as adopted in Europe, due to data availability) and is characterized by a larger number of contacts established by children with respect to adults (<italic>η</italic>
= 0.32 vs. <italic>η</italic>
= 0.79). Our predictions indicate that, for the same assortativity values, there exists a range of values of the epidemic transmission potential that is compatible with epidemic containment in Mexico, whereas Europe may experience spatial propagation of the disease (see Figure <xref rid="Fig6" ref-type="fig">6</xref>
B for <italic>R</italic>
<sub>0</sub>
= 1.05). Remarkably, the increase of the fraction of children and the change in the contact ratio, with no change in the across-groups mixing pattern (i.e. same <italic>ε</italic>
), is able to drive the system in the Mexican scenario to the extinction phase, reducing of 94% the global invasion threshold obtained for Europe if we consider a pathogen with a transmission potential compatible with the seasonal estimates of the H1N1 pandemic in Europe during summer 2009 [<xref ref-type="bibr" rid="CR3">3</xref>
]. Such results would be particularly important when considering the country where a new epidemic may emerge, as predictions obtained from previous works considering hypothetical seeding countries would be hardly applicable if demographic features are different. National preparedness plans may be informed by country-specific recommendations through the present approach, and an extensive exploration of the model’s results in terms of classes of demographic and contact pattern features would help in providing more general insights on classes of seeding scenarios.</p>
<p>In the case of <italic>R</italic>
<sub>0</sub>
= 1.4, i.e. the lower bound of the estimate of the reproductive number for Mexico based on the early outbreak of the 2009 H1N1 pandemic [<xref ref-type="bibr" rid="CR2">2</xref>
], our model predict that the country would be above the critical value, thus in agreement with the spatial spread that was observed.</p>
</sec>
<sec id="Sec12"><title>Effect of pre-existing immunity and travel reductions</title>
<p>If we consider pre-existing immunity in the population, calculating the fraction of individuals in the adult age class that corresponds to the estimated values from serological data available after the 2009 H1N1 pandemic [<xref ref-type="bibr" rid="CR40">40</xref>
–<xref ref-type="bibr" rid="CR42">42</xref>
], we find that immunity reduces the condition for the global invasion threshold, as we could expect since a fraction of the population is now modeled to be fully protected against the virus. Figure <xref rid="Fig7" ref-type="fig">7</xref>
addresses the comparison between the two cases, immunity and no-immunity, by showing the two corresponding critical curves <italic>R</italic>
<sub>*</sub>
= 1 in the <italic>α</italic>
, <italic>ε</italic>
plane for Europe and Mexico (panels A and B, respectively). The effect of pre-existing immunity in reducing the probability of a global epidemic spread is shown by the increase in the value of <italic>ε</italic>
corresponding to the invasion condition <italic>R</italic>
<sub>*</sub>
= 1; a larger mixing across age classes is therefore needed for the pathogen to spatially propagate in case a fraction of the older age class is immune, whereas a more assortative population is predicted to be able to contain the emerging epidemic. The effect is very small for the case of Mexico, while for Europe it is more visible given that the considered European population is, on average, older than the Mexican one, and thus a larger fraction of the adult population is assumed to be immune in Europe with respect to Mexico (9.6% in Europe vs. 4.4% in Mexico). The points (<italic>α</italic>
, <italic>ε</italic>
) parameterized with the average European data and with Mexican data (corresponding dots in the Figure <xref rid="Fig7" ref-type="fig">7</xref>
) both lie in the spatial invasion phase when we consider <italic>R</italic>
<sub>0</sub>
≥ 1.2 in Europe and <italic>R</italic>
<sub>0</sub>
≥ 1.4 in Mexico.<fig id="Fig7"><label>Figure 7</label>
<caption><p><bold>Case with immunity.</bold>
Threshold condition <italic>R</italic>
<sub>*</sub>
= 1 as a function of the across-groups mixing <italic>ɛ</italic>
and of the children fraction <italic>α</italic>
for Europe (panel <bold>A</bold>
) and Mexico (panel <bold>B</bold>
): comparison of the no-immunity case with the case of pre-existing immunity and of travel reduction, modeled by setting <italic>w</italic>
<sub>0</sub>
= 0.5, consistently with the empirically observed drop to/from Mexico during the early stage of the 2009 H1N1 pandemic [<xref ref-type="bibr" rid="CR37">37</xref>
]. Here we consider: <italic>R</italic>
<sub>0</sub>
= 1.2 in Europe and <italic>R</italic>
<sub>0</sub>
= 1.4 in Mexico, i.e. the lower bound of the reproductive number estimated for the country from the initial outbreak data [<xref ref-type="bibr" rid="CR2">2</xref>
]. All travelers are adults (<italic>r</italic>
= 0). The three lines, continuous red, dashed red and continuous blue, correspond to pre-existing immunity, no-immunity and travel reduction, respectively. Global epidemic invasion region is above each critical curve. The patterned gray area refers to the region of parameter values that do not satisfy the consistency relation <italic>ɛ</italic>
< min {<italic>α</italic>
, <italic>η</italic>
(1 − <italic>α</italic>
)}.</p>
</caption>
<graphic xlink:href="12879_2012_Article_2399_Fig7_HTML" id="d29e5291"></graphic>
</fig>
</p>
<p>As an additional factor, we also consider the effect resulting from the application of travel reductions. We simulate the travel controls applied by some countries in addition to the self-reaction of the population avoiding travel to the affected area that was observed during the early stage of the 2009 H1N1 pandemic [<xref ref-type="bibr" rid="CR39">39</xref>
], by setting <italic>w</italic>
<sub>0</sub>
= 0.5, i.e. a uniform reduction along all travel connections, independently of the age of travelers. Such reductions reduce the phase space of parameters leading to global invasion, as expected, however they would not be able to lead to a containment of the disease once the model is fitted to the Mexican and European data and to the H1N1 pandemic scenario, confirming previous findings [<xref ref-type="bibr" rid="CR28">28</xref>
, <xref ref-type="bibr" rid="CR39">39</xref>
, <xref ref-type="bibr" rid="CR46">46</xref>
, <xref ref-type="bibr" rid="CR57">57</xref>
–<xref ref-type="bibr" rid="CR60">60</xref>
].</p>
</sec>
<sec id="Sec13"><title>Limitations of the study</title>
<p>Our study is based on a multi-host stochastic metapopulation model that considers several simplifying assumptions that we discuss in this subsection.</p>
<p>The partition of the population into children and adult classes is clearly very schematic, especially if we consider that finer level classifications are available for demographic data at the global scale, and for contacts data, though in a very small set of countries that are the ones considered in this study. Similar considerations arise in the case of additional heterogeneities to be included in the model, as for instance differing patterns of transmission between households, schools, and workplace settings. A higher level of structuring of the population into classes is expected to decrease the epidemic sizes per comparable groups of classes, and to increase the probability of extinction of the epidemic, with respect to our predictions (see for instance the discussions in [<xref ref-type="bibr" rid="CR35">35</xref>
, <xref ref-type="bibr" rid="CR55">55</xref>
, <xref ref-type="bibr" rid="CR56">56</xref>
] and references therein). The inclusion of such features would prevent an analytical treatment of the model and therefore push the study towards an agent-based approach [<xref ref-type="bibr" rid="CR61">61</xref>
, <xref ref-type="bibr" rid="CR64">64</xref>
, <xref ref-type="bibr" rid="CR65">65</xref>
], for which numerical simulations would represent the only available methodology.</p>
<p>Another assumption concerns the exponentially distributed infectious period. More realistic descriptions of the infectious period – including constant, gamma-distributed or data-driven infectious periods – were found to alter the model results by reducing the probability of extinction. Such findings were however obtained in a single population model with homogeneous mixing [<xref ref-type="bibr" rid="CR66">66</xref>
] and in a two-population model coupled by mobility, where, on the other hand, individuals were allowed only one trip, i.e. to change only one time their subpopulation [<xref ref-type="bibr" rid="CR67">67</xref>
]. Such modeling approaches and corresponding results are not applicable to our case, and a systematic understanding of the impact of the infectious period distribution on the probability of extinction in a metapopulation model is still missing.</p>
<p>Our analytical approach also assumes that the importation or the emergence of an infectious disease is highly localized at the beginning of the outbreak, so that it is possible to approximate the spatial spreading process in terms of a branching process evolving in a set of subpopulations not yet affected by the disease. This is similar to the approximations used to calculate the basic reproductive number in a fully susceptible population, and it is required to treat the model analytically. While this assumption can generally be considered as a good approximation to describe the early phase of an outbreak, more complicated seeding events may occur that would require numerical approaches able to explicitly take into account the initial conditions and assess the epidemic risks.</p>
<p>Our model is fit to demographic statistics of a set of countries and it is informed with H1N1 epidemic estimates to provide quantitative information on the risk for the pandemic invasion in such countries. However, in all our predictions we assume the same population structure (<italic>α</italic>
), contacts and mixing profiles (<italic>ε</italic>
, <italic>η</italic>
), and travel behavior (<italic>r</italic>
) across all subpopulations of the metapopulation system, as informed by the data for a given country. We consider this approximation a reasonable one for regions characterized by population features that are quite uniform across space, for instance if we consider the subpopulations within a given country; already at the European level we noticed how small variations in demographic features and mixing patterns may be responsible for diverse outcomes regarding extinction or invasion, and additional layers of heterogeneities need to be considered when variations among populations within the system are larger (e.g. regarding travel behavior per age class, as shown in Figure <xref rid="Fig1" ref-type="fig">1</xref>
B). This could be achieved by considering subpopulation-dependent variables {<italic>α</italic>
<sub><italic>i</italic>
</sub>
, <italic>ε</italic>
<sub><italic>i</italic>
</sub>
, <italic>η</italic>
<sub><italic>i</italic>
</sub>
, <italic>r</italic>
<sub><italic>i</italic>
</sub>
} (with <italic>i</italic>
indicating the subpopulation), however preventing an analytical treatment to solve the system due to the additional complications considered, thus requiring the use of numerical approaches.</p>
<p>Recent work relying on large-scale transmission models has explored the ability of these approaches to predict the timing of spread of the 2009 A/H1N1 influenza pandemic around the world [<xref ref-type="bibr" rid="CR68">68</xref>
]. Differently from these numerical approaches that can describe the geotemporal propagation of the infectious disease in the population, our model does not provide any temporal information on the epidemic unfolding in that all dynamical aspects are synthetically summarized in a branching process leading to the condition for global invasion. On the other hand, given the relevance of demographic profiles, mixing patterns and age-specific travel resulting from the present study, it would be important to further extend large-scale spatial transmission computational models to include such features. While computationally feasible, the main limitation nowadays is represented by the availability of mixing data and travel behavior for a large set of countries, given a global level objective. This further supports the need to have multiple modeling frameworks that can complement each other in providing important information to characterize an emerging epidemic and its associated risks and impact.</p>
<p>Finally, we note that changes in time of population behavior as a response to the ongoing outbreak cannot be dynamically incorporated in the model. These may refer for example to changes in the contact patterns due to self-awareness or changes at the community level due to the implementation of intervention strategies to control the epidemic [<xref ref-type="bibr" rid="CR49">49</xref>
, <xref ref-type="bibr" rid="CR69">69</xref>
–<xref ref-type="bibr" rid="CR73">73</xref>
]. On the other hand, these scenarios can be separately studied with the model, assuming each of these features to be constant in time, in order to assess the effect of such changes on the corresponding risk of a major epidemic. As possible application examples we provided model predictions for different mixing patterns related to school terms and school holidays during H1N1 pandemic, as well as uniform travel restrictions that may result from self-impositions, national guidelines or travel bans. Future studies will focus on other historical epidemics, like e.g. the case of the 2002–2003 SARS outbreak, in order to explore which mechanisms, among the ones included in this approach and related to the applied interventions or to individual’s self-adaptation, hindered a fully global transmission of the disease.</p>
</sec>
</sec>
<sec id="Sec14"><title>Conclusions</title>
<p>The 2009 H1N1 pandemic represents an example of the important role that age classes have on the local and global spread of the disease during its early stage; the local outbreaks being mainly driven by children leading to epidemics in schools, whereas the adults were mainly responsible for the international dissemination by means of air travel. We introduced and solved a multi-host stochastic metapopulation model to quantify these aspects and characterize the conditions of the population partition and heterogeneous travel behavior that lead to the pandemic global invasion. Notwithstanding the high level of assortativity observed in contact patterns data by age, that increase the probability of pandemic extinction, the model explains the spread at the global level observed in the 2009 H1N1 pandemic as induced by the interplay between the heterogeneity of the air mobility network structure, favoring the spread, and the population partition. A major epidemic is always achieved for <italic>R</italic>
<sub>0</sub>
≥ 1.2 even in the case children are assumed not to travel, when the model is parameterized with European countries data and statistics. Results are also consistent with the occurrence of sporadic outbreaks in continental Europe during summer 2009 and widespread transmission in the UK, once the model is informed with the substantial reduction in transmission associated to school holidays.</p>
<p>Despite the presence of various other epidemiological factors that may influence the epidemic outcome, our results suggest that the variations of demographic and mixing profiles across countries are an important source of heterogeneity in the epidemic outcome. This applies in particular to the contacts ratio that is observed to vary significantly and to have a large impact on the invasion potential. Such findings calls for the need to develop further studies in order to identify the social factors that affect this parameter and design targeted interventions, such as work-related measures, that may lower it, thus reducing the risk of an outbreak.</p>
<p>Given the availability of data regarding demographic, mixing and travel profiles, the model results can be used to assess the risks of a given outbreak scenario in a specific country for a newly emerging pathogen. Collecting data on population partitions and mixing matrices or developing alternative methods to estimate the contact patterns based on the demographic information available [<xref ref-type="bibr" rid="CR25">25</xref>
, <xref ref-type="bibr" rid="CR74">74</xref>
] is therefore important to make this approach applicable to a larger range of scenarios, as shown with the example of the US integrating synthetic contact information [<xref ref-type="bibr" rid="CR25">25</xref>
].</p>
<p>Though based on simplifying assumption, the model is able to account for the heterogeneities in the spatial distribution of the population, in the mixing patterns and in the travel behavior, and provide a solution to assess the risk of a major epidemic. We considered a definition for the children age class up to 15 or 18 years old, justified by the available data and statistics, however the approach is transparent to this choice and analogous results to the ones presented can be reached by informing the model with a different definition of classes, as long as statistics informing the groups-specific parameters <italic>α</italic>
, ε, <italic>η</italic>
, and <italic>r</italic>
are available. The approach represents a general framework that can be applicable to other case studies and host population partitions that do not depend on age, such as for instance mixing patterns and travel behaviors depending on socio-economic aspects, or contact profiles and mobility within specific settings where classes correspond to professional roles or conditions of individuals (e.g. health-care workers and patients in hospitals).</p>
</sec>
<sec sec-type="supplementary-material"><title>Electronic supplementary material</title>
<sec id="Sec15"><p><supplementary-material content-type="local-data" id="MOESM1"><media xlink:href="12879_2012_2399_MOESM1_ESM.pdf"><caption><p>Additional file 1: Age-specific contacts and travel patterns in the spatial spread of 2009 H1N1 influenza pandemic.(PDF 5 MB)</p>
</caption>
</media>
</supplementary-material>
</p>
</sec>
</sec>
</body>
<back><app-group><app id="App1"><sec id="Sec16"><title>Authors’ original submitted files for images</title>
<p>Below are the links to the authors’ original submitted files for images.<media position="anchor" xlink:href="12879_2012_2399_MOESM2_ESM.tiff" id="MOESM2"><caption><p>Authors’ original file for figure 1</p>
</caption>
</media>
<media position="anchor" xlink:href="12879_2012_2399_MOESM3_ESM.tiff" id="MOESM3"><caption><p>Authors’ original file for figure 2</p>
</caption>
</media>
<media position="anchor" xlink:href="12879_2012_2399_MOESM4_ESM.tiff" id="MOESM4"><caption><p>Authors’ original file for figure 3</p>
</caption>
</media>
<media position="anchor" xlink:href="12879_2012_2399_MOESM5_ESM.tiff" id="MOESM5"><caption><p>Authors’ original file for figure 4</p>
</caption>
</media>
<media position="anchor" xlink:href="12879_2012_2399_MOESM6_ESM.tiff" id="MOESM6"><caption><p>Authors’ original file for figure 5</p>
</caption>
</media>
<media position="anchor" xlink:href="12879_2012_2399_MOESM7_ESM.tiff" id="MOESM7"><caption><p>Authors’ original file for figure 6</p>
</caption>
</media>
<media position="anchor" xlink:href="12879_2012_2399_MOESM8_ESM.tiff" id="MOESM8"><caption><p>Authors’ original file for figure 7</p>
</caption>
</media>
</p>
</sec>
</app>
</app-group>
<fn-group><fn><p>Andrea Apolloni, Chiara Poletto contributed equally to this work.</p>
</fn>
<fn><p><bold>Competing interests</bold>
</p>
<p>All authors declare that no competing interest exists.</p>
</fn>
<fn><p><bold>Authors’ contributions</bold>
</p>
<p>AA, CP, and VC have all contributed to conceive, design and carry out the study and draft the manuscript. All authors approved the final version of the manuscript.</p>
</fn>
</fn-group>
<ack><p>The authors would like to thank Niel Hens for providing us the daily contact data of the POLYMOD study and for his useful comments; Caterina Rizzo for the age profile of the laboratory confirmed cases in Italy during summer 2009; Stephen Eubank for providing us the synthetic contact data for Portland; and Jose J Ramasco and Pablo Jensen for useful discussions. This work has been partially funded by the ERC Ideas contract no. ERC-2007-Stg204863 (EPIFOR) and the EC-Health contract no. 278433 (PREDEMICS) to VC and CP; the EU-FP7 contract no. 231807 (EPIWORK) to AA; the ANR contract no. ANR-12-MONU-0018 (HARMSFLU) to VC.</p>
</ack>
<ref-list id="Bib1"><title>References</title>
<ref id="CR1"><label>1.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Khan</surname>
<given-names>K</given-names>
</name>
</person-group>
<article-title>Spread of a novel influenza a (H1N1) virus via global airline transportation</article-title>
<source>New Engl J Med</source>
<year>2009</year>
<volume>361</volume>
<issue>2</issue>
<fpage>212</fpage>
<lpage>214</lpage>
<pub-id pub-id-type="doi">10.1056/NEJMc0904559</pub-id>
<pub-id pub-id-type="pmid">19564630</pub-id>
</element-citation>
</ref>
<ref id="CR2"><label>2.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><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>
</person-group>
<article-title>Pandemic potential of a strain of influenza a (H1N1): early findings</article-title>
<source>Science</source>
<year>2009</year>
<volume>324</volume>
<fpage>1557</fpage>
<lpage>1561</lpage>
<pub-id pub-id-type="doi">10.1126/science.1176062</pub-id>
<pub-id pub-id-type="pmid">19433588</pub-id>
</element-citation>
</ref>
<ref id="CR3"><label>3.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Balcan</surname>
<given-names>D</given-names>
</name>
<name><surname>Hu</surname>
<given-names>H</given-names>
</name>
<name><surname>Goncalves</surname>
<given-names>B</given-names>
</name>
<name><surname>Bajardi</surname>
<given-names>P</given-names>
</name>
<name><surname>Poletto</surname>
<given-names>C</given-names>
</name>
<name><surname>Ramasco</surname>
<given-names>JJ</given-names>
</name>
<name><surname>Paolotti</surname>
<given-names>D</given-names>
</name>
<name><surname>Perra</surname>
<given-names>N</given-names>
</name>
<name><surname>Tizzoni</surname>
<given-names>M</given-names>
</name>
<name><surname>Van Den Broeck</surname>
<given-names>W</given-names>
</name>
<name><surname>Colizza</surname>
<given-names>V</given-names>
</name>
<name><surname>Vespignani</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>Seasonal transmission potential and activity peaks of the new influenza A(H1N1): a Monte Carlo likelihood analysis based on human mobility</article-title>
<source>BMC Med</source>
<year>2009</year>
<volume>7</volume>
<fpage>45</fpage>
<pub-id pub-id-type="doi">10.1186/1741-7015-7-45</pub-id>
<pub-id pub-id-type="pmid">19744314</pub-id>
</element-citation>
</ref>
<ref id="CR4"><label>4.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nishiura</surname>
<given-names>H</given-names>
</name>
<name><surname>Castillo-Chavez</surname>
<given-names>C</given-names>
</name>
<name><surname>Safan</surname>
<given-names>M</given-names>
</name>
<name><surname>Chowell</surname>
<given-names>G</given-names>
</name>
</person-group>
<article-title>Transmission potential of the new influenza A (H1N1) virus and its age-specificity in Japan</article-title>
<source>Euro Surveill</source>
<year>2009</year>
<volume>14</volume>
<fpage>19227</fpage>
<pub-id pub-id-type="pmid">19497256</pub-id>
</element-citation>
</ref>
<ref id="CR5"><label>5.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><collab>Novel swine-origin influenza A (H1N1) virus investigation team</collab>
</person-group>
<article-title>Emergence of a novel swine-origin influenza A (H1N1) virus in humans</article-title>
<source>N Engl J Med</source>
<year>2009</year>
<volume>360</volume>
<fpage>2605</fpage>
<lpage>2615</lpage>
<pub-id pub-id-type="doi">10.1056/NEJMoa0903810</pub-id>
<pub-id pub-id-type="pmid">19423869</pub-id>
</element-citation>
</ref>
<ref id="CR6"><label>6.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cutler</surname>
<given-names>J</given-names>
</name>
<name><surname>Schleihauf</surname>
<given-names>E</given-names>
</name>
<name><surname>Hatchette</surname>
<given-names>TF</given-names>
</name>
</person-group>
<article-title>Investigation of the first cases of human-to-human infection with the new swine-origin influenza A (H1N1) virus in Canada</article-title>
<source>CMAJ</source>
<year>2009</year>
<volume>181</volume>
<fpage>159</fpage>
<lpage>163</lpage>
<pub-id pub-id-type="doi">10.1503/cmaj.090859</pub-id>
<pub-id pub-id-type="pmid">19620268</pub-id>
</element-citation>
</ref>
<ref id="CR7"><label>7.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Miller</surname>
<given-names>E</given-names>
</name>
<name><surname>Hoschler</surname>
<given-names>K</given-names>
</name>
<name><surname>Stanford</surname>
<given-names>E</given-names>
</name>
<name><surname>Andrews</surname>
<given-names>N</given-names>
</name>
</person-group>
<article-title>Incidence of 2009 pandemic influenza A H1N1 infection in England: a cross-sectional serological study</article-title>
<source>Lancet</source>
<year>2010</year>
<volume>375</volume>
<fpage>1100</fpage>
<lpage>1108</lpage>
<pub-id pub-id-type="doi">10.1016/S0140-6736(09)62126-7</pub-id>
<pub-id pub-id-type="pmid">20096450</pub-id>
</element-citation>
</ref>
<ref id="CR8"><label>8.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hahné</surname>
<given-names>S</given-names>
</name>
<name><surname>Donker</surname>
<given-names>T</given-names>
</name>
<name><surname>Meijer</surname>
<given-names>A</given-names>
</name>
<name><surname>Timen</surname>
<given-names>A</given-names>
</name>
<name><surname>Van Steenbergen</surname>
<given-names>J</given-names>
</name>
<name><surname>Osterhaus</surname>
<given-names>A</given-names>
</name>
<name><surname>Van Der Sande</surname>
<given-names>M</given-names>
</name>
<name><surname>Koopmans</surname>
<given-names>M</given-names>
</name>
<name><surname>Wallinga</surname>
<given-names>J</given-names>
</name>
<name><surname>Coutinho</surname>
<given-names>R</given-names>
</name>
</person-group>
<article-title>The Dutch New Influenza A(H1N1)v Investigation Team. Epidemiology and control of influenza A(H1N1)v in the Netherlands: the first 115 cases</article-title>
<source>Eurosurveillance</source>
<year>2009</year>
<volume>14</volume>
<issue>27</issue>
<fpage>2</fpage>
<lpage>5</lpage>
</element-citation>
</ref>
<ref id="CR9"><label>9.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><collab>Belgian Working Group on influenza A(H1N1)v</collab>
</person-group>
<article-title>Influenza A(H1N1) virus infections in Belgium, May-June 2009</article-title>
<source>Eurosurveillance</source>
<year>2009</year>
<volume>14</volume>
<issue>28</issue>
<fpage>2</fpage>
<lpage>4</lpage>
</element-citation>
</ref>
<ref id="CR10"><label>10.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><collab>Health Protection Agency and Health Protection Scotland new influenza A(H1N1) investigation teams</collab>
</person-group>
<article-title>Epidemiology of the new influenza A(H1N1) in the United Kingdom, April-May 2009</article-title>
<source>Eurosurveillance</source>
<year>2009</year>
<volume>14</volume>
<issue>19</issue>
<fpage>8</fpage>
<lpage>9</lpage>
<pub-id pub-id-type="pmid">19317973</pub-id>
</element-citation>
</ref>
<ref id="CR11"><label>11.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><collab>Influenza A(H1N1)v investigation teams</collab>
</person-group>
<article-title>Modified surveillance of influenza A(H1N1)v virus infections in France</article-title>
<source>Eurosurveillance</source>
<year>2009</year>
<volume>14</volume>
<issue>29</issue>
<fpage>2</fpage>
<lpage>5</lpage>
</element-citation>
</ref>
<ref id="CR12"><label>12.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rizzo</surname>
<given-names>C</given-names>
</name>
<name><surname>Declich</surname>
<given-names>S</given-names>
</name>
<name><surname>Bella</surname>
<given-names>A</given-names>
</name>
<name><surname>Caporali</surname>
<given-names>MG</given-names>
</name>
<name><surname>Lana</surname>
<given-names>S</given-names>
</name>
<name><surname>Pompa</surname>
<given-names>MG</given-names>
</name>
<name><surname>Vellucci</surname>
<given-names>L</given-names>
</name>
<name><surname>Salmaso</surname>
<given-names>S</given-names>
</name>
</person-group>
<article-title>Enhanced epidemiological surveillance of influenza A(H1N1)v in Italy</article-title>
<source>Eurosurveillance</source>
<year>2009</year>
<volume>14</volume>
<issue>27</issue>
<fpage>6</fpage>
<lpage>9</lpage>
</element-citation>
</ref>
<ref id="CR13"><label>13.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nishiura</surname>
<given-names>H</given-names>
</name>
</person-group>
<article-title>Travel and age of influenza A (H1N1) 2009 virus infection</article-title>
<source>J Travel Med</source>
<year>2010</year>
<volume>17</volume>
<issue>4</issue>
<fpage>269</fpage>
<lpage>270</lpage>
<pub-id pub-id-type="doi">10.1111/j.1708-8305.2010.00418.x</pub-id>
<pub-id pub-id-type="pmid">20636601</pub-id>
</element-citation>
</ref>
<ref id="CR14"><label>14.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nishiura</surname>
<given-names>H</given-names>
</name>
<name><surname>Cook</surname>
<given-names>AR</given-names>
</name>
<name><surname>Cowling</surname>
<given-names>BJ</given-names>
</name>
</person-group>
<article-title>Assortativity and the probability of epidemic extinction: A case study of pandemic influenza A (H1N1-2009)</article-title>
<source>Interdiscipl Perspect Infect Dis</source>
<year>2011</year>
<volume>201</volume>
<issue>1</issue>
<fpage>194507</fpage>
</element-citation>
</ref>
<ref id="CR15"><label>15.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Barrat</surname>
<given-names>A</given-names>
</name>
<name><surname>Barthelemy</surname>
<given-names>M</given-names>
</name>
<name><surname>Pastor-Satorras</surname>
<given-names>R</given-names>
</name>
<name><surname>Vespignani</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>The architecture of complex weighted networks</article-title>
<source>Proc Natl Acad Sci USA</source>
<year>2004</year>
<volume>101</volume>
<fpage>3747</fpage>
<lpage>3752</lpage>
<pub-id pub-id-type="doi">10.1073/pnas.0400087101</pub-id>
<pub-id pub-id-type="pmid">15007165</pub-id>
</element-citation>
</ref>
<ref id="CR16"><label>16.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Guimera</surname>
<given-names>R</given-names>
</name>
<name><surname>Mossa</surname>
<given-names>S</given-names>
</name>
<name><surname>Turtschi</surname>
<given-names>A</given-names>
</name>
<name><surname>Amaral</surname>
<given-names>LAN</given-names>
</name>
</person-group>
<article-title>The worldwide air transportation network: Anomalous centrality, community structure, and cities global role</article-title>
<source>Proc Natl Acad Sci USA</source>
<year>2005</year>
<volume>102</volume>
<fpage>7794</fpage>
<lpage>7799</lpage>
<pub-id pub-id-type="doi">10.1073/pnas.0407994102</pub-id>
<pub-id pub-id-type="pmid">15911778</pub-id>
</element-citation>
</ref>
<ref id="CR17"><label>17.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>DeMontis</surname>
<given-names>A</given-names>
</name>
<name><surname>Barthelemy</surname>
<given-names>M</given-names>
</name>
<name><surname>Chessa</surname>
<given-names>A</given-names>
</name>
<name><surname>Vespignani</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>The structure of interurban traffic: a weighted network analysis</article-title>
<source>Environ Plann B</source>
<year>2007</year>
<volume>34</volume>
<issue>5</issue>
<fpage>905</fpage>
<lpage>924</lpage>
<pub-id pub-id-type="doi">10.1068/b32128</pub-id>
</element-citation>
</ref>
<ref id="CR18"><label>18.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chowell</surname>
<given-names>G</given-names>
</name>
<name><surname>Hyman</surname>
<given-names>JM</given-names>
</name>
<name><surname>Eubank</surname>
<given-names>S</given-names>
</name>
<name><surname>Castillo-Chavez</surname>
<given-names>C</given-names>
</name>
</person-group>
<article-title>Scaling laws for the movement of people between locations in a large city</article-title>
<source>Phys Rev E</source>
<year>2003</year>
<volume>68</volume>
<fpage>066102</fpage>
<pub-id pub-id-type="doi">10.1103/PhysRevE.68.066102</pub-id>
</element-citation>
</ref>
<ref id="CR19"><label>19.</label>
<element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Barrett</surname>
<given-names>CL</given-names>
</name>
</person-group>
<source>TRANSIMS: Transportation Analysis Simulation System</source>
<year>2000</year>
<publisher-loc>Los Alamos National Laboratory</publisher-loc>
<publisher-name>Technical Report LA-UR-00-1725</publisher-name>
</element-citation>
</ref>
<ref id="CR20"><label>20.</label>
<mixed-citation publication-type="other">Eurostat. <ext-link ext-link-type="uri" xlink:href="http://epp.eurostat.ec.europa.eu/portal/page/portal/eurostat/home/">http://epp.eurostat.ec.europa.eu/portal/page/portal/eurostat/home/</ext-link>
,</mixed-citation>
</ref>
<ref id="CR21"><label>21.</label>
<mixed-citation publication-type="other">United Nations: Department of Economics and Social Affairs, Population Division, Population Estimates and Projections Sections. <ext-link ext-link-type="uri" xlink:href="http://esa.un.org/unpd/wpp/Excel-Data/population.htm">http://esa.un.org/unpd/wpp/Excel-Data/population.htm</ext-link>
,</mixed-citation>
</ref>
<ref id="CR22"><label>22.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mossong</surname>
<given-names>J</given-names>
</name>
<name><surname>Hens</surname>
<given-names>N</given-names>
</name>
<name><surname>Jit</surname>
<given-names>M</given-names>
</name>
<name><surname>Beutels</surname>
<given-names>P</given-names>
</name>
<name><surname>Auranen</surname>
<given-names>K</given-names>
</name>
</person-group>
<article-title>Social Contacts and Mixing Patterns Relevant to the Spread of Infectious Diseases</article-title>
<source>PLoS Med</source>
<year>2008</year>
<volume>5</volume>
<issue>3</issue>
<fpage>e74</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pmed.0050074</pub-id>
<pub-id pub-id-type="pmid">18366252</pub-id>
</element-citation>
</ref>
<ref id="CR23"><label>23.</label>
<element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Anderson</surname>
<given-names>RM</given-names>
</name>
<name><surname>May</surname>
<given-names>RM</given-names>
</name>
</person-group>
<source>Infectious Diseases of Humans Dynamics and Control</source>
<year>1992</year>
<publisher-loc>Oxford</publisher-loc>
<publisher-name>Oxford University Press</publisher-name>
</element-citation>
</ref>
<ref id="CR24"><label>24.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wallinga</surname>
<given-names>J</given-names>
</name>
<name><surname>Teunis</surname>
<given-names>P</given-names>
</name>
<name><surname>Kretzschmar</surname>
<given-names>M</given-names>
</name>
</person-group>
<article-title>Using data on social contacts to estimate age-specific transmission parameters for respiratory-spread infectious agents</article-title>
<source>Am J Epidemiol</source>
<year>2006</year>
<volume>164</volume>
<fpage>936</fpage>
<lpage>944</lpage>
<pub-id pub-id-type="doi">10.1093/aje/kwj317</pub-id>
<pub-id pub-id-type="pmid">16968863</pub-id>
</element-citation>
</ref>
<ref id="CR25"><label>25.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Del Valle</surname>
<given-names>SY</given-names>
</name>
<name><surname>Hyman</surname>
<given-names>J</given-names>
</name>
<name><surname>Hethcote</surname>
<given-names>HW</given-names>
</name>
<name><surname>Eubank</surname>
<given-names>SG</given-names>
</name>
</person-group>
<article-title>Mixing patterns between age groups in social networks</article-title>
<source>Soc Network</source>
<year>2007</year>
<volume>29</volume>
<fpage>539</fpage>
<lpage>554</lpage>
<pub-id pub-id-type="doi">10.1016/j.socnet.2007.04.005</pub-id>
</element-citation>
</ref>
<ref id="CR26"><label>26.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kiss</surname>
<given-names>I</given-names>
</name>
<name><surname>Simon</surname>
<given-names>PL</given-names>
</name>
<name><surname>Kao</surname>
<given-names>R</given-names>
</name>
</person-group>
<article-title>A contact network-based formulation of a preferential mixing model</article-title>
<source>Bull Math Biol</source>
<year>2009</year>
<volume>71</volume>
<issue>4</issue>
<fpage>888</fpage>
<lpage>905</lpage>
<pub-id pub-id-type="doi">10.1007/s11538-008-9386-2</pub-id>
<pub-id pub-id-type="pmid">19214639</pub-id>
</element-citation>
</ref>
<ref id="CR27"><label>27.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Goeyvaerts</surname>
<given-names>N</given-names>
</name>
<name><surname>Hens</surname>
<given-names>N</given-names>
</name>
<name><surname>Aerts</surname>
<given-names>M</given-names>
</name>
<name><surname>Beutels</surname>
<given-names>P</given-names>
</name>
</person-group>
<article-title>Model structure analysis to estimate basic immunological processes and maternal risk for parvovirus B19</article-title>
<source>Biostatistics</source>
<year>2011</year>
<volume>12</volume>
<issue>2</issue>
<fpage>283</fpage>
<lpage>302</lpage>
<pub-id pub-id-type="doi">10.1093/biostatistics/kxq059</pub-id>
<pub-id pub-id-type="pmid">20841333</pub-id>
</element-citation>
</ref>
<ref id="CR28"><label>28.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Colizza</surname>
<given-names>V</given-names>
</name>
<name><surname>Vespignani</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>Epidemic modelling in metapopulation systems with heterogeneous coupling pattern: Theory and simulations</article-title>
<source>J Theor Biol</source>
<year>2008</year>
<volume>251</volume>
<fpage>450</fpage>
<lpage>467</lpage>
<pub-id pub-id-type="doi">10.1016/j.jtbi.2007.11.028</pub-id>
<pub-id pub-id-type="pmid">18222487</pub-id>
</element-citation>
</ref>
<ref id="CR29"><label>29.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pastor-Satorras</surname>
<given-names>R</given-names>
</name>
<name><surname>Vespignani</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>Epidemic spreading in scale-free networks</article-title>
<source>Phys Rev Lett</source>
<year>2001</year>
<volume>86</volume>
<fpage>3200</fpage>
<lpage>3203</lpage>
<pub-id pub-id-type="doi">10.1103/PhysRevLett.86.3200</pub-id>
<pub-id pub-id-type="pmid">11290142</pub-id>
</element-citation>
</ref>
<ref id="CR30"><label>30.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Colizza</surname>
<given-names>V</given-names>
</name>
<name><surname>Pastor-Satorras</surname>
<given-names>R</given-names>
</name>
<name><surname>Vespignani</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>Reaction–diffusion processes and metapopulation models in heterogeneous networks</article-title>
<source>Nat Phys</source>
<year>2007</year>
<volume>3</volume>
<fpage>276</fpage>
<lpage>282</lpage>
<pub-id pub-id-type="doi">10.1038/nphys560</pub-id>
</element-citation>
</ref>
<ref id="CR31"><label>31.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Goeyvaerts</surname>
<given-names>N</given-names>
</name>
<name><surname>Hens</surname>
<given-names>N</given-names>
</name>
<name><surname>Ogunjimi</surname>
<given-names>B</given-names>
</name>
<name><surname>Aerts</surname>
<given-names>M</given-names>
</name>
<name><surname>Shkedy</surname>
<given-names>Z</given-names>
</name>
<name><surname>Van Damme</surname>
<given-names>P</given-names>
</name>
<name><surname>Beutels</surname>
<given-names>P</given-names>
</name>
</person-group>
<article-title>Estimating infectious disease parameters from data on social contacts and serological status</article-title>
<source>J R Stat Soc</source>
<year>2010</year>
<volume>59</volume>
<fpage>255</fpage>
<pub-id pub-id-type="doi">10.1111/j.1467-9876.2009.00693.x</pub-id>
</element-citation>
</ref>
<ref id="CR32"><label>32.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Goldstein</surname>
<given-names>E</given-names>
</name>
<name><surname>Apolloni</surname>
<given-names>A</given-names>
</name>
<name><surname>Lewis</surname>
<given-names>B</given-names>
</name>
<name><surname>Miller</surname>
<given-names>JC</given-names>
</name>
<name><surname>Macauley</surname>
<given-names>M</given-names>
</name>
<name><surname>Eubank</surname>
<given-names>S</given-names>
</name>
<name><surname>Lipsitch</surname>
<given-names>MJ</given-names>
</name>
<name><surname>Wallinga</surname>
<given-names>J</given-names>
</name>
</person-group>
<article-title>Distribution of vaccine/antivirals and the “least spread line” in a stratified population</article-title>
<source>J R Soc Interface</source>
<year>2010</year>
<volume>7</volume>
<fpage>755</fpage>
<lpage>764</lpage>
<pub-id pub-id-type="doi">10.1098/rsif.2009.0393</pub-id>
<pub-id pub-id-type="pmid">19828505</pub-id>
</element-citation>
</ref>
<ref id="CR33"><label>33.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Diekmann</surname>
<given-names>O</given-names>
</name>
<name><surname>Heesterbeek</surname>
<given-names>JAP</given-names>
</name>
<name><surname>Roberts</surname>
<given-names>MG</given-names>
</name>
</person-group>
<article-title>The construction of next-generation matrices for compartmental epidemic models</article-title>
<source>J R Soc Interface</source>
<year>2010</year>
<volume>7</volume>
<fpage>873</fpage>
<lpage>885</lpage>
<pub-id pub-id-type="doi">10.1098/rsif.2009.0386</pub-id>
<pub-id pub-id-type="pmid">19892718</pub-id>
</element-citation>
</ref>
<ref id="CR34"><label>34.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Leung</surname>
<given-names>GM</given-names>
</name>
<name><surname>Hedley</surname>
<given-names>AJ</given-names>
</name>
<name><surname>Ho</surname>
<given-names>LM</given-names>
</name>
<name><surname>Chau</surname>
<given-names>P</given-names>
</name>
<name><surname>Wong</surname>
<given-names>IO</given-names>
</name>
<name><surname>Thach</surname>
<given-names>TQ</given-names>
</name>
</person-group>
<article-title>The epidemiology of severe acute respiratory syndrome in the 2003 Hong Kong epidemic: an analysis of all 1755 patients</article-title>
<source>Ann Intern Med</source>
<year>2004</year>
<volume>141</volume>
<fpage>662</fpage>
<lpage>673</lpage>
<pub-id pub-id-type="doi">10.7326/0003-4819-141-9-200411020-00006</pub-id>
<pub-id pub-id-type="pmid">15520422</pub-id>
</element-citation>
</ref>
<ref id="CR35"><label>35.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nishiura</surname>
<given-names>H</given-names>
</name>
<name><surname>Chowell</surname>
<given-names>G</given-names>
</name>
<name><surname>Safan</surname>
<given-names>M</given-names>
</name>
<name><surname>Castillo-Chavez</surname>
<given-names>C</given-names>
</name>
</person-group>
<article-title>Pros and cons of estimating the reproduction number from early epidemic growth rate of influenza A (H1N1) 2009</article-title>
<source>Theor Biol Med Model</source>
<year>2010</year>
<volume>7</volume>
<issue>1</issue>
<fpage>1</fpage>
<pub-id pub-id-type="doi">10.1186/1742-4682-7-1</pub-id>
<pub-id pub-id-type="pmid">20056004</pub-id>
</element-citation>
</ref>
<ref id="CR36"><label>36.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Eames</surname>
<given-names>KTD</given-names>
</name>
<name><surname>Tilston</surname>
<given-names>NL</given-names>
</name>
<name><surname>Brooks-Pollock</surname>
<given-names>E</given-names>
</name>
<name><surname>Edmunds</surname>
<given-names>WJ</given-names>
</name>
</person-group>
<article-title>Measured Dynamic Social Contact Patterns Explain the Spread of H1N1v Influenza</article-title>
<source>PLoS Comput Biol</source>
<year>2012</year>
<volume>8</volume>
<issue>3</issue>
<fpage>e1002425</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pcbi.1002425</pub-id>
<pub-id pub-id-type="pmid">22412366</pub-id>
</element-citation>
</ref>
<ref id="CR37"><label>37.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hens</surname>
<given-names>N</given-names>
</name>
<name><surname>Ayele</surname>
<given-names>GM</given-names>
</name>
<name><surname>Goeyvaerts</surname>
<given-names>N</given-names>
</name>
<name><surname>Aerts</surname>
<given-names>M</given-names>
</name>
<name><surname>Mossong</surname>
<given-names>J</given-names>
</name>
<name><surname>Edmunds</surname>
<given-names>JW</given-names>
</name>
<name><surname>Beutels</surname>
<given-names>P</given-names>
</name>
</person-group>
<article-title>Estimating the impact of school closure on social mixing behaviour and the transmission of close contact infections in eight European countries</article-title>
<source>BMC Infect Dis</source>
<year>2009</year>
<volume>9</volume>
<fpage>187</fpage>
<pub-id pub-id-type="doi">10.1186/1471-2334-9-187</pub-id>
<pub-id pub-id-type="pmid">19943919</pub-id>
</element-citation>
</ref>
<ref id="CR38"><label>38.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Diekmann</surname>
<given-names>O</given-names>
</name>
<name><surname>Heesterbeek</surname>
<given-names>JAP</given-names>
</name>
<name><surname>Metz</surname>
<given-names>JAJ</given-names>
</name>
</person-group>
<article-title>On the definition and the computation of the basic reproduction ratio r0 in models for infectious diseases in heterogeneous populations</article-title>
<source>J Math Biol</source>
<year>1990</year>
<volume>1990</volume>
<issue>28</issue>
<fpage>365</fpage>
<lpage>382</lpage>
</element-citation>
</ref>
<ref id="CR39"><label>39.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bajardi</surname>
<given-names>P</given-names>
</name>
<name><surname>Poletto</surname>
<given-names>C</given-names>
</name>
<name><surname>Ramasco</surname>
<given-names>JJ</given-names>
</name>
<name><surname>Tizzoni</surname>
<given-names>M</given-names>
</name>
<name><surname>Colizza</surname>
<given-names>V</given-names>
</name>
<name><surname>Vespignani</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>Human Mobility Networks, Travel Restrictions, and the Global Spread of 2009 H1N1 Pandemic</article-title>
<source>PLoS One</source>
<year>2011</year>
<volume>6</volume>
<issue>1</issue>
<fpage>e16591</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pone.0016591</pub-id>
<pub-id pub-id-type="pmid">21304943</pub-id>
</element-citation>
</ref>
<ref id="CR40"><label>40.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ikonen</surname>
<given-names>N</given-names>
</name>
<name><surname>Strengell</surname>
<given-names>M</given-names>
</name>
<name><surname>Kinnunen</surname>
<given-names>L</given-names>
</name>
<name><surname>Osterlund</surname>
<given-names>P</given-names>
</name>
<name><surname>Pirhonen</surname>
<given-names>J</given-names>
</name>
<name><surname>Broman</surname>
<given-names>M</given-names>
</name>
<name><surname>Davidkin</surname>
<given-names>I</given-names>
</name>
<name><surname>Ziegler</surname>
<given-names>T</given-names>
</name>
<name><surname>Julkunen</surname>
<given-names>I</given-names>
</name>
</person-group>
<article-title>High frequency of cross-reacting antibodies against 2009 pandemic influenza A(H1N1) virus among the elderly in Finland</article-title>
<source>Euro Surveill</source>
<year>2010</year>
<volume>15</volume>
<fpage>19478</fpage>
<pub-id pub-id-type="pmid">20144443</pub-id>
</element-citation>
</ref>
<ref id="CR41"><label>41.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><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>
<name><surname>Sun</surname>
<given-names>H</given-names>
</name>
<name><surname>Liu</surname>
<given-names>F</given-names>
</name>
<name><surname>Dong</surname>
<given-names>L</given-names>
</name>
<name><surname>DeVos</surname>
<given-names>JR</given-names>
</name>
<name><surname>Gargiullo</surname>
<given-names>PM</given-names>
</name>
<name><surname>Brammer</surname>
<given-names>TL</given-names>
</name>
<name><surname>Cox</surname>
<given-names>NJ</given-names>
</name>
<name><surname>Tumpey</surname>
<given-names>TM</given-names>
</name>
<name><surname>Katz</surname>
<given-names>JM</given-names>
</name>
</person-group>
<article-title>Cross-reactive antibody responses to the 2009 pandemic H1N1 influenza virus</article-title>
<source>N Engl J Med</source>
<year>2009</year>
<volume>361</volume>
<fpage>1945</fpage>
<lpage>1952</lpage>
<pub-id pub-id-type="doi">10.1056/NEJMoa0906453</pub-id>
<pub-id pub-id-type="pmid">19745214</pub-id>
</element-citation>
</ref>
<ref id="CR42"><label>42.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Allwinn</surname>
<given-names>R</given-names>
</name>
<name><surname>Geiler</surname>
<given-names>J</given-names>
</name>
<name><surname>Berger</surname>
<given-names>A</given-names>
</name>
<name><surname>Cinatl</surname>
<given-names>J</given-names>
</name>
<name><surname>Doer</surname>
<given-names>HW</given-names>
</name>
</person-group>
<article-title>Determination of serum antibodies against swine-origin influenza A virus H1N1/09 by immunofluorescence, haemagglutination inhibition, and by neutralization tests: how is the prevalence rate of protecting antibodies in humans?</article-title>
<source>Med Microbiol Immunol</source>
<year>2010</year>
<volume>199</volume>
<fpage>117</fpage>
<lpage>121</lpage>
<pub-id pub-id-type="doi">10.1007/s00430-010-0143-4</pub-id>
<pub-id pub-id-type="pmid">20162304</pub-id>
</element-citation>
</ref>
<ref id="CR43"><label>43.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ball</surname>
<given-names>F</given-names>
</name>
<name><surname>Mollison</surname>
<given-names>D</given-names>
</name>
<name><surname>Scalia-Tomba</surname>
<given-names>G</given-names>
</name>
</person-group>
<article-title>Epidemics with two levels of mixing</article-title>
<source>Ann Appl Probab</source>
<year>1997</year>
<volume>7</volume>
<fpage>46</fpage>
<lpage>89</lpage>
<pub-id pub-id-type="doi">10.1214/aoap/1034625252</pub-id>
</element-citation>
</ref>
<ref id="CR44"><label>44.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cross</surname>
<given-names>P</given-names>
</name>
<name><surname>Lloyd-Smith</surname>
<given-names>JO</given-names>
</name>
<name><surname>Johnson</surname>
<given-names>PLF</given-names>
</name>
<name><surname>Wayne</surname>
<given-names>MG</given-names>
</name>
</person-group>
<article-title>Duelling timescales of host movement and disease recovery determine invasion of disease in structured populations</article-title>
<source>Ecol Lett</source>
<year>2005</year>
<volume>8</volume>
<fpage>587</fpage>
<lpage>595</lpage>
<pub-id pub-id-type="doi">10.1111/j.1461-0248.2005.00760.x</pub-id>
</element-citation>
</ref>
<ref id="CR45"><label>45.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cross</surname>
<given-names>P</given-names>
</name>
<name><surname>Johnson</surname>
<given-names>PLF</given-names>
</name>
<name><surname>Lloyd-Smith</surname>
<given-names>JO</given-names>
</name>
<name><surname>Wayne</surname>
<given-names>MG</given-names>
</name>
</person-group>
<article-title>Utility of R0 as a predictor of disease invasion in structured populations</article-title>
<source>J R Soc Interface</source>
<year>2007</year>
<volume>4</volume>
<fpage>315</fpage>
<lpage>324</lpage>
<pub-id pub-id-type="doi">10.1098/rsif.2006.0185</pub-id>
<pub-id pub-id-type="pmid">17251146</pub-id>
</element-citation>
</ref>
<ref id="CR46"><label>46.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Colizza</surname>
<given-names>V</given-names>
</name>
<name><surname>Vespignani</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>Invasion threshold in heterogeneous metapopulation networks</article-title>
<source>Phys Rev Lett</source>
<year>2007</year>
<volume>99</volume>
<fpage>148701</fpage>
<pub-id pub-id-type="doi">10.1103/PhysRevLett.99.148701</pub-id>
<pub-id pub-id-type="pmid">17930732</pub-id>
</element-citation>
</ref>
<ref id="CR47"><label>47.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Balcan</surname>
<given-names>D</given-names>
</name>
<name><surname>Vespignani</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>Phase transitions in contagion processes mediated by recurrent mobility patterns</article-title>
<source>Nat Phys</source>
<year>2011</year>
<volume>7</volume>
<fpage>581</fpage>
<lpage>586</lpage>
<pub-id pub-id-type="doi">10.1038/nphys1944</pub-id>
<pub-id pub-id-type="pmid">21799702</pub-id>
</element-citation>
</ref>
<ref id="CR48"><label>48.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Belik</surname>
<given-names>V</given-names>
</name>
<name><surname>Geisel</surname>
<given-names>T</given-names>
</name>
<name><surname>Brockmann</surname>
<given-names>D</given-names>
</name>
</person-group>
<article-title>Natural human mobility patterns and spatial spread of infectious diseases</article-title>
<source>Phys Rev X</source>
<year>2011</year>
<volume>1</volume>
<fpage>011001</fpage>
</element-citation>
</ref>
<ref id="CR49"><label>49.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Meloni</surname>
<given-names>S</given-names>
</name>
<name><surname>Perra</surname>
<given-names>N</given-names>
</name>
<name><surname>Arenas</surname>
<given-names>A</given-names>
</name>
<name><surname>Gomes</surname>
<given-names>S</given-names>
</name>
<name><surname>Moreno</surname>
<given-names>Y</given-names>
</name>
<name><surname>Vespignani</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>Modeling human mobility responses to the large-scale spreading of infectious diseases</article-title>
<source>Sci Rep</source>
<year>2011</year>
<volume>1</volume>
<fpage>62</fpage>
<pub-id pub-id-type="doi">10.1038/srep00062</pub-id>
<pub-id pub-id-type="pmid">22355581</pub-id>
</element-citation>
</ref>
<ref id="CR50"><label>50.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Poletto</surname>
<given-names>C</given-names>
</name>
<name><surname>Tizzoni</surname>
<given-names>M</given-names>
</name>
<name><surname>Colizza</surname>
<given-names>V</given-names>
</name>
</person-group>
<article-title>Heterogeneous length of stay of hosts’ movements and spatial epidemic spread</article-title>
<source>Sci Rep</source>
<year>2012</year>
<volume>2</volume>
<fpage>476</fpage>
<pub-id pub-id-type="doi">10.1038/srep00476</pub-id>
<pub-id pub-id-type="pmid">22741060</pub-id>
</element-citation>
</ref>
<ref id="CR51"><label>51.</label>
<element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Harris</surname>
<given-names>TE</given-names>
</name>
</person-group>
<source>The Theory of Branching Processes</source>
<year>1989</year>
<publisher-loc>Mineola</publisher-loc>
<publisher-name>Dover Publications</publisher-name>
</element-citation>
</ref>
<ref id="CR52"><label>52.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ball</surname>
<given-names>F</given-names>
</name>
<name><surname>Clancy</surname>
<given-names>D</given-names>
</name>
</person-group>
<article-title>The final size and severity of a generalised stochastic multitype epidemic model</article-title>
<source>Adv Appl Probab</source>
<year>1993</year>
<volume>25</volume>
<issue>4</issue>
<fpage>721</fpage>
<lpage>736</lpage>
<pub-id pub-id-type="doi">10.2307/1427788</pub-id>
</element-citation>
</ref>
<ref id="CR53"><label>53.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Flasche</surname>
<given-names>S</given-names>
</name>
<name><surname>Hens</surname>
<given-names>N</given-names>
</name>
<name><surname>Boëlle</surname>
<given-names>PY</given-names>
</name>
<name><surname>Mossong</surname>
<given-names>J</given-names>
</name>
<name><surname>van Ballegooijen</surname>
<given-names>WM</given-names>
</name>
<name><surname>Nunes</surname>
<given-names>B</given-names>
</name>
<name><surname>Rizzo</surname>
<given-names>C</given-names>
</name>
<name><surname>Popovici</surname>
<given-names>F</given-names>
</name>
<name><surname>Santa-Olalla</surname>
<given-names>P</given-names>
</name>
<name><surname>Hrubá</surname>
<given-names>F</given-names>
</name>
<name><surname>Parmakova</surname>
<given-names>K</given-names>
</name>
<name><surname>Baguelin</surname>
<given-names>M</given-names>
</name>
<name><surname>van Hoek</surname>
<given-names>AJ</given-names>
</name>
<name><surname>Desenclos</surname>
<given-names>JC</given-names>
</name>
<name><surname>Bernillon</surname>
<given-names>P</given-names>
</name>
<name><surname>Cámara</surname>
<given-names>AL</given-names>
</name>
<name><surname>Wallinga</surname>
<given-names>J</given-names>
</name>
<name><surname>Asikainen</surname>
<given-names>T</given-names>
</name>
<name><surname>White</surname>
<given-names>PJ</given-names>
</name>
<name><surname>Edmunds</surname>
<given-names>WJ</given-names>
</name>
</person-group>
<article-title>Different transmission patterns in the early stages of the influenza A(H1N1)v pandemic: a comparative analysis of 12 European countries</article-title>
<source>Epidemics</source>
<year>2011</year>
<volume>3</volume>
<issue>2</issue>
<fpage>125</fpage>
<lpage>133</lpage>
<pub-id pub-id-type="doi">10.1016/j.epidem.2011.03.005</pub-id>
<pub-id pub-id-type="pmid">21624784</pub-id>
</element-citation>
</ref>
<ref id="CR54"><label>54.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mills</surname>
<given-names>CE</given-names>
</name>
<name><surname>Robins</surname>
<given-names>JM</given-names>
</name>
<name><surname>Lipsitch</surname>
<given-names>M</given-names>
</name>
</person-group>
<article-title>Transmissibility of the 1918 pandemic influenza</article-title>
<source>Nature</source>
<year>2004</year>
<volume>432</volume>
<fpage>904</fpage>
<lpage>906</lpage>
<pub-id pub-id-type="doi">10.1038/nature03063</pub-id>
<pub-id pub-id-type="pmid">15602562</pub-id>
</element-citation>
</ref>
<ref id="CR55"><label>55.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chao</surname>
<given-names>DL</given-names>
</name>
<name><surname>Halloran</surname>
<given-names>ME</given-names>
</name>
<name><surname>Obenchain</surname>
<given-names>VJ</given-names>
</name>
<name><surname>Longini</surname>
<given-names>IM</given-names>
</name>
</person-group>
<article-title>FluTE a publicly available stochastic influenza epidemic simulation model</article-title>
<source>PLoS Comput Biol</source>
<year>2010</year>
<volume>6</volume>
<fpage>e1000656</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pcbi.1000656</pub-id>
<pub-id pub-id-type="pmid">20126529</pub-id>
</element-citation>
</ref>
<ref id="CR56"><label>56.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Watts</surname>
<given-names>DJ</given-names>
</name>
<name><surname>Muhamad</surname>
<given-names>R</given-names>
</name>
<name><surname>Medina</surname>
<given-names>DC</given-names>
</name>
<name><surname>Dodds</surname>
<given-names>PS</given-names>
</name>
</person-group>
<article-title>Multiscale, resurgent epidemics in a hierarchical metapopulation model</article-title>
<source>Proc Natl Acad Sci USA</source>
<year>2005</year>
<volume>102</volume>
<fpage>11157</fpage>
<lpage>11162</lpage>
<pub-id pub-id-type="doi">10.1073/pnas.0501226102</pub-id>
<pub-id pub-id-type="pmid">16055564</pub-id>
</element-citation>
</ref>
<ref id="CR57"><label>57.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Epstein</surname>
<given-names>JM</given-names>
</name>
<name><surname>Goedecke</surname>
<given-names>DM</given-names>
</name>
<name><surname>Yu</surname>
<given-names>F</given-names>
</name>
<name><surname>Morris</surname>
<given-names>RJ</given-names>
</name>
<name><surname>Wagener</surname>
<given-names>DK</given-names>
</name>
<name><surname>Bobashev</surname>
<given-names>GV</given-names>
</name>
</person-group>
<article-title>Controlling Pandemic Flu: The Value of International Air Travel Restrictions</article-title>
<source>PLoS One</source>
<year>2007</year>
<volume>2</volume>
<fpage>1</fpage>
<lpage>11</lpage>
<pub-id pub-id-type="doi">10.1371/journal.pone.0000401</pub-id>
</element-citation>
</ref>
<ref id="CR58"><label>58.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cooper</surname>
<given-names>BS</given-names>
</name>
<name><surname>Pitman</surname>
<given-names>RJ</given-names>
</name>
<name><surname>Edmunds</surname>
<given-names>WJ</given-names>
</name>
<name><surname>Gay</surname>
<given-names>NJ</given-names>
</name>
</person-group>
<article-title>Delaying the international spread of pandemic influenza</article-title>
<source>PLoS Med</source>
<year>2006</year>
<volume>3</volume>
<fpage>845</fpage>
<lpage>855</lpage>
<pub-id pub-id-type="doi">10.1371/journal.pmed.0030212</pub-id>
</element-citation>
</ref>
<ref id="CR59"><label>59.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hollingsworth</surname>
<given-names>TD</given-names>
</name>
<name><surname>Ferguson</surname>
<given-names>NM</given-names>
</name>
<name><surname>Anderson</surname>
<given-names>RM</given-names>
</name>
</person-group>
<article-title>Will travel restrictions control the international spread of pandemic influenza?</article-title>
<source>Nat Med</source>
<year>2006</year>
<volume>12</volume>
<fpage>497</fpage>
<lpage>499</lpage>
<pub-id pub-id-type="doi">10.1038/nm0506-497</pub-id>
<pub-id pub-id-type="pmid">16675989</pub-id>
</element-citation>
</ref>
<ref id="CR60"><label>60.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lam</surname>
<given-names>EHY</given-names>
</name>
<name><surname>Cowling</surname>
<given-names>BJ</given-names>
</name>
<name><surname>Cook</surname>
<given-names>AR</given-names>
</name>
<name><surname>Wong</surname>
<given-names>JYT</given-names>
</name>
<name><surname>Lau</surname>
<given-names>MSY</given-names>
</name>
<name><surname>Nishiura</surname>
<given-names>H</given-names>
</name>
</person-group>
<article-title>The feasibility of age-specific travel restrictions during influenza pandemics</article-title>
<source>Theor Biol Med Model</source>
<year>2011</year>
<volume>8</volume>
<fpage>44</fpage>
<lpage>57</lpage>
<pub-id pub-id-type="doi">10.1186/1742-4682-8-44</pub-id>
<pub-id pub-id-type="pmid">22078655</pub-id>
</element-citation>
</ref>
<ref id="CR61"><label>61.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Merler</surname>
<given-names>S</given-names>
</name>
<name><surname>Ajelli</surname>
<given-names>M</given-names>
</name>
</person-group>
<article-title>The role of population heterogeneity and human mobility in the spread of pandemic influenza</article-title>
<source>Proc Roy Soc Biol Sci</source>
<year>2010</year>
<volume>277</volume>
<fpage>557</fpage>
<lpage>565</lpage>
<pub-id pub-id-type="doi">10.1098/rspb.2009.1605</pub-id>
</element-citation>
</ref>
<ref id="CR62"><label>62.</label>
<element-citation publication-type="book"><person-group person-group-type="author"><collab>European Centre for Disease Prevention and Control (ECDC)</collab>
</person-group>
<source>Guide to public health measures to reduce the impact of influenza epidemics in Europe</source>
<year>2009</year>
<publisher-loc>Stockholm</publisher-loc>
<publisher-name>The ECDC Menu’ ECDC Technical report</publisher-name>
</element-citation>
</ref>
<ref id="CR63"><label>63.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cowling</surname>
<given-names>BJ</given-names>
</name>
<name><surname>Lau</surname>
<given-names>LLH</given-names>
</name>
<name><surname>Wu</surname>
<given-names>P</given-names>
</name>
<name><surname>Wong</surname>
<given-names>HWC</given-names>
</name>
<name><surname>Fang</surname>
<given-names>VJ</given-names>
</name>
<name><surname>Riley</surname>
<given-names>S</given-names>
</name>
<name><surname>Nishiura</surname>
<given-names>H</given-names>
</name>
</person-group>
<article-title>Entry screening to delay local transmission of 2009 pandemic influenza A(H1N1)</article-title>
<source>BMC Infect Dis</source>
<year>2010</year>
<volume>10</volume>
<fpage>82</fpage>
<pub-id pub-id-type="doi">10.1186/1471-2334-10-82</pub-id>
<pub-id pub-id-type="pmid">20353566</pub-id>
</element-citation>
</ref>
<ref id="CR64"><label>64.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Eubank</surname>
<given-names>S</given-names>
</name>
<name><surname>Guclu</surname>
<given-names>H</given-names>
</name>
<name><surname>Anil Kumar</surname>
<given-names>VS</given-names>
</name>
<name><surname>Marathe</surname>
<given-names>MV</given-names>
</name>
<name><surname>Srinivasan</surname>
<given-names>A</given-names>
</name>
<name><surname>Toroczkai</surname>
<given-names>Z</given-names>
</name>
<name><surname>Wang</surname>
<given-names>N</given-names>
</name>
</person-group>
<article-title>Modelling disease outbreaks in realistic urban social networks</article-title>
<source>Nature</source>
<year>2004</year>
<volume>429</volume>
<fpage>180</fpage>
<lpage>184</lpage>
<pub-id pub-id-type="doi">10.1038/nature02541</pub-id>
<pub-id pub-id-type="pmid">15141212</pub-id>
</element-citation>
</ref>
<ref id="CR65"><label>65.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Halloran</surname>
<given-names>ME</given-names>
</name>
<name><surname>Ferguson</surname>
<given-names>NM</given-names>
</name>
<name><surname>Eubank</surname>
<given-names>S</given-names>
</name>
<name><surname>Longini</surname>
<given-names>IM</given-names>
</name>
<name><surname>Cummings</surname>
<given-names>DAT</given-names>
</name>
<name><surname>Lewis</surname>
<given-names>B</given-names>
</name>
<name><surname>Xu</surname>
<given-names>S</given-names>
</name>
<name><surname>Fraser</surname>
<given-names>C</given-names>
</name>
<name><surname>Vullikanti</surname>
<given-names>A</given-names>
</name>
<name><surname>Germann</surname>
<given-names>TC</given-names>
</name>
<name><surname>Wagener</surname>
<given-names>D</given-names>
</name>
<name><surname>Beckman</surname>
<given-names>R</given-names>
</name>
<name><surname>Kadau</surname>
<given-names>K</given-names>
</name>
<name><surname>Macken</surname>
<given-names>BDS</given-names>
</name>
<name><surname>Cooley</surname>
<given-names>P</given-names>
</name>
</person-group>
<article-title>Modeling targeted layered containment of an influenza pandemic in the United States</article-title>
<source>Proc Natl Acad Sci USA</source>
<year>2008</year>
<volume>105</volume>
<fpage>4639</fpage>
<lpage>4644</lpage>
<pub-id pub-id-type="doi">10.1073/pnas.0706849105</pub-id>
<pub-id pub-id-type="pmid">18332436</pub-id>
</element-citation>
</ref>
<ref id="CR66"><label>66.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Keeling</surname>
<given-names>MJ</given-names>
</name>
<name><surname>Grenfell</surname>
<given-names>BT</given-names>
</name>
</person-group>
<article-title>Effect of variability in infectious period on the persistence and spatial spread of infectious diseases</article-title>
<source>Math Biosci</source>
<year>1998</year>
<volume>147</volume>
<fpage>207</fpage>
<lpage>226</lpage>
<pub-id pub-id-type="doi">10.1016/S0025-5564(97)00101-6</pub-id>
<pub-id pub-id-type="pmid">9433063</pub-id>
</element-citation>
</ref>
<ref id="CR67"><label>67.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Vergu</surname>
<given-names>E</given-names>
</name>
<name><surname>Busson</surname>
<given-names>H</given-names>
</name>
<name><surname>Ezanno</surname>
<given-names>P</given-names>
</name>
</person-group>
<article-title>Impact of the Infection Period Distribution on the Epidemic Spread in a Metapopulation Model</article-title>
<source>PLoS One</source>
<year>2010</year>
<volume>5</volume>
<issue>2</issue>
<fpage>e9371</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pone.0009371</pub-id>
<pub-id pub-id-type="pmid">20195473</pub-id>
</element-citation>
</ref>
<ref id="CR68"><label>68.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tizzoni</surname>
<given-names>M</given-names>
</name>
<name><surname>Bajardi</surname>
<given-names>P</given-names>
</name>
<name><surname>Poletto</surname>
<given-names>C</given-names>
</name>
<name><surname>Ramasco</surname>
<given-names>JJ</given-names>
</name>
<name><surname>Balcan</surname>
<given-names>D</given-names>
</name>
<name><surname>Goncalves</surname>
<given-names>B</given-names>
</name>
<name><surname>Perra</surname>
<given-names>N</given-names>
</name>
<name><surname>Colizza</surname>
<given-names>V</given-names>
</name>
<name><surname>Vespignani</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>Real-time numerical forecast of global epidemic spreading: case study of 2009 A/H1N1pdm</article-title>
<source>BMC Med</source>
<year>2012</year>
<volume>10</volume>
<fpage>165</fpage>
<pub-id pub-id-type="doi">10.1186/1741-7015-10-165</pub-id>
<pub-id pub-id-type="pmid">23237460</pub-id>
</element-citation>
</ref>
<ref id="CR69"><label>69.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Funk</surname>
<given-names>S</given-names>
</name>
<name><surname>Salathe</surname>
<given-names>M</given-names>
</name>
<name><surname>Jansen</surname>
<given-names>VAA</given-names>
</name>
</person-group>
<article-title>Modeling the influence of human behavior on the spread of infectious diseases: a review</article-title>
<source>J Roy Soc Interface</source>
<year>2010</year>
<volume>7</volume>
<issue>50</issue>
<fpage>1247</fpage>
<lpage>1256</lpage>
<pub-id pub-id-type="doi">10.1098/rsif.2010.0142</pub-id>
<pub-id pub-id-type="pmid">20504800</pub-id>
</element-citation>
</ref>
<ref id="CR70"><label>70.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Salathé</surname>
<given-names>M</given-names>
</name>
<name><surname>Khandelwal</surname>
<given-names>S</given-names>
</name>
</person-group>
<article-title>Assessing Vaccination Sentiments with Online Social Media: Implications for Infectious Disease Dynamics and Control</article-title>
<source>PLoS Comput Biol</source>
<year>2011</year>
<volume>7</volume>
<issue>10</issue>
<fpage>e1002199</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pcbi.1002199</pub-id>
<pub-id pub-id-type="pmid">22022249</pub-id>
</element-citation>
</ref>
<ref id="CR71"><label>71.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Perra</surname>
<given-names>N</given-names>
</name>
<name><surname>Balcan</surname>
<given-names>D</given-names>
</name>
<name><surname>Goncalves</surname>
<given-names>B</given-names>
</name>
<name><surname>Vespignani</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>Towards a characterization of behavior-disease models</article-title>
<source>PLoS One</source>
<year>2011</year>
<volume>6</volume>
<issue>8</issue>
<fpage>e23084</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pone.0023084</pub-id>
<pub-id pub-id-type="pmid">21826228</pub-id>
</element-citation>
</ref>
<ref id="CR72"><label>72.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bagnoli</surname>
<given-names>F</given-names>
</name>
<name><surname>Lio</surname>
<given-names>P</given-names>
</name>
<name><surname>Sguanci</surname>
<given-names>L</given-names>
</name>
</person-group>
<article-title>Risk perception in epidemic modeling</article-title>
<source>Phys Rev E</source>
<year>2007</year>
<volume>76</volume>
<fpage>061904</fpage>
<pub-id pub-id-type="doi">10.1103/PhysRevE.76.061904</pub-id>
</element-citation>
</ref>
<ref id="CR73"><label>73.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Poletti</surname>
<given-names>P</given-names>
</name>
<name><surname>Caprile</surname>
<given-names>B</given-names>
</name>
<name><surname>Ajelli</surname>
<given-names>M</given-names>
</name>
<name><surname>Pugliese</surname>
<given-names>A</given-names>
</name>
<name><surname>Merler</surname>
<given-names>S</given-names>
</name>
</person-group>
<article-title>Spontaneous behavioural changes in response to epidemics</article-title>
<source>J Theor Biol</source>
<year>2009</year>
<volume>260</volume>
<fpage>31</fpage>
<pub-id pub-id-type="doi">10.1016/j.jtbi.2009.04.029</pub-id>
<pub-id pub-id-type="pmid">19446570</pub-id>
</element-citation>
</ref>
<ref id="CR74"><label>74.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Fumanelli</surname>
<given-names>L</given-names>
</name>
<name><surname>Ajelli</surname>
<given-names>M</given-names>
</name>
<name><surname>Manfredi</surname>
<given-names>P</given-names>
</name>
<name><surname>Vespignani</surname>
<given-names>A</given-names>
</name>
<name><surname>Merler</surname>
<given-names>S</given-names>
</name>
</person-group>
<article-title>Inferring the Structure of Social Contacts from Demographic Data in the Analysis of Infectious Diseases Spread</article-title>
<source>PLoS Comput Biol</source>
<year>2012</year>
<volume>8</volume>
<issue>9</issue>
<fpage>e1002673</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pcbi.1002673</pub-id>
<pub-id pub-id-type="pmid">23028275</pub-id>
</element-citation>
</ref>
<ref-list id="BSec1"><title>Pre-publication history</title>
<ref id="CR75"><mixed-citation publication-type="other">The pre-publication history for this paper can be accessed here:<ext-link ext-link-type="uri" xlink:href="http://www.biomedcentral.com/1471-2334/13/176/prepub">http://www.biomedcentral.com/1471-2334/13/176/prepub</ext-link>
</mixed-citation>
</ref>
</ref-list>
</ref-list>
</back>
</pmc>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Sante/explor/PandemieGrippaleV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 0001829 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 0001829 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Sante |area= PandemieGrippaleV1 |flux= Pmc |étape= Corpus |type= RBID |clé= |texte= }}
![]() | This area was generated with Dilib version V0.6.34. | ![]() |