Skip to main content

New H6 influenza virus reassortment strains isolated from Anser fabalis in Anhui Province, China



H6 subtype avian influenza viruses are globally distributed and, in recent years, have been isolated with increasing frequency from both domestic and wild bird species as well as infected humans. Many reports have examined the viruses in the context of poultry or several wild bird species, but there is less information regarding their presence in migratory birds.


Hemagglutination and hemagglutination inhibition tests were used to measure HA activity for different HA subtypes. Whole viral genomes were sequenced and analysed using DNAstar and MEGA 6 to understand their genetic evolution. Pathogenicity was evaluated using a mouse infection model.


We isolated 13 strains of H6 virus from faecal samples of migratory waterfowl in Anhui Province of China in 2014. Phylogenetic analysis showed gene reassortment between Eurasian and North American lineages. Five of the identified H6 strains had the ability to infect mice without adaptation.


Our findings suggest that regular surveillance of wild birds, especially migratory birds, is important for providing early warning and control of avian influenza outbreaks.


Avian influenza virus (AIV) is an important zoonotic pathogen [1] and can be classified into 16 hemagglutinin (HA) subtypes and 9 neuraminidase (NA) subtypes based on the antigenicity of these two surface glycoproteins. At least 136 species have been discovered in 26 different families of wild birds and harbour 144 subtypes of influenza A virus [14]. Most subtypes are found not only in common aquatic birds but also in lower mammals, such as bats, as has been demonstrated in recent reports [1, 57]. In Qinghai Lake, China in 2005, massive numbers of migrating birds died by H5N1, and this virus subsequently spread worldwide [8, 9]. In recent years, different influenza virus subtypes (such as H5, H6, H7, and H10) have emerged in succession in poultry, crossed the species barrier and caused frequent outbreaks of human infection [1, 1015]. The continuous outbreaks of avian influenza in recent decades have alarmed and increased focus on the role of wild birds, as major reservoirs, that maintain the persistence and variation of AIV, facilitating viral spread and disease occurrence.

H6 AIV was first isolated from turkeys in Massachusetts, USA in 1965 [10, 16]. Currently, H6 AIVs has a worldwide distribution, and strains of the virus have been detected in various animal species [17]. On June 21, 2013, the first case of a human infection by an avian-origin H6N1 influenza A virus (A/Taiwan/2/2013, Taiwan2) was reported by the Taiwan Centres for Disease Control [11, 18]. Human infection with H6 AIVs in mainland China has not been reported, but a serum antibody positive for the H6 virus was found in poultry workers [19]. More than 30% of the H6 AIVs circulating in poultry in China have enhanced affinity to human-like receptors (ɑ-2, 6 NeuAcGal). Some H6 subtypes can also infect mice without prior adaptation, and some can be efficiently transmitted among guinea pigs [17]. These events indicate that H6 AIV can cross species barriers and directly infect humans. As the natural host of AIVs, wild birds play an important role in virus reassortment as well as intra- and interspecies transmission [4, 5]. However, there is little research regarding the biological properties of H6 AIV strains from wild birds and their potential threat to mammals compared to the many studies of H6 AIV in poultry.

In the present study, we isolated 13 H6 AIV strains after surveying wild birds in the wetlands of the National Nature Reserves in Anhui Province, China. For all of these strains, each gene segment was sequenced to analyse phylogenetic and speculate its origin. We also preliminarily evaluated the potential threat of each strain to mice. Our study highlights the importance of performing regular surveillance of AIV strains in Anhui Province, China as an integral part of worldwide efforts to better understand AIV ecology and prevent the emergence of novel, potentially pandemic strains. Surveillance of avian influenza virus (AIVs) in wild birds for early warning, prevention, and control of viral outbreaks should be enhanced to reduce the risk of pathogen emergence from wildlife host reservoirs.


Virus isolation and identification

In spring and autumn of 2014, 2970 faecal samples were collected from wild birds in Anhui Province, China. Fresh, single faecal samples were collected with cotton swabs and placed into 5 ml EP tubes with 2 ml virus protection solution (phosphate-buffered solution (pH 7.2) supplemented with penicillin, streptomycin, and 10% glycerinum). The handheld end of each cotton swab was broken off, and the remainder was placed into the sample tube and covered tightly. The handheld ends of the swabs were inserted into the ground to designate collection locations to avoid repetition in sampling.

