Open Access

Positive selection in the hemagglutinin-neuraminidase gene of Newcastle disease virus and its effect on vaccine efficacy

  • Min Gu1,
  • Wujie Liu1,
  • Lijun Xu1,
  • Yongzhong Cao1,
  • Chunfeng Yao1,
  • Shunlin Hu1 and
  • Xiufan Liu1Email author
Contributed equally
Virology Journal20118:150

https://doi.org/10.1186/1743-422X-8-150

Received: 16 November 2010

Accepted: 31 March 2011

Published: 31 March 2011

Abstract

Background

To investigate the relationship between the selective pressure and the sequence variation of the hemagglutinin-neuraminidase (HN) protein, we performed the positive selection analysis by estimating the ratio of non-synonymous to synonymous substitutions with 132 complete HN gene sequences of Newcastle disease viruses (NDVs) isolated in China.

Results

The PAML software applying a maximum likelihood method was used for the analysis and three sites (residues 266, 347 and 540) in the HN protein were identified as being under positive selection. Codon 347 was located exactly in a recognized antigenic determinant (residues 345-353) and codon 266 in a predicted linear B-cell epitope. Substitutions at codon 540 contributed to the N-linked glycosylation potential of residue 538. To further evaluate the effect of positively selected sites on the vaccine efficacy, we constructed two recombinant fowlpox viruses rFPV-JS6HN and rFPV-LaSHN, expressing the HN proteins from a genotype VII field isolate Go/JS6/05 (with A266, K347 and A540) and vaccine strain La Sota (with V266, E347 and T540), respectively. Two groups of SPF chickens, 18 each, were vaccinated with the two recombinant fowlpox viruses and challenged by Go/JS6/05 at 3 weeks post-immunization. The results showed that rFPV-JS6HN could elicit more effective immunity against the prevalent virus infection than rFPV-LaSHN in terms of reducing virus shedding.

Conclusions

The analysis of positively selected codons and their effect on the vaccine efficacy indicated that the selective pressure on the HN protein can induce antigenic variation, and new vaccine to control the current ND epidemics should be developed.

Background

Newcastle disease (ND) is notorious for its devastations to the world poultry industry and listed as one of the notifiable terrestrial animal diseases by the World Organization for Animal Health (Office International des Epizooties). The causative agent, Newcastle disease virus (NDV), also known as avian paramyxovirus serotype 1, is a member of the family Paramyxoviridae [1]. The virus genome is a non-segmented, single-strand, negative sense RNA which codes for six major proteins including nucleocapsid protein (NP), phosphoprotein (P), matrix protein (M), fusion protein (F), hemagglutinin-neuraminidase (HN), and large RNA-directed RNA polymerase (L), in the order from the 3' to 5' terminus [2]. Since its emergence in fowls in 1926, NDV has undergone substantial genetic evolution and has developed into several distinct genotypes (I to IX) [3, 4]. Among these, genotype VII is considered to be responsible for the severe outbreaks in Western Europe [5], South Africa and Southern Europe [3], and East Asia [6, 7] in the 1990s. Presently, the genotype VII NDV is still prevalent in China [4, 810].

Although the cleavability of F protein is pivotal to NDV pathogenicity [11, 12], recent studies have shown that HN protein also contributes to tissue tropism and virulence [13]. HN is an important immunoprotective glycoprotein on the envelope of ND virions and responsible for essential viral functions, such as binding to sialic acid-containing cell receptors, facilitating the fusion activity of the F protein and removing sialic acid to release progeny virus particles [14]. Despite the critical role that HN protein plays in NDV immunity and pathogenesis, the positive selection pressure acting on HN during the viral evolution has not been well analyzed.

The ratio of non-synonymous (d N ) to synonymous (d S ) substitutions (ω = dN/dS) provides an important means for studying the selective pressure at the protein level, with ω = 1 denoting neutral mutations, ω < 1 purifying selection, and ω > 1 diversifying positive selection. As a high proportion of amino acids in many proteins is often largely invariable (with ω close to 0) due to strong structural and functional constraints, approaches conferring an average ω over all codons across the gene are not sensitive enough to detect positive selection [15]. The program PAML [16, 17], which applies a maximum likelihood (ML) criterion and a few simple models allowing for heterogeneous ω ratios among sites, has been considered an efficient integrated method to estimate positive selection and has been commonly used to study virus evolution [1821]. In this paper, the selective pressure on NDV HN protein was examined using 132 complete HN sequences (Chinese isolates), including 106 retrieved from GenBank (up to 14 April, 2009) and the other 26 obtained from field isolates. Based on the analysis, three codons of HN were identified under positive selection and their potential effect on the routine vaccine efficacy was then evaluated.

