Phylogeny, Pathogenicity, and Transmission of H5N1 Avian Influenza Viruses in Chickens.

We analyzed five H5N1 avian influenza viruses (AIVs) isolated from different birds in 2012 in China. Based on whole-genome sequences, we divided the viruses into four genotypes. The DKE26, GSE43, and DKE53 viruses belonged to Genotypes 1-3, respectively. The CKE93 and CKE96 viruses were classified into Genotype 4. Genotypes 1-3 correspond to the viruses containing the HA gene of clade 2.3.2, and Genotype 4 is the virus that bears the HA gene of clade 7.2. To better understand the pathogenicity and transmission of the viruses, we infected chickens with 103 EID50/0.1 ml GSE43 (clade 2.3.2) or CKE93 (clade 7.2) virus. Our results revealed that 6 of 7 specific-pathogen-free (SPF) chickens inoculated with GSE43 virus were dead before 7-day post-infection, but all the SPF chickens inoculated with CKE93 virus survived the infection. Both the GSE43 and CKE93 viruses replicated systemically in chickens. The virus titers of GSE43 virus in tested organs were obviously higher than those of CKE93 virus. Our results revealed that the pathogenicity and replication of GSE43 in chickens was much higher than those of CKE93. The GSE43 virus could transmit between chickens, but the CKE93 could not transmit between chickens by naïve contact. Therefore, different clades of H5N1 AIVs possessed variable pathogenicities and transmission abilities among chickens. Our study contributes to knowledge of pathogenic variations of prevalent H5N1 viruses.

Introduction

pan class="Species">Avian influenza viruses (AIVs) are single-stranded negative-sense RNA viruses belonging to the Orthomyxoviridae family (Webster et al., 1992). In 1996, the first n>n class="Species">H5N1 AIVs were found in sick geese in Guangdong, and the H5N1 viruses were divided into highly pathogenic avian influenza viruses (HPAIVs) and low pathogenic avian influenza viruses (LPAIVs) (Li et al., 2006). The H5N1 HPAIVs infected 18 humans in Hong Kong in 1997, six of whom died (Subbarao et al., 1998). Till February 14, 2017, the World Health Organization (WHO) reported that there have been 856 humans infected with H5N1 HPAIVs in 16 countries; 452 of these patients have died (WHO, 2017). Since the beginning of 2004, there had been significant outbreaks of H5N1 AIVs infection involving multiple farm flocks in more than 20 provinces in China (Chen, 2009). It was clear that H5N1 HPAIVs could cause enormous economic losses and posed a serious public health threat.

pan class="Species">H5N1 Hn>n class="Chemical">PAIVs could be perpetuated in birds or mammals. Although the natural reservoir for AIVs was thought to be waterfowls, some viruses were highly pathogenic to domestic poultry (Chen et al., 2006). Moreover, H5N1 HPAIVs could be transmitted from waterfowls to mammalian and domestic poultry (Webster et al., 1992). In China, the H5N1 HPAIVs transmitted between aquatic birds and domestic poultry had contributed to the genetic diversity of the circulated viruses in domestic poultry (Duan et al., 2008; Vijaykrishna et al., 2008). H5N1 HPAIVs were classified into 10 clades (0–9) and several second clades by H5 HA nomenclature (World Health Organization/World Organization for Animal Health/Food and Agriculture Organization (WHO/OIE/FAO) H5N1 Evolution Working Group, 2014). The H5N1 viruses including clade 2.2, 2.3.2, 2.3.4, 4, 7, and 9 had been co-circulating from 2005 to 2006 in China, while the circulated viruses in domestic birds and waterfowls during 2007–2009 in China belonged to clades 2.3.2, 2.3.4, and 7 (Jiang et al., 2010). Several previous studies suggested that the pathogenicity of clade 2.3.2 viruses to waterfowls were increasing (Sakoda et al., 2010). Since 2005, H5N1 HPAIVs of clade 7 have occurred in chickens in northern China. The clade 7.2 viruses broke out in China in 2008 (e.g., in the Ningxia and Jiangsu regions), and caused a significant number of chicken deaths (WHO, 2012). To understand the transmission and pathogenicity of the H5N1 viruses in chicken, we chose five viruses belonging to various clades isolated from different birds.

