Skip directly to site content Skip directly to page options Skip directly to A-Z link Skip directly to A-Z link Skip directly to A-Z link
Volume 28, Number 7—July 2022
Research

Epidemiologic, Clinical, and Genetic Characteristics of Human Infections with Influenza A(H5N6) Viruses, China

Author affiliations: National Institute for Viral Disease Control and Prevention, Chinese Center for Disease Control and Prevention, Beijing, China (W. Zhu, X. Li, J. Dong, H. Bo, J. Liu, Jiaying Yang, Y. Zhang, H. Wei, W. Huang, X. Zhao, T. Chen, Jing Yang, Z. Li, X. Zeng, C. Li. J. Tang, L. Xin, R. Gao, L. Liu, M. Tan, Y. Shu, L. Yang, D. Wang); Chinese Center for Disease Control and Prevention, Beijing (C. Li); World Health Organization Collaborating Centre for Reference and Research on Influenza, Beijing (W. Zhu, X. Li, J. Dong, H. Bo, J. Liu, Jiaying Yang, Y. Zhang, H. Wei, W. Huang, X. Zhao, T. Chen, Jing Yang, Z. Li, X. Zeng, J. Tang, L. Xin, R. Gao, L. Liu, M. Tan, Y. Shu, L. Yang, D. Wang); Key Laboratory for Medical Virology, National Health and Family Planning Commission, Beijing (W. Zhu, X. Li, J. Dong, H. Bo, J. Liu, Jiaying Yang, Y. Zhang, H. Wei, W. Huang, X. Zhao, T. Chen, Jing Yang, Z. Li, X. Zeng, J. Tang, L. Xin, R. Gao, L. Liu, M. Tan, Y. Shu, L. Yang, D. Wang); Sun Yat-sen University, Guangzhou, China (Jiaying Yang, Y. Shu)

Cite This Article

Abstract

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 2.3.4.4b; 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 2.3.4.4 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 (14). 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 2.3.4.4 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 2.3.4.4 H5 viruses spread westward and eastward from Asia to other continents, accompanied by multiple reassortment events between clade 2.3.4.4 H5 viruses and other AIVs circulating in wild and domestic birds (8,9). Clade 2.3.4.4 H5N6 AIV has been mainly endemic among birds in China and southeast Asia since 2013 (1012), 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 (1416). 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,1618).

Genetic and biologic characteristics indicated that H5N6 viruses were highly pathogenic in chickens (1921). Studies have shown that some H5N6 viruses of clade 2.3.4.4 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 2.3.4.4 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.

Methods

Study Cases

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.

Ethics

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.

Virus Isolation

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.

Statistical Analysis

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).

Results

Epidemiologic Findings

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.

Figure 1

Temporal distribution of 65 human infections with influenza A(H5N6) virus, by month, China, April 21, 2014–December 31, 2021.

Figure 1. Temporal distribution of 65 human infections with influenza A(H5N6) virus, by month, China, April 21, 2014–December 31, 2021.

Figure 2

Spatial distribution of 65 human infections with influenza A(H5N6) virus, China, by year, April 21, 2014–December 31, 2021.

Figure 2. Spatial distribution of 65 human infections with influenza A(H5N6) virus, China, by year, April 21, 2014–December 31, 2021.

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.

Clinical Findings

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).

Figure 3

Disease courses of human infections with influenza A (H5N6) virus, China, April 21, 2014–December 31, 2021. Intervals are given as median (ICU) days. ICU, intensive care unit; IQR, interquartile range.

Figure 3. Disease courses of human infections with influenza A (H5N6) virus, China, April 21, 2014–December 31, 2021. Intervals are given as median (ICU) days. ICU, intensive care unit; IQR, interquartile range.

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

Figure 4

Maximum clade credibility trees of hemagglutinin gene of influenza A(H5N6) viruses, China. Red indicates human-origin H5N6 viruses; blue triangles indicate H5Ny vaccine strains recommended by the World Health Organization; blue dot indicates the most recent common ancestry of clade 2.3.4.4b A(H5N6) human viruses in 2021. Posterior probabilities >0.9 are labeled on the branches.