Materials and methods

Viruses

Four pigeon isolates: NDV03-018, NDV03-044, NDV05-028 and NDV05-029 [22], were kindly provided by Dr. Zhiliang Wang (China Animal Health and Epidemiology Center). Two chicken isolates, QH-1/79 and QH-4/85 [23], were obtained from Dr. Dianjun Cao (Harbin Veterinary Research Institute, Chinese Academy of Agricultural Sciences). Twenty field strains were isolated from diseased chicken and goose flocks in China during 2005-2006. All of these viruses were subjected to three rounds of plaque-purification in chick embryo fibroblast (CEF) monolayers and subsequently propagated in 10-day-old specific pathogen free (SPF) chicken embryos. Infective allantoic fluid containing virus stocks was aliquoted and stored at -80°C before use.

RNA preparation, PCR, and sequencing

Viral RNAs were extracted directly from the allantoic fluid with the Trizol LS reagent (Invitrogen, Carlsbad, CA), following the manufacturer's instructions. Reverse transcription (RT) was conducted with random primers, and PCR was performed with a pair of primers (sense: 5'-CTTCACAACATCCGTTCTACC-3', antisense: 5'-ACCTTCCGAGTTTTATCATTCT-3') to amplify the full-length HN gene of NDV. The PCR products were purified with a DNA purification kit (QIAGEN, Hilden, Germany) and sequenced directly using the ABI PRISM BigDye Terminator v3.1 Cycle Sequencing kit (Applied Biosystems, Foster City, CA).

Sequence information and phylogenetic analysis

GenBank accession numbers assigned to the 26 strains characterized in the present study were as follows: FJ751918, FJ751919, FJ766528, EF666110, GQ338309-GQ338311, and EU044809-EU044827. In addition, the other 106 full-length HN sequences of NDV isolates from China (vaccine strains of La Sota and Mukteswar were included as they are used extensively in poultry flocks, while recombinant strains were excluded to ensure the accuracy of detecting positive selection at amino acid sites [24]) were retrieved directly from GenBank and their accession numbers were listed in Table 1. All the 132 HN sequences were edited and aligned with the Lasergene software (DNASTAR Inc., Madison, WI). The GTR (general time reversible) + I (invariable sites) + G (gamma distribution) evolutionary model was selected as the optimal nucleotide substitution model with the program Modeltest 3.7 [25]. Phylogenetic tree was then constructed by employing the ML method implemented in PAUP* version 4.0b [26] and neighbor-joining (NJ) method in MEGA version 4.0 [27]. The robustness of the statistical support for the tree branch was evaluated by 1000 bootstrap replicates. The online server, BepiPred 1.0 [28], was used to predict the position of linear B-cell epitopes of all the HN sequences.
Table 1

Accession numbers of the 106 complete HN gene sequences of NDVs isolated in China that were directly retrieved from GenBank (* sequences with lysine (K) at residue 347)

Accession Number

Accession Number

Accession Number

Accession Number

Accession Number

DQ060053

EU346660

FJ004152

DQ314571-DQ314572

DQ368684*

DQ228924

EU546165

FJ240169

DQ485229-DQ485231

DQ485270*

DQ228931

AF204872

AF456429-AF456434

DQ485262-DQ485268

DQ682452*

DQ485272

AY135171

EF211815-EF211818

DQ469832-DQ469833

DQ228925-DQ228930*

DQ486859

AY253912

FJ011444-FJ011447

DQ234586-DQ234587

DQ228932-DQ228935*

DQ659677

AY351959

FJ386392-FJ386396

DQ682447-DQ682451

DQ234584-DQ234585*

DQ858355

AY997298

DQ023148-DQ023155

DQ023156*

DQ234590-DQ234592*

EF201805

EF540730

DQ023554-DQ023557

DQ023558*

DQ469830-DQ469831*

EU649675

EF175144

DQ023559-DQ023560

DQ234579*

FJ217668-FJ217669*

EU481973

EF141104

DQ234580-DQ234583

DQ234588*

 

Positive selection detection