Bird species were identified based on morphological characteristics observed by binoculars before sampling and were further determined based on faecal shape and colour in addition to DNA barcoding [20]. The samples were vortexed, oscillation and finally centrifuged, and the collected supernatant was inoculated into 9-day-old specific-pathogen-free (SPF) chicken embryos. Allantoic fluid was harvested after 72 h of culture. The HA activity of the allantoic fluid was then evaluated using an HA test. Following this, influenza virus and HA subtypes were identified using a hemagglutination inhibition (HI) test [21] with H6 mono-antiserum and further verified using subtype-specific real-time PCR (RT-PCR) [22, 23]. Virus was standardized to a concentration of 4 HAU/25 μl prior to using the HI assay to identify HA subtype. Working in a PBC, 25 μL of standardized test virus for each unknown isolate was dispensed into a series of three wells (in triplicate) in a U-bottom microtitre plate. Additionally, 25 μL of a corresponding HA-subtype positive-control antigen was dispensed at a 4 HAU/25 μL concentration into positive-control wells. Then, 25 μl of the appropriate standardized antiserum was added to the first well of the HA-subtype series. The antiserum in the antigen wells was then serially diluted beginning with the first well (25 μL carry-back with the excess 25 μL from the final row being discarded). In this way, the serum was diluted into standardized antigen. Each subtype series was diluted as soon as possible after the addition of the antiserum for that series. After this step, 25 μL of liquid remained in each well. The plate was then covered and incubated at room temperature for 25 min. Following this, 25 μL of 1% chicken erythrocyte suspension was added to each well, and the plate was submitted to gentle shaking/agitation. The erythrocyte solution was mixed periodically during this step to ensure an even suspension of erythrocytes during the dispensing process. The plate was then covered with microtitre plate-sealing tape (the PBS was removed from the plates after sealing with the tape) and incubated at room temperature until a distinct button formed in the positive-control wells, which typically required 20–30 min. The assay plates were initially observed after approximately 20 min of incubation and checked frequently thereafter for evidence of hemagglutination. Because some isolates can elute (i.e., detach from erythrocytes) in as little as 30 min, the time window for evaluation of assay results is short in some cases [21]. NA subtypes were directly analysed using subtype-specific RT-PCR and sequencing analysis. Viral RNA was extracted from the allantoic fluid samples that were positive in the HA test using TRIzol reagent. Next, cDNA was synthesized by reverse transcription with the Uni12 primer. Viral genomes were PCR-amplified using primers complementary to the conserved promoter and non-coding region of each gene segment (Table 1). The PCR reaction contained 1 μl cDNA, 1 μl forward primer and reverse primer, 5 μl 10× Taq buffer (TAKARA, Japan), 4 μl 2.5 mM dNTPs (TAKARA, Japan), 1 μl Ex Taq (TAKARA, Japan) and 37 μl RNase-free water for a final volume of 50 μl. A single PCR program was used for all primers with the following conditions: initial denaturation at 95 °C for 10 min; 30 cycles of 95 °C for 30 s, 56 °C for 30 s, and 72 °C for 1.5 min; and extension at 72 °C for 10 min, after which the reaction was stored at 4 °C. The PCR products were purified using a PCR purification kit (TianGen, China) and sequenced on an Applied Biosystems DNA analyser using specific primers.

Table 1 Primers in RT-PCR

Genetic and phylogenetic analyses

Nucleotide sequences were edited using the SeqMan module of the DNASTAR package, and multiple sequence alignments were compiled using Clustal W [2224]. Phylogenetic analyses were performed with the neighbour-joining method with maximum likelihood trees using MEGA 6.0 software [25, 26]. Bootstrap values of 1,000 were used. The phylogenetic analyses were based on the following coding sequences (nucleotides): polymerase PB2 (PB2), 1 to 2280; polymerase PB1 (PB1), 1 to 2274; polymerase PBA (PA), 1 to 2151; hemagglutinin (HA), 1 to 1701; nucleoprotein (NP), 1 to 1497; neuraminidase (NA), 1 to 1410; matrix protein 1 (M), 1 to 982; and non-structural protein 1 (NS), 1 to 844.

Mouse studies

The 50% egg infection dose (EID50) was calculated according to the protocol of Reed and Muench. Six-week-old female BALB/c mice (Experimental Animal Centre of Vital River, Beijing, China) were randomly divided into four infection groups (8 mice/group) and one control group (5 mice/group). Animals in the infection groups were mildly anesthetized with dry ice and infected by intranasal inoculation with a dose of 106 EID50 (50 μl/mouse). Three of the eight mice were randomly euthanized on day 3 post inoculation (p.i.) for titration of virus from the lungs, nasal turbinate, kidneys, spleen, and brain. The remaining five mice were monitored daily for 2 weeks for changes in body weight and mortality.


Influenza virus isolation and subtype identification

In total, 13 H6 AIV strains were isolated from 2970 faecal samples collected from wild birds in Anhui Province, China, in 2014. The strains consisted of two NA subtypes, H6N1 (3 strains) and H6N2 (10 strains), and the host for all of these viruses was Anser fabalis (bean goose) (Table 2).

Table 2 Information of the 13 avian influenza viruses

Phylogenetic analysis of whole genes

The HA genes of all 13 H6 isolates, including 3 H6N1 viruses and 10 H6N2 viruses, belonged to the Eurasian lineage and were further divided into 2 groups according to nucleic acid homology (>95%). The homology of the 13 viral strains ranged between 93.1 and 100%. Group 1 consisted of only one virus, AH/L221, which exhibited the highest sequence identity with A/MDK/Vietnam/LBM455/2013(H6N2) in a GenBank search as well as high sequence identity with A/duck/eastern China/34/2005(H6N1) isolated from poultry. The 12 strains in group 2 exhibited high sequence identity with A/QA/Korea/CN20/2009(H6N1). Among poultry influenza viruses, these 12 viruses exhibited high sequence identity with A/duck/eastern China/11/2008(H6N1) (Fig. 1a). The three N1 NA genes formed only one group, and their homology was greater than 98% (Fig. 1b). The ten N2 NA genes exhibited greater diversity than the HA genes and formed 3 groups. Three N2 genes had a North American lineage, and seven other NA genes had a Eurasian lineage. The homology of the nucleotide sequences ranged between 86.5 and 99.8% (Fig. 1c).

