Characterization of emerging Newcastle disease virus isolates in China

Newcastle disease (ND) is a devastating worldwide disease of poultry characterized by increased respiration, circulatory disturbances, hemorrhagic enteritis, and nervous signs. Sequence analysis shows several amino acid residue substitutions at neutralizing epitopes on the F and HN proteins of recent Shaanxi strains. Both Cross protection and cross serum neutralization tests revealed that the traditional vaccine strains were unable to provide full protection for the flocks. To better understand the epidemiology of Newcastle disease outbreak, a portion of the F gene and the full-length HN gene were amplified from Shaanxi isolates by reverse transcription-polymerase chain reaction (RT-PCR) and then conducted sequence and phylogenetic analyzes. In pathogenicity analysis, both high intra-cerebral pathogenicity index (ICPI) and mean death time (MDT) tests of chicken embryo were carried out. Furthermore, a cross-protection experiment in which specific-pathogen-free chickens vaccinated with a LaSota vaccine strain were challenged by the recent Shaanxi strain was also performed. Nine Newcastle disease (ND) virus (NDV) isolates which were recovered from ND outbreaks in chicken flocks in China were genotypically and pathotypically characterized. Amino acid sequence analysis revealed that all the recent Shaanxi-isolated NDVs have 112R-R-Q-K-R-F117 for the C-terminus of the F2 protein and exhibit high ICPI and MDT of chicken embryos, suggesting that they were all classified as velogenic type of NDVs. Phylogenetic analysis of these isolates showed that they belong to subgenotype VIId that have been implicated in the recent outbreaks in northwestern China. The percentage of amino acid sequence identity of F protein between recent Shaanxi stains and five vaccine strains was in the range of 81.9 %–88.1 %, while the percentage of amino acid sequence identity of HN protein between recent Shaanxi strains and vaccine strains was in the range of 87.4 %–91.2 %. Furthermore, a number of amino acid residue substitutions at neutralizing epitopes on the F and HN proteins of these isolates were observed, which may lead to the change of antibody recognition and neutralization capacity. A cross-protection experiment indicated that specific-pathogen-free chickens vaccinated with a LaSota vaccine strain was not capable of providing full protection for the flocks that were challenged by the recent Shaanxi strain. Taken together, our findings reveal that recent Shannxi NDVstrains exhibit antigenic variations that could be responsible for recent outbreaks of NDVs in northwestern China.


