Emerg Infect DisEmerging Infect. DisEIDEmerging Infectious Diseases1080-60401080-6059Centers for Disease Control and Prevention23968540381092513-033610.3201/eid1909.130336ResearchResearchAntigenic and Molecular Characterization of Avian Influenza A(H9N2) Viruses, BangladeshAvian Influenza (H9N2), BangladeshShanmuganathamKarthikFeerozMohammed M.Jones-EngelLisaSmithGavin J.D.FourmentMathieuWalkerDavidMcClenaghanLauraAlamS.M. RabiulHasanM. KamrulSeilerPatrickFranksJohnDannerAngieBarmanSubrataMcKenziePamelaKraussScottWebbyRichard J.WebsterRobert G.St. Jude Children’s Research Hospital, Memphis, Tennessee, USA (K. Shanmuganatham, D. Walker, L. McClenaghan, P. Seiler, J. Franks, A. Danner, S. Barman, P. McKenzie, S. Krauss, R.J. Webby, R.G. Webster); Jahangirnagar University, Dhaka, Bangladesh (M.M. Feeroz, S.M.R. Alam, M.K. Hasan); University of Washington, Seattle, Washington, USA (L. Jones-Engel); Duke-NUS Graduate Medical School, Singapore (G.J.D. Smith, M. Fourment); Duke Global Health Institute, Durham, North Carolina, USA (G.J.D. Smith)Address for correspondence: Robert G. Webster, St. Jude Children’s Research Hospital, Department of Infectious Diseases, MS 330, 262 Danny Thomas Pl, Memphis, TN 38105, USA; email: robert.webster@stjude.org9201319913931402

Human infection with avian influenza A(H9N2) virus was identified in Bangladesh in 2011. Surveillance for influenza viruses in apparently healthy poultry in live-bird markets in Bangladesh during 2008–2011 showed that subtype H9N2 viruses are isolated year-round, whereas highly pathogenic subtype H5N1 viruses are co-isolated with subtype H9N2 primarily during the winter months. Phylogenetic analysis of the subtype H9N2 viruses showed that they are reassortants possessing 3 gene segments related to subtype H7N3; the remaining gene segments were from the subtype H9N2 G1 clade. We detected no reassortment with subtype H5N1 viruses. Serologic analyses of subtype H9N2 viruses from chickens revealed antigenic conservation, whereas analyses of viruses from quail showed antigenic drift. Molecular analysis showed that multiple mammalian-specific mutations have become fixed in the subtype H9N2 viruses, including changes in the hemagglutinin, matrix, and polymerase proteins. Our results indicate that these viruses could mutate to be transmissible from birds to mammals, including humans.

Keywords: influenzaavian influenza A(H9N2) virusesvirusesH9N2H5N1H7N3molecular characterizationreassortantsBangladeshzoonosestransmissionavian influenzarespiratory infections

Initially infecting poultry, avian influenza A(H9N2) viruses have been sporadically identified in pigs and humans, which suggests that some of these viruses have adapted to bind mammalian host receptors or have acquired mutations that increase mammalian receptor specificity (13). Human infection with avian influenza A(H9N2) virus was initially identified in Hong Kong and China in 1999 (4); in 2011, infection with this subtype was reported for a patient in Bangladesh (5). Detection of these viruses in humans outside of China highlights the necessity and urgency for comprehensive surveillance because of the viruses’ expanding host range.

Phylogenetically, avian influenza A(H9N2) viruses can be grouped into 3 distinct sublineages represented by their prototype strains: A/Qa/HK/G1/97 (G1-like), A/Dk/HK/Y280/97 (Y280-like), and A/Ck/Korea/38349-p96323/96 (Korean-like) (1,6,7). Genetic and antigenic analyses of subtype H9N2 isolates from the past 2 decades have shown that these viruses are gradually evolving from the Eurasian lineage into several distinct sublineages and are becoming established in domestic poultry (718). Phylogenetic analyses of subtype H9N2 viruses isolated in China and the Middle East have shown that these viruses have undergone reassortment with other subtypes to generate multiple novel genotypes consisting of gene segments from different lineages (7,11,13,19,20).

Worldwide, Bangladesh is among countries with the highest numbers of reported outbreaks of highly pathogenic avian influenza (HPAI) (H5N1) (21). Since an initial outbreak in February 2007, Bangladesh has reported 550 outbreaks of infection with HPAI (H5N1) virus (493 at commercial farms and 57 among backyard poultry) (2224). In Bangladesh, live-bird markets are the most common outlets for purchase of poultry and poultry meat; an estimated 95% of poultry meat and eggs sold in the country are sold at these markets (25). Previous surveillance conducted at live-bird markets in Bangladesh found that avian influenza virus (AIV) is prevalent (23%); the low pathogenicity H9N2 subtype predominated, but other subtypes were isolated, including H5N1, H1N2, H1N3, H3N6, H4N2, and H10N7 (26). However, this surveillance report did not include information about the molecular properties of circulating subtype H9N2 viruses.

We reviewed surveillance data and conducted molecular and genetic analyses of influenza A(H9N2) viruses circulating among poultry in Bangladesh. Our study had 3 primary goals: 1) characterize the antigenic and molecular properties of subtype H9N2 isolates; 2) define genetic and phylogenetic relationships between the genes identified in these viruses and those of other AIVs; and 3) determine whether these viruses have acquired genomic changes that could facilitate transmission from avian to mammalian hosts.

Materials and MethodsSample Collection and Screening

Surveillance for influenza viruses in poultry began in Bangladesh in November 2008. Each month, trained personnel collected 200–600 samples from apparently healthy domestic live birds (chickens, quail, pigeons, ducks, and turkeys) at retail markets, a pet market, chicken layer and duck farms, and wild birds (50–200 samples from each location). Samples consisted of oropharyngeal, cloacal, and environmental (fecal matter and water samples from cages and fecal digesters) swab specimens.

During November 2008–August 2011, a total of 17,438 samples (3,078 oropharyngeal, 3,377 cloacal, 10,983 environmental) were collected (Table 1). Screening for AIVs was performed as described (26). Briefly, all samples were subjected to real-time reverse transcription PCR by using influenza A–specific primer and probes; samples with positive results were injected into egg to confirm the presence or absence of virus. A total of 734 AIV (H9N2) isolates were extracted from the 17,438 samples collected (Table 1).