Fig. 1
figure 1

Phylogenetic analysis based on maximum likelihood of surface genes H6 (a), N1 (b) and N2 (c) of AIVs isolated in Anhui Province, 2014. Phylogenetic trees were generated with the MAGE 6.0 software package. The phylogenetic trees for HA (1A) were rooted to A/Turkey/Canada/1/1963 (H6N8). The phylogenetic tree for the N1 genes were rooted to A/Brevig_Mission/1/18(H1N1) and N2 genes were rooted to turkey/Massachussetts/3740/65(H6N2). The genomic sequences of the viruses listed in black and blue were downloaded from available databases; the viruses listed in red, yellow, and/or purple was sequenced in this study. The sequences shown in cyan and italics in Fig 1c were obtained from a paper by Wang GJ. Abbreviations: CK, chicken; DK, duck; GS, goose; SW, swine; ML, mallard; GT, green-winged teal; PD, Pacific black duck; WDK, wild duck; MLDK, mallard duck; WB, wild bird; EN, environment

To better understand the evolution of the H6 AIV strains in this study, 6 internal genes from the 13 viruses were phylogenetically analysed as a whole. These genes exhibited more diversity than the surface genes and were further divided into eight genotypes (Table 3). The nucleotide and amino acid homologies for the 13 H6 isolates in this study were 86.8–100% and 96.1–99.9% for the PB2 gene, 95.6–100% and 99.1–100% for the PB1 gene, 90.3–100% and 97.9–100% for the PA gene, 93.5–100% and 96.4–100% for the NP gene, 96.6–100% and 98.5–100% for the M gene, and 71.2–100% and 63.7–100% for the NS gene, respectively. All of the internal genes belong to a Eurasian lineage except for the PA gene of strain AH/L180 and the NS gene of strain AH/L93, which belong to the North American lineage. The genes were further classified into various subgroups within each lineage. The PB2 gene originated from poultry and wild birds in countries outside of China (i.e., Vietnam, Korea, and Republic of Georgia) and was divided into three groups (Fig. 2a). The PB1 gene, which showed few differences in nucleotide sequences among different strains, was classified into one group (Fig. 2b), while the PA gene was classified into three groups (Fig. 2c). The NP genes from the 13 strains were divided into two groups (Fig. 2d): three originated from ducks in Vietnam and the other ten originated from mallard in Republic of Georgia. The M genes of all 13 stains had high identity and were placed into one group (Fig. 2e). In addition, the NS phylogenetic tree showed that the NS genes of the H6 AIVs were clearly divided into two genetic lineages, termed allele A and allele B, and most of the H6 AIVs contained Eurasian allele B, although one virus from strain AH/L93 had allele A. The NS gene was classified into two groups (Fig. 2f).

Table 3 The genotypes of 13 H6 AIVs
Fig. 2
figure 2

Phylogenetic analysis by maximum likelihood of inner genes of H6 subtype AIVs isolated in Anhui Province, 2014: (a) PB2, (b) PB1, (c) PA, (d) NP, (e) M, and (f) NS. The colours of the virus names listed in the PB2, PB1, PA, NP, M, and NS trees match those used in the genotype table. The phylogenetic trees of inner genes were rooted to A/Turkey/Canada/1/1963 (H6N8). The nucleotide sequences of the viruses listed in black were downloaded from GenBank; the viruses listed in red, yellow, or green were sequenced in this study. The sequences shown in cyan and italics were obtained from a paper by Wang GJ. Abbreviations: CK, chicken; DK, duck; GS, goose; SW, swine; WDK, wild duck; WB, wild bird; WGS, wild goose; EN, environment

Preliminary analyses of the possible origins of the strains were conducted. There were three Genotype 1 viruses, which were tri-reassortment of H6N2, H5N1 and H3N6. There were three Genotype 2 viruses, which were tri-reassortment of H6N2, H9N2 and H3N6. There were two Genotype 3 viruses, which were bi-reassortment of H6N2 and H9N2. There was only one virus in Genotype 4 and one in Genotype 5: both were tri-reassortment of H6N2, H9N2 and H5N1. There was only one virus in Genotype 6, which was a bi-reassortment of H6N2 and H9N2. There was only one virus in Genotype 7 and one in Genotype 8: both were tri-reassortments of H6N2, H6N1 and H3N6 (Fig. 4).

Molecular characterization of viral genes