Background
Newcastle disease is a highly contagious disease, which caused by the virulent strains of Newcastle disease virus (NDV), and causes great economic loss to the global poultry industry. In addition to poultry species, at least 250 species of birds can be naturally or experimentally infected with NDV [1]. The major signs of the disease include respiratory distress, diarrhea, circulatory disturbances, and central nervous system impairment [1][2][3]. Although NDV has only 1 serotype, ND outbreaks still occur frequently in flocks vaccinated with NDV live vaccines (e.g. Lasota and Clone 30) and inactivated vaccines in China. Immune failure is considered to be related to the genetic variation of NDV epidemic strains [3][4][5][6].
On the basis of restriction fragment length polymorphism (RFLP) analysis and partial nucleotide sequencing of the cleavage recognition sequence of the F gene, NDV strains are grouped into at least 10 genogroups, among which genogroup VII can be further divided into several sub-genotypes [5,7,8]. The currently circulating NDVs in the worldwide exhibit multiple genotypes and high genetic diversity. According to pathogenicity assays in specific pathogen free (SPF) chickens, NDVs have been categorized into velogenic, mesogenic, and lentogenic pathotypes. The F protein mediates the fusion of viral and cellular membranes during penetration and spread between infected and adjacent cells. The virulence of NDV can be predicted by analyzing the cleavage site sequence of F protein [1,9,10]. By comparing the amino acid residues of the C-terminal HN extension of NDV, they were divided into at least 3 groups, including the numbers of amino acid residues of 571, 577 and 616 [11,12]. The viral HN protein mediates attachment to sialic acid-containing receptor(s) through its neuraminidase (NA) activity. HN is a type II membrane glycoprotein, which exists on the surface of virions and infected cells as a tetrameric spike. The ectodomain of the HN glycoprotein consists of a membrane-proximal, stalk-like segment supporting a terminal globular domain. The antigenic, receptor recognition, and NA active sites all reside in the latter [13]. Previous studies have suggested that there are at least five antigenic sites related to epitopes on the HN protein of NDV, which are related to antibody recognition [13,14].
In this study, nine NDV strains have been isolated from NDV-infected chickens in Shaanxi province in China. To investigate the epidemiology of Newcastle disease outbreak, a portion of the F gene and the fulllength HN gene were amplified from these isolates by reverse transcription-polymerase chain reaction (RT-PCR) and then conducted sequence analysis. The nucleotide sequences of the F and HN genes of NDVs from recent Shaanxi isolates were compared with published sequences to clarify the evolution, epidemiology and their pathotypes of NDV. Our findings reveal that all Shaanxi NDV strains belong to genotype VIId. Sequence analysis showed several amino acid residue substitutions at neutralizing epitopes on the F and HN proteins of recent isolated Shaanxi strains. Furthermore, cross protection and cross serum neutralization tests suggested that the traditional vaccine strains such as LaSota were unable to provide full protection for the flocks.  (Table 1). All Shaanxi strains of NDV with MDT of 38-50 h and with ICPI of 1.67-1.78 (Table 1) were classified as virulent NDVs. By using the primers F1 and F2, the fragments of the F gene were amplified from all Shaanxi NDV isolates with the expected size of 957 bp in length. By using the primer pairs (HN1-HN2 and HN3-HN4), 2 PCR products (1002 and 1338 bp) covered the full-length HN gene (1716 bp) could be amplified from all Shaanxi NDV isolates (data not shown).

Nucleotide and amino acid sequence analysis of the F gene of NDVs
All the recent Shaanxi strains contain the amino acid sequences 112 R-R-Q-K-R-F 117 for the C-terminus of the F2 protein, indicating that they are velogenic type of NDV. Furthermore, all the recent Shaanxi strains possess 101 K and 121 V (Table 2), a characteristic of genotype VII viruses reported previously [15].
The nucleotide and predicted amino acid sequences of the F gene of NDVs were aligned. The aligned sequences were then used in pairwise comparisons and the percentage of nucleotide and amino acid sequence identities were determined.
The percentage of nucleotide sequence identity of F gene among these Shaanxi strains was in the range of 98.9 %-99.6 %, while the percentage of amino acid sequence identity was in the range of 98.4 %-99.5 %. Sequence comparisons of Shaanxi strains with five vaccine strains (Herts/33, B1, Clone 30, LaSota, and Mukteswar) homologues revealed that there exist extensive sequence divergence, with 80.7 %-87.2 % and 81.9 %-88.1 % identities of nucleotide and amino acid sequence, respectively. The Shaanxi strains exhibited higher sequence identity with Herts/33 vaccine strain. The percentage of nucleotide sequence identity was in the range of 86.6 %-87.2 %, while the percentage of amino acid sequence identity was in the range of 87.6 %-88.1 %. Furthermore, Shaanxi strains Comparison of the deduced amino acid sequences of the F protein of recent Shaanxi strains with the traditional vaccine strains revealed that there were many amino acid residue substitutions at F protein (Table 2). Interestingly, compared to the past pandemic strains in China during the past 10 years [16][17][18][19][20], a number of amino acid residue substitutions at N-terminal of F protein were found in Shaanxi NDV strains. Examples include amino acid residue substitutions of all NDV at positions 13 (M to L and M to P), 14 (M to V), 19 (I to V or T to V), 22 (V to A, V to T and V to I), 28 (P to L, P to T and P to S), and 51 (N to S) ( Table 2). It is worth to note that two of amino acid residues at positions 14(V) and 51(S) were identified as specific for recent Shaanxi NDV strains in compared to the past pandemic strains in China (Table 2). Taken together, our results of the present work provide evidences suggesting that extensive sequence divergence occur among Shaanxi strains, past pandemic strains in China, and vaccine strains.

