Open Access

Molecular characterization of partial fusion gene and C-terminus extension length of haemagglutinin-neuraminidase gene of recently isolated Newcastle disease virus isolates in Malaysia

Virology Journal20107:183

https://doi.org/10.1186/1743-422X-7-183

Received: 29 December 2009

Accepted: 8 August 2010

Published: 8 August 2010

Abstract

Background

Newcastle disease (ND), caused by Newcastle disease virus (NDV), is a highly contagious disease of birds and has been one of the major causes of economic losses in the poultry industry. Despite routine vaccination programs, sporadic cases have occasionally occurred in the country and remain a constant threat to commercial poultry. Hence, the present study was aimed to characterize NDV isolates obtained from clinical cases in various locations of Malaysia between 2004 and 2007 based on sequence and phylogenetic analysis of partial F gene and C-terminus extension length of HN gene.

Results

The coding region of eleven NDV isolates fusion (F) gene and carboxyl terminal region of haemagglutinin-neuraminidase (HN) gene including extensions were amplified by reverse transcriptase PCR and directly sequenced. All the isolates have shown to have non-synonymous to synonymous base substitution rate ranging between 0.081 - 0.264 demonstrating presence of negative selection. Analysis based on F gene showed the characterized isolates possess three different types of protease cleavage site motifs; namely 112RRQKRF117, 112RRRKRF117 and 112GRQGRL117 and appear to show maximum identities with isolates in the region such as cockatoo/14698/90 (Indonesia), Ch/2000 (China), local isolate AF2240 indicating the high similarity of isolates circulating in the South East Asian countries. Meanwhile, one of the isolates resembles commonly used lentogenic vaccine strains. On further characterization of the HN gene, Malaysian isolates had C-terminus extensions of 0, 6 and 11 amino acids. Analysis of the phylogenetic tree revealed that the existence of three genetic groups; namely, genotype II, VII and VIII.

Conclusions

The study concluded that the occurrence of three types of NDV genotypes and presence of varied carboxyl terminus extension lengths among Malaysian isolates incriminated for sporadic cases.

Background

Newcastle disease (ND) is a highly contagious disease of birds and has been regarded throughout the world as one of the most important diseases of poultry and other birds [1], in which infection with the extremely virulent viruses may result in sudden, high mortality with comparatively few clinical signs. The causative agent, NDV, is avian Paramyxovirus under the Avulavirus and has a negative-sense, single-stranded RNA genome [2]. So far, NDV strains with genomic sizes of 15,186, 15192 and 15198 nucleotides which codes for at least six proteins including nucleoprotein (N), phosphoprotein (P), matrix (M) protein, fusion (F) protein, haemagglutinin-neuraminidase (HN) protein and RNA polymerase (L) [24] have been identified. Among the six major proteins, the two interactive surface glycoproteins, the F and the HN proteins, are involved in cell surface attachment and cell membrane fusion [3, 5].

The molecular basis for NDV pathogenicity has been shown to be dependent on the F protein cleavage site amino acid sequence which is clearly described by OIE [6] molecular definition of virulent NDV stating that any AMV-1 virus that has three basic amino acids, either lysine (K) or arginine (R), at the fusion protein cleavage site between residues 113 and 116 at the C-terminus of the F2, as well as phenylalanine at residue 117 of F1 and carboxyl terminus amino acid extension length which varies due to the varying location of termination codons within the HN protein, resulting in the expression of HN proteins with varying amino acid lengths. A longer reading frame consisting of HN0 precursor of 616 amino acid residue is expressed only by avirulent NDV strains and biologically active HN proteins of 571 and 577 amino acid residues are expressed by virulent and lentogenic viruses, respectively [7].

Three different NDV genotypes, II, III, and IV, were involved in the first panzootic of ND and were restricted to the specific geographic region; South East Asia in which the outbreak began. In the late 1960 s, NDV genotypes V and VI emerged and caused the second and third panzootics, respectively. After that, two novel NDV genotypes, VII and VIII, were found in Asia, Southern Africa, and a number of European countries [812]. Genotype VII was mainly responsible for recent outbreaks in the neighbouring countries of Taiwan and China [8, 1315] constituting the fourth panzootic of NDV.

Intensive vaccine programs have been implemented in Malaysia, but ND outbreaks and sporadic cases have occasionally occurred, even in well-vaccinated farms. A major epidemic of ND has occurred in Peninsular Malaysia from 2000-2001 peaking with 84 outbreaks and 525981 cases in 2001 [16] which cause substantial losses. Isolates of low virulence, HitchnerB1 and LaSota are the most common type of vaccines being used in the world including Malaysia. Other vaccines used in Malaysia include S, Ulster 2C, NDV-6/10 and enteric vaccine strain VG-GA [17]. Despite intensive vaccination programs with live vaccines, NDV remains a constant threat to the commercial poultry. Here, we describe molecular characterization of F and C-terminus extension length of HN protein genes of recently isolated Malaysian isolates and their phylogenetic relationship of among NDV isolates derived from other countries or regions. Thus, characterization of these recent isolates may help to gain invaluable information about the pathogenicity and epidemiological relationships.

Results

Nucleotide and predicted amino acid sequence analysis

