MEEGID-02813; No of Pages 11 Infection, Genetics and Evolution 44 (2016) xxx–xxx

Contents lists available at ScienceDirect

Infection, Genetics and Evolution journal homepage: www.elsevier.com/locate/meegid

Research paper

Prevalence and diversity of H9N2 avian influenza in chickens of Northern Vietnam, 2014 Duong Mai Thuy a, Thomas P. Peacock b,c, Vu Thi Ngoc Bich d, Thomas Fabrizio e, Dang Nguyen Hoang f, Nguyen Dang Tho g, Nguyen Thi Diep f, Minh Nguyen d, Le Nguyen Minh Hoa d, Hau Thi Thu Trang d, Marc Choisy d,g, Ken Inui h, Scott Newman h, Nguyen vu Trung i, Rogier van Doorn d,j, Thanh Long To a, Munir Iqbal b, Juliet E. Bryant d,j,⁎ a

National Center for Veterinary Diagnostics, Department of Animal Health, Hanoi, Vietnam Avian Viral Diseases programme, The Pirbright Institute, Woking, UK c St Mary's Campus, Imperial College London, London, UK d Oxford University Clinical Research Unit and Wellcome Trust Major Overseas Programme, Hanoi, Vietnam e St Jude's Center for Excellence in Influenza Research and Surveillance, Memphis, TN, USA f Division of Epidemiology, Department of Animal Health, Hanoi, Vietnam g MIVEGEC (UM1-UM2-CNRS 5290-IRD 224), Centre de Recherche IRD, Montpellier, France h Food and Agriculture Organization of the United Nations, Hanoi, Vietnam i National Hospital Tropical Diseases, Hanoi, Vietnam j Center for Tropical Medicine, Nuffield Department of Clinical Medicine, University of Oxford, Oxford, UK b

a r t i c l e

i n f o

Article history: Received 18 May 2016 Received in revised form 17 June 2016 Accepted 19 June 2016 Available online 20 June 2016 Keywords: Avian influenza Vietnam H9N2 Poultry Chicken H5 Zoonotic

a b s t r a c t Despite their classification as low pathogenicity avian influenza viruses (LPAIV), A/H9N2 viruses cause significant losses in poultry in many countries throughout Asia, the Middle East and North Africa. To date, poultry surveillance in Vietnam has focused on detection of influenza H5 viruses, and there is limited understanding of influenza H9 epidemiology and transmission dynamics. We determined prevalence and diversity of influenza A viruses in chickens from live bird markets (LBM) of 7 northern Vietnamese provinces, using pooled oropharyngeal swabs collected from October to December 2014. Screening by real time RT-PCR revealed 1207/4900 (24.6%) of pooled swabs to be influenza A virus positive; overall prevalence estimates after accounting for pooling (5 swabs/pools) were 5.8% (CI 5.4–6.0). Subtyping was performed on 468 pooled swabs with M gene Ct b 26. No influenza H7 was detected; 422 (90.1%) were H9 positive; and 22 (4.7%) were H5 positive. There was no evidence was of interaction between H9 and H5 virus detection rates. We sequenced 17 whole genomes of A/H9N2, 2 of A/H5N6, and 11 partial genomes. All H9N2 viruses had internal genes that clustered with genotype 57 and were closely related to Chinese human isolates of A/H7N9 and A/H10N8. Using a nucleotide divergence cutoff of 98%, we identified 9 distinct H9 genotypes. Phylogenetic analysis suggested multiple introductions of H9 viruses to northern Vietnam rather than in-situ transmission. Further investigations of H9 prevalence and diversity in other regions of Vietnam are warranted to assess H9 endemicity elsewhere in the country. © 2016 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).

1. Introduction Vietnam has a well-developed system for the active surveillance of highly pathogenic avian influenza (HPAI) strains within live bird markets (LBMs). From 2010 to 2015, surveillance focused primarily on

Abbreviations: LPAI, low pathogenicity avian influenza; HPAI, highly pathogenic avian influenza; FAO, Food and Agriculture Organization of the United Nations; LBM, live bird market. ⁎ Corresponding author at: Oxford University Clinical research Unit and Wellcome Trust Major Overseas Programme, Hanoi, Vietnam. E-mail address: [email protected] (J.E. Bryant).

monitoring genetic diversity and ongoing evolution of HPAI H5 subtype viruses within duck populations (Nguyen et al., 2014). Starting in 2013, additional surveillance activity was initiated in provinces along the Chinese border to promote early detection of H7N9 viruses, and was focused largely on collections from chickens and environmental samples. To date, this second surveillance activity has screened N33,480 pooled chicken oropharyngeal swabs and environmental samples, however no A/H7N9 viruses have been identified in Vietnam (FAO, 2014). Low pathogenicity avian influenza (LPAI) H9N2 is the most prevalent influenza subtype circulating endemically in chickens worldwide. H9N2 viruses are endemic in poultry across Asia, the Middle East and

http://dx.doi.org/10.1016/j.meegid.2016.06.038 1567-1348/© 2016 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).

Please cite this article as: Thuy, D.M., et al., Prevalence and diversity of H9N2 avian influenza in chickens of Northern Vietnam, 2014, Infection, Genetics and Evolution (2016), http://dx.doi.org/10.1016/j.meegid.2016.06.038

2

D.M. Thuy et al. / Infection, Genetics and Evolution 44 (2016) xxx–xxx

Northern Africa (Fusaro et al., 2011), and despite their status as LPAI, poultry outbreaks of H9N2 are associated with significant economic losses, largely due to reduced egg production, reduced feed conversion efficiencies, and highly lethal bacterial or viral co-infections (Li et al., 2005; Parvin et al., 2014). Epidemiological evidence indicates that the chicken adapted H9N2 genotype 57 (G57), also known as genotype S, contains an internal gene cassette (polymerase genes, nucleoprotein, matrix and non-structural genes) that confers increased host range and zoonotic transmission potential to non-H9N2 influenza strains, with recent examples including the zoonotic H7N9 and H10N8 outbreaks in China (H. Chen et al., 2014; Gu et al., 2014; Pu et al., 2014). Additionally H9N2 viruses pose a zoonotic risk in their own right, having caused sporadic human infections in China, Hong Kong, Bangladesh and Egypt (ICDDRB, 2011; Butt et al., 2005; Huang et al., 2015a, 2015b; WHO, 2015a, 2015b). It is suggested that the prevalence and variation of H9N2 influenza virus in farmed poultry could provide an important early warning of the emergence of novel reassortants with pandemic potential (Liu et al., 2014a, 2014b). Despite increasingly widespread use of H9 vaccines in neighboring China (Pu et al., 2014; Wei et al., 2016; P. Zhang et al., 2008), the poultry sector in Vietnam has not yet adopted any vaccination against H9N2 viruses, largely because poultry production remains dominated by smallholder farms, with little vertical integration and relatively low investment in veterinary vaccines. The only licensed avian influenza vaccines for poultry in Vietnam are inactivated virus vaccines against the H5 subtype (Nguyen et al., 2014). Research and surveillance activities in Vietnam from 2003 to 2013 previously identified a number of LPAI subtypes circulating in duck populations (Hotta et al., 2012; Jadhao et al., 2009; Kim et al., 2013; D. C. Nguyen et al., 2005; Nomura et al., 2012; Masatoshi Okamatsu et al., 2013; Takakuwa et al., 2013), including all 3 major antigenically divergent lineages of LPAI H9N2 (i.e. the Chinese BJ94/Y280-like lineage, the Middle Eastern G1-like lineage and the predominantly Korean Y439-like lineage) (Jadhao et al., 2009; Kim et al., 2013; Nomura et al., 2012). In contrast, LPAI transmission within chickens has not been well studied. Here we determine the prevalence and diversity of LPAI in chickens from northern Vietnam, and evaluate the feasibility of direct sequencing of AIV from pooled surveillance swabs.