Sequence analysis of NDV isolates of HN gene of NDV
Based on the C-terminal HN extension of NDV, 4 groups were classified, including amino acid residues of 571, 577, 580, and 616 [14]. In the current study, we found that recent Shaanxi strains have HN gene encoding only 571 amino acid residues of HN glycoprotein. The percentage of nucleotide sequence identity of HN gene among recent Shaanxi strains was in the range of 99.0 %-99.7 %, while the percentage of amino acid sequence identity was in the range of 98.5 %-99.6 %. Furthermore, the percentage of nucleotide sequence identity of HN gene between Shaanxi and vaccine strains was in the range of 81.4 %-84.8 %, while the percentage of amino acid sequence identity was in the range of 87.5 %-91.2 %. Sequence analysis of Shaanxi and vaccine strains exhibited that recent Shaanxi strains of NDV had the lowest percentage of sequence identity with LaSota vaccine strain. The percentage of nucleotide and amino acid sequence identity between recent Shaanxi and LaSota vaccine strains was in the range of 81.4 %-81.6 % and 87.5 %-88.0 %, respectively.
Amino acid residue substitutions at neutralizing epitopes on the HN glycoprotein among Shannxi isolates of NDV Analysis of antigenic sites on the HN glycoprotein of recent Shaanxi NDV strains revealed that there are many amino acid residue substitutions at neutralizing epitopes on HN protein ( Table 3). As compared to vaccine strains (Herts/33, B1, Clone 30, LaSota, and Mukteswar), many amino acid residue substitutions at positions 347E → R (sites 1 and 14), 350Y → H (sites 1 and 14), 494G → D (site 1 2), 514I → V (sites 2 and 12), 518S → N (sites 2 and 12), 519S → L (sites 2 and 12) were found in recent Shaanxi strains of NDV (Table 3). Furthermore, in comparison to the past pandemic strains in China during the past 10 years, several amino acid residue substitutions at neutralizing epitopes on HN protein were uncovered in recent Shaanxi strains. Examples include amino acid residue substitutions at positions 347 (E to R) with JSD0812 whose host was duck and 347 (E to R) with sh09 whose host were chickens as well as 518S → N (sites 2 and 12) and 519S → L (sites 2 and 12) (Table 3). Additionally, compared to the recent Taiwanese NDV isolates, several amino acid residue substitutions at neutralizing epitopes on the HN protein at positions 347G → R, 348H → Q, 350Y → H, 352 T → I, 518S → N and 519S → L were also uncovered ( Table 3). It is interesting to note that 4 amino acid residues at positions 347(R), 350(H), 518 (N), and 519 (L) were conserved in recent Shaanxi NDV isolates and identified as a specific marker for these isolates (Table 3).
Previous reports showed that LaSota vaccine can be fully protected against challenge by strains from genetic groups VIb, VIg, VIId and IX [16,21]. As compared to APMV-1/ chicken/NL/152608/93 [21], we found several substitutions in F and HN proteins of shaanxi NDV isolates, such as in the amino acid sequences of F gene with             [14]. However, similar mutation was not found in the Shaanxi strain. Additionally, examination of the critical amino acid residues 401E, 416R, and 526Y for receptor binding [22] in the HN glycoprotein of recent Shaanxi strains of NDV were also conserved.

Phylogenetic analysis of F and HN genes of NDV
Based on the nucleotide sequences (1-389 bp) of the F gene among NDV isolates, a phylogenetic tree was constructed (Fig. 1). The reference strains cover traditional vaccine strains and different pandemic VII genotype NDV strains in China ( Fig. 1; Table 1). All recent Shaanxi strains were phylogenetically closely related and belonged to genetype VIId (Fig. 1). Furthermore, all recent Shaanxi strains were not phylogenetically related to vaccine strains, such as Mukteswar, Clone 30, Lasota, B1 and Herts/33. Furthermore, all recent Shaanxi strains were phylogenetically related to HA-14-07-Ch and JS/5/06/ Go NDV strains that were isolated from the Eastern China.
Based on nucleotide sequences of the full-length HN gene among NDV isolates, a phylogenetic tree was created (Fig. 2). The results in Fig. 2 showed that all recent Shaanxi strains were phylogenetically closely related and had only 571 amino acid residues. Besides, all these strains were phylogenetically related to chicken/China/ JSX1/2010, HA-14-07-Ch and JS/5/06/Go isolated from the Eastern China. However, all recent Shaanxi strains were not phylogenetically related to vaccine strains. In addition, JS/1/97/Go strain was closely related to the vaccine strain Mukteswar and Herts/33 in the phylogenetic tree based on a variable portion (nt 1-389) of the F gene, but JS/1/97/Go strain was closely related to HA-14-07-Ch in the phylogenetic tree creased on the basis of the full-length HN gene.