To estimate the selective constraints on the HN protein, the codeml program of the PAML package (version 4) [16] was utilized to calculate the site-to-site variation in ω. Two nested site-specific models, consisting of a neutral model that does not allow positive selection (ω≤1) and an alternative model that permits positive selection (ω > 1), were compared. As recommended [15], the following models were used: M0 (one-ratio) v. M3 (discrete) and M7 (beta) v. M8 (beta & ω). M0 assumes a constant ω for all codons whereas M3 allows for discrete classes of sites with different ω ratios. M7 supposes a beta distribution with 10 categories of ω over sites, each corresponding to a unique ω value that is always less than 1 while M8 has an extra category with ω > 1. Then the log likelihood values for each pair of the above nested models were compared by a likelihood ratio test (LRT) [15, 17], in order to assess whether the model allowing for positive selection is significantly more suitable for the data. Finally, the Bayes empirical Bayes (BEB) procedure [29] was used to infer the particular codons under positive selection and to calculate their posterior probabilities.

Animal experiment

To further investigate the effect of positively selected sites on vaccine efficacy against the prevalent NDVs, Go/JS6/05 (field NDV strain) and La Sota (the most widely used vaccine strain in China) were chosen to construct corresponding recombinant fowlpox viruses (rFPVs) expressing each HN gene based on the transfer vector pP12LS developed by Sun et al [30]. The expression was identified by indirect immunofluorescence assay (IFA) in secondary CEF cultures using anti-NDV polyclonal antibody as previously described [31], and the levels of HN expression were further compared between the two generated rFPVs by flow cytometry [32] on DF-1 cells (a stable cell line of CEF) at a multiplicity of infection (MOI) of 5. Subsequently, two groups of five-day-old SPF White Leghorn chickens (18 birds/group, Beijing Merial Vital Laboratory Animal Technology, Beijing, China) were immunized respectively with the above two rFPVs at a dose of 1 × 104 PFU. A third group was served as a mock-vaccinated control. Three weeks later, all chickens were challenged oculonasally with 100 μL of PBS-diluted allantoic fluid containing 1 × 105 EID50 of Go/JS6/05. Tracheal and cloacal swabs were collected on days 3, 5 and 7 post-challenge (p.c.). Furthermore, six chickens from each vaccinated group were sacrificed humanely on day 5 p.c., and tissue samples including liver, brain, spleen, kidney, trachea and lung were collected. The swabs were immersed in PBS with antibiotics (8000 U/mL ampicilin, 5 mg/mL streptomycin and kanamycin, pH 7.2), and stored at -80°C until analyzed. The recovery of the challenged virus in these swabs or organ samples was confirmed by inoculation into embryonated chicken eggs. All animal work was approved by the Jiangsu Administrative Committee for Laboratory Animals (Permission number: SYXK-SU-2007-0005).

Results

Sequence analysis

The full coding region of the 132 HN genes analyzed in this study exhibited diverse phylogenetic phenotypes in Chinese poultry flocks, containing six of the nine recognized genotypes (II, III, VI, VII, VIII and IX), with the overwhelming majority (113/132) belonged to genotype VII (Additional file 1, Table S1). In addition, the phylogenetic tree (Figure 1) showed that our 26 sequences belonged to four different genotypes as follows: III (Ch/JS7/05 and Go/JS9/05), VI (NDV05-028 and NDV05-029), VIII (QH-1/79 and QH-4/85) and the remaining 20 fell into genotype VII.
Figure 1

Phylogenetic tree of the 26 NDV field isolates based on the full coding region of the HN gene. Reference sequences representing recognized genotypes were denoted with their taxon names underlined. The tree was generated using the neighbor-joining method implemented in program MEGA 4.0. The Class I NDV strain DE-R49/99 was served as an outgroup. The scale indicates the branch length based on the number of nucleotide substitutions per site and the numbers at nodes indicate bootstrap values calculated by resampling of 1000 replicates.

Detection of positive selection

An ML method implemented in the software package PAML was used for the identification of any positive selection on the HN protein of NDV. The log likelihood differences between M0 and M3, as well as M7 and M8 were found to be significant. Models that permitted positive selection showed a better fit to the data and contained a class of codons with the non-synonymous to synonymous substitutions ratio greater than one (ω > 1), indicating the existence of positive selection (Table 2). A further analysis on the amino acids most likely responsible for the detected non-neutral pattern revealed that three codons, 266, 347 and 540, were under positive selection identified by both M3 and M8, with the posterior probability over 95% for residue 266 and 347, and over 90% for residue 540. As the M0-M3 comparison is more a test of heterogeneity in the ω value among sites and not actually a test of positive selection [33], only the results obtained by M8 were investigated further.
Table 2

Parameter estimates and likelihood values corresponding to the results from different models

Model

lnL

dN/dS

Estimated parameters

Dfa

2lb

Positively selected codonsc

M0 (one-ratio)

-13384.287

0. 225

ω = 0. 225

4

 

None

M3 (discrete)

-13246. 011

0. 228

p0= 0.656, p1= 0.310, p2= 0.034