Materials and methods

Viruses

The five pan class="Species">H5N1 Hn>n class="Chemical">PAIVs A/Duck/China/E26/2012(H5N1) (DKE26), A/Goose/China/E43/2012(H5N1) (GSE43), A/Duck/China/E53/2012(H5N1) (DKE53), A/Chicken/China/E93/2012(H5N1) (CKE93), and A/Chicken/China/E96/2012(H5N1) (CKE96) were obtained from swabs of different poultry in markets in 2012. Nine- to ten-days-old specific-pathogen-free (SPF) embryonated chicken eggs were used to purify and propagate the swabs by three rounds of limiting dilution method (Jiao et al., 2014; Yuan et al., 2014). The viruses were collected from allantoic fluids of multiple eggs. Values of 50% egg infective doses (EID50) was calculated using the Behrens-Reed-Muench cumulant method (Thakur and Fezio, 1981). Virus isolation and purification were done in biosafety level 3 (BSL-3) facilities.

Sequence analysis

The genome of the viruses was sequenced. pan class="Chemical">Trizol LS Reagent (Invitrogen Life Technologies, Carlsbad, CA, USA) was used to extract viral RNA from allantoic fluid and reverse transcription was performed with MLV (Invitrogen Life Technologies, Carlsbad, CA, USA). PCR amplification used n>n class="Chemical">pfx (Invitrogen Life Technologies, Carlsbad, CA, USA). QIAquick PCR purification kit (Qiagen, Valencia, CA, USA) was used to purify the PCR products and automatic ABI Prism 3730 genetic analyzer (Applied Biosystems, Foster City, CA, USA) was used to sequence. Lasergene 7.1 (DNASTAR, Madison, WI, USA) was used to compile and edit the DNA sequences. The phylogenetic tree was generated by the distance-based neighbor-joining method using MEGA 5 (Sinauer Associates, Inc., Sunderland, MA, USA). Bootstrap analysis with 1,000 replicates was used to assess the reliability of the trees. The genetic distance and horizontal distances were proportional. The nucleotide sequences in this study are available from NCBI GenBank (MF116309–MF116348).

Pathogenicity and transmission

Four-week-old pan class="Disease">SPF White Leghorn n>n class="Species">chickens (n = 20) were allocated into two groups of 10 animals per group and housed in isolator cages. Seven chickens of the GSE43 or CKE93 group were inoculated intranasally with 103 EID50/0.1 ml of GSE43 or CKE93 virus, respectively (Yuan et al., 2014). Contact group cohabited with those animals inoculated with the GSE43 or CKE93 virus were three chickens inoculated with 0.1 ml phosphate buffered saline (PBS). Significance of differences in survival was tested by Log-Rank analysis. Clinical symptoms of the chickens were observed for 14 days. Viral replication in the lungs, kidneys, brain, heart, spleen and liver of three infected chickens in each group were detected on day post-infection (DPI) 3. So were done on dead chickens. On 3, 5, 7, 9, and 11 DPI, cloacal and oropharyngeal swabs were collected from chickens and suspended in 1 ml PBS. All of the tissues and swabs were collected and titrated for virus infectivity in eggs. In accordance with the biosafety committee of South China Agriculture University protocols, the animal infected experiments were done in ABSL-3 facilities. The animal handling was done by the experimental animal administration and the ethics committee of South China Agriculture University guidelines.

Results

Phylogenetic analysis of the H5N1 viruses