Relatedness (R) values of Shaanxi and LaSota strains
The neutralizing antibody (NA) titers of LaSota and NDV/Chicken/TC/1/2011 vaccination hyperimmune sera were determined by neutralization test. While NDV LaSota was used as an antigen, the NA geometric mean titer (GMT) to both LaSota and NDV/Chicken/TC/1/ 2011 sera were 128 and 148, respectively. Furthermore, the NA GMT to both LaSota and NDV/Chicken/TC/1/ 2011 vaccination sera were 67.2 and 241.8, respectively when all nine recent Shannxi strains were used as antigens.
Relatedness (R) values, determined by cross VN, revealed that all the recent Shaanxi isolates are closely related and belonged to a single serotype. However, the cross VN test revealed that R values obtained between NDV/Chicken/TC/1/2011 and LaSota vaccine strain was 0.54, indicating that they were minor subtype difference.

Hemagglutination inhibition assay
The HI assay was employed to evaluate antibody titers using LaSota antigens. As illustrated in Table 4, the HI titer was undetectable in un-inoculated and negative control group chickens. The chickens exhibited HI titer that received Oil-TC/1 were 7.1 log 2 , received Oil-Las were 7.2 log 2 , received Live-Las were 7.0 log 2 , received  Oil-TC/1and Live-Las were 8.1 log 2 , received Oil-Las and Live-Las were 8.7 log 2 . The cross-HI assay was employed to evaluate the HI antibody titers using different antigens between NDV/ Chicken/TC/1/2011 and LaSota. Anti-NDV/Chicken/ TC/1/2011 hyperimmune serum was used as an antibody, the HI titer were 9Log 2 when using the homologous antigens and7Log 2 when using the vaccine strain LaSota. When anti-NDV LaSota hyperimmune serum was used as an antibody, the HI titer were 10Log 2 that using the homologous antigens and 7Log 2 that using NDV/Chicken/TC/1/2011.

Neutralization test and cross-protectivity experiments
As shown in Table 5, neutralization titers of chicken immunized with different vaccine combinations while NDV/Chicken/TC/1/2011 strain was used as an antigen in neutralization test. The neutralization titer was undetectable in both un-inoculated and negative control group chickens. Our results reveal that chickens having neutralization titers of 64 could be protected against challenge with the NDV/Chicken/TC/1/2011 strain (Table 5).
To evaluate the protective efficacy of different vaccine combinations, the cross protectivity of each group vaccine was assessed in 3 week after the last vaccine (Table 6). Chickens in the negative control groups had no clinical signs during the course of the experiments. All birds in the un-inoculated control group displayed severe depression and typical lesions from day 3 to 5 pc and 100 % mortality was observed on day 7 pc. Necropsy of these un-inoculated control birds revealed gross lesions consistent with a virulent NDV infection including hemorrhages and edema in the conjunctiva of the lower eyelid, petechial hemorrhages in the thymus, and multifocal hemorrhage of the proventriculus. No obvious clinical signs or mortality was observed in birds receiving either Oil-TC/1 or Oil-TC/1and Live-Las vaccines, and conferred protection from morbidity and death in all chickens subjected to the NDV/Chicken/TC/1/2011 strain challenge (Table 6).
Obvious clinical signs were seen in part of other 3 vaccinated groups by 3 days post-inoculation (dpi). They displayed depressed, consumed less food and water. Two birds died on 5 dpi in Oil-Las group during the experiment, and the mortality rate was 17 %. Three birds died on 5 and 7 dpi in LaSota attenuated vaccines group during the experiment, and the mortality rate was 25 % (Table 6). One birds died on 7 dpi in the Oil-Las and Live-Las vaccines group during the experiment, the mortality rate was 9 %. However, both Oil-Las and LaSota attenuated vaccines conferred protection 83 % and 75 % of chickens, respectively. Furthermore, combination of the Oil-Las and Live-Las vaccine conferred protection to 92 % of chickens challenged with the NDV/Chicken/TC/ 1/2011 strain ( Table 6).

