Skip directly to site content Skip directly to page options Skip directly to A-Z link Skip directly to A-Z link Skip directly to A-Z link
Volume 25, Number 1—January 2019

Avian Influenza A(H9N2) Virus in Poultry Worker, Pakistan, 2015

Muzaffar Ali, Tahir Yaqub, Nadia Mukhtar, Muhammad Imran, Aamir Ghafoor, Muhammad Furqan Shahid, Muhammad Naeem, Munir Iqbal, Gavin J.D. SmithComments to Author , and Yvonne C.F. Su
Author affiliations: University of Veterinary and Animal Sciences, Lahore, Pakistan (M. Ali, T. Yaqub, M. Imran, A. Ghafoor, M.F. Shahid); Health Care Department, Government of Punjab, Lahore (N. Mukhtar); Bahauddin Zakariya University, Multan, Pakistan (M. Naeem); The Pirbright Institute, Compton Laboratory, Newbury, UK (M. Iqbal); Duke University, Durham, North Carolina, USA (G.J.D. Smith); Duke-National University Singapore Medical School, Singapore (G.J.D. Smith, Y.C.F. Su)

Cite This Article


Avian influenza A(H9N2) virus isolated from a poultry worker in Pakistan in 2015 was closely related to viruses detected in poultry farms. Observed mutations in the hemagglutinin related to receptor-binding affinity and antigenicity could affect cross-reactivity with prepandemic H9N2 vaccine strains.

Influenza A(H9N2) virus circulates in domestic poultry, and outbreaks have been recorded since the early 1990s in China (1). In Pakistan, H9N2 virus was first detected in a poultry outbreak in 1998; subsequent outbreaks have led to increased genetic diversifications of distinct viral lineages descended from H9N2 G1 lineage viruses (2). Serologic studies of H9 virus among persons of different occupations in Pakistan who had direct exposure to poultry (e.g., poultry workers, vaccinators, veterinarians) have shown high rates of seropositivity (30%–85%) (35). Although no human infection with H9N2 virus has been reported from Pakistan, sporadic clinical cases of H9N2 virus infection in humans have been reported in China (6), Hong Kong (7), and Bangladesh (8). We report the isolation of H9N2 virus from a poultry worker during avian influenza virus surveillance in Pakistan.

The Study

During January 2015–June 2016, avian influenza virus (AIV) surveillance was conducted in poultry farms throughout 19 districts of Punjab Province, Pakistan (9). In addition, after obtaining written informed consent, we collected 117 nasal swab specimens from male poultry workers 25–35 years of age. The Institutional Ethical Review Board at the Institute of Public Health (Lahore, Pakistan) reviewed and approved the study protocol.

Sterile swabs were used to take nasal swab samples from humans; the swabs were placed in sterile tubes containing 2 mL of viral transport media. To prevent cross-contamination between samples, human and chicken samples were collected in different zip-sealable plastic bags and transported on ice to the laboratory. Human and chicken samples were also processed and cultured separately on different dates. Individual human samples were inoculated in 9-day-old embryonated chicken eggs and amnio-allantoic fluid (AAF) harvested after 48 h incubation. Harvested AAF was first tested by hemagglutination assay, and positive AAF screened for H5, H7, H9, and Newcastle disease virus (NDV) using hemagglutination inhibition (HI) assay as previously described (9).

We detected 1 H9-positive sample collected from a poultry worker in Narowal District, Punjab, where only 1 (1.1%) of 88 chicken samples was H9N2-positive (9), and the flock from the same farm was H9 negative. The 36-year-old worker did not display major signs of influenza-like illness, which is typical of the mild to no symptoms shown in H9N2 virus infection (10). No human samples were positive for H5 or H7 virus. The virus isolate was confirmed as H9N2 by reverse transcription PCR, and the hemagglutinin (HA) and neuraminidase (NA) genes were sequenced as described (9). We sequenced the HA and NA genes of 8 additional viruses isolated from chickens in different districts of Punjab (Table 1).