The genome of the viruses was sequenced to determine the molecular characteristic of them. We compared their sequences and those of the representative pan class="Species">H5N1 viruses got from NCBI GenBank (Supplementary Table 5). Based on antigenic characteristics by the WHO, the HA genes of DKE26, n>n class="Chemical">GSE43, and DKE53 were classified to clade 2.3.2. The DKE53 could further be divided into clade 2.3.2.1.A. The DKE26 could further belong to 2.3.2.1.C. The GSE43 could further be classified into 2.3.2.1.B. The HA gene of other two viruses (CKE93 and CKE96) belonged to clade 7.2 (Figure 1). Compared with the GSGD1/96 virus nucleotide, the DKE26, GSE43, DKE53, CKE93, and CKE96 nucleotide similarities were 92.1, 91.3, 92.8, 91.8, and 91.6, respectively. We divided the HA genes of our isolates phylogenetically into three groups. The HA genes of the DKE26 and DKE53 viruses belonged to the MG-like group, and the HA gene of the GSE43 virus was classified into the GX-like group. The HA genes of the CKE93 and CKE96 viruses were clustered into the SX-like group. The HA genes similarity in groups were over 97%, and between the groups were <96%.

Figure 1

Phylogenetic analysis of HA. HA(A): nucleotides (nt) 29–1,732. Triangles indicate viruses characterized in this study, other virus sequence were downloaded from GenBank. QH, Qinghai; VN, Vietnam; IDN, Indonesia; GY, Guiyang; MG, Mongolia; GX, Guangxi; HD, Huadong; SX, Shanxi.

We retrieved complete NA sequences from NCBI GenBank to compare to the viruses containing the N1 gene. The NA genes of these viruses were divided into three groups, and the similarity among the three groups was <97%. The NA gene of the DKE26 virus belonged to the MG-like group. The NA genes of the pan class="Chemical">GSE43 and DKE53 viruses were divided into the GX-like group. The NA genes of the CKE93 and CKE96 viruses were clustered to the SX-like group (Figure 2).

Figure 2

Phylogenetic analysis of NA. NA: nt19–1,380. Except our isolate, other virus sequences were downloaded from GenBank. QH, Qinghai; VN, Vietnam; IDN, Indonesia; GY, Guiyang; MG, Mongolia; GX, Guangxi; HD, Huadong; SX, Shanxi.

The n class="Chemical">PA genes of these viruses were clustered into three groups. The MG-like group included of DKE26 virus, and the GX-like group contained n>n class="Chemical">GSE43 and DKE53 viruses. The CKE93 and CKE96 viruses were classified into the SX-like group. The PA genes of MG-like group viruses were <94% similarity with the GX-like group and the SX-like group viruses (Figure 3).

Figure 3

Phylogenetic analysis of PA. PA: nt 25–2,175. Except our isolate, other virus sequences were downloaded from GenBank. QH, Qinghai; VN, Vietnam; IDN, Indonesia; GY, Guiyang; MG, Mongolia; GX, Guangxi; HD, Huadong; SX, Shanxi.

We classified the pan class="Gene">PB1 genes of these n>n class="Species">H5N1 viruses into three groups. The PB1 gene of the DKE26 virus clustered to the MG-like group. The PB1 genes of the GSE43 and DKE53 viruses belonged to the GX-like group. The PB1 genes of the CKE93 and CKE96 viruses were classified into the SX-like group. The PB1 genes of SX-like group viruses were <95% similarity with other three groups, and the similarity among the SX-like group and the MG-like group ranged between 96 and 98% (Figure 4).

Figure 4

Phylogenetic analysis of PB1. PB1: nt 25–2,298. Except our isolate, other virus sequences were downloaded from GenBank. QH, Qinghai; VN, Vietnam; IDN, Indonesia; GY, Guiyang; MG, Mongolia; GX, Guangxi; HD, Huadong; SX, Shanxi.

Phylogenetic analysis of pan class="Gene">PB1. n>n class="Gene">PB1: nt 25–2,298. Except our isolate, other virus sequences were downloaded from GenBank. QH, Qinghai; VN, Vietnam; IDN, Indonesia; GY, Guiyang; MG, Mongolia; GX, Guangxi; HD, Huadong; SX, Shanxi.

The PB2 genes of these viruses consisted of three groups. The MG-like group included the PB2 gene of the DKE26 virus. The GX-like group contained the PB2 gene of the pan class="Chemical">GSE43 and DKE53 viruses. The PB2 genes of the n>n class="Chemical">CKE93 and CKE96 viruses belonged to the SX-like group. The PB2 genes of the viruses among the three groups shared <93% similarity with one other (Figure 5).