Frequency of isolation of challenge virus in different vaccine groups
Virus shedding was examined in all experimental chickens on day 0, 2, 4, 6, and 9 pc. As shown in Table 7, viruses could be detected in both oropharyngeal and cloacal swabs in challenge un-inoculated control group from day 2 pc. No virus was isolated from any bird in the negative control group. On day 4 pc, there were 83.3 %, 75 % and 66.7 % positive oropharyngeal swabs could be detected respectively from the Oil-Las vaccines group, Live-Las vaccines group and the vaccines group with Oil-Las and Live-Las. In this study, chickens immunized with an Oil-TC/1 vaccine have lower viruses shedding in day 2 pc from oropharyngeal samples than those immunized with a Live-Las vaccine (Fig. 3a). The virus titers of oropharyngeal swabs from both Oil-TC/1-and Oil-TC/1 plus Live-Las-inoculated chickens were significantly reduced compared to the titers of three other vaccinated groups as well as the un-inoculated control (Fig. 3b).

Discussion
There have been at least four major outbreaks of Newcastle disease worldwide so far and every epidemic had its specific genotypes [4,23]. The NDVs currently circulating worldwide exhibit multiple genotypes. Outbreaks of Newcastle disease in vaccinated chicken flocks have been reported previously in China and even in Taiwan [2,3,5,24]. Velogenic NDV strains were isolated from egg layer flocks with NDV vaccine immunization failure in China [24,25]. To date, the genotype VII of NDVs is predominant in the domestic poultry of Asia. The subgenotypes VIIc and VIId of NDVs existed in China and Far East countries [20,24,26], while sub-genotype VIIe of NDVs was found in Taiwan [5]. Earlier studies have suggested that high genetic diversity and multiple genotypes of NDVs in China were found [17,19,24,27]. Sub-genotypes VIIc andVIId of NDVs were the major (See figure on previous page.) Fig. 1 Phylogenetic tree of the nucleotide sequence of NDV strains based on a variable portion (nt 1-389) of the F gene. The accession numbers of the sequences derived in this study and those derived from GenBank were shown in Table 1. The tree was constructed using MEGA 4 neighborjoining method with 1000 replicates of bootstrap. Genotype, sub-genotype and recently Chinese isolates are indicated on the right. ▲: Shaanxi NDV isolates epidemic strains in China in recent years and could evolve into other NDV genotypes [24].
Analysis of the deduced amino acid sequences of the F protein of recent Shaanxi strains revealed that all these isolates possessed the amino acid sequences 112 R-R-Q-K-R-F 117 for the C-terminus of the F2 protein, indicating that they were all highly pathogenic viruses. The result was consistent with the pathogenicity assay. Additionally, all the recent Shaanxi strains had 101 K and 121 V, a characteristic of genotype VII viruses [15]. Furthermore, phylogenetic analysis indicated that all Shaanxi strains were classified as sub-genotype VIId. Analysis of deduced amino acid sequences of the F protein of recent Shannxi and traditional vaccine strains revealed many amino acid residue substitutions. It is worth to note that many amino acid residue substitutions on F protein of recent Shannxi NDV strains were also found as compared to the pandemic strains in China during the past 10 years. The findings suggest that these mutations on the F proteins may change antigenicity of recent Shaanxi NDV strains.
Earlier studies have demonstrated that the viral HN protein mediates attachment to sialic acid-containing receptor(s) via its neuraminidase (NA) activity, the apparently opposing activity of release of sialic acid from soluble and membrane-associated glycoconjugates [28][29][30][31]. It was also suggested that there are at least 5 antigenic sites related to neutralizing epitopes on the HN protein of NDV, such as residues 193 to 201 (site 23), residues 345 to 355 (sites 1 and 14) and a C-terminal domain composed of residues 494, 513 to 521, and 569 (sites 12 and 2). Furthermore, monoclonal antibodies (MAbs) to 3 overlapping sites (2, 12, and 23) on HN glycoprotein suppressed the NA activity of the virus [32,33]. MAbs to sites (1 and 14) inhibited HA predominantly by preventing viral attachment to chick cells [32,34]. As compared to vaccine and the past pandemic strains in China, recent Shaanxi NDV strains reveal a number of amino acid residue substitutions at neutralizing epitopes on HN protein, which suggest that these mutations may lead to antigenic change and have effect on viral attachment to the receptor on the cell surface [33,35]. This divergence further elucidates that amino acid residue substitutions may lead to the change of antibody recognition capabilities and types. Apparently, some new immune response-escaping antigenic variants were developed and resulted in ND outbreaks in Shannxi province in China. In this work, an inactivated oil mulsion (Oil-TC/1) vaccine prepared from a recent Shannxi isolate (NDV/Chicken/TC/1/2011), conferred protection from morbidity and death in all chicks subjected to homologous challenge whereas the inactivated oil mulsion Lasota vaccine (Oil-Las) and LaSota attenuated vaccine (Live-Las) conferred clinical protection to only 83 % and 75 % of chicks, respectively challenged with NDV/Chicken/TC/1/2011. It can be speculated that immunization failure may be caused by the difference of antigenic sites between vaccine strains and field strains.
In an earlier study, Liu and his colleagues also showed that SPF chickens vaccinated with LaSota vaccines could be fully protected against challenge by strains from genetic groups VIb, VIg, VIId, and IX [16]. Recently, Cornax et al. reported LaSota vaccine could provide full protection for NDV genotype VIId [21]. Conversely, several reports in using LaSota vaccine to immunize chicken in China or Taiwan [5,25] demonstrated that LaSota vaccine could not provide full protection for genotype VII of NDV. In the present study, LaSota vaccine was not  Table 1. The tree was constructed using MEGA 4 neighborjoining method with 1000 replicates of bootstrap. Genotype, sub-genotype and recently Chinese isolates are indicated on the right. ▲: Shaanxi NDV isolates capable of providing full protection for the flocks. As compared to APMV-1/chicken/NL/152608/93, XJ-2/97, and JS/1/97/Go NDV strains of genotypes VIId and VIIb [16,21], we found several substitutions in both F and HN proteins of shaanxi NDV isolates, thereby altering antigenicity of recent shaanxi NDV strains. Vaccination is the most important way to prevent NDV infections [16,24]. In the present study, the inactivated vaccine prepared from a recent Shaanxi strain NDV/Chicken/TC/1/2011 protected all the chickens from morbidity and mortality against its homologous challenge. However, the group immunized with LaSota strain inactivated and attenuated LaSota vaccines had 17 % and 25 % mortality rate, respectively. Our results are argreement with previous reports in using LaSota vaccine to immunize chicken in China or Taiwan [5,25]. Furthermore, no obvious clinical signs or mortality was observed in birds receiving either Oil-NDV/Chicken/ TC/1/2011 or Oil-NDV/Chicken/TC/1/2011 plus attenuated LaSota vaccines, and conferred protection from morbidity and death in all chickens subjected to NDV/ Chicken/TC/1/2011 strain challenge.