Thumbnail of Evolutionary relationships of the influenza A virus H9-HA gene isolated from avian and human hosts, Pakistan, 1998–2016. The phylogeny was generated using the uncorrelated lognormal relaxed molecular clock, the SRD06 codon position model, HKY85 plus gamma substitution model, and a Gaussian Markov random field (GMRF) Bayesian skyride in BEAST version 1.8.4. Two independent runs of 100 million Markov chain Monte Carlo generations were performed. Horizontal node bars represent the 95%

Figure. Evolutionary relationships of the influenza A virus H9-HA gene isolated from avian and human hosts, Pakistan, 1998–2016. The phylogeny was generated using the uncorrelated lognormal relaxed molecular clock, the SRD06 codon...

We analyzed 2,751 H9NX and 8,059 HXN2 avian and human virus sequences (collected during 1963–2017) obtained from GenBank and the GISAID database ( and reconstructed large H9-HA and N2-NA phylogenies using maximum-likelihood (ML) analysis. We then subsampled the datasets to a final dataset of 145 sequences for each H9-HA and N2-NA. Reconstruction of temporal phylogenies was performed as previously described (11). H9-HA (Figure; Appendix Figure 1) and N2-NA (Appendix Figure 2) phylogenies show that, within the G1 lineage, the human H9N2 isolate (A/Pakistan/486/2015) formed a strongly supported monophyletic group with chicken H9N2 viruses collected in Pakistan from 2015 to 2016. The HA and NA genes of A/Pakistan/486/2015 (H9N2) showed close genetic resemblance (99.8%–99.9% nt identity) with 2015–2016 Pakistan chicken viruses, most likely indicating direct cross-species transmission from poultry to human. These 2015–2016 Pakistan viruses are further grouped within sublineage B2 (2), which also includes viruses isolated during 2008–2014 from Afghanistan, Iran, and Pakistan. The mean times to the most recent common ancestor (tMRCAs) of the HA and NA genes for the 2015–2016 Pakistan strains were both estimated as late 2010 (Table 2), and sublineage B2 mean tMRCAs were estimated as 2005–2006. Together, these results indicate the circulation of unsampled virus diversity in poultry during the past decade, hindering our ability to investigate virus transmission and evolution of H9N2 virus in Pakistan and highlighting the need for systematic surveillance.

All 2015–2016 Pakistan lineage viruses possessed a leucine (L) residue at position 226 (H3 numbering) in the HA receptor-binding site (RBS), which is associated with greater affinity for α-2,6 binding. However, HA residue 228 (H3 numbering) in the RBS was glutamine (G), known for avian-like receptor specificity. All 2015–2016 Pakistan lineage viruses possessed a low pathogenic motif (KSSR/G) in the HA cleavage site. To further characterize the recent evolution of H9N2 viruses in Pakistan, we mapped the HA and NA amino acid substitutions onto the phylogenetic trees and used Mixed Effects Model of Evolution (MEME) analysis to detect amino acids under positive selection, as previously described (11). Three amino acid substitutions were present at the ancestral node of the HA of the 2015–2016 Pakistan lineage, T103I, L150N, and V393I (mature H9 protein numbering), although only residue 150 (H3 numbering 160) was under positive selection. One additional substitution, T186A (H3 numbering 196), was present in the human isolate A/Pakistan/486/2015 (H9N2). Both of these mutations are situated in the antigenic region of the H9-HA globular head and correspond to H3N2 epitope B (12). The biologic function of mutations at the remaining residues is not known. In addition, the HA phylogeny also indicated that 2015–2016 Pakistan viruses had diverged from lineages containing the 3 H9N2 G1 lineage candidate vaccine viruses (A/Hong Kong/1073/99, A/Hong Kong/33982/2009, and A/Bangladesh/994/2011) (Appendix Figure 1). These vaccine candidates may therefore not provide efficient protection against the avian 2015–2016 Pakistan lineage, although this assumption must be confirmed by antigenic assays.

No deletion in the NA gene, associated with aquatic to terrestrial host adaptation and increased replication in ferrets (13), was observed in the 2015–2016 Pakistan lineage viruses. Numerous substitutions were observed in the 2015–2016 Pakistan lineage viruses (Appendix Figure 2). Of these, the Q39R, K47E, and I62T mutations are functionally key residues in the NA stalk that may be associated with host adaptation and virus virulence (13). The mutations at 372 and 401 residues may be responsible for hemadsorption activity of the NA (14). The functional importance of the remaining mutations is unknown, including V263I that was under positive selection.


Our detection and isolation of H9N2 virus from a poultry worker in Pakistan highlights the potential for cross-species transmission of H9 viruses in the country. The World Health Organization considers avian H9N2 viruses a consistent pandemic threat because they are widespread in poultry and cause sporadic infection in humans. H9N2 viruses have been central to the generation of other viruses of pandemic concern and have contributed the internal genes to both H5 and H7 viruses in China (15). Although the subtype is relatively well studied in China, investigation in other countries is generally limited in scope. Within Pakistan, H9N2 viruses in chickens have circulated endemically for at least a decade, yet systematic surveillance is lacking.

Our results show continued diversification of H9N2 viruses in Pakistan; viruses isolated during 2015–2016 formed a distinct clade to earlier viruses from Afghanistan, Iran, and Pakistan isolated during 2008–2014. Dating analysis further estimated the tMRCA of the 2015–2016 Pakistan viruses as late 2010, indicating at least 5 years of unsampled virus diversity that circulated in poultry. We also observed mutations in HA related to changes in receptor-binding affinity and antigenicity that could affect cross-reactivity with the World Health Organization–recommended prepandemic H9N2 vaccine strains. None of the 3 G1 candidate vaccine viruses are closely related to strains from Pakistan.

Phylogenetic relationships indicate H9N2 virus transmission across South Asia and the Middle East, where the persistence and circulation of AIV are poorly understood. Increased surveillance in wild bird populations, poultry farms and markets, and occupationally exposed workers is needed in these regions to identify the emergence of antigenic variants and to maintain up-to-date H9 vaccine candidates.

Mr. Ali is a PhD student at the University of Veterinary and Animal Sciences, Lahore, Pakistan, and performed part of this work while an International Research Support Initiative Program fellow at Duke-National University Singapore Medical School. His research focuses on the evolution of avian influenza viruses.



We thank the Influenza Laboratory at the Department of Microbiology, University of Veterinary and Animal Sciences Lahore, Pakistan, for providing facilities, funded by Biotechnology and Biological Sciences Research Council, UK (project no. BB/L018853/1).

The study was supported by the Duke-National University Singapore Signature Research Programme funded by the Ministry of Health, Singapore, and by contract HHSN272201400006C from the National Institute of Allergy and Infectious Disease, National Institutes of Health, US Department of Health and Human Services. M.A. was supported by a fellowship from the International Research Support Initiative Program administered by the Higher Education Commission of Pakistan.



  1. Li  C, Yu  K, Tian  G, Yu  D, Liu  L, Jing  B, et al. Evolution of H9N2 influenza viruses from domestic poultry in Mainland China. Virology. 2005;340:7083. DOIPubMedGoogle Scholar
  2. Lee  D-H, Swayne  DE, Sharma  P, Rehmani  SF, Wajid  A, Suarez  DL, et al. H9N2 low pathogenic avian influenza in Pakistan (2012-2015). Vet Rec Open. 2016;3:e000171. DOIPubMedGoogle Scholar
  3. Ahad  A, Rabbani  M, Yaqub  T, Younus  M, Mahmood  A, Zubair  M, et al. Serosurveillance to H9 and H7 avian influenza virus among poultry workers in Punjab Province, Pakistan. Pak Vet J. 2013;33:10712.
  4. Rasheed  M, Rehmani  SF, Iqbal  M, Ahmad  A, Akhtar  F, Akhtar  R, et al. Seropositivity to avian influenza virus subtype H9N2 among human population of selected districts of Punjab, Pakistan. Journal of Infection and Molecular Biology. 2013;1:324.
  5. Ahad  A, Thornton  RN, Rabbani  M, Yaqub  T, Younus  M, Muhammad  K, et al. Risk factors for H7 and H9 infection in commercial poultry farm workers in provinces within Pakistan. Prev Vet Med. 2014;117:6104. DOIPubMedGoogle Scholar
  6. Pan  Y, Cui  S, Sun  Y, Zhang  X, Ma  C, Shi  W, et al. Human infection with H9N2 avian influenza in northern China. Clin Microbiol Infect. 2018;24:3213. DOIPubMedGoogle Scholar
  7. Butt  KM, Smith  GJ, Chen  H, Zhang  LJ, Leung  YH, Xu  KM, et al. Human infection with an avian H9N2 influenza A virus in Hong Kong in 2003. J Clin Microbiol. 2005;43:57607. DOIPubMedGoogle Scholar
  8. Food and Agriculture Organization. EMPRES-I Global Animal Disease Information System [cited 2018 Apr 3].
  9. Ali  M, Yaqub  T, Mukhtar  N, Imran  M, Ghafoor  A, Shahid  MF, et al. Prevalence and phylogenetics of H9N2 in backyard and commercial poultry in Pakistan. Avian Dis. 2018. In press.
  10. Khan  SU, Anderson  BD, Heil  GL, Liang  S, Gray  GC. A systematic review and meta-analysis of the seroprevalence of influenza A(H9N2) infection among humans. J Infect Dis. 2015;212:5629. DOIPubMedGoogle Scholar
  11. Su  YC, Bahl  J, Joseph  U, Butt  KM, Peck  HA, Koay  ES, et al. Phylodynamics of H1N1/2009 influenza reveals the transition from host adaptation to immune-driven selection. Nat Commun. 2015;6:7952. DOIPubMedGoogle Scholar
  12. Wan  Z, Ye  J, Xu  L, Shao  H, Jin  W, Qian  K, et al. Antigenic mapping of the hemagglutinin of an H9N2 avian influenza virus reveals novel critical amino acid positions in antigenic sites. J Virol. 2014;88:3898901. DOIPubMedGoogle Scholar
  13. Diederich  S, Berhane  Y, Embury-Hyatt  C, Hisanaga  T, Handel  K, Cottam-Birt  C, et al. Hemagglutinin-neuraminidase balance influences the virulence phenotype of a recombinant H5N3 influenza A virus possessing a polybasic HA0 cleavage site. J Virol. 2015;89:1072434. DOIPubMedGoogle Scholar
  14. Uhlendorff  J, Matrosovich  T, Klenk  H-D, Matrosovich  M. Functional significance of the hemadsorption activity of influenza virus neuraminidase and its alteration in pandemic viruses. Arch Virol. 2009;154:94557. DOIPubMedGoogle Scholar
  15. Zhu  H, Lam  TT-Y, Smith  DK, Guan  Y. Emergence and development of H7N9 influenza viruses in China. Curr Opin Virol. 2016;16:10613. DOIPubMedGoogle Scholar




Cite This Article

DOI: 10.3201/eid2501.180618

Original Publication Date: December 04, 2018

Table of Contents – Volume 25, Number 1—January 2019

EID Search Options
presentation_01 Advanced Article Search – Search articles by author and/or keyword.
presentation_01 Articles by Country Search – Search articles by the topic country.
presentation_01 Article Type Search – Search articles by article type and issue.



Please use the form below to submit correspondence to the authors or contact them at the following address:

Gavin J.D. Smith, Duke-National University Singapore Medical School, Programme in Emerging Infectious Diseases, 8 College Rd, Singapore 169857, Singapore

Send To

10000 character(s) remaining.


Page created: December 18, 2018
Page updated: December 18, 2018
Page reviewed: December 18, 2018
The conclusions, findings, and opinions expressed by authors contributing to this journal do not necessarily reflect the official position of the U.S. Department of Health and Human Services, the Public Health Service, the Centers for Disease Control and Prevention, or the authors' affiliated institutions. Use of trade names is for identification only and does not imply endorsement by any of the groups named above.