The isolated H6 strains, except for AH/L221, contained seven potential glycosylation sites according to analysis performed at; these included sites at positions 26 to 28, 27 to 29, 39 to 41, 306 to 308, and 311 to 313 on HA1 and at positions 498 to 500 and 557 to 559 on HA2 (H3 numbering, which is used throughout the manuscript). Strain AH/L221 contained eight potential glycosylation sites, including a glycosylation site on NNT at positions 182 to 184. All viruses except for strain AH/L144 contained the PQIETR/G motif at the cleavage site between HA1 and HA2. The cleavage site on AH/L144 contained the PLIETR/G motif, which includes an amino acid substitution from Q to L. None of the strains contained consecutive basic amino acids within the motif, and all conformed to the characteristics of low-pathogenicity AIVs. The amino acid substitutions A138S, E190V, P186L [27], Q226L, and G228S, as well as the absence of glycosylation at positions 158 to 160 in HA, have been reported to increase the affinity of influenza virus for human-type receptors [28]. The above-mentioned changes were not detected in any of the 13 H6 avian viruses sequenced. Necklace deletion in the NA gene confers enhanced virus lethality in mice. In the present study, none of the 13 H6 AIVs studied had a necklace deletion in the NA gene [29]. Some amino acid substitutions, including T271A, E627K, and D701N in PB2, contribute to increased virulence and transmission of influenza viruses in mammals. The H6 viruses sequenced in this study did not contain any of these changes either. No Y436H substitutions in the PB1 protein or T515A substitutions in the PA protein were found, which suggests that the studied viruses have low pathogenicity for mammalian and avian hosts. Additionally, no amino acid substitutions were found in the M2 transmembrane domain, suggesting that the viruses are sensitive to M2 ion channel inhibitors [30].

The N66S substitution in the PB1-F2 protein, which has been associated with the increased virulence of the 1918 pandemic virus and the high pathogenicity of the AI H5N1 virus in mice and ferrets [31, 32], was not found in any of the H6 strains examined in this study. Additionally, none of the viruses harboured the S31N substitution in the M2 protein, indicating that all are sensitive to amantadine inhibitors [33]. The virulence of influenza virus in humans is related to resistance to the antiviral effects of cytokines, such as interferon (IFN). In particular, the D92E mutation in the NS1 protein increases resistance to IFN; this mutation was not observed in any of the strains evaluated in this study [34].

Replication of H6 viruses in mice

To investigate the replication ability and virulence of the isolated H6 strains in mammals [35], five H6 strains were selected according to their NA gene subtype and phylogenetic characteristics of HA and NA and used to infect mice. Following infection, none of the viruses were detected in the spleen, kidneys, or brain of the mice. Four of the five viruses caused the mice to lose weight (Fig. 3a). In virus proliferation tests in chicken embryos, one strain of virus, strain AH/S45 (H6N2), was not detectable in any organs, whereas the four other strains were detected in the lungs (mean titre 2.4 to 6.3 log10 EID50) and nasal turbinate (mean titre 0.6 to 1.3 log10 EID50) (Fig. 3b). These results indicate that the majority of H6 influenza viruses that were isolated from wild birds can replicate in the respiratory systems of mice without preadaptation.

Fig. 3
figure 3

Changes in body weight. Body weight (a) and mortality (b) in BALB/c mice inoculated with strains A and B of H6 viruses. a Body weights and survival rates of mice were observed over 14 days after infection. b Lungs and trachea were collected at 3 days post infection (p.i.), and virus replication levels were measured by EID50 values in specific-pathogen-free (SPF) eggs


Based on body shape, physiological characteristics, and other features, wild birds have unique flight patterns, which has led to their wide geographical spread, strong regional activity and extended seasonal migration. Wild birds migrate from breeding to wintering places each year along unique migration routes. Wild aquatic birds are the natural host of influenza virus and can act as a reservoir for the transfer of genes between species, forming transient “genome constellations” that are continually reshuffled by reassortment [36, 37]. The birds’ migration behaviours expand the distribution of the virus and provide conditions for the cross-border spread of viruses. Many reports have analysed the biological and pathogenic characteristics of influenza virus in poultry in China; however, studies of influenza virus in wild birds are relatively rare.

In recent years, several AIV subtypes have been circulating and evolving in southern China; these include the H5, H9, and H6 subtypes as well as the newly emerged H7N9 and H10N8 viruses. This has led to frequent avian influenza outbreaks in poultry and humans. Anhui Province is located in south-eastern China (29°41′to 34°38′N, 114°54′to 119°37′E). Lying in the middle and lower valleys of the Yangtze and HuaiHe rivers, Anhui is abundant in wetlands, with wetland areas accounting for 21% of the total area of the province. These areas provide important stopover and wintering grounds for birds migrating along the East Asia-Australia migratory flyway. In the fall and winter of each year, a large number of migratory birds gather in the wetlands, resulting in extremely high population densities of the same or different species, which likely favours the spread and viral gene reassortment of the influenza virus. As an important agricultural province, Anhui is also home to a great number of poultry raised under different farming practices. These poultry often come into close contact with wild birds, creating the opportunity for cross-species AIV transmission.