Conclusions
In conclusion, this study provides evidences, suggesting that recent Shaanxi strains are immune response-escaping antigenic variants that were responsible for these new NDV outbreaks in northwestern China. Therefore, a ND inactivated vaccine must be prepared from new local strains of velogenic NDV.

Virus isolation and animals
Forty Newcastle diseases suspected lungs; tracheal rings of the dead chickens (either Hy-Line Brown or Roman White chickens) from nine different poultry farms were collected from both Tongchuan city and Yaoxian county of Shaanxi province in China during 2011. The ages of these chickens range from 9 to 38-week-old. All samples were propagated in the allantoic sacs of 10-day-old SPF chicken embrocated eggs at 37°C for 3-7 days. Hemagglutination (HA) tests and hemagglutinationinhibition (HI) tests were conducted as described previously [36]. Nine NDV strains were purified by plaque purification three times. Each virus was coded as follows: virus name/isolation host/isolation place abbreviation/ case number/year isolated. For example, NDV/Chicken/ TC/1/2011 indicates that the sample is case no. 1 of a NDV strain isolated from the 2011 outbreak in Tongchuan, Shaanxi province. Intracerebral pathogenicity index (ICPI) test and mean death time (MDT) of chicken embryos test were performed using the OIE recommended method. All NDV isolates and reference isolates are shown in Table 1.
Seventy-eight of 5 days old specific pathogen-free (SPF) White Leghorn chickens were supplied by Green Square Biological Engineering Company, Yangling, china. Table 5 Neutralization titres of chicken immunized with different vaccine combinations using NDV/Chicken/TC/1/2011 strain as an antigen All birds were housed in isolators in the laboratory animal facility under negative pressure with food and water provided ad libitum. At the age of 3 weeks, the birds were inoculated with vaccines. All animal experimental procedures were approved by the Ethical and Animal Welfare Committee of Shannxi Province and were in compliance with the China law on animal experiments.