2. Methods 2.1. Sample collections, molecular screening, and virus isolations Samples used in this study were collected at 10 time points from October to December 2014, representing one discrete period of surveillance activity in the FAO-supported programme for early detection of H7N9 viruses. This surveillance activity focused on sampling LBM in northern provinces believed to play a key role in cross-border poultry trade with China. The sampling frame was based on a stratified sampling design, and entailed collections in 7 northern provinces, with a primary focus on market trade routes along the Vietnam/China border. Five district markets were sampled at weekly intervals for approximately 10 weeks between October 30 and December 31 of 2014. At each sampling timepoint, oropharyngeal (OP) swabs of 35 chickens were sampled individually from 7 different vendors per market, and aggregated into pools (5 swabs per pool). The pooled OP swabs were stored in 2 mL viral transport medium (PBS pH 7.2 supplemented with antibiotics, gentamicin, and buffering agents), transported by cold chain to the National Center for Veterinary Diagnostics (NCVD) and stored at 80 °C prior to further processing. RNA of pooled swab samples were extracted using a tacoTM Nucleic Acid Automatic Extraction System (GeneReach Biotechnology Corp., Taiwan). Detection of influenza A and subtype-specific monoplex real time RTPCR were performed using standardized protocols (Supplementary T2). Allantoic fluid from inoculated embryonated chicken eggs were clarified by centrifugation for