Sources and subtypes of avian influenza viruses isolated in Bangladesh, 2008–2011
SitePrimary species from which samples were obtained*Primary species from which sample was influenza positiveTotal no. samplesSubtypes
H1N2, H1N3, H3N6, H3N8, H4N2, H10N7H5N1H9N2
Live bird market-1
Quail
Quail
821

20
37
Live bird market-2 Chickens, ducks, pigeons
Chickens1,691121385

Duck



2
Live bird market-3ChickensChickens1,877738360


Ducks



2
Live bird markets-4
Chickens
Chickens
46


12
Live bird markets-5ChickensChickens2,6771420205


Duck



2
Pet market
Pet and wild birds
Quail
1,859
3
1
20
Farm-1DucksDucks10024
Farm-2ChickensChickens3,94529
Farm-3Chickens40
Farm-4Ducks100
Farm-5Ducks100
Farm-6
Ducks

901



Wild birds, environmental-1Wild birds836
Wild birds, environmental -2Wild birds471
Wild ducks, environmental-1Water1,724
Wild ducks, environmental-2
Ducks

250



Total 17,4386292734

*Other species that were sampled at some locations included turkeys, wild birds, and exotic birds.

Virus Isolation and Propagation

To determine the genetic and evolutionary diversity of the AIV (H9N2) viruses circulating in Bangladesh, we selected 44 of the 734 isolates for further analysis (Table 2); these isolates were representative of the location, time, species, and sample types from which they were isolated. With the exception of 1 isolate from a chicken on a farm, all isolates examined were from chicken, duck, and quail samples collected at retail markets. The viruses were propagated in 10-day-old embryonated chicken eggs, and initial subtyping was done by sequencing the hemagglutinin (HA) and neuraminidase (NA) genes, as described (26).

Avian influenza A(H9N2) virus isolates from animals in live bird markets, Bangladesh, 2008–2011*
Isolate Host common nameSample type Mixtures†
A/Env/BD/907/2009 (quail)QuailFNDV
A/Env/BD/1041/2009 (duck)DuckFH5
A/Dk/BD/1231/2009DuckOPH5
A/Ck/BD/2075/2009ChickenOPH5
A/Pigeon/BD/4303/2009PigeonOPNDV
A/Env/BD/5144/2009 (quail)QuailFH5
A/Ck/BD/5209/2009ChickenOPNDV
A/Env/BD/5745/2010 (duck)DuckFH5
A/Env/BD/8202/2010 (chicken)ChickenF
A/Ck/BD/8411/2010ChickenOP
A/Ck/BD/8413/2010ChickenOP
A/Ck/BD/8415/2010ChickenOP
A/Env/BD/8463/2010 (chicken)ChickenW
A/Env/BD/8465/2010 (chicken)ChickenW
A/Ck/BD/8725/2010ChickenOP
A/Ck/BD/8731/2010ChickenOP
A/Ck/BD/8996/2010ChickenCNDV
A/Ck/BD/9029/2010ChickenOPNDV
A/Env/BD/9306/2010 (parrot)ParrotFNDV
A/Ck/BD/9334/2010ChickenOP
A/Env/BD/9350/2010 (chicken)ChickenWNDV
A/Env/BD/9457/2010 (chicken)ChickenF
A/Env/BD/10234/2011 (chicken)ChickenW
A/Env/BD/10306/2011 (quail)QuailFNDV
A/Env/BD/10307/2011 (quail)QuailFNDV
A/Env/BD/10313/2011 (quail)QuailF
A/Env/BD/10316/2011QuailF
A/Ck/BD/10401/2011ChickenOP
A/Ck/BD/10411/2011ChickenOP
A/Ck/BD/10450/2011ChickenCNDV
A/Ck/BD/10897/2011ChickenOP
A/Ck/BD/11154/2011ChickenCNDV
A/Env/BD/11173/2011 (chicken)ChickenF
A/Ck/BD/11309/2011ChickenOP
A/Ck/BD/11315/2011ChickenOP
A/Env/BD/11597/2011 (chicken)ChickenW
A/Env/BD/12068/2011 (pigeon)PigeonFNDV
A/Env/BD/12077/2011 (turkey)TurkeyFNDV
A/Env/BD/12093/2011 (quail)QuailF
A/Env/BD/12103/2011 (quail)QuailF
A/Env/BD/12116/2011 (quail)QuailF
A/Env/BD/12119/2011 (quail)QuailF
A/Ck/BD/13916/2011ChickenOPH5
A/Ck/BD/13962/2011ChickenOPH5

*F, fecal; NDV, Newcastle disease virus; OP, oropharyngeal; W, water; C, cloacal.
†Some hosts infected with influenza (H9N2) viruses were co-infected with H5N1 or NDV; thus, mixtures of isolates were obtained from some birds. Species in parentheses indicate that the samples were collected from the cages where these particular birds were caged.

Hemagglutination Inhibition Assay

Polyclonal serum samples were obtained from ferrets that had been inoculated with 1 of the following: prototype viruses from different influenza (H9N2) lineages; influenza (H9N2) isolates from chickens or quail collected in Bangladesh; or a prototype human isolate provided by the US Centers for Disease Control and Prevention. We then conducted hemagglutination inhibition (HI) assays for these serum samples by using 0.5% chicken erythrocytes, as described (27).

Phylogenetic Analysis and Molecular Characterization

Viral RNA extraction and reverse transcription PCR were performed as described (28,29). DNA sequencing using specific primers was completed in the Hartwell Center at St. Jude Children’s Research Hospital (Memphis, TN, USA). Samples were analyzed by using the BigDye Terminator v3.1 Cycle Sequencing Kit (Invitrogen, Carlsbad, CA, USA) on 3730XL DNA analyzers (Applied Biosystems, Foster City, CA, USA), according to the manufacturers’ recommendations. DNA sequences for all genes were edited, compiled, assembled, and analyzed by using SeqMan in Lasergene 8 (DNASTAR, Madison, WI, USA). Nucleotide sequences were compared with sequences available in GenBank. Multiple nucleotide sequence alignment and the alignment of the deduced amino acids of all of the gene segments were conducted by using ClustalW in BioEdit 7.09.0 (www.mbio.ncsu.edu/bioedit/bioedit.html). The phylogeny of each gene was investigated within the maximum-likelihood framework by using PhyML under the general time reversible substitution model with gamma rate heterogeneity (30). The robustness of the grouping was assessed by using 100 bootstrapping replicates (30). MegAlign in Lasergene 8 (DNASTAR) was used to determine the percentage of nucleotide sequence similarities, and the NetNGlyc 1.0 Server (31) was used to predict the glycosylation sites in the HA and NA genes. The sequences we obtained (Technical Appendix, Table 1) were deposited into GenBank under accession nos. KC757782–KC758115.

