Time-space analysis of highly pathogenic avian influenza H5N2 outbreak in the US
© The Author(s). 2016
Received: 28 March 2016
Accepted: 23 August 2016
Published: 30 August 2016
In early 2015, highly pathogenic avian influenza H5N2 caused outbreaks in commercial poultry farms in Minnesota and neighboring states where more than 48 million birds were affected. To date, the origin and transmission pathways of HPAI H5N2 have not been conclusively established.
In this study, we analyzed forty-six samples from turkeys and their environment that were collected at different time-points of the outbreak to identify origins and within outbreak evolutionary changes. We performed de-novo whole genome sequencing from primary samples and the most recent common ancestors of the PB2, PA, HA5, M and NS segments were traced back to Japanese HPAI H5N8 isolates. These segments appeared to have diverged from the ancestor around June and November 2014.
The time to most recent common ancestor analysis for PB1, NP and NA2 segments suggest two likely possibilities of reassortant HPAI H5N2 origin - either a reassortment in Alaska area or multiple reassortments with North American low pathogenic avian influenza strains, before the HPAI H5N2 outbreak strain emerged. Within the outbreak, viruses clustered into two and three subgroups suggesting high substitution rates of 0.702x10-2 - 1.665x10-2 (subs/site/year), over the 5-month outbreak period.
Data are suggestive of a fast evolving HPAI strain within an outbreak that should be taken into consideration in developing appropriate control strategies in the future.
KeywordsEvolutionary H5N2 Highly pathogenic avian influenza Minnesota
Avian influenza virus (AIV) is an enveloped virus that contains eight negative sense single strand RNA segments which encode at least ten functional proteins . Hemagglutinin (HA) and neuraminidase (NA) are surface proteins that are used to classify AIV into subtypes. At present, 16 HA and 9 NA subtypes have been identified in aquatic wild birds which are a natural reservoir of these viruses . AIV can be divided into low pathogenic avian influenza viruses (LPAI) and highly pathogenic avian influenza viruses (HPAI) by its ability to cause disease in poultry which is identified by intravenous pathogenicity index (IVPI) test or by it’s possession of poly-basic amino acid feature at the HA cleavage site .
Outbreaks HPAI H5 in Guangdong, China in 1996 have likely origins from goose Guangdong lineage HPAI viruses . Since 1997, HPAI H5N1 dispersed in more than 50 countries in Africa, Asia, Europe, and the Middle East. In 2003, a re-emergence of HPAI H5N1 in China and South East Asia was reported. Since its emergence, the H5 gene has continued to evolve and has been arbitrarily classified in the OIE nomenclature system into clades identified from 0 to 9 . The latest clade, 126.96.36.199, was designated in January 2015 and replaced the provisional clade 188.8.131.52. The first isolation of H5 clade 184.108.40.206 viruses was from domestic mallard ducks (Anas platyrhynchos) in China in 2008 . The H5 clade 220.127.116.11 has demonstrated an ability to reassort with multiple neuraminidase subtypes including N1, N2, N3, N5, N6 or N8 .
In January 2014, HPAI H5N8 caused outbreaks in 161 commercial poultry flocks in South Korea and led to culling of 14,000,000 birds . In early November 2014, HPAI H5N8 was detected in a turkey flock in Germany followed by an outbreak in a duck farm in England and a chicken farm in the Netherlands. Wild birds may play an important role in generating novel reassortant subtypes of clade 18.104.22.168 and carrying viruses across continents . In late November 2014, HPAI H5N2 caused outbreaks in turkey and broiler breeder flock in British Columbia, Canada. In the United States, the first case of HPAI H5 was detected in a captive gyrfalcon (Falco rusticolus) in Washington in early December and subsequently the first detection of the reassortant HPAI H5N2 was reported in a northern pintail duck (Anas acuta) in Washington. In December 2014 and January 2015, HPAI H5 viruses associated with various NA subtypes were detected in backyard poultry flocks in Oregon and Washington, USA . In January 2015, an HPAI H5N8 virus was detected in a commercial turkey flock in California . On March 2nd, 2015, the first case of HPAI H5N2 in the Midwestern USA was confirmed in Pope County, Minnesota. On June 17th, 2015, last confirmed case in this outbreak was reported in Iowa. In the Midwest region, more than 200 confirmed cases were reported and more than 48,000,000 birds were affected .
In this study, we investigated the origin and within outbreak evolution of this new emergent HPAI H5N2 that caused outbreaks in Midwestern states.
Samples were collected from turkeys and the environment at locations involved in outbreaks of HP clade 22.214.171.124 H5N2. Sampling was performed over a 106-day period during the outbreak. A total of 46 tissue or oropharyngeal or environmental samples were collected. Samples were divided into three phases representing an approximate 35-day interval each - early phase (March 4th–April 7th) [n = 18], mid phase (April 8th–May 12th) [n = 19] and late phase (May 13th–June 17th) [n = 9].
Whole genome sequencing
Forty-six samples from turkeys, drinker biofilm, air and environment submitted to the University of Minnesota Mid Central Research and Outreach Center (Willmar, MN) were used in this study. RNA was extracted with the MagMAX™-96 Viral RNA Isolation Kit (Ambion) using a magnetic particle processor (Kingfisher, model 700) according to manufacturer’s instructions. All RNA samples were tested by RT-PCR for matrix gene as described . The positive RNA samples were subjected to amplify all eight segments of virus simultaneously using a one-step RT-PCR as described . Briefly, SuperScript® III one-Step RT-PCR system with Platinum® Taq DNA polymerase (Invitrogen™; CA, USA) was used with 1.6 mM and 0.2 μM final concentration of Magnesium and primers, respectively. PCR products were imaged by electrophoresis in 1.5 % agarose gel and subsequently purified by QIAquick PCR Purification Kit (Qiagen®; Hilden, Germany). Purified PCR products were submitted to University of Minnesota Genomics Center (UMGC) for illumina paired-end 250 cycles sequencing by using Nextera XT DNA kit for library generation. Sequences were assembled by mapping all reads to a reference and whole- or partial-genome sequences were extracted via CLC genomics workbench module available on Minnesota supercomputing institute (MSI) resources at the University of Minnesota.
The reference nucleotide sequences were obtained from Influenza Research Database (http://www.fludb.org/) and GISAID (http://gisaid.org/) in May 2016. Reference sequences were selected to represent previous and recent avian influenza strains from varying geographic areas including North America and Eurasia. The following approach was applied to select reference sequences: 1) include an ancestral strain of clade 2.3.4; 2) use of sequences of isolates selected from 3 well characterized H5 (126.96.36.199 lineage) outbreaks in Asia, Europe and North America; and 3) use of a double selection criteria to identify the top 30 hits to the current outbreak isolates by BLAST and develop a tree. Subsequently, these closely related sequences were reanalyzed by BLAST to expand the reference database and used to reconstruct phylogeny. LPAI North American strains that were isolated during 2012–2015 were also included. Reference and sample sequences from the current study were aligned using Muscle v.3.6  subsequently any extra sequences beyond and after start and stop codon were trimmed. Maximum clade credibility (MCC) tree of each gene segments were generated by BEAST 1.8 with Bayesian Markov Chain Monte Carlo (BMCMC) algorithm. Strict clock model with coalescent constant population and HKY with gamma 4 substitution was used as model parameters [16, 17].
Mean substitution rate was estimated by Bayesian coalescent with constant population size  and strict clock model was applied. The Bayesian MCMC chain lengths were 10,000,000 generations with sampling every 10,000 generations and the effective sample size (ESS) value was assessed by using Tracer (v1.6.0) (Molecular evolution, phylogenetics and epidemiology, Edinburgh, Scotland, UK) . Every gene segment analysis had ESS value greater than 200 to suggesting minimal standard error. The resulting tree of each iteration was summarized for a representative clustering pattern by using a tree annotator with 10 % discarding of the chains as burn-in and the resulting maximum clade credibility tree was visualized with FigTree software (v1.4.2) (Molecular evolution, phylogenetics and epidemiology, Edinburgh, Scotland, UK).
First case of this series of outbreaks in the Midwest area was confirmed on March 4th, 2015 in Pope county, Minnesota and last case on June 17th, 2015 in Iowa . A 106-day period was divided into three phases representing an approximate 35-day interval - early phase (March 4th–April 7th), mid phase (April 8th–May 12th) and late phase (May 13th–June 17th). Eighteen, nineteen and nine samples were collected during early, mid and late phases of the outbreak, respectively (Table 1).
Origin of HPAI H5N2
Within outbreak evolution
Mean nucleotide substitution rate of the H5N2 epidemic
Mean Substitution rate (×10−2)
Substitution rate 95 % HPD (×10−2)
Genetic analysis for amino acid variability
Genetic analysis of the HA5 segment indicated multiple basic amino acid at HA cleavage site. Samples from early and mid phase of the outbreak contain the motif PLRERRRKR/GLF that is a characteristic of H5 clade 2.3.4, while samples collected during late phases of the outbreak carried the PQRERRRKR/GLF motif Glutamine (Q) is the major amino acid of other H5 clades . Amino acid substitutions were analyzed by comparing with A/Northern pintail/Washington/40964/2014 (A/NP/WA/14), the first isolate of HPAI H5N2 EA/NA in the US. At position 141 (H5 numbering system) on the HA gene, A/NP/WA/14 and sub-cluster I contained serine (S) while 17 samples, sub-cluster II and III, contain proline (P). S141P mutations are related to surface protein stability and reduced host immune response . The deduced amino acid sequence of the HA stalk region of all samples were identical.
All samples contained S at position 31 on M2 that has been shown to confer reduced susceptibility to amantadine and rimantadine antiviral drugs. However, analysis of resistance to neuraminidase inhibitors (NAI) at position 119, 151, 222, 224, 276, 292 and 371 on NA2 gene segment  showed that all samples in this study contain the conserved amino acid at every position which are expected to maintain NAI susceptibility. The E627K and D701N substitutions in the PB2 gene segment that are expected to increase virulence in mammals were examined. The results showed that all samples contain E and D in respective amino acid sites suggesting low mammalian virulence in this group of viruses.
A study addressing the spread of HPAI H5N8 from South Korea to Europe and Japan have suggested that these outbreaks originated from a single source population . Dispersion of the virus to different geographic localities coincides with wild bird migratory seasons strongly suggesting their important role in the transmission of this pathogen and has been suggested as a mechanism by which H5N8 entered South Korea . The most likely flyway involved in the recent spread of HPAI EA H5N8 to Japan was the East Asian Australian (EAA) Flyway . The EAA flyway extends from north-eastern Asia and western Alaska to the southern extent of Australia and New Zealand and Mongolia, western China and eastern India . On the North America continent, Alaska is an area where there is major overlap between migratory birds from the EAA and Pacific Flyways  but there exists no evidence of HP H5 clade 188.8.131.52 in Alaska. The existence of PB1, NP and NA2 ancestor that related to HP H5N2 in Alaska arise the possibility that HPAI H5N2 EA/NA originated in a location that coincides with Alaska and subsequently spread to British Columbia, Canada and Washington, USA. However, the mechanism and timing of viral transfer from the Pacific Flyway to the Mississippi Flyway in NA has not been established.
The ancestral strain and the overlapping of interval between TMRCA estimates of PB1, NP and NA2 suggest two hypothesis of HPAI EA H5N2 origin. First, it is possible that reassortant event occurred in Alaska. This finding is concordant with the report of avian influenza surveillance in Alaska during spring and summer 2015 that found reassortant LPAI between EA and NAm . Second, the outbreak viral genome constellation likely occurred via multiple reassortment events with an LPAI before addition of genes of HPAI EA H5N2 origin due to the earlier divergence time of PB1, NP and NA2. Furthermore, HA5 clade 184.108.40.206 has the potential to reassort with multiple NA subtypes including NA1, NA2, NA5, NA6 and NA8 . For example, reassortant H5N1 (EA/NA H5N1) containing PB2, HA5, NP and M from EA H5N8 and the rest of the gene segments (PB1, PA, NA1 and NS) from AIV North American lineage was detected in a green wing teal and an American wigeon in December 2014, with no subsequent case reports. This finding indicates that EA/NA H5N1 was likely unable to persist in the population while EA/NA H5N2 provided a compatible constellation for maintenance. Alternately, this may indicate that there was insufficient surveillance to detect the fate of that reassortment.
The nucleotide substitution rate of all segments within the outbreak are greater than 5.000 × 10−3 substitution/site/year which are higher than the rates estimated in long term analysis . These results indicated that viruses were evolving rapidly during the current outbreak. Host jump from wild birds to domestic poultry including turkeys and chickens is a factor that provoked this rapid evolution. Therefore, identification of bridge of interspecies transmission is an important process to prevent the further outbreaks. While the HA gene had an extremely high substitution rate, the HA stalk region of all samples were identical at the amino acid level. The HA stalk region, PA and NP proteins were highly conserved, opening avenues for a universal subunit vaccine concept and development.
At present, no human case has been reported due to HPAI H5N2 EA/NA and it is considered as low health risk to public health . However, three human cases with severe respiratory disease caused by HPAI H5N6 HA5 clade 220.127.116.11 have been reported in China . These cases confirmed that HA5 clade 18.104.22.168 is capable of causing human infection and thus ongoing monitoring of viral evolution should be performed.
Wild bird surveillance in United States has been performed annually during September and November, which is a migratory season and only fecal or cloacal swabs are traditionally collected. After HPAI H5 clade 22.214.171.124 was detected in December 2014, USDA/APHIS wildlife services responded by enhancing surveillance in Pacific Flyway and detected more viruses. Furthermore, multiple introductions of this clade of viruses cannot be ruled out. Therefore, there is a need for more comprehensive wild bird surveillance effort to address a complex wild-domestic animal-human interface to capture variations in subtype specific organ predilection and varying patterns of shedding by different host species.
The authors would like to thank Chulalongkorn University for a Dutsadi Phiphat scholarship to NN and TRF senior scholar to AA (RTA5680006). We also thank Drs. Toremorrell, Murtaugh, Marthaler and Mor for samples and critical discussions on our analysis. This project was funded by Agricultural Experiment Station grants expended as Rapid Agricultural Response fund awarded to CC.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Lamb BA, Krug RM. Orthomyxoviridae: The viruses and Their Replication. in: Fundamental virology. In: Fields BN, Knipe PM, Howley PM, editors. Fundamental virology. Third editionth ed. Philadelphia: Lippincott-Raven; 1996. p. 605–47.Google Scholar
- Webster RG, Bean WJ, Gorman OT, Chambers TM, Kawaoka Y. Evolution and ecology of influenza A viruses. Microbiol Rev. 1992;56:152–79.PubMedPubMed CentralGoogle Scholar
- Horimoto T, Kawaoka Y. Reverse genetics provides direct evidence for a correlation of hemagglutinin cleavability and virulence of an avian influenza A virus. J Virol. 1994;68:3120–8.PubMedPubMed CentralGoogle Scholar
- World Health Organization Global Influenza Program Surveillance N. Evolution of H5N1 avian influenza viruses in Asia. Emerg Infect Dis. 2005;11:1515–21.View ArticleGoogle Scholar
- Updated unified nomenclature system for the highly pathogenic H5N1 avian influenza viruses. http://www.who.int/influenza/gisrs_laboratory/h5n1_nomenclature/en/
- Gu M, Liu W, Cao Y, Peng D, Wang X, Wan H, Zhao G, Xu Q, Zhang W, Song Q, et al. Novel reassortant highly pathogenic avian influenza (H5N5) viruses in domestic ducks, China. Emerg Infect Dis. 2011;17:1060–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Ip HS, Torchetti MK, Crespo R, Kohrs P, DeBruyn P, Mansfield KG, Baszler T, Badcoe L, Bodenstein B, Shearn-Bochsler V, et al. Novel Eurasian highly pathogenic avian influenza A H5 viruses in wild birds, Washington, USA, 2014. Emerg Infect Dis. 2015;21:886–90.View ArticlePubMedPubMed CentralGoogle Scholar
- Jeong J, Kang HM, Lee EK, Song BM, Kwon YK, Kim HR, Choi KS, Kim JY, Lee HJ, Moon OK, et al. Highly pathogenic avian influenza virus (H5N8) in domestic poultry and its relationship with migratory birds in South Korea during 2014. Vet Microbiol. 2014;173:249–57.View ArticlePubMedGoogle Scholar
- Verhagen JH, van der Jeugd HP, Nolet BA, Slaterus R, Kharitonov SP, de Vries PP, Vuong O, Majoor F, Kuiken T, Fouchier RA. Wild bird surveillance around outbreaks of highly pathogenic avian influenza A(H5N8) virus in the Netherlands, 2014, within the context of global flyways. Euro Surveill. 2015;20.Google Scholar
- Highly Pathogenic H5 Avian Influenza Confirmed in Wild Birds in Washington State H5N2 Found in Northern Pintail Ducks & H5N8 Found in Captive Gyrfalcons. http://www.usda.gov/wps/portal/usda/usdahome?contentid=2014/12/0273.xml
- Highly Pathogenic H5N8 Avian Influenza Confirmed in Commercial Turkey Flock in California. https://www.aphis.usda.gov/aphis/newsroom/news/SA_By_Date/SA_2015/SA_01/CT_hpai_california
- Avian Influenza Disease. http://www.aphis.usda.gov/wps/portal/aphis/ourfocus/animalhealth/sa_animal_disease_information/sa_avian_health/
- Spackman E, Senne DA, Bulaga LL, Myers TJ, Perdue ML, Garber LP, Lohman K, Daum LT, Suarez DL. Development of real-time RT-PCR for the detection of avian influenza virus. Avian Dis. 2003;47:1079–82.View ArticlePubMedGoogle Scholar
- Zhou B, Donnelly ME, Scholes DT, St George K, Hatta M, Kawaoka Y, Wentworth DE. Single-reaction genomic amplification accelerates sequencing and vaccine production for classical and Swine origin human influenza a viruses. J Virol. 2009;83:10309–13.View ArticlePubMedPubMed CentralGoogle Scholar
- Edgar RC. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics. 2004;5:113.View ArticlePubMedPubMed CentralGoogle Scholar
- Drummond AJ, Nicholls GK, Rodrigo AG, Solomon W. Estimating mutation parameters, population history and genealogy simultaneously from temporally spaced sequence data. Genetics. 2002;161:1307–20.PubMedPubMed CentralGoogle Scholar
- Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73.View ArticlePubMedPubMed CentralGoogle Scholar
- Tracer v1.6. http://beast.bio.ed.ac.uk/Tracer
- Ozawa M, Matsuu A, Tokorozaki K, Horie M, Masatani T, Nakagawa H, Okuya K, Kawabata T, Toda S: Genetic diversity of highly pathogenic H5N8 avian influenza viruses at a single overwintering site of migratory birds in Japan, 2014/15. Euro Surveill. 2015;20.Google Scholar
- Influenza A Cleavage Sites. http://www.offlu.net/fileadmin/home/en/resource-centre/pdf/Influenza_A_Cleavage_Sites_Aug2015.pdf
- Cattoli G, Milani A, Temperton N, Zecchin B, Buratin A, Molesti E, Aly MM, Arafa A, Capua I. Antigenic drift in H5N1 avian influenza virus in poultry is driven by mutations in major antigenic sites of the hemagglutinin molecule analogous to those for human influenza virus. J Virol. 2011;85:8718–24.View ArticlePubMedPubMed CentralGoogle Scholar
- Eshaghi A, Shalhoub S, Rosenfeld P, Li A, Higgins RR, Stogios PJ, Savchenko A, Bastien N, Li Y, Rotstein C, Gubbay JB. Multiple influenza A (H3N2) mutations conferring resistance to neuraminidase inhibitors in a bone marrow transplant recipient. Antimicrob Agents Chemother. 2014;58:7188–97.View ArticlePubMedPubMed CentralGoogle Scholar
- Dalby AR, Iqbal M. The European and Japanese outbreaks of H5N8 derive from a single source population providing evidence for the dispersal along the long distance bird migratory flyways. Peer J. 2015;3:e934.View ArticlePubMedPubMed CentralGoogle Scholar
- Hill SC, Lee YJ, Song BM, Kang HM, Lee EK, Hanna A, Gilbert M, Brown IH, Pybus OG. Wild waterfowl migration and domestic duck density shape the epidemiology of highly pathogenic H5N8 influenza in the Republic of Korea. Infect Genet Evol. 2015;34:267–77.View ArticlePubMedPubMed CentralGoogle Scholar
- Bamford M, Watkins D, Bancroft W, Tischler G, Wahl J. Migratory Shorebirds of the East Asian - Australian Flyway: Population Estimates and Internationally Important Sites, Wetlands International - Oceania. 2008.Google Scholar
- Ramey AM, Pearce JM, Reeves AB, Poulson RL, Dobson J, Lefferts B, Spragens K, Stallknecht DE. Surveillance for Eurasian-origin and intercontinental reassortant highly pathogenic influenza A viruses in Alaska, spring and summer 2015. Virol J. 2016;13:55.View ArticlePubMedPubMed CentralGoogle Scholar
- Torchetti MK, Killian ML, Dusek RJ, Pedersen JC, Hines N, Bodenstein B, White CL, Ip HS. Novel H5 Clade 126.96.36.199 Reassortant (H5N1) Virus from a Green-Winged Teal in Washington, USA. Genome Announc. 2015;3:e00195-15.Google Scholar
- Vijaykrishna D, Bahl J, Riley S, Duan L, Zhang JX, Chen H, Peiris JS, Smith GJ, Guan Y. Evolutionary dynamics and emergence of panzootic H5N1 influenza viruses. PLoS Pathog. 2008;4:e1000161.View ArticlePubMedPubMed CentralGoogle Scholar
- New CDC Laboratory Study Suggests U.S. H5 Bird Flu Viruses Currently Pose Low Risk to People. http://www.cdc.gov/flu/news/bird-viruses-risk.htm
- Pan M, Gao R, Lv Q, Huang S, Zhou Z, Yang L, Li X, Zhao X, Zou X, Tong W, et al. Human infection with a novel highly pathogenic avian influenza A (H5N6) virus: Virological and clinical findings. J Infect. 2016;72(1):52–9.View ArticlePubMedGoogle Scholar