Disclaimer: Early release articles are not considered as final versions. Any changes will be reflected in the online version in the month the article is officially released.
Volume 28, Number 7—July 2022
Epidemiologic, Clinical, and Genetic Characteristics of Human Infections with Influenza A(H5N6) Viruses, China
The recent rise in the frequency of influenza A(H5N6) infections in China has raised serious concerns about whether the risk for human infection has increased. We surveyed epidemiologic, clinical, and genetic data of human infections with A(H5N6) viruses. Severe disease occurred in 93.8% of cases, and the fatality rate was 55.4%. Median patient age was 51 years. Most H5N6 hemagglutinin (HA) genes in human isolates in 2021 originated from subclade 22.214.171.124b; we estimated the time to most recent common ancestor as June 16, 2020. A total of 13 genotypes with HA genes from multiple subclades in clade 126.96.36.199 were identified in human isolates. Of note, 4 new genotypes detected in 2021 were the major causes of increased H5N6 virus infections. Mammalian-adapted mutations were found in HA and internal genes. Although we found no evidence of human-to-human transmission, continuous evolution of H5N6 viruses may increase the risk for human infections.
Various subtypes of avian influenza viruses (AIVs) circulate globally in wild birds. Domestic poultry are also susceptible to these viruses and, occasionally, AIVs have become able to cross the species barrier to infect humans (1–4). In China, situations in which human infections with AIVs have occurred can be complicated; many subtypes of AIVs have been reported, including influenza A(H5N1), A(H5N6), A(H7N4), A(H7N9), A(H9N2), A(H10N8), and A(H10N3) (4,5). Highly pathogenic avian influenza (HPAI) A(H5) viruses have continually caused worldwide outbreaks in both wild birds and poultry, with some spillover to humans, most notably 863 HPAI A(H5N1) cases, 456 of which were fatal (6).
In April 2014, the first human infection with HPAI A(H5N6) virus was reported in Sichuan Province, China (1). Since then, human cases have been continuously documented in China. By the end of 2021, 66 cases had been documented globally, 36 of which were fatal. Of those, 65 cases were in China, and the remaining case was reported in Laos.
Since its emergence in Guangdong Province in 1996, HPAI H5 type viruses became endemic in birds in China and other regions and developed into distinct clades (7). The clade 188.8.131.52 H5 viruses were first reported in migratory birds in eastern China in 2013, followed by outbreaks in both wild and domestic birds in South Korea at the beginning of 2014. Thereafter, clade 184.108.40.206 H5 viruses spread westward and eastward from Asia to other continents, accompanied by multiple reassortment events between clade 220.127.116.11 H5 viruses and other AIVs circulating in wild and domestic birds (8,9). Clade 18.104.22.168 H5N6 AIV has been mainly endemic among birds in China and southeast Asia since 2013 (10–12), gradually replacing H5N1 as a dominant AIV subtype in poultry across southern China (13). Similar to H7N9 AIV, reassortments among H5N6 and H9N2 AIVs occurred dynamically (14–16). Furthermore, many other subtypes of AIVs were documented to donate their internal genes to H5N6, including AIVs from wild birds (11,17,18). Thus, H5N6 genotypes diversified both in poultry and humans (11,16–18).
Genetic and biologic characteristics indicated that H5N6 viruses were highly pathogenic in chickens (19–21). Studies have shown that some H5N6 viruses of clade 22.214.171.124 possessed the ability to bind both avian-origin and human-origin sialic acid receptors and had the ability to attach to human tracheal epithelial and alveolar tissues (22). Although H5N6 viruses were not as pathogenic in mice and ferrets as their parental clade 2.3.4 H5N1 viruses, they exhibited greater transmissibility than H5N1 viruses in a ferret model (22,23). Thus, clade 126.96.36.199 H5N6 viruses demonstrated higher potential to transmit among humans.
During the COVID-19 pandemic, seasonal influenza infections notably decreased worldwide compared with previous flu seasons (24). However, human infections with zoonotic influenza viruses do not appear to have decreased. Of note, a recent rise in the frequency of H5N6 cases was observed (25). For the purposes of preparing for future possible pandemics, investigating the epidemiologic, clinical, and genetic characteristics of H5N6 AIVs that infect humans is critical.
In China, all laboratory-confirmed influenza A(H5N6) cases are reported through a national surveillance system. Patients with respiratory tract specimens that test positive for H5 and N6 by real-time reverse transcription PCR are confirmed as H5N6 infections. Demographic, epidemiologic, and basic clinical data on influenza A(H5N6) cases are collected on standardized forms. We included information regarding age; sex; place of residence; symptoms at illness onset; comorbidities associated with increased risk for influenza complications; dates of illness onset, hospital admission, death or discharge, and clinical treatments; and potential exposures to domestic or retail animals and visits to live poultry markets in our analysis.
The National Health Commission of the People’s Republic of China determined that the collection of data on each H5N6 case was part of a continuing public health investigation of an emerging outbreak. Thus, the ethics approval was exempt from official review by our institutional review board.
We used original samples from human cases for H5/N6 subtyping by real-time reverse transcription PCR (26). We selected positive samples for virus isolation. We conducted the isolation of H5N6 viruses in a BioSafety Level 3 laboratory by inoculating 0.2 mL of original sample into 9- to 11-day specific pathogen–free embryonated chicken eggs. After samples incubated for 36–48 hours at 37°C, we harvested allantoic fluids.
RNA Extraction and Genome Sequencing
We extracted RNA from the original samples and isolated viruses by using the QIAamp viral RNA mini kit (QIAGEN, https://www.qiagen.com) and performed sequencing as follows. As previously described (27), we subjected extracted RNA to reverse transcription and amplification. We then implemented whole-genome sequencing of FluA on the MiSeq or Miniseq high-throughput sequencing platform (Illumina, https://www.illumina.com) (28). We predominately conducted data analysis and genome sequence acquisition by using a pipeline established in our laboratory. We trimmed low-quality reads, sampled the filtered reads, and de novo assembled using Velvet version 1.2.10 (https://guix.gnu.org/en/packages/velvet-1.2.10) and Newbler Assembler version 2.5. We blasted contigs against a database generated by CD-HIT that clusters all FluA sequences collected from the GISAID EpiFlu database (http://www.gisaid.org) and National Center for Biotechnology Information (NCBI) Influenza Virus Database. We selected sequences with the highest similarity as references and used bowtie2 version 2.1.0 (https://sourceforge.net/projects/bowtie-bio/files/bowtie2/2.1.0) for read mapping. We obtained FluA genome sequences by extracting the consensus sequences from the mapping results, with a coverage depth of at least 30× at each site on the 8 segments. We submitted the genome sequences of influenza A(H5N6) viruses determined in this study to GISAID (Appendix Table 1).
Sequence Alignment and Phylogenetic Analysis
We collected all H5N6 human isolate sequences, including sequences from GISAID, and downloaded the top 100 sequences of avian viruses with high similarity to each representative sequence through BLAST (https://blast.ncbi.nlm.nih.gov/Blast.cgi) from GISAID. We used CD-HIT to reduce the sequence redundancy of phylogenetic analysis and performed sequence alignments by using MAFFT software version 6.857b (https://mafft.cbrc.jp/alignment/software). We constructed a maximum-likelihood phylogenetic tree for the nucleotide sequences of each gene of selected influenza viruses under the GTRGAMMAR model with 1,000 bootstrap replicates, using RAxML (29). To estimate the time to most recent common ancestor (tMRCA) of hemagglutinin (HA) and neuraminidase (NA) genes of the H5N6 virus, we selected nonredundant subdatasets to run time-measured Bayesian Markov chain Monte Carlo analysis by BEAST v1.10.4 (30). We used the SRD06 substitution model (31) and the uncorrelated relaxed molecular clock model and set the Bayesian skygrid coalescent as the tree prior. We ran the Bayesian Markov chain Monte Carlo for up to 1 × 108 steps, with samples for each 10,000 steps to achieve convergence. We used tracer version 1.6 (https://bioweb.pasteur.fr/packages/pack@Tracer@v1.6) to examine effective sample size values >200.
We summarized continuous variables as either means + SD or medians with interquartile ranges (IQRs). For categorical variables, we calculated the percentages of patients in each category. We used parametric tests to analyze normally distributed variables and nonparametric tests to analyze non–normally distributed variables. We used an unpaired Student t test, Wilcoxon rank-sum test, χ2 test, or Fisher exact test, as appropriate, to compare the epidemiologic and clinical characteristics of subgroups of patients who were infected with influenza A(H5N6) virus before or during 2021. We calculated 95% CIs for the means of normally distributed data and the risk estimates by using the binomial distribution. We considered a p value of <0.05 statistically significant. We performed all analyses in SAS version 9.4 (SAS Institute, Inc., https://www.sas.com).
During April 21, 2014–December 31, 2021, a total of 65 cases of human influenza A(H5N6) infection were reported in China; illness onset dates occurred during April 13, 2014–December 31, 2021. The case-fatality rate (CFR) was 55.4% (36/65). Most cases were reported through pneumonia surveillance and were identified by Chinese National Influenza Surveillance Network laboratories; however, 15 cases in 2020 and 2021 were first identified through third-party sequencing agencies and reported to Chinese National Influenza Surveillance Network laboratories for confirmation.
The number of human cases of H5N6 infection reported for each year from 2014 to 2021 was 2 for 2014, 6 for 2015, 9 for 2016, 2 for 2017, 4 for 2018, 1 for 2019, 5 for 2020, and 36 for 2021 (Figure 1). These cases were distributed across 14 provinces in China; most (64/65) cases were detected in southern China (Figure 2). The H5N6 cases were predominantly in adults; median patient age was 51 (IQR 36–57) years (Table 1). More than half (45/65, 69.2%) of H5N6 cases occurred in persons 18–59 years of age. Men accounted for 50.8% (33/65) of H5N6 infections (Table 1).
Poultry exposure is the main risk factor for human infection with AIVs. Nearly half (31/65) of the H5N6 cases occurred in residents of rural areas, where birds were raised. Among the 65 cases, information regarding history of poultry exposure was available for 61 (93.8%) persons; the most common exposures were visiting a live poultry market (38, 62.3%) and exposure to backyard poultry (28, 45.9%). Two cases reported in August 2021 were in a husband and wife, both of whom reported poultry exposure. All close contacts were investigated, and no influenza-like illness symptoms were observed for any household contacts. However, human-to-human transmission cannot be ruled out for the infections in the husband and wife.
Of the 65 case-patients, 31 reported comorbidities (Table 1), 24 were otherwise healthy, and 10 did not report information on other health conditions. The prevalence of coronary heart disease and cancer in influenza H5N6 case-patients was higher than in the general population in China (Appendix Table 2).
The most commonly reported symptoms at onset of illness were fever (73.4%) and cough (59.4%) (Table 2), indicating that H5N6 infection is difficult to distinguish from other respiratory illnesses, including seasonal influenza viruses, in terms of early symptoms. The median onset-to-admission interval was 4 days (IQR 2–6 days), and the median onset-to-laboratory confirmation delay was 9 days (IQR 7–12 days).
Antiviral treatments targeting NA were given to 60% of patients. However, the median time from illness onset to antiviral treatment was 5 days (IQR 3–9 days) (Figure 3). For patients in whom time from illness onset to initiation of antiviral treatment was <5 days, CFR was 38.1% (8/21), lower than in patients for whom that interval was >5 days (66.7% [12/18]). For case-patients who received no antiviral treatment, CFR was 65.2% (15/23). Although not statistically significant (p = 0.13), differences did exist in the CFR value among these groups, indicating that early initiation of antiviral treatments could increase the survival rate of H5N6 patients.
Severe disease occurred in 93.8% of cases, and complications developed in 86.2% of patients (Table 1). A total of 51 case-patients required admission to the intensive care unit for treatment; 41 persons underwent mechanical ventilation, and 14 persons required extracorporeal membrane oxygenation (Table 1). We estimated that the hospitalization fatality risks were 59.0% (95% CI 46.7%–71.4%). No significant difference was observed in terms of mean age between fatal cases (45.2 + 19.9 years) and survivors (45.1 + 20.2 years).
Comparative Epidemiology of H5N6 Infections
To identify any possible changes to risk for H5N6 infection before and during 2021 in China, we comparatively analyzed epidemiologic and clinical characteristics of human infections. We observed no significant difference in fatality rate, sex distribution, or the course of disease (Table 3; Appendix Figure 1). However, differences between the seasonal distributions of H5N6 infections were evident (Figure 1). In previous years, few cases were reported during June–September, but an increase in H5N6 infections occurred during summer 2021. In addition, we observed a marked difference in case-patients’ geographic location in 2021 (Table 3). Approximately 34.5% of H5N6 case-patients before 2021 lived in rural areas; during 2021, ≈58.3% of case-patients lived in rural areas (p<0.05). Furthermore, the average age of patients who were infected with or died of H5N6 viruses was much lower before 2021 than during 2021. The median age of those infected with H5N6 viruses was 40 years (IQR 25–50 years) before 2021 versus 54 years (IQR 49.5–60.5 years) during 2021 (p<0.001) (Table 3). The mean age of persons who died of H5N6 infection was 37.7 years (95% CI 27.9–47.4, SD 19.6 years) before 2021 versus 52.8 years (95% CI 44.0–61.6,; SD 17.7 years) during 2021 (p<0.05). Of note, third-party sequencing agencies started to play a role in diagnosing H5N6 viruses after the COVID-19 pandemic began. More than one third of cases in 2021 were first detected by hospitals sending samples from patients with pneumonia to third-party agencies.
Evolutionary Relationships of HA and NA Genes of H5N6 Viruses Isolated from Humans
To elucidate the evolutionary pattern of divergence, we analyzed the genetic relationships of the HA and NA genes of 42 human H5N6 viruses with other clade 188.8.131.52 H5N6 genomes available in public databases. According to World Health Organization (WHO) nomenclature of HA genes in clade 184.108.40.206, HPAI H5N6 viruses can be divided into 8 subclades, a–h (Figure 4). All human isolates of H5N6 viruses clustered with vaccine strains recommended by the WHO. The first human virus isolated in 2014, A/Sichuan/26221/2014(H5N6), belonged to subclade 220.127.116.11a; 5 viruses isolated in 2015 were grouped into subclade 18.104.22.168d; and viruses isolated in 2016 belonged to subclade 22.214.171.124g. Subclade 126.96.36.199h viruses were isolated from 2015 through 2021. Of note, most H5N6 human viruses in 2021 (18/20) grouped into genetic subclade 188.8.131.52b. The median tMRCA among the HA genes of clade 184.108.40.206b H5N6 human viruses in 2021 was estimated to be June 16, 2020 (95% highest posterior density March 29, 2020–August 23, 2020).
Phylogenetic analysis of the N6 genes revealed 2 subclades (Figure 5). Almost all viruses fell into the subgroup with an 11-residue deletion at position 59–69 in the NA stalk. Only 1 virus, A/Sichuan/26221/2014(H5N6), had full-length NA protein. The tMRCA among the N6 genes of clade 220.127.116.11b H5N6 human viruses from 2021 was estimated to be January 19, 2015 (95% highest posterior density June 16, 2014–June 17, 2015). These results indicated that clade 18.104.22.168b H5N6 viruses were generated by multiple reassortment events.
Classification and Temporal Distribution of H5N6 Genotypes Isolated from Humans
The phylogenetic analysis exhibited dynamic donation of the 8 gene segments to H5N6 viruses, including those from the Eurasian gene pool, H9N2 AIVs, and so on (Appendix Figure 2). Each gene of the human isolates showed a highly homologous sequence origin from poultry viruses, wild bird viruses, or both. On the basis of different clade combinations of 6 internal genes, we identified a total of 13 H5N6 genotypes, termed 2014A, 2014B, 2015A, 2015B, 2016A, 2016B, 2017A, 2018A, 2020A, 2021A, 2021B, 2021C, and 2021D (Figure 6). Several of these genotypes (e.g., 2014A and 2014B) were only identified in 1 case. In contrast, genotype 2015A was first detected in 2015 but continuously caused human infections during 2017–2021. Genotype 2015B contained a set of internal genes derived from H9N2 viruses in poultry and infected at least 8 persons during 2015–2016. However, no cases infected with this genotype were detected in the years after. Genotype 2020A emerged in 2020 and had 6 internal genes originating from the Eurasian gene pool and clade 22.214.171.124c H5N1 viruses. In 2021, a total of 4 new genotypes (2021A, 2021B, 2021C, and 2021D) acquiring genes from clade 126.96.36.199b H5N8 viruses, H5N1 viruses, and the Eurasian gene pool were identified. Five human isolates from Sichuan and Chongqing provinces belonged to genotype 2021A, and 9 human isolates from Hunan, Guangxi, and Zhejiang provinces were classified into genotype 2021B. Genotype 2021C and 2021D contained 1 virus each. Partial sequences were obtained from another 3 human viruses detected in 2021; their available sequences were closely related to those of the new genotype viruses. These findings might indicate the cross-species advantages of these newly emerged H5N6 genotypes to infect humans.
Key Amino Acid Mutations Occurring in H5N6 Viruses Isolated from Humans
We analyzed molecular substitutions associated with increased virulence and transmissibility in mammals and reduced susceptibility to antiviral drugs. All H5N6 human viruses obtained a 5–amino acid insertion in the HA cleavage site, except for A/Anhui/33162/2016(H5N6), which had 1 more amino acid insertion. Substitution Q226L (H3 numbering) in HA protein, which had previously been reportedly associated with a switch in receptor specificity from avian-type (α2–3Gal) to human-type (α2–6Gal) (34–36), was detected from 1 virus, A/Hunan/09285/2021(H5N6). Substitution S227R in HA protein, which could also alter receptor specificity (37), was detected in 24 human viruses. Another receptor-changing substitution, T192I (38), was newly detected in clade 188.8.131.52b viruses from 2021 cases (Table 4; Appendix Table 3).
Several mammalian-adapted mutations have occurred in PB2 protein of H5N6 viruses (Appendix Table 3). The substitutions E627K and D701N in PB2 protein, which were associated with increased polymerase activities or enhanced virulence in mice, occurred in 12/42 (E627K) and 5/42 (D701N) human viruses. A high proportion of these substitutions were observed in genotype 2015B viruses; 6/8 viruses acquired the E627K substitution and 1/8 viruses acquired the D701N substitution. These phenomena were frequently observed when an AIV transmitted to a mammalian host (39–42). These findings further confirmed that the infection source of H5N6 cases was poultry populations.
All genotype H5N6 viruses had 31S in M2 protein, except for genotype 2015B viruses. The genotype 2015B viruses contained 8 viruses, and exclusively had 31N in M2 protein, indicating their reduced susceptibility to amantadine (43). Mutations associated with drug resistance in the NA protein (E119D) only occurred in 1 human virus A/Anhui/33162/2016(H5N6). Mutation of I38V in PA, which might affect susceptibility to endonuclease inhibitors, was found in A/Jiangsu/32888/2018(H5N6). In addition, residues in the M1 and NS1 proteins of several H5N6 viruses showed changes that might alter virulence in mice (Appendix Table 3).
Among AIVs, viruses of subtypes H5 and H7 have been of particular concern because of their high rates of death. The hospitalization fatality risk for H5N6 infections (59.0%) was slightly lower than that for H5N1 infections (70.0%) but higher than that for H7N9 (35.0%) (44). In addition, the median age of patients with H5N6 virus (51 years) was older than the median age for patients with H5N1 virus (26 years) but younger than the age for patients with H7N9 virus (62 years) (44).
On the basis of epidemiologic investigations, 93.8% of influenza H5N6 case-patients were confirmed to have poultry exposure history. The contamination of live poultry markets and backyard birds, as well as the practice of processing poultry without personal protection, could be ongoing exposure sources for influenza A(H5N6) virus. While H5N6 viruses continue to circulate in poultry, human infections will undoubtedly continue.
WHO has recommended early antiviral therapy, ideally within 48 hours of symptom onset, for suspected or confirmed influenza patients (45). Our study found that early initiation of antiviral treatments could reduce the fatality rate in H5N6 patients to some extent. However, symptoms at the onset of influenza H5N6 infections were clinically similar to those of other respiratory pathogen infections. Thus, increased sensitivity of diagnostic systems is needed to improve case identification and initiate timely antiviral treatment.
During the COVID-19 pandemic, diagnostic capacity for respiratory illnesses among human health systems in China, including hospitals, Centers for Disease Control and Prevention at different levels, and third-party agencies, was increased. In our study, clinical samples from 15 H5N6 cases during the COVID-19 pandemic (2 in 2020 and 13 in 2021) were first identified by third party agencies. This additional diagnostic capacity contributed to the detection of H5N6 cases reported in 2021.
Genetic analysis in our study detected at least 13 types of reassortant H5N6 viruses in infected humans in China (Figure 6). Origins of internal genes were dramatically diversified, indicating the advanced genetic compatibility of H5N6 viruses with other AIVs. In 2021, a total of 4 new H5N6 genotype viruses emerged and accounted for almost all of the H5N6 human viruses based on the available full genome. Of note, different from previously isolated viruses, the HA gene of almost all H5N6 human viruses in 2021 belonged to genetic clade 184.108.40.206b and was closely related to that of the first H5N8 isolate, which caused infection in a patient in Russia (46). Additional mammal-adapted mutations, including Q226L and T192I in the HA protein, which could increase the viral affinity for human cells, were also detected (Appendix Table 3), indicating the viral adaptation process from birds to humans.
In summary, although we observed a rise in the number of influenza A(H5N6) infections in 2021, the disease course and CFR were comparable to previously detected H5N6 cases. Antiviral drugs remain effective if used early. However, new genotype viruses and mammal-adapted substitutions emerged. Moreover, reports from OFFLU (https://www.offlu.org) and WHO have documented that clade 220.127.116.11h and 18.104.22.168b H5N6 viruses and clade 22.214.171.124b H5N8 viruses have been detected in poultry and wild birds in China. Considering the continuous viral circulation in birds and incidence of human infection, more H5N6 variants and genotypes with further advantages in humans might emerge. Increased attention to such emerging viruses is vital for public health and pandemic preparedness.
Dr. Zhu works in the Chinese National Influenza Center, National Institute for Viral Disease Control and Prevention, China Center for Disease Control and Prevention. Her research interests include evolutionary analysis and pathogenicity mechanism of influenza viruses.
We thank the Chinese National Influenza Surveillance Network and the authors and laboratories for sharing the influenza A(H5N6) sequences in GISAID’s EpiFlu Database. We thank the following provincial Centers for Disease Control and Prevention for their data sharing on H5N6 cases: Anhui, Beijing, Chongqing, Fujian, Guangdong, Guangxi, Guizhou, Hubei, Hunan, Jiangsu, Jiangxi, Sichuan, Yunnan, and Zhejiang.
This study was supported by the National Nature Science Foundation of China (81961128002, 81971941, and 31970643). The funding bodies had no role in study design, data collection and analysis, preparation of the manuscript, or the decision to publish.
- Pan M, Gao R, Lv Q, Huang S, Zhou Z, Yang L, et al. Human infection with a novel, highly pathogenic avian influenza A (H5N6) virus: Virological and clinical findings. J Infect. 2016;72:52–9.
- Chen H, Yuan H, Gao R, Zhang J, Wang D, Xiong Y, et al. Clinical and epidemiological characteristics of a fatal case of avian influenza A H10N8 virus infection: a descriptive study. Lancet. 2014;383:714–21.
- Gao R, Cao B, Hu Y, Feng Z, Wang D, Hu W, et al. Human infection with a novel avian-origin influenza A (H7N9) virus. N Engl J Med. 2013;368:1888–97.
- Wang D, Zhu W, Yang L, Shu Y. The epidemiology, virology, and pathogenicity of human infections with avian influenza viruses. Cold Spring Harb Perspect Med. 2021;11:
- Yang L, Zhao X, Li X, Bo H, Li D, Liu J, et al. Case report for human infection with a highly pathogenic avian influenza A(H5N6) virus in Beijing, China 2019. Biosafety and Health. 2020;2:49–52.
- World Health Organization. Cumulative number of confirmed human cases for avian influenza A(H5N1) reported to WHO, 2003–2021, 1 October 2021 [cited 2021 Nov 19]. https://www.who.int/publications/m/item/cumulative-number-of-confirmed-human-cases-for-avian-influenza-a(h5n1)-reported-to-who-2003-2021-1-october-2021
- Li KS, Guan Y, Wang J, Smith GJ, Xu KM, Duan L, et al. Genesis of a highly pathogenic and potentially pandemic H5N1 influenza virus in eastern Asia. Nature. 2004;430:209–13.
- Jeong J, Kang HM, Lee EK, Song BM, Kwon YK, Kim HR, 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.
- Lee DH, Torchetti MK, Winker K, Ip HS, Song CS, Swayne DE. Intercontinental spread of Asian-origin H5N8 to North America through Beringia by migratory birds. J Virol. 2015;89:6521–4.
- Butler J, Stewart CR, Layton DS, Phommachanh P, Harper J, Payne J, et al. Novel reassortant H5N6 influenza A virus from the Lao People’s Democratic Republic is highly pathogenic in chickens. PLoS One. 2016;11:
- Yang L, Zhu W, Li X, Bo H, Zhang Y, Zou S, et al. Genesis and dissemination of highly pathogenic H5N6 avian influenza viruses. J Virol. 2017;91:e02199–16.
- Chu DH, Okamatsu M, Matsuno K, Hiono T, Ogasawara K, Nguyen LT, et al. Genetic and antigenic characterization of H5, H6 and H9 avian influenza viruses circulating in live bird markets with intervention in the center part of Vietnam. Vet Microbiol. 2016;192:194–203.
- Bi Y, Chen Q, Wang Q, Chen J, Jin T, Wong G, et al. Genesis, evolution and prevalence of H5N6 avian influenza viruses in China. Cell Host Microbe. 2016;20:810–21.
- Xu W, Li X, Bai T, Zhao X, Zhao X, Zhang Y, et al. A fatal case of infection with a further reassortant, highly pathogenic avian influenza (HPAI) H5N6 virus in Yunnan, China. Infect Genet Evol. 2016;40:63–6.
- Yang L, Zhu W, Li X, Bo H, Zhang Y, Zou S, et al. Genesis and dissemination of highly pathogenic H5N6 avian influenza viruses. J Virol. 2017;91:e02199–16.
- Cui Y, Li Y, Li M, Zhao L, Wang D, Tian J, et al. Evolution and extensive reassortment of H5 influenza viruses isolated from wild birds in China over the past decade. Emerg Microbes Infect. 2020;9:1793–803.
- Zhang Y, Chen M, Huang Y, Zhu W, Yang L, Gao L, et al. Human infections with novel reassortant H5N6 avian influenza viruses in China. Emerg Microbes Infect. 2017;6:
- Sun W, Li J, Hu J, Jiang D, Xing C, Zhan T, et al. Genetic analysis and biological characteristics of different internal gene origin H5N6 reassortment avian influenza virus in China in 2016. Vet Microbiol. 2018;219:200–11.
- Chen Z, Wang Z, Zhao X, Guan Y, Xue Q, Li J, et al. Pathogenicity of different H5N6 highly pathogenic avian influenza virus strains and host immune responses in chickens. Vet Microbiol. 2020;246:
- Sun H, Pu J, Hu J, Liu L, Xu G, Gao GF, et al. Characterization of clade 126.96.36.199 highly pathogenic H5 avian influenza viruses in ducks and chickens. Vet Microbiol. 2016;182:116–22.
- Li Y, Li M, Li Y, Tian J, Bai X, Yang C, et al. Outbreaks of highly pathogenic avian influenza (H5N6) virus subclade 188.8.131.52h in Swans, Xinjiang, Western China, 2020. Emerg Infect Dis. 2020;26:2956–60.
- Sun H, Pu J, Wei Y, Sun Y, Hu J, Liu L, et al. Highly pathogenic avian influenza H5N6 Viruses exhibit enhanced affinity for human type sialic acid receptor and in-contact transmission in model ferrets. J Virol. 2016;90:6235–43.
- Bi Y, Liu H, Xiong C, Di Liu , Shi W, Li M, et al. Novel avian influenza A (H5N6) viruses isolated in migratory waterfowl before the first human case reported in China, 2014. Sci Rep. 2016;6:29888.
- Koutsakos M, Wheatley AK, Laurie K, Kent SJ, Rockman S. Influenza lineage extinction during the COVID-19 pandemic? Nat Rev Microbiol. 2021;19:741–2.
- World Health Organization. Assessment of risk associated with highly pathogenic avian influenza A(H5N6) virus [cited 2021 Nov 19]. https://cdn.who.int/media/docs/default-source/influenza/avian-and-other-zoonotic-influenza/a(h5n6)-risk-assessment.pdf
- Yu L, Wang Z, Chen Y, Ding W, Jia H, Chan JF, et al. Clinical, virological, and histopathological manifestations of fatal human infections by avian influenza A(H7N9) virus. Clin Infect Dis. 2013;57:1449–57.
- Zhou B, Donnelly ME, Scholes DT, St George K, Hatta M, Kawaoka Y, et al. Single-reaction genomic amplification accelerates sequencing and vaccine production for classical and Swine origin human influenza a viruses. J Virol. 2009;83:10309–13.
- McGinnis J, Laplante J, Shudt M, George KS. Next generation sequencing for whole genome analysis and surveillance of influenza A viruses. [Erratum in: J Clin Virol. 2017;93:65]. J Clin Virol. 2016;79:44–50.
- Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30:1312–3.
- Drummond AJ, Rambaut A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007;7:214.
- Drummond AJ, Rambaut A, Shapiro B, Pybus OG. Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol. 2005;22:1185–92.
- Fiore AE, Shay DK, Broder K, Iskander JK, Uyeki TM, Mootrey G, et al.; Centers for Disease Control and Prevention (CDC); Advisory Committee on Immunization Practices (ACIP). Prevention and control of influenza: recommendations of the Advisory Committee on Immunization Practices (ACIP), 2008. MMWR Recomm Rep. 2008;57(RR-7):1–60.
- Influenza HE. In: Epidemiology and prevention of vaccine-preventable diseases, 14th edition. 2021 [cited 2022 Apr 21]. https://www.cdc.gov/vaccines/pubs/pinkbook/flu.html#impact
- Shi Y, Zhang W, Wang F, Qi J, Wu Y, Song H, et al. Structures and receptor binding of hemagglutinins from human-infecting H7N9 influenza viruses. Science. 2013;342:243–7.
- Xu R, de Vries RP, Zhu X, Nycholat CM, McBride R, Yu W, et al. Preferential recognition of avian-like receptors in human influenza A H7N9 viruses. Science. 2013;342:1230–5.
- Ramos I, Krammer F, Hai R, Aguilera D, Bernal-Rubio D, Steel J, et al. H7N9 influenza viruses interact preferentially with α2,3-linked sialic acids and bind weakly to α2,6-linked sialic acids. J Gen Virol. 2013;94:2417–23.
- Guo H, de Vries E, McBride R, Dekkers J, Peng W, Bouwman KM, et al. Highly pathogenic influenza A(H5Nx) viruses with altered H5 receptor-binding specificity. Emerg Infect Dis. 2017;23:220–31.
- Yang ZY, Wei CJ, Kong WP, Wu L, Xu L, Smith DF, et al. Immunization by avian H5 influenza hemagglutinin mutants with altered receptor binding specificity. Science. 2007;317:825–8.
- Steel J, Lowen AC, Mubareka S, Palese P. Transmission of influenza virus in a mammalian host is increased by PB2 amino acids 627K or 627E/701N. PLoS Pathog. 2009;5:
- Li Z, Chen H, Jiao P, Deng G, Tian G, Li Y, et al. Molecular basis of replication of duck H5N1 influenza viruses in a mammalian mouse model. J Virol. 2005;79:12058–64.
- Subbarao EK, London W, Murphy BR. A single amino acid in the PB2 gene of influenza A virus is a determinant of host range. J Virol. 1993;67:1761–4.
- Gabriel G, Dauber B, Wolff T, Planz O, Klenk HD, Stech J. The viral polymerase mediates adaptation of an avian influenza virus to a mammalian host. Proc Natl Acad Sci U S A. 2005;102:18590–5.
- Lan Y, Zhang Y, Dong L, Wang D, Huang W, Xin L, et al. A comprehensive surveillance of adamantane resistance among human influenza A virus isolated from mainland China between 1956 and 2009. Antivir Ther. 2010;15:853–9.
- Cowling BJ, Jin L, Lau EH, Liao Q, Wu P, Jiang H, et al. Comparative epidemiology of human infections with avian influenza A H7N9 and H5N1 viruses in China: a population-based study of laboratory-confirmed cases. Lancet. 2013;382:129–37.
- Fiore AE, Fry A, Shay D, Gubareva L, Bresee JS, Uyeki TM; Centers for Disease Control and Prevention (CDC). Antiviral agents for the treatment and chemoprophylaxis of influenza --- recommendations of the Advisory Committee on Immunization Practices (ACIP). MMWR Recomm Rep. 2011;60:1–24.
- Pyankova OG, Susloparov IM, Moiseeva AA, Kolosova NP, Onkhonova GS, Danilenko AV, et al. Isolation of clade 184.108.40.206b A(H5N8), a highly pathogenic avian influenza virus, from a worker during an outbreak on a poultry farm, Russia, December 2020. Euro Surveill. 2021;26:26.
Suggested citation for this article: Zhu W, Li X, Dong J, Bo H, Liu J, Yang J, et al. Epidemiologic, clinical, and genetic characteristics of human infections with influenza A(H5N6) viruses, China. Emerg Infect Dis. 2022 Jul [date cited]. https://doi.org/10.3201/eid2807.212482
Original Publication Date: April 27, 2022
1These authors contributed equally to this article.