Figure 5

Phylogenetic analysis of PB2. PB2: nt 28–230. Except our isolate, other virus sequences were downloaded from GenBank. QH, Qinghai; VN, Vietnam; IDN, Indonesia; GY, Guiyang; MG, Mongolia; GX, Guangxi; HD, Huadong; SX, Shanxi.

We divided the NP genes of these viruses into two groups. The NP genes of the DKE26, pan class="Chemical">GSE43, and DKE53 viruses derived from GX-like viruses. The NP genes of the n>n class="Chemical">CKE93 and CKE96 viruses were clustered into GY-like viruses. The NP genes of the viruses in the GX-like group shared <95% similarity with those viruses in the GY-like group (Figure 6).

Figure 6

Phylogenetic analysis of NP. NP: nt46–1,542. Except our isolate, other virus sequences were downloaded from GenBank. QH, Qinghai; VN, Vietnam; IDN, Indonesia; GY, Guiyang; MG, Mongolia; GX, Guangxi; HD, Huadong; SX, Shanxi.

The M and NS genes of the viruses were divided into three groups. The M and NS genes of the DKE26 virus belonged to the MG-like group. The M and NS genes of the pan class="Chemical">GSE43 and DKE53 viruses were clustered into the GX-like group. The M and NS genes of the n>n class="Chemical">CKE93 and CKE96 viruses were classified into the SX-like viruses, respectively (Figures 7, 8).

Figure 7

Phylogenetic analysis of M. M: nt26–795. Except our isolate, other virus sequences were downloaded from GenBank. QH, Qinghai; VN, Vietnam; IDN, Indonesia; GY, Guiyang; MG, Mongolia; GX, Guangxi; HD, Huadong; SX, Shanxi.

Figure 8

Phylogenetic analysis of NS. NS: nt27–678. Except our isolate, other virus sequences were downloaded from GenBank. QH, Qinghai; VN, Vietnam; IDN, Indonesia; GY, Guiyang; MG, Mongolia; GX, Guangxi; HD, Huadong; SX, Shanxi.

Based on genomic diversity, we here classified the viruses into four genotypes (Table 1). The pan class="Species">H5N1 viruses of Genotypes 1–3 were isolated from Waterfowl, and Genotype 4 contained viruses from n>n class="Species">chickens. Genotypes 1–3 contained the HA gene of clade 2.3.2, and Genotype 4 contained the virus that bore the HA gene of clade 7.2, which was originally detected in chickens in China.

Table 1

Genotypic evolution of the H5N1 viruses isolated in China in 2012.

Virus	Group of each gene segment in the phylogenetic treea								Genotype	Species
	HA	NA	PB2	PB1	PA	NP	M	NS
A/DK/CN/E26/2012	MG	MG	MG	MG	MG	GX	MG	MG	1	DK
A/GS/CN/E43/2012	GX	GX	GX	GX	GX	GX	GX	GX	2	GS
A/DK/CN/E53/2012	MG	GX	GX	GX	GX	GX	GX	GX	3	DK
A/CK/CN/E93/2012	SX	SX	SX	HD	SX	GY	SX	SX	4	CK
A/CK/CN/E96/2012	SX	SX	SX	HD	SX	GY	SX	SX	4	CK

The eight gene segments are indicated at the top of each bar. The number in each bar shows the group of genes indicated in Figures .

Genotypic evolution of the pan class="Species">H5N1 viruses isolated in China in 2012.

Molecular characterization of the H5N1 viruses

A series of pan class="Chemical">basic amino acids at the cleavage site of the HA (–RRRKR↓G–) were found in all the five viruses, which was characteristic of Hn>n class="Chemical">PAIVs (Nobusawa et al., 1991; Belser et al., 2009). The amino acid residues Q226 and G228 in HA showed that the viruses preferentially bound to the AIVs receptor (Ha et al., 2001). We found three potential N-linked glycosylation sites in HA1 (27, 39, and 301 or 302) and two in HA2 (499 or 500 and 558 or 559) in the five viruses. An extra potential N-linked glycosylation site in HA1 (181) was found in DKE26, GSE43, and DKE53 viruses, and another two extra potential N-linked glycosylation site in HA1 (156 and 209) were found in the GSE43 virus. Five extra potential N-linked glycosylation sites in HA1 (88, 142, 155, 178, and 251) were found in CKE93 and CKE96 viruses (Table 2).