10 min at 4 °C, and presence of virus was determined by hemagglutination assay as previously described (WHO, 2011). 2.2. Next generation sequencing (NGS) We used NGS protocols for influenza sequencing adopted from St Jude's Center for Excellence in influenza Research & Surveillance (SJCEIRS). Briefly, fresh nucleic acid extracts were prepared and samples were pre-amplified to enrich for influenza genomes using MBTUni12/ 13 primers as previously described (Zaraket et al., 2015; Zhou et al., 2009). DNA libraries were prepared using Nextera XT DNA Library Prep Kits (Illumina) with 96 dual-index barcodes according to manufacturer's instructions. Pooled libraries were sequenced on a MiSeq Illumina, using the v2 reagents and 300 cycles. Sequence reads were filtered for quality scores, and the high quality reads aligned to a custom reference library (kindly provided by SJCEIRS) and by de novo assembly. Genomic sequences and variant frequencies were generated in CLC genomics workbench v. 8.5.1 (http://www.clcbio.com/). 2.3. Phylogenetic and genotypic analysis Data used for phylogenetic analyses comprised the novel whole genome sequences generated by this study, as well as additional publicly available sequences of closely related viruses identified by Blastn analysis of NCBI and GISAID databases. Full datasets for each gene segment were aligned using the MUSCLE algorithm in Mega software version 6.06 (www.megasoftware.net). Alignments were trimmed to yield the maximal sequence coverage per segment as follows: 260-2289 for PB2; 78-2242 for PB1; 232-1981 for PA; 160-1599 for HA; 49-1507 for NP; 24-1386 for NA; 26-996 for MP; and 56-859 for NS. Phylogenetic trees were created using the neighbor-joining algorithm and the Tajima-Nei model, with 1000 bootstrap replications. Bootstrap values of b 70 were excluded from trees. Genotyping was performed using 2% nucleotide difference cutoff to indicate genetically divergent segments. 2.4. Statistical analysis Prevalence of influenza type A, H5 and H9 virus infection was estimated from the pooled samples at the market and province levels, using a maximum likelihood modeling approach described in Supplementary Appendix A. We used an influenza A Matrix gene Ct threshold of 35 as the standard cut-off for positivity. For subtype specific prevalence of H5 and H9, the threshold for positivity was set at Ct ≤38. Evidence for interaction between H9 and H5 from pooled samples was tested by a likelihood ratio test as described in the Supplementary Appendix B. 2.5. Nucleotide accession numbers The complete and partial genome sequences of the viruses sequenced in this study were submitted to the GenBank database (Table S2). 3. Results 3.1. Prevalence of influenza A in market chickens A total of 1207 of 2450 pooled swabs (49%) screened positive for influenza A by matrix gene real time RT-PCR. Using a maximum likelihood modeling approach that accounts for pooling (Suppl Appendix A), this detection level corresponded to an overall influenza A prevalence of 5.45% (95% Confidence Interval [CI] 5.4–6.0%). Fig. 1(a) presents estimated prevalence of influenza A per weekly sampling round (x-axis), depicted by province (rows) and by market (individual graphs). None of the 1207 influenza positive pools were found positive for H7 subtype. Further molecular subtyping to detect H5 and H9 viruses was conducted

Please cite this article as: Thuy, D.M., et al., Prevalence and diversity of H9N2 avian influenza in chickens of Northern Vietnam, 2014, Infection, Genetics and Evolution (2016), http://dx.doi.org/10.1016/j.meegid.2016.06.038

D.M. Thuy et al. / Infection, Genetics and Evolution 44 (2016) xxx–xxx

3

Fig. 1. (a) Influenza A prevalence estimates in 5 district markets each for 7 northern provinces from Oct.–Dec. 2014; (b) map showing provinces surveyed by the H7F programme. Yellow = Lao Cai province; orange = Ha Giang province; red = Cao Bang province; green = Ha Noi province; purple = Bac Giang province.

Table 1 H5 and H9 subtyping by RTPCR of H7F surveillance swabs, Oct.–Dec. 2014. AIV+ AIV+ swabs screened (%)⁎ Bac Giang Bac Ninh Cao Bang Ha Giang Ha Noi

230

97 (20.7)

220

82 (17.5)

254

93 (19.9)

167

68 (14.5)

85

34 (7.3)

Lao Cai

246

94 (20.1)

Hưng Yên Total

5 1207

0 468 (100)

H5+ (%)

H9+ (%)

Dual H5+ H9+

H5+ only

H5–H9−

2 (0.4) 0

1

1

34

0

0

2

1

0

1

3

3

5

0

0

0

13 (2.8) 0

62 (63.9%) 80 (97.6%) 92 (98.9%) 60 (88.2%) 34 (100%) 94 (100%) 0

13

0

0

0

0

0

22 (4.7)

422 (90.1)

18 (3.8)

4 42 (8.9) (0.85)

1 (0.2) 6 (1.3) 0

⁎ Only swabs with AIV matrix gene RTPCR Ct b 26 were selected for screening.

on all pooled swabs with an influenza A matrix gene Ct b 26 (n = 468) (Table 1). Detection rates for H9 and H5 subtype varied dramatically between provinces and individual markets, with H9 accounting for 64– 100% of influenza A positive pools; H5 was detected in 4.7% (22/468) of positive pools; 3.8% (18/468) pools screened positive for both H5 and H9 subtypes; and 8.9% (42/468) were classified as ‘subtype unknown’. The maximum likelihood prevalence (accounting for pooling) for H9 and H5 was 3.7% (95%CI 3.3–4.1%) and 1.09% (95%CI 0.7–1.6%), respectively. We found no evidence of significant interaction between H5 and H9 when using a maximum likelihood modeling approach (Suppl Appendix B). Virus isolation in embryonated eggs was attempted on selected representative H9 positive pools with Ct b 24 (n = 25), and all H5 positive pools and pools the were ‘unknown subtype’, regardless of Ct value (n = 22 and 42, respectively). We isolated 25 H9 viruses, however none of the H5 positive pools or the pools of ‘unknown subtype’ yielded a virus isolate.

Please cite this article as: Thuy, D.M., et al., Prevalence and diversity of H9N2 avian influenza in chickens of Northern Vietnam, 2014, Infection, Genetics and Evolution (2016), http://dx.doi.org/10.1016/j.meegid.2016.06.038

4

D.M. Thuy et al. / Infection, Genetics and Evolution 44 (2016) xxx–xxx

3.2. Phylogenetic analysis and genotyping We utilized next generation sequencing using the Illumina Nextera XT Sample Preparation Kit and the MiSeq platform as previously described (Zaraket et al., 2015) using both direct sequencing of pooled surveillance swabs, as well as egg-grown virus isolates. Amplification enrichment for influenza genomes was performed on a total of 96 samples (19 egg isolates, 77 pooled surveillance swabs), from which 55

were selected for Nextera XT library preparations, and 50 were included in the run. Next generation sequencing identified 17 H9N2 complete genomes, 2 H5N6 complete genomes, 11 partial genomes. The HPAI H5 genomes detected from LBM will be presented within a separate publication and are only briefly discussed here. To determine relationships among the novel A/H9N2 sequences, we performed Blastn on each gene segment to identify closely related viruses within public datasets, and included a number of additional references. All H9 HA genes fell

Fig. 2. Neighbor joining phylogenetic trees of Vietnamese H9N2 viruses (n = 17) and reference strains (n = 83 strains). Red = virus sequences generated in this study; Green = genotype 57 viruses as defined by Pu et al., PNAS 2015; Blue = H9 sequences from non-H9N2 subtype viruses; Pink = recent genotype 57 human H9N2 isolates from China; Black = additional representative H9N2 viruses available from Genbank. Bootstrap values below 70% excluded.

Please cite this article as: Thuy, D.M., et al., Prevalence and diversity of H9N2 avian influenza in chickens of Northern Vietnam, 2014, Infection, Genetics and Evolution (2016), http://dx.doi.org/10.1016/j.meegid.2016.06.038

D.M. Thuy et al. / Infection, Genetics and Evolution 44 (2016) xxx–xxx

5

into the BJ94/Y280-like lineage, the predominant H9 lineage found in China. Segment 6 (NA) exhibited the most diversity with 15.7% nucleotide differences between strains, and fell into 2 distinct clusters (Fig. 2). In contrast, the other segments showed reduced diversity, with maximum pairwise nucleotide differences of 6.0%, 4.0%, 2.9%, 5.0%, 4.8%, 1.6% and 6.0% for segments 1–5, 7 and 8 respectively. Each of the 17 whole genomes of H9N2 possessed the G57-like internal gene cassette, similar to zoonotic H7N9 and H10N8 viruses in mainland China. For example, the genes most closely related to the PB2 and PB1 segments of H7F-14-BN4-437 belonged to the human H10N8 and H7N9 isolates, A/Jiangxi/IPB13/2013(H10N8) and A/Guangdong/04/ 2013(H7N9) with percentage nucleotide homologies of 99.6% and 99.5% respectively. Additionally H7F-14-CB4-2, H7F-14-CB4-29, H7F14-CB4-31 and H7F-14-CB4-34 possessed PB1 and NS genes closely related to human and avian H7N9 isolates from 2013 with homologies of over 99.9%, and their PA, NP and MP genes were related to H10N8 viruses with homologies N99.3%. Using an N98% nucleotide difference cutoff for each segment, we determined patterns of clustering at the sub-genotype level, and identified 9 distinct sub-genotypes, all of which appeared to be reassortants between closely related H9N2 viruses (Fig. 3). (For reference, the comparable divergence level used for clade designations in H5 viruses is 1.5%) (WHO/OIE/FAO, 2012). Previous investigations of nucleotide substitution rates of H9N2 viruses have suggested that 2% divergence signifies approximately 4.5 to 9 years of evolution (Fusaro et al., 2011). Of the 18 surveillance swabs initially identified by RT-PCR as containing both H5 and H9, only 8 were subjected to MiSeq (the other 10 samples were excluded during quality control steps), and 1 yielded sequence results indicating presence of both subtypes (H7F-CB4-35). Additionally, there was one sample that tested positive for H5 but not H9 on initial RTPCR screening (H7F-HG4-604), however it yielded H9HA, H5HA, H9-like PB2, PA, NP, NA, MP and NS genes, and an H5-like PB1. The status of these samples as intersubtypic reassortants cannot be established as they represented pooled swabs from multiple birds. Of the samples designated as ‘influenza positive, unknown subtype’ by RT-PCR, 3 of 9 samples sequenced by MiSeq yielded partial sequence data of H9N2, and the remaining did not yield sequences of sufficient coverage for analysis. 3.3. Molecular characteristics of the novel H9N2 genomes 3.3.1. Hemagglutinin (HA) The 19 H9 sequences generated retained the classical BJ94/Y280-like ‘low-pathogenic’ dibasic RSSR cleavage motif that has previously shown to allow cleavage by epithelial proteases such as matriptase and is thought to facilitate systemic infections in poultry (Baron et al., 2013). The receptor binding site of all sequences had leucine at position 216 (226 in H3 numbering), a residue associated with increased binding to human-like α2,6 linked sialic acids in multiple subtypes including H9N2, and consequently associated with mammalian adaptation and transmission (Sorrell et al., 2009). This is a change from previously detected Vietnamese H9N2 viruses which tended to have a glutamine at this position which preferentially binds to the ‘avian’ α2.3 linked sialic acid (Jadhao et al., 2009). Position 217 (227, H3 numbering) has been implicated in both receptor binding and mammalian adaptation in a number of subtypes, however the significance of methionine at this position (as observed in the H9 Vietnamese sequences here) remains unknown (Sang et al., 2015; Stevens et al., 2006). Additionally position 180 (190 H3 numbering) is implicated in receptor binding specificity in the closely related H1 hemagglutinins (Matrosovich et al., 2001; Matrosovich et al., 2000) the Vietnamese viruses have either alanine or valine, which are very common in recent Chinese H9s, or threonine, which is more rarely found, at this position. The HA sequences also contained a number of previously described markers with the potential to confer mammalian transmission potential, as well as the chicken airborne transmission marker K363 (Rudneva et al., 2005; Kaverin et al.,

Fig. 2 (continued).

2005; Sun et al., 2013; Zhong et al., 2014a, 2014b) (Table S1). All viruses contained conserved N-linked glycosylation sequence (N-X-S/T, where X =/=P) at positions 11–13, 123–125, 280–282, 287–289 and 295– 297, and more than half of the sequences had an additional glycosylation sequence at position 200–202 (Table S1). Finally, by analyzing known H9N2 antigenic residues from the literature (Kaverin et al., 2004; Okamatsu et al., 2008; Peacock et al., 2016; Ping et al., 2008; Wan et al., 2014; Zhu et al., 2015), we predict there are at least two antigenically variant sub-groups circulating in northern Vietnam: those antigenically similar to H7F-14-CB4-2 and those similar to HF7-LC4-26 which vary from each other at 11 antigenically important residues (Fig. 4).

Please cite this article as: Thuy, D.M., et al., Prevalence and diversity of H9N2 avian influenza in chickens of Northern Vietnam, 2014, Infection, Genetics and Evolution (2016), http://dx.doi.org/10.1016/j.meegid.2016.06.038

6

D.M. Thuy et al. / Infection, Genetics and Evolution 44 (2016) xxx–xxx

Fig. 2 (continued).

N1 of HPAI H5N1 viruses, has been shown experimentally to increase the pathogenicity of H9N2 viruses in chickens and mice, and has become predominant in Chinese H9N2 circulation (94.8% of viruses since 2013) (Sun et al., 2013). Although the Vietnamese NA genes had a number of additional molecular markers of virulence and airborne transmission in mammals and poultry (Lv et al., 2015; Z. Zhang et al., 2011) (Table S1), they did not have any known molecular markers of oseltamivir resistance, and thus probably remain susceptible to the drug (Gubareva, 2004). Fig. 2 (continued).

3.3.2. Neuraminidase (NA) The N2 genes of the Vietnamese viruses fell into two distinct clades (Fig. 2) within the G57 lineage, clade 0 and 1 (Pu et al., 2014). These distinct clades vary in regards to a 3 amino acid deletion in the NA stalk (clade 1 NAs have the deletion, clade 0 do not). This deletion, although minor compared to the 20+ amino acid deletions found recently in the

3.3.3. NS1 and NEP Similar to divergence of NA, segment 8 also segregated phylogenetically into two groups within the G57-like lineage, with a major group represented by H7F-14-CB4-2 and a minor group represented by H7F14-BN4-315 (Fig. 2). Interestingly, the minor group contained a NS1 C-terminal elongation from 217 to 237 amino acids. This elongation has been shown to influence inflammation in poultry and mammalian experimental models, and to increase transmissibility among poultry (Kong et al., 2015). Similar to the Chinese human H9N2 isolates, the NS1 proteins of these virus contain a large number of mammalian

Please cite this article as: Thuy, D.M., et al., Prevalence and diversity of H9N2 avian influenza in chickens of Northern Vietnam, 2014, Infection, Genetics and Evolution (2016), http://dx.doi.org/10.1016/j.meegid.2016.06.038

D.M. Thuy et al. / Infection, Genetics and Evolution 44 (2016) xxx–xxx

7

markers (Dankar et al., 2011; Ping et al., 2011; Z. Zhang et al., 2011)

Fig. 2 (continued).

Fig. 2 (continued).

(Table S1).

3.3.4. Accessory proteins (PB1-F2, PA-X) In addition to the 10 well characterized and almost universally expressed core proteins of influenza, some virus strains also code for accessory proteins such as PB1-F2 and PA-X, which are thought to be

important in modulating virulence (Chen et al., 2001; Gao, 2015; Jagger et al., 2012). All sequenced PA segments coded for a complete PA-X protein, shown to increase mammalian replication and pathogenicity in H9N2 viruses (Gao, 2015). Length of the PB1-F2 protein, however, varied between different isolates. Of the H9N2-like PB1s, 18 isolates contained a full length PB1-F2 (90aa), 5 had a minor truncation (76aa), and a single isolate had a major truncation (27aa). 76aa long PB1-F2s from H9N2 viruses has previously been implied to increase the virulence of H7N9 viruses in mice but are thought not to interfere with function (Gibbs et al., 2003; Su et al., 2015). PB1-F2 functional knockouts, similar to the 27aa variant seen, have also been shown to increase virulence of avian influenza viruses (Gibbs et al., 2003). The PB1F2 of the mixed isolate H7F-HG4–604 was H5-like and had a length of 57 amino acids, which like the 27aa variants is often considered a functional knockout (Gibbs et al., 2003). All PB1-F2s of over 76aa also had several molecular markers associated with mammalian virulence and predisposition to secondary bacterial pneumonia (Alymova et al., 2014) (Table S1).

Please cite this article as: Thuy, D.M., et al., Prevalence and diversity of H9N2 avian influenza in chickens of Northern Vietnam, 2014, Infection, Genetics and Evolution (2016), http://dx.doi.org/10.1016/j.meegid.2016.06.038

8

D.M. Thuy et al. / Infection, Genetics and Evolution 44 (2016) xxx–xxx

Fig. 2 (continued).

3.3.5. Polymerase subunits and nucleoprotein The polymerase subunits of each of the H9N2 isolates, PB2, PB1, PA, and NP, displayed multiple mammalian adaptation markers (Chen et al., 2006; Mei et al., 2016; Park et al., 2015; Sun et al., 2013; Xiao et al., 2016; Yao et al., 2013) (Table S1). In all sequenced isolates the PA subunit had a marker of airborne transmission in chickens, L672 (Zhong et al., 2014a, 2014b). However, none of the viruses had any of the polymerase markers most well associated with mammalian adaptation of avian influenza viruses, i.e. PB2 E627K or D701N (Herfst et al., 2012; Imai et al., 2012; Mehle and Doudna, 2009; Subbarao et al., 1993).

Fig. 2 (continued).

3.3.6. Matrix proteins The M1 and M2 proteins, encoded as alternative splice products of segment 7 of the influenza genome, similarly to the other segments, contained a number of known markers of mammalian adaptation (Song et al., 2009; Yao et al., 2013; Zhang et al., 2011) (Table S1). All viruses had an asparagine (N) at amino acid position 31 of M2, a mutation associated with resistance to antivirals amantadine and rimantadine (Belshe et al., 1988).

Please cite this article as: Thuy, D.M., et al., Prevalence and diversity of H9N2 avian influenza in chickens of Northern Vietnam, 2014, Infection, Genetics and Evolution (2016), http://dx.doi.org/10.1016/j.meegid.2016.06.038

D.M. Thuy et al. / Infection, Genetics and Evolution 44 (2016) xxx–xxx

9

Fig. 3. Genotype constellations for H9N2. Segments of different colours represent N2% nucleotide difference. Different genotypes were assigned an arbitrary ‘Vietnam genotype #’ (VNG#).

4. Discussion Our results confirm the hyperendemicity of avian influenza viruses in chickens sampled in LBM of northern Vietnam, where influenza H9N2 viruses are by far the dominant subtype in circulation. Our detection rates are comparable to those reported from China, where overall prevalence of AIV in chickens has ranged from 21 to 36% in Hubei, 34– 46% in Zhejiang province (Chen et al., 2016), and 15.7% in Hunan (Huang et al., 2015a), and where H9N2 is the dominant AIV subtype detected throughout all provinces (Shi et al., 2014). Although a plethora of AIV subtypes are known to circulate in wild and domestic Vietnamese ducks, including H3, H4, H5, H6, H7 and H10 viruses (Hotta et al., 2012; Jadhao et al., 2009; D. C. Nguyen et al., 2005; Nomura et al., 2012; Okamatsu et al., 2013; Takakuwa et al., 2013), our deep sequencing data confirm that the only detectable AIV subtypes circulating endemically within chickens are H5 and H9. Additionally, although no Chinese-like H7N9 viruses have presently been found in Vietnam, the H9N2 viruses sequenced here shared the H7N9-like G57/genotype S internal gene cassette (Gu et al., 2014; Pu et al., 2014). This internal gene

Fig. 4. Structural diagram showing antigenically important amino acid substitutions observed in H9HA. Made using pymol, amino acids mapped onto the HA structure of A/ Swine/Hong Kong/9/98 (H9N2) PDBID:1JSD.

cassette is thought to contribute to the pathogenicity of H7N9 viruses in mammals, to confer zoonotic and poultry enzootic potential to other subtypes (e.g. H10N8), and is at least partly responsible for the extremely high transmissibility of H9N2 viruses in chickens (Pu et al., 2014; Zhang et al., 2014). The H9N2 viruses circulating in northern Vietnam also closely resembled recent H9N2 strains isolated from human clinical cases in China (Huang et al., 2015b). Like its more temperate neighbor, China, Vietnam has seasonal fluctuations in AIV prevalence, with peak incidence in HPAI the winter months of November to February (Durand et al., 2015; Nguyen et al., 2014), and a similar pattern for LPAI (Pu et al., 2014). Although the current surveillance system provides an excellent opportunity to monitor diversity of co-circulating HPAI H5 and LPAI H9N2, the process of combining swabs from multiple birds (i.e. pooling) complicates the inference of co-infection rates (and definitely decrease its power), and complicates analysis of intersubtypic reassortment. Among the deep sequencing data generated directly from surveillance swabs, although we could not assess intersubtypic reassortment between H5 and H9N2 viruses, the genomic analysis suggested high levels of intrasubtypic reassortment between different subclades of H9N2 (Liu et al., 2014a, 2014b; Sun et al., 2010). The genomic data demonstrated that most gene segments were closely related to Chinese viruses, which likely reflects proximity of our sampling sites to the Chinese-Vietnam border. Importation of live chickens from China has been illegal since May 2013; hence, these findings suggest that despite regulations to curtail illegal traffic in live poultry, there remains potential for cross-border spread. The pattern of phylogenetic clustering of the H9N2 sequenced here revealed two distinct clusters for each gene segment (or sometime 3 clusters, as in the case of the HA phylogeny), suggesting there may have been separate introduction events. However, H9N2 viruses are known to be highly transmissible between chickens and wild terrestrial birds (e.g. sparrows) (Iqbal et al., 2013), and thus cross-border wild bird migration may also have played a role in geographical spread of H9 viruses. Interestingly, the N2 sequences generated here were predominantly from the BJ94-like lineage, and, unlike all other gene segments, were most similar to Vietnamese N2s that have been circulating in Vietnam for at least a decade suggesting the possibility of underlying continual, in situ evolution of H9N2 viruses within Vietnam. Given the evidence that prevalence and diversity of AIV in neighboring China continues to increase at alarming rates (Chen, 2016; Huang 2015), and given that Vietnam ranks as one of the countries with highest clade diversity for HPAI H5 (Shi et al., 2014), the public health

Please cite this article as: Thuy, D.M., et al., Prevalence and diversity of H9N2 avian influenza in chickens of Northern Vietnam, 2014, Infection, Genetics and Evolution (2016), http://dx.doi.org/10.1016/j.meegid.2016.06.038

10

D.M. Thuy et al. / Infection, Genetics and Evolution 44 (2016) xxx–xxx

concerns regarding potential adaptation of AIV to human transmission remain as pressing today as ever. Yet the numbers of recorded human cases of HPAI H5 in Vietnam are extremely low (total of 127 cases since 2003, and no cases recorded for 2015) (WHO, 2015b). High levels of ‘silent’ circulation of HPAI H5 in ducks have been well documented for many years (Nguyen et al., 2014), and are known to play a key role in endemic transmission ecology. In contrast, however, ‘silent’ shedding of HPAI H5 from apparently healthy chickens may be a relatively recent phenomenon. Asymptomatic shedding of HPAI H5 from chickens could be explained by multiple factors, with at least four possible scenarios: 1) the lack of symptoms may indicate subtle shifts in virulence associated with reassortant genomes (i.e. emergence of attenuated phenotypes); 2) the lack of symptoms could be attributed to use of poultry H5 vaccines that successfully mitigate disease symptoms of H5 infections, but fail to completely block virus replication and shedding; 3) the lack of symptoms may be related to prior heterologous immunity to H9N2, which has been shown to reduce HPAI disease severity but does not impact transmission (Khalenkov et al., 2009; Seo and Webster, 2001); or 4) the chickens that are detected as H5 positive may have become infected after arrival in the market, in which case the lack of symptoms (and the relatively low viral loads) could reflect the stage of infection (early phase incubation). Further examination of these questions would be facilitated by more extensive and systematic testing for both H5 and H9 within the surveillance system, and improved data collection for epidemiological risk factors (i.e. trace-back data on vaccination status, time to sale in market), and modeling approaches to assess overall transmission dynamics of H5 and H9 within the local context of Vietnam. Regarding the human health implications of LPAI H9N2 circulation, previous studies of human seroprevalence conducted in northern Vietnam in 2001 documented HI titers to H9 and presence of H9 neutralizing antibodies in 3–4% of people tested, but with no significant differences between the general population and those with occupational exposures to poultry (Uyeki et al., 2012). We speculate that overall H9 prevalence has increased over the last 15 years, thus warranting revisiting the question of whether H9N2 exposures among people may influence heterosubtypic population immunity and overall influenza transmission dynamics. In conclusion, our findings suggest that circulation of H9N2 in Vietnamese chickens is sufficiently widespread to be influencing the overall ecology and emergence potential of AIV, and we emphasize the need for further research to better understand possible interactions between cocirculating LPAI and HPAI subtypes. We fully support the call for continued vigilance in monitoring ongoing AIV evolution.

