Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Epidemiological and genetic characterization of pH1N1 and H3N2 influenza viruses circulated in MENA region during 2009–2017</title><author><name sortKey="Al Khatib, Hebah A" sort="Al Khatib, Hebah A" uniqKey="Al Khatib H" first="Hebah A." last="Al Khatib">Hebah A. Al Khatib</name><affiliation><nlm:aff id="Aff1">Life Science division, College of Science and Engineering, Hamad Ben Khalifah University, Doha, 34110 Qatar</nlm:aff></affiliation></author><author><name sortKey="Al Thani, Asmaa A" sort="Al Thani, Asmaa A" uniqKey="Al Thani A" first="Asmaa A." last="Al Thani">Asmaa A. Al Thani</name><affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 0634 1084</institution-id><institution-id institution-id-type="GRID">grid.412603.2</institution-id><institution>Biomedical Research Center,</institution><institution>Qatar University,</institution></institution-wrap>Doha, 2713 Qatar</nlm:aff></affiliation></author><author><name sortKey="Gallouzi, Imed" sort="Gallouzi, Imed" uniqKey="Gallouzi I" first="Imed" last="Gallouzi">Imed Gallouzi</name><affiliation><nlm:aff id="Aff1">Life Science division, College of Science and Engineering, Hamad Ben Khalifah University, Doha, 34110 Qatar</nlm:aff></affiliation><affiliation><nlm:aff id="Aff3"><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 1936 8649</institution-id><institution-id institution-id-type="GRID">grid.14709.3b</institution-id><institution>Biochemistry Department and Goodman Cancer Center, 3655 Promenade Sir William Osler,</institution><institution>McGill University,</institution></institution-wrap>Montreal, Quebec, H3G1Y6 Canada</nlm:aff></affiliation></author><author><name sortKey="Yassine, Hadi M" sort="Yassine, Hadi M" uniqKey="Yassine H" first="Hadi M." last="Yassine">Hadi M. Yassine</name><affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 0634 1084</institution-id><institution-id institution-id-type="GRID">grid.412603.2</institution-id><institution>Biomedical Research Center,</institution><institution>Qatar University,</institution></institution-wrap>Doha, 2713 Qatar</nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">30971204</idno><idno type="pmc">6458790</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC6458790</idno><idno type="RBID">PMC:6458790</idno><idno type="doi">10.1186/s12879-019-3930-6</idno><date when="2019">2019</date><idno type="wicri:Area/Pmc/Corpus">000165</idno><idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000165</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Epidemiological and genetic characterization of pH1N1 and H3N2 influenza viruses circulated in MENA region during 2009–2017</title><author><name sortKey="Al Khatib, Hebah A" sort="Al Khatib, Hebah A" uniqKey="Al Khatib H" first="Hebah A." last="Al Khatib">Hebah A. Al Khatib</name><affiliation><nlm:aff id="Aff1">Life Science division, College of Science and Engineering, Hamad Ben Khalifah University, Doha, 34110 Qatar</nlm:aff></affiliation></author><author><name sortKey="Al Thani, Asmaa A" sort="Al Thani, Asmaa A" uniqKey="Al Thani A" first="Asmaa A." last="Al Thani">Asmaa A. Al Thani</name><affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 0634 1084</institution-id><institution-id institution-id-type="GRID">grid.412603.2</institution-id><institution>Biomedical Research Center,</institution><institution>Qatar University,</institution></institution-wrap>Doha, 2713 Qatar</nlm:aff></affiliation></author><author><name sortKey="Gallouzi, Imed" sort="Gallouzi, Imed" uniqKey="Gallouzi I" first="Imed" last="Gallouzi">Imed Gallouzi</name><affiliation><nlm:aff id="Aff1">Life Science division, College of Science and Engineering, Hamad Ben Khalifah University, Doha, 34110 Qatar</nlm:aff></affiliation><affiliation><nlm:aff id="Aff3"><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 1936 8649</institution-id><institution-id institution-id-type="GRID">grid.14709.3b</institution-id><institution>Biochemistry Department and Goodman Cancer Center, 3655 Promenade Sir William Osler,</institution><institution>McGill University,</institution></institution-wrap>Montreal, Quebec, H3G1Y6 Canada</nlm:aff></affiliation></author><author><name sortKey="Yassine, Hadi M" sort="Yassine, Hadi M" uniqKey="Yassine H" first="Hadi M." last="Yassine">Hadi M. Yassine</name><affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 0634 1084</institution-id><institution-id institution-id-type="GRID">grid.412603.2</institution-id><institution>Biomedical Research Center,</institution><institution>Qatar University,</institution></institution-wrap>Doha, 2713 Qatar</nlm:aff></affiliation></author></analytic><series><title level="j">BMC Infectious Diseases</title><idno type="eISSN">1471-2334</idno><imprint><date when="2019">2019</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><sec><title>Background</title><p id="Par1">Influenza surveillance is necessary for detection of emerging variants of epidemiologic and clinical significance. This study describes the epidemiology of influenza types A and B, and molecular characteristics of surface glycoproteins (hemagglutinin [HA] and neuraminidase [NA]) of influenza A subtypes: pH1N1 and H3N2 circulated in Arabian Gulf, Levant and North Africa regions during 2009–2017.</p></sec><sec><title>Methods</title><p id="Par2">Analysis of phylogenetics and evolution of HA and NA genes was done using full HA and NA sequences (<italic>n</italic> = 1229) downloaded from Influenza Research Database (IRD).</p></sec><sec><title>Results</title><p id="Par3">In total, 130,354 influenza positive cases were reported to WHO during study period. Of these, 50.8% were pH1N1 positive, 15.9% were H3N2 positives and 17.2% were influenza B positive. With few exceptions, all three regions were showing the typical seasonal influenza peak similar to that reported in Northern hemisphere (December–March). However, influenza activity started earlier (October) in both Gulf and North Africa while commenced later during November in Levant countries. The molecular analysis of the HA genes (influenza A subtypes) revealed similar mutations to those reported worldwide. Generally, amino acid substitutions were most frequently found in head domain in H1N1 pandemic viruses, while localized mainly in the stem region in H3N2 viruses. Expectedly, seasons with high pH1N1 influenza activity was associated with a relatively higher number of substitutions in the head domain of the HA in pH1N1 subtype. Furthermore, nucleotide variations were lower at the antigenic sites of pH1N1 viruses compared to H3N2 viruses, which experienced higher variability at the antigenic sites, reflecting the increased immunological pressure because of longer circulation and continuous vaccine changes. Analysis of NA gene of pH1N1 viruses revealed sporadic detections of oseltamivir-resistance mutation, H275Y, in 4% of reported sequences, however, none of NAI resistance mutations were found in the NA of H3N2 viruses.</p></sec><sec><title>Conclusions</title><p id="Par4">Molecular characterization of H1N1 and H3N2 viruses over 9 years revealed significant differences with regard to position and function of characterized substitutions. While pH1N1 virus substitutions were mainly found in HA head domain, H3N2 virus substitutions were mostly found in HA stem domain. Additionally, more fixed substitutions were encountered in H3N2 virus compared to larger number of non-fixed substitutions in pH1N1.</p></sec><sec><title>Electronic supplementary material</title><p>The online version of this article (10.1186/s12879-019-3930-6) contains supplementary material, which is available to authorized users.</p></sec></div></front><back><div1 type="bibliography"><listBibl><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Hay, Aj" uniqKey="Hay A">AJ Hay</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Molinari, Na" uniqKey="Molinari N">NA Molinari</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Schafer, Jr" uniqKey="Schafer J">JR Schafer</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Webster, Rg" uniqKey="Webster R">RG Webster</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Cohen, J" uniqKey="Cohen J">J. Cohen</name></author><author><name sortKey="Enserink, M" uniqKey="Enserink M">M. Enserink</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Cohen, Al" uniqKey="Cohen A">AL Cohen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Martinez, O" uniqKey="Martinez O">O Martinez</name></author><author><name sortKey="Tsibane, T" uniqKey="Tsibane T">T Tsibane</name></author><author><name sortKey="Basler, Cf" uniqKey="Basler C">CF Basler</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gerhard, W" uniqKey="Gerhard W">W Gerhard</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wiley, Dc" uniqKey="Wiley D">DC Wiley</name></author><author><name sortKey="Skehel, Jj" uniqKey="Skehel J">JJ Skehel</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Igarashi, M" uniqKey="Igarashi M">M Igarashi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gallagher, P" uniqKey="Gallagher P">P Gallagher</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sun, S" uniqKey="Sun S">S Sun</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Chen, Z" uniqKey="Chen Z">Z Chen</name></author><author><name sortKey="Zhou, H" uniqKey="Zhou H">H Zhou</name></author><author><name sortKey="Jin, H" uniqKey="Jin H">H Jin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Koel, Bf" uniqKey="Koel B">BF Koel</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Haq, Z" uniqKey="Haq Z">Z Haq</name></author><author><name sortKey="Mahjour, J" uniqKey="Mahjour J">J Mahjour</name></author><author><name sortKey="Khan, W" uniqKey="Khan W">W Khan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Baghi, Hossein Bannazadeh" uniqKey="Baghi H">Hossein Bannazadeh Baghi</name></author><author><name sortKey="Soroush, Mohammad Hossein" uniqKey="Soroush M">Mohammad Hossein Soroush</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Buliva, E" uniqKey="Buliva E">E Buliva</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Soebiyanto, Rp" uniqKey="Soebiyanto R">RP Soebiyanto</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Khan, W" uniqKey="Khan W">W Khan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tamerius, Jd" uniqKey="Tamerius J">JD Tamerius</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Nguyen, Ht" uniqKey="Nguyen H">HT Nguyen</name></author><author><name sortKey="Fry, Am" uniqKey="Fry A">AM Fry</name></author><author><name sortKey="Gubareva, Lv" uniqKey="Gubareva L">LV Gubareva</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tamura, K" uniqKey="Tamura K">K Tamura</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yang, H" uniqKey="Yang H">H Yang</name></author><author><name sortKey="Carney, P" uniqKey="Carney P">P Carney</name></author><author><name sortKey="Stevens, J" uniqKey="Stevens J">J Stevens</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Huang, J W" uniqKey="Huang J">J-W Huang</name></author><author><name sortKey="Yang, J M" uniqKey="Yang J">J-M Yang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Choi, Y" uniqKey="Choi Y">Y Choi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Drummond, Aj" uniqKey="Drummond A">AJ Drummond</name></author><author><name sortKey="Rambaut, A" uniqKey="Rambaut A">A Rambaut</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kosakovsky Pond, Sl" uniqKey="Kosakovsky Pond S">SL Kosakovsky Pond</name></author><author><name sortKey="Frost, Sd" uniqKey="Frost S">SD Frost</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Viboud, C" uniqKey="Viboud C">C Viboud</name></author><author><name sortKey="Alonso, Wj" uniqKey="Alonso W">WJ Alonso</name></author><author><name sortKey="Simonsen, L" uniqKey="Simonsen L">L Simonsen</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Chutinimitkul, S" uniqKey="Chutinimitkul S">S Chutinimitkul</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Schulze, It" uniqKey="Schulze I">IT Schulze</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sriwilaijaroen, N" uniqKey="Sriwilaijaroen N">N Sriwilaijaroen</name></author><author><name sortKey="Suzuki, Y" uniqKey="Suzuki Y">Y Suzuki</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Air, Gm" uniqKey="Air G">GM Air</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Maurer Stroh, S" uniqKey="Maurer Stroh S">S Maurer-Stroh</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Graham, M" uniqKey="Graham M">M Graham</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Liu, Ss" uniqKey="Liu S">SS Liu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Butler, J" uniqKey="Butler J">J Butler</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bragstad, K" uniqKey="Bragstad K">K Bragstad</name></author><author><name sortKey="Nielsen, Lp" uniqKey="Nielsen L">LP Nielsen</name></author><author><name sortKey="Fomsgaard, A" uniqKey="Fomsgaard A">A Fomsgaard</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Caini, S" uniqKey="Caini S">S Caini</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Simonsen, L" uniqKey="Simonsen L">L Simonsen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Simonsen, L" uniqKey="Simonsen L">L Simonsen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bedford, T" uniqKey="Bedford T">T Bedford</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Smith, Dj" uniqKey="Smith D">DJ Smith</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Klein, Ey" uniqKey="Klein E">EY Klein</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Liu, Y" uniqKey="Liu Y">Y Liu</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Gamblin, Sj" uniqKey="Gamblin S">SJ Gamblin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Maines, Tr" uniqKey="Maines T">TR Maines</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="De Vries, Rp" uniqKey="De Vries R">RP de Vries</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jayaraman, A" uniqKey="Jayaraman A">A Jayaraman</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Luksza, M" uniqKey="Luksza M">M Luksza</name></author><author><name sortKey="Lassig, M" uniqKey="Lassig M">M Lassig</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wei, Cj" uniqKey="Wei C">CJ Wei</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ison, Mg" uniqKey="Ison M">MG Ison</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Meijer, A" uniqKey="Meijer A">A Meijer</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Storms, Ad" uniqKey="Storms A">AD Storms</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gubareva, Larisa V" uniqKey="Gubareva L">Larisa V. Gubareva</name></author><author><name sortKey="Besselaar, Terry G" uniqKey="Besselaar T">Terry G. Besselaar</name></author><author><name sortKey="Daniels, Rod S" uniqKey="Daniels R">Rod S. Daniels</name></author><author><name sortKey="Fry, Alicia" uniqKey="Fry A">Alicia Fry</name></author><author><name sortKey="Gregory, Vicki" uniqKey="Gregory V">Vicki Gregory</name></author><author><name sortKey="Huang, Weijuan" uniqKey="Huang W">Weijuan Huang</name></author><author><name sortKey="Hurt, Aeron C" uniqKey="Hurt A">Aeron C. Hurt</name></author><author><name sortKey="Jorquera, Patricia A" uniqKey="Jorquera P">Patricia A. Jorquera</name></author><author><name sortKey="Lackenby, Angie" uniqKey="Lackenby A">Angie Lackenby</name></author><author><name sortKey="Leang, Sook Kwan" uniqKey="Leang S">Sook-Kwan Leang</name></author><author><name sortKey="Lo, Janice" uniqKey="Lo J">Janice Lo</name></author><author><name sortKey="Pereyaslov, Dmitriy" uniqKey="Pereyaslov D">Dmitriy Pereyaslov</name></author><author><name sortKey="Rebelo De Andrade, Helena" uniqKey="Rebelo De Andrade H">Helena Rebelo-de-Andrade</name></author><author><name sortKey="Siqueira, Marilda M" uniqKey="Siqueira M">Marilda M. Siqueira</name></author><author><name sortKey="Takashita, Emi" uniqKey="Takashita E">Emi Takashita</name></author><author><name sortKey="Odagiri, Takato" uniqKey="Odagiri T">Takato Odagiri</name></author><author><name sortKey="Wang, Dayan" uniqKey="Wang D">Dayan Wang</name></author><author><name sortKey="Zhang, Wenqing" uniqKey="Zhang W">Wenqing Zhang</name></author><author><name sortKey="Meijer, Adam" uniqKey="Meijer A">Adam Meijer</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hurt, Ac" uniqKey="Hurt A">AC Hurt</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Takashita, E" uniqKey="Takashita E">E Takashita</name></author><author><name sortKey="Ejima, M" uniqKey="Ejima M">M Ejima</name></author><author><name sortKey="Itoh, R" uniqKey="Itoh R">R Itoh</name></author><author><name sortKey="Miura, M" uniqKey="Miura M">M Miura</name></author><author><name sortKey="Ohnishi, A" uniqKey="Ohnishi A">A Ohnishi</name></author><author><name sortKey="Nishimura, H" uniqKey="Nishimura H">H Nishimura</name></author><author><name sortKey="Odagiri, T" uniqKey="Odagiri T">T Odagiri</name></author><author><name sortKey="Tashiro, M" uniqKey="Tashiro M">M Tashiro</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chen, Lf" uniqKey="Chen L">LF Chen</name></author></analytic></biblStruct><biblStruct></biblStruct></listBibl></div1></back></TEI><pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">BMC Infect Dis</journal-id><journal-id journal-id-type="iso-abbrev">BMC Infect. Dis</journal-id><journal-title-group><journal-title>BMC Infectious Diseases</journal-title></journal-title-group><issn pub-type="epub">1471-2334</issn><publisher><publisher-name>BioMed Central</publisher-name><publisher-loc>London</publisher-loc></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">30971204</article-id><article-id pub-id-type="pmc">6458790</article-id><article-id pub-id-type="publisher-id">3930</article-id><article-id pub-id-type="doi">10.1186/s12879-019-3930-6</article-id><article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject></subj-group></article-categories><title-group><article-title>Epidemiological and genetic characterization of pH1N1 and H3N2 influenza viruses circulated in MENA region during 2009–2017</article-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Al Khatib</surname><given-names>Hebah A.</given-names></name><address><email>halkhatib@qf.org.qa</email></address><xref ref-type="aff" rid="Aff1">1</xref></contrib><contrib contrib-type="author"><name><surname>Al Thani</surname><given-names>Asmaa A.</given-names></name><address><email>aaja@qu.edu.qa</email></address><xref ref-type="aff" rid="Aff2">2</xref></contrib><contrib contrib-type="author"><name><surname>Gallouzi</surname><given-names>Imed</given-names></name><address><email>igallouzi@qf.org.qa</email></address><xref ref-type="aff" rid="Aff1">1</xref><xref ref-type="aff" rid="Aff3">3</xref></contrib><contrib contrib-type="author" corresp="yes"><contrib-id contrib-id-type="orcid">http://orcid.org/0000-0001-7592-2788</contrib-id><name><surname>Yassine</surname><given-names>Hadi M.</given-names></name><address><phone>00 974 44036819</phone><email>Hyassine@qu.edu.qa</email></address><xref ref-type="aff" rid="Aff2">2</xref></contrib><aff id="Aff1"><label>1</label>Life Science division, College of Science and Engineering, Hamad Ben Khalifah University, Doha, 34110 Qatar</aff><aff id="Aff2"><label>2</label><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 0634 1084</institution-id><institution-id institution-id-type="GRID">grid.412603.2</institution-id><institution>Biomedical Research Center,</institution><institution>Qatar University,</institution></institution-wrap>Doha, 2713 Qatar</aff><aff id="Aff3"><label>3</label><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 1936 8649</institution-id><institution-id institution-id-type="GRID">grid.14709.3b</institution-id><institution>Biochemistry Department and Goodman Cancer Center, 3655 Promenade Sir William Osler,</institution><institution>McGill University,</institution></institution-wrap>Montreal, Quebec, H3G1Y6 Canada</aff></contrib-group><pub-date pub-type="epub"><day>11</day><month>4</month><year>2019</year></pub-date><pub-date pub-type="pmc-release"><day>11</day><month>4</month><year>2019</year></pub-date><pub-date pub-type="collection"><year>2019</year></pub-date><volume>19</volume><elocation-id>314</elocation-id><history><date date-type="received"><day>13</day><month>6</month><year>2018</year></date><date date-type="accepted"><day>20</day><month>3</month><year>2019</year></date></history><permissions><copyright-statement>© The Author(s). 2019</copyright-statement><license license-type="OpenAccess"><license-p><bold>Open Access</bold>This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">http://creativecommons.org/licenses/by/4.0/</ext-link>), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/publicdomain/zero/1.0/">http://creativecommons.org/publicdomain/zero/1.0/</ext-link>) applies to the data made available in this article, unless otherwise stated.</license-p></license></permissions><abstract id="Abs1"><sec><title>Background</title><p id="Par1">Influenza surveillance is necessary for detection of emerging variants of epidemiologic and clinical significance. This study describes the epidemiology of influenza types A and B, and molecular characteristics of surface glycoproteins (hemagglutinin [HA] and neuraminidase [NA]) of influenza A subtypes: pH1N1 and H3N2 circulated in Arabian Gulf, Levant and North Africa regions during 2009–2017.</p></sec><sec><title>Methods</title><p id="Par2">Analysis of phylogenetics and evolution of HA and NA genes was done using full HA and NA sequences (<italic>n</italic> = 1229) downloaded from Influenza Research Database (IRD).</p></sec><sec><title>Results</title><p id="Par3">In total, 130,354 influenza positive cases were reported to WHO during study period. Of these, 50.8% were pH1N1 positive, 15.9% were H3N2 positives and 17.2% were influenza B positive. With few exceptions, all three regions were showing the typical seasonal influenza peak similar to that reported in Northern hemisphere (December–March). However, influenza activity started earlier (October) in both Gulf and North Africa while commenced later during November in Levant countries. The molecular analysis of the HA genes (influenza A subtypes) revealed similar mutations to those reported worldwide. Generally, amino acid substitutions were most frequently found in head domain in H1N1 pandemic viruses, while localized mainly in the stem region in H3N2 viruses. Expectedly, seasons with high pH1N1 influenza activity was associated with a relatively higher number of substitutions in the head domain of the HA in pH1N1 subtype. Furthermore, nucleotide variations were lower at the antigenic sites of pH1N1 viruses compared to H3N2 viruses, which experienced higher variability at the antigenic sites, reflecting the increased immunological pressure because of longer circulation and continuous vaccine changes. Analysis of NA gene of pH1N1 viruses revealed sporadic detections of oseltamivir-resistance mutation, H275Y, in 4% of reported sequences, however, none of NAI resistance mutations were found in the NA of H3N2 viruses.</p></sec><sec><title>Conclusions</title><p id="Par4">Molecular characterization of H1N1 and H3N2 viruses over 9 years revealed significant differences with regard to position and function of characterized substitutions. While pH1N1 virus substitutions were mainly found in HA head domain, H3N2 virus substitutions were mostly found in HA stem domain. Additionally, more fixed substitutions were encountered in H3N2 virus compared to larger number of non-fixed substitutions in pH1N1.</p></sec><sec><title>Electronic supplementary material</title><p>The online version of this article (10.1186/s12879-019-3930-6) contains supplementary material, which is available to authorized users.</p></sec></abstract><kwd-group xml:lang="en"><title>Keywords</title><kwd>Influenza epidemics</kwd><kwd>pH1N1</kwd><kwd>H3N2</kwd><kwd>Phylogeny</kwd><kwd>Hemagglutinin</kwd><kwd>Neuraminidase</kwd><kwd>Molecular evolution</kwd></kwd-group><custom-meta-group><custom-meta><meta-name>issue-copyright-statement</meta-name><meta-value>© The Author(s) 2019</meta-value></custom-meta></custom-meta-group></article-meta></front><body><sec id="Sec1"><title>Background</title><p id="Par11">Recently published data indicates that up to 650,000 deaths are reported annually due to respiratory illness associated with seasonal influenza [<xref ref-type="bibr" rid="CR1">1</xref>]. Influenza viruses evolve rapidly, escaping natural or vaccine-induced immune response by accumulating mutations especially within surface glycoproteins genes: hemagglutinin (HA) and neuraminidase (NA) [<xref ref-type="bibr" rid="CR2">2</xref>]. This antigenic drift results in annual influenza epidemics in humans, a major cause of morbidity and mortality and consequently increased burden on health services [<xref ref-type="bibr" rid="CR3">3</xref>]. In addition to antigenic drift, influenza viruses evolve by genetic reassortment, leading in some occasions to life threatening pandemics. Three influenza A pandemics were reported in last century: 1918 H1N1, 1957 H2N2 and 1968 H3N2 [<xref ref-type="bibr" rid="CR4">4</xref>, <xref ref-type="bibr" rid="CR5">5</xref>]. In June 2009, WHO announced the first pandemic in the twenty-first century, which was caused by swine-origin H1N1 (pH1N1) and resulted in more than 250,000 deaths worldwide [<xref ref-type="bibr" rid="CR6">6</xref>, <xref ref-type="bibr" rid="CR7">7</xref>]. However, the pH1N1 virus replaced the seasonal H1N1 virus, and is currently co-circulating with H3N2 and influenza B viruses.</p><p id="Par12">Antigenic drift mutations are typically localized at the antigenic sites surrounding the receptor binding site (RBS) of the HA protein. Therefore, in addition to contributing to HA antigenicity, they may affect receptor preference and virulence [<xref ref-type="bibr" rid="CR8">8</xref>]. Four antigenic sites (Ca, Cb, Sa, Sb) have been identified for H1N1 viruses compared to five sites (A-E) for H3N2 viruses [<xref ref-type="bibr" rid="CR9">9</xref>, <xref ref-type="bibr" rid="CR10">10</xref>]. Changes in HA and NA proteins may sometimes lead to acquisition of glycosylation sites, which is believed to efficiently generate antigenic variants and play critical role in virus evolution [<xref ref-type="bibr" rid="CR11">11</xref>, <xref ref-type="bibr" rid="CR12">12</xref>] by masking epitopes, stabilizing HA structures and altering receptor binding specificity [<xref ref-type="bibr" rid="CR13">13</xref>].</p><p id="Par13">Each season, the World Health Organization (WHO) select virus strains for optimal vaccine performance with the recommendation for the northern hemisphere announced in February and that for the southern hemisphere in September [<xref ref-type="bibr" rid="CR14">14</xref>]. Following 2009 pandemic, unprecedented surveillance programs have been launched to monitor influenza viruses at the molecular and epidemiological levels [<xref ref-type="bibr" rid="CR15">15</xref>, <xref ref-type="bibr" rid="CR16">16</xref>].</p><p id="Par14">Countries in the Middle East and North Africa (MENA) region are continuously considered as one of the hotspot regions for emerging and reemerging infectious diseases [<xref ref-type="bibr" rid="CR17">17</xref>]. A major reason is the location of MENA countries in the pathway of three major migratory bird flyways: east Africa–west Asia, Black Sea–Mediterranean, and east Atlantic [<xref ref-type="bibr" rid="CR18">18</xref>]. Additional risk factors that contribute in the emergence and rapid spread of epidemic diseases in the MENA region include prolonged humanitarian crises which subsequently result in massive population mobility, weak surveillance systems and limited laboratory diagnostic capacity. Further, some countries in this region, such as the Gulf States, are hubs for expatriates that constitute more than half of its population, many of which arrive from East and South-East Asia [<xref ref-type="bibr" rid="CR19">19</xref>]. Globally shared factors also include climate change and increased human–animal interaction [<xref ref-type="bibr" rid="CR20">20</xref>].</p><p id="Par15">Until recently, little was known about epidemiology of influenza in MENA region. However, many countries have developed a systemic influenza surveillance programs after 2009 pandemic [<xref ref-type="bibr" rid="CR21">21</xref>]. Currently, fourteen countries in MENA region are reporting to WHO virological surveillance database ‘FluNet; on weekly basis. However, there is a relatively limited number of genetic sequences available for influenza virus from MENA countries. Moreover, most of publications characterizing the molecular evolution of influenza virus in MENA originated from very few countries, mainly Egypt and Iran [<xref ref-type="bibr" rid="CR22">22</xref>] . In the current study, we report on the epidemiology of co-circulating influenza A and B viruses in MENA since April 2009 until December 2017. Most of countries in this region are characterized by a subtropical climate [<xref ref-type="bibr" rid="CR23">23</xref>]. However, climatic parameters such as temperature and humidity are variable among different regions of MENA and hence could affect influenza seasonality patterns. Therefore, we evaluated the epidemiology of influenza viruses among the three areas of MENA region: Gulf States, Levant, and North Africa. Additionally, we analyzed the genetic characteristics and molecular evolution of HA and NA genes of circulating pH1N1 and H3N2influenza A subtypes.</p></sec><sec id="Sec2"><title>Methods</title><sec id="Sec3"><title>Data collection</title><p id="Par16">For this study, we have used virological surveillance data available on WHO FluNet database (<ext-link ext-link-type="uri" xlink:href="http://www.who.int/flunet">http://www.who.int/flunet</ext-link>) for investigating temporal and geographical distribution of influenza types A and B in MENA countries [<xref ref-type="bibr" rid="CR24">24</xref>]. All cases reported from sixteen counties since April 2009 to December 2017 (total of 563,087 cases) were included in this study, however, no data was available from Saudi Arabia, Yemen or Libya. For further analysis on molecular evolution of influenza A subtypes (pH1N1 and H3N2), we fetched all HA and NA sequences of both subtypes depicted in Influenza Research Database (IRD) (<ext-link ext-link-type="uri" xlink:href="http://www.fludb.org">http://www.fludb.org</ext-link>) by MENA countries during 2009–2017. A total of 1226 sequences were analyzed: 512 H1 sequences, 239 H3 sequences, 343 N1 sequences and 132 N2 sequences were used for sequence analysis. Influenza vaccine strains and representative strains published on the WHO collaborating center website [<xref ref-type="bibr" rid="CR25">25</xref>] were used as references. The corresponding HA and NA sequences of were obtained from NCBI and GISAID EpiFlu Database (global initiative on sharing all influenza data (<ext-link ext-link-type="uri" xlink:href="http://www.gisaid.org">http://www.gisaid.org</ext-link>). The accession numbers for vaccine stains of HA and NA genes used in this study are summarized in Additional file <xref rid="MOESM1" ref-type="media">1</xref>: Table S1. We acknowledge the authors and submitting laboratories to all databases used in this study.</p></sec><sec id="Sec4"><title>Sequence and phylogenetic analysis</title><p id="Par17">For molecular evolution analysis, sequences were compiled and edited using the DNASTAR Lasergene sequence analysis software (DNASTAR Inc. version 13.0). Multiple alignments for pH1N1 and H3N2 HA and NA sequences were carried out using ClustalW Algorithm. Herein, H1 numbering was used for H1 sequence analysis and H3 numbering scheme was used for H3 analysis (excluding signal peptide). N1 and N2 numbering were also applied for N1 and N2 sequence analysis respectively [<xref ref-type="bibr" rid="CR26">26</xref>]. Phylogenetic trees were constructed in MEGA7, using best-fitting nucleotide substitution model, maximum likelihood method with the General time reversible model (GTR) and gamma distribution [<xref ref-type="bibr" rid="CR27">27</xref>] . The statistical significance of the tree topology was evaluated by bootstrapping with 1000 iterations. Bootstrap probabilities higher than 50 were reported at the main branches to show support values. Clustering patterns were determined using representative influenza reference strains published on the website of the WHO Collaborating Centre, London, UK (WHO- Influenza-Centre, 2009–2016) and the corresponding HA sequences were obtained from GISAID (global initiative on sharing all influenza data) and IRD (Influenza Research Database). In addition to investigating genetic drift of HA and NA genes, we looked at possible variations in glycosylation patterns of HA and NA proteins. Potential N-linked glycosylation sites were predicted using web-based NetNGlyc 1.0 Server (<ext-link ext-link-type="uri" xlink:href="http://www.cbs.dtu.dk/services/NetNGlyc/">http://www.cbs.dtu.dk/services/NetNGlyc/</ext-link>) using default settings. Only threshold values of > 0·5 were considered as potential glycosylation sites. Only scores crossing the default threshold of 0.5 were considered positive for potential glycosylation sites.</p></sec><sec id="Sec5"><title>Visualization of amino acid substitutions and glycosylation sites in HA protein using 3D models of influenza a subtypes</title><p id="Par18">The identified substitutions within antigenic sites were depicted in 3D structure models using CLC Genomic Workbench. The 3D crystal structures were downloaded from protein data bank. We used 3D model of A/California/04/2009 (PDB 3LZG) as a representative of pH1N1 viruses and 2005 human H3N2 virus (PDB 2YP7) as structure model of H3N2 viruses. The antigenic sites of HAs of pH1N1 and H3N2 viruses were assigned according to Igarashi et al. (2010) and Yang et al. (2011) respectively [<xref ref-type="bibr" rid="CR11">11</xref>, <xref ref-type="bibr" rid="CR28">28</xref>, <xref ref-type="bibr" rid="CR29">29</xref>].</p></sec><sec id="Sec6"><title>Prediction of mutations effect on HA function</title><p id="Par19">For further analysis of reported mutations, we used PROVEAN software tool to predict the possible effect of identified amino acid substitution on the biological function of a HA protein (<ext-link ext-link-type="uri" xlink:href="http://provean.jcvi.org">http://provean.jcvi.org</ext-link>). This program determines the degree of amino acid conservation in comparison to other published sequences and provides a score defining the possible effect of the substitution on protein function. A default cutoff score of less than − 2.5 indicates high probability of deleterious mutations than can affect protein’s function [<xref ref-type="bibr" rid="CR30">30</xref>].</p></sec><sec id="Sec7"><title>Estimation of evolution rate and positive selection</title><p id="Par20">The Evolutionary rate of HA and NA genes was estimated using the Bayesian Markov Chain Monte Carlo (MCMC) method as implemented in the Bayesian Evolutionary Analysis Sampling Trees (BEAST) software package v1.8.4 (<ext-link ext-link-type="uri" xlink:href="http://beast.bio.ed.ac.uk">http://beast.bio.ed.ac.uk</ext-link>) [<xref ref-type="bibr" rid="CR31">31</xref>]. Dates of virus isolation (included in sequence names) were used to calibrate the molecular clock. For nucleotide substitution model selection, we used the recommended HKY model with gamma distributed rates. A strict molecular clock model was used with a chain length of 100,000,000-250,000,000 (relative to number of sequences) logged every 5000 states. Multiple runs were performed and combined when the number of sequences analyzed exceeded 250 sequences (H1 and N1). Results were evaluated by co-efficient of variance and 95% highest posterior density (HPD) values using the tracer v1.6 program [<xref ref-type="bibr" rid="CR31">31</xref>]. Mutations rates were estimated for both HA and NA genes of pH1N1 and H3N2 viruses. Gene- and Site-specific selection pressures for HA and NA genes were estimated as the ratio of non-synonymous (dN) to synonymous (dS) nucleotide substitutions per site per year. Gene selective pressure was estimated by Tajima test of neutrality implemented in MEGA7. Positive selection was considered when D value is found to be positive (greater than 0). The test compares the average number of nucleotide differences between pairs of sampled sequences (referred to as pairwise difference-) and the total number of polymorphic sites (segregating site- S) in the sampled DNA sequences. The difference in the expectations for these two variables (which can be positive or negative) defines the Tajima’s D test statistic. A negative Tajima’s D indicates purifying selection while positive D signifies selection. Codon selective pressure was determined using conservative single-likelihood ancestor counting (SLAC) and fixed-effects likelihood (FEL) methods to identify codon sites under positive selection in selected HA and NA gene sequences. In all cases, positively selected sites were defined when <italic>p</italic> value of dN/dS ratio is less than 0.05 using the substitution model selected by website. Both methods are available at the DataMonkey online version of the Hypothesis Testing Using Phylogenies (HYPhy) package (<ext-link ext-link-type="uri" xlink:href="http://www.datamonkey.org">http://www.datamonkey.org</ext-link>) [<xref ref-type="bibr" rid="CR32">32</xref>].</p></sec><sec id="Sec8"><title>Statistical analysis</title><p id="Par21">One-Way ANOVA followed by Dunn’s multiple comparison test was performed using GraphPad in Prism version7 software to perform statistical analysis on influenza activity of the three influenza viruses among Gulf, Levant and North Africa regions.</p></sec></sec><sec id="Sec9"><title>Results</title><sec id="Sec10"><title>Temporal distribution of influenza</title><p id="Par22">According to WHO, MENA region fall within West Asia and North Africa influenza transmission zones [<xref ref-type="bibr" rid="CR33">33</xref>]. Based on this classification, countries within the same geographical group show similar influenza transmission patterns. However, the spread and geographical location of MENA countries between tropical and temperate zones may affect seasonality patterns among MENA regions [<xref ref-type="bibr" rid="CR34">34</xref>]. More specifically, climatic factors such as humidity, temperature and rainfalls have been shown to affect the shape of influenza seasonality in temperate zone compared to a less defined seasonality pattern in tropical zone [<xref ref-type="bibr" rid="CR23">23</xref>]. Herein, we divided MENA region into three sub-regions: Arabian Gulf, Levant and North Africa based on differences in climate and geographical location [<xref ref-type="bibr" rid="CR35">35</xref>] .</p><p id="Par23">A total number of 563,087 samples were reported to WHO FluNet from sixteen countries in MENA region during January 2009 and December 2017 (Table <xref rid="Tab1" ref-type="table">1</xref>). Overall, influenza positive cases accounted for 23% of all respiratory infections, with pH1N1 being the dominant virus (50.8%) followed by influenza B and H3N2 with 17.2 and 15.9% respectively (Fig. <xref rid="Fig1" ref-type="fig">1</xref>a). The remaining influenza cases were either not subtyped or subtyped as H5 (Fig. <xref rid="Fig1" ref-type="fig">1</xref>b). Unsurprisingly, the largest number of influenza positive cases were reported during 2009 H1N1 pandemic (33,064 cases) while the smallest number of positive cases were reported in 2012 (4918 cases) (Fig. <xref rid="Fig1" ref-type="fig">1</xref>b). During study period, pH1N1 viruses dominated in all years except in 2012 during which H3N2 viruses dominated (Table <xref rid="Tab1" ref-type="table">1</xref>; Fig. <xref rid="Fig1" ref-type="fig">1</xref>b). After the 2009 pandemic, pH1N1 cases increased rapidly reaching up to 65.3% of influenza positive cases in 2009 thereafter peaked to highest amplitude in 2010 (72%). During 2011, however, the percentage of pH1N1 positive cases declined to 37.5%, while prevalence of H3N2 and influenza B increased gradually accounting to 30.9 and 22% of positive cases respectively (Fig. <xref rid="Fig1" ref-type="fig">1</xref>b). In 2013, influenza A was still the dominant type, with H3N2 being the dominant subtype representing 48.3% of the case, and only 13% were documented as pH1N1. Major activities were reported for influenza B during 2011–2012 reaching around 30%. Such activity was stable in the subsequent years. Interestingly, pH1N1 and H3N2 were showing opposite activities during following years, although pH1N1 continued to be the dominant virus (Fig. <xref rid="Fig1" ref-type="fig">1</xref>b).<table-wrap id="Tab1"><label>Table 1</label><caption><p>Total number of influenza A positive samples and number of HA/NA sequences analyzed in this study. Influenza surveillance data submitted from 14 countries in MENA region (2009–2017) were obtained from WHO FluNet online database. For evolution analysis, complete HA and NA sequences from pH1N1 and H3N2 viruses were downloaded from Influenza Research Database (IRD)-NCBI. Accession numbers of all vaccine strains are provided in Additional file <xref rid="MOESM1" ref-type="media">1</xref>: Table S1</p></caption><table frame="hsides" rules="groups"><thead><tr><th rowspan="3">Year</th><th rowspan="3">Dominant Subtype</th><th colspan="4">pH1N1</th><th colspan="4">H3N2</th></tr><tr><th rowspan="2">Number of positive samples (WHO)</th><th colspan="2">Number of sequences analyzed (NCBI)</th><th rowspan="2">Vaccine strain</th><th rowspan="2">Number of positive samples (WHO)</th><th colspan="2">Number of sequences analyzed (NCBI)</th><th rowspan="2">Vaccine strain</th></tr><tr><th>H1</th><th>N1</th><th>H3</th><th>N2</th></tr></thead><tbody><tr><td><bold>2009</bold></td><td>pH1N1</td><td>21,586</td><td>132</td><td>97</td><td>A/California/07/2009</td><td>1093</td><td>22</td><td>5</td><td>A/Brisbane/10/2007</td></tr><tr><td><bold>2010</bold></td><td>pH1N1</td><td>4136</td><td>68</td><td>42</td><td>A/California/07/2009</td><td>280</td><td>1</td><td>1</td><td>A/Brisbane/10/2007</td></tr><tr><td><bold>2011</bold></td><td>pH1N1</td><td>3717</td><td>36</td><td>36</td><td>A/California/07/2009</td><td>2183</td><td>22</td><td>22</td><td>A/Perth/16/2009</td></tr><tr><td><bold>2012</bold></td><td>H3N2</td><td>660</td><td>27</td><td>29</td><td>A/California/07/2009</td><td>2373</td><td>37</td><td>32</td><td>A/Perth/16//2009</td></tr><tr><td><bold>2013</bold></td><td>pH1N1</td><td>4395</td><td>33</td><td>32</td><td>A/California/07/2009</td><td>1568</td><td>13</td><td>8</td><td>A/Victoria//361/2011</td></tr><tr><td><bold>2014</bold></td><td>pH1N1</td><td>3559</td><td>25</td><td>25</td><td>A/California/07/2009</td><td>2883</td><td>40</td><td>33</td><td>A/Texas/50/2012</td></tr><tr><td><bold>2015</bold></td><td>pH1N1</td><td>9742</td><td>157</td><td>79</td><td>A/California/07/2009</td><td>1688</td><td>47</td><td>23</td><td>A/Switzerland/9715293/2013</td></tr><tr><td><bold>2016</bold></td><td>pH1N1</td><td>8320</td><td>60</td><td>2</td><td>A/California/07/2009</td><td>4784</td><td>47</td><td>8</td><td>A/Hong Kong/4801/2014</td></tr><tr><td><bold>2017</bold></td><td>pH1N1</td><td>10,143</td><td>–</td><td>–</td><td>A/Michigan/45/2015</td><td>3891</td><td>10</td><td>–</td><td>A/Hong Kong/4801/2014</td></tr><tr><td><bold>Total</bold></td><td></td><td>66,228</td><td>538</td><td>342</td><td></td><td>20,743</td><td>239</td><td>132</td><td></td></tr></tbody></table></table-wrap><fig id="Fig1"><label>Fig. 1</label><caption><p><bold>a</bold> Percentage of influenza types/subtypes reported to FluNet database from MENA (2009–2017). The bar chart represents the contribution of each type/subtype to the overall percentage of positive samples. <bold>b</bold> Overall trends of influenza activity in MENA region during 2009–2017 as reported to FluNet database. Total number of influenza positive samples received per year is presented by gray bars. Annual percentages of reported cases positive for influenza A subtypes pH1N1 (red), H3N2 (blue) and influenza B viruses are also presented</p></caption><graphic xlink:href="12879_2019_3930_Fig1_HTML" id="MO1"></graphic></fig></p><p id="Par24">The mean proportions of influenza positive cases reported during 2009–2017 had shown no significant differences among Gulf, Levant and North Africa regions (<italic>P</italic> value 0.39). More specifically, the mean proportion of pH1N1 and influenza B cases did not differ significantly (<italic>P</italic> value 0.85 and 0.988 respectively) among regions under the study. However, a significant difference was found in percentage of H3N2 cases (<italic>P</italic> value 0.0158) between Gulf and North Africa. Although the overall epidemiological trend of influenza A and B viruses were similar among the three regions, minor variations were observed after 2011 (Fig. <xref rid="Fig2" ref-type="fig">2</xref>). Although the activity of pH1N1 was the lowest among other viruses during 2012 in Levant and North Africa, all three viruses were reported to have similar prevalence in Gulf during 2012. During following years, variable influenza activity was seen of the three viruses among the three regions (Fig. <xref rid="Fig2" ref-type="fig">2</xref>a). In Levant, pH1N1 and H3N2 viruses were showing alternative activities with pH1N1 dominating in 2013, 2015 and 2016 while H3N2 dominating in 2012, 2014 and 2017 (Fig. <xref rid="Fig2" ref-type="fig">2</xref>b). On the other hand, pH1N1 remained the dominant virus during 2013–2017 in both Gulf and North Africa region (Fig. <xref rid="Fig2" ref-type="fig">2</xref>b and d). Influenza B showed similar activity trends to that of H3N2 in Gulf and Levant except in 2013 and 2015, during which, influenza B was showing similar trend to that of pH1N1. However, opposite activity patterns of the influenza B and H3N2 viruses were observed in North Africa during the same period (2013–2017) (Fig. <xref rid="Fig2" ref-type="fig">2</xref>c). Of note, pH1N1 positive samples dropped by half in Levant (18%) and North Africa (21%) regions in 2017 while increased to 59% in Gulf (Fig. <xref rid="Fig2" ref-type="fig">2</xref>b-d).<fig id="Fig2"><label>Fig. 2</label><caption><p>Influenza activity patterns of co-circulating influenza A subtypes (pH1N1 and H3N2) and influenza B in MENA, 2009–2017. <bold>a</bold> Heat map showing the levels of influenza activity of each type/subtype in each region (<italic>n</italic> = 3) from April 2009 to December 2017 period. The color scale denotes the percentage of influenza positive samples relative to corresponding color shade. Annual activity trends of influenza viruses: pH1N1 (red), H3N2 (blue) and influenza B in Gulf (<bold>b</bold>), Levant (<bold>c</bold>) and North Africa (<bold>d</bold>) regions since April 2009 to December 2017. Bars represent the total number of samples collected in each year while doted pattern in each bar denotes number of influenza positive samples in each specific year</p></caption><graphic xlink:href="12879_2019_3930_Fig2_HTML" id="MO2"></graphic></fig></p><p id="Par25">The well-established trend of influenza A was similar to that of Northern Hemisphere with a single consistent peak in all epidemics of all three regions (Additional file <xref rid="MOESM2" ref-type="media">2</xref>: Figure S1) [<xref ref-type="bibr" rid="CR36">36</xref>]. The peak period, however, commenced in October in Gulf and North Africa while started later in December in Levant. During epidemics, influenza-positive cases rapidly increased from October and peaked in December–January in Gulf and North Africa. In contrast, highest influenza activity was reported during January–February in Levant. Additionally, the duration of influenza activity was similar in Gulf and North Africa; in which the activity continued for a longer period (October–April) compared to a sharper peak in Levant (December–March). Interestingly, the typical seasonal influenza peak was altered during specific epidemics. In North Africa, for example, the influenza activity peaked during February–March in 2013 and 2014 epidemics instead of the classical peak in December–January. Moreover, a late secondary peak with high amplitude (77%) was observed in Levant during February–May in 2015 following an early primary peak (amplitude of 17%) in December. Finally, a stable influenza activity was noticed during December to April of 2014 and 2015 in Gulf with no distinctive primary peak observed (Additional file <xref rid="MOESM2" ref-type="media">2</xref>: Figure S1).</p></sec><sec id="Sec11"><title>Molecular evolution analysis of HA and NA genes</title><sec id="Sec12"><title>Evolution rates and positive selection</title><p id="Par26">The mean evolution rate of the H1 and H3 genes of pH1N1 and H3N2 viruses analyzed in this study amounted to 1.91 × 10<sup>− 3</sup> and 1.47 × 10<sup>− 3</sup> substitutions per site per year respectively. For NA genes, values of 1.92 × 10<sup>− 3</sup> and 3.09 × 10<sup>− 3</sup> substitutions per site per year were obtained for the N1 and N2 respectively (Table <xref rid="Tab2" ref-type="table">2</xref>). Selective pressure acting on HA and NA genes showed that both genes of pH1N1 and H3N2 were under negative selection (<italic>P</italic> value of pH1N1: HA = − 2.13 and NA = − 2.19) and (P value of H3N2: HA = − 1.92 and NA = − 1.95). Besides, the positive selection at codon level was estimated by assessing the overall number of non-synonymous (dN) to synonymous substitutions (dS). Analysis showed that the overall dN/dS values of the coding HA regions of pH1N1 and H3N2 were 0.32 and 0.28, respectively (Table <xref rid="Tab2" ref-type="table">2</xref>). Two sites in H1 gene were found to be under positive selection by two computational models (SLAC and FEL) with significant statistical <italic>p</italic>-values (< 0.05 for both). The two amino acids which were found to be prone to adaptive evolution were 239 and 240 in Ca antigenic site. These two residues are located in an overlapping position of antigenic and receptor binding site [<xref ref-type="bibr" rid="CR29">29</xref>]. In particular, amino acid in 239 position, is a key determinant for virus binding to sialylated glycan receptors [<xref ref-type="bibr" rid="CR37">37</xref>]. Additionally, amino acids in 202 and 468 positions corresponding to receptor binding site and stem respectively, were also found to be under positive selection using FEL method. Two codons in the signal peptide (positions 10 and 13) and three codons in stem were additionally identified to be under diversifying positive selection. On the other hand, applying FEL and SLAC methods on H3 genes resulted in different outputs regarding residues under positive selection. While 239 position was found to be under positive selection using SLAC method, position 16 and 69 were determined to be positively selected by FEL method (<italic>p</italic> value < 0.05). None of residues of NA genes of both subtypes were reported to be under positive pressure.<table-wrap id="Tab2"><label>Table 2</label><caption><p>Estimation of evolution rate and positive selection for HA and NA genes of pH1N1 and H3N2 viruses. Evolution rates were estimated using the Bayesian Markov Chain Monte Carlo method in the Bayesian Evolutionary Analysis Sampling Trees (BEAST) package (v1.8.4). HPD represents lower and upper limits of highest posterior density that contains 95% of the sampled values. The selective pressure at codon level was also calculated using the most conservative ‘single-likelihood ancestor counting’ algorithm (SLAC) and the ‘fixed effects likelihood’ (FEL) methods available at the DataMonkey web server of the HYPHY software. (N) refers to the number of positively selected amino acids and (AA) refers to amino acids that were found to be under positive selection. Positively selected sites were defined when dN/dS p value is less than 0.05 using the best fit model as recommended by website</p></caption><table frame="hsides" rules="groups"><thead><tr><th rowspan="3"></th><th rowspan="3">Number of analysed sequences</th><th colspan="2">Nucleotide substitution rate (×10<sup>−3</sup>)</th><th colspan="6">Positive selection</th></tr><tr><th rowspan="2">Mean</th><th rowspan="2">95% HPD</th><th colspan="3">SLAC</th><th colspan="3">FEL</th></tr><tr><th>Mean (dN/dS)</th><th>N</th><th>AA</th><th>Mean (dN/dS)</th><th>N</th><th>AA</th></tr></thead><tbody><tr><td>H1</td><td>512</td><td>1.91</td><td>1.68–2.1</td><td>0.336</td><td>3</td><td>10–239-240</td><td>0.311</td><td>8</td><td>13–202–239-240-305-468-543-546</td></tr><tr><td>H3</td><td>239</td><td>1.47</td><td>1.27–1.68</td><td>0.299</td><td>1</td><td>239</td><td>0.267</td><td>2</td><td>19–69</td></tr><tr><td>N1</td><td>342</td><td>1.93</td><td>1.7–2</td><td>0.352</td><td>–</td><td>–</td><td>0.311</td><td>–</td><td>–</td></tr><tr><td>N2</td><td>132</td><td>3.09</td><td>2.5–3.1</td><td>0.238</td><td>–</td><td>–</td><td>0.208</td><td>–</td><td>–</td></tr></tbody></table></table-wrap></p></sec><sec id="Sec13"><title>Phylogenetic tree analysis</title><p id="Par27">To assess the evolution pH1N1 and H3N2 viruses during the same period, HA sequences of circulating strains were used to construct phylogenetic trees along with the vaccine and reference sequences. The phylogenetic trees of HA genes of both subtypes formed phylogenetic clusters representative of each season between 2009 and 2017. The phylogenetic analysis of HA sequences revealed clustering patterns similar to those reported in the annual and interim reports of WHO since 2009 until 2016 [<xref ref-type="bibr" rid="CR25">25</xref>]. Analysis of pH1N1 viruses revealed seven genetic groups that evolved between 2009 and 2017 (Fig. <xref rid="Fig3" ref-type="fig">3</xref>). In 2009, the majority of viruses belonged to the HA clade 1 represented by the vaccine strain A/California/07/2009, however, minority of viruses (7.6%) were also found to belong to clade 4. In contrast, pH1N1 viruses isolated during 2010 were showing more diversity by clustering within four genetic clades: clade 1 (37%), clade 3 (7%), clade 5 (10%) and clade 7 (26%). In 2011, viruses within clades 5 and 7 were predominant (50 and 28% respectively) in addition to clades 2, 3 and 6B which represents 27% of viruses in 2011. Following 2011, all pH1N1 viruses belonged to clade 6 and its related sub-clades (6A, 6B and 6C). Clade 6 is characterized by carrying amino acid substitutions: D114N, S202 T, S220 T and K300E in HA1 and E391K, S468 N and E516K in HA2. In 2012, pH1N1 viruses were mainly found within subclades 6A (44%) and 6C (26%). Viruses belonging to both clades continued to circulate during 2013, however, minority of viruses were also found to belong to sub-clade 6B which has the additional substitutions: K180Q and A273T in HA1 and E516K in HA2. The HA sequences of pH1N1 viruses isolated in 2014 clustered mainly within clade 6B (80%) and sub-clade 6B.2 (12%) which is characterized by E508G substitution in HA2. In 2015, clade 6B.1 became the predominant genetic group, while only 5% of viruses belonged to 6B.2 sub-clade. Viruses within 6B.1 group characterized by S101 N, S179 N and I233T in HA1 as compared to parental clade 6B. All Viruses collected in 2016 and 2017 (only 3 sequences found in IRD) fell into genetic clade 6B.1.<fig id="Fig3"><label>Fig. 3</label><caption><p>Phylogenetic trees of HA nucleotide sequences of pH1N1 and H3N2 viruses. HA sequences of pH1N1 (<italic>n</italic> = 512) and H3N2 (<italic>n</italic> = 239) reported from MENA region between 2009 and 2017 were compared to vaccine strains recommended by WHO (denoted with green triangles) and the reference strains of known clades as reported in WHO Influenza Center London, indicated by navy blue circles [25]. Phylogenic trees were generated using maximum likelihood method by GTR + G model using in MEGA7 software. Bootstrap values of 1000 replicates and only values larger than 50 are indicated at the nodes. At the major nodes are the signature amino acid changes in different colors according to antigenic site: Red for Sa (pH1N1) and A (H3N2), blue for Ca (pH1N1) and B (H3N2), green for C (H3N2). Scale bar represents approximately 0.005 nucleotide difference between close relatives</p></caption><graphic xlink:href="12879_2019_3930_Fig3_HTML" id="MO3"></graphic></fig></p><p id="Par28">On the other hand, phylogenetic analysis of H3 gene clustered viruses into three genetic groups in which vaccine strain A/Perth/16/2009 represent clade 1 (Fig. <xref rid="Fig3" ref-type="fig">3</xref>). H3N2 viruses in 2009 clustered with A/Perth/16/2009, the vaccine strain recommended during 2010–2012. During 2011 and 2012, H3N2 strains belonged to sub-clades of genetic group 3 (3A, 3B and 3C) which share amino acid substitutions: N145S, V223I in HA1 and D487N in HA2. Of the genetic group B, H3N2 viruses belonging to sub-clade 3B dominated in 2011 (41%) and in 2012 (54%). Moreover, 32% of sequences in 2012 fell into clade 5 which is characterized by D53N, Y94H and E280A. H3N2 strains from 2013 and 2014 grouped mainly into 3C.3 and 3C.2 sub-clades. In 2013, 38% of viruses were found within 3C.3 related subclades (3C.3a and 3C.3b), 23% within 3C.2 and 23% within 3A. Similarly, majority of viruses sequenced in 2014 belonged to sub-clade 3C.3 (65%) compared to 12.5% belonged to 3C.2 sub-clade. All H3N2 viruses of 2015 and 2016 belonged to sub-clade 3C. 2 as defined by N145S substitution. Of these, more than 90% were grouped within 3C.2a sub-clade that is characterized by additional substitutions such as L3I, N144S, K160 T and N225D. By the end of 2015, a novel subgroup of viruses emerged within 3C.2a subclade, 3C.2a1, with majority of viruses belonging to 3C.2a1 subclades carrying N121K substitution. All H3N2 viruses isolated from MENA region in 2017 were found to belong to 3C.2a1 subclade.</p></sec></sec><sec id="Sec14"><title>Amino acid sequence variation in HA and NA proteins</title><p id="Par29">Here, amino acid sequences for HA and NA proteins of pH1N1 and H3N2 viruses circulated in MENA region during 2009–2017 were compared to vaccine strains of respective years. Genetic characterization was performed by sequencing analysis of 512 HA and 343 NA sequences of pH1N1; and 239 HA and 132 NA sequences of H3N2 viruses. Full length HA and NA proteins were used for alignment and subsequent analysis: the length of HA and NA genes were 566 and 469 amino acids for both viruses respectively.</p></sec><sec id="Sec15"><title>Genetic characterization of pH1N1 viruses</title><sec id="Sec16"><title>Molecular evolution of H1 protein</title><p id="Par30">Alignment of 512 HA amino acid sequences of pH1N1 viruses revealed multiple amino acid changes in comparison to vaccine strain A/California/07/2009. Overall, amino acid sequences of full HA protein were highly conserved showing > 98% similarity to Cal09 virus in all pH1N1 viruses isolated since pandemic in 2009 until 2017. Comparison of HA antigenic sites of 2016 pH1N1 with vaccine strain is shown in Table <xref rid="Tab3" ref-type="table">3</xref>. Two substitutions at 179 and 180 positions were found in Sa antigenic site that showed 84.6% similarity to Cal09 vaccine virus. Sb and Ca sites, on the other hand, showed 91.7 and 94.7% similarities to A/California/07/2009 virus respectively during 2016 while no substitutions were detected in Cb site (Table <xref rid="Tab3" ref-type="table">3</xref>).<table-wrap id="Tab3"><label>Table 3</label><caption><p>Identification of amino acid substitutions in H1 antigenic sites of pH1N1 viruses. Antigenic sites of HA sequences were compared to vaccine strains A/California/07/2009. Amino acids in <italic>italics</italic> indicate amino acids that differ from vaccine sequence. Bold amino acids denote amino acid changes found in vaccine strain recommended for subsequent season (A/Michigan/45/2015). A/California/07/2009 was used in influenza vaccine composition since the emergence of pandemic H1N1 in 2009 before being replaced by A/Michigan/45/2015 in 2016. Similarity of antigenic sites between vaccine strains and circulating viruses is shown as percentages in the last column</p></caption><table frame="hsides" rules="groups"><thead><tr><th colspan="5">H1</th><th rowspan="3" colspan="3">Similarity to pH1N1 vaccine strain (%)</th></tr><tr><th rowspan="2"></th><th colspan="2">Sa</th><th>Sb</th><th>Ca</th></tr><tr><th>179</th><th>180</th><th>202</th><th>220</th></tr></thead><tbody><tr><td>California/2009</td><td>S</td><td>K</td><td>S</td><td>S</td><td>Sa</td><td>Sb</td><td>Ca</td></tr><tr><td>2016</td><td><italic>N</italic></td><td><italic>Q</italic></td><td><italic>T</italic></td><td><italic>T</italic></td><td>84.6</td><td>91.7</td><td>94.7</td></tr><tr><td>Michigan/2015<sup>a</sup></td><td>N</td><td>Q</td><td>T</td><td>T</td><td colspan="3"></td></tr></tbody></table><table-wrap-foot><p><sup>a</sup>A/Michigan/45/2015 vaccine was recommended to be used in 2016–2017 season, however, only 3 sequences were depicted in NCBI-IRD from MENA region during 2017 therefore sequences isolated in 2016 were used for comparison purpose</p></table-wrap-foot></table-wrap></p><p id="Par31">In total, 26 substitutions were reported in HA protein; 15 in HA globular head and 6 in the stem region. Seven of identified substitutions were found in antigenic sites including S179 N and K180Q in Sa site; S202 T in Sb site and A203T, S220 T, R222K, D239G/E/N and Q240R in Ca site (Table <xref rid="Tab3" ref-type="table">3</xref>). Substitutions were also detected sporadically near antigenic sites such as Q240R and G187E/R substitutions near the Ca site and A214T substitution near to the Sb site. Majority of amino acid residues within RBS were highly conserved showing only four substitutions: V149E and A151T in 130-loop; S202 T and A203T in the 190-helix and D239N/E/G in the 220-loop (Table <xref rid="Tab4" ref-type="table">4</xref>). Additionally, two substitutions were also observed near RBS: S200P and I233T that have been identified to alter receptor binding affinity. Ten of the amino acid substitutions were successfully introduced into the viral population in 2015, and were identified in majority of viruses (> 80%) circulating in subsequent seasons. Half of these substitutions are located in HA head (S179 N, K180Q, S220 T, I233T and A273T) while the other half are located in HA stem region (D114N, K300E, E391K, S468 N and E516K).<table-wrap id="Tab4"><label>Table 4</label><caption><p>Amino acid substitutions identified in HAs of pH1N1 viruses sequenced during 2009–2017. Amino acid substitutions of pH1N1 are reported in comparison to A/California/07/2009. Bordered cells indicate substitutions in receptor binding site. Deleterious mutations were predicted by online server PROVEAN and are indicated in bold</p></caption><graphic xlink:href="12879_2019_3930_Tab4_HTML" id="MO4"></graphic><table-wrap-foot><p>*Substitution Q240R was found in 16 viruses isolated exclusively from Iran</p></table-wrap-foot></table-wrap></p><p id="Par32">In 2009, the majority of viruses (68%) emerged with S220 T substitution at Ca site. Additional substitutions were also reported with moderate frequencies such as D239E (21%), P314S (17%) and E391K (12%) while only minority of viruses were carrying D114N (3%), A203T (2%) and R222K (4%) substitutions in respective season (Table <xref rid="Tab4" ref-type="table">4</xref>). Of these substitutions, S220 T and E391K as well as D114N were the genetic signature of subsequently circulating viruses. Viruses isolated in 2010 were characterized by the emergence of S160G (23%), S202 T (26%), A214T (26%) substitutions. Additionally, a minor group of viruses had N104K (10%), S200P (7%) and N458K (6%) substitutions. The later substitution (N458K) was predicted by PROVEAN to exert deleterious effect on HA1 function (score − 4.98). Further, two of previously circulating variants, D239E (22%) and P314S (23%), were lost from viruses isolated after 2010. Two of detected substitutions: G105 W (12%) and S138I (6%) were introduced in the viral population and were only found in viruses isolated in 2010 (Fig. <xref rid="Fig4" ref-type="fig">4</xref>a). During 2011, no additional variants were encountered in pH1N1 viruses, instead, previously reported mutations kept circulating with variable frequencies. Forty-seven percent of the strains had R222K substitution at the Ca site and 33% had A214T substitution. S200P substitution was also observed in minority of viruses (5%). Importantly, 8 and 3% of circulating pH1N1 strains in 2011were detected with D239G and N458K mutations respectively, both of which were predicted to be have deleterious effect on HA1 by Provean. S160G, S202 T and N458K were not detected in subsequent influenza seasons in MENA region. Despite the minor pH1N1 activity during 2012, new viruses evolved carrying additional substitutions: E516K and V149A (11 and 50% respectively) in addition to K300E which was encountered in more than 50% of viruses. Along with E516K, K300E substitution in stem domain became fixed in subsequent seasons. All viruses isolated during 2012 were also carrying D114N, S220 T, E391K and S468 N substitutions. During 2013 season, all viruses were characterized by S220 T/E391K/E516K changes (Fig. <xref rid="Fig4" ref-type="fig">4</xref>a). Two additional substitutions, K180Q and A273T, was detected in in 9 and 6% of viruses respectively during 2013 season, however both became characteristic for the subsequent pH1N1 viruses isolated during following seasons (Fig. <xref rid="Fig4" ref-type="fig">4</xref>a). The deleterious mutation D239G was also detected at low frequency (6%) in addition to previously circulated mutations such as S200P (3%), A203T (12%) and A214T (12%). As 2014 year started, five substitutions were found in 100% of viruses isolated in subsequent seasons: D114N, S220 T, E391K, S468 N and E516K. Furthermore, 16% of viruses in this season had E508G substitution in the HA2 subunit (Table <xref rid="Tab4" ref-type="table">4</xref>). In 2015, the majority of viruses were characterized by two novel variants, S179 N (87%) and I233T (87%), which gradually replaced original amino acids and became dominant variants in 2016. Viruses circulated during 2016 season were more stable compared to previous seasons with no novel variants detected. Instead, almost all viruses isolated in this season (95–100%) were carrying fixed mutations that has been established in previous seasons such as D114N, S179 N, K180Q and I233T in HA1 subunit and E391K, N468S and E516K in HA2 subunit (Fig. <xref rid="Fig4" ref-type="fig">4</xref>a).<fig id="Fig4"><label>Fig. 4</label><caption><p>Summary of amino acid substitutions identified in HA proteins of (A) pH1N1 viruses and (B) H3N2 viruses between 2009 and 2017 in MENA. Substitutions are colored with respect to their localization in the antigenic sites. For pH1N1: Sa is colored red and Ca is colored blue. For H3N2: A is colored red; B is colored blue; C is colored green and D is colored orange. Sporadic and fixed substitutions (bold) are indicated by the dashed and continuous lines respectively. Fixed substitutions that have been found in > 80% of viruses are listed in parallel (on the right) to respective year. Only one HA sequence has been found in IRD representing H3N2 virus from MENA during 2010 (red star)</p></caption><graphic xlink:href="12879_2019_3930_Fig4_HTML" id="MO5"></graphic></fig></p><p id="Par33">Glycosylations of HA protein can also alter antigenic properties of the virus by masking antigenic sites and preventing antibodies recognition [<xref ref-type="bibr" rid="CR38">38</xref>]. Seven potential glycosylation sites were found in the HA protein of pH1N1 strains, including N28, N40, N104, N304, N498 and N557 which are conserved in all H1 sequences (Table <xref rid="Tab5" ref-type="table">5</xref>) [<xref ref-type="bibr" rid="CR39">39</xref>]. Interestingly, the appearance of N179 residue at Sa site in 2015 has led to acquisition of a novel potential glycosylation site at head of HA protein (Fig. <xref rid="Fig5" ref-type="fig">5</xref>). Potential glycosylation at N179 site located in HA globular head was first detected in 2% of HA sequences during 2011. However, 85 and 99% of HA sequences have acquired N179 glycosylation sequon during 2015 and 2016, respectively (Table <xref rid="Tab5" ref-type="table">5</xref>). An additional glycosylation site, N489, was detected in sequences isolated from Iran during 2013 (18%) and 2016 (2%) (Table <xref rid="Tab5" ref-type="table">5</xref>). Of note, glycosylation sequons at N28, N104, N498 and N557 were lost from viruses isolate from Iran at different time periods. In 2015, for example, N28 and N557 sequons were missing from 1.5 and 5% of HA sequences, respectively. N104 glycosylation was also missing from 2% of HA sequences, all of which have been isolated from Iran during 2010, 2013 and 2016.<table-wrap id="Tab5"><label>Table 5</label><caption><p>Dynamics of HA and NA glycosylation for pH1N1 viruses isolated from MENA region during 2009–2017. In total, glycosylation patterns of 512 HA and 343 NA sequences of pH1N1 virus were analyzed. Numbers between brackets represent the percentage of sequences expressing glycosylation sequon in corresponding year while conserved glycosylation sites are indicated in bold</p></caption><table frame="hsides" rules="groups"><thead><tr><th colspan="13">Localization of potential glycosylation sequons in HA and NA of pH1N1 viruses</th></tr></thead><tbody><tr><td colspan="13">HA</td></tr><tr><td></td><td></td><td colspan="11">Amino acid position</td></tr><tr><td>Year</td><td>Number of sequences analyzed</td><td colspan="2">Stem</td><td colspan="3">Head</td><td colspan="6">Stem</td></tr><tr><td>2009</td><td>132</td><td><bold>28</bold></td><td><bold>40</bold></td><td colspan="2"><bold>104</bold></td><td></td><td colspan="2"><bold>304</bold></td><td></td><td><bold>498</bold></td><td colspan="2"><bold>557</bold></td></tr><tr><td>2010</td><td>68</td><td><bold>28</bold></td><td><bold>40</bold></td><td colspan="2"><bold>104</bold></td><td></td><td colspan="2"><bold>304</bold></td><td></td><td><bold>498</bold></td><td colspan="2"><bold>557</bold></td></tr><tr><td>2011</td><td>36</td><td><bold>28</bold></td><td><bold>40</bold></td><td colspan="2"><bold>104</bold></td><td>179 (2%)</td><td colspan="2"><bold>304</bold></td><td></td><td><bold>498</bold></td><td colspan="2"><bold>557</bold></td></tr><tr><td>2012</td><td>27</td><td><bold>28</bold></td><td><bold>40</bold></td><td colspan="2"><bold>104</bold></td><td></td><td colspan="2"><bold>304</bold></td><td></td><td><bold>498</bold></td><td colspan="2"><bold>557</bold></td></tr><tr><td>2013</td><td>33</td><td><bold>28</bold></td><td><bold>40</bold></td><td colspan="2"><bold>104</bold></td><td></td><td colspan="2"><bold>304</bold></td><td>489* (18%)</td><td><bold>498</bold></td><td colspan="2"><bold>557</bold></td></tr><tr><td>2014</td><td>25</td><td><bold>28</bold></td><td><bold>40</bold></td><td colspan="2"><bold>104</bold></td><td></td><td colspan="2"><bold>304</bold></td><td></td><td><bold>498</bold></td><td colspan="2"><bold>557</bold></td></tr><tr><td>2015</td><td>157</td><td><bold>28</bold></td><td><bold>40</bold></td><td colspan="2"><bold>104</bold></td><td>179 (85%)</td><td colspan="2"><bold>304</bold></td><td></td><td><bold>498</bold></td><td colspan="2"><bold>557</bold></td></tr><tr><td>2016</td><td>60</td><td><bold>28</bold></td><td><bold>40</bold></td><td colspan="2"><bold>104</bold></td><td>179 (99%)</td><td colspan="2"><bold>304</bold></td><td>489* (2%)</td><td><bold>498</bold></td><td colspan="2"><bold>557</bold></td></tr><tr><td colspan="13">NA</td></tr><tr><td>Year</td><td>Number of sequences analyzed</td><td colspan="11">Amino acid position</td></tr><tr><td>2009</td><td>97</td><td></td><td></td><td char="(" align="char">44 (11%)</td><td><bold>50</bold></td><td><bold>58</bold></td><td char="." align="char"><bold>63</bold></td><td><bold>68</bold></td><td><bold>88</bold></td><td><bold>146</bold></td><td><bold>235</bold></td><td><bold>386</bold></td></tr><tr><td>2010</td><td>42</td><td></td><td>42 (36%)</td><td></td><td><bold>50</bold></td><td><bold>58</bold></td><td char="." align="char"><bold>63</bold></td><td><bold>68</bold></td><td><bold>88</bold></td><td><bold>146</bold></td><td><bold>235</bold></td><td><bold>386</bold></td></tr><tr><td>2011</td><td>36</td><td></td><td>42 (100%)</td><td></td><td><bold>50</bold></td><td><bold>58</bold></td><td char="." align="char"><bold>63</bold></td><td><bold>68</bold></td><td><bold>88</bold></td><td><bold>146</bold></td><td><bold>235</bold></td><td><bold>386*</bold></td></tr><tr><td>2012</td><td>29</td><td>28 (7%)</td><td>42 (100%)</td><td></td><td><bold>50</bold></td><td><bold>58</bold></td><td char="." align="char"><bold>63</bold></td><td><bold>68</bold></td><td><bold>88</bold></td><td><bold>146</bold></td><td><bold>235</bold></td><td><bold>386</bold></td></tr><tr><td>2013</td><td>32</td><td></td><td>42 (100%)</td><td></td><td><bold>50</bold></td><td><bold>58</bold></td><td char="." align="char"><bold>63</bold></td><td><bold>68</bold></td><td><bold>88</bold></td><td><bold>146</bold></td><td><bold>235</bold></td><td><bold>386</bold></td></tr><tr><td>2014</td><td>25</td><td></td><td>42 (100%)</td><td></td><td><bold>50</bold></td><td><bold>58</bold></td><td char="." align="char"><bold>63</bold></td><td><bold>68</bold></td><td><bold>88</bold></td><td><bold>146</bold></td><td><bold>235</bold></td><td><bold>386</bold></td></tr><tr><td>2015</td><td>79</td><td></td><td>42 (100%)</td><td></td><td><bold>50</bold></td><td><bold>58</bold></td><td char="." align="char"><bold>63</bold></td><td><bold>68</bold></td><td><bold>88</bold></td><td><bold>146</bold></td><td><bold>235</bold></td><td><bold>386</bold></td></tr><tr><td>2016</td><td>3</td><td></td><td>42 (100%)</td><td></td><td><bold>50</bold></td><td><bold>58</bold></td><td char="." align="char"><bold>63</bold></td><td><bold>68</bold></td><td><bold>88</bold></td><td><bold>146</bold></td><td><bold>235</bold></td><td></td></tr></tbody></table><table-wrap-foot><p>HA Glycosylation: glycosylation site, N489, was detected in sequences isolated from Iran. Glycosylation sequons at N28, N104, N498 and N557 were lost from viruses isolate from Iran at different time periods. In 2015, for example, N28 and N557 sequons were missing from 1.5 and 5% of HA sequences, respectively. NA Glycosylation: N386 glycosylation sequon was missing from one sequence in 2011, however, during following years more NA sequences were found to lack 386 glycosylation (31% of sequences in 2012; 80% in 2014; 77% in 2015) before disappearing completely from NA sequences in 2016</p></table-wrap-foot></table-wrap><fig id="Fig5"><label>Fig. 5</label><caption><p>The glycosylation sites of HA proteins of pH1N1 and H3N2 viruses. The 3D structures of HA monomers from pH1N1 (PDB 3LZG) and H3N2 (PDB 2YP7) are displayed using CLC genomic workbench version 11. The 3D structures show glycosylation sites with respect to their localization in the HA protein: conserved glycosylation sequon are colored green and acquired glycosylation sequon are colored in red</p></caption><graphic xlink:href="12879_2019_3930_Fig5_HTML" id="MO6"></graphic></fig></p></sec><sec id="Sec17"><title>Molecular evolution of N1 protein</title><p id="Par34">A total of 343 NA protein sequences of pH1N1 viruses were analyzed in comparison to NA of vaccine strain Cal09 (KC781784). The catalytic sites represented by the following residues: 118R, 119E, 151D, 152R, 179 W, 223I, 225R, 277E, 293R, 368R, and 402Y as well as framework residues which support the catalytic residues: 156R, 180S, 228E, 247S, 278E, and 295 N [<xref ref-type="bibr" rid="CR40">40</xref>, <xref ref-type="bibr" rid="CR41">41</xref>] were highly conserved in all pH1N1 viruses. However, I117M (3%), N270K (10%) and N369K (47%) substitutions were found near the catalytic sites, whereas N248D (93%) and N270K (10%) substitutions were encountered close to the framework sites (Additional file <xref rid="MOESM3" ref-type="media">3</xref>: Table S2). Fourteen distinct substitutions were observed at NA antigenic sites denoted by NA residues 83–143, 156–190, 252–303, 330, 332, 340–345, 368, 370, 387–395, 400, 431–435, and 448–468 [<xref ref-type="bibr" rid="CR42">42</xref>]. Most common substitutions in antigenic sites were V106I (57.5%) followed by K432E (23.6%). The latter was detected in majority pH1N1 viruses in 2015. Additional substitutions in antigenic sites were N68 T (5%); H275Y (4%); N341D (3%) and E462K (6%) (Summarized in Additional file <xref rid="MOESM3" ref-type="media">3</xref>: Table S2). In addition to K432E substitutions mentioned earlier, five of emerged substitutions predominated in 2015: N200S, V241I, N248D, I314M, N386K and N369K. Of these, V241I, N369K, N386K and K432E, alone or in conjunction with oseltamivir resistance variant, H275Y, have been reported to significantly affect binding affinity of oseltamivir to neuraminidase, rendering neuraminidase less susceptible [<xref ref-type="bibr" rid="CR43">43</xref>, <xref ref-type="bibr" rid="CR44">44</xref>].</p><p id="Par35">Interestingly, the presence of I46T mutation in 11% of NA sequences in 2009 resulted in the acquisition of potential glycosylation at N44 residue. In total, eight glycosylation sites were found in N1: N50, N58, N63, N68, N88, N146, N235 and N386. Of these, N68 glycosylation was lost form 16% of pH1N1 viruses isolated from Lebanon during 2009. Similarly, N368 glycosylation sequon was missing from one sequence in 2011, however, during following years more NA sequences were found to lack 386 glycosylation (31% of sequences in 2012; 80% in 2014; 77% in 2015). The potential N-linked glycosylation sites in NA of pH1N1 viruses (2009–2017) are shown in Table <xref rid="Tab5" ref-type="table">5</xref>.</p></sec></sec><sec id="Sec18"><title>Genetic characterization of H3N2 viruses</title><sec id="Sec19"><title>Molecular evolution of H3 protein</title><p id="Par36">In order to track the molecular evolution of H3N2 viruses circulated in MENA region during 2009–2017, the potential genetic drift in HA and NA genes were investigated. Complete 239 HA protein sequences were analyzed and compared against vaccine strains of the respective years. The overall HA amino acid identities among the H3N2 strains isolated since 2009 compared to the respective vaccine strains were > 98%. The similarities of analyzed HA sequences were estimated to range from 99.1% in 2009 to 98.6% in 2017. Antigenic epitope A, was the most variable compared to vaccine strains reaching 77% identity to 2017 vaccine strain Hong Kong/2014, with A3 site being the least conserved compared to A1 and A2 sites (Table <xref rid="Tab6" ref-type="table">6</xref>). B antigenic epitope had 89.5% similarity to vaccine strains in all years except for viruses circulated in 2012 which had 100% similarity to Perth/2009 vaccine strain. Antigenic epitope C of 2016 viruses was only 66.6% conserved compared to vaccine strain (Switzerland/2013) while it had 100% similarities to vaccine strains in other years (2010–2015-2017). Similarities in D epitope ranged from 88% in 2012 with two substitutions, to 100% in subsequent years. The most conserved epitope of all was the E site, with no substitutions detected throughout the study period. Of interest, variants at two antigenic sites changed over time. Specifically, amino acid residue 144 of the antigenic epitope A3 which was asparagine in 2009 strains, was then mutated to Serine in all viruses isolated in 2015 and then to either Lysine (30%) and Serine (60%) of sequences isolated in 2017. A similar variation pattern was observed for amino acid 159 of B1 epitope. This amino acid residue was phenylalanine between 2009 and 2013, and then mutated to Tyrosine or Serine in 30% of HA sequences in 2014, and finally to Tyrosine in 2017 (Table <xref rid="Tab6" ref-type="table">6</xref>).<table-wrap id="Tab6"><label>Table 6</label><caption><p>Identification of amino acid substitutions in H3 antigenic sites of H3N2 viruses. Antigenic sites of H3 protein were compared to vaccine strains Brisbane/2007<bold>.</bold> Here, we are presenting changes in antigenic sites of viruses isolated following any modification in H3N2 vaccine strain (2012, 2014, 2015, 2016 and 2017). Dash (−) indicates amino acid similar to vaccine sequence while amino acids different than vaccine strain are indicated in <italic>italics</italic>. Bold amino acids denote amino acid changes found in vaccine strains recommended for subsequent season in comparison to Brisbane/2007. Similarity of antigenic sites between vaccine strain and circulating viruses is shown as percentages in the last column</p></caption><table frame="hsides" rules="groups"><thead><tr><th colspan="28">H3</th><th colspan="5">Similarity to vaccine strains (%)</th></tr><tr><th></th><th colspan="2">A1</th><th>A2</th><th colspan="6">A3</th><th colspan="8">B1</th><th colspan="2">B2</th><th colspan="2">C1</th><th colspan="2">C2</th><th>E1</th><th colspan="3">D</th><th rowspan="2" colspan="5"></th></tr><tr><th></th><th colspan="2">121</th><th>135</th><th colspan="2">142</th><th colspan="2">144</th><th colspan="2">145</th><th colspan="2">157</th><th colspan="2">158</th><th colspan="3">159</th><th>160</th><th>189</th><th>194</th><th colspan="2">53</th><th colspan="2">278</th><th>62</th><th colspan="2">212</th><th>215</th></tr></thead><tbody><tr><td>Brisbane/2007</td><td colspan="2">N</td><td>T</td><td colspan="2">R</td><td colspan="2">N</td><td colspan="2">N</td><td colspan="2">L</td><td colspan="2">K</td><td colspan="3">F</td><td>K</td><td>N</td><td>L</td><td colspan="2">D</td><td colspan="2">N</td><td>E</td><td colspan="2">T</td><td>P</td><td><bold>A</bold></td><td><bold>B</bold></td><td><bold>C</bold></td><td><bold>E</bold></td><td><bold>D</bold></td></tr><tr><td>2009</td><td colspan="2">–</td><td>–</td><td colspan="2">–</td><td colspan="2">–</td><td colspan="2">–</td><td colspan="2">–</td><td>K</td><td><italic>N</italic></td><td colspan="3">–</td><td>–</td><td><italic>K</italic></td><td>–</td><td colspan="2">–</td><td colspan="2">–</td><td>–</td><td>–</td><td><italic>A</italic></td><td>–</td><td>100</td><td>89.5</td><td>100</td><td>100</td><td>97</td></tr><tr><td>Perth/2009</td><td colspan="2">N</td><td>T</td><td colspan="2">R</td><td colspan="2"><bold>K</bold></td><td colspan="2">N</td><td colspan="2">L</td><td colspan="2"><bold>N</bold></td><td colspan="3">F</td><td>K</td><td><bold>K</bold></td><td>L</td><td colspan="2">D</td><td colspan="2">N</td><td><bold>K</bold></td><td colspan="2">T</td><td><bold>S</bold></td><td colspan="5"></td></tr><tr><td>2012</td><td colspan="2">–</td><td>–</td><td colspan="2">–</td><td>K</td><td><italic>N</italic></td><td>N</td><td><italic>S</italic></td><td colspan="2">–</td><td colspan="2">–</td><td colspan="3">–</td><td>–</td><td>–</td><td>–</td><td>D</td><td><italic>A</italic></td><td colspan="2"><italic>P</italic></td><td>–</td><td colspan="2"><italic>A</italic></td><td><italic>P</italic></td><td><italic>90</italic></td><td><italic>100</italic></td><td><italic>84</italic></td><td>100</td><td><italic>88</italic></td></tr><tr><td>Victoria/2011</td><td colspan="2">N</td><td>T</td><td colspan="2">R</td><td colspan="2">N</td><td colspan="2">N</td><td colspan="2">L</td><td colspan="2"><bold>N</bold></td><td colspan="3">F</td><td>K</td><td><bold>K</bold></td><td>L</td><td colspan="2">D</td><td colspan="2">N</td><td>E</td><td colspan="2"><bold>A</bold></td><td>P</td><td colspan="5"></td></tr><tr><td>2014</td><td colspan="2">–</td><td>–</td><td colspan="2"><italic>G</italic></td><td colspan="2">–</td><td colspan="2"><italic>S</italic></td><td>L</td><td><italic>S</italic></td><td colspan="2">–</td><td>F</td><td><italic>Y</italic></td><td><italic>S</italic></td><td>–</td><td>–</td><td>–</td><td colspan="2">–</td><td colspan="2">–</td><td>–</td><td colspan="2">–</td><td>–</td><td>90</td><td>89.5</td><td>84</td><td>100</td><td>100</td></tr><tr><td>Texas/2012</td><td colspan="2">N</td><td>T</td><td colspan="2">R</td><td colspan="2">N</td><td colspan="2">N</td><td colspan="2">L</td><td colspan="2"><bold>N</bold></td><td colspan="3">F</td><td>K</td><td><bold>K</bold></td><td>L</td><td colspan="2">D</td><td colspan="2"><bold>K</bold></td><td>E</td><td colspan="2"><bold>A</bold></td><td>P</td><td colspan="5"></td></tr><tr><td>2015</td><td colspan="2">–</td><td>–</td><td>R</td><td><italic>G</italic></td><td colspan="2"><italic>S</italic></td><td colspan="2"><italic>S</italic></td><td colspan="2">–</td><td colspan="2">–</td><td colspan="3"><italic>Y</italic></td><td><italic>T</italic></td><td>–</td><td>–</td><td colspan="2">–</td><td colspan="2">–</td><td>–</td><td colspan="2">–</td><td>–</td><td>86.3</td><td>89.5</td><td>100</td><td>100</td><td>100</td></tr><tr><td>Switzerland/2013</td><td colspan="2">N</td><td>T</td><td colspan="2"><bold>G</bold></td><td colspan="2">N</td><td colspan="2"><bold>S</bold></td><td colspan="2">L</td><td colspan="2"><bold>N</bold></td><td colspan="3"><bold>S</bold></td><td>K</td><td><bold>K</bold></td><td>L</td><td colspan="2">D</td><td colspan="2"><bold>K</bold></td><td>E</td><td colspan="2"><bold>A</bold></td><td>P</td><td colspan="5"></td></tr><tr><td>2016</td><td>N</td><td><italic>K</italic></td><td>–</td><td>G</td><td><italic>R</italic></td><td><italic>S</italic></td><td><italic>K</italic></td><td colspan="2">–</td><td colspan="2">–</td><td colspan="2">–</td><td><italic>F</italic></td><td colspan="2"><italic>Y</italic></td><td>–</td><td>–</td><td>–</td><td>–</td><td><italic>N</italic></td><td>K</td><td><italic>N</italic></td><td>–</td><td colspan="2">–</td><td>–</td><td>82</td><td>89.5</td><td>66.6</td><td>100</td><td>100</td></tr><tr><td>Hong Kong/2014</td><td colspan="2">N</td><td>T</td><td colspan="2">R</td><td colspan="2"><bold>S</bold></td><td colspan="2"><bold>S</bold></td><td colspan="2">L</td><td colspan="2"><bold>N</bold></td><td colspan="3"><bold>Y</bold></td><td>K</td><td><bold>K</bold></td><td><bold>P</bold></td><td colspan="2">D</td><td colspan="2"><bold>K</bold></td><td>E</td><td colspan="2"><bold>A</bold></td><td>P</td><td colspan="5"></td></tr><tr><td>2017</td><td colspan="2"><italic>K</italic></td><td><italic>K</italic></td><td>R</td><td><italic>G</italic></td><td>S</td><td><italic>K</italic></td><td colspan="2">–</td><td colspan="2">–</td><td colspan="2">–</td><td colspan="3">–</td><td><italic>T</italic></td><td>–</td><td><italic>L</italic></td><td colspan="2">–</td><td colspan="2">–</td><td>–</td><td colspan="2">–</td><td>–</td><td>77</td><td>89.5</td><td>100</td><td>100</td><td>100</td></tr></tbody></table></table-wrap></p><p id="Par37">A total of 32 amino acid substitutions were encountered in H3 proteins, of which 13 substitutions were located in HA head and 19 in the stem region (Table <xref rid="Tab7" ref-type="table">7</xref>). The antigenic sites had the largest proportion of substitutions found in HA head domain with four substitutions in each A and B epitopes, two in C epitope and one in E epitope. Three substitutions were observed within the HA RBS: one in 190-helix (N189K) and two in 220-loop (V223I and N225D). All of which became dominant variants in viruses circulated in 2016.Two more substitutions were also observed close to RBS: the sporadically detected R142G/K substitution (A site) and fixed A198S substitution. Overall, 16 of detected substitutions were fixed in the trunk of the evolutionary tree and were found in all H3N2 viruses in 2017 (Fig. <xref rid="Fig4" ref-type="fig">4</xref>b). Three additional substitutions were also found in > 50% of H3N2 viruses in 2017. Herein, the amino acids were numbered using H3 numbering scheme exclusive of the signal peptide.<table-wrap id="Tab7"><label>Table 7</label><caption><p>Amino acid substitutions identified in the HA of H3N2 viruses sequenced during 2009–2017. Amino acid substitutions were identified in comparison to A/Brisbane/10/2007. Bordered amino acids indicate substitutions in receptor binding site. Deleterious mutations were predicted by online server PROVEAN and are indicated in bold. Only one HA sequence was deposited from MENA countries during 2010 hence this year was excluded from mutation and glycosylation analysis</p></caption><table frame="hsides" rules="groups"><thead><tr><th>Reference strain</th><th></th><th colspan="2">AA substitution</th><th>Antigenic Site</th><th>2009 (%)</th><th>2011 (%)</th><th>2012 (%)</th><th>2013 (%)</th><th>2014 (%)</th><th>2015 (%)</th><th>2016 (%)</th><th>2017 (%)</th></tr></thead><tbody><tr><td>Brisbane/2007</td><td rowspan="31">HA1</td><td colspan="2">L3I</td><td></td><td></td><td></td><td></td><td>7</td><td>5</td><td>85</td><td>64</td><td>100</td></tr><tr><td rowspan="35"></td><td colspan="2">Q33R</td><td></td><td></td><td></td><td char="." align="char">9</td><td>19</td><td>80</td><td>93</td><td>68</td><td>100</td></tr><tr><td colspan="2">S45 N</td><td></td><td></td><td char="." align="char">23</td><td char="." align="char">5</td><td>69</td><td>90</td><td>96</td><td>83</td><td>100</td></tr><tr><td colspan="2">T48I</td><td></td><td></td><td char="." align="char">23</td><td char="." align="char">5</td><td>54</td><td>87.5</td><td>91</td><td>72</td><td>100</td></tr><tr><td colspan="2">D53N</td><td><bold>C1</bold></td><td></td><td char="." align="char">9</td><td char="." align="char">32</td><td>30.7</td><td></td><td>4</td><td>4</td><td>20</td></tr><tr><td colspan="2">E50R</td><td></td><td></td><td></td><td></td><td></td><td>10</td><td></td><td></td><td></td></tr><tr><td colspan="2">Q75H</td><td></td><td></td><td></td><td></td><td></td><td>10</td><td></td><td></td><td></td></tr><tr><td colspan="2">N81K</td><td><bold>E2</bold></td><td></td><td></td><td char="." align="char">8</td><td></td><td></td><td></td><td>6</td><td></td></tr><tr><td colspan="2">Y94H</td><td></td><td></td><td char="." align="char">9</td><td char="." align="char">38</td><td></td><td>12.5</td><td></td><td>4</td><td></td></tr><tr><td colspan="2">N121K</td><td><bold>A1</bold></td><td></td><td></td><td></td><td></td><td></td><td>4</td><td>25.5</td><td>80</td></tr><tr><td colspan="2">T128A</td><td></td><td></td><td></td><td></td><td>31</td><td>72.5</td><td>6</td><td>12.7</td><td></td></tr><tr><td rowspan="2">R142</td><td>G</td><td rowspan="2"><bold>A3</bold></td><td></td><td char="." align="char">9</td><td char="." align="char">2</td><td>30</td><td>72</td><td>10</td><td>12.7</td><td>30</td></tr><tr><td>K</td><td></td><td></td><td char="." align="char">13.5</td><td></td><td></td><td></td><td>12.7</td><td>10</td></tr><tr><td rowspan="2">N144</td><td>S</td><td rowspan="2"><bold>A3</bold></td><td></td><td></td><td></td><td></td><td>5</td><td>87</td><td>40</td><td>60</td></tr><tr><td>K</td><td></td><td></td><td></td><td></td><td></td><td></td><td>21</td><td>30</td></tr><tr><td colspan="2">N145S</td><td><bold>A3</bold></td><td></td><td char="." align="char">59</td><td char="." align="char">48.5</td><td>100</td><td>100</td><td>100</td><td>100</td><td>100</td></tr><tr><td colspan="2">L157S</td><td><bold>B1</bold></td><td char="." align="char">4</td><td char="." align="char">4</td><td></td><td>15</td><td>17.5</td><td>4</td><td>6</td><td></td></tr><tr><td rowspan="2">F159</td><td>Y</td><td rowspan="2"><bold>B1</bold></td><td></td><td></td><td></td><td></td><td>15</td><td>87</td><td>62</td><td>100</td></tr><tr><td>S</td><td></td><td></td><td></td><td></td><td>15</td><td></td><td>8</td><td></td></tr><tr><td colspan="2">K160 T</td><td><bold>B1</bold></td><td></td><td></td><td></td><td></td><td></td><td>83</td><td>53</td><td>90</td></tr><tr><td colspan="2">N171K</td><td></td><td></td><td></td><td></td><td></td><td></td><td>6</td><td>19</td><td>60</td></tr><tr><td colspan="2">N189K</td><td><bold>B2</bold></td><td char="." align="char">73</td><td char="." align="char">100</td><td char="." align="char">100</td><td>100</td><td>100</td><td>100</td><td>100</td><td>100</td></tr><tr><td colspan="2">A198S</td><td></td><td></td><td char="." align="char">77</td><td char="." align="char">43</td><td>61.5</td><td>100</td><td>91</td><td>83</td><td>100</td></tr><tr><td colspan="2">T212A</td><td><bold>D</bold></td><td char="." align="char">63.6</td><td char="." align="char">95</td><td char="." align="char">100</td><td>84</td><td>100</td><td>91</td><td>85</td><td>100</td></tr><tr><td colspan="2">V223I</td><td></td><td></td><td char="." align="char">91</td><td char="." align="char">46</td><td>84</td><td>100</td><td>100</td><td>85</td><td>100</td></tr><tr><td colspan="2">N225D</td><td></td><td></td><td></td><td></td><td></td><td>20</td><td>89</td><td>68</td><td>100</td></tr><tr><td colspan="2">R261Q</td><td></td><td></td><td></td><td></td><td>15</td><td></td><td></td><td>8</td><td>20</td></tr><tr><td colspan="2">N278K</td><td><bold>C2</bold></td><td></td><td char="." align="char">9</td><td></td><td>54</td><td>97</td><td>93</td><td>74</td><td>100</td></tr><tr><td colspan="2">E280A</td><td></td><td></td><td char="." align="char">4</td><td char="." align="char">38</td><td></td><td></td><td></td><td>4</td><td></td></tr><tr><td colspan="2">Q311H</td><td></td><td></td><td></td><td></td><td></td><td>7</td><td>85</td><td>59</td><td>60</td></tr><tr><td colspan="2">N312S</td><td></td><td char="." align="char">82</td><td char="." align="char">91</td><td char="." align="char">62</td><td>69</td><td>100</td><td>100</td><td>91</td><td>100</td></tr><tr><td rowspan="5">HA2</td><td rowspan="2">V347</td><td>M</td><td></td><td></td><td char="." align="char">9</td><td></td><td>23</td><td>15</td><td></td><td></td><td></td></tr><tr><td>K</td><td></td><td></td><td char="." align="char">22</td><td></td><td></td><td>40</td><td>4</td><td></td><td></td></tr><tr><td colspan="2"><bold>G484E</bold></td><td></td><td></td><td></td><td></td><td></td><td></td><td>2</td><td>27.5</td><td>50</td></tr><tr><td colspan="2">D487N</td><td></td><td></td><td char="." align="char">59</td><td char="." align="char">43</td><td>23</td><td></td><td></td><td></td><td></td></tr><tr><td colspan="2">D489N</td><td></td><td char="." align="char">9</td><td char="." align="char">13</td><td></td><td>15</td><td>12</td><td>87</td><td>64</td><td>100</td></tr><tr><td colspan="4">Total number of substitutions/years</td><td></td><td char="." align="char"><bold>5</bold></td><td char="." align="char"><bold>17</bold></td><td char="." align="char"><bold>16</bold></td><td><bold>19</bold></td><td><bold>23</bold></td><td><bold>25</bold></td><td><bold>28</bold></td><td><bold>23</bold></td></tr></tbody></table></table-wrap></p><p id="Par38">After 2009 pandemic, a few H3N2 viruses circulated carrying five amino acid substitutions compared to Brisbane/2007 vaccine strain. Most prevalent substitutions detected were N312S (82%) followed by N189K (73%) at B epitope and T212A (63.6%). Minority of viruses were also possessing L157S (4%) at B site and D489N (9%) substitutions. Unlike L157S, which was detected sporadically over years, D489N variant gradually replaced the parental variant in subsequent seasons before becoming the dominant variant in 2017(Fig. <xref rid="Fig4" ref-type="fig">4</xref>b). Only one HA sequence of H3N2 virus was deposited in IRD from MENA region during 2010 and was carrying N189K substitution. Despite the limited activity of H3N2 during 2011, various substitutions were introduced to H3N2 virus population. Besides previously reported N189K, T212A and N312S which became dominant variants (> 90%) in H3N2 viruses, 12 new substitutions were additionally identified during 2011. Among these, five substitutions were found at antigenic sites: D53N (9%) and N278K (9%) at C epitope; R142G (9%) and N145S (59%) at A epitope and L157S at B epitope. Further, six changes were introduced in 2011 and gradually predominated variants detected in later seasons: S45 N (23%), T48I (23%), A198S (77%) and V223I (91%) in HA1 subunit as well as D489N (13%) in HA2 subunit (Table <xref rid="Tab7" ref-type="table">7</xref>). In 2012, only two novel substitutions were detected; N81K (8%) at E epitope and Q33R (9%). In contrast to N81K that was reported sporadically over years, the frequency of Q33R increased gradually in subsequent years before becoming characteristic of all H3N2 viruses in 2017. All viruses had the N189K and T212A, while approximately 50% of them were carrying N145S, A198S, V223I and D487N substitutions by end of 2012. Other sporadic substitutions were also detected such as R142G/K (15.5%) and D53N (32%) found in antigenic epitopes A and C, respectively (Table <xref rid="Tab7" ref-type="table">7</xref>). The vaccine strain of the 2012–2013 season (Victoria/2011) was characterized by S45 N, T48I, N189K, A198S, T212A, V223I and N312S substitutions. Predominantly circulating viruses in 2013 carried these substitutions and additionally L3I (7%), T128A (31%) and R261Q (15%) substitutions. Unlike T128A and R261Q, L3I accumulated in H3N2 virus population and become a dominant variant in 2017. Of note, T128A mutation near A epitope of 17.6% of long circulated H3N2 viruses had been previously reported to be associated with high mortality rate among children [<xref ref-type="bibr" rid="CR45">45</xref>]. Moreover, previously reported N278K and N145S were found in 54 and 100% of viruses, respectively. In 2014, majority of H3N2 viruses (> 75%) were already possessing 11 of previously identified substitutions (Fig. <xref rid="Fig4" ref-type="fig">4</xref>b). Circulated viruses were also characterized by the appearance of four novel substitutions; N144S (5%) F159Y/S (30%), N225D (20%) and Q311H (7%). All of these variants became fixed in the backbone of H3 phylogenetic tree with F159Y and N225D being found in all viruses isolated in 2017. Also, 10% of viruses emerged with E50R and Q75H substitutions were exclusively found in 2014. During 2015, circulating viruses were showing 97.2% similarity to respective vaccine strain (Texas/2012) (Table <xref rid="Tab7" ref-type="table">7</xref>). Additional four novel substitutions were acquired: N121K (4%) at A epitope, K160 T (83%) at B epitope, N171K (6%) and G484E (2%). The frequencies of these substitutions increased dramatically in subsequent years reaching up to 80% for N121K and 90% for K160 T in 2017. In most cases, H3N2 viruses in 2015 possessed 17 fixed changes that have accumulated since 2009 (Fig. <xref rid="Fig4" ref-type="fig">4</xref>b). Towards 2016, viruses have become more stable carrying previously established substitutions and showing no novel mutation. In total, 28 amino acid substitutions were reported, however, 14 of them were already expressed in vaccine strain recommended for 2016 (Switzerland/2013). Likewise, no novel mutations were observed in HA protein during 2017. All H3N2 viruses were carrying 20 of fixed substitutions.</p><p id="Par39">Overall, H3N2 viruses acquired fixed mutations rapidly. Since 2009, the number of fixed mutations increased dramatically from one fixed substitution in 2009 to 15 in 2017 (Fig. <xref rid="Fig4" ref-type="fig">4</xref>b). By the end of 2015, the prevalence of fixed mutations was > 80%. Interestingly, the fixed mutation, G484E, which appeared first in 2015 (2%) was predicted to exert a deleterious effect on HA function (provean score − 4.94). Sporadic substitutions were also introduced several times in the virus population since 2009, e.g. D53N, Y94H, T128A, L157S, R261Q and E280A in HA1 subunit and V347 M/K and D487N in HA2 subunit (Table <xref rid="Tab7" ref-type="table">7</xref>). Other sporadic substitutions are listed in Table <xref rid="Tab7" ref-type="table">7</xref>.</p><p id="Par40">In addition to genetic drift, HA glycosylation pattern of H3N2 viruses has also changed since 2009. Nine putative N-glycosylation sites were found in all H3N2 isolates at 8, 22, 63, 126, 133, 165, 246, 285 and 483 (Fig. <xref rid="Fig5" ref-type="fig">5</xref>). However, the emergence of S45 N substitution in 2011 and K160 T substitution in 2015 was associated with acquisition of two additional glycosylation sequons at N45 (HA stem) and N158 (HA head) residues (Table <xref rid="Tab8" ref-type="table">8</xref>). Glycosylation sequon at N45 and N158 were found in 23 and 5% of viruses in 2011 and 2015 respectively. However, number of viruses expressing these two sequons has increased gradually until found in all H3N2 viruses in 2017 (Table <xref rid="Tab8" ref-type="table">8</xref>).<table-wrap id="Tab8"><label>Table 8</label><caption><p>Dynamics of HA and NA glycosylation for H3N2 viruses isolated from MENA region during 2009–2017. In total, glycosylation patterns of 239 HA and 132 NA sequences of H3N2 virus were analyzed. Numbers between brackets represent the percentage of sequences expressing glycosylation sequon in corresponding year while conserved glycosylation sites are indicated in bold</p></caption><table frame="hsides" rules="groups"><thead><tr><th colspan="13">Localization of potential glycosylation sequons in HA and NA of H3N2 viruses</th></tr></thead><tbody><tr><td colspan="13">HA</td></tr><tr><td></td><td></td><td colspan="11">Amino acid position</td></tr><tr><td>Year</td><td>Number of sequences analyzed</td><td colspan="5">Stem</td><td colspan="4">Head</td><td colspan="2">Stem</td></tr><tr><td>2009</td><td>22</td><td><bold>8</bold></td><td><bold>22</bold></td><td></td><td><bold>63</bold></td><td><bold>126</bold></td><td><bold>133</bold></td><td></td><td><bold>165</bold></td><td><bold>246</bold></td><td><bold>385</bold></td><td><bold>483</bold></td></tr><tr><td>2011</td><td>22</td><td><bold>8</bold></td><td><bold>22</bold></td><td>45(23%)</td><td><bold>63</bold></td><td><bold>126</bold></td><td><bold>133</bold></td><td></td><td><bold>165</bold></td><td><bold>246</bold></td><td><bold>385</bold></td><td><bold>483</bold></td></tr><tr><td>2012</td><td>37</td><td><bold>8</bold></td><td><bold>22</bold></td><td>45(5%)</td><td><bold>63</bold></td><td><bold>126</bold></td><td><bold>133</bold></td><td></td><td><bold>165</bold></td><td><bold>246</bold></td><td><bold>385</bold></td><td><bold>483</bold></td></tr><tr><td>2013</td><td>13</td><td><bold>8</bold></td><td><bold>22</bold></td><td>45(69%)</td><td><bold>63</bold></td><td><bold>126</bold></td><td><bold>133</bold></td><td></td><td><bold>165</bold></td><td><bold>246</bold></td><td><bold>385</bold></td><td><bold>483</bold></td></tr><tr><td>2014</td><td>40</td><td><bold>8</bold></td><td><bold>22</bold></td><td>45(92%)</td><td><bold>63</bold></td><td><bold>126</bold></td><td><bold>133</bold></td><td>158(5%)</td><td><bold>165</bold></td><td><bold>246</bold></td><td><bold>385</bold></td><td><bold>483</bold></td></tr><tr><td>2015</td><td>47</td><td><bold>8</bold></td><td><bold>22</bold></td><td>45(100%)</td><td><bold>63</bold></td><td><bold>126</bold></td><td><bold>133</bold></td><td>158(81%)</td><td><bold>165</bold></td><td><bold>246</bold></td><td><bold>385</bold></td><td><bold>483</bold></td></tr><tr><td>2016</td><td>47</td><td><bold>8</bold></td><td><bold>22</bold></td><td>45(80%)</td><td><bold>63</bold></td><td><bold>126</bold></td><td><bold>133</bold></td><td>158(53%)</td><td><bold>165</bold></td><td><bold>246</bold></td><td><bold>385</bold></td><td><bold>483</bold></td></tr><tr><td>2017</td><td>10</td><td><bold>8</bold></td><td><bold>22</bold></td><td>45(100%)</td><td><bold>63</bold></td><td><bold>126</bold></td><td><bold>133</bold></td><td>158(100%)</td><td><bold>165</bold></td><td><bold>246</bold></td><td><bold>385</bold></td><td><bold>483</bold></td></tr><tr><td colspan="13">NA</td></tr><tr><td>Year</td><td>Number of sequences analyzed</td><td colspan="11">Amino acid position</td></tr><tr><td>2009</td><td>5</td><td colspan="3"><bold>61</bold></td><td><bold>70</bold></td><td colspan="2"><bold>86</bold></td><td colspan="3"><bold>146</bold></td><td><bold>234</bold></td><td></td></tr><tr><td>2011</td><td>22</td><td colspan="3"><bold>61</bold></td><td><bold>70</bold></td><td colspan="2"><bold>86</bold></td><td colspan="3"><bold>146</bold></td><td><bold>234</bold></td><td>367 (90.6%)</td></tr><tr><td>2012</td><td>32</td><td colspan="3"><bold>61</bold></td><td><bold>70</bold></td><td colspan="2"><bold>86</bold></td><td colspan="3"><bold>146</bold></td><td><bold>234</bold></td><td>367 (100%)</td></tr><tr><td>2013</td><td>8</td><td colspan="3"><bold>61</bold></td><td><bold>70</bold></td><td colspan="2"><bold>86</bold></td><td colspan="3"><bold>146</bold></td><td><bold>234</bold></td><td>367 (100%)</td></tr><tr><td>2014</td><td>33</td><td colspan="3"><bold>61</bold></td><td><bold>70</bold></td><td colspan="2"><bold>86</bold></td><td colspan="3"><bold>146</bold></td><td><bold>234</bold></td><td>367 (100%)</td></tr><tr><td>2015</td><td>23</td><td colspan="3"><bold>61</bold></td><td><bold>70</bold></td><td colspan="2"><bold>86</bold></td><td colspan="3"><bold>146</bold></td><td><bold>234</bold></td><td>367 (100%)</td></tr><tr><td>2016</td><td>8</td><td colspan="3"><bold>61</bold></td><td><bold>70</bold></td><td colspan="2"><bold>86</bold></td><td colspan="3"><bold>146</bold></td><td><bold>234</bold></td><td>367 (100%)</td></tr></tbody></table></table-wrap></p></sec><sec id="Sec20"><title>Molecular evolution of N2 protein</title><p id="Par41">The NA gene was highly variable with total of 20 mutations detected in H3N2 strains included in this study compared to A/Brisbane/2007 virus (Additional file <xref rid="MOESM4" ref-type="media">4</xref>: Table S3). The catalytic and framework sites that support the catalytic residues within NA were highly conserved and are generally stable over time. However, eight substitutions were found near NA active site: Y155F, I176M, E221D, S245 N, S247 T, S367 N, K369 T and N402D. With the exception of Y155F and I176M, these substitutions were acquired by majority of H3N2 viruses (75%) in 2016. Importantly, none of the analyzed strains had the oseltamivir resistance sites H274Y, E119I or R292K [<xref ref-type="bibr" rid="CR26">26</xref>]. Six potential glycosylation sites were predicted in the NA gene of 132 H3N2 strains, including N61, N70, N86, N146, N234 and N367. The appearance of two substitutions S367 N and K369 T has led to acquisition of a new glycosylation site at 367 position was found in all N2 sequences isolated after 2010 (Table <xref rid="Tab7" ref-type="table">7</xref>).</p></sec><sec id="Sec21"><title>Comparison of mutation accumulation at head and stem of HA protein in pH1N1 and H3N2 subtypes</title><p id="Par42">When comparing the localization of substitutions in the HA protein of both subtypes, we found that most of substitutions in H3N2 were located in the stem (<italic>n</italic> = 19) whereas in pH1N1, we found the larger number of substitutions in the head region (<italic>n</italic> = 16) (Fig. <xref rid="Fig6" ref-type="fig">6</xref>a). Similarly, number of substitutions in antigenic site were much higher in H3N2 viruses (<italic>n</italic> = 12) compared to only 4 substitutions in pH1N1 viruses (Fig. <xref rid="Fig6" ref-type="fig">6</xref>a). Of the substitutions detected in antigenic sites during 2009, 3 were fixed in pH1N1 virus population compared to 7 fixed substitutions in H3N2 viruses in 2016 (Fig. <xref rid="Fig6" ref-type="fig">6</xref>b). Analysis of total number of fixed substitutions versus non-fixed substitutions in both subtypes revealed that the non-fixed substitutions predominated in pH1N1 (<italic>n</italic> = 18 out of 26) while fixed substitutions predominated in H3N2 (<italic>n</italic> = 20 out of 32). However, number of fixed substitutions in HA head and stem domains of both subtypes was equal. Interestingly, higher number of substitutions in HA head domain of pH1N1 virus was associated with higher number of pH1N1 positive cases during 2009–2011 and 2013–2016 (Fig. <xref rid="Fig7" ref-type="fig">7</xref>a). However, the higher number of substitutions in HA stem domain of pH1N1 was associated with a lower pH1N1 activity during 2012 and 2017. In contrast, the higher number of substitutions in HA stem was seen in H3N2 viruses during 2009–2016 was not associated with increased influenza activity (Fig. <xref rid="Fig7" ref-type="fig">7</xref>b).<fig id="Fig6"><label>Fig. 6</label><caption><p>Amino acid (AA) substitutions in HAs of pH1N1 and H3N2 viruses circulated between 2009 and 2017. <bold>a</bold> Number of amino acid substitutions in HA head and stem domains of pH1N1 and H3N2 viruses. Amino acids are colored according to their location in head (red), stem (blue) and antigenic sites (dark red). <bold>b</bold> Amino acid substitutions identified in the antigenic sites of HAs of pH1N1 and H3N2 viruses in 2009 and 2016. The 3D crystal structures of A/California/04/2009 (PDB 3LZG) and for 2005 human H3N2 Virus (PDB 2YP7) are displayed using CLC genomic workbench version 11. The antigenic sites of H1 are colored: Sa (red), Sb (green), Ca (blue) and Cb (orange). The antigenic sites of H3 are colored: A (red), B (blue), C (green), D (orange) and E (purple). The inner site of H3 structure (2016) is shown to visualize D53N substitution in C1 epitope. Fixed substitutions developed during the evolution of the influenza viruses are labeled in bold. AA substitutions in RBS are shown in <italic>italics</italic>. AA substitutions that are associated with changes of potential glycosylation sites are framed</p></caption><graphic xlink:href="12879_2019_3930_Fig6_HTML" id="MO7"></graphic></fig><fig id="Fig7"><label>Fig. 7</label><caption><p>Amino acid substitutions in the HA genes of (<bold>a</bold>) pH1N1 and (<bold>b</bold>) A(H3N2) viruses circulated during 2009–2017. Amino acid substitutions in the head domain of the HA protein are labeled in red and amino acid (AA) substitutions in the stem domain are displayed in blue. The black curve shows the number of positive samples as reported to FluNet (WHO) during the same period. The green arrows denote the introduction of a new vaccine strain</p></caption><graphic xlink:href="12879_2019_3930_Fig7_HTML" id="MO8"></graphic></fig></p></sec></sec></sec><sec id="Sec22"><title>Discussion</title><p id="Par43">We compared the relative contribution of influenza types (A and B) and subtypes (pH1N1 and H3N2) to influenza epidemics and investigated seasonality patterns in Arabian Gulf, Levant and North Africa regions during 2009–2017. Influenza surveillance data were retrieved from the global web-based tool ‘FluNet database’ which was launched first in 1997 by WHO for influenza surveillance [<xref ref-type="bibr" rid="CR24">24</xref>]. A preliminary check of data available revealed that reported influenza cases were unevenly distributed across countries, in addition to lack of data from three countries (Saudi Arabia, Libya and Yemen), thus a total number of 17 countries were involved in current study. A study comparing influenza seasonality patterns among MENA countries has been published earlier in 2018 by Caini and his colleagues [<xref ref-type="bibr" rid="CR46">46</xref>]. In the current study, we compared influenza seasonality based on WHO classification of influenza transmission zones and divided MENA countries in three main regions: Arabian Gulf (5 countries) and Levant (6 countries) which belong to Western Asia transmission zone and North Africa (5 countries) which is classified as North Africa transmission zone [<xref ref-type="bibr" rid="CR33">33</xref>].</p><p id="Par44">Influenza A was dominant in all years, with the pH1N1 accounting for the majority of cases except in 2012, during which H3N2 dominated. Influenza A and B accounted for 79.5 and 20.5% of all influenza cases during the study period, respectively. Generally, seasonality of influenza epidemics in all regions under study was in line with the northern hemisphere with a typical peak placed between December and March. Although Levant and Arabian Gulf lies in West Asia transmission zone (based on WHO classification), seasonality trends in Gulf were similar to that of North Africa rather than Levant. In all years, influenza activity started during October at high levels in Gulf and North Africa compared to a weak, late commence of the season in Levant, albeit the earlier winter seasons in the later. Furthermore, influenza positive cases were gradually increasing/decreasing in Gulf and North Africa, whereas a sharp increase/decrease of activity was reported in Levant (Additional file <xref rid="MOESM2" ref-type="media">2</xref>: Figure S1). The less defined seasonality pattern in Gulf countries could be partially explained by the tropical-like climate of these countries which is generally characterized by a wide primary peak in winter in addition to secondary smaller peaks seen throughout the year [<xref ref-type="bibr" rid="CR34">34</xref>]. In contrast, the location of Levant countries (Turkey-Lebanon-Syria-Jordan-Cyprus) in temperate regions was demonstrated by a clear seasonal variations of influenza epidemics, with a marked peak in cold winter months [<xref ref-type="bibr" rid="CR46">46</xref>–<xref ref-type="bibr" rid="CR48">48</xref>].</p><p id="Par45">In addition to monitoring activity of circulating influenza viruses, WHO utilizes surveillance data in studying virus evolution and assessing the antigenic characteristics of circulating viruses before selecting the most appropriate influenza vaccine strains for both hemispheres. Therefore, the molecular evolution of surface glycoproteins (HA and NA) of influenza A subtypes were also investigated in current study. All HA and NA sequences of pH1N1 and H3N2 viruses depicted in IRD since 2009 (total of 1260 sequences) were analyzed. Due to the limited number of sequences available, sequence analysis was done per MENA as whole rather than independent three regions. Phylogenetic analysis of both genes revealed that influenza viruses clustering in the phylogenetic tree is mainly chronological for both subtypes. Additionally, the phylogenetic analysis of the HA sequences revealed similar evolution of mutations circulating worldwide and hence similar clustering with globally circulating strains. In accordance to published data [<xref ref-type="bibr" rid="CR25">25</xref>], majority of pH1N1 viruses after 2014 clustered within the genetic clade 6B and related sub-clade, 6B.1, while in H3N2 viruses’ sub-clades 3C.2 and related sub-clade 3C.2a were the predominating. In contrast to our findings, though, pH1N1 viruses belonging to group 6 related sub-clades (6A, 6B and 6C) were detected globally during 2013–2014 compared to an earlier appearance (2011–2012) in MENA countries. Moreover, group 6C pH1N1 viruses were not found after 2013 in MENA countries, however, were found in viruses isolated in other countries during 2015. Likewise, phylogenetic analysis of H3N2 viruses has revealed some differences compared to globally reported data [<xref ref-type="bibr" rid="CR25">25</xref>]. Although global reports indicate that the majority of H3N2 viruses circulated during 2011–2012 fell into sub-clade 3C, viruses were mainly clustering within subclades 3A and 3B in MENA countries during that period.</p><p id="Par46">For further evolution analysis of the pH1N1 and H3N2 HAs, we determined the positions of identified substitutions with respect to their localization in HA protein structure. Generally, more substitutions were detected in H3N2 viruses (<italic>n</italic> = 32) compared to pH1N1 viruses (<italic>n</italic> = 24) during study period. In accordance with our findings, analysis of antigenic dynamics of both subtypes’ HAs revealed that H3N2 viruses evolve faster than human H1 subtype [<xref ref-type="bibr" rid="CR49">49</xref>]. The average amino acid mutation rate of the HA protein for H3N2 virus was 3.6 per year compared to 2.45 per year for H1N1 [<xref ref-type="bibr" rid="CR50">50</xref>]. One reason could be the continuous circulation of H3N2 virus in human population since 1968, resulting in higher immune pressure as compared to the relatively new pH1N1 virus [<xref ref-type="bibr" rid="CR51">51</xref>]. Substitutions were most frequently found in the head domain in pH1N1 viruses, including eight changes in the antigenic sites and five in the RBS, while interestingly, changes were mainly localized in stem domain in H3N2 viruses. Yet, most of substitutions identified in head domain of H3N2 viruses were localized in or near major antigenic sites (<italic>n</italic> = 8). Of note, changes of the antigenic HA composition of H3N2 viruses were accompanied by a six-time exchange of the H3N2 vaccine strains which exert more pressure on H3 antigenic sites compared to only two changes in pH1N1 vaccine strains during study period (Additional file <xref rid="MOESM1" ref-type="media">1</xref>: Table S1). Overall, the H3N2 strains displayed more diversity and faster accumulation of mutations in antigenic epitopes compared to pH1N1 strains [<xref ref-type="bibr" rid="CR49">49</xref>]. HA sequences from the H3N2 strains possessed two amino acid changes in one epitope (B; L157S and N189K) during the 2009. In 2011, this number increased to seven changes in three epitopes: two in epitope A, three in epitope B, and two in epitope C (Fig. <xref rid="Fig4" ref-type="fig">4</xref>b). In following years, more antigenic drift occurred due to changes in H3 antigenic site, however, most of these changes were found in recommended vaccine strain of subsequent seasons. The continuous change in H3N2 vaccine strains could also explain H3N2 lower activity despite the larger number of accumulated substitutions (Fig. <xref rid="Fig7" ref-type="fig">7</xref>b). Reports measuring influenza vaccine effectiveness have shown that the overall vaccine effectiveness (VE) was highest during 2010–2011 season (60%), while lowest during the 2014–2015 season (19%) [<xref ref-type="bibr" rid="CR52">52</xref>]. This was associated with change in vaccine strain of H3N2 from A/Victoria/361/2011 to A/Texas/50/2012. Thus, vaccine composition was reformulated again for 2015–2016 season to contain A/Switzerland/2013 and was shown to be more effective (48%) [<xref ref-type="bibr" rid="CR52">52</xref>]. In 2015, H3N2 viruses were carrying four substitutions in antigenic sites that were not found in A/Texas/2012 strain including: N144S and N145S in A3 site, and Y159F and K160 T in B1 site. Moreover, 87% of pH1N1 viruses had acquired S179 N mutation in 2015, which was associated with acquisition of potential glycosylation site. Together, these changes in H3N2 and pH1N1 viruses isolated in 2015 might -at least- partially explain the lower vaccine effectiveness during earlier months of 2015.</p><p id="Par47">The pandemic strain, pH1N1, originated from swine reservoir and circulated in human population for much shorter time (~ 9 years) compared to H3N2 viruses that has been circulating since 1968. Therefore, pH1N1 viruses are still considered to be in adaptive phase in terms of recognition of new receptor as well as avoiding immune response of the new host. Evidence of such adaptation process was revealed while analyzing mutations acquired by different pH1N1 viruses since 2009. Positive selection test of HA gene in pH1N1 viruses revealed that two amino acids in Ca antigenic site (D239G/E/N and Q240R) showed evidence of adaptive evolution. Owing to its presence in HA antigenic and RBS, variants of 239 position, in particular, has been found to be a key determinant of virus capability to bind epithelial cells in upper and/or lower respiratory tract [<xref ref-type="bibr" rid="CR53">53</xref>]. D239G variant has demonstrated the ability to bind α2–6 receptors in addition to increased binding affinity for α2–3 receptors expressed by pneumocytes in lower respiratory tract, and hence resulting in enhanced illness [<xref ref-type="bibr" rid="CR37">37</xref>]. 225G and 225 N variants were only detected in specimens obtained from patients with severe clinical symptoms, whereas D225E variants were detected in both severe and mild cases with similar frequencies [<xref ref-type="bibr" rid="CR53">53</xref>, <xref ref-type="bibr" rid="CR54">54</xref>]. In this study, D239 variants occurred in 10.3% of pH1N1 viruses. In contrast to D239E (8%) which appeared during 2009–2010, D239G and D239N substitutions formed no genetic cluster and reported sporadically suggesting that there was no spreading of viruses carrying this variant. Of note, variants of D239 position were also reported in seasonal strains circulated before 2009 [<xref ref-type="bibr" rid="CR39">39</xref>]. Receptor binding site of HA protein represented by 130-loop (148–152), 190-helix (201–208) and 220-loop (235–242) (H1 numbering) is one of the major determinants of virus transmission ability [<xref ref-type="bibr" rid="CR55">55</xref>]. Therefore, amino acid residues in and around the receptor binding site were predicted to change during adaptation process to human receptors. Those amino acids include: S200, S202, D204, A214, I233 and E241 [<xref ref-type="bibr" rid="CR56">56</xref>]. Four of these sites were found to be mutated in pH1N1 viruses included in this study. Of these, S200P, S202 and A214T were detected sporadically while I233T was introduced successfully in more than 80% of viruses after 2014. In 2013, de Vries and his colleagues has reported that S200P and S202 T substitutions would increase the receptor-binding avidity whereas the presence of A214T substitution was linked to decrease binding avidity [<xref ref-type="bibr" rid="CR57">57</xref>]. Interestingly, sequences expressing A214T variant were always associated with S200P and S202 T, which could be considered as compensatory process. With respect to I233T mutation, the presence of non-polar Isoleucine in 233 position has been found to disrupt the positioning of residues in the RBS from making a stable network of interactions [<xref ref-type="bibr" rid="CR58">58</xref>]. The substitution the Isoleucine with the polar amino acid Lysine, has been shown to significantly increase HA affinity to human receptor [<xref ref-type="bibr" rid="CR58">58</xref>]. Similarly, the displacement of the non-polar Isoleucine residue at 233 position with a polar amino acid, Threonine, could generate a stable ionic interaction between polar T233 residue in RBS and the acidic E241 residue in 220 loop, consequently increasing binding affinity to human receptor in comparison to earlier pH1N1 isolates [<xref ref-type="bibr" rid="CR58">58</xref>]. In addition to mutations in RBS, limited number of mutations in HA stem domain were also acquired in an early stage of pH1N1 evolution and thereafter introduced successfully in virus population. Among these mutations are E391K and E516K, which are found in the vicinity of the fusion peptide and transmembrane domain respectively. Such mutations are usually considered to be under negative selection because they might alter HA stability or functions [<xref ref-type="bibr" rid="CR59">59</xref>]. However, for a newly emerged virus such as pH1N1, such mutations may be prone to positive selection because it has to adapt to the new host. It is noteworthy, both subtype viruses acquired fixed mutations generally rapidly, however, more fixed mutations were accumulated in H3N2 compared to more non-fixed substitution in pH1N1, another evidence that the latter is still in adaptation process.</p><p id="Par48">In addition to antigenic drift, glycans on the globular head of HA can shield antigenic sites from neutralizing antibodies, but at the same time can affect HA binding to host cell receptors [<xref ref-type="bibr" rid="CR38">38</xref>]. Compared to vaccine strains, two amino acid mutations involving S179 N (Sa epitope) and K160 T (B epitope) in globular head of H1 and H3, respectively, have resulted in acquisition of putative glycosylation sequons. Glycosylation at the antigenic sites of HA head could be associated with increased viral pathogenicity by hiding antigenic sites from immune recognition. Interestingly, S179 N glycosylation was also reported in seasonal H1N1 viruses circulated during 1940–1948 before being replaced by N177 glycosylation in subsequent seasons until 2009 [<xref ref-type="bibr" rid="CR60">60</xref>]. Additional glycosylation sequons in positions 42 and 367 have been also acquired by N1 and N2 respectively. Unexpectedly, mutation rate of N2 (3.09 × 10<sup>− 3</sup>) was found to be higher than that of H3 (1.47 × 10<sup>− 3</sup>), however, this could be related to unusual number of mutations in N2 sequences which were not reported elsewhere. That in addition to low number of N2 sequences which were mainly deposited from single country (Iran).</p><p id="Par49">Systematic monitoring for influenza viruses is also important for evaluating the efficacy of the influenza NA inhibitors. Currently, neuraminidase inhibitors (NAIs) such as oseltamivir and zanamivir are the main drugs for treatment of influenza infections [<xref ref-type="bibr" rid="CR61">61</xref>]. The effectiveness of oseltamivir in treating influenza infections was threatened by the predominance of oseltamivir resistance among seasonal H1N1 viruses 2009–2010, even in countries where oseltamivir had not been used [<xref ref-type="bibr" rid="CR62">62</xref>]. However, seasonal virus was subsequently displaced by the pH1N1 virus, which is largely sensitive to oseltamivir [<xref ref-type="bibr" rid="CR63">63</xref>]. It is noteworthy, oseltamivir-resistant strains were rarely reported during 2009–2010 [<xref ref-type="bibr" rid="CR64">64</xref>]. In following years, though, a notable increase in the proportion of oseltamivir-resistant pH1N1 viruses were reported worldwide amongst patients with or without NAI treatment [<xref ref-type="bibr" rid="CR65">65</xref>, <xref ref-type="bibr" rid="CR66">66</xref>]. In accordance with our findings, most of reported oseltamivir-resistant cases exhibited geographical and temporal clustering. Such cases were reported during 2011 in New England, Australia [<xref ref-type="bibr" rid="CR67">67</xref>] and during 2013 in Sapporo, Japan [<xref ref-type="bibr" rid="CR68">68</xref>]. Here, all oseltamivir-resistant strains were reported were clustering in Oman and Iran during 2013 and 2015, respectively, suggesting efficient transmission of oseltamivir-resistant viruses among humans [<xref ref-type="bibr" rid="CR69">69</xref>]. This raises the concern that prevalence of oseltamivir-resistant pH1N1 viruses may increase in the future as the case of the previously circulating seasonal H1N1 viruses. Interestingly, two of NA fixed mutations reported in current study, V241I and N369K, have been found to offer robust fitness on pH1N1 viruses carrying H275Y mutation by enhancing NA surface expression and its enzymatic activity [<xref ref-type="bibr" rid="CR44">44</xref>]. The presence of such permissive mutations in all pH1N1 viruses suggesting that they are more permissive to the acquisition of H275Y and hence increase the risk of oseltamivir-resistant strains to spread globally. Here, we found oseltamivir-resistance mutation, H275Y, in only 4% of pH1N1 virus while no resistance mutations were reported in H3N2 viruses besides both subtypes have retained susceptibility to zanamivir. Unfortunately, we were unable to find published resources or reports regarding the prevalence of NAIs use in MENA countries.</p><p id="Par50">There are limitations in this study. Surveillance data was not available from Saudi Arabia, Yemen or Libya. The lack of influenza surveillance data from Saudi Arabia especially during Umra and Haj seasons could underestimate the actual prevalence of influenza infection in MENA region. Nevertheless, several studies investigating the burden of respiratory infection during Hajj season were published [<xref ref-type="bibr" rid="CR70">70</xref>]. Moreover, there is a relatively small number of HA and NA genes depicted in Influenza Research Database form MENA countries, especially for H3N2 viruses. This mandates further molecular characterization of circulating viruses in the region. Analysis of small number of sequences could preclude exploring the actual genetic diversity of influenza viruses.</p></sec><sec id="Sec23"><title>Conclusion</title><p id="Par51">In the context of the national influenza surveillance, we estimated the epidemiology of influenza types A and B in Arabian Gulf, Levant and North Africa regions since 2009. The pandemic strain pH1N1 was the dominant subtype in all years except in 2012 which was dominated by H3N2. Although WHO includes both Gulf and Levant in Western Asia, we found that seasonality patterns of Gulf and North Africa regions are similar compared to Levant region. We have also characterized the genetic evolution of influenza A subtypes: pH1N1 and H3N2 viruses during same period. Comparison of identified mutations of pH1N1 and H3N2 viruses encountered over 9 years revealed significant differences with regard to position and function of identified substitutions. pH1N1 viruses accumulated more substitutions in the head domain whereas H3N2 viruses had most of the reported substitutions in the stem domain. Moreover, H3N2 viruses acquired more fixed substitutions that were successfully introduced in viral population compared to larger number of non-fixed substitutions in pH1N1. More importantly, amino acids in the antigenic sites of H3N2 viruses were showing more variability than those of pH1N1 viruses, probably indicating greater immunological pressure. Interestingly, the increase of substitutions in HA globular head of pH1N1 virus was seen in strong influenza activity while that was not the case in H3N2 viruses despite the higher number of accumulated mutations. A possible explanation could be that the H3N2 viruses have been circulating in human population since 1968 resulting in circulation of different variants of virus compared to the relatively novel pH1N1. On the other hand, molecular analysis of the neuraminidase genes revealed the clustering of Y275H oseltamivir resistant mutation in pH1N1 viruses isolated from Oman and Iran, however, neither NAs carried any of substitutions associated with reduced susceptibility to zanamivir. Major challenges are yet to be faced by MENA countries in order to contribute to global surveillance of influenza virus by improving their surveillance systems for detection and response to all public health events.</p></sec><sec sec-type="supplementary-material"><title>Additional files</title><sec id="Sec24"><p><supplementary-material content-type="local-data" id="MOESM1"><media xlink:href="12879_2019_3930_MOESM1_ESM.docx"><label>Additional file 1:</label><caption><p><bold>Table S1.</bold> Accession numbers of pH1N1 and H3N2 vaccine and representative strains as obtained from NCBI and GISAID. Vaccine strains were used as references when analyzing HA and NA sequences. Representative stains were used to identify the clades and sub-clades in phylogenetic tree. (DOCX 31 kb)</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="MOESM2"><media xlink:href="12879_2019_3930_MOESM2_ESM.docx"><label>Additional file 2:</label><caption><p><bold>Figure S1.</bold> Seasonality patterns of influenza viruses in MENA during 2009–2017. Here, we are presenting influenza activity from September to April of each year in Arabian Gulf, Levant and North Africa regions. Monthly trend of each subtype is presented by bars: red bars represents percentage of pH1N1 positive cases, Blue represents H3N2 positive cases and green represents Flu B positive cases. (DOCX 159 kb)</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="MOESM3"><media xlink:href="12879_2019_3930_MOESM3_ESM.docx"><label>Additional file 3:</label><caption><p><bold>Table S2.</bold> Accumulation of amino acid substitutions in N1 protein of pH1N1 viruses of during 2009–2017: Amino acid substitutions were identified relative to A/California/07/2009 vaccine strain. The last column shows the overall prevalence of each substitution throughout study period (2009–2017). Substitutions associated with NAIs activity are indicated in bold. N1 numbering was used for reporting substitutions. (DOCX 34 kb)</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="MOESM4"><media xlink:href="12879_2019_3930_MOESM4_ESM.docx"><label>Additional file 4:</label><caption><p><bold>Table S3.</bold> Accumulation of amino acid substitutions in NA protein of in N2 protein of H3N2 viruses of during 2009–2017: Amino acid substitutions were identified relative to A/Brisbane/2007 vaccine strain. The last column shows the overall prevalence of each substitution in all NA sequences included in study (2009–2017). N2 numbering was used for reporting substitutions. (DOCX 34 kb)</p></caption></media></supplementary-material></p></sec></sec></body><back><glossary><title>Abbreviations</title><def-list><def-item><term>HA</term><def><p id="Par5">Hemagglutinin</p></def></def-item><def-item><term>IRD</term><def><p id="Par6">Influenza Research Database</p></def></def-item><def-item><term>MENA</term><def><p id="Par7">Middle East and North Africa</p></def></def-item><def-item><term>NA</term><def><p id="Par8">Neuraminidase</p></def></def-item><def-item><term>NAIs</term><def><p id="Par9">Neuraminidase inhibitors</p></def></def-item><def-item><term>RBS</term><def><p id="Par10">Receptor Binding Site</p></def></def-item></def-list></glossary><ack><p>We highly appreciate all submitting and originating laboratories of the surveillance data (FluNet-WHO) and sequences (GISAID and IRD) we used in this study.</p><sec id="FPar1"><title>Funding</title><p id="Par52">This study was supported by funds from Qatar University, grant number QUCG-BRC-2018\2019–1. The publication of this article was funded by the Qatar National Library.</p></sec><sec id="FPar2" sec-type="data-availability"><title>Availability of data and materials</title><p id="Par53">Virological surveillance data is available on WHO FluNet database (<ext-link ext-link-type="uri" xlink:href="http://www.who.int/flunet">http://www.who.int/flunet</ext-link>) while HA and NA sequences of both subtypes is available on Influenza Research Database (IRD) (<ext-link ext-link-type="uri" xlink:href="http://www.fludb.org">http://www.fludb.org</ext-link>).</p></sec></ack><notes notes-type="author-contribution"><title>Authors’ contributions</title><p>HMY developed the theory, HAA, AAA and HMY designed the experiments, HAA performed the computations, HMY and IG validated the analysis, HAA wrote the first draft of the paper, all authors read and approved the manuscript.</p></notes><notes notes-type="COI-statement"><sec id="FPar3"><title>Ethics approval and consent to participate</title><p>Not applicable.</p></sec><sec id="FPar4"><title>Consent for publication</title><p>Not applicable.</p></sec><sec id="FPar5"><title>Competing interests</title><p>Hadi Yassine is a member of the editorial board of the journal. Otherwise, the authors declare they have no conflict of interest.</p></sec><sec id="FPar6"><title>Publisher’s Note</title><p>Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.</p></sec></notes><ref-list id="Bib1"><title>References</title><ref id="CR1"><label>1.</label><mixed-citation publication-type="other">World Health Organization- Influenza fact sheet. Influenza (Seasonal) Fact sheet. 2017 [cited 2018 February]; Available from: <ext-link ext-link-type="uri" xlink:href="http://www.who.int/news-room/detail/14-12-2017-up-to-650-000-people-die-of-respiratory-diseases-linked-to-seasonal-flu-each-year">http://www.who.int/news-room/detail/14-12-2017-up-to-650-000-people-die-of-respiratory-diseases-linked-to-seasonal-flu-each-year</ext-link></mixed-citation></ref><ref id="CR2"><label>2.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hay</surname><given-names>AJ</given-names></name><etal></etal></person-group><article-title>The evolution of human influenza viruses</article-title><source>Philos Trans R Soc Lond Ser B Biol Sci</source><year>2001</year><volume>356</volume><issue>1416</issue><fpage>1861</fpage><lpage>1870</lpage><pub-id pub-id-type="doi">10.1098/rstb.2001.0999</pub-id><pub-id pub-id-type="pmid">11779385</pub-id></element-citation></ref><ref id="CR3"><label>3.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Molinari</surname><given-names>NA</given-names></name><etal></etal></person-group><article-title>The annual impact of seasonal influenza in the US: measuring disease burden and costs</article-title><source>Vaccine</source><year>2007</year><volume>25</volume><issue>27</issue><fpage>5086</fpage><lpage>5096</lpage><pub-id pub-id-type="doi">10.1016/j.vaccine.2007.03.046</pub-id><pub-id pub-id-type="pmid">17544181</pub-id></element-citation></ref><ref id="CR4"><label>4.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Schafer</surname><given-names>JR</given-names></name><etal></etal></person-group><article-title>Origin of the pandemic 1957 H2 influenza a virus and the persistence of its possible progenitors in the avian reservoir</article-title><source>Virology</source><year>1993</year><volume>194</volume><issue>2</issue><fpage>781</fpage><lpage>788</lpage><pub-id pub-id-type="doi">10.1006/viro.1993.1319</pub-id><pub-id pub-id-type="pmid">7684877</pub-id></element-citation></ref><ref id="CR5"><label>5.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Webster</surname><given-names>RG</given-names></name><etal></etal></person-group><article-title>Evolution and ecology of influenza a viruses</article-title><source>Microbiol Rev</source><year>1992</year><volume>56</volume><issue>1</issue><fpage>152</fpage><lpage>179</lpage><pub-id pub-id-type="pmid">1579108</pub-id></element-citation></ref><ref id="CR6"><label>6.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cohen</surname><given-names>J.</given-names></name><name><surname>Enserink</surname><given-names>M.</given-names></name></person-group><article-title>After Delays, WHO Agrees: The 2009 Pandemic Has Begun</article-title><source>Science</source><year>2009</year><volume>324</volume><issue>5934</issue><fpage>1496</fpage><lpage>1497</lpage><pub-id pub-id-type="doi">10.1126/science.324_1496</pub-id><pub-id pub-id-type="pmid">19541958</pub-id></element-citation></ref><ref id="CR7"><label>7.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cohen</surname><given-names>AL</given-names></name><etal></etal></person-group><article-title>Epidemiology of influenza virus types and subtypes in South Africa, 2009-2012</article-title><source>Emerg Infect Dis</source><year>2014</year><volume>20</volume><issue>7</issue><fpage>1162</fpage><lpage>1169</lpage><pub-id pub-id-type="doi">10.3201/eid2007.131869</pub-id><pub-id pub-id-type="pmid">24960314</pub-id></element-citation></ref><ref id="CR8"><label>8.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Martinez</surname><given-names>O</given-names></name><name><surname>Tsibane</surname><given-names>T</given-names></name><name><surname>Basler</surname><given-names>CF</given-names></name></person-group><article-title>Neutralizing anti-influenza virus monoclonal antibodies: therapeutics and tools for discovery</article-title><source>Int Rev Immunol</source><year>2009</year><volume>28</volume><issue>1</issue><fpage>69</fpage><lpage>92</lpage><pub-id pub-id-type="doi">10.1080/08830180802593540</pub-id><pub-id pub-id-type="pmid">19241254</pub-id></element-citation></ref><ref id="CR9"><label>9.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gerhard</surname><given-names>W</given-names></name><etal></etal></person-group><article-title>Antigenic structure of influenza virus haemagglutinin defined by hybridoma antibodies</article-title><source>Nature</source><year>1981</year><volume>290</volume><issue>5808</issue><fpage>713</fpage><lpage>717</lpage><pub-id pub-id-type="doi">10.1038/290713a0</pub-id><pub-id pub-id-type="pmid">6163993</pub-id></element-citation></ref><ref id="CR10"><label>10.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wiley</surname><given-names>DC</given-names></name><name><surname>Skehel</surname><given-names>JJ</given-names></name></person-group><article-title>The structure and function of the hemagglutinin membrane glycoprotein of influenza virus</article-title><source>Annu Rev Biochem</source><year>1987</year><volume>56</volume><fpage>365</fpage><lpage>394</lpage><pub-id pub-id-type="doi">10.1146/annurev.bi.56.070187.002053</pub-id><pub-id pub-id-type="pmid">3304138</pub-id></element-citation></ref><ref id="CR11"><label>11.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Igarashi</surname><given-names>M</given-names></name><etal></etal></person-group><article-title>Predicting the antigenic structure of the pandemic (H1N1) 2009 influenza virus hemagglutinin</article-title><source>PLoS One</source><year>2010</year><volume>5</volume><issue>1</issue><fpage>e8553</fpage><pub-id pub-id-type="doi">10.1371/journal.pone.0008553</pub-id><pub-id pub-id-type="pmid">20049332</pub-id></element-citation></ref><ref id="CR12"><label>12.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gallagher</surname><given-names>P</given-names></name><etal></etal></person-group><article-title>Addition of carbohydrate side chains at novel sites on influenza virus hemagglutinin can modulate the folding, transport, and activity of the molecule</article-title><source>J Cell Biol</source><year>1988</year><volume>107</volume><issue>6 Pt 1</issue><fpage>2059</fpage><lpage>2073</lpage><pub-id pub-id-type="doi">10.1083/jcb.107.6.2059</pub-id><pub-id pub-id-type="pmid">2461945</pub-id></element-citation></ref><ref id="CR13"><label>13.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sun</surname><given-names>S</given-names></name><etal></etal></person-group><article-title>Prediction of biological functions on glycosylation site migrations in human influenza H1N1 viruses</article-title><source>PLoS One</source><year>2012</year><volume>7</volume><issue>2</issue><fpage>e32119</fpage><pub-id pub-id-type="doi">10.1371/journal.pone.0032119</pub-id><pub-id pub-id-type="pmid">22355413</pub-id></element-citation></ref><ref id="CR14"><label>14.</label><mixed-citation publication-type="other">World Health Organization-Vaccines. WHO recommendations on the composition of influenza virus vaccines. 2018 [cited 2018 January]; Available from: <ext-link ext-link-type="uri" xlink:href="http://www.who.int/influenza/vaccines/virus/recommendations/en/">http://www.who.int/influenza/vaccines/virus/recommendations/en/</ext-link>.</mixed-citation></ref><ref id="CR15"><label>15.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname><given-names>Z</given-names></name><name><surname>Zhou</surname><given-names>H</given-names></name><name><surname>Jin</surname><given-names>H</given-names></name></person-group><article-title>The impact of key amino acid substitutions in the hemagglutinin of influenza a (H3N2) viruses on vaccine production and antibody response</article-title><source>Vaccine</source><year>2010</year><volume>28</volume><issue>24</issue><fpage>4079</fpage><lpage>4085</lpage><pub-id pub-id-type="doi">10.1016/j.vaccine.2010.03.078</pub-id><pub-id pub-id-type="pmid">20399830</pub-id></element-citation></ref><ref id="CR16"><label>16.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Koel</surname><given-names>BF</given-names></name><etal></etal></person-group><article-title>Identification of amino acid substitutions supporting antigenic change of influenza a(H1N1)pdm09 viruses</article-title><source>J Virol</source><year>2015</year><volume>89</volume><issue>7</issue><fpage>3763</fpage><lpage>3775</lpage><pub-id pub-id-type="doi">10.1128/JVI.02962-14</pub-id><pub-id pub-id-type="pmid">25609810</pub-id></element-citation></ref><ref id="CR17"><label>17.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Haq</surname><given-names>Z</given-names></name><name><surname>Mahjour</surname><given-names>J</given-names></name><name><surname>Khan</surname><given-names>W</given-names></name></person-group><article-title>Communicable diseases in the eastern Mediterranean region: prevention and control 2010-2011</article-title><source>East Mediterr Health J</source><year>2013</year><volume>19</volume><issue>10</issue><fpage>888</fpage><lpage>891</lpage><pub-id pub-id-type="doi">10.26719/2013.19.10.888</pub-id><pub-id pub-id-type="pmid">24313154</pub-id></element-citation></ref><ref id="CR18"><label>18.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Baghi</surname><given-names>Hossein Bannazadeh</given-names></name><name><surname>Soroush</surname><given-names>Mohammad Hossein</given-names></name></person-group><article-title>Influenza: the role of the Middle East and north Africa?</article-title><source>The Lancet Infectious Diseases</source><year>2018</year><volume>18</volume><issue>5</issue><fpage>498</fpage><pub-id pub-id-type="doi">10.1016/S1473-3099(18)30213-5</pub-id><pub-id pub-id-type="pmid">29695361</pub-id></element-citation></ref><ref id="CR19"><label>19.</label><mixed-citation publication-type="other">Gulf labour markets, m.a.p.p. Total population and percentage of nationals and foreign nationals in GCC countries 2016 [cited 2018 May 2018]; Available from: <ext-link ext-link-type="uri" xlink:href="http://gulfmigration.eu/gcc-total-population-percentage-nationals-foreign-nationals-gcc-countries-national-statistics-2010-2016-numbers/">http://gulfmigration.eu/gcc-total-population-percentage-nationals-foreign-nationals-gcc-countries-national-statistics-2010-2016-numbers/</ext-link>.</mixed-citation></ref><ref id="CR20"><label>20.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Buliva</surname><given-names>E</given-names></name><etal></etal></person-group><article-title>Emerging and reemerging diseases in the World Health Organization (WHO) eastern Mediterranean region-Progress, challenges, and WHO initiatives</article-title><source>Front Public Health</source><year>2017</year><volume>5</volume><fpage>276</fpage><pub-id pub-id-type="doi">10.3389/fpubh.2017.00276</pub-id><pub-id pub-id-type="pmid">29098145</pub-id></element-citation></ref><ref id="CR21"><label>21.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Soebiyanto</surname><given-names>RP</given-names></name><etal></etal></person-group><article-title>Associations between meteorological parameters and influenza activity in Berlin (Germany), Ljubljana (Slovenia), castile and Leon (Spain) and Israeli districts</article-title><source>PLoS One</source><year>2015</year><volume>10</volume><issue>8</issue><fpage>e0134701</fpage><pub-id pub-id-type="doi">10.1371/journal.pone.0134701</pub-id><pub-id pub-id-type="pmid">26309214</pub-id></element-citation></ref><ref id="CR22"><label>22.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Khan</surname><given-names>W</given-names></name><etal></etal></person-group><article-title>Influenza Research in the eastern Mediterranean region: a review</article-title><source>Oman Med J</source><year>2017</year><volume>32</volume><issue>5</issue><fpage>359</fpage><lpage>364</lpage><pub-id pub-id-type="doi">10.5001/omj.2017.70</pub-id><pub-id pub-id-type="pmid">29026466</pub-id></element-citation></ref><ref id="CR23"><label>23.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tamerius</surname><given-names>JD</given-names></name><etal></etal></person-group><article-title>Environmental predictors of seasonal influenza epidemics across temperate and tropical climates</article-title><source>PLoS Pathog</source><year>2013</year><volume>9</volume><issue>3</issue><fpage>e1003194</fpage><pub-id pub-id-type="doi">10.1371/journal.ppat.1003194</pub-id><pub-id pub-id-type="pmid">23505366</pub-id></element-citation></ref><ref id="CR24"><label>24.</label><mixed-citation publication-type="other">World Health Organization-FluNet. FluNet database. 2018 [cited 2018 January]; Available from: <ext-link ext-link-type="uri" xlink:href="http://www.who.int/influenza/gisrs_laboratory/flunet/en/">http://www.who.int/influenza/gisrs_laboratory/flunet/en/</ext-link>.</mixed-citation></ref><ref id="CR25"><label>25.</label><mixed-citation publication-type="other">World Health Organization- London Collaberationg center. Influenza Collaberating Center-London-Annual and Interim Reports. 2018 [cited 2018 February]; Available from: <ext-link ext-link-type="uri" xlink:href="http://www.nimr.mrc.ac.uk/who-influenza-centre/annual-and-interim-reports/">http://www.nimr.mrc.ac.uk/who-influenza-centre/annual-and-interim-reports/</ext-link>.</mixed-citation></ref><ref id="CR26"><label>26.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nguyen</surname><given-names>HT</given-names></name><name><surname>Fry</surname><given-names>AM</given-names></name><name><surname>Gubareva</surname><given-names>LV</given-names></name></person-group><article-title>Neuraminidase inhibitor resistance in influenza viruses and laboratory testing methods</article-title><source>Antivir Ther</source><year>2012</year><volume>17</volume><issue>1 Pt B</issue><fpage>159</fpage><lpage>173</lpage><pub-id pub-id-type="doi">10.3851/IMP2067</pub-id><pub-id pub-id-type="pmid">22311680</pub-id></element-citation></ref><ref id="CR27"><label>27.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tamura</surname><given-names>K</given-names></name><etal></etal></person-group><article-title>MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods</article-title><source>Mol Biol Evol</source><year>2011</year><volume>28</volume><issue>10</issue><fpage>2731</fpage><lpage>2739</lpage><pub-id pub-id-type="doi">10.1093/molbev/msr121</pub-id><pub-id pub-id-type="pmid">21546353</pub-id></element-citation></ref><ref id="CR28"><label>28.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yang</surname><given-names>H</given-names></name><name><surname>Carney</surname><given-names>P</given-names></name><name><surname>Stevens</surname><given-names>J</given-names></name></person-group><article-title>Structure and receptor binding properties of a pandemic H1N1 virus hemagglutinin</article-title><source>PLoS Curr</source><year>2010</year><volume>2</volume><fpage>Rrn1152</fpage><pub-id pub-id-type="doi">10.1371/currents.RRN1152</pub-id><pub-id pub-id-type="pmid">20352039</pub-id></element-citation></ref><ref id="CR29"><label>29.</label><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Huang</surname><given-names>J-W</given-names></name><name><surname>Yang</surname><given-names>J-M</given-names></name></person-group><source>Changed epitopes drive the antigenic drift for influenza A (H3N2) viruses</source><year>2011</year></element-citation></ref><ref id="CR30"><label>30.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Choi</surname><given-names>Y</given-names></name><etal></etal></person-group><article-title>Predicting the functional effect of amino acid substitutions and indels</article-title><source>PLoS One</source><year>2012</year><volume>7</volume><issue>10</issue><fpage>e46688</fpage><pub-id pub-id-type="doi">10.1371/journal.pone.0046688</pub-id><pub-id pub-id-type="pmid">23056405</pub-id></element-citation></ref><ref id="CR31"><label>31.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Drummond</surname><given-names>AJ</given-names></name><name><surname>Rambaut</surname><given-names>A</given-names></name></person-group><article-title>BEAST: Bayesian evolutionary analysis by sampling trees</article-title><source>BMC Evol Biol</source><year>2007</year><volume>7</volume><fpage>214</fpage><pub-id pub-id-type="doi">10.1186/1471-2148-7-214</pub-id><pub-id pub-id-type="pmid">17996036</pub-id></element-citation></ref><ref id="CR32"><label>32.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kosakovsky Pond</surname><given-names>SL</given-names></name><name><surname>Frost</surname><given-names>SD</given-names></name></person-group><article-title>Not so different after all: a comparison of methods for detecting amino acid sites under selection</article-title><source>Mol Biol Evol</source><year>2005</year><volume>22</volume><issue>5</issue><fpage>1208</fpage><lpage>1222</lpage><pub-id pub-id-type="doi">10.1093/molbev/msi105</pub-id><pub-id pub-id-type="pmid">15703242</pub-id></element-citation></ref><ref id="CR33"><label>33.</label><mixed-citation publication-type="other">World Health Organization- Influenza transmission zones. Influenza transmission zones. 2018 [cited 2017 December]; Available from: <ext-link ext-link-type="uri" xlink:href="http://www.who.int/influenza/surveillance_monitoring/updates/EN_GIP_Influenza_transmission_zones">http://www.who.int/influenza/surveillance_monitoring/updates/EN_GIP_Influenza_transmission_zones</ext-link>.</mixed-citation></ref><ref id="CR34"><label>34.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Viboud</surname><given-names>C</given-names></name><name><surname>Alonso</surname><given-names>WJ</given-names></name><name><surname>Simonsen</surname><given-names>L</given-names></name></person-group><article-title>Influenza in tropical regions</article-title><source>PLoS Med</source><year>2006</year><volume>3</volume><issue>4</issue><fpage>e89</fpage><pub-id pub-id-type="doi">10.1371/journal.pmed.0030089</pub-id><pub-id pub-id-type="pmid">16509764</pub-id></element-citation></ref><ref id="CR35"><label>35.</label><mixed-citation publication-type="other">Tanarhte MHP, Lelieveld J. Intercomparison of temperature and precipitation data sets based on observations in the Mediterranean and the Middle East. J Geophysical Res. 2012:117.</mixed-citation></ref><ref id="CR36"><label>36.</label><mixed-citation publication-type="other">World Health Organisation, Influenza Updat. 2009-2016, WHO.</mixed-citation></ref><ref id="CR37"><label>37.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chutinimitkul</surname><given-names>S</given-names></name><etal></etal></person-group><article-title>Virulence-associated substitution D222G in the hemagglutinin of 2009 pandemic influenza a(H1N1) virus affects receptor binding</article-title><source>J Virol</source><year>2010</year><volume>84</volume><issue>22</issue><fpage>11802</fpage><lpage>11813</lpage><pub-id pub-id-type="doi">10.1128/JVI.01136-10</pub-id><pub-id pub-id-type="pmid">20844044</pub-id></element-citation></ref><ref id="CR38"><label>38.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Schulze</surname><given-names>IT</given-names></name></person-group><article-title>Effects of glycosylation on the properties and functions of influenza virus hemagglutinin</article-title><source>J Infect Dis</source><year>1997</year><volume>176</volume><issue>Suppl 1</issue><fpage>S24</fpage><lpage>S28</lpage><pub-id pub-id-type="doi">10.1086/514170</pub-id><pub-id pub-id-type="pmid">9240690</pub-id></element-citation></ref><ref id="CR39"><label>39.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sriwilaijaroen</surname><given-names>N</given-names></name><name><surname>Suzuki</surname><given-names>Y</given-names></name></person-group><article-title>Molecular basis of the structure and function of H1 hemagglutinin of influenza virus</article-title><source>Proc Jpn Acad Ser B Phys Biol Sci</source><year>2012</year><volume>88</volume><issue>6</issue><fpage>226</fpage><lpage>249</lpage><pub-id pub-id-type="doi">10.2183/pjab.88.226</pub-id></element-citation></ref><ref id="CR40"><label>40.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Air</surname><given-names>GM</given-names></name></person-group><article-title>Influenza neuraminidase</article-title><source>Influenza Other Respir Viruses</source><year>2012</year><volume>6</volume><issue>4</issue><fpage>245</fpage><lpage>256</lpage><pub-id pub-id-type="doi">10.1111/j.1750-2659.2011.00304.x</pub-id><pub-id pub-id-type="pmid">22085243</pub-id></element-citation></ref><ref id="CR41"><label>41.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Maurer-Stroh</surname><given-names>S</given-names></name><etal></etal></person-group><article-title>A new common mutation in the hemagglutinin of the 2009 (H1N1) influenza a virus</article-title><source>PLoS Curr</source><year>2010</year><volume>2</volume><fpage>Rrn1162</fpage><pub-id pub-id-type="doi">10.1371/currents.RRN1162</pub-id><pub-id pub-id-type="pmid">20535229</pub-id></element-citation></ref><ref id="CR42"><label>42.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Graham</surname><given-names>M</given-names></name><etal></etal></person-group><article-title>Nationwide molecular surveillance of pandemic H1N1 influenza a virus genomes: Canada, 2009</article-title><source>PLoS One</source><year>2011</year><volume>6</volume><issue>1</issue><fpage>e16087</fpage><pub-id pub-id-type="doi">10.1371/journal.pone.0016087</pub-id><pub-id pub-id-type="pmid">21249207</pub-id></element-citation></ref><ref id="CR43"><label>43.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname><given-names>SS</given-names></name><etal></etal></person-group><article-title>Susceptibility of influenza a(H1N1)/pdm2009, seasonal a(H3N2) and B viruses to oseltamivir in Guangdong, China between 2009 and 2014</article-title><source>Sci Rep</source><year>2017</year><volume>7</volume><issue>1</issue><fpage>8488</fpage><pub-id pub-id-type="doi">10.1038/s41598-017-08282-6</pub-id><pub-id pub-id-type="pmid">28814737</pub-id></element-citation></ref><ref id="CR44"><label>44.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Butler</surname><given-names>J</given-names></name><etal></etal></person-group><article-title>Estimating the fitness advantage conferred by permissive neuraminidase mutations in recent oseltamivir-resistant a(H1N1)pdm09 influenza viruses</article-title><source>PLoS Pathog</source><year>2014</year><volume>10</volume><issue>4</issue><fpage>e1004065</fpage><pub-id pub-id-type="doi">10.1371/journal.ppat.1004065</pub-id><pub-id pub-id-type="pmid">24699865</pub-id></element-citation></ref><ref id="CR45"><label>45.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bragstad</surname><given-names>K</given-names></name><name><surname>Nielsen</surname><given-names>LP</given-names></name><name><surname>Fomsgaard</surname><given-names>A</given-names></name></person-group><article-title>The evolution of human influenza a viruses from 1999 to 2006: a complete genome study</article-title><source>Virol J</source><year>2008</year><volume>5</volume><fpage>40</fpage><pub-id pub-id-type="doi">10.1186/1743-422X-5-40</pub-id><pub-id pub-id-type="pmid">18325125</pub-id></element-citation></ref><ref id="CR46"><label>46.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Caini</surname><given-names>S</given-names></name><etal></etal></person-group><article-title>Epidemiology of seasonal influenza in the Middle East and North Africa regions, 2010-2016: circulating influenza a and B viruses and spatial timing of epidemics</article-title><source>Influenza Other Respir Viruses</source><year>2018</year><volume>12</volume><issue>3</issue><fpage>344</fpage><lpage>352</lpage><pub-id pub-id-type="doi">10.1111/irv.12544</pub-id><pub-id pub-id-type="pmid">29405575</pub-id></element-citation></ref><ref id="CR47"><label>47.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Simonsen</surname><given-names>L</given-names></name></person-group><article-title>The global impact of influenza on morbidity and mortality</article-title><source>Vaccine</source><year>1999</year><volume>17</volume><issue>Suppl 1</issue><fpage>S3</fpage><lpage>10</lpage><pub-id pub-id-type="doi">10.1016/S0264-410X(99)00099-7</pub-id><pub-id pub-id-type="pmid">10471173</pub-id></element-citation></ref><ref id="CR48"><label>48.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Simonsen</surname><given-names>L</given-names></name><etal></etal></person-group><article-title>Pandemic versus epidemic influenza mortality: a pattern of changing age distribution</article-title><source>J Infect Dis</source><year>1998</year><volume>178</volume><issue>1</issue><fpage>53</fpage><lpage>60</lpage><pub-id pub-id-type="doi">10.1086/515616</pub-id><pub-id pub-id-type="pmid">9652423</pub-id></element-citation></ref><ref id="CR49"><label>49.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bedford</surname><given-names>T</given-names></name><etal></etal></person-group><article-title>Integrating influenza antigenic dynamics with molecular evolution</article-title><source>Elife</source><year>2014</year><volume>3</volume><fpage>e01914</fpage><pub-id pub-id-type="doi">10.7554/eLife.01914</pub-id><pub-id pub-id-type="pmid">24497547</pub-id></element-citation></ref><ref id="CR50"><label>50.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Smith</surname><given-names>DJ</given-names></name><etal></etal></person-group><article-title>Mapping the antigenic and genetic evolution of influenza virus</article-title><source>Science</source><year>2004</year><volume>305</volume><issue>5682</issue><fpage>371</fpage><lpage>376</lpage><pub-id pub-id-type="doi">10.1126/science.1097211</pub-id><pub-id pub-id-type="pmid">15218094</pub-id></element-citation></ref><ref id="CR51"><label>51.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Klein</surname><given-names>EY</given-names></name><etal></etal></person-group><article-title>Influenza a H1N1 pandemic strain evolution--divergence and the potential for antigenic drift variants</article-title><source>PLoS One</source><year>2014</year><volume>9</volume><issue>4</issue><fpage>e93632</fpage><pub-id pub-id-type="doi">10.1371/journal.pone.0093632</pub-id><pub-id pub-id-type="pmid">24699432</pub-id></element-citation></ref><ref id="CR52"><label>52.</label><mixed-citation publication-type="other">CDC, Seasonal Influenza Vaccine Effectiveness, 2004-2018. 2004-2018, Centers for Disease Control and Prevention.</mixed-citation></ref><ref id="CR53"><label>53.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname><given-names>Y</given-names></name><etal></etal></person-group><article-title>Altered receptor specificity and cell tropism of D222G hemagglutinin mutants isolated from fatal cases of pandemic a(H1N1) 2009 influenza virus</article-title><source>J Virol</source><year>2010</year><volume>84</volume><issue>22</issue><fpage>12069</fpage><lpage>12074</lpage><pub-id pub-id-type="doi">10.1128/JVI.01639-10</pub-id><pub-id pub-id-type="pmid">20826688</pub-id></element-citation></ref><ref id="CR54"><label>54.</label><mixed-citation publication-type="other">Kilander A, et al. Observed association between the HA1 mutation D222G in the 2009 pandemic influenza a(H1N1) virus and severe clinical outcome, Norway 2009-2010. Euro Surveill. 2010;15(9).</mixed-citation></ref><ref id="CR55"><label>55.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gamblin</surname><given-names>SJ</given-names></name><etal></etal></person-group><article-title>The structure and receptor binding properties of the 1918 influenza hemagglutinin</article-title><source>Science</source><year>2004</year><volume>303</volume><issue>5665</issue><fpage>1838</fpage><lpage>1842</lpage><pub-id pub-id-type="doi">10.1126/science.1093155</pub-id><pub-id pub-id-type="pmid">14764886</pub-id></element-citation></ref><ref id="CR56"><label>56.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Maines</surname><given-names>TR</given-names></name><etal></etal></person-group><article-title>Transmission and pathogenesis of swine-origin 2009 a(H1N1) influenza viruses in ferrets and mice</article-title><source>Science</source><year>2009</year><volume>325</volume><issue>5939</issue><fpage>484</fpage><lpage>487</lpage><pub-id pub-id-type="doi">10.1126/science.1177238</pub-id><pub-id pub-id-type="pmid">19574347</pub-id></element-citation></ref><ref id="CR57"><label>57.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>de Vries</surname><given-names>RP</given-names></name><etal></etal></person-group><article-title>Evolution of the hemagglutinin protein of the new pandemic H1N1 influenza virus: maintaining optimal receptor binding by compensatory substitutions</article-title><source>J Virol</source><year>2013</year><volume>87</volume><issue>24</issue><fpage>13868</fpage><lpage>13877</lpage><pub-id pub-id-type="doi">10.1128/JVI.01955-13</pub-id><pub-id pub-id-type="pmid">24109242</pub-id></element-citation></ref><ref id="CR58"><label>58.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jayaraman</surname><given-names>A</given-names></name><etal></etal></person-group><article-title>A single base-pair change in 2009 H1N1 hemagglutinin increases human receptor affinity and leads to efficient airborne viral transmission in ferrets</article-title><source>PLoS One</source><year>2011</year><volume>6</volume><issue>3</issue><fpage>e17616</fpage><pub-id pub-id-type="doi">10.1371/journal.pone.0017616</pub-id><pub-id pub-id-type="pmid">21407805</pub-id></element-citation></ref><ref id="CR59"><label>59.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Luksza</surname><given-names>M</given-names></name><name><surname>Lassig</surname><given-names>M</given-names></name></person-group><article-title>A predictive fitness model for influenza</article-title><source>Nature</source><year>2014</year><volume>507</volume><issue>7490</issue><fpage>57</fpage><lpage>61</lpage><pub-id pub-id-type="doi">10.1038/nature13087</pub-id><pub-id pub-id-type="pmid">24572367</pub-id></element-citation></ref><ref id="CR60"><label>60.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wei</surname><given-names>CJ</given-names></name><etal></etal></person-group><article-title>Cross-neutralization of 1918 and 2009 influenza viruses: role of glycans in viral evolution and vaccine design</article-title><source>Sci Transl Med</source><year>2010</year><volume>2</volume><issue>24</issue><fpage>24ra21</fpage><pub-id pub-id-type="doi">10.1126/scitranslmed.3000799</pub-id></element-citation></ref><ref id="CR61"><label>61.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ison</surname><given-names>MG</given-names></name></person-group><article-title>Antivirals and resistance: influenza virus</article-title><source>Curr Opin Virol</source><year>2011</year><volume>1</volume><issue>6</issue><fpage>563</fpage><lpage>573</lpage><pub-id pub-id-type="doi">10.1016/j.coviro.2011.09.002</pub-id><pub-id pub-id-type="pmid">22440914</pub-id></element-citation></ref><ref id="CR62"><label>62.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Meijer</surname><given-names>A</given-names></name><etal></etal></person-group><article-title>Oseltamivir-resistant influenza virus a (H1N1), Europe, 2007-08 season</article-title><source>Emerg Infect Dis</source><year>2009</year><volume>15</volume><issue>4</issue><fpage>552</fpage><lpage>560</lpage><pub-id pub-id-type="doi">10.3201/eid1504.181280</pub-id><pub-id pub-id-type="pmid">19331731</pub-id></element-citation></ref><ref id="CR63"><label>63.</label><mixed-citation publication-type="other">World Health Organization. Positive samples update. World Health Organization. Number of specimens positive for influenza by subtypes from week no. 17 (2009) to 21 (2010). 2010 [cited 2018; Available from: <ext-link ext-link-type="uri" xlink:href="http://www.who.int/csr/disease/swine%20u/%20Virologicaldata2010_06_11.pdf">http://www.who.int/csr/disease/swine u/ Virologicaldata2010_06_11.pdf</ext-link>.</mixed-citation></ref><ref id="CR64"><label>64.</label><mixed-citation publication-type="other">Lackenby, A., et al., Continued emergence and changing epidemiology of oseltamivir-resistant influenza a(H1N1)2009 virus, United Kingdom, winter 2010/11. Euro Surveill. 2011;16(5).</mixed-citation></ref><ref id="CR65"><label>65.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Storms</surname><given-names>AD</given-names></name><etal></etal></person-group><article-title>Oseltamivir-resistant pandemic (H1N1) 2009 virus infections, United States, 2010-11</article-title><source>Emerg Infect Dis</source><year>2012</year><volume>18</volume><issue>2</issue><fpage>308</fpage><lpage>311</lpage><pub-id pub-id-type="doi">10.3201/eid1802.111466</pub-id><pub-id pub-id-type="pmid">22305467</pub-id></element-citation></ref><ref id="CR66"><label>66.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gubareva</surname><given-names>Larisa V.</given-names></name><name><surname>Besselaar</surname><given-names>Terry G.</given-names></name><name><surname>Daniels</surname><given-names>Rod S.</given-names></name><name><surname>Fry</surname><given-names>Alicia</given-names></name><name><surname>Gregory</surname><given-names>Vicki</given-names></name><name><surname>Huang</surname><given-names>Weijuan</given-names></name><name><surname>Hurt</surname><given-names>Aeron C.</given-names></name><name><surname>Jorquera</surname><given-names>Patricia A.</given-names></name><name><surname>Lackenby</surname><given-names>Angie</given-names></name><name><surname>Leang</surname><given-names>Sook-Kwan</given-names></name><name><surname>Lo</surname><given-names>Janice</given-names></name><name><surname>Pereyaslov</surname><given-names>Dmitriy</given-names></name><name><surname>Rebelo-de-Andrade</surname><given-names>Helena</given-names></name><name><surname>Siqueira</surname><given-names>Marilda M.</given-names></name><name><surname>Takashita</surname><given-names>Emi</given-names></name><name><surname>Odagiri</surname><given-names>Takato</given-names></name><name><surname>Wang</surname><given-names>Dayan</given-names></name><name><surname>Zhang</surname><given-names>Wenqing</given-names></name><name><surname>Meijer</surname><given-names>Adam</given-names></name></person-group><article-title>Global update on the susceptibility of human influenza viruses to neuraminidase inhibitors, 2015–2016</article-title><source>Antiviral Research</source><year>2017</year><volume>146</volume><fpage>12</fpage><lpage>20</lpage><pub-id pub-id-type="doi">10.1016/j.antiviral.2017.08.004</pub-id><pub-id pub-id-type="pmid">28802866</pub-id></element-citation></ref><ref id="CR67"><label>67.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hurt</surname><given-names>AC</given-names></name><etal></etal></person-group><article-title>Community transmission of oseltamivir-resistant a(H1N1)pdm09 influenza</article-title><source>N Engl J Med</source><year>2011</year><volume>365</volume><issue>26</issue><fpage>2541</fpage><lpage>2542</lpage><pub-id pub-id-type="doi">10.1056/NEJMc1111078</pub-id><pub-id pub-id-type="pmid">22204735</pub-id></element-citation></ref><ref id="CR68"><label>68.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Takashita</surname><given-names>E</given-names></name><name><surname>Ejima</surname><given-names>M</given-names></name><name><surname>Itoh</surname><given-names>R</given-names></name><name><surname>Miura</surname><given-names>M</given-names></name><name><surname>Ohnishi</surname><given-names>A</given-names></name><name><surname>Nishimura</surname><given-names>H</given-names></name><name><surname>Odagiri</surname><given-names>T</given-names></name><name><surname>Tashiro</surname><given-names>M</given-names></name></person-group><article-title>A community cluster of influenza A(H1N1)pdm09 virus exhibiting cross-resistance to oseltamivir and peramivir in Japan, November to December 2013</article-title><source>Eurosurveillance</source><year>2014</year><volume>19</volume><issue>1</issue><fpage>20666</fpage><pub-id pub-id-type="doi">10.2807/1560-7917.ES2014.19.1.20666</pub-id><pub-id pub-id-type="pmid">24434172</pub-id></element-citation></ref><ref id="CR69"><label>69.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname><given-names>LF</given-names></name><etal></etal></person-group><article-title>Cluster of oseltamivir-resistant 2009 pandemic influenza a (H1N1) virus infections on a hospital ward among immunocompromised patients--North Carolina, 2009</article-title><source>J Infect Dis</source><year>2011</year><volume>203</volume><issue>6</issue><fpage>838</fpage><lpage>846</lpage><pub-id pub-id-type="doi">10.1093/infdis/jiq124</pub-id><pub-id pub-id-type="pmid">21343149</pub-id></element-citation></ref><ref id="CR70"><label>70.</label><mixed-citation publication-type="other">Al-Tawfiq JA, Benkouiten S, Memish ZA. A systematic review of emerging respiratory viruses at the hajj and possible coinfection with Streptococcus pneumoniae. Travel Med Infect Dis. 2018.</mixed-citation></ref></ref-list></back></pmc></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Sante/explor/H2N2V1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000165 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000165 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Sante
|area= H2N2V1
|flux= Pmc
|étape= Corpus
|type= RBID
|clé=
|texte=
}}
| This area was generated with Dilib version V0.6.33. |  |