In 2 years of surveillance work in Anhui, we obtained 31 AIV strains from 2970 faecal samples. Among these strains, eight subtype combinations, including H1N1, H1N2, H3N3, H3N8, H6N1, H6N2, H9N2, and H11N9, were identified using HA-HI tests and PCR. The H6 subtype accounted for the largest proportion of those identified, being found in approximately 41.9% (13/31) of the samples. H6N2 accounted for the majority of these H6 AIVs (10/13). Overall, the goal of this work was to obtain baseline information regarding the epidemiological and virological characteristics of H6 virus strains that exist in wild waterfowl. Toward this aim, we assessed 13 strains of H6 virus isolated from faecal samples of bean goose (Anser fabalis) in Anhui Province, China in 2014.

These 13 H6 strains originated from complicated reassortment between viruses found in poultry and wild birds from China and other countries (Fig. 4). Interestingly, all of the HA genes originated from poultry in China. Except for strain AH/L221, which originated from the A/duck/Eastern China/34/2005(H6N1)-like gene pool, the other 12 H6 strains all originated from the A/duck/Hunan/2110/2006(H6N2)-like gene pool. The ten N2 strains all exhibited high sequence identity with wild birds from different countries and were reassortment between H9N2 and H6N2 in bean goose, duck, and mallard. Three N1 genes were found in H6N1 isolates from poultry in their native countries. The internal genes also showed complicated reassortment from varied sources, including poultry and different wild bird species from China and other countries. This finding further shows that wild birds, as the natural host of AIV, play an important role in gene reassortment or rearrangement of viruses, which can subsequently be transmitted to domestic poultry. A former study showed that the gene flow and reassortment of the Eurasian lineage in the North American gene pool dramatically changed the evolutionary dynamics of influenza virus in natural reservoir hosts [38]. The HA, PA and NS genes in the strains examined in this study had high identity with A/bean goose/Korea/220/2011(H9N2), which belongs to the North American lineage [38]. China and Korea are both along the same bird migration pathway, which increases the gene flux of influenza virus between the two countries and may accelerate evolution driven by North American-Eurasian reassortment. During seasonal migration for breeding and wintering, migratory birds can distribute genes to local poultry or wild birds and can acquire additional genes [39]. This invasion of viruses from wild bird reservoirs can increase the risk of emergence of highly pathogenic avian influenza (HPAI) viruses and threaten food security in addition to increasing the chance of virus infection in humans.

Fig. 4
figure 4

A simplified schematic showing the putative genomic composition of the novel H6N1 and H6N2 reassortant avian influenza viruses identified in this study. The eight gene segments (from top to bottom) in each virus are PB2, PB1, PA, HA, NP, NA, M, and NS. Each colour represents a separate source background. Eight genotypes are listed here. The six different colours represent separate virus backgrounds. The illustration is based on nucleotide-distance comparisons and phylogenetic analyses

Gillim-Ross et al. [40] described 14 strains of H6 virus from poultry that could infect and replicate in the lungs of mice and showed that some of these strains are lethal. Wang et al. also reported that 37 of 38 tested H6 viruses from poultry could replicate efficiently in the lungs of mice [17]; however, none were lethal. In the present study, experimental infection in mice revealed that all but one of the studied H6 strains could replicate in the lungs and turbinate of mice without adaptation. Some of the viruses also induced weight loss in mice. These data support the idea that H6 AIVs from wild birds can replicate in BALB/c mice without adaptation and that mammals can be directly infected with wild-bird AIV. Following infection, the AH/S45 strain was not found in any organs, but this strain has high sequence identity with strain AH/S20: these two strains have no reported differences in amino acid sequence. Further research is required to determine whether differences in amino acids that influence the viral replication ability exist between these two viruses. Some studies have reported that red blood cells from turkeys, sheep, and horses possess a-2, 3 receptors; however, red blood cells in humans with Type O blood and those from guinea pigs possess a-2, 6 receptors. Interestingly, in the present study, the 13 H6 AIV strains studied all exhibited high HA activity in red blood cells from horses and sheep but low activity in red blood cells from Type O humans and from guinea pigs. These results coincide with the properties bestowed by amino acids 226 and 228 in the HA gene of the virus, which together form the characteristic poultry receptor-binding site on the virus [27].

H6 AIV has a worldwide distribution, similar to H5N1 and H9N2 [41, 42], and has caused considerable losses in poultry farming. In one report, out of a total of 15,689 serum specimens gathered from 22 provinces in mainland China, approximately 63 specimens were positive for H6 AIV [12]. These cases of H6-infected humans serve as a warning. As the natural host of influenza virus, wild birds, especially migratory species, play an indispensable role in the spread and reassortment of AIV and provide unknown opportunities for the mutation and emergence of novel influenza viruses [43]. To fully understand the ecology, biology, and potential hazards of influenza virus, regular surveillance of all influenza virus hosts, and especially wild birds and poultry, is necessary. The data generated from such efforts must be combined to accurately assess the potential threat. Notably, the present study provides only a general overview of the H6 virus strains that are currently circulating among Anseriformes in Anhui, China. More intensive sampling and sequencing of AIVs from wild birds is needed to help refine our understanding of host specificity and viral evolution [44], which will also reveal the processes driving viral adaptation and maintenance in alternative hosts.