ω0= 0.069, ω1= 0.446, ω2 = 1.313 d

 

276.552 (p < 0.01)

266, 347, 540

M7 (beta)

-13253. 972

0. 227

p = 0.486, q = 1.635

2

 

Not allowed

M8 (beta& ω)

-13247. 102

0. 228

p0= 0.968, p1= 0.032

p = 0.662, q = 2.730, ω = 1.312

 

13.740 (p < 0.01)

266, 347, 540

a Df: degree of freedom.

b 2l: Log likelihood difference between nested models using the χ2-test.

c Codons estimated to be under positive selection at the 90% level. Numbers in bold denote codons with posterior probabilities over 95%.

d Underlined parameters are evidence of positive selection.

Amino acid variations of positively selected codons

As shown in Table 3, each of the three positively selected codons identified by M8 exhibited diversity in amino acid substitutions which could induce the variations in hydrophobicity or charge. In particular, codon 266 contained the most abundant amino acid substitutions. According to program BepiPred 1.0, a fragment of the HN sequence covering residues 255 to 265, sometimes extending to residue 267 or even 270, was detected as a potential B-cell epitope. And codon 266 was involved in the above predicted epitope of 119 HN sequences, including all the 113 genotype VII viruses (Additional file 1, Table S1).
Table 3

Amino acid variations of the positively selected codons identified by M8

Codon

Residue

Substitution frequencya

Hydrophobicityb

Charge

site ω( ± S.E.)c

266

A

91

Y

Neutral

1. 477 ± 0. 125

 

T

23

N

Neutral

 
 

V

10

Y

Neutral

 
 

I

4

Y

Neutral

 
 

P

2

Y

Neutral

 
 

D

1

N

Negative

 
 

S

1

N

Neutral

 

347

E

62

N

Negative

1. 485 ± 0. 100

 

G

41

Y

Neutral

 
 

K

28

N

Positive

 
 

N

1

N

Neutral

 

540

A

94

Y

Neutral

1. 438 ± 0. 202

 

V

23

Y

Neutral

 
 

T

15

N

Neutral

 

a Numbers represent the occurrence times of each specific residue among the 132 HN sequences.

b Y stands for hydrophobic and N for hydrophilic.

c Estimated ω value at single site with standard error.

Generation of rFPVs expressing HN genes

Two rFPVs, rFPV-LaSHN and rFPV-JS6HN respectively expressing the HN proteins of vaccine strain La Sota (with V266, E347 and T540) and genotype VII field isolate JS-06/05 (with A266, K347 and A540) were generated by homologous recombination between the corresponding transfer plasmid and the wild-type parental FPV. Fluorescence was readily observed in CEF cells transfected with either rFPV-LaSHN (Figure 2A) or rFPV-JS6HN (Figure 2B), demonstrating the expression of the inserted HN gene. The flow cytometric analysis indicated that there was no significant difference in the expression levels of HN between rFPV-LaSHN and rFPV-JS6HN on DF-1 cells at 24 hours post-infection (data not shown).
Figure 2

Identification of HN protein expressed in the recombinant fowlpox virus (rFPV) by indirect immunofluoresence assay (IFA). Chick embryo fibroblast (CEF) cells were infected with rFPV-LaSHN (A), rFPV-JS6HN (B), or wild-type parental FPV (C) at a multiplicity of infection (MOI) of 0.01. Sixty hours post-infection, IFA was performed using polyclonal anti-NDV chicken serum and FITC-conjugated rabbit anti-chicken IgG as primary and secondary antibody, respectively. The fluorescence was observed under an inverted fluorescence microscope.

Protective efficacies of rFPV-JS6HN and rFPV-LaSHN

On day 21 after immunization, effective antibody responses were induced in both rFPV vaccinated groups, with higher hemagglutinin-inhibition (HI) titers to the homologous antigen, and serum from birds inoculated with rFPV-LaSHN reacted poorly with Go/JS6/05 (Table 4). On day 5 p.c., virus replication in different visceral organs of rFPV-JS6HN or rFPV-LaSHN vaccinated chickens was examined by inoculation into 10-day-old embryonated chicken eggs. As shown in Table 4, the frequencies of virus isolation in the brain, spleen and lung from rFPV-LaSHN group were higher than that of the rFPV-JS6HN, though no statistically significant difference was observed.
Table 4

HI titers to NDV at day pre-challenge and virus detection in visceral organs from different rFPV immunized groups on day 5 p.c.

Vaccine group

HI titer ± standard deviationa

Virus detection (positive/total)c

 

Go/JS6/05

La Sota

Liver

Brain

Spleen

Kidney

Trachea