Figure 4. Maximum clade credibility trees of hemagglutinin gene of influenza A(H5N6) viruses, China. Red indicates human-origin H5N6 viruses; blue triangles indicate H5Ny vaccine strains recommended by the World Health Organization; blue...

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 2.3.4.4 H5N6 genomes available in public databases. According to World Health Organization (WHO) nomenclature of HA genes in clade 2.3.4.4, 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 2.3.4.4a; 5 viruses isolated in 2015 were grouped into subclade 2.3.4.4d; and viruses isolated in 2016 belonged to subclade 2.3.4.4g. Subclade 2.3.4.4h viruses were isolated from 2015 through 2021. Of note, most H5N6 human viruses in 2021 (18/20) grouped into genetic subclade 2.3.4.4b. The median tMRCA among the HA genes of clade 2.3.4.4b H5N6 human viruses in 2021 was estimated to be June 16, 2020 (95% highest posterior density March 29, 2020–August 23, 2020).

Figure 5

Maximum clade credibility tree of neuraminidase gene of influenza A(H5N6) viruses, China. Red indicates human-origin H5N6 viruses; blue triangles indicate vaccine strains recommended by the World Health Organization; blue dot indicates the most recent common ancestry of clade 2.3.4.4b H5N6 human viruses in 2021. Posterior probabilities >0.9 are labeled on the branches.

Figure 5. Maximum clade credibility tree of neuraminidase gene of influenza A(H5N6) viruses, China. Red indicates human-origin H5N6 viruses; blue triangles indicate vaccine strains recommended by the World Health Organization; blue dot...

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 2.3.4.4b 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 2.3.4.4b H5N6 viruses were generated by multiple reassortment events.

Classification and Temporal Distribution of H5N6 Genotypes Isolated from Humans

Figure 6

Diversity and prevalence of influenza A(H5N6) viruses isolated from humans, China. Circles represent the corresponding virus genotypes and their times of isolation. Gene segments are ordered as polymerase basic 2, polymerase basic 1, polymerase acidic, hemagglutinin, nucleoprotein, neuraminidase, matrix, and nonstructural from top to bottom within circles. A total of 13 genotypes are listed; the genotype name is shown under each circle. Names of each human H5N6 virus are listed besides the genotype to which they belong. To illustrate the history of reassortant events, segments in descendant viruses are colored according to their corresponding source viruses on the bottom line.

Figure 6. Diversity and prevalence of influenza A(H5N6) viruses isolated from humans, China. Circles represent the corresponding virus genotypes and their times of isolation. Gene segments are ordered as polymerase basic 2,...

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 2.3.2.1c H5N1 viruses. In 2021, a total of 4 new genotypes (2021A, 2021B, 2021C, and 2021D) acquiring genes from clade 2.3.4.4b 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) (3436), was detected from 2 viruses, A/ Sichuan-Luzhou/LZ20211439-Q1/2021(H5N6) and A/Hunan/09285/2021(H5N6). Substitution S227R in HA protein, which could also alter receptor specificity (37), was detected in 32 human viruses. Another receptor-changing substitution, T192I (38), was newly detected in clade 2.3.4.4b 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 (3942). 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).

Discussion

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 2.3.4.4b 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 2.3.4.4h and 2.3.4.4b H5N6 viruses and clade 2.3.4.4b 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.

Top

Acknowledgments

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.

Top