ResultsPrevalence of Avian Influenza in Bangladesh

During November 2008–August 2011, a total of 17,438 samples were collected from live-bird markets, farms, and wild and environmental sources (26). In the markets and farms, samples were predominantly collected from chicken, quail, and ducks. The environmental samples were collected from cages in the markets in which only 1 species was housed (e.g., only chickens or only ducks). The wild bird samples were from a lake that was a feeding station for wild migratory ducks.

The predominant subtypes identified in these samples were low pathogenicity H9N2 (n = 734) and highly pathogenic H5N1 (n = 92) (Table 1). Other influenza subtypes identified, at low percentages and only during 2008–2009, were H1N2, H1N3, H3N6, H3N8, H4N2, and H10N7 (Table 1). None of the wild-bird samples tested and only 35 (0.87%) of 4,045 farm samples had positive test results for avian influenza (Table 1). The subtype H9N2 and H5N1 viruses were isolated from domestic poultry, mainly chickens and quail that were sampled in live-bird or pet markets (Table 1).

The number of subtype H5N1 viruses isolated from field samples increased during the surveillance period (Figure 1, panel B). The numbers of other influenza subtypes identified were low during 2009–2010; during 2010–2011, subtypes H9N2 and H5N1 were most commonly isolated.

Influenza virus surveillance in Bangladesh, November 2008–September 2012, showing the total number of avian influenza viruses isolated per month. A) Low pathogenicity subtype H9N2 viruses; B) highly pathogenic subtype H5N1 viruses.

AIVs and Newcastle disease virus (NDV) are the 2 most common viruses that infect poultry in Bangladesh; therefore, all of the influenza-positive samples that were subtyped and used in this study were screened again for H9N2, H5N1, and NDV to determine whether a single host species was co-infected with AIV and NDV. In a substantial proportion of the isolates screened, we identified mixtures of either subtype H9N2 and NDV or subtypes H9N2 and H5N1 (Table 2). All of the samples containing a mixture of viruses were obtained from apparently healthy birds. The highest incidence of influenza (H5N1) virus isolation was during the winter months (January–March), whereas subtype H9N2 was isolated year-round (Figure 1, panel A)

Antigenic Properties

Antigenic analysis showed that most chicken H9N2 viruses circulating in Bangladesh were antigenically homogenous and similar to the human H9N2 virus (A/BD/0994/11) isolated from Bangladesh (Table 3, Appendix). Antigenic variation was seen in subtype H9N2 viruses isolated from quail; specifically, A/Env/BD/907/2009 (H9N2) (quail) was distinguishable from other contemporary quail isolates when postinfection ferret serum of A/Env/BD/907/2009 (H9N2) (quail) and chicken H9N2 viruses was used (Table 3, Appendix). The chicken H9N2 viruses had moderate cross-reaction to antiserum of A/Env/BD/907/2009 (H9N2) (quail), indicating slight antigenic variation from subtype H9N2 viruses isolated from quail (Table 3, Appendix). The low level of cross-reactivity seen across the quail viruses and between the quail and chicken viruses indicates that the HAs of the subtype H9N2 viruses maintained in quail were antigenically distinct from those of subtype H9N2 viruses maintained in chickens. In addition to the quail isolate, A/Env/BD/12068/2011 (H9N2) (pigeon) reacted with 2 serum samples, the polyclonal serum against A/Env/BD/907/2009 (H9N2) (quail) and that against A/Ck/Pk/(NARC-2434)/06 (H9N2). The HI titers of those reactions were 320 and 160, respectively.

Antigenic analyses of avian influenza A(H9N2) viruses isolated from poultry in Bangladesh during 2008–2011 compared with reference strains*
Viruses and designationsHI titers† obtained with reference ferret hyperimmune serum directed against
G1 clade
Y280 Clade
G1 clade
Bangladesh subtype H9N2 viruses
Qa/HK/G1/97HK/1073/97Ck/HK/G9/97Ck/Bei/1/94Ck/Pk (Swabi)/NARC/2434/06BD/0994/11Env/BD/600/08Env/BD/5473/09Env/BD/907/09
Reference viruses and designation
A/Qa/HK/G1/97160320
A/Ck/HK/G9/9740806403201601608080
A/Ck/Beijing/1/9480
A/HK/1073/97 (human)80320
A/Ck/Pk(Swabi)/NARC–2434/06801603203202,5602,5602,5601,280320
A/BD/0994/2011 (human)801601603201,2802,5601,280640160
A/Env/BD/600/2008 (chicken)801601603201,2802,5601,2801,280160
A/Env/BD/907/2009 (quail)40804040808040160
A/Env/BD/5473/2009 (chicken)
80
160

320
320

2,560

2,560
2,560
1,280
640
Bangladesh isolates
From chickens, n = 27‡
A/Ck/BD/2075/2009801603206402,5601,280640160
A/Ck/BD/5209/2009401603203201,2802,5602,560640320
A/Env/BD/8202/2010 (chicken)801603203201,2801,2801,280640320
A/Env/BD/8465/2010 (chicken)1603206403202,560>5,1202,5602,560320
A/Ck/BD/8725/2010801603203201,2802,5601,280640160
A/Ck/BD/9334/201040801601606401,280640320160
A/Env/BD/10234/2011 (chicken)401603203201,280>5,1202,5601,280320
A/Ck/BD/10450/2011801603203206401,2806401,280160
A/Ck/BD/11309/2011401603203201,2802,5601,2801,280320
A/Ck/BD/13962/2011401601603203206401,28032080
From quail, n = 10‡         
A/Env/BD/5144/ 2009 (quail)4040401601601608040
A/Env/BD /10306/2011 (quail)401604040404080
A/Env/BD/10316/2011 (quail)160320808016080808080
A/Env/BD/12103/ 2011 (quail)40404040404040
A/Env/BD/12116/ 2011 (quail)
A/Env/BD/12119/2011 (quail)80160808080160808080
From other poultry species, n = 7‡         
A/Env/BD/1041/2009 (duck)
A/Dk/BD/1231/2009401601603201,2802,5601,280640320
A/Pi/BD/4303/2009401601601606402,5601,280640160
A/Env/BD/12068/2011 (pigeon)16016040808080404080
A/Env/BD/9306/2010 (parrot)801601603206402,5601,280640320
A/Env/BD/12077/2011 (turkey)401603203201,280>5,1202,560640320

