Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Global migration of influenza A viruses in swine</title><author><name sortKey="Nelson, Martha I" sort="Nelson, Martha I" uniqKey="Nelson M" first="Martha I." last="Nelson">Martha I. Nelson</name><affiliation><nlm:aff id="A1">Division of International Epidemiology and Population Studies, Fogarty International Center, National Institutes of Health, Bethesda, MD 20892 USA</nlm:aff></affiliation></author><author><name sortKey="Viboud, Cecile" sort="Viboud, Cecile" uniqKey="Viboud C" first="Cécile" last="Viboud">Cécile Viboud</name><affiliation><nlm:aff id="A1">Division of International Epidemiology and Population Studies, Fogarty International Center, National Institutes of Health, Bethesda, MD 20892 USA</nlm:aff></affiliation></author><author><name sortKey="Vincent, Amy L" sort="Vincent, Amy L" uniqKey="Vincent A" first="Amy L." last="Vincent">Amy L. Vincent</name><affiliation><nlm:aff id="A2">Virus and Prion Diseases of Livestock Research Unit, National Animal Disease Center, USDA-ARS, Ames, IA 50010 USA</nlm:aff></affiliation></author><author><name sortKey="Culhane, Marie R" sort="Culhane, Marie R" uniqKey="Culhane M" first="Marie R." last="Culhane">Marie R. Culhane</name><affiliation><nlm:aff id="A3">University of Minnesota Veterinary Diagnostic Laboratory, St. Paul, MN 55108 USA</nlm:aff></affiliation></author><author><name sortKey="Detmer, Susan E" sort="Detmer, Susan E" uniqKey="Detmer S" first="Susan E." last="Detmer">Susan E. Detmer</name><affiliation><nlm:aff id="A4">Western College of Veterinary Medicine, University of Saskatchewan, Saskatchewan, Canada</nlm:aff></affiliation></author><author><name sortKey="Wentworth, David E" sort="Wentworth, David E" uniqKey="Wentworth D" first="David E." last="Wentworth">David E. Wentworth</name><affiliation><nlm:aff id="A5">J Craig Venter Institute, Rockville, MD USA 20850 USA</nlm:aff></affiliation></author><author><name sortKey="Rambaut, Andrew" sort="Rambaut, Andrew" uniqKey="Rambaut A" first="Andrew" last="Rambaut">Andrew Rambaut</name><affiliation><nlm:aff id="A1">Division of International Epidemiology and Population Studies, Fogarty International Center, National Institutes of Health, Bethesda, MD 20892 USA</nlm:aff></affiliation><affiliation><nlm:aff id="A6">Institute of Evolutionary Biology, University of Edinburgh, Ashworth Laboratories, Edinburgh, EH9 FLT United Kingdom</nlm:aff></affiliation><affiliation><nlm:aff id="A7">Centre for Immunology, Infection and Evolution, University of Edinburgh, Ashworth Laboratories, Edinburgh, EH9 FLT United Kingdom</nlm:aff></affiliation></author><author><name sortKey="Suchard, Marc A" sort="Suchard, Marc A" uniqKey="Suchard M" first="Marc A." last="Suchard">Marc A. Suchard</name><affiliation><nlm:aff id="A8">Departments of Biomathematics and Human Genetics, David Geffen School of Medicine at UCLA, and Department of Biostatistics, UCLA School of Public Health, Los Angeles, CA 90095 USA</nlm:aff></affiliation></author><author><name sortKey="Holmes, Edward C" sort="Holmes, Edward C" uniqKey="Holmes E" first="Edward C." last="Holmes">Edward C. Holmes</name><affiliation><nlm:aff id="A9">Marie Bashir Institute for Infectious Diseases and Biosecurity, Charles Perkins Centre, School of Biological Sciences and Sydney Medical School, University of Sydney, NSW 2006, Australia</nlm:aff></affiliation></author><author><name sortKey="Lemey, Philippe" sort="Lemey, Philippe" uniqKey="Lemey P" first="Philippe" last="Lemey">Philippe Lemey</name><affiliation><nlm:aff id="A10">Department of Microbiology and Immunology, Rega Institute, KU Leuven – University of Leuven, Leuven, Belgium</nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">25813399</idno><idno type="pmc">4380236</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4380236</idno><idno type="RBID">PMC:4380236</idno><idno type="doi">10.1038/ncomms7696</idno><date when="2015">2015</date><idno type="wicri:Area/Pmc/Corpus">000E59</idno><idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000E59</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Global migration of influenza A viruses in swine</title><author><name sortKey="Nelson, Martha I" sort="Nelson, Martha I" uniqKey="Nelson M" first="Martha I." last="Nelson">Martha I. Nelson</name><affiliation><nlm:aff id="A1">Division of International Epidemiology and Population Studies, Fogarty International Center, National Institutes of Health, Bethesda, MD 20892 USA</nlm:aff></affiliation></author><author><name sortKey="Viboud, Cecile" sort="Viboud, Cecile" uniqKey="Viboud C" first="Cécile" last="Viboud">Cécile Viboud</name><affiliation><nlm:aff id="A1">Division of International Epidemiology and Population Studies, Fogarty International Center, National Institutes of Health, Bethesda, MD 20892 USA</nlm:aff></affiliation></author><author><name sortKey="Vincent, Amy L" sort="Vincent, Amy L" uniqKey="Vincent A" first="Amy L." last="Vincent">Amy L. Vincent</name><affiliation><nlm:aff id="A2">Virus and Prion Diseases of Livestock Research Unit, National Animal Disease Center, USDA-ARS, Ames, IA 50010 USA</nlm:aff></affiliation></author><author><name sortKey="Culhane, Marie R" sort="Culhane, Marie R" uniqKey="Culhane M" first="Marie R." last="Culhane">Marie R. Culhane</name><affiliation><nlm:aff id="A3">University of Minnesota Veterinary Diagnostic Laboratory, St. Paul, MN 55108 USA</nlm:aff></affiliation></author><author><name sortKey="Detmer, Susan E" sort="Detmer, Susan E" uniqKey="Detmer S" first="Susan E." last="Detmer">Susan E. Detmer</name><affiliation><nlm:aff id="A4">Western College of Veterinary Medicine, University of Saskatchewan, Saskatchewan, Canada</nlm:aff></affiliation></author><author><name sortKey="Wentworth, David E" sort="Wentworth, David E" uniqKey="Wentworth D" first="David E." last="Wentworth">David E. Wentworth</name><affiliation><nlm:aff id="A5">J Craig Venter Institute, Rockville, MD USA 20850 USA</nlm:aff></affiliation></author><author><name sortKey="Rambaut, Andrew" sort="Rambaut, Andrew" uniqKey="Rambaut A" first="Andrew" last="Rambaut">Andrew Rambaut</name><affiliation><nlm:aff id="A1">Division of International Epidemiology and Population Studies, Fogarty International Center, National Institutes of Health, Bethesda, MD 20892 USA</nlm:aff></affiliation><affiliation><nlm:aff id="A6">Institute of Evolutionary Biology, University of Edinburgh, Ashworth Laboratories, Edinburgh, EH9 FLT United Kingdom</nlm:aff></affiliation><affiliation><nlm:aff id="A7">Centre for Immunology, Infection and Evolution, University of Edinburgh, Ashworth Laboratories, Edinburgh, EH9 FLT United Kingdom</nlm:aff></affiliation></author><author><name sortKey="Suchard, Marc A" sort="Suchard, Marc A" uniqKey="Suchard M" first="Marc A." last="Suchard">Marc A. Suchard</name><affiliation><nlm:aff id="A8">Departments of Biomathematics and Human Genetics, David Geffen School of Medicine at UCLA, and Department of Biostatistics, UCLA School of Public Health, Los Angeles, CA 90095 USA</nlm:aff></affiliation></author><author><name sortKey="Holmes, Edward C" sort="Holmes, Edward C" uniqKey="Holmes E" first="Edward C." last="Holmes">Edward C. Holmes</name><affiliation><nlm:aff id="A9">Marie Bashir Institute for Infectious Diseases and Biosecurity, Charles Perkins Centre, School of Biological Sciences and Sydney Medical School, University of Sydney, NSW 2006, Australia</nlm:aff></affiliation></author><author><name sortKey="Lemey, Philippe" sort="Lemey, Philippe" uniqKey="Lemey P" first="Philippe" last="Lemey">Philippe Lemey</name><affiliation><nlm:aff id="A10">Department of Microbiology and Immunology, Rega Institute, KU Leuven – University of Leuven, Leuven, Belgium</nlm:aff></affiliation></author></analytic><series><title level="j">Nature communications</title><idno type="eISSN">2041-1723</idno><imprint><date when="2015">2015</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><p id="P1">The complex and unresolved evolutionary origins of the 2009 H1N1 influenza pandemic exposed major gaps in our knowledge of the global spatial ecology and evolution of influenza A viruses in swine (swIAVs). Here we undertake an expansive phylogenetic analysis of swIAV sequence data and demonstrate that the global live swine trade strongly predicts the spatial dissemination of swIAVs, with Europe and North America acting as sources of viruses in Asian countries. In contrast, China has the world’s largest swine population but is not a major exporter of live swine, and is not an important source of swIAVs in neighboring Asian countries or globally. A meta-population simulation model incorporating trade data predicts that the global ecology of swIAVs is more complex than previously thought, and the US and China’s large swine populations are unlikely to be representative of swIAV diversity in their respective geographic regions, requiring independent surveillance efforts throughout Latin America and Asia.</p></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Garten, Rj" uniqKey="Garten R">RJ Garten</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Smith, Gjd" uniqKey="Smith G">GJD Smith</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lemey, P" uniqKey="Lemey P">P Lemey</name></author><author><name sortKey="Suchard, M" uniqKey="Suchard M">M Suchard</name></author><author><name sortKey="Rambaut, A" uniqKey="Rambaut A">A Rambaut</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Koen, Js" uniqKey="Koen J">JS Koen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhou, Nn" uniqKey="Zhou N">NN Zhou</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Karasin, A" uniqKey="Karasin A">A Karasin</name></author><author><name sortKey="Carman, S" uniqKey="Carman S">S Carman</name></author><author><name sortKey="Olsen, C" uniqKey="Olsen C">C Olsen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Vincent, Al" uniqKey="Vincent A">AL Vincent</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ottis, K" uniqKey="Ottis K">K Ottis</name></author><author><name sortKey="Sidoli, L" uniqKey="Sidoli L">L Sidoli</name></author><author><name sortKey="Bachmann, P" uniqKey="Bachmann P">P Bachmann</name></author><author><name sortKey="Webster, Rg" uniqKey="Webster R">RG Webster</name></author><author><name sortKey="Kaplan, Mm" uniqKey="Kaplan M">MM Kaplan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Vijaykrishna, D" uniqKey="Vijaykrishna D">D Vijaykrishna</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Takemae, N" uniqKey="Takemae N">N Takemae</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Poonsuk, S" uniqKey="Poonsuk S">S Poonsuk</name></author><author><name sortKey="Sangthong, P" uniqKey="Sangthong P">P Sangthong</name></author><author><name sortKey="Petcharat, N" uniqKey="Petcharat N">N Petcharat</name></author><author><name sortKey="Lekcharoensuk, P" uniqKey="Lekcharoensuk P">P Lekcharoensuk</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Song, D" uniqKey="Song D">D Song</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nelson, Mi" uniqKey="Nelson M">MI Nelson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zell, R" uniqKey="Zell R">R Zell</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nelson, Mi" uniqKey="Nelson M">MI Nelson</name></author><author><name sortKey="Vincent, Al" uniqKey="Vincent A">AL Vincent</name></author><author><name sortKey="Kitikoon, P" uniqKey="Kitikoon P">P Kitikoon</name></author><author><name sortKey="Holmes, Ec" uniqKey="Holmes E">EC Holmes</name></author><author><name sortKey="Gramer, Mr" uniqKey="Gramer M">MR Gramer</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ducatez, Mf" uniqKey="Ducatez M">MF Ducatez</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nelson, Mi" uniqKey="Nelson M">MI Nelson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shi, M" uniqKey="Shi M">M Shi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Anderson, Tk" uniqKey="Anderson T">TK Anderson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Minin, Vn" uniqKey="Minin V">VN Minin</name></author><author><name sortKey="Suchard, Ma" uniqKey="Suchard M">MA Suchard</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lemey, P" uniqKey="Lemey P">P Lemey</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bielejec, F" uniqKey="Bielejec F">F Bielejec</name></author><author><name sortKey="Lemey, P" uniqKey="Lemey P">P Lemey</name></author><author><name sortKey="Baele, G" uniqKey="Baele G">G Baele</name></author><author><name sortKey="Rambaut, A" uniqKey="Rambaut A">A Rambaut</name></author><author><name sortKey="Suchard, Ma" uniqKey="Suchard M">MA Suchard</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Russell, Ca" uniqKey="Russell C">CA Russell</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nelson, Mi" uniqKey="Nelson M">MI Nelson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nelson, Mi" uniqKey="Nelson M">MI Nelson</name></author><author><name sortKey="Vincent, Al" uniqKey="Vincent A">AL Vincent</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Holyoake, Pk" uniqKey="Holyoake P">PK Holyoake</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Deng, Y" uniqKey="Deng Y">Y Deng</name></author><author><name sortKey="Iannello, P" uniqKey="Iannello P">P Iannello</name></author><author><name sortKey="Smith, I" uniqKey="Smith I">I Smith</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Biondo, N" uniqKey="Biondo N">N Biondo</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Schaefer, R" uniqKey="Schaefer R">R Schaefer</name></author><author><name sortKey="Zanella, Jrc" uniqKey="Zanella J">JRC Zanella</name></author><author><name sortKey="Brentano, L" uniqKey="Brentano L">L Brentano</name></author><author><name sortKey="Vincent, Al" uniqKey="Vincent A">AL Vincent</name></author><author><name sortKey="Ritterbusch, Ga" uniqKey="Ritterbusch G">GA Ritterbusch</name></author><author><name sortKey="Silveira, S" uniqKey="Silveira S">S Silveira</name></author><author><name sortKey="Caron, L" uniqKey="Caron L">L Caron</name></author><author><name sortKey="M, N" uniqKey="M N">N M</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nagarajan, K" uniqKey="Nagarajan K">K Nagarajan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Njabo, Ky" uniqKey="Njabo K">KY Njabo</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Escalera Zamudio, M" uniqKey="Escalera Zamudio M">M Escalera-Zamudio</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Meseko, Ca" uniqKey="Meseko C">CA Meseko</name></author><author><name sortKey="Odaibo, Gn" uniqKey="Odaibo G">GN Odaibo</name></author><author><name sortKey="Olaleye, Do" uniqKey="Olaleye D">DO Olaleye</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Perera, Hkk" uniqKey="Perera H">HKK Perera</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hofshagen, M" uniqKey="Hofshagen M">M Hofshagen</name></author><author><name sortKey="Gjerset, B" uniqKey="Gjerset B">B Gjerset</name></author><author><name sortKey="Er, C" uniqKey="Er C">C Er</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nokireki, T" uniqKey="Nokireki T">T Nokireki</name></author><author><name sortKey="Laine, T" uniqKey="Laine T">T Laine</name></author><author><name sortKey="London, L" uniqKey="London L">L London</name></author><author><name sortKey="Ikonen, N" uniqKey="Ikonen N">N Ikonen</name></author><author><name sortKey="Huovilainen, A" uniqKey="Huovilainen A">A Huovilainen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Banyai, K" uniqKey="Banyai K">K Bányai</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lam, Tty" uniqKey="Lam T">TTY Lam</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Cappuccio, Ja" uniqKey="Cappuccio J">Ja Cappuccio</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Schaefer, R" uniqKey="Schaefer R">R Schaefer</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhou, B" uniqKey="Zhou B">B Zhou</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bao, Y" uniqKey="Bao Y">Y Bao</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Edgar, Rc" uniqKey="Edgar R">RC Edgar</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nelson, Mi" uniqKey="Nelson M">MI Nelson</name></author><author><name sortKey="Gramer, Mr" uniqKey="Gramer M">MR Gramer</name></author><author><name sortKey="Vincent, Al" uniqKey="Vincent A">AL Vincent</name></author><author><name sortKey="Holmes, Ec" uniqKey="Holmes E">EC Holmes</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Drummond, Aj" uniqKey="Drummond A">AJ Drummond</name></author><author><name sortKey="Suchard, Ma" uniqKey="Suchard M">MA Suchard</name></author><author><name sortKey="Xie, D" uniqKey="Xie D">D Xie</name></author><author><name sortKey="Rambaut, A" uniqKey="Rambaut A">A Rambaut</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Suchard, Ma" uniqKey="Suchard M">MA Suchard</name></author><author><name sortKey="Rambaut, A" uniqKey="Rambaut A">A Rambaut</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Viboud, C" uniqKey="Viboud C">C Viboud</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Romagosa, A" uniqKey="Romagosa A">A Romagosa</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wong, Kk" uniqKey="Wong K">KK Wong</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Robinson, Tp" uniqKey="Robinson T">TP Robinson</name></author></analytic></biblStruct></listBibl></div1></back></TEI><pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<pmc-dir>properties manuscript</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-journal-id">101528555</journal-id><journal-id journal-id-type="pubmed-jr-id">37539</journal-id><journal-id journal-id-type="nlm-ta">Nat Commun</journal-id><journal-id journal-id-type="iso-abbrev">Nat Commun</journal-id><journal-title-group><journal-title>Nature communications</journal-title></journal-title-group><issn pub-type="epub">2041-1723</issn></journal-meta><article-meta><article-id pub-id-type="pmid">25813399</article-id><article-id pub-id-type="pmc">4380236</article-id><article-id pub-id-type="doi">10.1038/ncomms7696</article-id><article-id pub-id-type="manuscript">NIHMS666251</article-id><article-categories><subj-group subj-group-type="heading"><subject>Article</subject></subj-group></article-categories><title-group><article-title>Global migration of influenza A viruses in swine</article-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Nelson</surname><given-names>Martha I.</given-names></name><xref ref-type="aff" rid="A1">1</xref><xref rid="FN1" ref-type="author-notes">*</xref></contrib><contrib contrib-type="author"><name><surname>Viboud</surname><given-names>Cécile</given-names></name><xref ref-type="aff" rid="A1">1</xref></contrib><contrib contrib-type="author"><name><surname>Vincent</surname><given-names>Amy L.</given-names></name><xref ref-type="aff" rid="A2">2</xref></contrib><contrib contrib-type="author"><name><surname>Culhane</surname><given-names>Marie R.</given-names></name><xref ref-type="aff" rid="A3">3</xref></contrib><contrib contrib-type="author"><name><surname>Detmer</surname><given-names>Susan E.</given-names></name><xref ref-type="aff" rid="A4">4</xref></contrib><contrib contrib-type="author"><name><surname>Wentworth</surname><given-names>David E.</given-names></name><xref ref-type="aff" rid="A5">5</xref><xref rid="FN2" ref-type="author-notes">#</xref></contrib><contrib contrib-type="author"><name><surname>Rambaut</surname><given-names>Andrew</given-names></name><xref ref-type="aff" rid="A1">1</xref><xref ref-type="aff" rid="A6">6</xref><xref ref-type="aff" rid="A7">7</xref></contrib><contrib contrib-type="author"><name><surname>Suchard</surname><given-names>Marc A.</given-names></name><xref ref-type="aff" rid="A8">8</xref></contrib><contrib contrib-type="author"><name><surname>Holmes</surname><given-names>Edward C.</given-names></name><xref ref-type="aff" rid="A9">9</xref></contrib><contrib contrib-type="author"><name><surname>Lemey</surname><given-names>Philippe</given-names></name><xref ref-type="aff" rid="A10">10</xref></contrib></contrib-group><aff id="A1"><label>1</label>Division of International Epidemiology and Population Studies, Fogarty International Center, National Institutes of Health, Bethesda, MD 20892 USA</aff><aff id="A2"><label>2</label>Virus and Prion Diseases of Livestock Research Unit, National Animal Disease Center, USDA-ARS, Ames, IA 50010 USA</aff><aff id="A3"><label>3</label>University of Minnesota Veterinary Diagnostic Laboratory, St. Paul, MN 55108 USA</aff><aff id="A4"><label>4</label>Western College of Veterinary Medicine, University of Saskatchewan, Saskatchewan, Canada</aff><aff id="A5"><label>5</label>J Craig Venter Institute, Rockville, MD USA 20850 USA</aff><aff id="A6"><label>6</label>Institute of Evolutionary Biology, University of Edinburgh, Ashworth Laboratories, Edinburgh, EH9 FLT United Kingdom</aff><aff id="A7"><label>7</label>Centre for Immunology, Infection and Evolution, University of Edinburgh, Ashworth Laboratories, Edinburgh, EH9 FLT United Kingdom</aff><aff id="A8"><label>8</label>Departments of Biomathematics and Human Genetics, David Geffen School of Medicine at UCLA, and Department of Biostatistics, UCLA School of Public Health, Los Angeles, CA 90095 USA</aff><aff id="A9"><label>9</label>Marie Bashir Institute for Infectious Diseases and Biosecurity, Charles Perkins Centre, School of Biological Sciences and Sydney Medical School, University of Sydney, NSW 2006, Australia</aff><aff id="A10"><label>10</label>Department of Microbiology and Immunology, Rega Institute, KU Leuven – University of Leuven, Leuven, Belgium</aff><author-notes><corresp id="FN1"><label>*</label>Corresponding author: Division of International Epidemiology and Population Studies, Fogarty International Center, National Institutes of Health, 16 Center Drive, Building 16, Room 202, Bethesda, MD 20892 USA. Tel: 1 301-402-5203; Fax: 1 301-496-8496; <email>nelsonma@mail.nih.gov</email></corresp><fn id="FN2" fn-type="present-address"><label>#</label><p>Current address: Influenza Division, National Center for Immunization and Respiratory Diseases, Centers for Disease Control and Prevention, Atlanta, Georgia, United States of America</p></fn></author-notes><pub-date pub-type="nihms-submitted"><day>24</day><month>2</month><year>2015</year></pub-date><pub-date pub-type="epub"><day>27</day><month>3</month><year>2015</year></pub-date><pub-date pub-type="collection"><year>2015</year></pub-date><pub-date pub-type="pmc-release"><day>27</day><month>9</month><year>2015</year></pub-date><volume>6</volume><fpage>6696</fpage><lpage>6696</lpage><pmc-comment>elocation-id from pubmed: 10.1038/ncomms7696</pmc-comment>
<permissions><license xlink:href="http://www.nature.com/authors/editorial_policies/license.html#terms"><license-p>Users may view, print, copy, and download text and data-mine the content in such documents, for the purposes of academic research, subject always to the full Conditions of use:<ext-link ext-link-type="uri" xlink:href="http://www.nature.com/authors/editorial_policies/license.html#terms">http://www.nature.com/authors/editorial_policies/license.html#terms</ext-link></license-p></license></permissions><abstract><p id="P1">The complex and unresolved evolutionary origins of the 2009 H1N1 influenza pandemic exposed major gaps in our knowledge of the global spatial ecology and evolution of influenza A viruses in swine (swIAVs). Here we undertake an expansive phylogenetic analysis of swIAV sequence data and demonstrate that the global live swine trade strongly predicts the spatial dissemination of swIAVs, with Europe and North America acting as sources of viruses in Asian countries. In contrast, China has the world’s largest swine population but is not a major exporter of live swine, and is not an important source of swIAVs in neighboring Asian countries or globally. A meta-population simulation model incorporating trade data predicts that the global ecology of swIAVs is more complex than previously thought, and the US and China’s large swine populations are unlikely to be representative of swIAV diversity in their respective geographic regions, requiring independent surveillance efforts throughout Latin America and Asia.</p></abstract></article-meta></front><body><sec sec-type="intro" id="S1"><title>Introduction</title><p id="P2">In early 2009, a novel reassortant H1N1 influenza A virus with gene segments from two swine virus (swIAV) lineages emerged in humans, initiating the first influenza pandemic of the 21<sup>st</sup> century. The virus had a complex genetic composition that had not been previously detected in swine, with six genome segments of North American triple reassortant swine virus origin (PB2, PB1, PA, HA (H1), NP, and NS) and two genome segments of Eurasian avian-like swine virus origin (NA (N1) and MP)<sup><xref rid="R1" ref-type="bibr">1</xref></sup>. Evolutionary analysis of this novel North American/Eurasian reassortant virus indicated that these segments had circulated undetected in swine for at least eight years<sup><xref rid="R2" ref-type="bibr">2</xref></sup>. The first human outbreak of the pandemic H1N1 virus (pH1N1) occurred in Mexico, and the extent of viral genetic diversity observed in Mexico supports the hypothesis that the virus first emerged there in humans<sup><xref rid="R3" ref-type="bibr">3</xref></sup>. However, efforts to detect the last common ancestor of the pH1N1 virus in Mexican swine populations have not been successful to date, and the opaque evolutionary history of the pandemic virus in swine highlights the gaps in our understanding of swIAV dynamics at a global scale.</p><p id="P3">In general, influenza viruses in swine are spatially separated into distinct North American and European swIAV lineages, although viruses of North American and European origin both circulate in Asia. Multiple viral lineages co-circulate in North American swine, including (i) ‘classical’ swine viruses, which descend from the 1918 H1N1 pandemic<sup><xref rid="R4" ref-type="bibr">4</xref></sup>; (ii) ‘triple reassortant’ swine viruses, which emerged in the mid-1990s with a combination of human, swine, and avian segments<sup><xref rid="R5" ref-type="bibr">5</xref></sup>; and (iii) ‘delta’ (δ) viruses that are closely related to human seasonal H1 viruses from the early 2000s<sup><xref rid="R6" ref-type="bibr">6</xref>,<xref rid="R7" ref-type="bibr">7</xref></sup>. The main European swIAV lineages include ‘avian-like’ H1N1 viruses that jumped from birds to swine in the 1970s, human-origin H1N1 viruses from the 1980s, and human-origin H3N2 viruses that are antigenically described as A/Port Chalmers/1/1973-like<sup><xref rid="R8" ref-type="bibr">8</xref></sup>. Multiple North American and European-origin swIAV lineages have both been identified in Asian countries<sup><xref rid="R9" ref-type="bibr">9</xref>–<xref rid="R12" ref-type="bibr">12</xref></sup>. Due to high levels of co-infection, segmental reassortment occurs frequently in swine, such that they are an important reservoir host for influenza virus genetic diversity<sup><xref rid="R9" ref-type="bibr">9</xref>,<xref rid="R11" ref-type="bibr">11</xref>,<xref rid="R13" ref-type="bibr">13</xref>–<xref rid="R16" ref-type="bibr">16</xref></sup>.</p><p id="P4">Live transport is routine in swine farming, and in the United States the transport of millions of swine from Southern to Midwestern regions for end-stage production appears to drive the strongly directional dissemination of swIAVs from Southern US states with high hog production (e.g., North Carolina, Texas, and Oklahoma) to the traditional center of swine farming located in the Midwestern ‘corn belt’<sup><xref rid="R17" ref-type="bibr">17</xref></sup>. Large numbers of swine also enter the United States from Canada, which has been implicated in the dissemination of other important swine pathogens, including Porcine Reproductive and Respiratory Syndrome Virus (PRRSV)<sup><xref rid="R18" ref-type="bibr">18</xref></sup>. Intercontinental trade of live swine also occurs, for end-stage production or to acquire female breeding pigs for genetic improvement of swine reproduction or growth traits. Globally, the largest swine population is found in China, where over 450 million hogs reside (<xref rid="F1" ref-type="fig">Fig. 1</xref>). Large swine populations also are found in the United States (> 60 million hogs), Brazil (> 30 million hogs), Vietnam (> 20 million hogs), Germany (> 20 million hogs), and Spain (> 20 million hogs).</p><p id="P5">Despite the global nature of both swine farming and swIAV circulation, the patterns and dynamics of the worldwide spread of this economically important virus are unknown. To characterize the phylodynamics and phylogeography of swIAVs at a global scale, here we conduct a phylogenetic analysis of 785 whole-genome swIAV sequences collected from ten countries/regions representing four continents, the largest study of its kind undertaken to date. To assess the drivers of viral migration, we compare the phylogeographic patterns with empirical data on live swine trade and swine population sizes. Based on these findings, we build a meta-population model to simulate the spatial dissemination of swIAVs at a global scale and identify regions at high risk for co-invasion of divergent lineages, increased total genetic diversity, and emergence of viruses with pandemic potential.</p></sec><sec sec-type="results" id="S2"><title>Results</title><sec id="S3"><title>Global migration of swIAVs</title><p id="P6">Phylogenetic analysis revealed that long-distance movement of influenza A viruses between countries and continents has occurred continuously in swine since the 1970s (summarized in <xref rid="F2" ref-type="fig">Fig. 2</xref>). Our estimate of 18 international viral migration events is a minimum based on the currently available influenza virus sequence data and certainly underestimates the true number. This lower-bound estimate is based on discrete monophyletic groups (defined by country) that are supported by high posterior probabilities (> 80), reflecting migration events that led to successful onward transmission. The estimate does not include the much higher number of international viral movements between the United States and Canada or between countries in Europe, which are each considered as meta-populations in our analysis. The estimate also does not consider viral migration events for which only one sequence is available, any viruses that do not form a well-supported clade, or which only partial sequence data was available.</p><p id="P7">Although global surveillance and sequencing of swIAVs has increased dramatically in the last five years, our time-scaled MCC phylogenies indicate that most intercontinental viral migration events occurred prior to this increase in surveillance (representative phylogeny of the NA segment is presented in <xref rid="F3" ref-type="fig">Fig. 3</xref>; phylogenies for other seven segments are available in <xref ref-type="supplementary-material" rid="SD1">Supplementary Figs. 1–7</xref>). Eight of the 18 viral migration events identified in our study were evident on the phylogenies inferred for all eight viral genome segments, indicative of onward transmission of the full viral genome in the new location (introductions 1, 3, 4, 5, 8, 9, 10, 15). The consistency of spatial-temporal inferences across these eight segments strengthens inferences of when and where each of these introductions occurred (<xref ref-type="supplementary-material" rid="SD1">Supplementary Table 1</xref>). Ten viral migration events could only be identified by a subset of genome segments, as at least one segment has been replaced in intervening years by reassortment (introductions 2, 6, 7, 11, 12, 13, 14, 16, 17, and 18). There is no evidence in these data that either of the δ-1 or δ-2 virus lineages that emerged in North American swine in the early 2000s has transmitted to swine on any other continents, despite high rates of detection in US swine populations of δ-1 viruses<sup><xref rid="R19" ref-type="bibr">19</xref></sup>.</p><p id="P8">A consistent global spatial dynamic was observed for swIAVs during 1970–2013, based on both a conservative measure of strongly supported monophyletic groups (<xref rid="F3" ref-type="fig">Fig. 3</xref> and <xref ref-type="supplementary-material" rid="SD1">Supplementary Figs. 1–7</xref>) as well as ‘Markov jump’ counts<sup><xref rid="R20" ref-type="bibr">20</xref></sup> of the expected number of location state transitions along the branches of the tree. ‘Markov jump’ counts provide a quantitative measure of gene flow between regions that includes singletons and clusters that may have poor phylogenetic resolution (<xref rid="F4" ref-type="fig">Fig. 4</xref>, <xref ref-type="supplementary-material" rid="SD1">Supplementary Table 2</xref>). Overall, North America (in this case referring to the United States and Canada) and Europe represent independent viral source populations for the Asian countries sampled in our study: China, Japan, South Korea, Thailand, and Vietnam. In contrast, only a single swIAV migration event was observed between North America and Europe (introduction 5).</p><sec id="S4"><title>Spatial dynamics of swIAVs in North America</title><p id="P9">Bi-directional viral migration between Canada and the United States is so frequent (reflected by the extremely high number of Markov jump counts, <xref rid="F4" ref-type="fig">Fig. 4</xref>) that Canada and the US were considered as a single meta-population, similar to Europe (<xref rid="F2" ref-type="fig">Fig. 2</xref>). The higher availability of swIAV sequence data from US swine makes it particularly difficult to distinguish whether an introduction was specifically of US or Canadian swine origin, and the origin of such introductions is more appropriately characterized as ‘North American’. Using newly generated sequence data from five swIAVs of the H3N2 subtype that were collected in Mexico during 2010–2011 (A/sw/Mexico/SG1442/2010, A/sw/Mexico/SG1444/2011, A/sw/Mexico/SG1447/2011, A/sw/Mexico/SG1448/2011, and A/sw/Mexico/SG1449/2011, accession numbers available in <xref ref-type="supplementary-material" rid="SD1">Supplementary Table 3</xref>), our phylogenetic analysis provides evidence of a single introduction of a H3N2 triple reassortant swIAV from the United States into Mexico that occurred between mid-2005 and mid-2006 (introduction 11, <xref rid="F2" ref-type="fig">Fig. 2</xref>). At the time of sampling, all five Mexican swIAVs had acquired at least one pH1N1 segment of human origin via reassortment, evidence that pH1N1 viruses also have circulated in Mexico’s swine population.</p></sec></sec><sec id="S5"><title>Spatial dynamics of swIAVs in Asia</title><p id="P10">Given China’s large swine population and long-term surveillance, it may have been expected that this country (encompassing mainland China, Hong Kong, Taiwan, and Macau) would be an important source of swIAV diversity for neighboring Asian countries in this study. However, since 1970 there have been 16 swIAV introductions of European or North American origin into Asia, compared to only two swIAV migration events between Asian countries, and only one definitive introduction of a swIAV from China into another country (introduction 16). Overall, the genetic diversity of swIAVs in Asia derives from five swIAV introductions of European origin and eleven swIAV introductions of North American origin. Six viral introductions from Europe and North America were observed in Thailand, and five introductions were observed in China, including the earliest intercontinental swIAV migration detected to date (introduction 1, <xref rid="F2" ref-type="fig">Fig. 2</xref>).</p><p id="P11">In contrast to the frequent exchange of swIAVs across European country borders and across the US-Canadian border, only two swIAV migration events were observed between any two Asian countries. A pair of H1 segments collected in South Korea in 2013 are positioned within a clade of avian-origin Eurasian viruses from China, suggesting China-to-Korea migration (introduction 16, <xref rid="F2" ref-type="fig">Fig. 2</xref> and <xref ref-type="supplementary-material" rid="SD1">Supplementary Fig. 4</xref>). Closely related North American triple reassortant viruses also were identified in China and Vietnam, suggestive of viral migration between these Asian countries. However, the location state probability for the node representing the common ancestor of the Chinese and Vietnamese clades is too low (ranging from 0.46 to 0.65 across the eight genome segments) to determine whether the North American virus was first introduced to Vietnam and disseminated to China, or vice versa (<xref rid="F2" ref-type="fig">Fig. 2</xref> and <xref ref-type="supplementary-material" rid="SD1">Supplementary Figs. 1–7</xref>).</p><p id="P12">Although much of the swIAV diversity in Asia appears to have emerged in the last two decades, the phylogenies suggest long-term circulation of swIAVs in Thailand and Japan, either via imports from North American or European swine in the 1970s and 1980s (<xref rid="F2" ref-type="fig">Fig. 2</xref>) or direct introductions from humans as early as the 1960s (<xref rid="F3" ref-type="fig">Fig. 3c</xref>). Long branch lengths and the lack of historical swIAV data from Asian countries make it difficult to infer with confidence the spatial history of older viral lineages in Asia. The relative lack of swIAV surveillance in Thailand prior to 2000 particularly complicates estimates of the timing and spatial pathway of the multiple viral introductions from North America and Europe that likely occurred during the 1980s and 1990s (e.g., introductions 6, 14, 17, and 18). At this time, there is little evidence of viral dissemination from Japan or Thailand to other Asian countries in our study, despite many decades of potential swIAV circulation. However, the relatively long branches may not reflect a single direct transition between the origin and (final) destination location, and may conceal additional spatial movements during the elapsed time. The evolutionary history of swIAVs in Thailand is also made more complex by the frequency of reassortment involving multiple clades, poorly supported clusters, and singletons. Whereas avian-like Eurasian viruses in China are monophyletic and result from a single introduction from Europe (introduction 15), the Thai viruses from this lineage are monophyletic only in the PB2, NP, and N1 trees. Our estimate of two introductions of avian-like Eurasian viruses into Thailand is therefore conservative and likely underestimates the true number (introductions 13 and 14, <xref rid="F2" ref-type="fig">Figs. 2</xref> and <xref rid="F3" ref-type="fig">3</xref>, <xref ref-type="supplementary-material" rid="SD1">Supplementary Figs. 1–7</xref>). Singletons and unsupported clusters also were observed among South Korean viruses, complicating the estimation of the number of viral migration events into South Korea as well.</p></sec><sec id="S6"><title>Importance of live swine trade in the global dispersal of swIAVs</title><p id="P13">We used a generalized linear model (GLM) extension of phylogeographic inference<sup><xref rid="R21" ref-type="bibr">21</xref></sup> to identify the putative drivers of swIAV migration events inferred from the genetic data. This approach considers single introductions and clusters that may have poor resolution, the uncertainty of which is accommodated by the analysis. The Bayesian model averaging approach found consistent and strong evidence that the amount of asymmetric live swine trade (measured in USD for the years 1996–2012) from one country to another is the dominant driver of the dispersal of IAVs in swine globally. This is reflected by the maximal estimated inclusion probability of live swine trade for all six internal gene segments (probability of 1, results for PB2, PB1, PA, and NP are presented in <xref rid="F5" ref-type="fig">Fig. 5</xref>; results for the shorter MP and NS segments are available in <xref ref-type="supplementary-material" rid="SD1">Supplementary Fig. 8</xref>; the analysis could not be performed on the HA or NA due to the high frequency of human viruses in these phylogenies). Accordingly, viral migration, measured by ‘Markov jump’ counts, was positively correlated with live swine trade volume (USD) (rho = 0.52, p = 1.5 × 10<sup>−7</sup>, Spearman correlation, <xref ref-type="supplementary-material" rid="SD1">Supplementary Fig. 9</xref>).</p><p id="P14">Other potential predictors, including swine population size, contributed little to the observed global spatial patterns, except for the number of sequences of the country of destination (average probability of ~0.5 across the six internal segments). Not only is the effect of live swine trade consistently robust to the inclusion of sample size, the high inclusion probability of live swine trade in cases where the effect of sample size is particularly low (e.g., the PA segment) indicates that it is also independent of sample size. The conditional effect size (the size of the effect conditional on the predictor being included in the model) ranges between about 3 to 5 on a log scale, implying that viral lineage movement probability is several orders of magnitude higher for connections with the highest swine trade compared to connections without trade (<xref rid="F5" ref-type="fig">Fig. 5</xref>).</p><p id="P15">Our GLM analysis could be affected by regional differences in the early establishment of swIAVs, with swIAVs having been potentially seeded later in Asia and decreasing the probability of viral export from Asia. To further explore this hypothesis, we conducted a sensitivity analysis focused on two more recent periods that correspond approximately to the emergence in China of classical North American swIAVs (1990–2013) and avian-origin swIAVs from Europe (2000–2013), using an ‘epoch’ extension of the diffusion model<sup><xref rid="R22" ref-type="bibr">22</xref></sup>. We find that the volume of live swine is still the only well-supported predictor of viral migration for both periods (<xref ref-type="supplementary-material" rid="SD1">Supplementary Table 4</xref>), including when swIAVs are thought to be endemic at high levels in Asia as well as Europe and North America. In addition to the relatively low number of highly supported migration events during 2000–2013 (<xref rid="F2" ref-type="fig">Fig. 2</xref>), transitions involving singleton viruses or clades that do not have high bootstrap support also contribute to the signal over this restricted time period, including migrations across the relatively porous US-Canada border and between the US/Canada and South Korea.</p></sec><sec id="S7"><title>Predicted spatial dissemination of swIAVs</title><p id="P16">To explore how the network of global live swine trade may drive movements of swIAVs beyond the ten countries for which whole-genome sequence data was available, we used data on pairwise live swine trade between 146 countries to simulate the patterns of viral dissemination under different epidemiological scenarios. <xref rid="F6" ref-type="fig">Fig. 6</xref> explores the predicted spatiotemporal spread of a new swIAV lineage that hypothetically originates in swine in five countries with large swine populations: Canada, China, France, Mexico, and the United States. These predictions are largely consistent with the spatial movements observed in the genetic data, with a high probability of viral export from the US and Canada into Asia, and from Europe to Asia, whereas epidemics originating from China have low probabilities of onward dissemination to other countries (<xref ref-type="supplementary-material" rid="SD1">Supplementary Table 5</xref>). In addition, we identified predicted connections with countries not sampled in our study, including swIAV migration from the USA and Canada to many countries in Latin America, as well as Russia, Kazakhstan, Malaysia, Thailand, and Singapore. Interestingly, our model suggests that a virus seeded in Mexican swine is comparatively less likely to disseminate to swine in other countries.</p><p id="P17">We also used our model to estimate the probability of co-invasion of European and North American swIAVs lineages, illustrating the potential for reassortment between lineages of European and North American descents, of the kind that generated the 2009 pH1N1 virus. Overall, co-invasion is strongly regionalized, with the highest probability in East and South-East Asian countries, particularly China, South Korea, and Russia (<xref rid="F6" ref-type="fig">Fig. 6</xref>, <xref ref-type="supplementary-material" rid="SD1">Supplementary Table 5</xref>). Conversely, South Asia, the Middle East, Africa, and Australia exhibited a low probability of invasion by each of these lineages. Mid-level probabilities were found in regions with a high probability of invasion by only one of the two lineages (i.e., the North American lineage in the Americas, and the Eurasian lineage in Europe). Interestingly, these simulations reveal a low probability of co-invasion in Mexico, where pH1N1 first emerged in humans, owing to the low probability of invasion by a European swine virus in Latin America.</p></sec></sec><sec sec-type="discussion" id="S8"><title>Discussion</title><p id="P18">The unknown origins of the swine virus that begot the 2009 H1N1 pandemic underscores the importance of understanding how influenza A viruses evolve in swine at a global scale, including regions where swIAV surveillance is lacking. Our expansive phylogenetic analysis of global swIAV sequence data demonstrates the importance of the asymmetrical nature of the global live swine trade on the global ecology and evolution of swIAVs. Using a phylogeographic GLM approach to assess the strength of specific predictors, we determine that the size of a country’s swine population is not a major factor in the rate of viral export to other countries. As a notable case in point, China, which hosts the world’s largest swine population, does not appear to be a major source of the viral diversity observed in other Asian countries. Rather, Japan, Thailand, Vietnam, and South Korea independently imported novel viruses from Europe and North America, most likely via long-distance live swine trade.</p><p id="P19">The reported pattern of swIAV dissemination is a reverse of a model proposed for the global spread of A/H3N2 seasonal influenza viruses in humans, in which a highly connected network of South-East Asian countries, including China, acts as a key source of viruses for Europe, North America, and other continents<sup><xref rid="R21" ref-type="bibr">21</xref>,<xref rid="R23" ref-type="bibr">23</xref></sup>, reflecting differences in the disease and mobility patterns of humans and swine. These findings have important implications for swIAV surveillance strategies, as the relatively low levels of viral gene flow between Asian countries means that no single country in Asia can serve as a proxy for the region, including China’s large swine herds. The extent of viral genetic diversity in Thailand highlights the importance of enhancing surveillance throughout South-East Asia, including countries not sampled in our study, such as Malaysia, Indonesia, Singapore, Laos, and Cambodia, as well as undersampled countries such as Thailand, Vietnam, Japan, and South Korea. Furthermore, Russia emerged as a hotspot for invasion and co-invasion of divergent lineages in our simulations, and yet Russia has no publicly available whole-genome swIAV sequences.</p><p id="P20">The limited number of sequences from Asian countries other than China (particularly via Hong Kong, the final destination of large numbers of hogs from mainland China) reduces our ability to detect viral migration events within Asia, particularly those that do not transmit onward in swine for many years. However, the high number of viruses identified in Asia that were of North American and European origin indicates that sample bias alone cannot explain the lack of viral exchange observed between Asian countries. Analysis of larger, less-constrained data sets including all available HA swIAV sequences from Asia identified several additional viral migration events from Europe and North American swine into Asia, but only limited evidence for one additional putative connection between two Asian countries (<xref ref-type="supplementary-material" rid="SD1">Supplementary Figs. 10–11</xref>). However, all inferences of spatial connections must be interpreted within the context of the many countries that are unsampled and under-sampled, and long branches may conceal additional spatial movements between the origin and (final) destination location.</p><p id="P21">It is important to note that our study focused only on the international dissemination of swIAVs, and did not consider the probability of initial emergence of an epidemic within a country, which is likely to be influenced by numerous local factors related to national swine farming practices, including the size and density of farms, movements of pigs within countries, and opportunities for interspecies transmission. As demonstrated previously, the dynamics of outbreaks within a large country like the US can be complex, with different regions acting as source and sink populations for viral diversity<sup><xref rid="R17" ref-type="bibr">17</xref></sup>. Novel IAVs of human origin have emerged repeatedly in swine in countries in North America, South America, Asia, and Europe, suggesting that swine populations in these regions can sustain new epidemics<sup><xref rid="R24" ref-type="bibr">24</xref></sup>. The extent of viral export from a country of origin is a product of both the national prevalence of circulating swIAVs and the volume of live swine export. In this study, we were unable to assess whether geographic differences in the prevalence of swIAVs affects large-scale viral migration, as population-level virological and serological data indicative of swIAV prevalence is available only from a limited number of study sites and time periods that are unlikely to be sufficiently representative for a global study. We therefore recognize that there are scenarios where live swine trade alone would not be a good predictor of viral migration. For example, if the major exporters of live swine (North America and Europe) did not have large endemic swIAV populations, then live swine trade alone would not be a good predictor of viral migration. This does not appear to be the case, as North America and Europe have long histories of endemic swIAV circulation and the highest volumes of outgoing international live swine trade. The apparent association between viral endemicity and trade export may not be a coincidence, as the features that enable countries in Europe and North America to export high volumes of live swine (i.e., large-scale commercial swine production) also are likely to be conducive to sustained endemic swIAV transmission. Finally, our analysis did not consider swine influenza vaccine use, which is highly heterogeneous within and between countries, and would be an important, albeit challenging, factor to integrate into future studies that consider the interaction between the national prevalence of swIAVs, the structure of swine industries, and the global spatial dynamics of trade and viral migration.</p><p id="P22">Although it is not possible at this time to incorporate empirical data on historical differences in swIAV prevalence by region, we were able to explore the interaction of influenza virus prevalence and trade volume using our simulation model. Of particular interest was the question of whether the low levels of viral export from Asia observed in our study could be an artifact of historically lower levels of endemic swIAV activity in Asia. As a consequence, viral exports from Asian countries might be expected to rise in the future as swIAVs become endemic at higher levels throughout the region. There was no support for this hypothesis, as our simulations predicted very low rates of viral export from China to other countries even when the prevalence of swIAVs in China is set unrealistically high (e.g., 58% of Chinese swine infected when R<sub>0</sub> = 1.5). Our simulations predicted low rates of swIAV export from Japan and Thailand under similar transmissibility scenarios (<xref ref-type="supplementary-material" rid="SD1">Supplementary Fig. 12</xref>), with the exception of viral dissemination from Thailand to Cambodia, as these countries are trade partners. These predictions are consistent with our observation from the genetic data that swIAVs could have circulated in Japan and Thailand for many decades without substantial onward dissemination to other Asian locations sampled in our study. In addition, a sensitivity analysis of the GLM model limited to the more recent 1990–2013 and 2000–2013 periods where swIAVs were established in Asia lends further support for the importance of trade. Overall, our findings suggest that viral exchange between Asian countries with low levels of trade is unlikely to increase in the future, regardless of the potential increases in endemic swIAV activity in the region as farming practices are modernized and swine farms become larger.</p><p id="P23">Despite the importance of swine trade in the global ecology of swIAVs, it should be noted that humans may be equally, if not more, important in disseminating IAV diversity to swine herds globally<sup><xref rid="R25" ref-type="bibr">25</xref></sup>. Even in the absence of international swine trade, swIAVs of human origin would likely still circulate in the majority of countries in our study, including in Asia. Quarantine and other restrictions in international trade may have the potential to reduce the genetic diversity of swIAVs, but are not likely to prevent swIAVs from circulating in a country’s swine population (Australia is a case in point<sup><xref rid="R26" ref-type="bibr">26</xref></sup>). The frequency of human-to-swine transmission has been even more apparent since the 2009 H1N1 pandemic, and humans have disseminated pH1N1 viruses to swine in numerous countries that had not previously reported IAV activity in swine, including Australia<sup><xref rid="R26" ref-type="bibr">26</xref>,<xref rid="R27" ref-type="bibr">27</xref></sup>, Brazil<sup><xref rid="R28" ref-type="bibr">28</xref>,<xref rid="R29" ref-type="bibr">29</xref></sup>, India<sup><xref rid="R30" ref-type="bibr">30</xref></sup>, Cameroon<sup><xref rid="R31" ref-type="bibr">31</xref></sup>, Mexico<sup><xref rid="R32" ref-type="bibr">32</xref></sup>, Nigeria<sup><xref rid="R33" ref-type="bibr">33</xref></sup>, Sri Lanka<sup><xref rid="R34" ref-type="bibr">34</xref></sup>, and several countries in Europe<sup><xref rid="R35" ref-type="bibr">35</xref>–<xref rid="R37" ref-type="bibr">37</xref></sup>.</p><p id="P24">Unfortunately, these new data did not advance our understanding of the evolution of the pH1N1 virus during the many years of undetected circulation in swine prior to 2009. Given that the human pandemic likely emerged in Mexico<sup><xref rid="R3" ref-type="bibr">3</xref></sup>, the most parsimonious explanation is that the pH1N1 virus transmitted from swine to humans in Mexico or a nearby locality. However, extremely little swIAV sequence data is available from swine in Mexico and other parts of Latin America, and Eurasian viruses have not been detected in any part of the Americas to date. Our simulation model provides a quantitative indicator of where the reassortment event that produced the pH1N1 virus in swine was most likely to have occurred, based on the probability of invasion with both North American and European viruses. Given limitations in our model, including the lack of information on within-country dynamics and the likelihood of initial viral emergence within seed countries, we consider the relative ranking of probabilities to be more important than their absolute values. To date, Asia is the only region where any reassortant viruses containing both North American and Eurasian virus segments have been detected<sup><xref rid="R38" ref-type="bibr">38</xref></sup>, consistent with our simulations, which show a high probability of co-invasion by both North American and Eurasian swine lineages in China, South-East Asia, and Russia. However, it remains unclear how a reassortant virus that most likely emerged in swine in Asia caused its first outbreak in humans in Mexico. Given the lack of south-to-north swine trade flows in the Americas, some swIAV lineages are likely to be exclusive to Latin America<sup><xref rid="R39" ref-type="bibr">39</xref>,<xref rid="R40" ref-type="bibr">40</xref></sup> and do not reach the United States or Canada. Strengthened surveillance in Latin America is warranted to gain a better understanding of swIAV diversity in the region.</p><p id="P25">Finally, we have focused our study on the global dynamics of influenza A viruses in swine, but our findings invite investigation of how trade, quarantine, and swine farming practices affect the spatial dynamics of other globally dispersed swine pathogens, such as PRRSV and the porcine epidemic diarrhea virus (PEDV) that emerged in US swine herds in 2013. Modeling studies rooted in pathogen sequence information, demographics, and mobility data have the power to inform global surveillance and control strategies for major animal and human disease threats.</p></sec><sec sec-type="materials|methods" id="S9"><title>Materials and Methods</title><sec id="S10"><title>Influenza virus sample preparation</title><p id="P26">Influenza A virus samples collected from swine via routine diagnostic submissions for the years 2002–2011 were randomly selected from the existing influenza virus archive at the University of Minnesota Veterinary Diagnostic Laboratory (UMVDL). These samples were chosen to best represent this time period and the three main geographical regions of US hog production: the Southeast region (US states of Alabama, Georgia, Kentucky, North Carolina, Tennessee, Virginia), South-central/west region (Arkansas, Colorado, New Mexico, Oklahoma, Texas), and Midwest region (Iowa, Illinois, Indiana, Kansas, Minnesota, Missouri, Nebraska). Samples from swine in Canada (2005–2011) and Mexico (2010–2011) also were selected from UMVDL, as available. Original specimen material (nasal swab supernatant or lung tissue homogenate stored at −80°C) were aliquotted from the archived samples and sent to the J. Craig Venter Institute (JCVI) in Rockville, MD for sequencing.</p></sec><sec id="S11"><title>Influenza virus genome sequencing</title><p id="P27">The complete genomes of 240 influenza viruses collected from North American swine were sequenced at JCVI. Viral RNA was isolated using the ZR 96 Viral RNA kit (Zymo Research Corporation, Irvine, CA, USA). The influenza A genomic RNA segments were simultaneously amplified from 3 μL of purified RNA using a multi-segment RT-PCR strategy<sup><xref rid="R41" ref-type="bibr">41</xref></sup>. The M-RTPCR amplicons were sequenced using Nextera Library construction using the MiSeq platform (Illumina, Inc., San Diego, CA, USA). Additionally, M-RTPCR amplicons were sheared for 7 minutes and Ion Torrent compatible barcoded adapters were ligated to create 200 base-pair libraries that were purified and sequenced using Ion Torrent (Life Technologies, Grand Island, NY, USA). All data sequenced for this study was submitted to the Influenza Virus Resource at the National Center for Biotechnology Information’s GenBank<sup><xref rid="R42" ref-type="bibr">42</xref></sup>, and accession codes are available in <xref ref-type="supplementary-material" rid="SD1">Supplementary Table 3</xref>.</p></sec><sec id="S12"><title>Phylogenetic analysis</title><p id="P28">In addition to the sequences generated for this study, whole-genome sequences from influenza A viruses collected in swine globally during 1960–2013 were downloaded from the Influenza Virus Resource at GenBank<sup><xref rid="R42" ref-type="bibr">42</xref></sup>. Viruses were removed that (a) had truncated sequences, (b) were of avian origin with no evidence of circulation in swine, (c) had unknown geographic origin, or (e) had evidence of lab errors (assessed by root-to-tip divergence using the program Path-O-Gen v1.3). Due to the disproportionately large number of swIAV sequences from the United States for the years 2009–2013, 100 of these were randomly sub-sampled.</p><p id="P29">Sequence alignments were constructed for each of the six internal gene segments (PB2, PB1, PA, NP, MP, and NS) and for the H1, H3, N1, and N2 antigenic segments separately using MUSCLE v3.8.31<sup><xref rid="R43" ref-type="bibr">43</xref></sup>, with manual correction in Se-Al v2.0 (available at <ext-link ext-link-type="uri" xlink:href="http://tree.bio.ed.ac.uk/software/seal/">http://tree.bio.ed.ac.uk/software/seal/</ext-link>). Phylogenetic trees were inferred using the neighbor-joining method available in PAUP v4.0b10 for each of the ten alignments (available at <ext-link ext-link-type="uri" xlink:href="http://paup.csit.fsu.edu/">http://paup.csit.fsu.edu/</ext-link>). For the PB2, PB1, and PA segments, each virus was categorized as belonging to one of the following lineages: (a) classical swine virus lineage, (b) triple reassortant (‘trig’) lineage, (c) avian-origin Eurasian swine lineage, (d) the pandemic H1N1 lineage that emerged in humans in 2009 and transmitted from humans to swine globally during 2009–2013 or (e) related to human seasonal influenza A viruses. For the NP, MP, and NS segments, each virus was categorized as (a) classical, (b) avian-origin Eurasian, (c) pandemic, or (d) human seasonal. No ‘trig’ category exists for the NP, MP, and NS segments because triple reassortant viruses contain classical virus NP, MP, and NS segments acquired through reassortment. For the H1 and N1 segments, each virus was categorized as (a) classical, (b) avian-origin Eurasian, (c) pandemic, or (d) human seasonal virus origin. All H3 and N2 segments belonged to the same category: human seasonal H3N2 virus-related.</p><p id="P30">The sequence alignment for each segment was further divided into each of these lineages. For the PB2, PB1, PA, NP, MP, NS, and N1 segments, very few (< 10) viruses were found to be of recent human seasonal virus origin, which is consistent with previous findings<sup><xref rid="R24" ref-type="bibr">24</xref></sup>, so these data sets were excluded from further analyses. The pandemic H1N1 viruses that have been recently transmitted from humans to swine since 2009 remain spatially structured by country (or continent), with little evidence of long-distance migration between continents<sup><xref rid="R44" ref-type="bibr">44</xref></sup>, and therefore sequences of pH1N1 origin were not included in further analyses. Similarly, no global migration was observed among N1 swIAV sequences that were closely related to human seasonal H1N1 viruses, and they were excluded from the study. In total, 22 segment- and lineage-specific data sets were included in the analysis (<xref ref-type="supplementary-material" rid="SD1">Supplementary Table 6; Supplementary Data 1</xref>). For the H3, N2, and H1 human-like (delta) lineages, human seasonal influenza virus H3, N2, and H1 sequences also were included as background. To reduce the impact of sample bias, additional phylogenies were inferred using all available full-length swIAV H1 and H3 sequence data from Asia, which included an additional 206 swIAVs from China, South Korea, Thailand, and Japan for the classical H1 segment, 191 swIAVs from China for the avian-origin Eurasian H1 segment, and 230 swIAVs from China, Thailand, Japan, Mongolia, South Korea, and Indonesia for the human seasonal virus origin H3 segment.</p><p id="P31">Phylogenetic relationships were inferred for each of the 22 data sets separately using the time-scaled Bayesian approach using MCMC available via the BEAST v1.8.0 package<sup><xref rid="R45" ref-type="bibr">45</xref></sup> and the high-performance computational capabilities of the Biowulf Linux cluster at the National Institutes of Health, Bethesda, MD (<ext-link ext-link-type="uri" xlink:href="http://biowulf.nih.gov">http://biowulf.nih.gov</ext-link>). A relaxed uncorrelated lognormal (UCLN) molecular clock was used, with a flexible Bayesian skyline plot (BSP) demographic model (10 piece-wise constant groups), and a general-time reversible (GTR) model of nucleotide substitution with gamma-distributed rate variation among sites. For viruses for which only the year of viral collection was available, the lack of tip date precision was accommodated by sampling uniformly across a one-year window from January 1<sup>st</sup> to December 31<sup>st</sup>. The MCMC chain was run separately three times for each of the data sets for at least 100 million iterations with sub-sampling every 10,000 iterations, using the BEAGLE library to improve computational performance<sup><xref rid="R46" ref-type="bibr">46</xref></sup>. All parameters reached convergence, as assessed visually using Tracer v.1.6, with statistical uncertainty reflected in values of the 95% highest posterior density (HPD). At least 10% of the chain was removed as burn-in, and runs for the same lineage and segment were combined using LogCombiner v1.8.0 and downsampled to generate a final posterior distribution of 1,000 trees that was used in subsequent analyses.</p><p id="P32">The phylogeographic analysis considered ten locations: Argentina, Canada, China, Europe, Japan, Mexico, Thailand, United States of America, South Korea, and Vietnam. All viruses from Europe (our study included swIAV data from Belgium, Czech Republic, Denmark, France, Germany, Italy, Netherlands, Poland, Spain, and the United Kingdom) were categorized into a single spatial category due to the high level of influenza virus mixing within Europe (<xref ref-type="supplementary-material" rid="SD1">Supplementary Fig. 13</xref>) and the relatively low level of sampling of individual countries. Similarly, we considered Hong Kong viruses to be part of China based on genetic similarities. The location state was specified for each viral sequence, allowing the expected number of location state transitions in the ancestral history conditional on the data observed at the tree tips to be estimated using ‘Markov jump’ counts<sup><xref rid="R20" ref-type="bibr">20</xref></sup>, which provided a quantitative measure of asymmetry in gene flow between regions (a representative XML file used in the analysis is provided in <xref ref-type="supplementary-material" rid="SD1">Supplementary Data 2</xref>). For computational efficiency the phylogeographic analysis was run using an empirical distribution of 1,000 trees<sup><xref rid="R21" ref-type="bibr">21</xref></sup>, allowing the MCMC chain to be run for 25 million iterations, sampling every 1,000. A Bayesian stochastic search variable selection (BSSVS) was employed to improve statistical efficiency for all data sets containing more than four location states. Maximum clade credibility (MCC) trees were summarized using TreeAnnotator v1.8.0 and the trees were visualized in FigTree v1.4.2.</p></sec><sec id="S13"><title>Testing predictors of global swIAV migration</title><p id="P33">A generalized linear model (GLM)<sup><xref rid="R21" ref-type="bibr">21</xref></sup> parameterization of the discrete phylogeographic diffusion model was employed to estimate the contribution of potential predictors to the migration patterns of swIAVs (a representative XML file used in the analysis is provided in <xref ref-type="supplementary-material" rid="SD1">Supplementary Data 3</xref>; the R code used to summarize the estimates is provided in <xref ref-type="supplementary-material" rid="SD1">Supplementary Data 4</xref>). First, the trade value (USD) for live swine trade between countries (asymmetric) for the years 1996–2012 was obtained from the United Nations’ Commodity Trade Statistics Database (available at <ext-link ext-link-type="uri" xlink:href="http://comtrade.un.org">http://comtrade.un.org</ext-link>, accessed March 20, 2014) (<xref ref-type="supplementary-material" rid="SD1">Supplementary Data 5</xref>). For the purposes of our study, we calculated the total trade value for each country. Data from all European countries was aggregated within the category ‘Europe’, and data from mainland China, Macao SAR, and Hong Kong SAR was aggregated within the category ‘China’. Second, estimates of the number of live swine by country and the total number of live swine imports and exports per country were obtained for a longer time period (1969–2010) from the Food and Agriculture Organization (FAO) of the United Nations Datasets repository (available at <ext-link ext-link-type="uri" xlink:href="http://data.fao.org/datasets">http://data.fao.org/datasets</ext-link>, accessed March 21, 2014) (<xref ref-type="supplementary-material" rid="SD1">Supplementary Table 7</xref>). Again, data were aggregated across years and for Europe as well as for China. Although there is variance in trade volumes between years, some of which may also reflect variance in reporting, consistent differences in trade volume were evident among countries across years (<xref ref-type="supplementary-material" rid="SD1">Supplementary Fig. 14</xref>). Further, we estimated the country-specific percentage change in the pig population over the study period (ratio of the numbers of live swine in 1969 versus 2010) as an additional putative predictor of swIAVs migration. All the predictors were log-transformed and standardized prior to their specification in the GLM parameterization. We performed the GLM analysis separately for each of the six internal gene segments (PB2, PB1, PA, NP, MP, and NS) and jointly for all three swIAV lineages (avian-origin Eurasian, triple reassortant, and classical) for each segment. We achieve this, for each segment, by sharing a single GLM-diffusion model across the independent evolutionary histories of the three viral lineages. For the NP, MP, and NS segments only two lineages were included, as the triple reassortant lineage is an extension of the classical lineage for these three segments. We also excluded Argentina from the GLM analysis because no viral migration was observed between Argentina and any other country in our study. Additionally, to explore the effect of regional differences in the early establishment of swIAVs we used an ‘epoch’ extension of the diffusion model<sup><xref rid="R22" ref-type="bibr">22</xref></sup> for two periods that correspond approximately to the emergence in China of classical North American swIAVs (1990–2013) and avian-origin swIAVs from Europe (2000–2013) (representative XML file available provided in <xref ref-type="supplementary-material" rid="SD1">Supplementary Data 6</xref>).</p></sec><sec id="S14"><title>Meta-population model simulating the global spread of swine influenza</title><p id="P34">Next, we built a meta-population model to explore the global spread of swIAVs and identify potential geographical hotspots for reassortment between viruses originating from different regions, which may generate novel viruses with pandemic potential. We built a meta-population model to explore the global spread of swIAVs and identify potential geographical hotspots for reassortment. We employed a stochastic patch-based SIR model adapted from an earlier model for the global diffusion of human influenza<sup><xref rid="R47" ref-type="bibr">47</xref></sup>. Each patch represents the swine population of an individual country, and patches are linked based on live swine trade movements. To calibrate the model, we obtained swine population sizes and pairwise between-country trade information for 146 countries reporting to FAO during 1969–2010 (<ext-link ext-link-type="uri" xlink:href="http://data.fao.org/datasets">http://data.fao.org/datasets</ext-link>; <xref ref-type="supplementary-material" rid="SD1">Supplementary Table 7, Supplementary Data 5</xref>), which coincides with the study period considered for phylogenetic analysis of swIAVs. We used the averages of swine population sizes and pairwise trade volumes throughout the study period for model simulations.</p><p id="P35">The influenza simulation model is as follows. Let <bold>S</bold>, <bold>I</bold>, and <bold>R</bold> denote vectors representing the number of susceptible, infected, and recovered swine at any time point in each of the 146 countries studied. We let <italic>μ</italic> = 1/5 denote the daily probability that an individual recovers (so that the infectious period is 5 days<sup><xref rid="R48" ref-type="bibr">48</xref>,<xref rid="R49" ref-type="bibr">49</xref></sup>), and <italic>β</italic><bold>I</bold> the daily per-capita rate of infection from infectious individuals within the same country. Here, <italic>β</italic> varies between countries; it is a vector scaled such that the effective reproduction number <italic>R</italic><sub>eff</sub> = <italic>β</italic><bold>N</bold> / <italic>μ</italic> is the same in all countries, where N is the vector of swine population sizes. We use an <italic>R<sub>0</sub></italic> of 1.5 in main analyses, consistent with limited information on swine influenza dynamics<sup><xref rid="R48" ref-type="bibr">48</xref>,<xref rid="R49" ref-type="bibr">49</xref></sup> The per-capita rate of contacts with infectious swine from other countries is given by <bold>G</bold>*<bold>I</bold>, where <bold>G</bold> represents the 146*146 matrix of between-country coupling.</p><p id="P36">To build <bold>G</bold>, we first create <bold>T</bold>, a 146*146 matrix with off-diagonal elements based on empirical live swine trade and zeros along the diagonal. We then rescale <bold>T</bold> by the estimated trade coefficients of the phylogeographic GLM model, as the GLM model suggests the relationship between trade and viral migration is not linear (but results are qualitatively similar with no rescaling). Following past work<sup><xref rid="R47" ref-type="bibr">47</xref></sup>, the rescaled matrix <bold>T’</bold> is then tuned by a free parameter, <italic>c</italic>, which governs the amount of international vs. domestic contacts between swine, while at the same time allowing conversion of empirical swine trade data (provided in $ amount by FAO) into actual population movements. Tuning parameter <italic>c</italic> allows obtaining realistic time course of infection, with global epidemics lasting between several months to several years, in line with (limited data available on) the global spread of past swine outbreaks. Our final coupling matrix is <bold>G</bold> = <italic>c</italic><bold>G</bold>’, where <italic>c</italic> is such that the maximum element of the coupling matrix is lower than 10<sup>−3</sup>.</p><p id="P37">We use a spatially-extended chain-binomial system to update the progression of the epidemic in each patch<sup><xref rid="R47" ref-type="bibr">47</xref></sup>:
<disp-formula id="FD1"><mml:math id="M1" display="block" overflow="scroll"><mml:mtable columnalign="left"><mml:mtr columnalign="left"><mml:mtd columnalign="left"><mml:mrow><mml:msub><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>S</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mstyle></mml:msub><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mo>=</mml:mo></mml:mstyle><mml:msub><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>S</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>t</mml:mi></mml:mstyle></mml:msub><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mo>–</mml:mo></mml:mstyle><mml:msub><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>W</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>t</mml:mi></mml:mstyle></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr columnalign="left"><mml:mtd columnalign="left"><mml:mrow><mml:msub><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>I</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mstyle></mml:msub><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mo> </mml:mo><mml:mo>=</mml:mo></mml:mstyle><mml:msub><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>I</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>t</mml:mi></mml:mstyle></mml:msub><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mo>+</mml:mo></mml:mstyle><mml:msub><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>W</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>t</mml:mi></mml:mstyle></mml:msub><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mo> </mml:mo><mml:mo>-</mml:mo></mml:mstyle><mml:msub><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>V</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>t</mml:mi></mml:mstyle></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr columnalign="left"><mml:mtd columnalign="left"><mml:mrow><mml:msub><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>R</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mstyle></mml:msub><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mo> </mml:mo><mml:mo>=</mml:mo></mml:mstyle><mml:msub><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>R</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>t</mml:mi></mml:mstyle></mml:msub><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mo>+</mml:mo></mml:mstyle><mml:msub><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>V</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>t</mml:mi></mml:mstyle></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr columnalign="left"><mml:mtd columnalign="left"><mml:mrow><mml:msub><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>W</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>t</mml:mi></mml:mstyle></mml:msub><mml:mo> </mml:mo><mml:mo> </mml:mo><mml:mo>=</mml:mo><mml:mtext> binomial</mml:mtext><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>S</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>t</mml:mi></mml:mstyle></mml:msub><mml:mo>,</mml:mo><mml:mspace width="0.2em"></mml:mspace><mml:mn>1</mml:mn><mml:mo>–</mml:mo><mml:mtext>e^</mml:mtext><mml:mo>−</mml:mo><mml:mi>β</mml:mi><mml:mo>∗</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>I</mml:mi><mml:mo>+</mml:mo><mml:mi>G</mml:mi><mml:mo>∗</mml:mo><mml:mi>I</mml:mi></mml:mstyle><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr columnalign="left"><mml:mtd columnalign="left"><mml:mrow><mml:msub><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>V</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>t</mml:mi></mml:mstyle></mml:msub><mml:mo> </mml:mo><mml:mo> </mml:mo><mml:mo>=</mml:mo><mml:mtext> binomial</mml:mtext><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>I</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>t</mml:mi></mml:mstyle></mml:msub><mml:mo>,</mml:mo><mml:mspace width="0.2em"></mml:mspace><mml:mi>μ</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>Here, <bold>W<sub>t</sub></bold> is the daily incidence of swine influenza (ie the number of new cases), and <bold>V<sub>t</sub></bold> is the number of new ‘recovereds’.</p><p id="P38">In simulations, the epidemic is initialized by infecting 5 swine in a predetermined seed country; we explored various scenarios with seeds in countries of the Americas, Europe and Asia. After the first infection occurs in each country, we draw a multinomial to determine the source of the infection, from a normalized vector <bold>G*I</bold>. For each scenario, involving a given source country (e.g., US, Canada, UK, France, China), we ran 1,000 simulations, allowed to run over a three year time period, and assess the probability of swine flu invasion in each non-source country, and its most likely source of infection. We conducted sensitivity analyses with higher and lower values of the free parameter <italic>c</italic> and <italic>R<sub>0</sub></italic>, which mostly affected the time course of the global epidemic and the synchronicity of epidemics across locations, but did not change dramatically the identification of hotspots countries for onward spread. To explore the probability of reassortment between viruses originating from North America and Europe, we ran 1000 independent simulations of epidemics starting in North America (USA and Canada), and in each the 10 European countries with largest swine populations, and compute the co-invasion probability in country <italic>i</italic>, following:
<disp-formula id="FD2"><mml:math id="M2" display="block" overflow="scroll"><mml:mi>P</mml:mi><mml:mo>.</mml:mo><mml:msub><mml:mi mathvariant="italic">coinv</mml:mi><mml:mrow><mml:mo>.</mml:mo><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mo stretchy="false">[</mml:mo><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:msub><mml:mi mathvariant="normal">Π</mml:mi><mml:mrow><mml:mi>j</mml:mi><mml:mi mathvariant="italic"> source</mml:mi><mml:mi mathvariant="italic"> countries </mml:mi><mml:mi mathvariant="italic">in</mml:mi><mml:mi mathvariant="italic"> North</mml:mi><mml:mi mathvariant="italic"> America</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>p</mml:mi><mml:mo>.</mml:mo><mml:msub><mml:mi mathvariant="italic">inv</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">]</mml:mo><mml:mo>∗</mml:mo><mml:mo stretchy="false">[</mml:mo><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:msub><mml:mi mathvariant="normal">Π</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mi mathvariant="italic"> source</mml:mi><mml:mi mathvariant="italic"> countries</mml:mi><mml:mi mathvariant="italic"> in</mml:mi><mml:mi mathvariant="italic"> Europe</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>p</mml:mi><mml:mo>.</mml:mo><mml:mi>i</mml:mi><mml:mi>n</mml:mi><mml:msub><mml:mi>v</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">]</mml:mo><mml:mo>,</mml:mo></mml:math></disp-formula>where <italic>p.inv<sub>i,k</sub></italic> is the probability that country <italic>i</italic> is invaded when the outbreak is seeded in country <italic>k</italic>.</p><p id="P39">Estimates were then used to produce risk maps for swIAV invasion and co-invasion using software available in the R package (<ext-link ext-link-type="uri" xlink:href="http://www.r-project.org/">http://www.r-project.org/</ext-link>).</p></sec></sec><sec sec-type="supplementary-material" id="S15"><title>Supplementary Material</title><supplementary-material content-type="local-data" id="SD1"><label>1</label><media xlink:href="NIHMS666251-supplement-1.pdf" orientation="portrait" xlink:type="simple" id="d37e1186" position="anchor"></media></supplementary-material></sec></body><back><ack id="S16"><p>The authors wish to thank Drs. Sagar Goyal and Devi Patnayak, along with Wendy Wiese, Becca Wheeldon, and Lotus Smasal at UMVDL for their virology work, and Rebecca Halpin, Nadia Fedorova, Xudong Lin, and Timothy B. Stockwell at the J. Craig Venter Institute for their assistance in viral sequencing. This manuscript was prepared while DEW was employed at the J. Craig Venter Institute. Partial support for this work was provided by a grant from the Saskatchewan Ministry of Agriculture contract number C70077. This work was supported in part by the Multinational Influenza Seasonal Mortality Study (MISMS), an on-going international collaborative effort to understand influenza epidemiology and evolution, led by the Fogarty International Center, National Institutes of Health, with funding from Office of Pandemics and Emerging Threats at the U.S. Department of Health and Human Services [AR, MIN, CV]. This project has also been funded in part with federal funds from the National Institute of Allergy and Infectious Diseases, National Institutes of Health, Department of Health and Human Services under the [contract number HHSN272200900007C]. PL, MAS, and AR acknowledge funding from the European Union Seventh Framework Programme [FP7/2007-2013] under Grant Agreement no. 278433-PREDEMICS and ERC Grant agreement no. 260864. ECH is supported by an NHMRC Australia Fellowship (AF30).</p></ack><fn-group><fn id="FN3" fn-type="con"><p><bold>Author contributions</bold></p><p>Generated the data used in the analysis: MRC SED DEW. Conceived, designed, and performed the analysis: MIN CV PL. Analyzed the data: MIN CV ECH PL. Wrote the paper: MIN CV ALV MRC SED DEW AR MAS ECH PL. The opinions expressed in this article are the author’s own and do not reflect the view of the Centers for Disease Control and Prevention, the Department of Health and Human Services, or the United States government.</p></fn><fn id="FN4" fn-type="conflict"><p><bold>Competing interests</bold></p><p>The authors have declared that no competing interests exist.</p></fn><fn id="FN5"><p><bold>Accession codes</bold></p><p>All sequences from swine influenza A viruses that were generated and used in this study have been deposited in the GenBank database with the accession codes AHB20822 to AHB20891, AHB20919 to AHB23256, AHB23664 to AHB23860, AHB24418 to AHB24833, and AHB82076 to AHB82086.</p></fn></fn-group><ref-list><ref id="R1"><label>1</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Garten</surname><given-names>RJ</given-names></name><etal></etal></person-group><article-title>Antigenic and genetic characteristics of swine-origin 2009 A(H1N1) influenza viruses circulating in humans</article-title><source>Science</source><volume>325</volume><fpage>197</fpage><lpage>201</lpage><year>2009</year><pub-id pub-id-type="pmid">19465683</pub-id></element-citation></ref><ref id="R2"><label>2</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Smith</surname><given-names>GJD</given-names></name><etal></etal></person-group><article-title>Origins and evolutionary genomics of the 2009 swine-origin H1N1 influenza A epidemic</article-title><source>Nature</source><volume>459</volume><fpage>1122</fpage><lpage>1125</lpage><year>2009</year><pub-id pub-id-type="pmid">19516283</pub-id></element-citation></ref><ref id="R3"><label>3</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lemey</surname><given-names>P</given-names></name><name><surname>Suchard</surname><given-names>M</given-names></name><name><surname>Rambaut</surname><given-names>A</given-names></name></person-group><article-title>Reconstructing the initial global spread of a human influenza pandemic: A Bayesian spatial-temporal model for the global spread of H1N1pdm</article-title><source>PLoS Curr</source><volume>1</volume><fpage>RRN1031</fpage><year>2009</year><pub-id pub-id-type="pmid">20029613</pub-id></element-citation></ref><ref id="R4"><label>4</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Koen</surname><given-names>JS</given-names></name></person-group><article-title>A practical method for field diagnosis of swine diseases</article-title><source>Am J Vet Med</source><volume>14</volume><fpage>468</fpage><lpage>470</lpage><year>1919</year></element-citation></ref><ref id="R5"><label>5</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhou</surname><given-names>NN</given-names></name><etal></etal></person-group><article-title>Genetic reassortment of avian, swine, and human influenza A viruses in American pigs</article-title><source>J Virol</source><volume>73</volume><fpage>8851</fpage><lpage>8856</lpage><year>1999</year><pub-id pub-id-type="pmid">10482643</pub-id></element-citation></ref><ref id="R6"><label>6</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Karasin</surname><given-names>A</given-names></name><name><surname>Carman</surname><given-names>S</given-names></name><name><surname>Olsen</surname><given-names>C</given-names></name></person-group><article-title>Identification of human H1N2 and human-swine reassortant H1N2 and H1N1 influenza A viruses among pigs in Ontario, Canada (2003 to 2005)</article-title><source>J Clin Microbiol</source><volume>44</volume><fpage>1123</fpage><lpage>1126</lpage><year>2006</year><pub-id pub-id-type="pmid">16517910</pub-id></element-citation></ref><ref id="R7"><label>7</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Vincent</surname><given-names>AL</given-names></name><etal></etal></person-group><article-title>Characterization of a newly emerged genetic cluster of H1N1 and H1N2 swine influenza virus in the United States</article-title><source>Virus Genes</source><volume>39</volume><fpage>176</fpage><lpage>85</lpage><year>2009</year><pub-id pub-id-type="pmid">19597980</pub-id></element-citation></ref><ref id="R8"><label>8</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ottis</surname><given-names>K</given-names></name><name><surname>Sidoli</surname><given-names>L</given-names></name><name><surname>Bachmann</surname><given-names>P</given-names></name><name><surname>Webster</surname><given-names>RG</given-names></name><name><surname>Kaplan</surname><given-names>MM</given-names></name></person-group><article-title>Human influenza A viruses in pigs: isolation of a H3N2 strain antigenically related to A/England/42/72 and evidence for continuous circulation of human viruses in the pig population</article-title><source>Arch Virol</source><volume>73</volume><fpage>103</fpage><lpage>108</lpage><year>1982</year><pub-id pub-id-type="pmid">6184032</pub-id></element-citation></ref><ref id="R9"><label>9</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Vijaykrishna</surname><given-names>D</given-names></name><etal></etal></person-group><article-title>Long-term evolution and transmission dynamics of swine influenza A virus</article-title><source>Nature</source><volume>473</volume><fpage>519</fpage><lpage>522</lpage><year>2011</year><pub-id pub-id-type="pmid">21614079</pub-id></element-citation></ref><ref id="R10"><label>10</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Takemae</surname><given-names>N</given-names></name><etal></etal></person-group><article-title>Antigenic variation of H1N1, H1N2 and H3N2 swine influenza viruses in Japan and Vietnam</article-title><source>Arch Virol</source><volume>158</volume><fpage>859</fpage><lpage>876</lpage><year>2013</year><pub-id pub-id-type="pmid">23435952</pub-id></element-citation></ref><ref id="R11"><label>11</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Poonsuk</surname><given-names>S</given-names></name><name><surname>Sangthong</surname><given-names>P</given-names></name><name><surname>Petcharat</surname><given-names>N</given-names></name><name><surname>Lekcharoensuk</surname><given-names>P</given-names></name></person-group><article-title>Genesis and genetic constellations of swine influenza viruses in Thailand</article-title><source>Vet Microbiol</source><volume>167</volume><fpage>314</fpage><lpage>26</lpage><year>2013</year><pub-id pub-id-type="pmid">24095146</pub-id></element-citation></ref><ref id="R12"><label>12</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Song</surname><given-names>D</given-names></name><etal></etal></person-group><article-title>A novel reassortant canine H3N1 influenza virus between pandemic H1N1 and canine H3N2 influenza viruses in Korea</article-title><source>J Gen Virol</source><volume>93</volume><fpage>551</fpage><lpage>554</lpage><year>2012</year><pub-id pub-id-type="pmid">22131311</pub-id></element-citation></ref><ref id="R13"><label>13</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nelson</surname><given-names>MI</given-names></name><etal></etal></person-group><article-title>Genomic reassortment of influenza A virus in North American swine, 1998–2011</article-title><source>J Gen Virol</source><volume>93</volume><fpage>2584</fpage><lpage>2589</lpage><year>2012</year><pub-id pub-id-type="pmid">22993190</pub-id></element-citation></ref><ref id="R14"><label>14</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zell</surname><given-names>R</given-names></name><etal></etal></person-group><article-title>Novel reassortant of swine influenza H1N2 virus in Germany</article-title><source>J Gen Virol</source><volume>89</volume><fpage>271</fpage><lpage>276</lpage><year>2008</year><pub-id pub-id-type="pmid">18089751</pub-id></element-citation></ref><ref id="R15"><label>15</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nelson</surname><given-names>MI</given-names></name><name><surname>Vincent</surname><given-names>AL</given-names></name><name><surname>Kitikoon</surname><given-names>P</given-names></name><name><surname>Holmes</surname><given-names>EC</given-names></name><name><surname>Gramer</surname><given-names>MR</given-names></name></person-group><article-title>Evolution of novel reassortant A/H3N2 influenza viruses in North American swine and humans, 2009–2011</article-title><source>J Virol</source><volume>86</volume><fpage>8872</fpage><lpage>8878</lpage><year>2012</year><pub-id pub-id-type="pmid">22696653</pub-id></element-citation></ref><ref id="R16"><label>16</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ducatez</surname><given-names>MF</given-names></name><etal></etal></person-group><article-title>Multiple reassortment between pandemic (H1N1) 2009 and endemic influenza viruses in pigs, United States</article-title><source>Emerg Infect Dis</source><volume>17</volume><fpage>1624</fpage><lpage>9</lpage><year>2011</year><pub-id pub-id-type="pmid">21892996</pub-id></element-citation></ref><ref id="R17"><label>17</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nelson</surname><given-names>MI</given-names></name><etal></etal></person-group><article-title>Spatial dynamics of human-origin H1 influenza A virus in North American swine</article-title><source>PLoS Path</source><volume>7</volume><fpage>e1002077</fpage><year>2011</year></element-citation></ref><ref id="R18"><label>18</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Shi</surname><given-names>M</given-names></name><etal></etal></person-group><article-title>The spread of type 2 Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) in North America: a phylogeographic approach</article-title><source>Virology</source><volume>447</volume><fpage>146</fpage><lpage>54</lpage><year>2013</year><pub-id pub-id-type="pmid">24210109</pub-id></element-citation></ref><ref id="R19"><label>19</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Anderson</surname><given-names>TK</given-names></name><etal></etal></person-group><article-title>Population dynamics of cocirculating swine influenza A viruses in the United States from 2009 to 2012</article-title><source>Influenza Other Respi Viruses</source><volume>7</volume><issue>Suppl 4</issue><fpage>42</fpage><lpage>51</lpage><year>2013</year></element-citation></ref><ref id="R20"><label>20</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Minin</surname><given-names>VN</given-names></name><name><surname>Suchard</surname><given-names>MA</given-names></name></person-group><article-title>Counting labeled transitions in continuous-time Markov models of evolution</article-title><source>J Math Biol</source><volume>56</volume><fpage>391</fpage><lpage>412</lpage><year>2008</year><pub-id pub-id-type="pmid">17874105</pub-id></element-citation></ref><ref id="R21"><label>21</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lemey</surname><given-names>P</given-names></name><etal></etal></person-group><article-title>Unifying Viral Genetics and Human Transportation Data to Predict the Global Transmission Dynamics of Human Influenza H3N2</article-title><source>PLoS Pathog</source><volume>10</volume><fpage>e1003932</fpage><year>2014</year><pub-id pub-id-type="pmid">24586153</pub-id></element-citation></ref><ref id="R22"><label>22</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bielejec</surname><given-names>F</given-names></name><name><surname>Lemey</surname><given-names>P</given-names></name><name><surname>Baele</surname><given-names>G</given-names></name><name><surname>Rambaut</surname><given-names>A</given-names></name><name><surname>Suchard</surname><given-names>MA</given-names></name></person-group><article-title>Inferring heterogeneous evolutionary processes through time: from sequence substitution to phylogeography</article-title><source>Syst Biol</source><volume>63</volume><fpage>493</fpage><lpage>504</lpage><year>2014</year><pub-id pub-id-type="pmid">24627184</pub-id></element-citation></ref><ref id="R23"><label>23</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Russell</surname><given-names>CA</given-names></name><etal></etal></person-group><article-title>The global circulation of seasonal influenza A (H3N2) viruses</article-title><source>Science</source><volume>320</volume><fpage>340</fpage><lpage>6</lpage><year>2008</year><pub-id pub-id-type="pmid">18420927</pub-id></element-citation></ref><ref id="R24"><label>24</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nelson</surname><given-names>MI</given-names></name><etal></etal></person-group><article-title>Introductions and evolution of human-origin seasonal influenza A viruses in multinational swine populations</article-title><source>J Virol</source><volume>88</volume><fpage>10110</fpage><lpage>9</lpage><year>2014</year><pub-id pub-id-type="pmid">24965467</pub-id></element-citation></ref><ref id="R25"><label>25</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nelson</surname><given-names>MI</given-names></name><name><surname>Vincent</surname><given-names>AL</given-names></name></person-group><article-title>Reverse zoonosis of influenza to swine: new perspectives on the human-animal interface</article-title><source>Trends Microbiol</source><year>2015</year><pub-id pub-id-type="doi">10.1016/j.tim.2014.12.002</pub-id></element-citation></ref><ref id="R26"><label>26</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Holyoake</surname><given-names>PK</given-names></name><etal></etal></person-group><article-title>The first identified case of pandemic H1N1 influenza in pigs in Australia</article-title><source>Aust Vet J</source><volume>89</volume><fpage>427</fpage><lpage>431</lpage><year>2011</year><pub-id pub-id-type="pmid">22008120</pub-id></element-citation></ref><ref id="R27"><label>27</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Deng</surname><given-names>Y</given-names></name><name><surname>Iannello</surname><given-names>P</given-names></name><name><surname>Smith</surname><given-names>I</given-names></name></person-group><article-title>Transmission of influenza A (H1N1) 2009 pandemic viruses in Australian swine</article-title><source>Influ Other Respi Viruses</source><volume>6</volume><fpage>42</fpage><lpage>47</lpage><year>2012</year></element-citation></ref><ref id="R28"><label>28</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Biondo</surname><given-names>N</given-names></name><etal></etal></person-group><article-title>Genomic analysis of influenza A virus from captive wild boars in Brazil reveals a human-like H1N2 influenza virus</article-title><source>Vet Microbiol</source><volume>168</volume><fpage>34</fpage><lpage>40</lpage><year>2014</year><pub-id pub-id-type="pmid">24238665</pub-id></element-citation></ref><ref id="R29"><label>29</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Schaefer</surname><given-names>R</given-names></name><name><surname>Zanella</surname><given-names>JRC</given-names></name><name><surname>Brentano</surname><given-names>L</given-names></name><name><surname>Vincent</surname><given-names>AL</given-names></name><name><surname>Ritterbusch</surname><given-names>GA</given-names></name><name><surname>Silveira</surname><given-names>S</given-names></name><name><surname>Caron</surname><given-names>L</given-names></name><name><surname>M</surname><given-names>N</given-names></name></person-group><article-title>Isolation characterization of a pandemic H1N1 influenza virus in pigs in Brazil</article-title><source>Pesqui Vet Bras</source><volume>31</volume><fpage>761</fpage><lpage>767</lpage><year>2011</year></element-citation></ref><ref id="R30"><label>30</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nagarajan</surname><given-names>K</given-names></name><etal></etal></person-group><article-title>Influenza A H1N1 virus in Indian pigs & its genetic relatedness with pandemic human influenza A 2009 H1N1</article-title><source>Indian J Med Res</source><volume>132</volume><fpage>160</fpage><lpage>167</lpage><year>2010</year><pub-id pub-id-type="pmid">20716816</pub-id></element-citation></ref><ref id="R31"><label>31</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Njabo</surname><given-names>KY</given-names></name><etal></etal></person-group><article-title>Pandemic A/H1N1/2009 influenza virus in swine, Cameroon, 2010</article-title><source>Vet Microbiol</source><volume>156</volume><fpage>189</fpage><lpage>192</lpage><year>2012</year><pub-id pub-id-type="pmid">21963416</pub-id></element-citation></ref><ref id="R32"><label>32</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Escalera-Zamudio</surname><given-names>M</given-names></name><etal></etal></person-group><article-title>Characterization of an influenza A virus in Mexican swine that is related to the A/H1N1/2009 pandemic clade</article-title><source>Virology</source><volume>433</volume><fpage>176</fpage><lpage>182</lpage><year>2012</year><pub-id pub-id-type="pmid">22921315</pub-id></element-citation></ref><ref id="R33"><label>33</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Meseko</surname><given-names>CA</given-names></name><name><surname>Odaibo</surname><given-names>GN</given-names></name><name><surname>Olaleye</surname><given-names>DO</given-names></name></person-group><article-title>Detection and isolation of 2009 pandemic influenza A/H1N1 virus in commercial piggery, Lagos Nigeria</article-title><source>Vet Microbiol</source><volume>168</volume><fpage>197</fpage><lpage>201</lpage><year>2014</year><pub-id pub-id-type="pmid">24321146</pub-id></element-citation></ref><ref id="R34"><label>34</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Perera</surname><given-names>HKK</given-names></name><etal></etal></person-group><article-title>Swine influenza in Sri Lanka</article-title><source>Emerg Infect Dis</source><volume>19</volume><fpage>481</fpage><lpage>484</lpage><year>2013</year><pub-id pub-id-type="pmid">23621918</pub-id></element-citation></ref><ref id="R35"><label>35</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hofshagen</surname><given-names>M</given-names></name><name><surname>Gjerset</surname><given-names>B</given-names></name><name><surname>Er</surname><given-names>C</given-names></name></person-group><article-title>Pandemic Influenza A(H1N1)v: Human to pig transmission in Norway</article-title><source>Euro Surveill</source><volume>14</volume><fpage>1</fpage><lpage>3</lpage><year>2009</year></element-citation></ref><ref id="R36"><label>36</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nokireki</surname><given-names>T</given-names></name><name><surname>Laine</surname><given-names>T</given-names></name><name><surname>London</surname><given-names>L</given-names></name><name><surname>Ikonen</surname><given-names>N</given-names></name><name><surname>Huovilainen</surname><given-names>A</given-names></name></person-group><article-title>The first detection of influenza in the Finnish pig population: a retrospective study</article-title><source>Acta Vet Scand</source><volume>55</volume><fpage>69</fpage><year>2013</year><pub-id pub-id-type="pmid">24047612</pub-id></element-citation></ref><ref id="R37"><label>37</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bányai</surname><given-names>K</given-names></name><etal></etal></person-group><article-title>Genome sequence of a monoreassortant H1N1 swine influenza virus isolated from a pig in Hungary</article-title><source>J Virol</source><volume>86</volume><fpage>13133</fpage><year>2012</year><pub-id pub-id-type="pmid">23118459</pub-id></element-citation></ref><ref id="R38"><label>38</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lam</surname><given-names>TTY</given-names></name><etal></etal></person-group><article-title>Reassortment events among swine influenza A viruses in China: implications for the origin of the 2009 influenza pandemic</article-title><source>J Virol</source><volume>85</volume><fpage>10279</fpage><lpage>10285</lpage><year>2011</year><pub-id pub-id-type="pmid">21795347</pub-id></element-citation></ref><ref id="R39"><label>39</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cappuccio</surname><given-names>Ja</given-names></name><etal></etal></person-group><article-title>Outbreak of swine influenza in Argentina reveals a non-contemporary human H3N2 virus highly transmissible among pigs</article-title><source>J Gen Virol</source><volume>92</volume><fpage>2871</fpage><lpage>2878</lpage><year>2011</year><pub-id pub-id-type="pmid">21849519</pub-id></element-citation></ref><ref id="R40"><label>40</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Schaefer</surname><given-names>R</given-names></name><etal></etal></person-group><article-title>A human-like H1N2 influenza virus detected during an outbreak of acute respiratory disease in swine in Brazil</article-title><source>Arch Virol</source><year>2014</year><pub-id pub-id-type="doi">10.1007/s00705-014-2223-z</pub-id></element-citation></ref><ref id="R41"><label>41</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhou</surname><given-names>B</given-names></name><etal></etal></person-group><article-title>Single-reaction genomic amplification accelerates sequencing and vaccine production for classical and Swine origin human influenza a viruses</article-title><source>J Virol</source><volume>83</volume><fpage>10309</fpage><lpage>13</lpage><year>2009</year><pub-id pub-id-type="pmid">19605485</pub-id></element-citation></ref><ref id="R42"><label>42</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bao</surname><given-names>Y</given-names></name><etal></etal></person-group><article-title>The influenza virus resource at the National Center for Biotechnology Information</article-title><source>J Virol</source><volume>82</volume><fpage>596</fpage><lpage>601</lpage><year>2008</year><pub-id pub-id-type="pmid">17942553</pub-id></element-citation></ref><ref id="R43"><label>43</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Edgar</surname><given-names>RC</given-names></name></person-group><article-title>MUSCLE: multiple sequence alignment with high accuracy and high throughput</article-title><source>Nucleic Acids Res</source><volume>32</volume><fpage>1792</fpage><lpage>7</lpage><year>2004</year><pub-id pub-id-type="pmid">15034147</pub-id></element-citation></ref><ref id="R44"><label>44</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nelson</surname><given-names>MI</given-names></name><name><surname>Gramer</surname><given-names>MR</given-names></name><name><surname>Vincent</surname><given-names>AL</given-names></name><name><surname>Holmes</surname><given-names>EC</given-names></name></person-group><article-title>Global transmission of influenza viruses from humans to swine</article-title><source>J Gen Virol</source><volume>93</volume><fpage>2195</fpage><lpage>2203</lpage><year>2012</year><pub-id pub-id-type="pmid">22791604</pub-id></element-citation></ref><ref id="R45"><label>45</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Drummond</surname><given-names>AJ</given-names></name><name><surname>Suchard</surname><given-names>MA</given-names></name><name><surname>Xie</surname><given-names>D</given-names></name><name><surname>Rambaut</surname><given-names>A</given-names></name></person-group><article-title>Bayesian phylogenetics with BEAUti and the BEAST 1.7</article-title><source>Mol Biol Evol</source><volume>29</volume><fpage>1969</fpage><lpage>73</lpage><year>2012</year><pub-id pub-id-type="pmid">22367748</pub-id></element-citation></ref><ref id="R46"><label>46</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Suchard</surname><given-names>MA</given-names></name><name><surname>Rambaut</surname><given-names>A</given-names></name></person-group><article-title>Many-core algorithms for statistical phylogenetics</article-title><source>Bioinformatics</source><volume>25</volume><fpage>1370</fpage><lpage>6</lpage><year>2009</year><pub-id pub-id-type="pmid">19369496</pub-id></element-citation></ref><ref id="R47"><label>47</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Viboud</surname><given-names>C</given-names></name><etal></etal></person-group><article-title>Synchrony, waves, and spatial hierarchies in the spread of influenza</article-title><source>Science</source><volume>312</volume><fpage>447</fpage><lpage>51</lpage><year>2006</year><pub-id pub-id-type="pmid">16574822</pub-id></element-citation></ref><ref id="R48"><label>48</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Romagosa</surname><given-names>A</given-names></name><etal></etal></person-group><article-title>Vaccination of influenza a virus decreases transmission rates in pigs</article-title><source>Vet Res</source><volume>42</volume><fpage>120</fpage><year>2011</year><pub-id pub-id-type="pmid">22185601</pub-id></element-citation></ref><ref id="R49"><label>49</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wong</surname><given-names>KK</given-names></name><etal></etal></person-group><article-title>Transmissibility of variant influenza from Swine to humans: a modeling approach</article-title><source>Clin Infect Dis</source><volume>57</volume><issue>Suppl 1</issue><fpage>S16</fpage><lpage>22</lpage><year>2013</year><pub-id pub-id-type="pmid">23794727</pub-id></element-citation></ref><ref id="R50"><label>50</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Robinson</surname><given-names>TP</given-names></name><etal></etal></person-group><article-title>Mapping the global distribution of livestock</article-title><source>PLoS One</source><volume>9</volume><fpage>e96084</fpage><year>2014</year><pub-id pub-id-type="pmid">24875496</pub-id></element-citation></ref></ref-list></back><floats-group><fig id="F1" orientation="portrait" position="float"><label>Fig. 1</label><caption><title>Modeled global swine distributions</title><p>Digital layers from Gridded Livestock of the World (GLW) (version 2.01)<sup><xref rid="R50" ref-type="bibr">50</xref></sup>, downloaded from the publically available Livestock Geo-Wiki database (<ext-link ext-link-type="uri" xlink:href="http://www.livestock.geo-wiki.org">http://www.livestock.geo-wiki.org</ext-link>) and manually edited in QGIS v.1.7.0. Swine densities are represented by black shading.</p></caption><graphic xlink:href="nihms666251f1"></graphic></fig><fig id="F2" orientation="portrait" position="float"><label>Fig. 2</label><caption><title>Inter-continental migration events of swIAVs</title><p>Circles represent the country of origin, based on the estimates summarized in the MCC tree, and are shaded accordingly. Lines represent the inferred time period of inter-continental transmission, within a level of uncertainty, inferred from the estimated date of ancestral nodes on the MCC tree. Triangles represent clades resulting from onwards transmission of the introduced viruses, are shaded by the country of destination, and extend as far forward in time as the most recently sampled virus. Numbers of introduction (1–18) correspond to the clade numbers on the phylogenies (<xref rid="F3" ref-type="fig">Fig. 3</xref> and <xref ref-type="supplementary-material" rid="SD1">Supplementary Figs. 1–7</xref>). The asterisks indicate that additional HA and NA swIAV sequence data was used to estimate the timing of introduction 18. Countries/regions are abbreviated as follows: CHN = China (including Hong Kong SAR and Taiwan), THA = Thailand, VNM = Vietnam, KOR = South Korea, JPN = Japan, MEX = Mexico, and EUR = Europe.</p></caption><graphic xlink:href="nihms666251f2"></graphic></fig><fig id="F3" orientation="portrait" position="float"><label>Fig. 3</label><caption><title>MCC trees of the NA lineages in swine</title><p>Time-scaled Bayesian MCC trees inferred for the NA segment for the three major swine virus lineages: (a) avian-origin Eurasian N1 swIAV lineage, the (b) classical N1 swIAV lineage, and the (c) multiple human seasonal virus-origin N2 swIAV lineages circulating in swine. Branches of human seasonal H3N2 influenza virus origin are shaded grey in (c), while branches associated with viruses from swine are shaded by country of origin: Argentina = brown; Canada = red; China (including Hong Kong SAR and Taiwan) = yellow; Europe = black; Japan = pink; Mexico = light blue; South Korea = green; Thailand = orange; USA = dark blue; Vietnam = purple. Posterior probabilities > 0.8 are included for key nodes, and international migration events that are supported by high posterior probabilities and long branch lengths are labeled according to <xref rid="F2" ref-type="fig">Fig. 2</xref>.</p></caption><graphic xlink:href="nihms666251f3"></graphic></fig><fig id="F4" orientation="portrait" position="float"><label>Fig. 4</label><caption><title>Heat-map of swIAV migration between locations</title><p>Countries are listed in order of increasing geographical distance from Argentina (ARG). MEX = Mexico, USA = United States, CAN = Canada, EUR = Europe, JPN = Japan, CHN = China (including Hong Kong SAR and Taiwan), KOR = South Korea, THA = Thailand, VNM = Vietnam. The intensity of the color (red = high; white = low) reflects the number of ‘Markov jump’ counts inferred over the totality of phylogenies (all segments, all lineages) from one location to another (asymmetrical). Markov jump counts measure the number of inferred location state transitions, modelled by a continuous-time Markov chain process, that occur along the branches of the phylogeny. For clarity the heat-map has been divided into four sections representing (a) viral migration events within the Americas and between the Americas and Europe; (b) migrations from the Americas/Europe to Asia; (c) migrations from Asia to the Americas/Europe; and (d) migrations between Asian countries.</p></caption><graphic xlink:href="nihms666251f4"></graphic></fig><fig id="F5" orientation="portrait" position="float"><label>Fig. 5</label><caption><title>The support and contribution of swIAV diffusion predictors among 9 countries</title><p>Twelve predictors were considered: geographical distance (km), volume of live swine trade, 1996–2012 (USD), swine population size for the years 1969–2010, the total number of imports of live swine during 1969–2010, the total number of swine exports during 1969–2010, the percent change in swine population size from 1969–2010, and the number of sequences available from a given country for our analysis. ‘O’ refers to the swine population of origin, and ‘d’ refers to the swine population of destination. Support for each predictor is represented by an inclusion probability that is estimated as the posterior expectation for the indicator variable associated with each predictor (<bold>E</bold>[<italic>δ</italic>]). The contribution of each predictor is represented by the mean and credible intervals of the GLM coefficients (<italic>β</italic>) on a log scale conditional on the predictor being included in the model (<italic>β</italic>|<italic>δ</italic>=1). See <xref ref-type="supplementary-material" rid="SD1">Supplementary Fig. 8</xref> for MP and NS results.</p></caption><graphic xlink:href="nihms666251f5"></graphic></fig><fig id="F6" orientation="portrait" position="float"><label>Fig. 6</label><caption><title>Maps of simulated spread of influenza viruses via live swine trade flows</title><p>Simulated spread of an influenza virus from 5 seed countries (shaded in black) to 146 countries for which live swine trade is available from the United Nations Commodity Trade Statistics Database (available at <ext-link ext-link-type="uri" xlink:href="http://comtrade.un.org">http://comtrade.un.org</ext-link>) (a–e). Probability of an outbreak in the invaded country is shaded from white (probability of 0) to red (probability of 1). The probability of co-invasion by both a virus seeded in North America (Canada and the United States) and Europe also is shaded from white (probability of 0) to red (probability of 1) (f). Arrows represent the direction of viral dissemination for countries with a probability of an outbreak > 0.25 (see <xref ref-type="supplementary-material" rid="SD1">Supplementary Table 5</xref> for a complete list of all outbreak probabilities by country).</p></caption><graphic xlink:href="nihms666251f6"></graphic></fig></floats-group></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 000E59 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000E59 | 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. |  |