Table 2

Cleavage site and potential glycosylation sites in HA of the five H5N1 HPAIVs.

Strains	Cleavage site	Potential glycosylation sites
	341–361/342–347	27	39	88	142	155	156	178	181	209	251	288/289	301/302	499/500	558/559
	−RRRKR/G–	NST	NVT	NVS	NTS	NPS	NSS	NYT	NNT	NPT	NDT	NCS	NSS	NGT	NGS
DKE26	+a	+	+	−b	−	−	−	−	+	−	−	−	+	+	+
GSE43	+	+	+	−	−	−	−	−	+	−	−	−	+	+	+
DKE53	+	+	+	−	−	−	+	−	+	+	−	−	+	+	+
CKE93	+	+	+	+	+	+	−	+	−	−	+	−	+	+	+
CKE96	+	+	+	+	+	+	−	+	−	−	+	−	+	+	+

The “+” means the amino acid sequences of glycosylation sites are same with list above.

The “−” means the glycosylation sites are lost.

Cleavage site and potential glycosylation sites in HA of the five pan class="Species">H5N1 Hn>n class="Chemical">PAIVs.

The NA genes of five viruses had 20 amino acid deletions at positions 49–68. The mutations pan class="Mutation">H274Y and n>n class="Mutation">N294S were not found in the viruses, which indicated the absence of antiviral drug-resistant residues (Gubareva et al., 2001).

pan class="Chemical">Amantadine and n>n class="Chemical">rimantadine target the M2 protein, and single mutations in the trans-membrane domain of M2 (e.g., residues L26F, V27A/T, A30T/V, S31N/R, and G34E) could confer resistance to these drugs (Suzuki et al., 2003). An S31N mutation was observed in the M2 gene of CKE93 and CKE96 viruses, suggesting that the CKE93 and CKE96 viruses were not sensitive to this class of antiviral drugs. No amino acid substitutions were found in other residues.

The pan class="Mutation">S314N mutation in the NP of n>n class="Species">H5N1 virus could cause a defect in nuclear localization at high temperature (Siboonnan et al., 2013). The mutation N319K in NP might influence pathogenicity (Gabriel et al., 2005), but no amino acid substitutions were found in these residues in the viruses.

The pan class="Gene">NS1 genes of the DKE26, n>n class="Chemical">GSE43, DKE53, CKE93, and CKE96 viruses had a five-amino-acid deletion at positions 80–84, which may contribute to increased virulence (Long et al., 2008). Previous studies had shown that P42S mutations may contribute to the virulence of H5N1 viruses in mice (Jiao et al., 2008); these mutations were observed the NS1 genes of the five viruses. An importation virulence factor of influenza A viruses was the PDZ-binding motif of NS1 (Liu et al., 2010). The sequences of the DKE26, DKE53, CKE93, and CKE96 viruses in the PDZ domain were ESEV, which also was found in the NS1 of 1918 pandemic virus and 1997–2003 H5N1 viruses (Jackson et al., 2008).

Previous study showed that the mutation pan class="Mutation">T515A in n>n class="Chemical">PA gene contributes to the viruses transmission in ducks (Hulse-Post et al., 2007). The amino acid 224P in PA increased the replication of the virus in duck embryo fibroblasts (Song et al., 2011). The absence of the S224P and T515A mutations in the PA protein suggests poor replication of those viruses in ducks. The mutation R185K and C241Y in PA gene enhanced the growth capability of the viruses in human cells (Fan et al., 2014; Yamaji et al., 2015). The presence of 185R and 241C in the PA protein showed low growth capability of the viruses in human cell.