Lung

rFPV-JS6HN

3.80 ± 0.92 b

3.30 ± 1.16

0/6

3/6

3/6

4/6

0/6

2/6

rFPV-LaSHN

1.50 ± 0.85

3.90 ± 0.57

0/6

5/6

6/6

4/6

0/6

5/6

a HI titers to Go/JS6/05 or La Sota were evaluated on day 21 post vaccination.

b titers to homologous antigens were noted in bold.

c Virus replication was detected on day 5 post challenge by inoculation into 10-day-old SPF embryonated chicken eggs.

Each of the rFPV-JS6HN and rFPV-LaSHN vaccinated chickens was fully protected against mortality after challenge, whereas unvaccinated birds died within 5 days p.c.. On day 3 p.c., virus shedding from the cloaca and trachea showed that both the rFPV vaccines remarkably decreased the level of virus excretion from cloaca, and that the rFPV-JS6HN group could significantly reduce the virus recovery rates from trachea when compared with the rFPV-LaSHN group. Moreover, on days 5 and 7 p.c., the frequencies of virus isolation from cloaca in rFPV-JS6HN vaccinated birds were much lower than that in rFPV-LaSHN immunized fowls (Table 5).
Table 5

Virus recovery of swab samples from different rFPV immunized groups challenged with NDV Go/JS6/05

Vaccine group

No. of birds shedding/total no. of birds on day after challenge:

Mortality

 

3

5

7

 
 

Trachea

Cloaca

Trachea

Cloaca

Trachea

Cloaca

 

rFPV-JS6HN

11/18B

4/18B

0/18A

1/18A

0/12A

0/12A

0/18A

rFPV-LaSHN

16/18AB

7/18B

1/18A

9/18B

0/12A

2/12A

0/18A

Unvaccinated

17/18A

14/18A

NS

NS

NS

NS

18/18B

A, B Different superscript letters denote significant difference between treatment groups within the column; Fisher's exact test, P < 0.05.

NS: No survivors.

Discussion

Positive selection is an evolutionary process that could drive the fixation of emerging advantageous mutations in the population with higher frequencies compared to the wild-type allele [34]. Therefore, identifying proteins or protein domains that experiencing adaptive selection will improve the understanding of their genomic functions and the recognition of genetic variation that leads to phenotypic diversity [35].

Observations from previous genetic and antigenic studies of viruses such as FMDV (foot-and-mouth disease virus) [18], HIV-1 (human immunodeficiency virus type 1) [19], RHDV (rabbit hemorrhagic disease virus) [20] and influenza B virus [21], have indicated that signatures of positive selection are generally functionally important and/or associated with antigenicity. To date, seven antigenic determinants that form a continuum on HN protein have been characterized by a panel of monoclonal antibodies (mAbs) against the HN of the velogenic Australia-Victoria/32 (AV) strain, including the amino acids positions 193, 194, 201, 263, 287, 321, 332, 333, 345, 347, 350, 353, 356, 494, 513, 514, 516, 521 and 569 [36, 37]. One of our positively selected sites, codon 347, was located exactly in those defined epitopes. In the present study, both La Sota and Mukteswar, which are widely used vaccine strains in China, have a glutamic acid (E) occupied codon 347, in contrast to that with a lysine (K) substitution resulting in the opposite residual charge exclusively in genotype VII viruses (accession numbers in italic in Table 1). Furthermore, as referred to recent work of Cho et al [38] and Hu et al (2009) [39], it is reasonable to postulate that the emergence of E347K substitution might be closely related to the host immune pressure. However, the codon 266 was not included in the aforementioned antigenic sites [36, 37], instead, it was involved in a predicted linear B-cell epitope that simultaneously held the epitope residue 263, suggesting that site 266 may lie in some antigenic regions yet to be recognized.

N-linked glycosylation, one of the most common forms of protein post-translational modifications, is known to be correlated with viral infectivity and immune escape [40]. There are six potential N-glycosylation sites (amino acids 119, 341, 433, 481, 508 and 538) in the HN protein of the AV strain [41]. In our analysis, positive selection was detected at codon 540, which comprised three different amino acids: alanine (A), valine (V) and threonine (T). Residue 538 was conserved with asparagine (N) in all the 132 HN sequences and tended to be a putative N-glycosylation site if T was present at site 540. However, the vast majority of prevalent strains owned A or V at codon 540 (Additional file 1, Table S1), which would deprive the possibility of N538 being glycosylated. The exact function of the resulted deglycosylation at site 540 remains unknown and needs to be further explored.