Viral RNA extraction for reverse transcription (RT)polymerase chain reaction (PCR)
Viral RNA was extracted using TIAN amp Virus DNA/ RNA kit according to manufacturer instructions. To perform RT-PCR, the NDV-RNA was extracted directly from the allantoic fluid of the NDV-inoculated embryos. After purification, the purified NDV genomic RNA was used to synthesize cDNA by reverse transcription. The primers for amplification of the F and HN genes of NDV were designed on the basis of LaSota sequences (accession no. JF950510). The primer sequences for amplification of F gene were F1 (5′-TGTAGTAACGGG AGACAAAG CAG-3′; identical to nucleotides 4669 to 4691) and F2 (5′-GAATAAATACCAG AGACATA GGGA-3′; complementary to nucleotides 5625 to 5601). The expected size of a PCR-amplified fragment is 957 bp in length. The upstream region of the HN gene of NDV isolates was amplified by PCR with the primers: HN1 (5′-AATTTCATCCCAGCGCCTAC-3′, identical to nucleotides 6893 to 6912) and HN2 (5′-ATCAT-CAAGCATCGTCCC-3′; complementary to nucleotides 7894 to 7877). Primers HN3 (5′-AGAGTCGTGCTG-GAGAAT-3′, identical to nucleotides 6431 to 6448) and HN4 (5′-CTACCTGTCCGAGTGAAA-3′; complementary to nucleotides 7768 to 7751) were used to amplify the downstream region of the HN gene. The expected upstream and downstream regions of the PCR products are 1002 bp and 1338 bp in length, respectively. These 2 regions overlap and cover the full-length HN gene. The synthesized cDNA was used as the template for PCR amplification. PCR reactions of F gene were subjected to 35 cycles consisting of denaturation for 50 s at 94°C, annealing for 50 s at 49°C, and extension for 1 min at 72°C and one final extension cycle at 72°C for 10 min. The annealing temperatures of HN1 and HN2 genes were 43°C and 46°C, respectively, and other procedures were same as described above. After completing PCR, 10 μL of the reaction mixture was loaded onto a 1.5 %

Hemagglutination inhibition assay
Blood samples were taken on 3 week after the last vaccination, and sera were isolated and the hemagglutination inhibition (HI) assay was performed using four haemagglutinating units of antigen. Sera of un-inoculated birds were collected and used as a negative control. All sera were inactivated at 56°Cfor 30 min and stored at −80°C until used. The geometric mean antibody titers were expressed on a log 2 scale. The antigen used in this experiment was LaSota vaccine strain that is used as a standard antigen in the HI test routinely used in ND diagnostic laboratories.