*Boldface indicates homologous serum. HI, hemagglutination inhibition. 
†HI tests were performed according to standard procedures (27).Titers are expressed as the reciprocal of the highest dilution of the last dilution that completely inhibited hemagglutination of 0.5% chicken erythrocytes.
‡Total number of influenza (H9N2) isolates with a similar HI pattern.

These antigenic analyses show that antigenic variation occurs among influenza (H9N2) viruses isolated from different host species. The AIV (H9N2) viruses from Bangladesh cross-reacted well to antiserum against G1 clade H9N2 viruses that were isolated in Bangladesh, Dubai, or Pakistan; however, they showed low reactivity to antiserum against A/Qa/HK/G1/97 (prototype G1), indicating that the viruses from Bangladesh were related to but very distinct from the G1 clade isolates. These results indicate antigenic diversity within the G1 clade.

Phylogenetic Relationships

The phylogenetic relationships of all 8 genes of the 44 representative AIV (H9N2) viruses were analyzed on the basis of their nucleotide sequences and by comparing them with nucleotide sequences of subtype H9N2 viruses belonging to the main lineages from Asia (14). All sequences from Bangladesh except 1 formed a monophyletic cluster. Phylogenetic analyses showed that all the subtype H9N2 viruses isolated from Bangladesh are part of the lineage represented by A/Qa/HK/G1/97 (G1), which is the most dominant lineage worldwide. The HA, NA, nucleoprotein (NP), matrix (M), and polymerase PB2 genes of the viruses from Bangladesh originated from the prototype G1 virus A/Qa/HK/G1/97 (H9N2). The nonstructural (NS) and polymerase PA and PB1 genes were closely related to those in the subtype H7N3 isolate from Pakistan (14).

Phylogenetically, the HA and NA genes of most of the subtype H9N2 viruses from Bangladesh tightly clustered with those of viruses from India and Pakistan (Figure 2; Technical Appendix Figure). These viruses had >95% sequence homology, grouped within the G1 lineage, and shared a common ancestor with A/Qa/Hk/G1/97. One duck isolate, A/Env/BD/1041/09 (H9N2) (duck), clustered with subtype H9N2 viruses from Korea and shared a common ancestor with A/Dk/HK/Y439/97 (Figure 2). This isolate also showed only 85% nt similarity with the other viruses from Bangladesh.

Phylogenetic relationships of hemagglutinin genes of avian influenza (H9N2) viruses (boldface) isolated in Bangladesh. Full-length DNA sequencing, starting from the first codon, was used. The phylogenetic trees were generated by PhyML (30) within the maximum-likelihood framework. Numbers above the branches indicate bootstrap values; only values >60 are shown. Boldface italics indicate prototype subtype H9N2 viruses from the Ck/Bei and G1 clades. Scale bar indicates distance between sequence pairs.

Three internal genes (NP, M, and PB2) of the subtype H9N2 viruses from Bangladesh showed a high level of nucleotide homology with the G1-like lineage. All 3 genes were phylogenetically closely related to viruses from India, Pakistan, or the Middle East (Technical Appendix Figure). The duck isolate A/Env/BD/1041/09 (H9N2) (duck) showed similar characteristics in the HA and NA trees, where it shared a common ancestor with a lineage from Korea. The other 3 internal genes (NS, PB1, and PA) of all the subtype H9N2 viruses from Bangladesh grouped into a single cluster, sharing high sequence homology with A/Ck/Karachi/NARC-100/04 (H7N3) (Figure 3; Technical Appendix Figure). These genes formed a distinct clade adjacent to the G1 and Ck/Bei lineages. Although subtype H5N1 viruses were isolated separately and in a mixture with subtype H9N2 viruses, phylogenetic analysis showed no reassortment between the isolates characterized in this study and the subtype H5N1 strains.

Phylogenetic relationships of nonstructural protein genes of avian influenza (H9N2) viruses (boldface) isolated in Bangladesh. Full-length DNA sequencing, starting from the first codon, was used. The phylogenetic trees were generated by PhyML (30) within the maximum-likelihood framework. Numbers above the branches indicate bootstrap values; only values >60 are shown. Boldface italics indicate prototype subtype H9N2 viruses from the Ck/Bei and G1 clades. Scale bar indicates distance between sequence pairs.

Molecular Characteristics

To identify the possible determinants of transmission of AIV (H9N2) from birds to humans, we aligned the amino acid sequences of all genes of the subtype H9N2 viruses from Bangladesh and compared them with those of representative subtype H9N2 viruses from different clades (Technical Appendix Tables 2, 3). When compared with the prototype G1 virus, the viruses from Bangladesh showed that they have evolved to acquire mammalian host–specific mutations throughout the genome (Technical Appendix Tables 2, 3). Comparing the amino acid sequences of all genes of the subtype H9N2 viruses from Bangladesh with the prototype subtype H9N2 viruses showed that certain amino acid substitutions throughout the viral genome have become fixed as the viruses have evolved (Figure 4). The receptor-binding site (RBS) of the virus HA influences the generation of human viruses from avian precursors (1). Within the RBS, 42 (95.5%) of 44 isolates had leucine (L) at position 226 (H3 numbering), whereas the other 2 viruses had glutamine (Q) at the same position. Of the 2 isolates with 226Q, 1 was a duck virus isolated in early 2009 and 1 was a quail virus isolated in 2011. All viruses with 226L were isolated during 2009–2011 and were found in all of the poultry species sampled. Amino acids at positions 183, 189, 190, and 226 (H3 numbering) are located within the RBS of the HA protein, and the combination of these 4 residues is essential for respiratory droplet transmission of a subtype H9N2 or H3N2 reassortants in ferrets (6). Of the 44 viruses analyzed, 42 had 2 (183H and 226L) of the 4 (183H, 189H, 190E, 226L) substitutions. Variations in the number of potential glycosylation sites within the HA are thought to be associated with the adaptation of duck viruses to poultry (7). Seven potential glycosylation sites (29, 105,141, 298, 305, 492, and 551) were found in the HAs of 43 (97.7%) of the 44 isolates; the N-X-T/S motif (in which X may be any amino acid except proline) and 2 potential glycosylation sites were lost at positions 208 and 218 (Technical Appendix Tables 2, 3). In the NA proteins of all of the analyzed viruses, we found no R292K substitution, which is associated with resistance to the neuraminidase inhibitor oseltamivir.

