A novel Clade 184.108.40.206c H5N1 reassortant virus caused several outbreaks in wild birds in some regions of China from late 2014 to 2015. Based on the genetic and phylogenetic analyses, the viruses possess a stable gene constellation with a Clade 220.127.116.11c HA, a H9N2-derived PB2 gene and the other six genes of Asian H5N1-origin. The Clade 18.104.22.168c H5N1 reassortants displayed a high genetic relationship to a human H5N1 strain (A/Alberta/01/2014). Further analysis showed that similar viruses have been circulating in wild birds in China, Russia, Dubai (Western Asia), Bulgaria and Romania (Europe), as well as domestic poultry in some regions of Africa. The affected areas include the Central Asian, East Asian-Australasian, West Asian-East African, and Black Sea/Mediterranean flyways. These results show that the novel Clade 22.214.171.124c reassortant viruses are circulating worldwide and may have gained a selective advantage in migratory birds, thus posing a serious threat to wild birds and potentially humans.
Citation: Yuhai Bi, Jianjun Chen, Zhenjie Zhang, Mingxin Li, Tianlong Cai, Kirill Sharshov, Ivan Susloparov, Alexander Shestopalov, Gary Wong, Yubang He, Zhi Xing, Jianqing Sun, Di Liu, Yingxia Liu, Lei Liu, Wenjun Liu, Fumin Lei, Weifeng Shi, George F. Gao. Highly pathogenic avian influenza H5N1 Clade 126.96.36.199c virus in migratory birds, 2014-2015[J]. VIROLOGICA SINICA, 2016, 31 (4): 300-305 https://doi.org/10.1007/s12250-016-3750-4
Received: 22 February, 2016; Accepted: 8 June 2016; Published: 30 June 2016
Copyright: © Wuhan Institute of Virology, CAS and Springer Science+Business Media Singapore 2016
Data Availability: All relevant data are within the paper and its Supporting Information files.
# These authors contributed equally to this article.
Highly pathogenic avian influenza (HPAI) H5N1 virus was initially isolated from China in 1996, and termed A/Goose/Guangdong/1/96(Gs/Gd). The descendants of Gs/Gd caused much concern during 1997, when an H5N1 virus outbreak in the Hong Kong Special Admi-nistration Region (HKSAR) of China resulted in 18 human infections and 6 deaths (Webster et al., 2006). Re-emergence of H5N1 was then reported during 2003 and continuing occurrences have been recorded ever since, including outbreaks in migratory birds (Liu et al., 2005; Wang et al., 2008; Hu et al., 2011; Sakoda et al., 2012; Naguib et al., 2015). As of 14 December 2015, transmission of Gs/Gd-lineage H5N1 viruses from avian hosts to humans has been reported in 16 countries worldwide, resulting in at least 844 laboratory-confirmed infections and 449 deaths (http://www.who.int/influenza/human_animal_interface/H5N1_cumulative_table_archives/en/). H5N1 viruses are now considered endemic amongst domestic birds in China, Southeast Asia and Africa, and have evolved into numerous Clades and sub-Clades (World Health Organization, 2014). In particular, Clade 188.8.131.52 H5N1 viruses were identified for the first time from a dead Chinese pond heron in HKSAR during 2004, and circulated in wild birds during 2009~2011(Hu et al., 2011; Sakoda et al., 2012). Since 2009, Clade 184.108.40.206 viruses have diversified into sub-Clades 220.127.116.11a, 18.104.22.168b and 22.214.171.124c (WHO, 2014; Bi et al., 2015a).
During 2014-2015, a novel Clade 126.96.36.199c H5N1 vi-rus caused several outbreaks amongst wild or domestic birds in Russia, China (Bi et al., 2015a), The United Arab Emirates (Western Asia) (Naguib et al., 2015), Bulgaria and Romania (Europe), and several African countries (FAO, 2015). In China, there have been at least three outbreaks in Sanmenxia Reservior, Inner Mongolia and Qinghai Lake, resulting in the deaths of thousands of wild birds. Therefore, a systematic analysis of the genetic evolution and relationships between these globally circulating viruses will help efforts to understand, prevent and control the spread of HPAI H5N1 virus.
Virus isolation and gene sequencing
Oropharyngeal and cloacal swabs, tissues and organs from sick or dead birds, and feces from apparently healthy birds were collected for virus isolation. Swabs and tissues were maintained at 4 ℃ in viral-transport medium and sample tubes, respectively. Before inoculation, swabs with viral-transport medium were thoroughly mixed. The tissues and organs were homoge-nized using a Qiagen TissueLyser Ⅱ machine (30 cycles/s; 4 min) in 1-mL cold PBS under sterile conditions. The solid debris in the swab and tissue samples was pelleted by centrifugation at 5, 000 ×g for 10 min. The supernatants were then inoculated in 10-day-old chicken embryos for 72 hours at 37 ℃. Viral RNA from HA titer-positive allantoic fluid was extracted using TRIZOL (Invitrogen) according to manufacturer instructions. All gene segments were amplified using Ex TaqTM DNA polymerase (Takara) with the improved primers (Bi et al., 2015b), and sequenced by the ABI 3730XL automatic DNA analyzer.
Genetic and phylogenetic analysis
The gene and amino acid sequences were assembled using DNAMAN (Version 7.0) and BioEdit (Version 188.8.131.52). Representative reference sequences were downloaded from the NCBI (http://www.ncbi.nlm.nih.gov) and GISAID (www.gisaid.org) databases via the online Blast search. Multiple sequence alignment was performed using Muscle (Edgar, 2004). Maximum-likelihood phylogenetic trees were inferred using the RAxML software under the GTRGAMMA model with 1000 bootstrap replicates (Stamatakis, 2014).
The GISAID accession numbers for the novel genome sequences of the H5N1 viruses are EPI669753-60, EPI669761-68, EPI669769-76, EPI669777-84, EPI669785-92, EPI669793-800, EPI669801-08, EPI669809-16, EPI669817-24, EPI669825-32, EPI669833-40, EPI669841-48 and EPI669849-56.
Disease outbreak, virus isolation and gene sequencing
In May 2015, approximately one hundred black-necked grebes (Podiceps nigricollis) were found dead in Inner Mongolia (Qehan Lake), China. Clinical symptoms, including neurological signs (tremor and opisthotonus) as well as diarrhea were observed in these birds. Various organs, oropharynx and rectal swabs were collected from four diseased or dead birds. Influenza A (H5N1) virus was detected and isolated from these samples, and a total of ten full-length genome sequences of the viruses were obtained.
On 8 July 2015, several great black-headed gulls (Larusi-chthyaetus) were found dead in a village near Qinghai Lake, located approximately 1250 km from Qehan Lake. As of 14 July, a total of more than 2300 birds were found dead in Qinghai Lake, where an HPAI H5N1 virus outbreak had resulted in the deaths of thousands of migra-tory birds in 2005(Liu et al., 2005). The majority of the affected birds in the 2015 outbreak were great black-headed gulls, and these birds were found on a single islet or river estuary in Qinghai Lake. Neurological symptoms were also observed in the infected birds. On 14 July, intestine samples were collected from 11 dead birds, including 8 great black-headed gulls, 2 bar-headed geese (Anser indicus) and 1 great cormorant (Phalacrocorax carbo). Three of the birds were detected to be positive for influenza A (H5N1) virus, and full-length genome sequences were obtained.
The most closely related sequences for the 13 novel ge-nome sequences were obtained from the GISAID database using BLAST. In addition, the reference sequences for the H5 subtype of HPAI were included for genetic grouping (World Health Organization, 2014). All of the sequences were aligned and phylogenetically analyzed. The HA phylogeny revealed that all of the novel sequences belonged to Clade 184.108.40.206c (Figure 1) and were closely related to sequences isolated from Sanmenxia, China, where HPAI H5N1 viruses caused the deaths of approximately one hundred migratory birds (mostly whooper swans) between late 2014 and early 2015(Bi et al., 2015a). Notably, the Chinese strains also clustered together with those from Russia, West Asia, Europe and Africa that were isolated during May 2015, winter of 2014(Naguib et al., 2015), and spring of 2015(Monne et al., 2015), respectively (Figure 1). For the remaining seven genes, sequences from the strains mentioned above were clustered together and formed an independent sub-Clade (Figure 2 and Supplementary Figures S1-S6), indicating that they belonged to the same genotype. This genotype consisted of a Clade 220.127.116.11c HA (Figure 1) and a PB2 gene derived from H9N2(Figure 2), with the remaining genes originating from other HPAI H5N1 viruses (Bi et al., 2015a; Monne et al., 2015). It should be noted that the prototype virus of this genotype (A/Alberta/01/2014) was first identified from a Canadian patient who traveled to Beijing, China in December 2013(Pabbaraju et al., 2014). However, Clade 18.104.22.168c viruses were subsequently identified from wild birds in the central regions of China, the western regions of Russia, Western Asia and Africa along the Central Asian flyway (Figure 3), suggesting that the novel Clade 22.214.171.124c reassortant viruses are being spread by wild birds along migratory flyways.
The amino acid motif at the HA cleavage site of the H5N1 isolates is RERRRKRG (Supplementary Table S1), indicating that they are highly pathogenic to chickens. In addition, they possess Q226 and G228(H3 numbering) at the receptor binding site, suggesting that they favor avian-like receptors (Stevens et al., 2006). However, they also contain the amino acid substitution T160A, which is thought to enhance the binding capacity to human-like receptors (Zhang et al., 2013). The amino acid deletion at the NA stalk (49-68) is considered a molecular marker for avian influenza virus from aquatic birds for adaptation to land poultry (Matsuoka et al., 2009; Zhou et al., 2009; Cauldwell et al., 2014). Therefore, the Clade 126.96.36.199c H5N1 viruses may have already adapted to poultry before their introduction into migratory birds. Amino acids Q591(Mok et al., 2011), E627(Subbarao et al., 1993; Hatta et al., 2001) and D701(Zhou et al., 2013) in the PB2 protein of the H5N1 viruses suggest that they have not yet adapted to mammalian hosts. No drug-resistance-associated mutations (H274Y in NA; S31N in M2(Hay et al., 1985; Pinto et al., 1992)) were observed among migratory H5N1 viruses, suggesting that these H5N1 isolates are still sensitive to NA and M2 inhibitors.
Considering that influenza A viruses undergo frequent genetic reassortment, the relatively stable genetic constellation of this genotype during its circulation over long periods of time and distances imply that this novel HPAI H5N1 virus might have gained certain selective advantages that allow for the stable passaging of its viral ge-nome. Clade 188.8.131.52c viruses were identified from wild birds in the central regions of China and the western regions of Russia along the Central Asian flyway (Figure 3). This suggests that the novel Clade 184.108.40.206c reassortant viruses are being spread by wild birds along migra-tory flyways. Therefore, countries located along the East Asian-Australian and Central Asian flyways should be aware of the potential introduction of this reassortant vi-rus, as migratory birds carrying Clade 220.127.116.11 H5N1 vi-ruses have already struck Japan and South Korea in the past (Sakoda et al., 2012). Additionally, other countries should also be on the alert for the potential introduction of this virus into their territories, as exemplified by the introduction of Clade 18.104.22.168 H5N1 viruses into North America and Europe during late 2014 (Jhung and Nelson, 2015; Verhagen et al., 2015). A successful example for outbreak control was the H5N1 outbreak amongst whooper swans in Sanmenxia during January 2015, in which high risk regions for virus transmission along the flyways (such as inner Mongolia) were alerted immediately (Bi et al., 2015a), which helped their preparations and in turn played a considerable role in controlling the outbreak. In addition, human infections with Clade 22.214.171.124c viruses constitute an important reminder of the public health risks from migratory birds carrying H5N1.
This study was supported by grants from the National Natural Science Foundation of China (31311120063, 81470096, 31570026, 31471253), the Ministry of Education and Science of the Russian Federation (Chinese-Russian project: RFMEFI61315X0045), the intramural special grant for influenza virus research from the Chinese Academy of Sciences (KJZD-EW-L09), the IDRC-APEIR program (106915-001), and Special Project of Ministry of Science and Technology (2013FY113500). GFG is a leading principal investigator of the NSFC Innovative Research Group (81321063). WFS was supported by the "Taishan Scholar" project of Shandong Province. G.W. is the recipient of a Banting Postdoctoral Fellowship from the Canadian Institutes of Health Research (CIHR) and the President's International Fellowship Initiative from the Chinese Academy of Sciences (CAS). We thank the data submitters from the GISAID and GenBank Flu databases for the AIV sequences.
COMPLIANCE WITH ETHICS GUIDELINES
All the authors declare that they have no conflict of interests. This article does not contain any studies with human or animal subjects performed by any of the authors.
YHB, JJC, WFS and GFG conceived and designed the research. YHB, JJC and WFS performed the experiments. ZJZ, MXL, KS, IS, AS, TLC, GW and YBH performed sample collections. DL, WJL and FML provided helpful suggestions about the study.
Supplementary table/figures are available on the websites of Virologica Sinica: www.virosin.org; link.springer.com/journal/12250.
Table S1. Comparison of key amino acid pattern of Clade 126.96.36.199c HPAI A (H5N1) viruses from different geographic regions
- . Bi Y, Zhang Z, Liu W, Yin Y, Hong J, Li X, Wang H, Wong G, Chen J, Li Y, Ru W, Gao R, Liu D, Liu Y, Zhou B, Gao GF, Shi W, Lei F. 2015a. Highly Pathogenic Avian Influenza A (H5N1) Virus Struck Migratory Birds in China in 2015. Sci Rep, 5: 12986.
- . Bi YH, Mei K, Shi WF, Liu D, Yu XL, Gao ZM, Zhao LH, Gao GF, Chen JJ, Chen QJ. 2015b. Two novel reassortants of avian influenza A (H5N6) virus in China. J Gen Virol, 96: 975-981.
- . Cauldwell AV, Long JS, Moncorgé O, Barclay WS. 2014. Viral determinants of influenza A virus host range. J Gen Virol, 95: 1193-1210.
- . Edgar RC. 2004. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res, 32: 1792-1797.
- . FAO. 2015. H5N1 HPAI spread in Nigeria and increased risk for neighbouring countries in West Africa. http://www.fao.org/3/a-i4561e.pdf.
- . Hatta M, Gao P, Halfmann P, Kawaoka Y. 2001. Molecular basis for high virulence of Hong Kong H5N1 influenza A viruses. Science, 293: 1840-1842.
- . Hay AJ, Wolstenholme AJ, Skehel JJ, Smith MH. 1985. The molecular basis of the specific anti-influenza action of amantadine. The EMBO J, 4: 3021-3024.
- . Hu X, Liu D, Wang M, Yang L, Zhu Q, Li L, Gao GF. 2011. Clade 2.3.2 avian influenza virus (H5N1), Qinghai Lake region, China, 2009-2010. Emerg Infect Dis, 17: 560-562.
- . Jhung MA, Nelson DI. 2015. Outbreaks of avian influenza A (H5N2), (H5N8), and (H5N1) among birds—United States, December 2014-January 2015. MMWR (Morb Mortal Wkly Rep), 64: 111.
- . Liu J, Xiao H, Lei F, Zhu Q, Qin K, Zhang XW, Zhang XL, Zhao D, Wang G, Feng Y, Ma J, Liu W, Wang J, Gao GF. 2005. Highly pathogenic H5N1 influenza virus infection in migratory birds. Science, 309: 1206.
- . Matsuoka Y, Swayne DE, Thomas C, Rameix-Welti M-A, Naffakh N, Warnes C, Altholtz M, Donis R, Subbarao K. 2009. Neuraminidase stalk length and additional glycosylation of the hemagglutinin influence the virulence of influenza H5N1 vi-ruses for mice. J Virol, 83: 4704-4708.
- . Mok CKP, Yen HL, Yu MYM, Yuen KM, Sia SF, Chan MCW, Qin G, Tu WW, Peiris JSM. 2011. Amino acid residues 253 and 591 of the PB2 protein of avian influenza virus A H9N2 contribute to mammalian pathogenesis. J Virol, 85: 9641-9645.
- . Monne I, Meseko C, Joannis T, Shittu I, Ahmed M, Tassoni L, Fusaro A, Cattoli G. 2015. Highly pathogenic avian influenza A (H5N1) Virus in poultry, Nigeria, 2015. Emerg Infect Dis, 21: 1275-1277.
- . Naguib MM, Kinne J, Chen H, Chan KH, Joseph S, Wong PC, Woo PC, Wernery R, Beer M, Wernery U, Harder TC. 2015. Outbreaks of highly pathogenic avian influenza H5N1 Clade 188.8.131.52c in hunting falcons and kept wild birds in Dubai implicate intercontinental virus spread. J Gen Virol, 96: 3212-3212.
- . Pabbaraju K, Tellier R, Wong S, Li Y, Bastien N, Tang JW, Drews SJ, Jang Y, Davis CT, Fonseca K, Tipples GA. 2014. Full-genome analysis of avian influenza A (H5N1) virus from a human, North America, 2013. Emerg Infect Dis, 20: 887-891.
- . Pinto LH, Holsinger LJ, Lamb RA. 1992. Influenza virus M2 protein has ion channel activity. Cell, 69: 517-528.
- . Sakoda Y, Ito H, Uchida Y, Okamatsu M, Yamamoto N, Soda K, Nomura N, Kuribayashi S, Shichinohe S, Sunden Y, Umemura T, Usui T, Ozaki H, Yamaguchi T, Murase T, Ito T, Saito T, Takada A, Kida H. 2012. Reintroduction of H5N1 highly pathogenic avian influenza virus by migratory water birds, causing poultry outbreaks in the 2010-2011 winter season in Japan. J Gen Virol, 93: 541-550.
- . Stamatakis A. 2014. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics, 30: 1312-1313.
- . Stevens J, Blixt O, Tumpey TM, Taubenberger JK, Paulson JC, Wilson IA. 2006. Structure and receptor specificity of the hemagglutinin from an H5N1 influenza virus. Science, 312: 404-410.
- . Subbarao EK, London W, Murphy BR. 1993. A single amino-acid in the Pb2-Gene of influenza a virus is a determinant of host range. J Virol, 67: 1761-1764.
- . Verhagen JH, Herfst S, Fouchier RAM. 2015. How a virus travels the world. Science, 347: 616-617.
- . Wang G, Zhan D, Li L, Lei F, Liu B, Liu D, Xiao H, Feng Y, Li J, Yang B, Yin Z, Song X, Zhu X, Cong Y, Pu J, Wang J, Liu J, Gao GF, Zhu Q. 2008. H5N1 avian influenza re-emergence of Lake Qinghai: phylogenetic and antigenic analyses of the newly isolated viruses and roles of migratory birds in virus circulation. J Gen Virol, 89: 697-702.
- . Webster RG, Peiris M, Chen H, Guan Y. 2006. H5N1 outbreaks and enzootic influenza. Emerg Infect Dis, 12: 3-8.
- . World Health Organization (WHO). 2014. Revised and updated nomenclature for highly pathogenic avian influenza A (H5N1) viruses. Influenza Other Respir Viruses, 8: 384-388.
- . Zhang W, Shi Y, Lu XS, Shu YL, Qi JX, Gao GF. 2013. An airborne transmissible avian influenza H5 hemagglutinin seen at the atomic level. Science, 340: 1463-1467.
- . Zhou B, Pearce MB, Li Y, Wang JR, Mason RJ, Tumpey TM, Wentworth DE. 2013. Asparagine substitution at PB2 residue 701 enhances the replication, pathogenicity, and transmission of the 2009 pandemic H1N1 influenza A virus. PloS One: 8.
- . Zhou H, Yu Z, Hu Y, Tu J, Zou W, Peng Y, Zhu J, Li Y, Zhang A, Yu Z. 2009. The special neuraminidase stalk-motif responsible for increased virulence and pathogenesis of H5N1 influenza A virus. PLoS One, 4: e6277.