Cross virus-neutralization (VN) assay
The serologic relatedness of LaSota and NDV/Chicken/ TC/1/2011 NDV strains as well as between NDV/ Chicken/TC/1/2011 strain and other eight Shannxi strains were determined using a virus-neutralization (VN) test. The cross-reactivity was examined and R value was determined as described previously [37]. R value is to determine the degree of antigenic similarity of two viruses, which can be used for the similarity analysis of the protein antigen in different virus strains. Ten groups of 30 specific-pathogen-free broilers each were first vaccinated with either inactivated vaccines (nine Shannxi NDV strains) or LaSota vaccine, and then boosted by either Shannxi NDV strains or LaSota vaccine. Serum was reacted against each other in a VN test. Briefly, Two-fold serial dilutions of sera were mixed with approximately 100 tissue culture infectious doses (TCID) 50 % for 60 min at 37°C in 96-well microtiter plates. The TCID 50 of nine Shannxi NDV strains and LaSota vaccine strain were measured by the method of Reed and Muench on chicken embryo fibroblast (CEF). After the hyperimmune serum of two viruses were 2-fold serially diluted, they were mixed with equal volume containing 100 TCID 50 two viruses and then reacted at room temperature for 1 h. The culture medium was discarded and virus-serum mixture was added to trace cell reaction plates containing monolayer cells as 0.1 ml/well. Each dilution was carried out for 8 repeats. Addition of maintenance medium after 37°C for 2 h, while virus control containing 100 TCID 50 , blank cell control and serum control were included. Calculation of the 50 % serum neutralization end after culturing 72 h was performed. All neutralization tests were performed twice and geometric mean titer (GMT) were calculated. The R-value was calculated as described previously [38]. R values higher than 0.70 prove antigenic identity, R values between 0.70 and 0.33 prove antigenic relatedness meaning minor subtype differences, and R values between 0.32 and 0.11 indicate loose relatedness meaning major subtype difference, whereas R values below 0.11 indicate no relatedness at all meaning serotype difference [37].

Cross-protectivity
Both NDV/Chicken/TC/1/2011 strains and vaccine strain LaSota were used to prepare monovalent oilemulsion vaccines as described previously [5]. Live vaccine of the LaSota strain from a commercial source (Green Square Biological Engineering Company, Yangling, china) was also used. As shown in Table 8, seventy-eight SPF White Leghorn chickens were randomly divided into seven groups. At the age of 3 weeks, the chickens were inoculated. Birds received LaSota attenuated vaccine viruses that were inoculated one dose commercial live-LaSota (Live-Las) via eye-drop and intra-nasal routes. Birds received inactivated vaccines that were injected subcutaneously with 0.4 ml of inactivated oil emulsion-NDV/ Chicken/TC/1/2011 strains (Oil-TC/1) and oil emulsion-LaSota (Oil-Las), respectively, while the control group was injected with phosphate-buffered saline (PBS). After 2 weeks, a booster dose of each vaccine was administered to the birds. The cross protectivity of each group vaccine was assessed in 3 week after booster. All birds were challenged through eye-drop and intra-nasal routes with 10 4 ELD 50 of NDV strain NDV/Chicken/TC/1/2011. Challenge-free birds were administrated with PBS via the same route as the negative control. Following challenge, birds were observed for clinical signs and death during 14 day postchallenge (pc). Moribund chickens were euthanized with intravenous sodium pentobarbital at a dose of 100 mg/kg and counted dead for the next day.
Necropsies were completed on selected birds to assess the presence of gross pathological lesions. Oropharyngeal and cloacal swabs were collected at 0, 2, 4, 6 and 9 days pc for virus isolation and titration as previously described [39] and virus titers were expressed as log 10 EID 50 /ml. The oropharyngeal positive samples were then quantified the viral loads as previously reported [4].

Statistical analysis
Statistical analysis of serology titers and virus titers were performed using IBM statistical package for social sciences (SPSS) statistical software. A probability (p) value < 0.05 was considered statistically significant.