The presence of pan class="Mutation">L13P and n>n class="Mutation">S678N mutations in the PB1 protein increases mammalian pathogenicity (Zell et al., 2007). Although the L13P mutation was present in the five isolates, the S678N mutation was not found, which indicates that those viruses did not acquire increased mammalian pathogenicity. A 90-residue PB1-F2 was observed in DKE26, DKE53, CKE93, and CKE96 viruses, but the GSE43 virus encoded a 57aa PB1-F2 protein.

The pan class="Mutation">E627K substitution was absent in the PB2 protein of the n>n class="Chemical">GSE43, DKE53, CKE93, and CKE96 viruses, but it was present in the DKE26 virus. The presence of 701D and 714S in the PB2 protein of the five viruses suggested that the virus might not be capable of increased virulence in mammalian species.

Pathogenicity of H5N1 HPAIVs in chickens

The pan class="Chemical">GSE43 of clade 2.3.2.1B and n>n class="Chemical">CKE93 of clade 7.2 viruses were obtained from swabs of gooses and chickens with apparent clinical symptoms. We focused on the two viruses to assess the pathogenicity of the H5N1 HPAIVs in chickens. SPF chickens were infected i.n. at 103 EID50. The Mean Death Time (MDT) of SPF chickens inoculated with GSE43 was 4.6 days (83.3% mortality) (Figure 9). No chicken inoculated with CKE93 virus died during the observation period (Figure 9). Torticollis and neurological symptoms were observed among SPF chickens infected with GSE43. Only one SPF chicken with slight torticollis was observed in the CKE93-inoculated group. Therefore, the two viruses had various pathogenicities in chickens.

Figure 9

Lethality of the A/GS/CN/E43/2012 and A/CK/CN/E93/2012 viruses in SPF chicken. Death patterns of the SPF chicken infected with A/GS/CN/E43/2012 and A/CK/CN/E93/2012 viruses with the doses of 103 EID50. Curves are significantly different (P < 0.05) by log-rank analysis.

Lethality of the A/GS/CN/E43/2012 and A/CK/CN/E93/2012 viruses in SPF pan class="Species">chicken. Death patterns of the SPF n>n class="Species">chicken infected with A/GS/CN/E43/2012 and A/CK/CN/E93/2012 viruses with the doses of 103 EID50. Curves are significantly different (P < 0.05) by log-rank analysis.

The pan class="Chemical">GSE43 and n>n class="Chemical">CKE93 viruses replicated systemically in chickens, which could be found from all of the tested organs—including the heart, liver, spleen, lungs, kidneys, and brain—3 DPI. The virus titers of GSE43 virus in tested organs were obviously higher than those of CKE93 virus. The mean virus titers of the GSE43 in the heart, liver, spleen, lungs, kidneys and brain were 6.63, 8.63, 8.63, 7.75, 8.63, and 8.63 log10EID50, respectively. The virus titer of CKE93 virus in the heart, liver, spleen, lungs, kidneys and brain were 3.67, 3.00, 3.67, 4.00, 3.75, and 3.67 log10EID50, respectively (Table 3). All in all, the replication of GSE43 in chickens was much higher than that of CKE93.

Table 3

Replication and lethality in chickens of the H5N1 viruses after inoculated intranasally.

Strains	Clade	Titer (log10EID50)		Virus replication on 3 DPI (log10EID50/0.1ml)[b] in
				Heart	Liver	Spleen	Lung	Kidney	Brain
GSE43	Clade 2.3.2.1.B	6.8	Infectedc	6.63	8.63	8.63	7.75	8.63	8.63
			Contactd	6.63	8.63	8.63	8.00 ± 1.41	7.81 ± 1.15	7.81 ± 1.15
CKE93	Clade 7.2	7.3	Infectedc	3.67	3.00	3.67	4.00	3.75	3.67
			Contactd	–	–	–	–	–	–

Four-week-old SPF chickens were inoculated intranasally (i.n.) with 10.

For statistical analysis, a value of 1.5 was assigned if the virus was not detected from the undiluted sample in three embryonated hen eggs. Virus titers are expressed as means standard deviation in log10EID.

Ducks inoculated with virus.

Naive contact ducks housed with those inoculated.

Replication and lethality in pan class="Species">chickens of the n>n class="Species">H5N1 viruses after inoculated intranasally.