Host range and pathogenicity determinants in avian influenza (H9N2) viruses isolated from different poultry species in Bangladesh during 2008–2011. Numbers in parentheses indicate number of viruses containing specific amino acid residues of the 44 virus isolates analyzed. Blue indicates the 11 genes that were assessed; red indicates the residues that are critical for influenza pathogenesis, enhanced replication in mammalian hosts, or those that are identical to residues present in human influenza viruses; green indicates unique substitutions in the viruses. HA, hemagglutinin; NA, neuraminidase; M, matrix; NS, nonstructural; NP, nucleoprotein; PA and PB, polymerase genes.

Sequence analysis of the internal gene products also identified mammalian host–specific markers (11,14) that were fixed in the M protein (M1 V15I, 43 isolates; M2 L55F, 44 isolates) and NP (11,14) (R214K, 37 isolates; E372D, 40 isolates; PB1-L13P, 40 isolates) (Figure 4). In the NS protein, all isolates were full length and possessed the PDZ-(PL) C-terminal motif (11,14). The NS1 protein also harbored the mammalian-specific glutamate (E) to lysine (K) substitution at position 227 in 43 of the 44 isolates tested (Figure 4).

Our comparison of the deduced amino acid sequences of the subtype H9N2 viruses from Bangladesh revealed that, in positions that were previously identified as important for host specificity, unique substitutions whose functions are unknown have become fixed. In the RBS of the HA protein, the Q227I substitution was in 35 of 44 isolates. This substitution was previously identified in subtype H9N2 viruses from Pakistan and the Middle East but was not found in any of the reference strains. The other unique substitutions were in the M2 protein (E16V and I28F, 43 isolates) and the PB1-F2 protein (R79L, 40 isolates) (Figure 4). Alignment of the M2 proteins showed that 41 of the 44 isolates had a substitution at position 31 (S31N), suggesting resistance to M2-blocker antiviral drugs such as amantadine.

The amino acid sequence of the HA1/HA2-connecting peptide of HA is a major determinant of pathogenicity in terrestrial poultry, and the pathogenicity of highly pathogenic subtype H5N1 and H7N3 viruses is influenced by the presence of a polybasic cleavage site in the connecting peptide (32). The amino acid sequences of the cleavage site of the HA1/HA2 junction in the subtype H9N2 viruses from Bangladesh possessed 4 cleavage motifs (Figure 4). Of the 44 isolates we tested, 28 had the KSKR/GLF motif (identified in all recent poultry isolates and a 2011 quail isolate); 14 isolates had the KSSR/GLF motif (identified in chickens, ducks, and older quail isolates); and 1 duck and 1 pigeon isolate had the ASDR/GLF and KASR/GLF motifs, respectively. Thus, the altered amino acid sequences in the connecting peptide of viruses we analyzed were distinct from the RSSR/GLF motif found in the connecting peptides of prototype subtype H9N2 viruses. This finding is clear evidence of genetic variation within isolates from different poultry species.

Discussion

We isolated a large number of low pathogenicity avian influenza A(H9N2) viruses and an increasing number of highly pathogenic avian influenza A(H5N1) viruses from apparently healthy birds in Bangladesh. During 2008–2011, subtype H9N2 viruses were isolated from multiple poultry species (predominantly chickens and quail) throughout the year, without seasonal prevalence; subtype H5N1 viruses were predominantly isolated from quail during the winter of 2011. Phylogenetic analyses revealed that the subtype H9N2 viruses are of the G1 clade; antigenically, the G1 clade consists of 2 branches distinguished by host species, but the Bangladesh viruses were divergent from the prototype G1 virus. Our comparison of the deduced amino acid sequence showed that the subtype H9N2 viruses have acquired mammalian host–specific mutations in their surface glycoproteins and internal genes.

We isolated subtype H9N2 and H5N1 viruses individually and in viral mixtures; the latter may result from co-infection of the same host species with multiple influenza subtypes. In our surveillance, we isolated influenza subtypes other these than mostly during 2008–2009, with very few other subtypes isolated after 2009. This finding suggests that subtypes H9N2 and H5N1 predominate in live-bird markets.

Phylogenetic analyses showed that, with 1 exception, the isolates we analyzed were homogenous and shared the common ancestor A/Qa/HK/G1/97 (H9N2). G1 was the only lineage identified. The viruses we analyzed shared close nucleotide homology with those of G1 lineage isolated from birds in Pakistan or India, which suggests that the Bangladesh viruses may have their origin in those countries (or vice versa). We speculate that, during the early 2000s, subtype H9N2 viruses circulating in Pakistan were introduced into Bangladesh through regional poultry or pet trade and then established themselves in domestic poultry. Although subtype H9N2 and H5N1 viruses were co-isolated, phylogenetic analysis did not identify any reassortment between the subtypes. However, the subtype H9N2 viruses still maintained the internal genes that were part of the subtype H5N1 reassortant previously implicated in human infection (7).

Results from our analyses support our hypothesis that the subtype H9N2 isolates from Bangladesh are related to those from Pakistan and have evolved into 2 distinct subpopulations based on host species. The chicken viruses we analyzed were antigenically homogenous and are related to human viruses but distinct from quail viruses. Furthermore, small changes in the quail viruses should be monitored closely because quail were the host species implicated during the emergence of the subtype H9N2 virus that jumped to humans in 1999 (7). In addition, quail were suggested to be the contributing host species of the H9N2/H5N1 reassortant of the subtype H5N1 virus that infected humans in 1997 (4)

Our analyses also shows that subtype H9N2 viruses circulating in Bangladesh show a trend toward accumulating molecular markers that favor interspecies transmission. Distinct mutations throughout the viral genomes have been implicated in the adaptation of viruses to a mammalian host (Table 4, Appendix; and Table 5, Appendix). In almost all of the Bangladesh viruses, the RBS of HA has the Q226L substitution, which has been implicated in human virus–like receptor specificity and is critical for replication and direct transmission of these viruses in ferrets (33,34).

HA, NA, and M gene substitutions that favor replication in mammals in avian influenza A(H9N2) viruses isolated from poultry in Bangladesh during 2008–2011*†
Residues and isolates
HA
NA
M1
M2 drug resistance
HA RBS

HA1/HA2

Glyco site
HB site
15
16
20
28
55
26
31
Residues
H9N2 numbering166191234399218372403
H3 numbering
158
183
226
391
