Author disclosure statement No competing financial interests exist.

Acknowledgements The authors would like to thank the many staff from Sub-Departments of Animal Health and the national veterinary laboratory network, for their dedication and continued support. The H7F surveillance activities were funded by Food and Agriculture Organization (FAO), and by the BBSRC grant BB/L018853/1, entitled Zoonoses of Emerging Livestock Systems.

Appendix A. Supplementary data Supplementary data to this article can be found online at http://dx. doi.org/10.1016/j.meegid.2016.06.038.

References Alymova, I.V., Samarasinghe, A., Vogel, P., Green, A.M., Weinlich, R., McCullers, J.a., 2014. A novel cytotoxic sequence contributes to influenza A viral protein PB1-F2 pathogenicity and predisposition to secondary bacterial infection. J. Virol. 88 (1), 503–515. http://dx.doi.org/10.1128/JVI.01373-13. Baron, J., Tarnow, C., Mayoli-Nussle, D., Schilling, E., Meyer, D., Hammami, M., ... BottcherFriebertshauser, E., 2013. Matriptase, HAT, and TMPRSS2 activate the hemagglutinin of H9N2 influenza A viruses. J. Virol. 87 (3), 1811–1820. http://dx.doi.org/10.1128/ JVI.02320-12. Belshe, R., Smith, M., Hall, C., Betts, R., Hay, A., 1988. Genetic basic of resistance to rimantadine emerging during treatment of influenza virus infection. J. Virol. 62, 1508–1512. Butt, K.M., Smith, G.J.D., Chen, H., Zhang, L.J., Leung, H.C., Xu, K.M., ... Guan, Y., 2005. Human infection with an avian H9N2 influenza A virus in Hong Kong in 2003 human infection with an avian H9N2 influenza A virus in Hong Kong in 2003. J. Clin. Microbiol. 43 (11), 5760–5767. http://dx.doi.org/10.1128/JCM.43.11.5760. Chen, W.S., Calvo, P.A., Malide, D., Gibbs, J., Schubert, U., Bacik, I., ... Yewdell, J.W., 2001. A novel influenza A virus mitochondrial protein that induces cell death. Nat. Med. 7 (12), 1306–1312. http://dx.doi.org/10.1038/nm1201-1306. Chen, G.-W., Chang, S.-C., Mok, C., Lo, Y.-L., Kung, Y.-N., Huang, J.-H., ... Shih, S.-R., 2006. Genomic signatures of human versus avian influenza A viruses. Emerg. Infect. Dis. 12 (9), 1353–1360. http://dx.doi.org/10.3201/eid1209.060276. Chen, H., Yuan, H., Gao, R., Zhang, J., Wang, D., Xiong, Y., ... Shu, Y., 2014. Clinical and epidemiological characteristics of a fatal case of avian influenza A H10N8 virus infection: a descriptive study. Lancet 383 (9918), 714–721. http://dx.doi.org/10.1016/ S0140-6736(14)60111-2. Chen, L., Chen, L., Lin, X., Tian, J., Wang, W., Wang, M., Holmes, E.C., 2016. Diversity and evolution of avian influenza viruses in live poultry markets , free range poultry and wild wetland birds in China. J. Gen. Virol. (epub ahead). Dankar, S.K., Wang, S., Ping, J., Forbes, N.E., Keleta, L., Li, Y., Brown, E.G., 2011. Influenza A virus NS1 gene mutations F103L and M106I increase replication and virulence. Virol. J. 8, 13. http://dx.doi.org/10.1186/1743-422X-8-13. Durand, L.O., Glew, P., Gross, D., Kasper, M., Trock, S., Kim, I.K., ... Azziz-Baumgartner, E., 2015. Timing of influenza A (H5N1) in poultry and humans and seasonal influenza activity worldwide, 2004–2013. Emerg. Infect. Dis. 21 (2), 202–208. http://dx.doi. org/10.3201/eid2102.140877. FAO, 2014. EMPRES animal health 360. EMPRES-Animal Health 360. 44 (Rome). Fusaro, A., Monne, I., Salviato, A., Schivo, A., Amarin, N.M., Ismail, M.M., ... Gonzalez, C., 2011. Phylogeography and Evolutionary History of Reassortant H9N2 Viruses With Potential Human Health Implications Phylogeography and Evolutionary History of Reassortant H9N2 Viruses With Potential Human Health Implications. http://dx. doi.org/10.1128/JVI.00219-11. Gao, H., 2015. Pa-X is a virulence factor in avian H9N2 influenza virus. J. Gen. Virol. C, 2–3. http://dx.doi.org/10.1099/jgv.0.000016. Gibbs, J.S., Malide, D., Hornung, F., Bennink, J.R., Yewdell, J.W., 2003. The influenza A virus PB1-F2 protein targets the inner mitochondrial membrane via a predicted basic amphipathic helix that disrupts mitochondrial function. J. Virol. 77 (13), 7214–7224. http://dx.doi.org/10.1128/JVI.77.13.7214-7224.2003. Gu, M., Chen, H., Li, Q., Huang, J., Zhao, M., Gu, X., ... Liu, X., 2014. Enzootic genotype S of H9N2 avian influenza viruses donates internal genes to emerging zoonotic influenza viruses in China. Vet. Microbiol. 174 (3–4), 309–315. http://dx.doi.org/10.1016/j. vetmic.2014.09.029. Gubareva, L.V., 2004. Molecular mechanisms of influenza virus resistance to neuraminidase inhibitors. Virus Res. 103 (1–2), 199–203. http://dx.doi.org/10.1016/j.virusres. 2004.02.034. Herfst, S., Schrauwen, E.J.a., Linster, M., Chutinimitkul, S., de Wit, E., Munster, V.J., ... Fouchier, R.a.M., 2012. Airborne transmission of influenza A/H5N1 virus between ferrets. Science 336 (6088), 1534–1541. http://dx.doi.org/10.1126/science.1213362 (New York, N.Y.). Hotta, K., Takakuwa, H., Le, Q.M.T., Phuong, S.L., Murase, T., Ono, E., ... Yamashiro, T., 2012. Isolation and characterization of H6N1 and H9N2 avian influenza viruses from Ducks in Hanoi, Vietnam. Virus Res. 163 (2), 448–453. http://dx.doi.org/10.1016/j.virusres. 2011.11.004. Huang, Y., Li, X., Zhang, H., Chen, B., Jiang, Y., Yang, L., ... Gao, L., 2015a. Human infection with an avian influenza A (H9N2) virus in the middle region of China. J. Med. Virol. 1648 (May), 1641–1648. http://dx.doi.org/10.1002/jmv. Huang, Y., Zhang, H., Li, X., Hu, S., Cai, L., Sun, Q., ... Gao, L., 2015b. Detection and genetic characteristics of H9N2 avian influenza viruses from live poultry markets in Hunan Province, China. PLoS One 10 (11), e0142584. http://dx.doi.org/10.1371/journal. pone.0142584. ICDDRB, 2011. Outbreak of mild respiratory disease caused by H5N1 and H9N2 infections among young children in Dhaka, Bangladesh. Health and Science Bulletin. Vol. 9. Imai, M., Watanabe, T., Hatta, M., Das, S.C., Ozawa, M., Shinya, K., ... Kawaoka, Y., 2012. Experimental adaptation of an influenza H5 HA confers respiratory droplet transmission to a reassortant H5 HA/H1N1 virus in ferrets. Nature 486 (7403), 420–428. http://dx.doi.org/10.1038/nature10831. Iqbal, M., Yaqub, T., Mukhtar, N., Shabbir, M.Z., McCauley, J.W., 2013. Infectivity and transmissibility of H9N2 avian influenza virus in chickens and wild terrestrial birds. Vet. Res. 44 (1), 1. http://dx.doi.org/10.1186/1297-9716-44-100. Jadhao, S.J., Nguyen, D.C., Uyeki, T.M., Shaw, M., Maines, T., Rowe, T., ... Katz, J.M., 2009. Genetic analysis of avian influenza A viruses isolated from domestic waterfowl in livebird markets of Hanoi, Vietnam, preceding fatal H5N1 human infections in 2004. Arch. Virol. 154 (8), 1249–1261. http://dx.doi.org/10.1007/s00705-009-0429-2. Jagger, B., Wise, H., Kash, J., Walters, K., Wills, N., 2012. An overlapping protein-coding region in influenza A protein that induces cell death. Science 337, 199–204. Kaverin, N.V., Rudneva, I.A., Ilyushina, N.A., Lipatov, A.S., Krauss, S., Webster, R.G., 2004. Structural differences among hemagglutinins of influenza A virus subtypes are