Four-week-old SPF pan class="Species">chickens were inoculated intranasally (i.n.) with 10.

pan class="Species">Ducks inoculated with virus.

Naive contact pan class="Species">ducks housed with those inoculated.

Oropharyngeal and cloacal swabs collected on 3, 5, 7, 9, and 11 pan class="Chemical">DPI were used to detect viruses shedding from the inoculated n>n class="Species">chickens. The GSE43 virus could not be tested from oropharyngeal swabs, but it was detected from cloacal swabs 3 DPI. The GSE43 virus was not detected after 3 DPI. The virus shedding of CKE93 were not detected from swabs of inoculated SPF chickens (Table 4). In conclusion, the virus shedding time of chickens inoculated with GSE43 was longer than that of chickens infected with CKE93.

Table 4

Virus titers in cloacal and oropharyngeal swabs from SPF chicken.

Strains		Days post-inoculation (log10EID50/0.1 ml) ± SDa
		3 Day		5 Day		7 Day		9 Day		11 Day
		Oropharyngeal swabs	Cloacal swabs	Oropharyngeal swabs	Cloacal swabs	Oropharyngeal swabs	Cloacal swabs	Oropharyngeal swabs	Cloacal swabs	Oropharyngeal swabs	Cloacal swabs
GSE43	Infectedb	0/3	1/3	0/1	0/1	0/1	0/1	0/1	0/1	–	–
	Contactc	1/3	0/3	1/2	1/2	0/1	0/1	0/1	0/1	0/1	0/1
CKE93	Infectedb	0/6	0/6	0/3	0/3	0/3	0/3	0/3	0/3	0/3	0/3
	Contactc	0/3	0/3	0/3	0/3	0/3	0/3	0/3	0/3	0/3	0/3

For statistical purposes, a value of 1.5 was assigned if virus was not detected from the undiluted sample in three embryonated hen's eggs (Sun et al., .

Chickens inoculated with virus.

Naive contact chickens housed with those inoculated.

Virus titers in cloacal and oropharyngeal swabs from SPF pan class="Species">chicken.

pan class="Species">Chickens inoculated with virus.

Naive contact pan class="Species">chickens housed with those inoculated.

Transmission of H5N1 HPAIVs in chickens

To understand the horizontal transmission of these two viruses, three SPF pan class="Species">chickens were inoculated intranasally with 0.1 ml n>n class="Chemical">PBS as a naïve-contact group. These animals were then housed with chickens inoculated with the GSE43 and CKE93 viruses. We noted that two SPF chickens in the naïve-contact group, housed with inoculated GSE43 chickens, died 3 and 5 DPI. The last naïve-contact chicken seroconverted 14 DPI. The GSE43 virus could be detected from oropharyngeal swabs of naïve-contact chickens 3 and 5 DPI and from cloacal swabs 5 DPI. The GSE43 virus also replicated systemically in naïve-contact chickens; the mean virus titers in the heart, liver, spleen, lungs, kidneys and brain were 6.63, 8.63, 8.63, 8.00 ± 1.41, 7.81 ± 1.15, and 7.81 ± 1.15 log10EID50, respectively. All of the three naïve-contact chickens housed with CKE93 survived and virus shedding was not detected from swabs. The CKE93 virus was not detected in any tested organs in naïve-contact SPF chickens (Table 3).

Therefore, we found that the pan class="Chemical">GSE43 virus could transmit between n>n class="Species">chickens by naïve contact. However, the CKE93 virus could not be transmitted to SPF chickens via naïve contact.

Discussion