Isolates
A/Env/BD/907/2009NHLKKSSRNoAWIVSFFLN
A/Env/BD1041/2009SHQKASDRNRTFSWVESILLS
A/Ck/BD/2075/2009NHLKKSSRNoAWIVSFFLS
A/Env/BD/8202/2010NHLKKSSRNoAWIVNFFFN
A/Ck/BD/8996/2010NHLKKSKRNoAWIVSFFLN
A/Ck/BD/9334/2010NHLKKSKRNoAWIVSFFLN
A/Env/BD/10234/2011NHLKKSKRNoAWIVSFFLN
A/Env/BD/10306/2011DHQKKSKRNoAWIVSFFLN
A/Ck/BD/10450/2011NHLKKSKRNoAWIVSFFLN
A/Ck/BD/10897/2011NHLKKSKRNoAWIVSFFLN
A/Ck/BD/11154/2011NHLKKSKRNoAWIVSFFLN
A/Env/BD/11173/2011NHLKKSKRNoAWIVSFFLN
A/Env/BD/11309/2011NHLKKSKRNoAWIVSFFLN
A/Env/BD/11315/2011NHLKKSKRNoAWIVSFFLN
A/Env/BD/12116/2011NHLKKSKRNoNDNDIVSFFLN
A/Env/BD/12119/2011NHLKKSKRNoAWIVSFFLN
A/Env/BD/13916/2011NHLKKSSRNoARIVSFFLN
A/Ck/BD/13962/2011NHLKKSKRNoAWIVSFFLN

*A total of 31 unique substitutions were identified in 18 H9 isolates. Isolates 8202 (M2 adamantine-like drug resistance mutation), 10306 (HA158D) and13916 (NA403R) contain unique substitution not seen in most isolates. Boldface indicates key amino acid substitutions implicated in mammalian transmission. HA, hemagglutinin; NA, neuraminidase; M, matrix; RBS, receptor-binding site; glyco, glycosylation; HB, hemabsorbing; ND, analysis not done because sequences were not available.
†Amino acids: A, alanine; D, aspartic acid; E, glutamic acid; H, histidine; I, isoleucine; K, lysine; L, leucine; N, asparagine; Q, glutamine; R, arginine; S, serine; V, valine.

NS, NP, PA, PB1, and PB2 gene substitutions that favor replication in mammals in avian influenza A(H9N2) viruses isolated from poultry in Bangladesh during 2008–2011*†
IsolateNS
NP
PA
PA-XPB1
PB1-F2
PB2
NS1 (227)PL motif3310921431337257409133363756668738281318
A/Env/BD/907/2009KKSEVVIKFDRSYesPVNNTKSTK
A/Env/BD1041/2009EESEVVIRFERSNDNDNDNDNDNDNDTR
A/Ck/BD/2075/2009KKSEVVIRFDRSPVTNTKLTK
A/Env/BD/8202/2010KKPEVVVKYDRSPVSNTKLIK
A/Ck/BD/8996/2010KKSEVVVKFDRNPVSNTKLIK
A/Ck/BD/9334/2010NDNDIIKFDRSPVNNTKLTK
A/Env/BD/10234/2011KKSEVVVKFDRSPVSNTKLIK
A/Env/BD/10306/2011KKSEVVVKFDRSPVNNIKLIK
A/Ck/BD/10450/2011KKSEVVVKFDRSPVNNTKLTK
A/Ck/BD/10897/2011KKSEVVIKFDRSPVSNTKLTK
A/Ck/BD/11154/2011KKSEVVVKFDRSPVNNTKLTK
A/Env/BD/11173/2011KKSEVVVKFDRSPVNNTKLIK
A/Env/BD/11309/2011KKSEVVVKFDRSPVNNTKLNDND
A/Env/BD/11315/2011KKSEVVVKFDRSPVNNTKLIK
A/Env/BD/12116/2011KKSEVVVKFDNDNDPVNNTKLIK
A/Env/BD/12119/2011KKSEVVVKFDNDNDPVNNTKLIK
A/Env/BD/13916/2011KKSEVVIKFDNDNDPINNTKLTK
A/Ck/BD/13962/2011KKSEVVVKFDNDNDPVNSTRLIK

*A total of 31 unique substitutions were identified in 18 H9 isolates. Isolates 907 (PB1-F2, 82S), 8202 (NP-313Y), 10306 (PB1-F2 68I), 10897 (PB1375S), 13916 (PB1336I), and 13962(PB1-F2 66S/82R) contain unique substitution not seen in most isolates. Boldface indicates key amino acid substitutions implicated in mammalian transmission. NS, nonstructural; NP, nucleoprotein; PA and PB, polymerase genes; ND, analysis not done because sequences were not available.
†Amino acids: D, aspartic acid; E, glutamic acid; I, isoleucine; K, lysine; L, leucine; N, asparagine; R, arginine; S, serine; T, threonine; V, valine.

Although the viruses we isolated have low pathogenicity and do not possess a polybasic cleavage site, they carry a motif similar to cleavage sites seen in the highly pathogenic H5N1 and H7N3 viruses. This result suggests that subtype H9N2 viruses have the potential to acquire a polybasic site and become highly pathogenic. In the HA of nearly all the isolates that we analyzed, 2 potential glycosylation sites were lost; this process was also seen in the adaptation of duck viruses to poultry (32,35). Because the duck isolate A/Env/BD/1041/2009 retained its glycosylation site, we hypothesize that the introduction of the Y439 lineage via ducks was unsuccessful because the G1 lineage was already dominant in the local poultry species.

In most of the viruses we analyzed, we observed residues in the avian–human signature positions becoming fixed (HA 226L and 183H, M1 V15I, M2 L55F, NS E227K, NP R214K, E372D, PB1-L13P); in certain isolates, the mutations were fixed by unique substitutions whose functions are unknown (HA Q227I, M2 E16V, I28F, PB1-F2 R79L). In addition, these viruses have acquired substitutions similar to those seen in highly pathogenic subtype H5N1 and H7N3 and pandemic subtype H1N1 and H3N2 viruses (36). Noticeably, in the NS gene’s PDZ motif, most subtype H9N2 isolates carried the KSEV motif seen in the 1918 pandemic subtype H1N1 virus (11,37). In the PA genes of 2 subtype H9N2 viruses, we identified R57Q and S409N substitutions, both of which were in the 1968 pandemic subtype H3N2 virus. These observations confirm that the Bangladesh subtype H9N2 viruses have accumulated molecular markers that influence host specificity and pathogenesis.