Please cite this article as: Thuy, D.M., et al., Prevalence and diversity of H9N2 avian influenza in chickens of Northern Vietnam, 2014, Infection, Genetics and Evolution (2016), http://dx.doi.org/10.1016/j.meegid.2016.06.038

D.M. Thuy et al. / Infection, Genetics and Evolution 44 (2016) xxx–xxx reflected in their antigenic architecture : analysis of H9 escape mutants. 78 (1), 240–249. http://dx.doi.org/10.1128/JVI.78.1.240. Khalenkov, A., Perk, S., Panshin, A., Golender, N., Robert, G., 2009. Modulation of the severity of highly pathogenic H5N1 influenza in chickens previously inoculated with Israeli H9N2 influenza viruses. Virology 383 (1), 32–38. http://dx.doi.org/10.1016/j.virol. 2008.09.026.Modulation. Kim, K., Choi, J., Kang, H., To, T. L., Nguyen, T.D., Song, B., ... Lee, Y., 2013. Geographical distribution of low pathogenic avian influenza viruses of domestic poultry in Vietnam and their genetic relevance with Asian isolates. Poult. Sci. 92, 2012–2023. Kong, W., Liu, L., Wang, Y., He, Q., Wu, S., 2015. C-terminal elongation of NS1 of H9N2 influenza virus induces a high level of inflammatory cytokines and increases transmission. J. Gen. Virol. 96, 259–268. Li, C., Yu, K., Tian, G., Yu, D., Liu, L., Jing, B., ... Chen, H., 2005. Evolution of H9N2 influenza viruses from domestic poultry in Mainland China. Virology 340 (1), 70–83. http://dx. doi.org/10.1016/j.virol.2005.06.025. Liu, D., Shi, W., Gao, G.F., 2014a. Poultry carrying H9N2 act as incubators for novel human avian influenza viruses. Lancet 383 (9920), 869. http://dx.doi.org/10.1016/S01406736(14)60386-X. Liu, J., Zhang, J., Huang, F., Zhang, Y., Luo, H., Zhang, H., 2014b. Complex reassortment of polymerase genes in Asian influenza A virus H7 and H9 subtypes. Infect. Genet. Evol. 23, 203–208. http://dx.doi.org/10.1016/j.meegid.2014.02.016. Lv, J., Wei, L., Yang, Y., Wang, B., Liang, W., Gao, Y., ... Chai, T., 2015. Amino acid substitutions in the neuraminidase protein of an H9N2 avian influenza virus affect its airborne transmission in chickens. Vet. Res. 46 (1). http://dx.doi.org/10.1186/s13567-014-0142-3. Matrosovich, M., Tuzikov, a., Bovin, N., Gambaryan, a., Klimov, a., Castrucci, M.R., ... Kawaoka, Y., 2000. Early alterations of the receptor-binding properties of H1, H2, and H3 avian influenza virus hemagglutinins after their introduction into mammals. J. Virol. 74 (18), 8502–8512. http://dx.doi.org/10.1128/JVI.74.18.8502-8512.2000. Matrosovich, M.N., Krauss, S., Webster, R.G., 2001. H9N2 influenza A viruses from poultry in Asia have human virus-like receptor specificity. Virology 281 (2), 156–162. http:// dx.doi.org/10.1006/viro.2000.0799. Mehle, A., Doudna, J.A., 2009. Adaptive strategies of the influenza virus polymerase for replication in humans. Proc. Natl. Acad. Sci. U. S. A. 106 (50), 21312–21316. http:// dx.doi.org/10.1073/pnas.0911915106. Mei, K., Liu, G., Chen, Z., Gao, Z., Zhao, L., Jin, T., ... Chen, Q., 2016. Deep sequencing reveals the viral adaptation process of environment-derived H10N8 in mice. Infect. Genet. Evol. 37, 8–13. http://dx.doi.org/10.1016/j.meegid.2015.10.016. Nguyen, D.C., Uyeki, T.M., Jadhao, S., Maines, T., Shaw, M., Matsuoka, Y., ... Katz, J.M., 2005. Isolation and characterization of avian influenza viruses, including highly pathogenic H5N1, from poultry in live bird markets in Hanoi, Vietnam, in 2001. J. Virol. 79 (7), 4201–4212. http://dx.doi.org/10.1128/JVI.79.7.4201-4212.2005. Nguyen, D.T., Bryant, J.E., Davis, C.T., Nguyen, L.V., Pham, L.T., Loth, L., ... E, F.D.V.P., 2014. Prevalence and distribution of avian influenza A (H5N1) virus clade variants in live bird markets of Vietnam, 2011–2013. Avian Dis. 58 (4), 599–608. Nomura, N., Sakoda, Y., Endo, M., Yoshida, H., Yamamoto, N., Okamatsu, M., ... Kida, H., 2012. Characterization of avian influenza viruses isolated from domestic ducks in Vietnam in 2009 and 2010. Arch. Virol. 157 (2), 247–257. http://dx.doi.org/10. 1007/s00705-011-1152-3. Okamatsu, M., Sakoda, Y., Kishida, N., Isoda, N., Kida, H., 2008. Antigenic structure of the hemagglutinin of H9N2 influenza viruses. Arch. Virol. 153 (12), 2189–2195. Okamatsu, M., Nishi, T., Nomura, N., Yamamoto, N., Sakoda, Y., Sakurai, K., ... Kida, H., 2013. The genetic and antigenic diversity of avian influenza viruses isolated from domestic ducks, muscovy ducks, and chickens in northern and southern Vietnam, 2010–2012. Virus Genes 47 (2), 317–329. http://dx.doi.org/10.1007/s11262-013-0954-7. Park, K.J., Song, M.-S., Kim, E.-H., Kwon, H., Baek, Y.H., Choi, E., ... Choi, Y.K., 2015. Molecular characterization of mammalian-adapted Korean-type avian H9N2 virus and evaluation of its virulence in mice. J. Microbiol. 53 (8), 570–577. http://dx.doi.org/ 10.1007/s12275-015-5329-4. Parvin, R., Heenemann, K., Halami, M.Y., Chowdhury, E.H., Islam, M.R., Vahlenkamp, T.W., 2014. Full-genome analysis of avian influenza virus H9N2 from Bangladesh reveals internal gene reassortments with two distinct highly pathogenic avian influenza viruses. Arch. Virol. 159 (7), 1651–1661. http://dx.doi.org/10.1007/s00705-014-1976-8. Peacock, T., Reddy, K., James, J., Adamiak, B., Barclay, W., Shelton, H., Iqbal, M., 2016. Antigenic mapping of an H9N2 avian influenza virus reveals two discrete antigenic sites and a novel mechanism of immune escape. Sci. Report. 6 (October 2015), 18745. http://dx.doi.org/10.1038/srep18745. Ping, J., Li, C., Deng, G., Jiang, Y., Tian, G., Zhang, S., ... Chen, H., 2008. Single-amino-acid mutation in the HA alters the recognition of H9N2 influenza virus by a monoclonal antibody. Biochem. Biophys. Res. Commun. 371 (1), 168–171. http://dx.doi.org/10. 1016/j.bbrc.2008.04.045. Ping, J., Keleta, L., Forbes, N.E., Dankar, S., Stecho, W., Tyler, S., ... Brown, E.G., 2011. Genomic and protein structural maps of adaptive evolution of human influenza A virus to increased virulence in the mouse. PLoS One 6 (6), e21740. http://dx.doi.org/10.1371/ journal.pone.0021740. Pu, J., Wang, S., Yin, Y., Zhang, G., Carter, R.A., Wang, J., ... Webster, R.G., 2014. Evolution of the H9N2 influenza genotype that facilitated the genesis of the novel H7N9 virus. PNAS http://dx.doi.org/10.1073/pnas.1422456112. Rudneva, I.a., Ilyushina, N.a., Timofeeva, T.a., Webster, R.G., Kaverin, N.V., 2005. Restoration of virulence of escape mutants of H5 and H9 influenza viruses by their readaptation to mice. J. Gen. Virol. 86 (Pt 10), 2831–2838. http://dx.doi.org/10.1099/vir.0.81185-0. Sang, X., Wang, A., Ding, J., Kong, H., Gao, X., Li, L., ... Xia, X., 2015. Adaptation of H9N2 AIV in guinea pigs enables efficient transmission by direct contact and inefficient transmission by respiratory droplets. Sci. Report. 5 (November), 15928. http://dx.doi. org/10.1038/srep15928.