Compared to vaccine strain La Sota, most genotype VII NDV isolates possessed different amino acids at the three identified positively selected sites. To further evaluate the effect of those sites on the vaccine protective efficacy, Go/JS6/05 was chosen together with La Sota for the recombinant fowlpox-virus construction. Before challenge, serum collected from the rFPV-LaSHN immunized chickens displayed lower HI titers to Go/JS6/05 than to La Sota (Table 4), which may suggest that substitutions at the positively selected sites are partially responsible for the antigenic variation between the two HN proteins. After challenge, virus shedding results showed that the rFPV-JS6HN could prevent the excretion of the challenged virus more efficiently than rFPV-LaSHN, indicating that the positively selected sites on the HN protein could affect the vaccine immune efficacy against the prevalent NDV infection.

Although an intensive vaccination program against ND has been executed in China in the last few decades, epidemic infections with velogenic genotype VII NDV in vaccinated birds are still frequently reported in recent years [4, 810]. Currently, the most extensively used vaccine strains, such as La Sota (genotype II) and Mukteswar (genotype III), were isolated and characterized in the 1940s and belonged to the "early" genotypes (I-IV and IX), which have evident amino acid sequence divergence from the "late" ones (V-VII), especially genotype VII [42]. The results in this study suggest that positive selection may play a role in the formation of such differentiation and even induce antigenic variations compared with the vaccine strains. Therefore, new vaccine to better control the ND epizootics of prevalent NDV strains carrying novel variations at identified positively selected sites should be developed to meet the challenge.

Notes

Declarations

Acknowledgements

This work was supported by National Natural Science Foundation of China (grant numbers: 30630048 and 30800814) and China Layer Industry Improvement program (grant number: nycytx-41-g07).

Authors’ Affiliations

(1)
Animal Infectious Disease Laboratory, School of Veterinary Medicine, Yangzhou University