In conclusion, our results show that the Bangladesh subtype H9N2 viruses are genetically similar to, but distinct from, the A/Qa/HK/G1/97 isolate, which has been previously implicated in human infection (7). However, the Bangladesh viruses have accumulated molecular characteristics required for infecting humans and are antigenically similar to the human subtype H9N2 virus isolated from a patient in Bangladesh. These viruses have a tendency to reassort with highly pathogenic subtype H5N1 (7) and H7 viruses (38,39), which is particularly concerning in light of the emergence of avian influenza A(H7N9) virus in China that contains 6 gene segments from H9N2 influenza viruses and causes lethal infection in humans (40). The emergence of this reassortant virus emphasizes the potential for, and the danger of, transmission of subtype H9N2 viruses to humans.

<supplementary-material content-type="loca-data" id="SD1"><caption><title>Technical Appendix

Genes of influenza A(H9N2) viruses from Bangladesh that were not sequenced or only partially sequenced, amino acid substitutions in hemagglutinin and neuraminidase gene products and in internal genes that determine host specificity of the viruses, and phylogenetic relations of key internal genes of the virus isolates.

Suggested citation for this article: Shanmuganatham K, Feeroz MM, Jones-Engel L, Smith GJD, Fourment M, Walker D, et al. Antigenic and molecular characterization of avian influenza A(H9N2) viruses, Bangladesh. Emerg Infect Dis [Internet]. 2013 Sep [date cited]. http://dx.doi.org/10.3201/eid1909.130336

Acknowledgments

We thank Jerry Parker and Richard Elia for maintaining the influenza database at St. Jude Children’s Research Hospital, James Knowles for providing administrative assistance, and Angela McArthur for scientific editing. We also thank the Centers for Disease Control and Prevention for providing the human influenza A(H9N2) virus from Bangladesh.

This work was supported by the National Institute of Allergy and Infectious Diseases, National Institutes of Health, Department of Health and Human Services (contract nos. N01-AI30071 and HHSN272200900007C); and the American Lebanese Syrian Associated Charities.

Dr Shanmugantham is a postdoctoral fellow working in Dr Webster’s laboratory at St. Jude Children’s Research Hospital. His primary research interest is the epidemiology, control, and prevention of zoonotic diseases.