11

Seo, S.H., Webster, R.G., 2001. Cross-reactive, cell-mediated immunity and protection of chickens from lethal H5N1 influenza virus infection in Hong Kong poultry markets. 75 (6), 2516–2525. http://dx.doi.org/10.1128/JVI.75.6.2516. Shi, W., Li, W., Li, X., Haywood, J., Ma, J., Gao, G.F., Liu, D., 2014. Phylogenetics of varied subtypes of avian influenza viruses in China: potential threat to humans. Protein Cell 5 (4), 253–257. http://dx.doi.org/10.1007/s13238-014-0036-1. Song, M.-S., Pascua, P.N.Q., Lee, J.H., Baek, Y.H., Lee, O.-J., Kim, C.-J., ... Choi, Y.K., 2009. The polymerase acidic protein gene of influenza a virus contributes to pathogenicity in a mouse model. J. Virol. 83 (23), 12325–12335. http://dx.doi.org/10.1128/JVI.01373-09. Sorrell, E.M., Wan, H., Araya, Y., Song, H., Perez, D.R., 2009. Minimal molecular constraints for respiratory droplet transmission of an avian-human H9N2 influenza A virus. PNAS 106 (18), 7565–7570. http://dx.doi.org/10.1073/pnas.0900877106. Stevens, J., Blixt, O., Tumpey, T.M., Taubenberger, J.K., Paulson, J.C., Wilson, I.a., 2006. Structure and receptor specificity of the hemagglutinin from an H5N1 influenza virus. Science 312 (5772), 404–410. http://dx.doi.org/10.1126/science.1124513 (New York, N.Y.). Su, W., Wang, C., Luo, J., Zhao, Y., Wu, Y., Chen, L., ... He, H., 2015. Testing the effect of internal genes derived from a wild-bird-origin H9N2 influenza A virus on the pathogenicity of an A/H7N9 virus. Cell Rep. 12 (11), 1831–1841. http://dx.doi.org/10.1016/j. celrep.2015.08.029. Subbarao, E.K., London, W., Murphy, B.R., 1993. A single amino acid in the PB2 gene of influenza A virus is a determinant of host range. J. Virol. 67 (4), 1761–1764. Sun, Y., Pu, J., Jiang, Z., Guan, T., Xia, Y., Xu, Q., ... Liu, J., 2010. Genotypic evolution and antigenic drift of H9N2 influenza viruses in China from 1994 to 2008. Vet. Microbiol. 146 (3–4), 215–225. http://dx.doi.org/10.1016/j.vetmic.2010.05.010. Sun, Y., Tan, Y., Wei, K., Sun, H., Shi, Y., Pu, J., ... Liu, J., 2013. Amino acid 316 of hemagglutinin and the neuraminidase stalk length influence virulence of H9N2 influenza virus in chickens and mice. J. Virol. 87 (5), 2963–2968. http://dx.doi.org/10.1128/JVI.02688-12. Takakuwa, H., Yamashiro, T., Le, M.Q., Phuong, L.S., Ozaki, H., Tsunekuni, R., ... Otsuki, K., 2013. The characterization of low pathogenic avian influenza viruses isolated from wild birds in northern Vietnam from 2006 to 2009. Comp. Immunol. Microbiol. Infect. Dis. 36 (6), 581–590. http://dx.doi.org/10.1016/j.cimid.2013.06.004. Uyeki, T.M., Nguyen, D.C., Rowe, T., Lu, X., Hu-Primmer, J., Huynh, L.P., ... Katz, J.M., 2012. Seroprevalence of antibodies to avian influenza A (H5) and A (H9) viruses among market poultry workers, Hanoi, Vietnam, 2001. PLoS One 7 (8), e43948. http://dx. doi.org/10.1371/journal.pone.0043948. Wan, Z., Ye, J., Xu, L., Shao, H., Jin, W., Qian, K., ... Qin, A., 2014. Antigenic mapping of the hemagglutinin of an H9N2 avian influenza virus reveals novel critical amino acid positions in antigenic sites. J. Virol. 88 (7), 3898–3901. http://dx.doi.org/10.1128/JVI.03440-13. Wei, Y., Xu, G., Zhang, G., Wen, C., Anwar, F., Wang, S., ... Pu, J., 2016. Antigenic evolution of H9N2 chicken influenza viruses isolated in China during 2009–2013 and selection of a candidate vaccine strain with broad cross-reactivity. Vet. Microbiol. 182, 1–7. http://dx.doi.org/10.1016/j.vetmic.2015.10.031. WHO, 2015a. Influenza at the human-animal interface, summary and assessment as of 3 March 2015. WHO Monthly Human Animal Interface Reports. WHO, 2015b. WHO. Cumulative Number of Confirmed Human Cases for Avian Influenza H5N1. (Retrieved from http://www.who.int/influenza/human_animal_interface/EN_ GIP_20151214cumulativeNumberH5N1cases.pdf). WHO Global Influenza Surveillance, 2011. Manual for the Laboratory Diagnosis and Virological Surveillance of Influenza. (Retrieved from http://whqlibdoc.who.int/ publications/2011/9789241548090_eng.pdf). Xiao, C., Ma, W., Sun, N., Huang, L., Li, Y., Zeng, Z., ... Liao, M., 2016. PB2-588V promotes the mammalian adaptation of H10N8, H7N9 and H9N2 avian influenza viruses. Sci. Report. 6 (August 2015), 19474. http://dx.doi.org/10.1038/srep19474. Yao, Y., Wang, H., Chen, Q., Zhang, H., Zhang, T., 2013. Characterization of low-pathogenic H6N6 avian influenza viruses in central China. Arch. Virol. 158, 367–377. Zaraket, H., Baranovich, T., Kaplan, B.S., Carter, R., Song, M.-S., Paulson, J.C., ... Webby, R.J., 2015. Mammalian adaptation of influenza A (H7N9) virus is limited by a narrow genetic bottleneck. Nat. Commun. 6, 6553. http://dx.doi.org/10.1038/ncomms7553. Zhang, P., Tang, Y., Liu, X., Peng, D., Liu, W., Liu, H., ... Liu, X., 2008. Characterization of H9N2 influenza viruses isolated from vaccinated flocks in an integrated broiler chicken operation in eastern China during a 5 year period (1998–2002). J. Gen. Virol. 89 (Pt 12), 3102–3112. http://dx.doi.org/10.1099/vir.0.2008/005652-0. Zhang, Z., Hu, S., Li, Z., Wang, X., Liu, M., Guo, Z., ... Jin, H., 2011. Multiple amino acid substitutions involved in enhanced pathogenicity of LPAI H9N2 in mice. Infect. Genet. Evol. 11 (7), 1790–1797. http://dx.doi.org/10.1016/j.meegid.2011.07.025. Zhang, T., Bi, Y., Tian, H., Li, X., Liu, D., Wu, Y., ... Xu, B., 2014. Human infection with influenza virus A (H10N8) from live poultry markets, China, 2014. Emerg. Infect. Dis. 20 (12), 2012–2015. Zhong, L., Wang, X., Li, Q., Liu, D., Chen, H., Zhao, M., ... Liu, X., 2014a. Molecular mechanism of the airborne transmissibility of H9N2 avian influenza A viruses in chickens. J. Virol. 88 (17), 9568–9578. http://dx.doi.org/10.1128/JVI.00943-14. Zhong, L., Wang, X., Li, Q., Liu, D., Zhao, M., Gu, X., ... Gu, M., 2014b. Molecular mechanism of the airborne transmissibility of H9N2 avian influenza A in chickens. J. Virol. 88 (77), 9568. http://dx.doi.org/10.1128/JVI.00943-14. Zhou, B., Donnelly, M.E., Scholes, D.T., St George, K., Hatta, M., Kawaoka, Y., Wentworth, D.E., 2009. Single-reaction genomic amplification accelerates sequencing and vaccine production for classical and swine origin human influenza a viruses. J. Virol. 83 (19), 10309–10313. http://dx.doi.org/10.1128/JVI.01109-09. Zhu, Y., Yang, D., Ren, Q., Yang, Y., Liu, X., Xu, X., ... Liu, X., 2015. Identification and characterization of a novel antigenic epitope in the hemagglutinin of the escape mutants of H9N2 avian influenza viruses. Vet. Microbiol. 178 (1–2), 144–149. http://dx.doi.org/ 10.1016/j.vetmic.2015.04.012.

Please cite this article as: Thuy, D.M., et al., Prevalence and diversity of H9N2 avian influenza in chickens of Northern Vietnam, 2014, Infection, Genetics and Evolution (2016), http://dx.doi.org/10.1016/j.meegid.2016.06.038