References

  1. Alexander DJ: Newcastle disease and other avian paramyxoviruses. Rev Sci Tech 2000, 19: 443-462.PubMedGoogle Scholar
  2. Kattenbelt JA, Stevens MP, Gould AR: Sequence variation in the Newcastle disease virus genome. Virus Res 2006, 116: 168-184. 10.1016/j.virusres.2005.10.001View ArticlePubMedGoogle Scholar
  3. Herczeg J, Wehmann E, Bragg RR, Travassos Dias PM, Hadjiev G, Werner O, Lomniczi B: Two novel genetic groups (VIIb and VIII) responsible for recent Newcastle disease outbreaks in Southern Africa, one (VIIb) of which reached Southern Europe. Arch Virol 1999, 144: 2087-2099. 10.1007/s007050050624View ArticlePubMedGoogle Scholar
  4. Liu XF, Wan HQ, Ni XX, Wu YT, Liu WB: Pathotypical and genotypical characterization of strains of Newcastle disease virus isolated from outbreaks in chicken and goose flocks in some regions of China during 1985-2001. Arch Virol 2003, 148: 1387-1403.PubMedGoogle Scholar
  5. Lomniczi B, Wehmann E, Herczeg J, Ballagi-Pordany A, Kaleta EF, Werner O, Meulemans G, Jorgensen PH, Mante AP, Gielkens AL, et al.: Newcastle disease outbreaks in recent years in western Europe were caused by an old (VI) and a novel genotype (VII). Arch Virol 1998, 143: 49-64. 10.1007/s007050050267View ArticlePubMedGoogle Scholar
  6. Yu L, Wang Z, Jiang Y, Chang L, Kwang J: Characterization of newly emerging Newcastle disease virus isolates from the People's Republic of China and Taiwan. J Clin Microbiol 2001, 39: 3512-3519. 10.1128/JCM.39.10.3512-3519.2001PubMed CentralView ArticlePubMedGoogle Scholar
  7. Mase M, Imai K, Sanada Y, Sanada N, Yuasa N, Imada T, Tsukamoto K, Yamaguchi S: Phylogenetic analysis of Newcastle disease virus genotypes isolated in Japan. J Clin Microbiol 2002, 40: 3826-3830. 10.1128/JCM.40.10.3826-3830.2002PubMed CentralView ArticlePubMedGoogle Scholar
  8. Qin ZM, Tan LT, Xu HY, Ma BC, Wang YL, Yuan XY, Liu WJ: Pathotypical characterization and molecular epidemiology of Newcastle disease virus isolates from different hosts in China from 1996 to 2005. J Clin Microbiol 2008, 46: 601-611. 10.1128/JCM.01356-07PubMed CentralView ArticlePubMedGoogle Scholar
  9. Liu H, Wang Z, Wu Y, Wu Y, Sun C, Zheng D, Xu T, Li J: Molecular characterization and phylogenetic analysis of new Newcastle disease virus isolates from the mainland of China. Res Vet Sci 2008, 85: 612-616. 10.1016/j.rvsc.2008.02.013View ArticlePubMedGoogle Scholar
  10. Liu H, Wang Z, Wu Y, Zheng D, Sun C, Bi D, Zuo Y, Xu T: Molecular epidemiological analysis of Newcastle disease virus isolated in China in 2005. J Virol Methods 2007, 140: 206-211. 10.1016/j.jviromet.2006.10.012View ArticlePubMedGoogle Scholar
  11. Panda A, Huang Z, Elankumaran S, Rockemann DD, Samal SK: Role of fusion protein cleavage site in the virulence of Newcastle disease virus. Microb Pathog 2004, 36: 1-10. 10.1016/j.micpath.2003.07.003View ArticlePubMedGoogle Scholar
  12. Peeters BP, de Leeuw OS, Koch G, Gielkens AL: Rescue of Newcastle disease virus from cloned cDNA: evidence that cleavability of the fusion protein is a major determinant for virulence. J Virol 1999, 73: 5001-5009.PubMed CentralPubMedGoogle Scholar
  13. Huang Z, Panda A, Elankumaran S, Govindarajan D, Rockemann DD, Samal SK: The hemagglutinin-neuraminidase protein of Newcastle disease virus determines tropism and virulence. J Virol 2004, 78: 4176-4184. 10.1128/JVI.78.8.4176-4184.2004PubMed CentralView ArticlePubMedGoogle Scholar
  14. Scheid A, Choppin PW: Isolation and purification of the envelope proteins of Newcastle disease virus. J Virol 1973, 11: 263-271.PubMed CentralPubMedGoogle Scholar
  15. Yang Z, Nielsen R, Goldman N, Pedersen AM: Codon-substitution models for heterogeneous selection pressure at amino acid sites. Genetics 2000, 155: 431-449.PubMed CentralPubMedGoogle Scholar
  16. Yang Z: PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol 2007, 24: 1586-1591. 10.1093/molbev/msm088View ArticlePubMedGoogle Scholar
  17. Nielsen R, Yang Z: Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene. Genetics 1998, 148: 929-936.PubMed CentralPubMedGoogle Scholar
  18. Haydon DT, Bastos AD, Knowles NJ, Samuel AR: Evidence for positive selection in foot-and-mouth disease virus capsid genes from field isolates. Genetics 2001, 157: 7-15.PubMed CentralPubMedGoogle Scholar
  19. Travers SA, O'Connell MJ, McCormack GP, McInerney JO: Evidence for heterogeneous selective pressures in the evolution of the env gene in different human immunodeficiency virus type 1 subtypes. J Virol 2005, 79: 1836-1841. 10.1128/JVI.79.3.1836-1841.2005PubMed CentralView ArticlePubMedGoogle Scholar
  20. Esteves PJ, Abrantes J, Carneiro M, Muller A, Thompson G, van der Loo W: Detection of positive selection in the major capsid protein VP60 of the rabbit haemorrhagic disease virus (RHDV). Virus Res 2008, 137: 253-256. 10.1016/j.virusres.2008.07.025View ArticlePubMedGoogle Scholar
  21. Shen J, Kirk BD, Ma J, Wang Q: Diversifying selective pressure on influenza B virus hemagglutinin. J Med Virol 2009, 81: 114-124. 10.1002/jmv.21335PubMed CentralView ArticlePubMedGoogle Scholar
  22. Liu H, Wang Z, Son C, Wang Y, Yu B, Zheng D, Sun C, Wu Y: Characterization of pigeon-origin Newcastle disease virus isolated in China. Avian Dis 2006, 50: 636-640. 10.1637/7618-042606R1.1View ArticlePubMedGoogle Scholar
  23. Liang R, Cao DJ, Li JQ, Chen J, Guo X, Zhuang FF, Duan MX: Newcastle disease outbreaks in western China were caused by the genotypes VIIa and VIII. Vet Microbiol 2002, 87: 193-203. 10.1016/S0378-1135(02)00050-0View ArticlePubMedGoogle Scholar
  24. Anisimova M, Nielsen R, Yang Z: Effect of recombination on the accuracy of the likelihood method for detecting positive selection at amino acid sites. Genetics 2003, 164: 1229-1236.PubMed CentralPubMedGoogle Scholar
  25. Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817View ArticlePubMedGoogle Scholar
  26. Swofford DL: PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods). Sunderland, Massachusetts: Sinauer Associates; 2002.Google Scholar
  27. Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol 2007, 24: 1596-1599. 10.1093/molbev/msm092View ArticlePubMedGoogle Scholar
  28. Larsen JE, Lund O, Nielsen M: Improved method for predicting linear B-cell epitopes. Immunome Res 2006, 2: 2. 10.1186/1745-7580-2-2PubMed CentralView ArticlePubMedGoogle Scholar
  29. Yang Z, Wong WS, Nielsen R: Bayes empirical bayes inference of amino acid sites under positive selection. Mol Biol Evol 2005, 22: 1107-1118. 10.1093/molbev/msi097View ArticlePubMedGoogle Scholar
  30. Sun L, Liu WJ, Chen SJ, Shi HY, Liu XF: Influence of nonessential region on protective efficacy of recombinant fowlpox viruses. Acta Microbiologica Sinica 2005, 45: 359-362. (In Chinese)PubMedGoogle Scholar
  31. Liu WJ, Sun L, Chen SJ, Xu ZL, Liu JB, Liu XF: Construction and protective efficacies of two recombinant fowlpox viruses expressing the HA gene of H9 subtype AIV. Acta Veterinaria et Zootechnica Sinica 2008, 39: 327-332. (In Chinese)Google Scholar
  32. Ferreira L, Munoz-Barroso I, Marcos F, Shnyrov VL, Villar E: Sialidase, receptor-binding and fusion-promotion activities of Newcastle disease virus haemagglutinin-neuraminidase glycoprotein: a mutational and kinetic study. J Gen Virol 2004, 85: 1981-1988. 10.1099/vir.0.79877-0View ArticlePubMedGoogle Scholar
  33. Wong WS, Yang Z, Goldman N, Nielsen R: Accuracy and power of statistical methods for detecting adaptive evolution in protein coding sequences and for identifying positively selected sites. Genetics 2004, 168: 1041-1051. 10.1534/genetics.104.031153PubMed CentralView ArticlePubMedGoogle Scholar
  34. Suzuki Y, Gojobori T: A method for detecting positive selection at single amino acid sites. Mol Biol Evol 1999, 16: 1315-1328.View ArticlePubMedGoogle Scholar
  35. Biswas S, Akey JM: Genomic insights into positive selection. Trends Genet 2006, 22: 437-446. 10.1016/j.tig.2006.06.005View ArticlePubMedGoogle Scholar
  36. Iorio RM, Bratt MA: Monoclonal antibodies to newcastle disease virus: delineation of four epitopes on the HN glycoprotein. J Virol 1983, 48: 440-450.PubMed CentralPubMedGoogle Scholar
  37. Iorio RM, Borgman JB, Glickman RL, Bratt MA: Genetic variation within a neutralizing domain on the haemagglutinin-neuraminidase glycoprotein of Newcastle disease virus. J Gen Virol 1986,67(Pt 7):1393-1403. 10.1099/0022-1317-67-7-1393View ArticlePubMedGoogle Scholar
  38. Cho SH, Kwon HJ, Kim TE, Kim JH, Yoo HS, Kim SJ: Variation of a newcastle disease virus hemagglutinin-neuraminidase linear epitope. J Clin Microbiol 2008, 46: 1541-1544. 10.1128/JCM.00187-08PubMed CentralView ArticlePubMedGoogle Scholar
  39. Hu S, Wang T, Liu Y, Meng C, Wang X, Wu Y, Liu X: Identification of a variable epitope on the Newcastle disease virus hemagglutinin-neuraminidase protein. Vet Microbiol 2010, 140: 92-97. 10.1016/j.vetmic.2009.07.029View ArticlePubMedGoogle Scholar
  40. Vigerust DJ, Shepherd VL: Virus glycosylation: role in virulence and immune interactions. Trends Microbiol 2007, 15: 211-218. 10.1016/j.tim.2007.03.003View ArticlePubMedGoogle Scholar
  41. McGinnes LW, Wilde A, Morrison TG: Nucleotide sequence of the gene encoding the Newcastle disease virus hemagglutinin-neuraminidase protein and comparisons of paramyxovirus hemagglutinin-neuraminidase protein sequences. Virus Res 1987, 7: 187-202. 10.1016/0168-1702(87)90027-XView ArticlePubMedGoogle Scholar
  42. Miller PJ, Decanini EL, Afonso CL: Newcastle disease: evolution of genotypes and the related diagnostic challenges. Infect Genet Evol 2010, 10: 26-35. 10.1016/j.meegid.2009.09.012View ArticlePubMedGoogle Scholar

Copyright

© Gu et al; licensee BioMed Central Ltd. 2011

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Advertisement