ReferencesButt KM, Smith GJ, Chen H, Zhang LJ, Leung YH, Xu KM, Human infection with an avian H9N2 influenza A virus in Hong Kong in 2003. J Clin Microbiol. 2005;43:57607 and. 10.1128/JCM.43.11.5760-5767.200516272514Lin YP, Shaw M, Gregory V, Cameron K, Lim W, Klimov A, Avian-to-human transmission of H9N2 subtype influenza A viruses: relationship between H9N2 and H5N1 human isolates. Proc Natl Acad Sci U S A. 2000;97:96548 and. 10.1073/pnas.16027069710920197Peiris JS, Guan Y, Markwell D, Ghose P, Webster RG, Shortridge KF. Cocirculation of avian H9N2 and contemporary “human” H3N2 influenza A viruses in pigs in southeastern China: potential for genetic reassortment? J Virol. 2001;75:967986 and. 10.1128/JVI.75.20.9679-9686.200111559800Peiris M, Yuen KY, Leung CW, Chan KH, Ip PLS, Lai RWM, Human infection with influenza H9N2. Lancet. 1999;354:9167 and. 10.1016/S0140-6736(99)03311-510489954Apurba Chakraborty Outbreak of mild respiratory disease caused by H5N1 and H9N2 infections among young children in Dhaka, Bangladesh. Health and Science Bulletin (ICDDR,B). 2011 [cited 2012 Dec 18]. http://www.icddrb.org/what-we-do/publications/cat_view/52-publications/10042-icddrb-periodicals/10048-health-and-science-bulletin-bangla-and-english/11089-vol-9-no-2-english-2011-/11092-outbreak-of-mild-respiratory-disease-caused-by-h5n1-and-h9n2-infections-among-young-children-in-dhaka-bangladesh-2011Guan Y, Shortridge KF, Krauss S, Chin PS, Dyrting KC, Ellis TM, H9N2 influenza viruses possessing H5N1-like internal genomes continue to circulate in poultry in southeastern China. J Virol. 2000;74:937280 and. 10.1128/JVI.74.20.9372-9380.200011000205Guan Y, Shortridge KF, Krauss S, Webster RG. Molecular characterization of H9N2 influenza viruses: were they the donors of the “internal” genes of H5N1 viruses in Hong Kong? Proc Natl Acad Sci U S A. 1999;96:93637 and. 10.1073/pnas.96.16.936310430948Aamir UB, Wernery U, Ilyushina N, Webster RG. Characterization of avian H9N2 influenza viruses from United Arab Emirates 2000 to 2003. Virology. 2007;361:4555 and. 10.1016/j.virol.2006.10.03717157891Alexander DJ. Report on avian influenza in the Eastern Hemisphere during 1997–2002. Avian Dis. 2003;47(Suppl):7927 and. 10.1637/0005-2086-47.s3.79214575066Cameron KR, Gregory V, Banks J, Brown IH, Alexander DJ, Hay AJ, H9N2 subtype influenza A viruses in poultry in Pakistan are closely related to the H9N2 viruses responsible for human infection in Hong Kong. Virology. 2000;278:3641 and. 10.1006/viro.2000.058511112478Fusaro A, Monne I, Salviato A, Valastro V, Schivo A, Amarin NM, Phylogeography and evolutionary history of reassortant H9N2 viruses with potential human health implications. J Virol. 2011;85:841321 10.1128/JVI.00219-1121680519Ge FF, Zhou JP, Liu J, Wang J, Zhang WY, Sheng LP, Genetic evolution of H9 subtype influenza viruses from live poultry markets in Shanghai, China. J Clin Microbiol. 2009;47:3294300 and. 10.1128/JCM.00355-0919656985Golender N, Panshin A, Banet-Noach C, Nagar S, Pokamunski S, Pirak M, Genetic characterization of avian influenza viruses isolated in Israel during 2000–2006. Virus Genes. 2008;37:28997 and. 10.1007/s11262-008-0272-718712589Iqbal M, Yaqub T, Reddy K, McCauley JW. Novel genotypes of H9N2 influenza A viruses isolated from poultry in Pakistan containing NS genes similar to highly pathogenic H7N3 and H5N1 viruses. PLoS ONE. 2009;4:e5788 and. 10.1371/journal.pone.000578819517011Lee YJ, Shin JY, Song MS, Lee YM, Choi JG, Lee EK, Continuing evolution of H9 influenza viruses in Korean poultry. Virology. 2007;359:31323 and. 10.1016/j.virol.2006.09.02517056087Li C, Yu K, Tian G, Yu D, Liu L, Jing B, Evolution of H9N2 influenza viruses from domestic poultry in Mainland China. Virology. 2005;340:7083 and. 10.1016/j.virol.2005.06.02516026813Toroghi R, Momayez R. Biological and molecular characterization of avian influenza virus (H9N2) isolates from Iran. Acta Virol. 2006;50:1638 .17131934Xu KM, Li KS, Smith GJ, Li JW, Tai H, Zhang JX, Evolution and molecular epidemiology of H9N2 influenza A viruses from quail in southern China, 2000 to 2005. J Virol. 2007;81:263545 and. 10.1128/JVI.02316-0617192315Guo YJ, Krauss S, Senne DA, Mo IP, Lo KS, Xiong XP, Characterization of the pathogenicity of members of the newly established H9N2 influenza virus lineages in Asia. Virology. 2000;267:27988 and. 10.1006/viro.1999.011510662623Li KS, Xu KM, Peiris JS, Poon LL, Yu KZ, Yuen KY, Characterization of H9 subtype influenza viruses from the ducks of southern China: a candidate for the next influenza pandemic in humans? J Virol. 2003;77:698894 and. 10.1128/JVI.77.12.6988-6994.200312768017Ahmed SS, Ersboll AK, Biswas PK, Christensen JP, Hannan AS, Toft N. Ecological determinants of highly pathogenic avian influenza (H5N1) outbreaks in Bangladesh. PLoS ONE. 2012;7:e33938 and. 10.1371/journal.pone.003393822470496International Centre for Diarrhoeal Disease Research Bangladesh. Backyard poultry raising practices in Bangladesh. Health and Science Bulletin (ICDDR,B). 2012 [cited 2012 Nov 26]. http://www.icddrb.org/what-we-do/publications/cat_view/52-publications/10042-icddrb-periodicals/10048-health-and-science-bulletin-bangla-and-english/14138-vol-10-no-3-english-2012/14144-backyard-poultry-raising-practices-in-bangladesh-implications-for-risk-of-avian-influenza-infection-in-humansWorld Health Organization Cumulative number of confirmed human cases of avian influenza A(H5N1) [cited 2012 Nov 26]. http://www.who.int/influenza/human_animal_ interface/H5N1_cumulative_table_archives/en/index.htmlExternalMondal SP, Tardif-Douglin D, Ryan-Silva R, Magnani R. Controlling highly pathogenic avian influenza, Bangladesh. Emerg Infect Dis. 2012;18:20835 and. 10.3201/eid1812.12063523171830Dolberg F. Poultry sector country review: Bangladesh. Food and Agricultural Organization of the United Nations. 2009 [cited 2012 Nov 26]. ftp://ftp.fao.org/docrep/fao/011/ai319e/ai319e00.pdfNegovetich NJ, Feeroz MM, Jones-Engel L, Walker D, Alam SM, Hasan K, Live bird markets of Bangladesh: H9N2 viruses and the near absence of highly pathogenic H5N1 influenza. PLoS ONE. 2011;6:e19311 and. 10.1371/journal.pone.001931121541296Palmer DFDW, Coleman MT, Schild GC. Advanced laboratory techniques for influenza diagnosis. Immunology series no. 6. Atlanta: Centers for Disease Control, US Department of Health, Education and Welfare; 1975Guan Y, Peiris JS, Lipatov AS, Ellis TM, Dyrting KC, Krauss S, Emergence of multiple genotypes of H5N1 avian influenza viruses in Hong Kong SAR. Proc Natl Acad Sci U S A. 2002;99:89505 and. 10.1073/pnas.13226899912077307Hoffmann E, Stech J, Guan Y, Webster RG, Perez DR. Universal primer set for the full-length amplification of all influenza A viruses. Arch Virol. 2001;146:227589 and. 10.1007/s00705017000211811679Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010;59:30721 and. 10.1093/sysbio/syq01020525638Gupta R, Jung E, Brunak S. Prediction of N-glycosylation sites in human proteins [cited 2012 Nov 26]. http://www.cbs.dtu.dk/services/NetNGlyc/Baigent SJ, McCauley JW. Influenza type A in humans, mammals and birds: determinants of virus virulence, host-range and interspecies transmission. Bioessays. 2003;25:65771 and. 10.1002/bies.1030312815721Wan H, Perez DR. Amino acid 226 in the hemagglutinin of H9N2 influenza viruses determines cell tropism and replication in human airway epithelial cells. J Virol. 2007;81:518191 and. 10.1128/JVI.02827-0617344280Wan H, Sorrell EM, Song H, Hossain MJ, Ramirez-Nieto G, Monne I, Replication and transmission of H9N2 influenza viruses in ferrets: evaluation of pandemic potential. PLoS ONE. 2008;3:e2923 and. 10.1371/journal.pone.000292318698430Matrosovich M, Zhou N, Kawaoka Y, Webster R. The surface glycoproteins of H5 influenza viruses isolated from humans, chickens, and wild aquatic birds have distinguishable properties. J Virol. 1999;73:114655 .9882316Chen GW, Chang SC, Mok CK, Lo YL, Kung YN, Huang JH, Genomic signatures of human versus avian influenza A viruses. Emerg Infect Dis. 2006;12:135360 and. 10.3201/eid1209.06027617073083Jackson D, Hossain MJ, Hickman D, Perez DR, Lamb RA. A new influenza virus virulence determinant: the NS1 protein four C-terminal residues modulate pathogenicity. Proc Natl Acad Sci U S A. 2008;105:43816 and. 10.1073/pnas.080048210518334632Abbas MA, Spackman E, Swayne DE, Ahmed Z, Sarmento L, Siddique N, Sequence and phylogenetic analysis of H7N3 avian influenza viruses isolated from poultry in Pakistan 1995–2004. Virol J. 2010;7:137 and. 10.1186/1743-422X-7-13720576101Kageyama T, Fujisaki S, Takashita E, Xu H, Yamada S, Uchida Y, Genetic analysis of novel avian A(H7N9) influenza viruses isolated from patients in China, February to April 2013. Euro Surveill. 2013;18:2045323594575Gao R, Cao B, Hu Y, Feng Z, Wang D, Hu W, Human infection with a novel avian-origin influenza A (H7N9) virus. N Engl J Med. 2013;368:188897 and. 10.1056/NEJMoa130445923577628