References

  1. 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:529. DOIPubMedGoogle Scholar
  2. 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:71421. DOIPubMedGoogle Scholar
  3. 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:188897. DOIPubMedGoogle Scholar
  4. 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:a038620. DOIPubMedGoogle Scholar
  5. 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:4952. DOIGoogle Scholar
  6. 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
  7. 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:20913. DOIPubMedGoogle Scholar
  8. 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:24957. DOIPubMedGoogle Scholar
  9. 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:65214. DOIPubMedGoogle Scholar
  10. 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:e0162375. DOIPubMedGoogle Scholar
  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:e0219916. DOIPubMedGoogle Scholar
  12. 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:194203. DOIPubMedGoogle Scholar
  13. 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:81021. DOIPubMedGoogle Scholar
  14. 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:636. DOIPubMedGoogle Scholar
  15. 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:e0219916. DOIPubMedGoogle Scholar
  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:1793803. DOIPubMedGoogle Scholar
  17. 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:e50. DOIPubMedGoogle Scholar
  18. 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:20011. DOIPubMedGoogle Scholar
  19. 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:108745. DOIPubMedGoogle Scholar
  20. Sun  H, Pu  J, Hu  J, Liu  L, Xu  G, Gao  GF, et al. Characterization of clade 2.3.4.4 highly pathogenic H5 avian influenza viruses in ducks and chickens. Vet Microbiol. 2016;182:11622. DOIPubMedGoogle Scholar
  21. Li  Y, Li  M, Li  Y, Tian  J, Bai  X, Yang  C, et al. Outbreaks of highly pathogenic avian influenza (H5N6) virus subclade 2.3.4.4h in Swans, Xinjiang, Western China, 2020. Emerg Infect Dis. 2020;26:295660. DOIPubMedGoogle Scholar
  22. 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:623543. DOIPubMedGoogle Scholar
  23. 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. DOIPubMedGoogle Scholar
  24. Koutsakos  M, Wheatley  AK, Laurie  K, Kent  SJ, Rockman  S. Influenza lineage extinction during the COVID-19 pandemic? Nat Rev Microbiol. 2021;19:7412. DOIPubMedGoogle Scholar
  25. 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
  26. 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:144957. DOIPubMedGoogle Scholar
  27. 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:1030913. DOIPubMedGoogle Scholar
  28. 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:4450. DOIPubMedGoogle Scholar
  29. Stamatakis  A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30:13123. DOIPubMedGoogle Scholar
  30. Drummond  AJ, Rambaut  A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007;7:214. DOIPubMedGoogle Scholar
  31. Drummond  AJ, Rambaut  A, Shapiro  B, Pybus  OG. Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol. 2005;22:118592. DOIPubMedGoogle Scholar
  32. 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):160.PubMedGoogle Scholar
  33. 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
  34. 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:2437. DOIPubMedGoogle Scholar
  35. 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:12305. DOIPubMedGoogle Scholar
  36. 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:241723. DOIPubMedGoogle Scholar
  37. 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:22031. DOIPubMedGoogle Scholar
  38. 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:8258. DOIPubMedGoogle Scholar
  39. 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:e1000252. DOIPubMedGoogle Scholar
  40. 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:1205864. DOIPubMedGoogle Scholar
  41. 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:17614. DOIPubMedGoogle Scholar
  42. 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:185905. DOIPubMedGoogle Scholar
  43. 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:8539. DOIPubMedGoogle Scholar
  44. 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:12937. DOIPubMedGoogle Scholar
  45. 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:124.PubMedGoogle Scholar
  46. Pyankova  OG, Susloparov  IM, Moiseeva  AA, Kolosova  NP, Onkhonova  GS, Danilenko  AV, et al. Isolation of clade 2.3.4.4b 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. DOIPubMedGoogle Scholar

Top

Figures
Tables

Top

Cite This Article

DOI: 10.3201/eid2807.212482

Original Publication Date: April 27, 2022

1These authors contributed equally to this article.

Table of Contents – Volume 28, Number 7—July 2022

EID Search Options
presentation_01 Advanced Article Search – Search articles by author and/or keyword.
presentation_01 Articles by Country Search – Search articles by the topic country.
presentation_01 Article Type Search – Search articles by article type and issue.

Top

Comments

Please use the form below to submit correspondence to the authors or contact them at the following address:

Lei Yang or Dayan Wang, National Institute for Viral Disease Control and Prevention, China CDC. 155 Changbai Rd, Changping District, Beijing, 102206, China

Send To

10000 character(s) remaining.

Top

Page created: April 16, 2022
Page updated: June 18, 2022
Page reviewed: June 18, 2022
The conclusions, findings, and opinions expressed by authors contributing to this journal do not necessarily reflect the official position of the U.S. Department of Health and Human Services, the Public Health Service, the Centers for Disease Control and Prevention, or the authors' affiliated institutions. Use of trade names is for identification only and does not imply endorsement by any of the groups named above.
file_external