Since late 2003, thousands of wild birds and millions of poultry death have been caused by pan class="Species">H5N1 Hn>n class="Chemical">PAI viruses, especially in southern China and Southeast Asia. Despite substantial efforts to control infection rates in poultry, H5N1 HPAI has broken out in some Chinese provinces. Because the wild birds migrate and the domestic poultry move, the viruses might spread at any time. In 2004, the first clade 2.3.2 virus was isolated from a dead Chinese pond heron in Hong Kong; it is now widely distributed in Asia and in Eastern Europe. Therefore, viruses of clade 2.3.2 have spread geographically and evolved genetically (WHO, 2012). Viruses of clade 2.3.2.1 have been circulating widely in China since 2010 and may have been responsible for a new wave of cross-continental spreading from Asia to Europe (Smith et al., 2009; Jiang et al., 2010; Li et al., 2010). This clade was already widely distributed in chickens, ducks, geese, and wild birds (WHO, 2015). Since 2003, the clade 7 H5N1 HPAIVs has been spreading widely in Northern China. The viruses of clade 7.2 have been spreading in China and Vietnam since 2008 (WHO, 2011). From 2010 to 2013, the clade 2.3.2.1 and 7.2 H5N1 viruses were dominant in Southern China and other clades had occasionally been found (WHO, 2013). The circulated viruses had been responsible for significantly damaging the poultry industry. Therefore, it is necessary to understand the phylogeny, pathogenicity and transmission of the viruses from clade 2.3.2 and 7.2.

The previous reported that the clade 2.3.2 AIVs had reasserted with H9 viruses and the reassortment between clades 2.3.2.1c, 2.3.2.1b, and 2.3.2.1a virus had been found (Marinova-Petkova et al., 2016; Nguyen et al., 2016). Our results showed The HA gene of the Genotypes 1 virus was same as that of Genotype 3, but the other seven genes were different from Genotypes 2 and 3. The NA, pan class="Chemical">PA, n>n class="Gene">PB1, PB2, NP, M, and NS genes of the Genotype 2 and 3 viruses were quite similar, but their HA genes exhibited significant diversity (Table 1). So the Genotype 3 virus may be a reassortant from Genotype 1 and 2 viruses. The SX-like viruses were first detected in chickens from the Shanxi province in Northern China, but the origin of these viruses remains unclear. However, genomic analyses confirmed that the SX-like viruses reassorting with new HA, NA, and PB1 genes were newly introduced into poultry in China (Li et al., 2010). Now the clade 7.2 viruses had derived gene segments from viruses of clade 2.3.4 or H9N2 (Liu et al., 2016). Our results showed that the HA, NA, M, NS, PA, and PB1 genes of Genotypes 4 viruses belonged to SX-like group viruses, but the NP and PB1 genes belonged to GY-like group viruses (clade 2.3.3) and HD-like group viruses, respectively. So the clade 7.2 viruses may had been reassorted with clade 2.3.3 viruses. This finding showed that multiple clades viruses cocirculating in waterfowl may promote the reassortment of the viruses.

Previous studies have shown that both clade 2.3.2 and 7.2 viruses were highly pathogenic to pan class="Species">chickens, and could transmit to the naive contact n>n class="Species">chickens (Yuan et al., 2014; Jiao et al., 2016). In our study, the GSE43 (clade 2.3.2) and CKE93 (clade 7.2) viruses from the different birds in 2012 were evaluated the pathogenicity of H5N1 HPAIVs in chickens. They could replicate systemically in chickens, but the virus titer of GSE43 virus was higher than that of the CKE93 virus. The pathogenicity of the CKE93 virus was much lower than that of the GSE43 virus. Especially, only the GSE43 virus could be transmitted between SPF chickens by naïve contact. Therefore, the H5N1 HPAIVs of clade 2.3.2 obtained from waterfowls were highly pathogenic to chickens and replicated systemically in chickens. And we found that the clade 7.2 virus could not be transmitted in chickens when the animals were inoculated with a low dose.

In conclusion, our findings revealed that clade 2.3.2 and 7.2 pan class="Species">H5N1 viruses had varying levels of pathogenicity and transmission in n>n class="Species">chickens. Therefore, we should intensify virological surveillance of the H5N1 viruses to understand the antigenic and pathogenic variations of prevalent H5N1 viruses.

Author contributions

JC and PJ designed this study, performed the experiment sand participated in the data collection and analysis. NQ, YG, LC, SW, KM, and YL assisted with animal experiment. JC and PJ drafted the manuscript. ML, HS, and ZQ participated in writing the discussion.

Conflict of interest statement

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