The H5, H6, and H9 subtypes of influenza A virus are found worldwide. Although it is known that H6 AIV strains found in poultry pose a potential threat to humans, few studies have surveyed this subtype in infected wild birds, and its threat to humans is therefore not well understood. In the present study, we performed a risk assessment of H6 influenza A virus transmission using phylogenetic and pathogenic analyses to better understand the potential threat of wild birds, especially migratory birds, transmitting influenza virus to poultry and/or humans. In doing so, we identified 13 strains of H6 AIV in Anseriformes. Some of these strains were reassortment from the Eurasian and North American lineages. Overall, we conclude that H6 AIV strains isolated from wild birds have the potential to infect mammals and humans.



Influenza A virus












Acidic polymerase


Polymerase basic 1


Polymerase basic 2


Reverse-transcription polymerase chain reaction


  1. Olsen B, Munster VJ, Wallensten A, Waldenström J, Osterhaus AD, Fouchier RA. Global patterns of influenza A virus in wild birds. Science. 2006;312:384–8. doi:10.1126/science.1122438.

    Article  CAS  PubMed  Google Scholar 

  2. Lewis NS, Verhagen JH, Javakhishvili Z, Russell CA, Lexmond P, Westgeest KB, et al. Influenza A virus evolution and spatio-temporal dynamics in Eurasian wild birds: a phylogenetic and phylogeographical study of whole-genome sequence data. J Gen Virol. 2015;96:2050–60. doi:10.1099/vir.0.000155.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Munster VJ, Baas C, Lexmond P, Waldenström J, Wallensten A, Fransson T, et al. Spatial, temporal, and species variation in prevalence of influenza A viruses in wild migratory birds. PLoS Pathog. 2007;3:e61. doi:10.1371/journal.ppat.0030061.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Stallknecht DE, Shane SM. Host range of avian influenza virus in free-living birds. Vet Res Commun. 1988;12:125–41. doi:10.1007/BF00362792.

    Article  CAS  PubMed  Google Scholar 

  5. Zhu X, Yu W, McBride R, Li Y, Chen LM, Donis RO, et al. Hemagglutinin homologue from H17N10 bat influenza virus exhibits divergent receptor-binding and pH-dependent fusion activities. Proc Natl Acad Sci U S A. 2013;110:1458–63. doi:10.1073/pnas.1218509110.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Fouchier RA, Munster VJ. Epidemiology of low pathogenic avian influenza viruses in wild birds. Rev Sci Tech. 2009;28:49–58. doi:10.20506/rst.28.1.1863.

  7. Tong S, Li Y, Rivailler P, Conrardy C, Castillo DA, Chen LM, et al. A distinct lineage of influenza A virus from bats. Proc Natl Acad Sci U S A. 2012;109:4269–74. doi:10.1073/pnas.1116200109.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Tong S, Zhu X, Li Y, Shi M, Zhang J, Bourgeois M, et al. New world bats harbor diverse influenza A viruses. PLoS Pathog. 2013;9:e1003657. doi:10.1371/journal.ppat.1003657.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Chen HL, Li YB, Li ZJ, Shi JZ, Shinya K, Deng GH, et al. Properties and dissemination of H5N1 viruses isolated during an influenza outbreak in migratory waterfowl in western China. J Virol. 2006;80:5976–83. doi:10.1128/JVI.00110-06.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Downie JC, Laver WG. Isolation of a type A influenza virus from an Australian pelagic bird. Virology. 1973;51:259–69. doi:10.1016/0042-6822(73)90426-1.

    Article  CAS  PubMed  Google Scholar 

  11. Shi W, Shi Y, Wu Y, Liu D, Gao GF. Origin and molecular characterization of the human-infecting H6N1 influenza virus in Taiwan. Protein Cell. 2013;4:846–53. doi:10.1007/s13238-013-3083-0.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Xin L, Bai T, Zhou JF, Chen YK, Li XD, Zhu WF, et al. Seropositivity for avian influenza H6 virus among humans, china. Emerg Infect Dis. 2015;21:1267–9. doi:10.3201/eid2107.150135.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Chen Z, Santos C, Aspelund A, Gillim-Ross L, Jin H, Kemble G, et al. Evaluation of live attenuated influenza A virus H6 vaccines in mice and ferrets. J Virol. 2009;83:65–72. doi:10.1128/JVI.01775-08.

    Article  CAS  PubMed  Google Scholar 

  14. 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. doi:10.1016/S0140-6736(14)60111-2.

    Article  PubMed  Google Scholar 

  15. Cui L, Liu D, Shi W, Pan J, Qi X, Li X, et al. Dynamic reassortments and genetic heterogeneity of the human-infecting influenza A (H7N9) virus. Nat Commun. 2014;5:3142. doi:10.1038/ncomms4142.

    Article  PubMed  Google Scholar 

  16. Liu J, Xiao H, Lei F, Zhu Q, Qin K, Zhang XW, et al. Highly pathogenic H5N1 influenza virus infection in migratory birds. Science. 2005;309:1206. doi:10.1126/science.1115273.

    Article  CAS  PubMed  Google Scholar 

  17. Wang G, Deng G, Shi J, Luo W, Zhang G, Zhang Q, et al. H6 influenza viruses pose a potential threat to human health. J Virol. 2014;88:3953–64. doi:10.1128/JVI.03292-13.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Wu H, Peng X, Peng X, Cheng L, Wu N. Molecular characterization of novel reassortant H6N2 subtype avian influenza viruses isolated from poultry in eastern China, in 2014. Infect Genet Evol. 2015;36:41–5. doi:10.1016/j.meegid.2015.08.043.

    Article  PubMed  Google Scholar 

  19. Yuan J, Zhang L, Kan X, Jiang L, Yang J, Guo Z, et al. Origin and molecular characteristics of a novel 2013 avian influenza A(H6N1) virus causing human infection in Taiwan. Clin Infect Dis. 2013;57:1367–8. doi:10.1093/cid/cit479.

    Article  CAS  PubMed  Google Scholar 

  20. Lee DH, Lee HJ, Lee YJ, Kang HM, Jeong OM, Kim MC, et al. DNA barcoding techniques for avian influenza virus surveillance in migratory bird habitats. J Wildl Dis. 2010;46:649–54. doi:10.7589/0090-3558-46.2.649.

    Article  PubMed  Google Scholar 

  21. Pedersen JC. Hemagglutination-inhibition test for avian influenza virus subtype identification and the detection and quantitation of serum antibodies to the avian influenza virus. Methods Mol Biol. 2008;436:53–66. doi:10.1007/978-1-59745-279-3_8.

    CAS  PubMed  Google Scholar 

  22. Liang L, Deng G, Shi J, Wang S, Zhang Q, Kong H, et al. Genetics, receptor binding, replication, and mammalian transmission of H4 avian influenza viruses isolated from live poultry markets in China. J Virol. 2015;90:1455–69. doi:10.1128/JVI.02692-15.

    Article  PubMed  Google Scholar 

  23. Yang H, Chen Y, Qiao C, He X, Zhou H, Sun Y, et al. Prevalence, genetics, and transmissibility in ferrets of Eurasian avian-like H1N1 swine influenza viruses. Proc Natl Acad Sci U S A. 2016;113:392–7. doi:10.1073/pnas.1522643113.

    Article  CAS  PubMed  Google Scholar 

  24. Li X, Fu Y, Yang J, Guo J, He J, Guo J, et al. Genetic and biological characterization of two novel reassortant H5N6 swine influenza viruses in mice and chickens. Infect Genet Evol. 2015;36:462–6. doi:10.1016/j.meegid.2015.08.017.

    Article  CAS  PubMed  Google Scholar 

  25. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9. doi:10.1093/molbev/mst197.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Zou S, Gao R, Zhang Y, Li X, Chen W, Bai T, et al. Molecular characterization of H6 subtype influenza viruses in southern China from 2009 to 2011. Emerg Microbes Infect. 2016;5:e73. doi:10.1038/emi.2016.71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Wang F, Qi J, Bi Y, Zhang W, Wang M, Zhang B, et al. Adaptation of avian influenza A (H6N1) virus from avian to human receptor-binding preference. EMBO J. 2015;34:1661–73. doi:10.15252/embj.201590960.

  28. Matrosovich MN, Gambaryan AS, Teneberg S, Piskarev VE, Yamnikova SS, Lvov DK, et al. Avian influenza A viruses differ from human viruses by recognition of sialyloligosaccharides and gangliosides and by a higher conservation of the HA receptor-binding site. Virology. 1997;233:224–34. doi:10.1006/viro.1997.8580.

    Article  CAS  PubMed  Google Scholar 

  29. Sun Y, Tan Y, Wei K, Sun H, Shi Y, Pu J, et al. Amino acid 316 of hemagglutinin and the neuraminidase stalk length influence virulence of H9N2 influenza virus in chickens and mice. J Virol. 2013;87:2963–8. doi:10.1128/JVI.02688-12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Hulse-Post DJ, Franks J, Boyd K, Salomon R, Hoffmann E, Yen HL, et al. Molecular changes in the polymerase genes (PA and PB1) associated with high pathogenicity of H5N1 influenza virus in Mallard ducks. J Virol. 2007;81:8515–24. doi:10.1128/JVI.00435-07.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Zamarin D, Ortigoza MB, Palese P. Influenza A virus PB1-F2 protein contributes to viral pathogenesis in mice. J Virol. 2006;80:7976–83. doi:10.1128/JVI.00415-06.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. McAuley JL, Hornung F, Boyd KL, Smith AM, McKeon R, Bennink J, et al. Expression of the 1918 influenza A virus PB1-F2 enhances the pathogenesis of viral and secondary bacterial pneumonia. Cell Host Microbe. 2007;2:240–9. doi:10.1016/j.chom.2007.09.001.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Deyde VM, Xu X, Bright RA, Shaw M, Smith CB, Zhang Y, et al. Surveillance of resistance to adamantanes among influenza A(H3N2) and A(H1N1) viruses isolated worldwide. J Infect Dis. 2007;196:249–57. doi:10.1086/518936.

    Article  CAS  PubMed  Google Scholar 

  34. Hall JS, Hallgrimsson GT, Suwannanarn K, Sreevatsen S, Ip HS, Magnusdottir E, et al. Avian influenza virus ecology in Iceland shorebirds: intercontinental reassortment and movement. Infect Genet Evol. 2014;28:130–6. doi:10.1016/j.meegid.2014.09.013.

    Article  PubMed  Google Scholar 

  35. Fislová T, Gocník M, Sládková T, Durmanová V, Rajcáni J, Varecková E, et al. Multiorgan distribution of human influenza A virus strains observed in a mouse model. Arch Virol. 2009;154:409–19. doi:10.1007/s00705-009-0318-8.

    Article  PubMed  Google Scholar 

  36. Yoon SW, Webby RJ, Webster RG. Evolution and ecology of influenza A viruses. Curr Top Microbiol Immunol. 2014;385:359–75. doi:10.1007/82_2014_396.

    PubMed  Google Scholar 

  37. Bahl J, Vijaykrishna D, Holmes EC, Smith GJ, Guan Y. Gene flow and competitive exclusion of avian influenza A virus in natural reservoir hosts. Virology. 2009;390:289–97. doi:10.1016/j.virol.2009.05.002.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Lee DH, Park JK, Yuk SS, Erdene-Ochir TO, Kwon JH, Lee JB, et al. Complete genome sequence of a natural reassortant H9N2 avian influenza virus found in bean goose (Anser fabalis): direct evidence for virus exchange between Korea and China via wild birds. Infect Genet Evol. 2014;26:250–4. doi:10.1016/j.meegid.2014.06.007.

    Article  CAS  PubMed  Google Scholar 

  39. Newman SH, Hill NJ, Spragens KA, Janies D, Voronkin IO, Prosser DJ, et al. Eco-virological approach for assessing the role of wild birds in the spread of avian influenza H5N1 along the Central Asian flyway. PLoS One. 2012;7:e30636. doi:10.1371/journal.pone.0030636.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Gillim-Ross L, Santos C, Chen Z, Aspelund A, Yang CF, Ye D, et al. Avian influenza h6 viruses productively infect and cause illness in mice and ferrets. J Virol. 2008;82:10854–63. doi:10.1128/JVI.01206-08.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Pepin KM, Wang J, Webb CT, Smith GJD, Poss M, Hudson PJ, et al. Multiannual patterns of influenza A transmission in Chinese live bird market systems. Influenza Other Respir Viruses. 2013;7:97–107. doi:10.1111/j.1750-2659.2012.00354.x.

    Article  PubMed  Google Scholar 

  42. Huang K, Bahl J, Fan XH, Vijaykrishna D, Cheung CL, Webby RJ, et al. Establishment of an H6N2 influenza virus lineage in domestic ducks in southern China. J Virol. 2010;84:6978–86. doi:10.1128/JVI.00256-10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Driskell EA, Jones CA, Stallknecht DE, Howerth EW, Tompkins SM. Avian influenza virus isolates from wild birds replicate and cause disease in a mouse model of infection. Virology. 2010;399:280–9. doi:10.1016/j.virol.2010.01.005.

    Article  CAS  PubMed  Google Scholar 

  44. Runstadler J, Hill N, Hussein ITM, Puryear W, Keogh M. Connecting the study of wild influenza with the potential for pandemic disease. Infect Genet Evol. 2013;17:162–87. doi:10.1016/j.meegid.2013.02.020.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Deng G, Shi J, Wang J, Kong H, Cui P, Zhang F, et al. Genetics, receptor binding, and virulence in mice of H10N8 influenza viruses isolated from ducks and chickens in live poultry markets in China. J Virol. 2015;89:6506–10. doi:10.1128/JVI.00017-15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