Sequence analysis of a region between nucleotide positions 47-435 of F gene, encompassing proteolytic cleavage site of F0 protein revealed 18-55 nucleotide and 4-21 amino acid substitutions in Malaysian isolates and the calculation of synonymous and non-synonymous substitution rate demonstrate that all isolates have a rate ranging between 0.081 - 0.264, which is a value less than one. As illustrated in Table 1, in all the NDV isolates examined, the rate of synonymous substitution (0.126 - 0.306) dominated the rate of non-synonymous substitution (0.014 - 0.032). The maximum amino acid substitution was observed at amino acid position 11 compared to consensus in which, 6 isolates (MB043/06, MB091/05, MB093/05, MB095/05, MB128/04 and MB85/05) had threonine (T) residue, and 3 isolates (MB016/07, MB064/05, MB061/06) had an alanine (A) residue. Isolates MB091/05, MB093/05, MB095/05 and MB128/04 shared unique amino acid substitutions with MB043/06 at amino acid position 13; proline (P) for leucine (L), C25F, T29A, K78R and Q114R (Table 2 and Table 3). The unique amino acid substitution suggests that these isolates might have a common origin despite being isolated in different areas in Malaysia. As shown in Table 4, the deduced amino acid sequence of the F0 protein cleavage site revealed that isolates MB043/06, MB091/05, MB093/05, MB095/05, and MB128/04 had an arginine (R) residue at position 114, which result in an 112RRRKRF117 motif and the remaining isolates MB047/05, MB064/05, MB076/05 and MB085/05) carry the residue 112RRQKRF117. Only one isolate (MB016/07) possessed the motif of 112GRQGRL117 common to avirulent isolates. Comparisons of nucleotide identities with published local and foreign isolates indicated that seven isolates (MB076/05, MB047/05, MB043/06, MB091/05, MB093/05, MB095/05 and MB128/04) were found to have maximum identity (91.2% to 96.4%) with Indonesian isolate (Cockatoo/Indonesia/14698/90), two isolates (MB064/05 and MB016/07) with strain Ch/2000 from China (94.8 to 96.6%) and the remaining isolates, MB085/05 and MB061/06 with the local Malaysian isolate AF2240 (97.4%) and vaccine strain LaSota (98.8%), respectively (data not shown). These findings indicated the high similarity of isolates circulating in the South East Asian countries but still having limited variation with isolates from different geographical areas.
Table 1

Number of nucleotide and amino acid substitutions and Ka/Ks ratio of in the current Malaysian NDV isolates.

NDV Isolates

No. of

nucleotide

substitution*

No. of

amino acid

substitution*

Non-synonymous

substitution rate (KA)

Synonymous

substitution

rate (KS)

Ratio

(KA/KS)

MB061/06

55

21

0.081

0.306

0.265

MB043/06

37

8

0.030

0.276

0.109

MB016/07

29

9

0.028

0.202

0.139

MB064/05

22

5

0.014

0.173

0.081

MB076/05

24

8

0.032

0.144

0.222

MB047/05

18

4

0.018

0.126

0.143

MB091/05

34

9

0.032

0.242

0.132

MB093/05

34

9

0.032

0.242

0.132

MB095/05

35

9

0.032

0.252

0.127

MB128/04

34

9

0.032

0.242

0.132

MB085/05

18

4

0.014

0.134

0.104

* compared to the consensus nucleotide and amino acid sequences.

Table 2

Amino acid residue substitution of F gene for NDV strains of different genotypes.

Isolate/genotype/subgenotype

11

12

13

14

16

17

18

19

20

22

25

26

28

29

30

36

Consensus

V

P

L

M

I

T

R

I

M

I

C

I

L

T

S

P

MB061/06II

A

.

M

.

T

I

.

V

A

V

.

.

P

A

N

.

MB043/06VII

T

.

P

.

.

.

.

.

.

.

F

.

.

A

.

.

MB016/07 VII

A

.

.

.

.

I

.

.

.

.

Y

V

.

.

.

.

MB064/05VII

A

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

MB076/05VII

.

.

.

I

T

.

.

.

.

.

.

.

.

.

.

S

MB047/05VII

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

MB091/05 VII

T

S

P

.

.

.

Q

.

.

.

F

.

.

A

.

.

MB093/05 VII

T

S

P

.

.

.

Q

.

.

.

F

.

.

A

.

.

MB095/05 VII

T

S

P

.

.

.

Q

.

.

.

F

.

.

A

.

.

MB128/04 VII

T

S

P

.

.

.

Q

.

.

.

F

.

.

A

.

.

MB085/05 VIII

T

.

.

.

.

.

.

.

T

.

.

.

.

.

.

.

AF2240 VIII

T

.

.

.

.

.

.

.

T

.

.

.

P

.

.

.

QH-1/79 VIII

T

.

S

.

.

.

.

.

.

.

.

.

.

.

G

.

QH-4/85VIII

T

.

S

.

.

.

.

.

.

.

.

.

.

.

G

.

V4 Queensland I

.

.

.

.

T

V

.

V

.

A

.

V

P

.

.

.

Ulster/67 I

.

.

.

.

T

V

.

V

A

E

.

V

P

.

.

.

LaSota II

A

.

M

T

T

I

.

V

A

V

.

.

P

A

N

.

Miyadera III

A

.

.

.

T

I

W

.

A

A

.

V

.

.

.

.

Herts/33 IV

A

.

P

.

.

I

.

.

V

T

.

.

.

.

.

.

Italien IV

.

.

.

.

.

I

.

.

A

T

.

.

.

.

.

.

CA1085/71 V

.

.

.

.

.

.

.

.

T

.

.

.

.

.

.

.

H-10/72 V

.

.

.

.

.

.

.

.

T

.

.

.

.

.

.

.

Iraq AG68 VI

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

Lebanon -70 VI

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

TX3503/04 VI

A

.

.

.

.

.

.

.

T

.

.

.

.

.

.

.

NDV05-027VI

A

.

.

.

.

.

.

.

T

V

.

.

.

.

.

.

Q-GB 506/97VI

P

.

.

.

.

.

.

.

V

.

S

V

.

.

.

.

DK-1/95VI

.

.

P

.

.

.

.

.

T

.

.

.

.

 

.

.

DE 143/95VII

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

Cockatoo/14698/90VII

.

.

.

.

.

.

.

.

.

.

.

.

S

.

.

.

ZA360/95VII

.

.

.

.

.

.

.

V

.

.

.

.

.

.

G

.

ZW 3422/95VII

.

.

.

.

.

.

.

V

.

.

.

.

.

.

G

.

NDV05-055 VII

A

.

.

I

.

.

.

.

.

 

.

.

.

A

.

.

Ch/2000VII

A

.

.

.

.

.

.

.

.

.

.

.

P

.

.

.

TW/2000 VII

A

.

.

T

.

.

.

.

.

.

.

.

.

.

.

.

TW/95-1VII

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

.

ZhJ-1/85

A

.

.

.

T

A

.

.

A

A

.

V

.

.

N

.

FJ-1/85

A

.

.

.

T

V

.

.

A

A

.

V

.

.

N

.

Genetic identity with consensus sequence is marked by dot (.).The bolded residues were residues discussed in the text.

Table 3

Amino acid residue substitution of F gene for NDV strains of different genotypes.

Isolate/genotype/Subgenotype

52

71

78

79

82

101

104

107

108

114

115

117

121

124

Consensus

I

K

K

A

E

R

G

S

T

Q

K

F

I

S

MB061/06 II

.

.

.

.

D

.

E

T

.

.

G

L

.

G

MB043/06 VII

.

.

R

.

.

K

.

.

M

R

.

.

.

.

MB016/07 VII

V

R

.

.

.

K

.

.

.

.

.

.

V

.

MB064/05 VII

V

R

.

.

.

K

.

.

.

.

.

.

V

.

MB076/05 VII

.

.

.

.

.

.

.

.

.

.

.

.

V

.

MB047/05 VII

.

.

.

.

.

K

.

A

.

.

.

.

V

.

MB091/05 VII

.

.

R

.

.

K

.

.

.

R

.

.

.

.

MB093/05 VII

.

.

R

.

.

K

.

.

.

R

.

.

.

.

MB095/05 VII

.

.

R

.

.

K

.

.

.

R

.

.

.

.

MB128/04 VII

.

.

R

.

.

K

.

.

.

R

.

.

.

.

MB085/05 VIII

.

.

.

T

.

.

.

T

.

.

.

.

.

.

AF2240 VIII

.

.

.

P

.

.

.

T

.

.

.

.

V

.

QH-1/79 VIII

.

.

R

T

.

.

.

T

.

.

.

.

V

.

QH-4/85 VIII

.

.

R

T

.

.

.

T

.

.

.

.

V

.

V4 Queensland I

.

.

.

.

.

.

E

T

.

.

G

L

.

G

Ulster/67 I

.

.

.

.

.

.

E

T

.

.

G

L

.

G

LaSota II

.

.

.

.

D

.

E

T

.

.

G

L

.

G

Miyadera III

.

.

.

.

.

.

E

T

.

.

R

.

.

.

Herts/33 IV

.

.

.

.

.

.

E

T

.

.

R

.

.

.

Italien IV

.

.

.

.

.

.

E

T

.

.

R

   

CA1085/71 V

.

.

.

.

.

.

.

T

.

.

.

.

.

.

H-10/72 V

.

.

.

.

.

.

.

T

.

.

.

.

.

.

Iraq AG68 VIa

.

.

.

.

.

.

.

.

.

.

.

.

.

.

Lebanon -70 VIa

.

.

.

.

.

.

.

.

.

.

.

.

.

.

TX3503/04 VIb

.

.

.

.

.

.

.

.

.

K

.

.

.

.

NDV05-027VIb

.

.

.

.

.

.

.

.

.

.

.

.

.

.

Q-GB 506/97VIc

.

R

.

.

 

.

.

.

.

.

.

.

.

.

DK-1/95VId

.

.

.

.

.

.

.

.

.

.

.

.

.

.

DE 143/95VIIa

.

.

.

.

.

K

.

.

.

.

.

.

V

.

Cockatoo/14698/90VIIa

.

.

.

.

.

K

.

.

.

.

.

.

V

.

ZA360/95VIIb

.

.

.

.

.

.

.

.

.

.

.

.

V

.

ZW 3422/95VIIb

.

.

.

.

.

.

.

.

.

.

.

.

V

.

NDV05-055 VIIc

.

.

.

.

.

K

.

.

.

.

.

.

V

.

Ch/2000VIId

V

R

.

.

.

K

.

.

.

.

.

.

V

.

TW/2000 VIId

V

.

.

.

.

K

.

.

.

.

.

.

V

.

TW/95-1VIIe

.

.

.

T

.

K

.

.

.

.

.

.

V

.

ZhJ-1/85 IX

.

.

.

.

.

.

E

T

.

.

.

.

.

.

FJ-1/85 IX

.

.

.

.

.

.

E

T

.

.

.

.

.

.

Table 4

List of Newcastle disease virus isolates collected in Malaysia between 2004 and 2007 and their deduced amino acid sequence at the cleavage site.

Isolate ɸ

Date of collection (DD/MM/YY)

F protein cleavage site*

Accession No. (F gene)

Origin

MB061/06

15/11/06

112GRQGR↓L117

GQ901891

Selangor

MB043/06

17/10/06

112RRRKR↓F117

GQ901896

Selangor

MB016/07

14/06/07

112RRQKR↓F117

GQ901894

Melaka

MB064/05

24/8/05

112RRQKR↓F117

GQ901893

Selangor

MB076/05

7/10/05

112RRQKR↓F117

GQ901892

Sabah

MB047/05

5/07/05

112RRQKR↓F117

GQ901895

Selangor

MB091/05

22/11/05

112RRRKR↓F117

GQ901897

Sabah

MB093/05

6/12/05

112RRRKR↓F117

GQ901898

Sabah

MB095/05

6/12/05

112RRRKR↓F117

GQ901899

Sabah

MB128/04

29/12/04

112RRRKR↓F117

GQ901900

Selangor

MB085/05

22/11/05

112RRQKR↓F117

GQ901901

Sabah

ɸ All the isolates were isolated from chicken except for MB061/06 isolated from parrot.

* Numbers indicate amino acid positions.

Carboxyl terminus of HN protein gene

Nucleotide sequencing and subsequent deduction of the amino acid sequence covering the C- terminus of the HN protein revealed differences in HN length of amino acid sequences. Isolates MB016/07, MB043/06, MB047/05, MB064/05, MB076/05, MB091/05, MB093/05, MB095/05, and MB128/04 had no amino acid extension length with a total HN length of 571 amino acids regardless of their cleavage site sequence profiles, whereas isolates MB085/05 and MB061/06 had HN C-terminus extension length of 11 aa and 6 aa, respectively, as shown in Table 5. Isolate MB085/05 had the same C-terminus aa extension length and composition as that of the local isolate AF2240 while isolate MB061/06 revealed similar aa and C-terminus extension length with that of common vaccine strains such as LaSota, VG/GA, and Herts/33 used in Malaysia.
Table 5

Deduced HN protein amino acid sequence and C-terminus extension lengths of Malaysian NDV isolates.

NDV isolate

Deduced amino acid sequence

C-terminus amino acid extension length**

HN***

Accession No.

MB061/06

KDDGVREARSG*

6

577

GQ922495

MB043/06

KDNRA*

0

571

GQ922500

MB016/07

KDDRV*

0

571

GQ922498

MB064/05

KDDRV*

0

571

GQ922497

MB076/05

KDDRV*

0

571

GQ922496

MB047/05

KDDRV*

0

571

GQ922499

MB091/05

KDDRV*

0

571

GQ922501

MB093/05

KDDRV*

0

571

GQ922502

MB095/05

KDDRV*

0

571

GQ922503

MB128/04

KDDRV*

0

571

GQ922504

MB085/05

KDDRLQEVRSGRLSQP*

11

581

GQ922505

*Indicates the stop codon

**No C-terminus amino acid extension length was detected in all the isolates except for MB061/06 and MB085/05 with 6 and 11 amino acid extension lengths, respectively.

***Predicted length of the HN protein in number of amino acids based on ORF analysis of the gene nucleotide sequence.

Phylogenetic analysis

Phylogenetic tree generated based on the variable portion of the F gene between nt 47 - 435 of the 3 isolates and other 35 local and foreign NDV isolates worldwide separated into 9 potential clusters corresponding to the different genotypes of NDV. All the Malaysian isolates formed six major clusters despite their geographical entity and within these, the presence of viruses belonging to three of the nine genotypes were identified (Figure 1). Isolates MB016/07and MB043/06 characterized in this study were assigned to genotype VII within they belonged to previously established subgenotype VIId (Figure 1). Isolate MB043/06 together with other five local strains segregated into a monophyletic group with 100% bootstrap value within genotype VII as shown in Figure 1. Considering the topology of the phylogenetic tree, these Malaysian isolates might have a common origin
Figure 1

Phylogenetic relationship among 38 NDV isolates on basis of F gene nucleotide sequences between position 47 and 435. Sequences were obtained either from the present study or GenBank. The phylogenetic tree was constructed using neighbour- joining method on MEGA 4. Malaysian isolates in the current study were indicated by the rotated black square.

Previous phylogenetic analysis performed based on variable region of F gene at nt 206-421(216 bp) classified isolate MB085/05 as genotype VII [18]. In this study, however, it was grouped together with previously characterized local NDV isolate AF2240 (isolated in the 1960s) and Chinese isolates (QH-1/79 and QH-4/85) in genotype VIII. The phylogenetic grouping of MB085/05 as genotype VIII is supported by unique V11→T, A79→T and S107→T substitutions although isolates grouped into genotype VII in previous studies lacked the aforementioned characteristic (Table 2 and Table 3). The phylogenetic tree also showed that the lentogenic isolate (MB061/06) was closely related to strain LaSota.

Discussion

Analysis of the ratio of synonymous and non-synonymous substitution rate in the current isolates demonstrated the presence of only purifying (negative) selection, despite the most variable portion of the F gene between nucleotide positions 47-435 was used in the analysis. Similar results were also observed previously on analysis of other entire NDV genes such as L (Ka/Ks = 0.064), M (Ka/Ks = 0.14 and P (Ka/Ks = 0.25) genes [1921] which corresponds to negative selection.

The amino acid sequence of the protease cleavage site reveals that all of the isolates except isolate MB061/06 maintained multiple basic amino acids motifs within the penta-amino acid sequence of the F0 cleavage signal between residues 113 and 116 and phenylalanine (F) on the residue 117 on the N-terminus. It is widely accepted that the number of basic amino acids immediately upstream to the F0 protein cleavage site determines viral pathogenicity which is clearly described by OIE [6]. The presence of these characteristic patterns of amino acid demonstrated that the isolates could be considered as virulent except isolate MB061/06. It is a paramount importance to note that, the F0 cleavage site of isolate MB043/06 isolated from Selangor, Peninsular Malaysia, is unusual, containing an arginine (R) for glutamine (Q) substitution at residue 114. In our previous studies, isolates MB091/05, MB093/05, MB095/05, and MB128/04 from different districts of Sabah, West Malaysia showed similar F0 cleavage site characteristics [18]. However, in recent years, emergence of similar results were reported, in South African genotype VIII viruses [11], in Taiwan [15] and from Eurasian collard dove and pigeons isolates containing a 112R-R-K-K-R116 , 112R-R-Q-K-R116, and 112R-R-R-K-R116 motif [2226]. Even though the contribution of arginine (R) at amino acid 114 in our isolates needs further studies, other studies indicated that arginine residue at different positions 113, 115 and 116 contribute for intracellular cleavage of virulent NDV fusion proteins [27].

Analysis of C-terminus extension length of HN protein gene revealed that the nine virulent NDV isolates shared 0 amino acid extension length with a total HN length of 571 amino acids regardless of their cleavage site sequence profiles (terminating in the sequence KDDRV and KDNRA), except MB085/05. In the present study, isolate MB085/05 was shown to have 11aa HN C-terminus extension length. This might have been due to a point mutations occurred at position 8164 (T→C) replacing the stop codon 8164TAA8166 with 8164CAA8166 codon and hence delayed terminations. This result may suggest a genetic relationship of MB085/05 sequences to the local virulent isolates with a genetic evolution that has led to the presence of these varied HN C- terminus extension lengths. In contrast, MB061/06 isolated from parrot was shown to have 6 amino acid C-terminus extension lengths, similar to that of common vaccine strains such as LaSota, VG/GA, and Herts/33 used in Malaysia. However, it is shorter than other avirulent strains such as Queensland V4 strain and Ulster/67, most of which have 45 amino acid extensions. Substitutions at positions 8149 (C→G) and 8150 (G→A) might have contributed for generation of the stop codon 8149TAG8151, terminating early and hence fewer amino acid extensions. It has been indicated that the shortest amino acid value (571 amino acids) was found in genotypes III-VIII consisting exclusively of viscerotropic velogenic strains [4]. The phylogenetic relationships based on the partial sequence of F protein gene of these isolates also showed that all the nine NDV isolates with shortest amino acid value (571) appeared to be closely related to viruses from genotype VII. According to Gould et al [28] and Kattenbelt et al [29], Australian isolates have HN extensions of 7, 9, 11, 14 or 45 amino acids and they indicated that viruses with 7 or 14 amino acid extensions were shown to be associated with summer respiratory disease.

The phylogenetic relationship result showed that the 3 Malaysian isolates recovered in 2006 and 2007; belong to genotype II and VII. Closer phylogenetic relationship of genotype VII viruses suggests that isolate MB043/06 and local viruses isolated from 2004-2005 belong to a subgenotype VIId with 100% bootstrap value, all of the members of which encodes the 112RRRKRF117 motif at the fusion cleavage site (Figure 1). In addition, we note that all those isolates share 6 common unique substitutions that distinguished them from all other subgenotypes/genotypes. According to Yu et al [30], genotype VIIb viruses evolved from genotype VIb viruses by producing a VII specific V121 for I substitution and then changed to VIIa and VIIc by producing a K101-for-R substitution and became VIId by producing additional V52-for-I and Y314-for-F substitutions. Exceptions were observed in the substitution of V at amino acid position 121 in our isolate but it did share the K substitution at amino acid residue 101 in which both K101 and V121 were unique features of genotype VII. Lien et al. [31] reports that Taiwanese subgenotype VIId viruses collected between 2003 and 2006 contained K71 and V52 in which similar features were observed in isolate MB016/07. Thus it is predicted that isolate MB016/07 can also be grouped under subgenotype VIId, which represents the most newly emerging NDV strains frequently occurring in the region especially from China and Taiwan. These findings supported by topology of phylogenetic tree suggest that it genetically related to genotype VII viruses but with the distinguishable subcluster.

Our study has also indicated the occurrence of other sub-genotype/genotype; VIIb and VIII viruses. Malaysian subgenotype VIIb virus in this study would have been more closely related to local and South East Asian countries' isolates because of their regional proximity. However, it has shared high sequence identities (91.7 - 92.5%) with Southern Africa subgenotype VIIb viruses than other Malaysian genotype VII viruses (87.4 - 90.4%). Previous studies also proposed that subgenotype VIIb strains isolated in Southern Africa might have originated in the Far East as their lack of diversity indicated that they cannot be indigenous to South Africa [8]. The finding in this work revealed that Malaysian genotype VII isolates are also phylogenetically related and share a common ancestor (supported by a high bootstrap value of 98) with Southern African isolates strengthening the hypothesis that genotype VII viruses were endemic in South East Asian countries. It has been documented that genotype VIII viruses were isolated from ND outbreaks in western China between 1979 and 1985, Southern Africa and are endemic in Southern Africa [8, 11, 13]. The presence of genotype VIII virus neighbouring country Singapore in early 1960 s and in Malaysia (AF2240) that gave rise to the recent genotype VIII virus (MB085/05) suggests that it was maintained in the region as endemic infection.

In Malaysia, NDV has been isolated frequently from chickens and most of the isolates were belong to genotype VII viruses [18, 32]. The MB061/06 isolate characterized in this study is the first genotype II NDV isolated from parrot in the peninsular Malaysia since it was reported in chicken in early 1980 s [33]. It is worth to point out that the grouping of isolate MB061/06 and LaSota under the same cluster in genotype II and the existence of similar C-terminus amino acid extensions length of HN protein suggested MB061/06 might be generated in nature from the LaSota vaccine strain.

Conclusions

Our study indicates the occurrence of genotypes; II, VII and VIII and the presence of varied C-terminus extension lengths within Malaysian isolates. The deduced amino acid sequence of the F0 protein cleavage site showed also a unique amino acid motif in one of the isolates incriminated for sporadic cases occurred in different areas of the country. Approaches towards full characterization of isolates with unique F0 cleavage signal amino acid sequence and 6 amino acid extensions should be continued and intensified for deeper molecular knowledge and better intervention strategies.

Methods

Isolates

A total of eleven (11) NDV isolates, which were previously detected positive by real time RT-PCR between 2004 and 2007 at Biologics Laboratory, Faculty of Veterinary Medicine, Universiti Putra Malaysia (UPM) were employed in this study (Table 1). The virus isolates were propagated in the allantoic cavity of 9 days old specific-pathogen-free (SPF) embryonated chicken eggs according to European Community directive 92/66/EC [34] and identified by haemaglutination tests. Infected allantoic fluid samples were clarified by centrifugation and supernatant was stored at -70°C for later analysis.

Viral RNA extraction

RNA extraction was performed from infected allantoic fluid using TRIzol® Reagent (Invitrogen, USA) in accordance with the manufacturer's instructions. The RNA pellet was re-suspended in 20 μl nuclease-free water (Promega, USA) after air dried for immediate use or kept in -80°C for later use.

Primers and reverse transcriptase-polymerase chain reaction (RT-PCR)

PCR amplification and sequencing were performed using previously described degenerative primers 5'-ATGGGC(C/T)CCAGA(C/T)CTTCTAC-3' (sense) and 5'-CTGCCACTGCTAGTTGTGATAATCC-3' (antisense) specific to fusion (F) protein gene [13] and HNNDV314 5'-ATATCCCGCAGTCGCATAAC-3'(sense) and HNNDV304 5'-TTTTTCTTAATCAAGTGACT-3'(antisense) specific to HN protein gene [35]. The primers generate an expected amplicon size of 535 bp (nt 47-535) fragment spanning the regions between nucleotides 47 and 581 of the fusion protein, which includes the F0 cleavage site and 320 bp products representing a fragment within HN protein gene, respectively. Standard RT-PCR was performed using Access® One- step RT-PCR kit (Promega, USA) in a reaction volume of 25 μl. The cycling parameters for F gene specific primers were 48°C for 45 min at reverse transcription (RT), and 35 cycles of 94°C for 2 min, 56°C for 2 min, and 72°C for 1 min, followed by 72°C for 10 min and for HN gene specific primers; 48 min at 45°C for RT, 95 °C for 5 min, followed by 35 cycles of 1 min at 95°C, 1 min at 51°C and 1 min at 72°C and a final extension step at 72°C for 10 min (MyCycler® Thermal Cycler, Bio-Rad, Hercules, CA, USA). The amplicons were separated by 1.5% agarose gel electrophoresis and visualized under ultraviolet after being stained with ethidium bromide.

Sequencing of PCR products

Specific bands for each gene of interest were excised and purified by using GENEALL™Gel SV kit (General Biosystem, Inc., Korea) following the manufacturer's instructions. The purified DNA encoded for F and HN regions were sequenced by direct sequencing in both directions. Sequencing reactions were performed using ABI PRISM® BigDye™ Terminator Cycle Sequencing Ready reaction kit v2.0 (Perkin Elmer, USA) in an automated DNA Sequencer (ABI PRISM® 377 DNA Sequencer). Each sample was sequenced three times to confirm consistency of the results.

Nucleotide and deduced amino acid sequence analysis

Percent nucleotide identity and sequence editing were carried out using BioEdit software package version 7.01 [36]. Nucleotide analysis, prediction of amino acid sequences, and alignments were performed by Molecular Evolutionary Genetics Analysis, version 4.01 (MEGA 4) [37]. Determination of synonymous and non-synonymous substitution rates was conducted using the Nei-Gojobori method in Mega 4 [37]. Nucleotide sequences of partial F and HN gene of studied isolates were deposited in the GenBank and accession numbers are shown in Tables 4 and 5, respectively.

Phylogenetic analysis

Phylogenetic analysis was carried out by comparing the variable portion of the F gene between 47 to 435 (389 bp) nucleotide sequences using the Clustal W multiple alignment method using MEGA 4 [37]. The phylogenetic tree was constructed using the neighbour-joining method after 1,000 bootstrap replicates. The sequence of the representative strains from each genotype and different geographical areas were retrieved from GenBank. Since no full-lengh sequence information at nt 47 to 435 was available for the remaining eight Malaysian isolates which were phylogenetically analysed in our previous study [18], sequences in this study were used. The accession numbers and country are as follows: AF2240 [AF048763], Malaysia; cockatoo/14698/90 [AY562985], Indonesia; Italien [EU293914], Italy; V4 Queensland [AF217084], Australia; Ulster/67 [AY562991], North Ireland; LaSota [DQ195265], USA; Miyadera [M18456], Japan; Herts/33 [AY741404], UK; CA1085/71 [AF001106], USA; H-10/72 [AF001107], Hungary; TX3503/04 [EU477190], USA; IraqAG68 [AF001108], Iraq; Ch/2000 [AF358788], China; TW/2000 [AF358786], Taiwan; TW/95-1 [AF083960], Taiwan; NDV05-055 [DQ439910], China; DK-1/95 [AF001129], Denmark; NDV05-027 [DQ439884], China; Q-GB 506/97 [AF109887], UK; QH-1/79 [AF378250], China; QH-4/85 [AF378252], China; ZhJ-1/85 [AF458023], China; FJ-1/85 [AF458009], China; ZA360/95 [AF109876], South Africa; ZW 3422/95 [AF109877], Zimbabwe; DE143/95 [AF109881], UK; Lebanon-70 [AF001110], Lebanon.

Declarations

Acknowledgements

This research was funded by grant number 01-02-04-009 BTK/ER/38 from the Ministry of Science, Technology and Innovation, Malaysia. Ayalew Berhanu was supported by a grant from the Netherlands Organization for International Cooperation in Higher Education (Nuffic).

Authors’ Affiliations

(1)
Faculty of Veterinary Medicine, Universiti Putra Malaysia
(2)
Faculty of Veterinary Medicine, University of Gondar
(3)
Institute of Bioscience, Universiti Putra Malaysia

References

  1. Aldous EW, Alexander DJ: Detection and differentiation of Newcastle disease virus (Avian paramyxovirus type 1). Avian Pathol. 2001, 30: 117-128. 10.1080/03079450120044515.PubMedView ArticleGoogle Scholar
  2. Mayo MA: A summary of taxonomic changes recently approved by ICTV. Arch Virol. 2002, 147: 1655-1656. 10.1007/s007050200039.PubMedView ArticleGoogle Scholar
  3. de Leeuw O, Peeters B: Complete nucleotide sequence of Newcastle disease virus: evidence for the existence of a new genus within the subfamily Paramyxovirinae. J Gen Virol. 1999, 80: 131-136.PubMedView ArticleGoogle Scholar
  4. Czeglédi A, Ujvari D, Somogyia E, Wehmanna E, Werner O, Lomniczi B: Third genome size category of avian paramyxovirus serotype 1 (Newcastle disease virus) and evolutionary implications. Virus Res. 2006, 120: 36-48. 10.1016/j.virusres.2005.11.009.PubMedView ArticleGoogle Scholar
  5. Lamb RA, Kolakofsky D: Paramyxoviridae: The viruses and their replication. Fields Virology. Edited by: Fields BN, Knipe DM, Howley PM. 1996, Philadelphia: Lippincott-Raven, 1177-1203. 3Google Scholar
  6. Office International des Epizooties (OIE): Manual of Diagnostic Tests and Vaccines for Terrestrial Animals. Paris. 2009, Chapter 2.3.14: 576-589. Accessed 10 September 2009, [http://www.oie.int/fr/normes/mmanual/A_summry.htm]Google Scholar
  7. Sakaguchi T, Toyoda T, Gotoh B, Inocencio NM, Kuma K, Miyata T, Negai Y: Newcastle disease virus evolution I. multiple lineages defined by sequence variability of the haemagglutinin-neuraminidase gene. Virol. 1989, 169: 260-272. 10.1016/0042-6822(89)90151-7.View ArticleGoogle Scholar
  8. 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/s007050050624.PubMedView ArticleGoogle Scholar
  9. Ke MG, Liu JH, Lin YMH, Chen J, Tsai SS, Chang CP: Molecular characterization of Newcastle disease viruses isolated from recent outbreaks in Taiwan. J Virol Methods. 2001, 97: 1-11. 10.1016/S0166-0934(01)00296-8.PubMedView ArticleGoogle Scholar
  10. Lomniczi B, Wehmann E, Herczeg J, Ballagi-Pordany A, Kaleta EF, Werner O, Meulemans G, Jorgensen PH, Mante AP, Gielkens AL, Capua I, Damoser J: 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/s007050050267.PubMedView ArticleGoogle Scholar
  11. Abolnik C, Horner RF, Bisschop SPR, Parker ME, Romito M, Viljoen GJ: A phylogenetic study of South African Newcastle disease virus strains isolated between 1990 and 2002 suggests epidemiological origins in the Far East. Arch Virol. 2004, 149: 603-619. 10.1007/s00705-003-0218-2.PubMedView ArticleGoogle Scholar
  12. 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.012.PubMedView ArticleGoogle Scholar
  13. 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-0.PubMedView ArticleGoogle Scholar
  14. Yang CY, Shieh HK, Lin YL, Chang PC: Newcastle disease virus isolated from recent outbreaks in Taiwan phylogenetically related to viruses (genotype VII) from recent outbreaks in Western Europe. Avian Dis. 1999, 43: 125-130. 10.2307/1592771.PubMedView ArticleGoogle Scholar
  15. Tsai H, Chang KH, Tseng CH, Frost KM, Manvell RJ, Alexander DJ: Antigenic and genotypical characterization of Newcastle disease viruses isolated in Taiwan between 1969 and 1996. Vet Microbiol. 2004, 104: 19-30. 10.1016/j.vetmic.2004.09.005.PubMedView ArticleGoogle Scholar
  16. OIE Multiannual Animal Disease Status (2004). Accessed 15 October 2007, [http://www.oie.int/hs2/report.asp]
  17. Aini I: Newcastle disease. Diseases of Poultry in Southeast Asia. Edited by: Zamri-Saad M. 2006, Serdang, Malaysia: Malindo Printers Sdn. Bhd, 225-233.Google Scholar
  18. Tan SW, Ideris A, Omar AR, Yesuff K, Hair-Bejo M: Sequence and phylogenetic analysis of Newcastle disease virus genotypes isolated in Malaysia between 2004 and 2005. Arch Virol. 2010, 155: 63-70. 10.1007/s00705-009-0540-4.PubMedView ArticleGoogle Scholar
  19. Seal BS, King DJ, Meinersmann RJ: Molecular evolution of the Newcastle disease virus matrix protein gene and phylogenetic relationships among the Paramyxoviridae. Virus Res. 2000, 66: 1-11. 10.1016/S0168-1702(99)00119-7.PubMedView ArticleGoogle Scholar
  20. Wise MG, Sellers HS, Alvarez R, Seal BS: RNA-dependent RNA polymerase gene analysis of worldwide Newcastle disease virus isolates representing different virulence types and their phylogenetic relationship with other members of the Paramyxoviridae. Virus Res. 2004, 104: 71-80. 10.1016/j.virusres.2004.01.034.PubMedView ArticleGoogle Scholar
  21. Locke DP, Sellers HS, Crawford JM, Schultz-Cherry S, King DJ, Meinersmann RJ, Seal BS: Newcastle disease virus phosphoprotein gene analysis and transcriptional editing in avian cells. Virus Res. 2000, 69: 55-68. 10.1016/S0168-1702(00)00175-1.PubMedView ArticleGoogle Scholar
  22. Terregino C, Cattoli G, Grossele B, Bertoli E, Tisato E, Capua I: Characterization of Newcastle disease virus isolates obtained from Eurasian collared doves (Streptopelia decaocto) in Italy. Avian Pathol. 2003, 32: 63-68. 10.1080/0307 945021000070732.PubMedView ArticleGoogle Scholar
  23. Ujvari D, Wehmann E, Kaleta EF, Werner O, Savic V, Nagy E, Czifra G, Lomniczi B: Phylogenetic analysis reveals extensive evolution of avian paramyxovirus type 1 strains of pigeons (Columba livia) and suggests multiple species transmission. Virus Res. 2003, 96: 63-73. 10.1016/S0168-1702(03)00173-4.PubMedView ArticleGoogle Scholar
  24. 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.2002.PubMedPubMed CentralView ArticleGoogle Scholar
  25. Mebatsion T, Koolen MJM, de Vaan LTC, de Haas N, Braber M, Römer-Oberdörfer A, van den Elzen P, van der Marel P: Newcastle disease virus (NDV) marker vaccine: an immunodominant epitope on the nucleoprotein gene of NDV can be deleted or replaced by a foreign epitope. J Virol. 2002, 76: 10138-10146. 10.1128/JVI.76.20.10138-10146.2002.PubMedPubMed CentralView ArticleGoogle Scholar
  26. Huovilainen A, Ek-Kommonen C, Manvell R, Kinnunen L: Phylogenetic analysis of avian paramyxovirus 1 strains isolated in Finland. Arch Virol. 2001, 146: 1775-1785. 10.1007/s007050170063.PubMedView ArticleGoogle Scholar
  27. Fuji Y, Sakaguchi T, Kiyotani K, Yoshida T: Comparison of substrate specificities against the fusion glycoprotein of virulent Newcastle disease virus between a chick embryo fibroblast processing protease and mammalian subtilisin-like proteases. Microbiol and Immunol. 1999, 43: 133-140.View ArticleGoogle Scholar
  28. Gould AR, Hansson E, Selleck K, Kattenbelt JA, MacKenzie M, Della-Porta AJ: Newcastle disease virus fusion and haemagglutinin-neuraminidase gene motifs as markers for viral lineage. Avian Pathol. 2003, 32: 361-373. 10.1080/0307945031000121112.PubMedView ArticleGoogle Scholar
  29. 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.001.PubMedView ArticleGoogle Scholar
  30. Yu L, Wang Z, Jiang Y, Chang L, Kwang J: Characterization of newly emerging Newcastle disease virus isolates from the People's Republic China and Taiwan. J Clin Microbiol. 2001, 39: 3512-3519. 10.1128/JCM.39.10.3512-3519.2001.PubMedPubMed CentralView ArticleGoogle Scholar
  31. Lien YY, Lee JW, Su HY, Tsai HJ, Tsai MC, Hsieh CY, Tsai SS: Phylogenetic characterization of Newcastle disease viruses isolated in Taiwan during 2003-2006. Vet Microbiol. 2007, 123: 194-202. 10.1016/j.vetmic.2007.03.006.PubMedView ArticleGoogle Scholar
  32. Maizan M, Sharifah SH, Ong GH, Fuad R, Suriani MN: Genetic and phylogenetic analysis of Newcastle disease virus isolated from recent outbreaks. Proceedings of the 2nd International Congress/13th VAM Congress and CVA-Australasia/Oceania Regional Symposium, Kuala Lumpur, Malaysia. 67-69.Google Scholar
  33. Aldous EW, Mynn JK, Banks J, Alexander DJ: A molecular epidemiological study of avian paramyxovirus type I (Newcastle disease virus) isolates by phylogenetic analysis of a partial nucleotide sequence of fusion protein gene. Avian Pathol. 32: 239-257. 10.1080/030794503100009783.Google Scholar
  34. CEC: Council Directive 92/66/EC of 14 July 1992 introducing community measures for the control of Newcastle disease. Offic J Eur Comm. 1992, L260: 1-20.Google Scholar
  35. Peroulis-Kourtis I, O'Riley K, Grix D, Condron RJ, Ainsworth C: Molecular characterisation of Victorian Newcastle disease virus isolates from 1976 to 1999. Aust Vet J. 2002, 80: 422-424. 10.1111/j.1751-0813.2002.tb11004.x.PubMedView ArticleGoogle Scholar
  36. Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser. 1999, 41: 95-98.Google Scholar
  37. Tamura K, Dudley J, Nei M, Kumar S: MEGA 4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol and Evolu. 2007, 24: 1596-1599. 10.1093/molbev/msm092.View ArticleGoogle Scholar

Copyright

© Berhanu et al; licensee BioMed Central Ltd. 2010

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