- Open Access
Molecular epidemiology of Newcastle disease virus isolates from vaccinated commercial poultry farms in non-epidemic areas of Japan
Virology Journal volume 10, Article number: 330 (2013)
Newcastle Disease (ND) is a highly contagious and economically devastating disease of poultry. At present, limited molecular epidemiological data are available regarding the causes of ND outbreaks in vaccinated commercial poultry farms. Knowing the genomic characteristics of Newcastle disease virus (NDV) infecting commercial poultry operations in spite of vaccination might give important insights on the infection dynamics of these viruses. In addition, molecular analyses at the subgenotype level and studies on the relationship of Japanese NDVs with other isolates from around the world are lacking. Therefore, in the present study, a molecular epidemiological investigation was conducted to characterize nine NDVs isolated from vaccinated commercial poultry flocks in five different Prefectures in non-epidemic areas of Japan between 1969 and 2002.
Nucleotide sequencing and phylogenetic studies were performed to characterize the complete fusion (F)-protein gene, 3-prime end of the nucleoprotein (NP)-gene and 5-prime end of the RNA dependent RNA polymerase (L)-gene. Sequence data were compared with 180 NDV strains from GenBank representing different NDV genotypes and subgenotypes from different regions of the world at different time periods. Deduced amino acids were analyzed for homologies, recombination and mutation. Recombination events were estimated using Recombination Detection Program (RDP) version 3.44. Phylogenetic trees were constructed to determine evolutionary relationships among strains.
Mean death time (MDT: 48-56 hr), Intracerebral Pathogenicity Index (ICPI: 1.7-1.9) and deduced amino acid sequences of the F0 proteolytic cleavage site (112RRQKR116) revealed that all nine field isolates were velogenic. Phylogenetic analysis showed that these isolates could be classified into two genetic lineages and three sublineages namely genotypes VIa (lineage 4a), VId (lineage 4d) and VIId (lineage 5d). No recombination events were observed but a point mutation in one of the neutralizing epitope of the F-protein was identified in the field isolates from Japan.
All field isolates from vaccinated commercial poultry in non-epidemic areas of Japan were part of much bigger outbreaks in provinces and regions and, in some cases, continents. In general, four ND panzootics occurred in Japan and that these outbreaks were mostly characterized by co-circulation of genetically distinct virus lineages due to involvements of infected wild birds. The point mutation identified in the field isolates from Japan may be due to escape from vaccine pressure. The identification of such mutation may be useful for future site-directed mutagenesis to understand the dynamics of NDV infection in vaccinated chickens.
Newcastle Disease (ND) is a highly contagious and economically devastating disease of poultry. It is caused by the Newcastle disease virus (NDV) [also known as avian paramyxovirus type-1 (APMV-1)] of the genus Avulavirus of the family Paramyxoviridae. NDV infects a wide range of domestic and wild bird species worldwide. Among animal viruses, it is one of the biggest contributors of economic losses to the world’s economy [1, 2].
NDV is an enveloped, non-segmented, single-stranded, negative-sense RNA virus with a helical morphology. Its genome has six open reading frames (ORF) in the order of 3′-NP-P-M-F-HN-L-5′. These genes encode for the following proteins: nucleoprotein (NP), phosphoprotein (P), matrix protein (M), fusion protein (F), hemagglutinin-neuraminidase (HN) and the RNA dependent RNA polymerase (L) respectively. During P-gene transcription, two additional non-structural proteins, the V and the W proteins, are also generated through RNA editing . Based on genomic size and the nucleotide sequences of the F and L genes, NDV strains can be categorized as class I or class II . Class I NDVs, which have a genomic size of 15,198 nucleotides , are occasionally isolated from wild aquatic birds and domestic poultry and are mostly avirulent to chickens. Class II NDVs comprise the majority of virulent NDV strains and some avirulent NDV strains . Class II NDVs are further subdivided into 11 genotypes (I-XI) [1, 5–10]. Early sublineages of Class II NDVs that occurred before the 1960s (genotypes I to IV) have a genomic size of 15,186 nucleotides, whereas late Class II NDV sublineages (genotypes VI to XI) have a genomic size of 15,192 nucleotides. Class II NDVs under genotype VI and VII are further subdivided into eight (a-h) subgenotypes [5–10]. Aldous et al. proposed the creation of lineages and sublineages in classifying NDVs to make it possible to rapidly type future virus isolates on the basis of their nucleotide sequence and make inferences about their origins. They proposed that NDVs could be divided into six broadly distinct groups (lineages 1 to 6), where lineages 3 and 4 were further subdivided into four sublineages (a to d) and lineage 5 was further subdivided into 5 sublineages (a to e). Genotypes I and II correspond to lineage 1 and 2, while genotype III corresponds to sublineage 3a; genotype IV to sublineage 3b; genotype V to sublineage 3c; genotype VIII to sublineage 3d; genotype VIa and VIe to sublineage 4a; genotypes VIb, VIc and VId to sublineages 4b, 4c and 4d; genotypes VIIa, VIIb, VIIc and VIId to sublineages 5a, 5b, 5c and 5d; 5e to previously characterized genotype VII NDVs from Taiwan and a quarantine isolate in UK that formed a separate cluster from other lineage 5 NDVs; and lineage 6 represents a new NDV genogroup. Recently a novel lineage, provisionally named lineage 7 was reported in West and Central Africa .
Different NDV strains vary greatly in pathogenicity [12–14]. NDV isolates can be broadly grouped into five pathotypes on the basis of clinical signs in infected chickens. ND may manifest as viscerotropic velogenic, neurotropic velogenic, mesogenic, lentogenic and asymptomatic enteric [14, 15]. Other factors, such as host species, host immune status and age, environmental stress, coinfection with other organisms, viral dose and route of exposure, may also influence the severity of the disease .
In Japan, ND was first reported during the first panzootic in the 1930s. This panzootic was caused by a genotype III NDV. After this time, large ND outbreaks were reported to occur until commercial vaccines became available in the late 1960s . Since then, sporadic outbreaks, mostly in small unvaccinated backyard flocks and pet birds, have been reported . In spite of vaccination, few sporadic outbreaks in vaccinated commercial poultry have also been observed .
At present, limited molecular epidemiological data are available regarding the causes of ND outbreaks in vaccinated poultry farms. Knowing the molecular characteristics of NDV strains affecting commercial poultry in spite of vaccination might give important insights on the possible origins and genetic nature of these viruses which may help in formulating more effective ND prevention and control strategies. In addition, no studies have been performed yet investigating the classification of Japanese NDVs at the subgenotype level and if recombination events occur in Japanese NDVs. Knowing the subgenotype classification of NDVs and occurrence of recombination events are essential since these may provide a more direct understanding on the epidemiological relationship of Japanese NDVs with other strains from different parts of the world, which may help further elucidate the mechanisms of global and transcontinental dynamics of transmission and spread of this disease. Therefore in the present study, field strains of NDVs with different geographical and temporal distribution patterns that were isolated from vaccinated commercial poultry flocks in non-epidemic areas of Japan were analyzed. Sequence data were extensively compared with 180 NDV strains from different parts of the world from different time periods.
Biological and pathotypical characterizations
All suspected NDV isolates yielded a 766-bp product in the nested PCR amplification step. This confirmed that all isolates belong to Avian Paramyxovirus type-1 viruses. All strains exhibited mean death time (MDT) of 48 to 56 hours in embryonated chicken eggs. Intracerebral pathogenicity index (ICPI) values ranged from 1.7 to 1.9 (Table 1). Nucleotide sequence analyses of the variable region of the F-gene (47–421 nt) showed that all isolates had multiple basic amino acids at its F0 proteolytic cleavage site (residues 112–117). Predicted amino acid sequence of the F0 cleavage site was 112RRQKR116 at the F2 protein and phenylalanine (116 F117) at the N-terminus of the F1 protein for all strains. These results indicated that all isolates were velogenic.
Genetic and phylogenetic characterizations
A total of 1662 nucleotides encoding for 553 amino acid residues were identified in the complete coding region of the F-gene of all field strains. Six potential N-glycosylation sites (Asn-X-Ser/Thr where X is any amino acid except proline or aspartic acid) located at positions 85 to 87, 191 to 193, 366 to 368, 447 to 449, 471 to 473 and 541 to 543 were recognized. Twelve cysteine residues located at positions 25, 76, 199, 338, 347, 362, 370, 394, 399, 401, 424 and 523 were identified. Comparison of glycosylation sites and cysteine residues showed no changes in the amino acid sequence in all field strains, which may indicate that these sites were highly conserved.
Analysis of seven neutralizing epitopes located at positions 72, 74, 75, 78, 79, 157 to 171 and 343 of the F-protein showed a K to R amino acid substitution at position 78 in seven of the nine strains (Table 2). These strains were JP/Ibaraki/SG106/99, JP/Ibaraki/IS5/02, JP/Chiba/BY103/01, JP/Ibaraki/IS2/02, JP/Miyagi/AGT/02, JP/Chiba/BY7/02 and JP/Fukushima/NYF-3/02. Furthermore, analysis of amino acid substitutions showed 13 point mutations in the variable region of the F-gene of the Japanese field strains (Table 2). Genotypic and subgenotypic-specific amino acid substitutions were also observed, which were consistent with the proposed theory of NDV evolution as reported previously .
Phylogenetic analyses of nine field strains and 180 NDV strains from GenBank were performed by using contiguous nucleotide sequences of the F-gene, NP-gene and L-gene. Reference strains from GenBank were selected as representatives of nine of the 11 ND genotypes (genotypes I to XI) representing isolates from different regions of the world (Figures 1, 2 and 3).
Field NDV strains were observed to belong to two distinct genotypic groups (genotype VI and VII) using phylogenetic analysis of the complete coding sequence of the F-gene. Early isolates such as JP/Osaka/2440/69 and JP/Ibaraki/SM87/87 were genotype VI, while field strains isolated from 1999 onwards were from genotype VII (Figures 1 and 4). Phylogenetic analysis using the 3-prime portion of NP-gene and 5-prime portion of L-gene yielded the same tree topology and phylogenetic groupings (Figures 2 and 3). Phylogenetic analyses at the subgenotype level using the variable region of the F-gene sequences revealed that JP/Osaka/2440/69 belongs to subgenotype VIa, JP/Ibaraki/SM87/87 to subgenotype VId and all the recent field isolates to subgenotype VIId (Figures 5 and 6).
F-gene nucleotide sequence of JP/Osaka/2440/69 was found to be closely related (98.5-99.2% nucleotide sequence identity) to isolates from the Middle East. The F-gene sequence of JP/Ibaraki/SM87/87 was closely related (96.3-100.0%) with isolates from Japan and China. All the other VIId isolates were highly similar (97.2-100.0%) to isolates from Japan, China and Taiwan and interestingly to a goose isolate from China (98.4-98.6%) (Table 3).
No intragenic nor intergenic recombination events (Unique events = 0; Recombination signals = 0) were observed involving the field isolates using all the described methods. Analysis of the over-all mean evolutionary distance among the Japanese field isolates showed rates of 3.9×10-2 [standard error (s.e) 0.01], 3.1×10-2 (s.e. 0.01) and 2.0×10-2 (s.e. 0.01) base substitution per site in the full F-gene, partial NP-gene and partial L-genes, respectively. In contrast, over-all mean evolutionary distance among the recent field strains (1999–2002) were 1.0×10-3 (s.e. 0.001) base substitution per site for the F and NP-genes and zero base substitution for the L-gene. In addition, 25 sites in the F-gene of all field strains were observed to be under negative selection (p-value < 0.05) (Table 4).
ND remains a serious threat to commercial poultry even though intensive vaccination programs are being applied. In Japan, occasional outbreaks have been reported in commercial poultry mostly due to improper vaccination, immunosuppression due to infectious and non-infectious causes, and challenge by more velogenic viruses . However, limited data are available regarding the genomic characteristics of NDVs occurring in vaccinated commercial poultry flocks. Knowing the genetic characteristics of wild strains of NDV affecting vaccinated poultry might give important insights on the possible origins, transmission mechanisms and infection routes of these viruses. Molecular and phylogenetic studies like this are important since these might lead to better understanding on how to prevent, control and manage future ND cases.
Molecular characterization of NDV strains mostly considered the F-gene with particular emphasis given on the variable region (47–421 nt) because it codes for a number of functionally important structures such as signal peptide [amino acid (aa) 1–31], cleavage activation sequence (aa 112–116), portion of the fusion inducing hydrophobic region (aa 117–142) and it is characterized by both variable and conserved regions [19, 20]. Nucleotide sequence of the F-gene fragment (nt 47–420) is regarded as standard criterion for genotyping . A molecular basis of pathogenicity has also been well established through sequence analysis of F-protein cleavage site. It was reported that the motiff 112R/K-R-Q-K/R-R116 at the C-terminus of the F2 protein and F (phenylalanine) at the N-terminus of the F1 protein (residue 117) are major determinants of viral virulence [1, 2, 18, 22–25]. A huge database of sequence data especially on F-gene sequences of NDVs isolated throughout the world has also been published and available for sequence comparison and phylogenetic studies .
Records of management and farm history showed that NDV strains used in this study originated from farms with diverse geographical, temporal and disease profiles. It is noteworthy to emphasize that these farms were thoroughly vaccinated against NDV but they were still infected with the disease. Moreover, deduced amino acid sequence of cleavage site of the F-gene of all field isolates revealed the motif 112R-R-Q-K-R-F117 indicating that all strains were velogenic. This was further confirmed by MDT and ICPI tests. These indicate that in spite of the regular use of inactivated and live vaccines, velogenic ND may still occur in vaccinated flocks. However, some affected birds showed only mild respiratory symptoms without significant mortalities and severe pathological lesions. In most flocks, only mild to moderate decrease in egg production was observed.
Seven major epitopes have been identified involving the fusion inhibition and neutralization of F-protein [19, 27, 28]. Individual amino acids at 72, 74, 75, 78, 79 and 343 and a stretch of amino acids from residues 157–171 were identified to be critical for both structures and functions of the F-gene. In this study, nucleotide substitution in one of the fusion inhibition and neutralizing epitope (p.K78R) was identified in all of the seven VIId strains. 13 point mutations were also identified in the variable region of the F-gene. Comparison with sequence data from reference strains (n = 180) showed that among these mutations, p.K4I were conserved only in NDV strains originating from Japan while p.L21P, p.I52V, p.K78R and p.R101K were conserved in strains originating from the Far East Asia (Japan, China and Taiwan). These substitutions maybe used as crude molecular markers of geographic origins of NDVs.
Analysis of genotypic and subgenotypic substitutions in the hypervariable region of the F-gene showed findings that were in conformity with the proposed theory of NDV genetic evolution . It was proposed that subgenotype VIa from the second pandemic probably evolved to VIc by production of a crucial p.S107T substitution and VId by production of p.S93T substitution; VIIb evolved from VIb via a VII-specific p.V121I substitution; VIIb evolved to become VIIa and VIIc through p.K101R substitution; and VIIc evolved to become VIId by the production of additional p.I52V and p.F314Y substitution .
A point mutation in the F-gene may have resulted to neutralizing epitope variants. However whether this mutation was part of adaptive mechanism of NDVs to evade the immune response to be able to infect vaccinated chickens is unclear or whether this mutation was actually the effect of selective immune pressure exerted on ND viral particles as a consequence of vaccination is also unknown. To understand how wild NDVs infect vaccinated chickens, this identified mutation may be useful for future site-directed mutagenesis studies.
Phylogenetic analyses on the field strains using the variable region of the F-gene (47–421 nt) revealed that JP/Osaka/2440/69 belongs to genotype VIa (Figure 5). As reported previously , genotype VIa was responsible for the second ND panzootic that started in the Middle East during the 1960s and then spread to most countries around the world as a result of enormous trade and importation of captive caged birds and technological advances in air transportation. Interestingly, this isolate shared 100% sequence identity with JP/Narashino/68, which was isolated from a Japanese Blue Magpie. It is interesting to note that nucleotide sequence identity of Iraq/AG-68 was 98.5% similar while strains Kuwait/256/68, Lebanon/70 and Israel/70 were 98.7-99.2% similar with these Japanese strains. It is possible that JP/Narashino/68 was a foreign strain that was introduced to Japan from wild birds. JP/Narashino/68 and/or its progenitor might have been then spread to domestic chickens in Japan, leading to the isolation of JP/Osaka/2440/69 (100% similar) (Table 2).
Phylogenetic analyses on JP/SM87/87 showed that this strain belongs to VId ND viruses (Figure 5). VId viruses together with VIb and VIc were responsible for the third panzootic, which were reported to be spread by pigeons. This strain shared 100% sequence identity with JP/Tochigi/85 and JP/Ibaraki/85, which may indicate that JP/SM87/87 was a product of the ongoing outbreak.
The seven remaining field isolates belong to VIId (Figure 6). Genotype VII is the most predominant NDV genotype that is responsible for most outbreaks in East Asian countries including Taiwan, Korea and China since the 1980s, constituting the fourth pandemic [1, 5–7, 18, 30, 31]. Also in Japan, the isolation of genotype VII viruses was reported previously . Therefore, this genotype has been the most predominant NDV in recent outbreaks in Japan.
The earliest VIId viruses on record infected chickens from South Korea in 1995 (Figure 6). These strains include Kr-279/95, Kr-146/95 and Kr-077/95. On the other hand, the earliest VIId NDVs that were reported from Japan were JP/Tokyo/96 from chickens and JP/Ibaraki-ph/97 from a pheasant (99.0% similar to one another), which might indicate that the two strains were part of an ongoing outbreak. Remarkably, JP/Tokyo/96 shared 99.0% sequence identity with GX-1/97, which was isolated from a chicken flock in Western China, FJ-2/99 from a fowl from China and GD/1/98/Go from a goose from China (Table 2). Remarkably, GD/1/98/Go was 98.4-98.6% similar with JP/Ibaraki/SG106/99 and all the other VIId field strains in this study. These findings may indicate that wild birds have played a role in the circulation of VIId viruses across the Far East Asian countries (Korea, Japan and China). A comparison of homologies with contemporary isolates showed that JP/Ibaraki/SG106/99 and all the other VIId field strains were highly similar (99-100%) with JP/Ibaraki/00. Interestingly, a p.K78R amino acid substitution in the F-protein of this strain was reported previously . It was observed that chickens that were challenged with JP/Ibaraki/00 survived a cross-protection test after vaccination with B1 strain, however it was noted that vaccination did not prevent infection and excretion of the virus . This result was partially correlated with the clinical profile of the infected flocks seen in this study. Although the infected flocks survived the infection, problems with production performance were observed.
Comparison among the recent Japanese field strains showed that these strains have high F-gene homologies (99.7-100%), which may indicate that these strains may have been epidemiologically related. Computation of the over-all mean evolutionary distance of the F-gene of these recent strains showed a substitution rate of 1.0×10-3 (s.e. 0.001) base substitution per site in the span of 3 years (1999–2002). In contrast, over-all F-gene mean evolutionary distance in all field strains (1969–2002) was 3.9×10-2 (s.e. 0.010) base substitution per site. Interestingly, the partial NP gene showed an almost same substitution rate (1.0×10-3 and 3.1×10-2) while partial L-gene had the lowest rate of substitution (zero and 2.0×10-2). However, direct comparison among these substitution rates is not feasible because of incomplete sequence data in NP and L-genes in this study. In other reports, it was shown that among the NDV proteins, the F and P-protein have the highest rate of change (0.78-1.98×10-3 and 0.78-2.32×10-3 substitution/site/year, respectively) while L and NP-proteins have the lowest rate of change (0.59-1.44×10-3 and 0.45-1.50×10-3 substitution/site/year, respectively) [32, 33]. In addition, analysis of evolutionary selection profiles of the Japanese field strains revealed 25 sites in the F-gene sequence that were under negative selection (p-value < 0.05) (Table 4). No positive selection sites were identified. This is in agreement with the findings of other authors that over-all NDV proteins are under strong purifying and negative selection pressures [32, 33].
Phylogenetic analyses using the nucleotide sequences of NP gene, L-gene, complete F-gene coding sequence and variable region of the F-gene resulted to almost similar tree topologies. Surprisingly phylogenetic analysis using NP gene resulted to a clearer differentiation among field strains. These may indicate that NP and L-genes may be alternative methods to characterize NDVs given that like the F-gene, these genes are also involved in the dynamics of viral virulence, play important functional roles in the NDV replication cycle; and are also characterized by regions with high conservation necessary to identify homologies among strains but also characterized by regions with high variations necessary to identify specific variations between strains. Phylogenetic analysis of NP and L-genes in conjunction with the F-gene may also help detect possible natural or artificial recombination events.
This investigation showed that all field isolates from vaccinated commercial poultry were part of much bigger outbreaks affecting not only provinces or regions but even entire continents. To determine how commercial farms are being infected with NDV, the epidemiology of NDV in the whole of Japan and parts of Far East Asia was analyzed. This study showed that Japanese poultry was affected by at least four pandemics and that outbreaks were mostly characterized by co-circulation of genetically distinct virus lineages that were consistent with the predominant virus genotype circulating in a particular time period. Moreover, no distinct transition was observed from each pandemic. Aside from involvement of local strains, ND outbreaks in Japan were mostly due to virus transmission from infected wild birds either by international bird trade or migration patterns.
A point mutation in one of the neutralizing epitopes of the F-protein resulting to occurrence of neutralizing epitope variants was identified. This identified mutation may be useful for future site-directed mutagenesis to understand the dynamics of NDV infection in vaccinated chickens.
Materials and methods
Nine NDV strains isolated from commercial poultry farms with different spatial (Osaka, Ibaraki, Chiba, Fukushima and Miyagi Prefectures) and temporal distribution patterns (1969, 1987, 1999, 2001 and 2002) were used to investigate the molecular epidemiological relationships of ND outbreaks in vaccinated commercial poultry flocks in Japan. These strains were isolated by one to two passages of pooled infected tissues in 10-day-old embryonated specific pathogen-free (SPF) chicken eggs. Infective allantoic fluids were harvested and kept in lyophilized form or in serum tubes and stored at −80°C until further use. All isolates used in this study were provided by Poultry Products Quality Control Co. Ltd. (Fukushima, Japan).
Farm history and clinical profile
Records of management and farm history were obtained to characterize the clinical profile of nine suspected NDV strains (Table 1). The oldest strain was from Osaka Prefecture in 1969 (JP/Osaka/2440/69). This isolate was recovered from six dead layer birds that were submitted to the Osaka Veterinary Municipal Office for diagnosis. The flock was vaccinated with the live B1 vaccine in drinking water at 5–7 days of age and killed ND vaccine (Sato strain) at 25 days. The affected flock was around 80 days of age when the disease occurred. Mortality was reported to be around 60-70% but no data were given regarding observed clinical signs and production performance. The second strain was from a layer farm in Ibaraki Prefecture isolated in 1987 (JP/Ibaraki/SM87/87). Total population of affected farm was 33,360 birds in 12 open-type houses of 2,780 birds each. The disease was reported in one of the houses. The flock was vaccinated with the live B1 spray at 10 days of age, killed Ishii/B1 at 45 days, live B1 spray again at 60 days and killed Ishii/B1 again at 120 days. The disease occurred two months after the last vaccination at 180 days of age. The disease was characterized by gasping and 25% drop in egg production. Mortality was less than three percent. The third strain was from Ibaraki Prefecture in 1999 that was isolated from dead spent hens sent for diagnosis (JP/Ibaraki/SG106/99). The isolate was from an unknown farm raising spent hens for liquid egg production. No detailed farm information was obtained. Chickens submitted for diagnosis were approximately more than 700 days old. Necropsy findings were proventriculitis, hemorrhagic lesions in duodenum and petechiae in lymphocytic tissues. The fourth strain was from a replacement pullet farm in Chiba Prefecture isolated in 2001 (JP/Chiba/BY103/01). The disease occurred in a flock of 21,000 birds. The flock was vaccinated with the live B1 strain in drinking water at 4 and 10 days of age, live B1 spray at 28 days and killed Ishii strain at 45 and 90 days of age. It was reported that the killed vaccines were injected by hired professional vaccination staff that travel from farm to farms. Six days after last vaccination, the disease occurred characterized by gasping, nervous symptoms, leg weakness, twisting of neck and greenish diarrhea. Mortality was around 10%. The fifth strain was from a layer farm in Ibaraki prefecture isolated in 2002 (JP/Chiba/BY7/02). The affected farm had a population of 125,000 birds (five houses of 25,000 birds) and affected flock was around 25,000 birds (one house). The flock was vaccinated with the live B1 strain in drinking water at 10 days of age, live B1 spray at 24 days, killed Ishii strain at 45 days, live B1 spray again at 60 days and killed Ishii strain again at 95 days. The disease occurred at 110 days of age characterized by gasping, swollen face and infectious bronchitis (IB)-like respiratory signs and 20% drop in egg production without significant mortalities. Necropsy findings were necrotic ovarian follicles and necrotic catarrhal inflammation of the intestines. The sixth strain was from a layer farm located in Ibaraki prefecture in 2002 (JP/Ibaraki/IS5/02). The affected farm had a population of approximately 200,000 birds and the affected flock was around 41,000 birds. The flock was vaccinated with the live B1 spray at 10 and 28 days old, killed Ishii strain at 45 days, live B1 spray again at 60 days and killed Ishii strain again at 90 days. The disease occurred at 336 days of age characterized by mild respiratory signs such as gasping, seven percent decrease in egg production with no marked mortalities. The seventh strain was recovered from dead birds from a replacement pullet farm in Ibaraki prefecture in 2002 (JP/Ibaraki/IS2/02) that were submitted for diagnosis. Total farm population was around 10,000 birds. The flock was vaccinated with the live B1 spray at 10 days of age. The disease occurred four days after vaccination characterized by severe depression. The eight strain was isolated from a layer farm in Fukushima prefecture (JP/Fukushima/NYF-3/02). The farm had a population of around 120,000 birds and affected flock was approximately 16,000 layers. The flock was vaccinated with the live B1 spray at 10 and 28 days of age, killed Ishii strain at 45 days, live B1 spray again at 60 days and killed Ishii strain again at 90 days. The disease occurred around 532 days of age characterized by mild gasping, increase in soft-shelled eggs and eight percent decrease in egg production. Mortalities were minimal and within production standards. The ninth strain was isolated from a layer farm in Miyagi prefecture in 2002 (JP/Miyagi/AGT/02). The flock was vaccinated with the live B1 in drinking water at 10 days of age, live B1 spray at 24 days, killed Ishii strain at 45 days, live B1 spray again at 60 days and killed Ishii strain at 95 days. It was informed that the flock was kept in multiple age houses (six different flocks in one house). Disease occurred at 250 days of age characterized by gasping, greenish white diarrhea and 70% decrease in egg production. No marked mortalities were observed.
Biological and pathotypical characterizations
Biological and pathotypical characterization of isolates were performed using mean death time (MDT) in 10-days-old embryonated SPF chicken eggs and intracerebral pathogenicity index (ICPI) in 1-day old chicks according to the protocols described in the OIE manual . Confirmation of pathotypes was performed by nucleotide sequence analysis of the F0 proteolytic cleavage site (residues 112–117). Pathotypical data for strain JP/Fukushima/NYF3/02 were obtained in our previous study.
Reverse transcription polymerase chain reaction (RT-PCR)
Nested RT-PCR was performed to confirm the identities of suspected NDV strains. In brief, isolates were propagated once in 10-day-old embryonated SPF eggs. Viral RNA from infected allantoic fluids was extracted directly by using a QIAamp® Viral RNA Mini Kit (Qiagen, West Sussex, UK). Viral RNA was transcribed to cDNA by using random hexamers and Primescript® Reverse Transcriptase (Takara Bio-Inc, Shiga, Japan). cDNA was amplified by PCR as described previously . A two-step nested PCR was performed to amplify the region comprising the 3′ end of the M-gene and the 5′ end of the F-gene using KOD dash® (Toyobo, Osaka, Japan), 5 uM of external and internal primers as described by Mase et al. (Table 5). Thermocycling conditions for the first and second PCR steps were as follows: prewarming at 94°C for 2 min (1 cycle), denaturation at 94°C for 30 sec, annealing at 50°C for 30 sec and extension at 72°C for 1 min. Amplification steps were performed for 35 cycles. The final extension was performed at 72°C for 30 sec (1 cycle).
Nucleotide sequence analysis
Confirmed NDV strains were subjected to additional RT-PCR amplification to characterize the open reading frame of the F-gene sequences (4504–6295 nt of whole NDV genome), 3-prime end of the NP-gene sequences (16–783 nt of whole NDV genome) and 5-prime end of the L-gene sequences (13995–14719 nt of whole NDV genome). In brief, RT-PCR was performed by using SapphireAmp® Fast PCR Master Mix (Takara Bio), 5 uM of forward and reverse primers (Table 5) and cDNAs that were transcribed previously. Thermocycling conditions consisted of initial denaturation at 95°C for 2 min followed by 35 cycles of 98°C for 10 sec (denaturation), 55°C for 10 sec (annealing), 72°C for 10 sec (extensions) and final extension at 72°C for 2 min. PCR products were analyzed by electrophoresis with 1.2% agarose gel and purified by using QIAquick® Gel Extraction Kit (Qiagen, Valencia, CA). The nucleotide sequences of PCR products were determined by Big Dye terminator cycle-sequencing kit version 3.1 (Applied Biosystems Inc., Foster City, CA) and an ABI Prism 3130 Genetic Analyzer (Applied Biosystems). DNA products were sequenced from both directions. The complete F-gene sequence for strain JP/Fukushima/NYF3/02 was obtained in our previous study.
Sequence assembly and editing were performed using CodonCode Aligner® (version 3.7.1, CodonCode Corporation, MA) and ClustalX® (version 2.1, Conway Institute UCD Dublin, Ireland). Deduced amino acid sequences were determined using Bioedit® software package version 126.96.36.199 . Confirmation of identity and homology were performed using BLAST http://www.ncbi.nlm.nih.gov.
To determine the molecular epidemiological relationships of field strains, 180 NDV strains isolated from different regions of the world at different time periods were obtained from GenBank. These reference strains were representatives of all the different NDV genotypes and subgenotypes. Phylogenetic and molecular evolutionary analyses were conducted using MEGA version 4 . Phylogenetic trees of the variable region of the F-gene sequences (47–421 nt), complete coding region of the F-gene sequences (1–1662 nt), coding region of the 3-prime end of NP-gene sequences (1–622 nt) and coding region of the 5-prime end of L-gene sequences (5629–6333 nt) were constructed by the neighbor-joining method with the maximum composite likelihood substitution model at 1000 bootstrap replicates.
Determination of recombination events, evolutionary distances and selection profile
Intragenic recombination events in the NDV nucleotide sequences were determined using RDP v3.44 program. . Seven different algorithms integrated in the program namely RDP, GeneConv, Bootscan, MaxChi, Chimaera, SiScan and 3Seq were applied to detect any putative recombination breakpoints and to estimate the occurrence of any recombination events within all the analyzed genes. Sequences with recombination events identified by at least two detection methods (p < 0.01) were considered as true recombinants. Intragenic recombination events within the F-gene were also determined by comparing the tree topology of phylogenetic analyses using the variable region (47–421 nt) and complete coding region (1–1662 nt) of the F-gene. Intergenic recombination events were determined by comparison of topologies of the generated F, NP and L-gene phylogenetic trees.
Evolutionary distances were calculated using MEGA version 4 using the Maximum Composite Likelihood method . Codon positions included were the 1st + 2nd + 3rd + noncoding. All positions containing gaps and missing data were eliminated from the dataset (complete deletion option). Estimates of s.e. were obtained by a bootstrap procedure of 1000 replicates. Analysis of evolutionary selection profile was performed using Datamonkey http://www.datamonkey.org/ following the Fixed- Effect Likelihood (FEL) method [Hasegawa, Kishino and Yano (HKY) model, p-value less than 0.05] [38, 39].
Nucleotide sequences used in this study were submitted to the DNA Databse of Japan (DDBJ) with the following accession numbers: JP/Osaka/2440/1969 F-gene sequence: AB853926; NP-gene: AB854728; L-gene: AB854729. JP/Ibaraki/SM87/1987 F-gene sequence: AB853928; NP-gene: AB854730; L-gene: AB854731. JP/Ibaraki/SG106/1999 F-gene sequence: AB853927; NP-gene: AB854732; L-gene: AB854733. JP/Ibaraki/BY103/2001 F-gene sequence: AB853931; NP-gene: AB854734; L-gene: AB854735. JP/Ibaraki/BY7/2002 F-gene sequence: AB853930; NP-gene: AB854736; L-gene: AB854737; JP/Ibaraki/IS5/2002 F-gene sequence: AB853933; NP-gene: AB854744; L-gene: AB854745. JP/Ibaraki/IS2/2002 F-gene sequence: AB853932; NP-gene: AB854742; L-gene: AB854743; JP/Fukushima/NYF-3/2002 F-gene sequence: AB853329; NP-gene: AB854740; L-gene: AB854741; JP/Miyagi/AGT/2002 F-gene sequence: AB853929; NP-gene: AB854738; L-gene: AB854739.
Animal experiments in this manuscript were critically reviewed and approved by the Animal Ethical Committee of Poultry Products Quality Control (Permit number: 2013-03A) and performed in accordance with the Japanese laws and international guidelines for the use of animals in research.
Avian paramyxovirus type-1
Complementary deoxyribonucleic acid
Fixed- effect likelihood
Hasegawa, Kishino and Yano
Intracerebral pathogenicity index
APMV1/chicken/Japan/Ibaraki/IS5/2002, JP/Ibaraki/SG106/99, APMV1/chicken/Japan/Ibaraki/SG106/1999
RNA dependent RNA polymerase
Mean death time
Newcastle disease virus
Recombination detection program
Reverse transcription polymerase chain reaction
Specific pathogen free.
Alexander DJ: Newcastle disease, other avian paramyxoviruses and pneumovirus infection. In Diseases of poultry. 12th edition. Edited by: Saif YM. Iowa: Blackwell Publishing; 2008:75-93.
Alexander DJ: Ecology and epidemiology of Newcastle disease. In Avian influenza and Newcastle disease. Edited by: Capua I, Alexander DJ. Milan: Springer; 2009:19-26.
Steward M, Vipond I, Millar N, Emmerson P: RNA editing in Newcastle disease virus. J Gen Virol 1993,74(12):2539-2547. 10.1099/0022-1317-74-12-2539
Liu XF, Wang X, Wu S, Hu S, Peng Y, Xue F, Liu X: Surveillance for avirulent Newcastle disease viruses in domestic ducks ( anas platyrhynchos and cairina moschata ) at live bird markets in eastern china and characterization of the viruses isolated. Avian Pathol 2009, 38: 377-391. 10.1080/03079450903183637
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
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
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.
Tsai H, Chang JKH, 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
Cattoli G, Fusaro A, Monne I, Molia S, Le Menach A, Maregeya B, Nchare A, Bangana I, Maina AG, Koffi JN, Thiam H, Bezeid OE, Salviato A, Nisi R, Terregino C, Capua I: Emergence of a new genetic lineage of Newcastle disease virus in west and central Africa – implications for diagnosis and control. Vet Microbiol 2010,142(3–4):168-176.
Maminiaina OF, Gil P, Briand FX, Albina E, Keita D, Andriamanivo H, Chevalier V, Lancelot R, Martinez D, Rakotondravao R, Rajaonarison J, Koko M, Andriantsimahavandy A, Jestin V, de Almeida R: Newcastle disease virus in Madagascar: identification of an original genotype possibly deriving from a died out ancestor of genotype IV. PLoS ONE 2010, 5: e13987. doi: 10.1371/journal.pone.0013987 10.1371/journal.pone.0013987
Aldous EW, Mynn JK, Banks J, Alexander DJ: A molecular epidemiological study of avian paramyxovirus type 1 (Newcastle disease virus) isolates by phylogenetic analysis of a partial nucleotide sequence of the fusion protein gene. Avian Pathol 2003, 32: 239-257. 10.1080/030794503100009783
Islam M, Ito T, Takakuwa H, Takada A, Itakura C, Kida H: Acquisition of pathogenicity of a Newcastle disease virus isolated from a Japanese quail by intracerebral passage in chickens. Jpn J Vet Res 1995, 42: 147-156.
Miller P, Decanini L, Afonso C: Newcastle disease: evolution of genotypes and the related diagnostic challenges. Infect Genet Evol 2010, 10: 26-35. 10.1016/j.meegid.2009.09.012
Tsunekuni R, Ito H, Kida H, Otsuki K, Ito T: Increase in the neuraminidase activity of a nonpathogenic Newcastle disease virus isolate during passaging in chickens. J Vet Med Sci 2010, 72: 453-457. 10.1292/jvms.09-0474
Alexander DJ: Newcastle disease and other paramyxovirus infections. In Diseases of poultry. 9th edition. Edited by: Calnek BW, Barnes HJ, Beard CW, Reid WM, Yoder HW. Iowa: Iowa State University Press; 1991:496-519.
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
Nakamura K, Ohtsu N, Nakamura T, Yamamoto Y, Yamada M, Mase M, Imai K: Pathologic and immunohistochemical studies of Newcastle disease in broiler chickens vaccinated with ND: severe nonpurulent encephalitis and necrotizing pancreatitis. Vet Pathol 2008,2008(45):928-933.
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.2001
Toyoda T, Gotoh B, Sakaguchi T, Kida H, Nagai Y: Identification of amino acids relevant to three antigenic determinants on the fusion protein of Newcastle disease virus that are involved in fusion inhibition and neutralization. J Virol 1988, 62: 4227-4430.
Yussof K, Tan WS: Newcastle disease virus: macromolecules and opportunities. Avian Pathol 2001, 30: 439-455. 10.1080/03079450120078626
Qin Z, Tan L, 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(2):601-611. 10.1128/JCM.01356-07
Nagai Y, Klenk HD, Rott R: Proteolytic cleavage of the viral glycoproteins and its significance for the virulence of Newcastle disease virus. Virology 1976,72(2):494-508. 10.1016/0042-6822(76)90178-1
Gotoh B, Ohnishi Y, Inocencio NM, Esaki E, Nakayama K, Barr PJ, Thomas G, Nagai Y: Mammalian subtilisin-related proteinases in cleavage activation of the paramyxovirus fusion glycoprotein: superiority of furin/PACE to PC2 or PC1/PC3. J Virol 1992,66(11):6391-6397.
Ogasawara T, Gotoh B, Suzuki H, Asaka J, Shimokata K, Rott R, Nagai Y: Expression of factor X and its significance for the determination of paramyxovirus tropism in the chick embryo. EMBO J 1992,11(2):467-472.
Mase M, Murayama K, Karino A, Inoue T: Analysis of the fusion protein gene of Newcastle disease viruses isolated in Japan. J Vet Med Sci 2011, 73: 47-54. 10.1292/jvms.10-0281
Jindal N, Chander Y, Chockalingam AK, de Abin M, Rdig PT, Goyal SM: Phyogenetic analysis of Newcastle disease viruses isolated from waterfowl in the upper Midwest region of the United States. Virol J 2009, 6: 191. 10.1186/1743-422X-6-191
Neyt C, Geliebter J, Slaoui M, Morales D, Meulemans G, Burny A: Mutations located on both F1 and F2 subunits of the Newcastle disease virus fusion protein confer resistance to neutralization with monoclonal antibodies. J Virol 1989, 63: 952-954.
Yusoff K, Nesbit M, Mc Cartney H, Meulemans G, Alexander DJ, Collins MS, Emmerson PT, Samson AC: Location of neutralizing epitopes on the fusion protein of Newcastle disease virus strain beaudette C. J Gen Virol 1989, 70: 3105-3109. 10.1099/0022-1317-70-11-3105
Alexander DJ: Historical aspects. In Newcastle disease. Edited by: Alexander DJ. Boston: Kluwer Academic Publishers; 1988:1-10.
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
Lien Y, Lee J, Su H, Tsai H, Tsai M, Hsieh C, Tsai S: 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
Chong Y, Padhi A, Hudson P, Poss M: The effect of vaccination on the evolution and population dynamics of avian paramyxovirus-1. PLoS Pathog 6(4): doi:10.1371/journal.ppat.1000872
Miller P, Mia Kim L, Ip H, Alfonso C: Evolutionary dynamics of Newcastle disease virus. Virology 2009, 391: 64-72. 10.1016/j.virol.2009.05.033
Office International des Epizooties (OIE): Newcastle disease. Manual of diagnostic tests and vaccines for terrestrial animals Chapter 2.1.14. 2010. [http://www.oie.int/fileadmin/Home/eng/Health_standards/tahm/2.03.14_NEWCASTLE_DIS.pdf Chapter 2.1.14. 2010. [
Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for windows 95/98/NT. Nucleic Acids Symp Ser 1999, 41: 95-98.
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/msm092
Martin DP, Lemey P, Lott M, Moulton V, Posada D, Lefeuvre P: RDP3: a flexible and fast computer program for analyzing recombination. Bioinformatics 2010, 26: 2462-2463. 10.1093/bioinformatics/btq467
Delport W, Poon A, Frost S, Kosakovsky PS: Datamonkey 2010: a suite of phylogenetic analysis tools for evolutionary biology. Bioinformatics 2010, 26: 2455-2457. 10.1093/bioinformatics/btq429
Kosakovsky Pond S, Frost S: Datamonkey: rapid detection of selective pressure on individual sites of codon alignments. Bioinformatics 2005, 21: 2531-2533. 10.1093/bioinformatics/bti320
Felsenstein J: Confidence limits on phylogenies: an approach using the bootstrap. Evolution 1985, 39: 783-791. 10.2307/2408678
Tamura K, Nei M, Kumar S: Prospects for inferring very large phylogenies by using the neighbor-joining method. Proc Natl Acad Sci U S A 2004, 101: 11030-11035. 10.1073/pnas.0404206101
We thank Mrs. Masako Katoh and the PPQC Scholarship Foundation for funding this research. We also thank Dr. Koichi Otsuki for the constructive review and valuable recommendations during the preparation of the manuscript and to Ms. Makiko Ishii, Ms. Misa Saito, Ms. Yukari Mutoh and Dr. Fletcher del Valle for the excellent technical assistance.
The authors declare that they have no competing interests.
HK, TI and HI conceived and designed the research. HI and DVU performed RT-PCR, nucleotide sequencing and generation of phylogenetic trees. HK, TI, HI and DVU analyzed the data: HK, KS and TS collected and provided the NDV samples: HK, TI, HI and DVU drafted the manuscript. All authors’ read and approved the final manuscript.
Authors’ original submitted files for images
About this article
Cite this article
Umali, D.V., Ito, H., Suzuki, T. et al. Molecular epidemiology of Newcastle disease virus isolates from vaccinated commercial poultry farms in non-epidemic areas of Japan. Virol J 10, 330 (2013). https://doi.org/10.1186/1743-422X-10-330
- Newcastle disease virus
- RNA dependent RNA polymerase
- Neutralizing epitope variant
- Phylogenetic analysis