We are grateful for the technical support and staff of related laboratories in the Natural Conservation & Management Station of Anhui Province.


This research was supported by the National Key National Key Research and Development Program of China(2016YFD0500201) and the Surveillance of Wildlife Diseases Department from the State Forestry Administration of China.

Availability of data and materials

Data supporting the conclusions of this article are presented in the manuscript.

Authors’ contributions

This study is a result of collaborative work. YG participated in data collection, conceiving the study design, and drafting the manuscript. QY and SC worked on isolating viruses. XW and ZF participated in sample collection. HChen provided equipment used for the study. JM worked on data analysis. HChai participated in discussions and modified the manuscript. GD provided equipment used for the study as well as funding support. YH participated in study design, made suggestions throughout the study, revised the manuscript and provided funding. All authors read and approved the final manuscript.

Competing interests

The authors declare that their work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Consent for publication

Not applicable.

Ethics approval

This study was carried out in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the Ministry of Science and Technology of the People’s Republic of China. General animal care was provided as required by the Institutional Animal Care and Use Committee (approval number BRDW-XBS-14 for mice) [17, 45].

Author information

Authors and Affiliations


Corresponding authors

Correspondence to Guohua Deng or Hualan Chen.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ge, Y., Chai, H., Fan, Z. et al. New H6 influenza virus reassortment strains isolated from Anser fabalis in Anhui Province, China. Virol J 14, 36 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: