- Research article
- Open Access
The deletion of an extra six nucleotides in the 5′ -untranslated region of the nucleoprotein gene of Newcastle disease virus NA-1 decreases virulence
BMC Veterinary Research volume 10, Article number: 964 (2014)
The virulent Newcastle disease virus (NDV) strain NA-1 (genotype VII) was isolated from an epizootic involving a flock of geese in Jilin Province, Northeast China, in 1999. Compared with the classical NDV strains, which have a genome size of 15,186 bp, the more recently isolated NDV strains, including that involved in the goose outbreak, have an extra six nucleotides in the 5′-untranslated region (UTR) of the nucleoprotein (NP) gene. This extra sequence, TCCCAC, is highly conserved and has been found in multiple NDV strains, including ZJ-1, WF00G, JSG0210, and NA-1. In the current study, an infectious clone from strain NA-1 was isolated and designated rNA-1. Subsequently, strain rNA-1 was mutated to delete the six-nucleotide insertion, producing strain rNA-1(−). Virulence of the recombinant virus was then assayed in chickens and geese.
The recombinant virus rNA-1(−) showed similar biological characteristics to the parental NA-1 strain in DF-1 chicken fibroblast cells. However, the virulence of rNA-1(−) in 9-day-old embryonated chicken eggs and 1-day-old specific pathogen-free (SPF) chickens was decreased compared with the rNA-1 control. Furthermore, the virulence of the recombinant strain was slightly decreased in 1-day-old SPF chickens when compared with that in 1-day-old geese.
Following deletion of six nucleotides in the 5′-UTR of the NP gene of NDV strain NA-1, the virulence of the rNA-1(−) recombinant strain was decreased in both chickens and geese. However, rNA-1(−) was more virulent in chickens than geese, as seen by the higher intracerebral pathogenicity index value.
Newcastle disease (ND), caused by Newcastle disease virus (NDV), is one of the most serious avian diseases worldwide. NDV, also known as avian paramyxovirus serotype-1 (APMV-1), is a member of the genus Avulavirus, within the Paramyxoviridae family ,, and is a single-stranded, negative-sense, non-segmented, enveloped RNA virus . The NDV genome is composed of six genes, encoding six structural proteins: nucleoprotein (NP), phosphoprotein (P), matrix (M), fusion (F), hemagglutinin-neuraminidase (HN), and RNA polymerase (L) . RNA editing of P produces the additional non-structural proteins V, and possibly, W . The genomic RNA associates with the NP, P, and L proteins to form the ribonucleoprotein complex, the machinery for RNA synthesis . Although all NDV strains belong to a single serotype , there is both genetic and antigenic diversity between the strains, and the detection of progressive changes in strains isolated over successive years indicates continuous evolution ,.
Based on standard pathotyping assays, NDV strains are classified into three main phenotypes: velogenic, mesogenic, and lentogenic . For decades, even the most virulent NDV strain was not considered to be a cause of serious disease in waterfowl. However, at least four NDV strains have been implicated in fatal diseases of geese in many provinces of South, East, and Northeast China since 1997 -. Several studies have shown that conventional live-attenuated NDV vaccines can protect birds against morbidity and mortality, but because of genotype diversity, do not prevent infection and virus shedding -. Untranslated regions (UTRs) of viral genomes play a role in the regulation of viral transcription and translation. For example, the long 3′-UTR of the M protein of measles virus and the 5′-UTR of the F protein of canine distemper virus play important roles in replication and virulence of these viruses ,.
Genome size in all NDV strains follows the ‘rule of six’ , which plays an important role in the replication of paramyxoviruses. Classical NDV strains which isolated before 1990s have a genome size of 15,186 bp ,; however, NDV strains isolated recently from geese have an extra region of six highly-conserved nucleotides (TCCCAC) located in the 5′-UTR of the NP gene after 1990s in China. This sequence has been found in NDV strains such as ZJ-1 (GenBank accession no. AF431744), WF00G (GenBank accession no. FJ754273), JSG0210 (GenBank accession no. JF340367), and NA-1, the strain used in this study .
The current study aimed to develop a rescue system that would allow genetic modification of the virulent NA-1 strain. The system was then used to examine the effects of mutations within the 5′-UTR of the NP gene. This rescue system, based on the intracellular synthesis of T7 RNA polymerase transcripts corresponding to the NDV genome, not only serves as a technical platform for the investigation of NA-1 gene function, but may also provide an opportunity to develop a live-attenuated vaccine candidate for geese.
Virus, plasmids, and cells
The genotype VII NDV virus designated NA-1 was isolated from geese in Jilin Province, Northeast China, in 1999. This virus has previously been sequenced  and the sequence deposited in GenBank (accession no. DQ659677). Purified viruses were propagated in 10-day-old embryonated specific-pathogen-free (SPF) eggs. The helper plasmids pCI-NP, pCI-P, and pCI-L, expressing NP, P, and L, respectively, were constructed and verified as to functionality. The pBR322 plasmid, containing T7 terminator and HdvRz sequence, was constructed by our lab. Baby hamster kidney (BHK-21) cells , expressing the T7 RNA polymerase, and DF-1 chicken fibroblast cells were grown in Dulbecco’s modified Eagle’s medium (DMEM) containing 10% fetal calf serum.
All birds were housed in individual ventilated cages supplied with filtered air under a 12-h light/dark cycle. SPF hygiene status was approved by Harbin Veterinary Research Institute. Food and water were available ad libitum.
Construction of a full-length NA-1 cDNA clone
RNA from NA-1 was extracted from allantoic fluid using TRIzol (Invitrogen) according to the manufacturer’s instructions. cDNA was prepared using a Superscript II reverse transcriptase kit (Invitrogen) and random hexamers, as per the manufacturer’s instructions. Six pairs of specific primers (sequences available on request) were designed to amplify the complete genomic cDNA complement. Resulting PCR amplicons, named A, B, C, D, E, and F, were purified and cloned into pMD18-T vector (TaKaRa). A full-length cDNA copy of the NA-1 genome was constructed using a modified pBR322 cloning vector (Figure 1), and designated pAF. Six segments of the NA-1 genome were amplified from pAF and sequenced to verify accurate assembly.
There is one RsrII site in the F protein segment of the NA-1 genome. To differentiate rescue virus from wild-type, the interior RsrII site was synonymously mutated from cggtccg to cggtctg using overlap PCR with primers 6FMP3 (5′-CCTCAGCCTCCATGTAAATCACaGACCGAATGACTGT-3′) and 6FMP4 (5′-CACATTGACACAGTCATTCGGTCtGTGATTTAC-3′). Another pair of primers, 6FP1 (5′-CGCGTCGACTCAGAGAACGAGAGGACAGAGATAC-3′) and 6FP2 (5′-CGCCGGTCCGACCAAACAGAGATTTGGTGAATG-3′), were designed to amplify the F segment of the NA-1 genome in a two-step process. The first part of the F segment was amplified using 6FP1 and 6FMP3 at 54°C. Then, primers 6FMP4 and 6FP2 were used to amplify the second part at 55°C.
Construction of a full-length cDNA clone with 5′-UTR deletion
The 5′-UTR of the NP gene was located in PCR fragment A described above. Using the plasmid constructed from segment A and the pMD18-T vector (named T-A), the six additional nucleotides (nucleotides 1648–1653) in the 5′-UTR of the NP gene were deleted using PCR. Primer pair 6AMP1 (5′-TCCGCCCAAAACCCACTCCCCGACCC-3′) and 6AMP2 (5′-GGGTCGGGGAGTGGGTTTTGGGCGGA-3′) was designed to amplify T-A and delete the six nucleotides. The resulting PCR product was digested using DpnI. Plasmids were transformed into competent cells DH5α, and clones containing the 6-bp deletion (named T-A(−)) were screened by PCR. To construct a full-length cDNA clone, T-A(−) was used to replace the same segment in pAF using SalI and AgeI. The resulting product was named pAF(−). Six segments of rNA-1 were amplified from pAF(−) and sequenced.
Recovery and identification of infectious NDV from cDNA
BHK-21 cells were grown overnight to 70–80% confluence in six-well culture dishes. These were cotransfected with 0.25 μg of pAF or pAF(−), 0.4 μg of pCI-NP, 0.2 μg of pCI-P, and 0.2 μg of pCI-L, using 8 μL of X-tremeGENE HP DNA Transfection Reagent (Roche). pAF, pAF(−), pCI-NP, pCI-P, and pCI-L were transfected into BHK-21 cells using liposome method separately as controls. Following incubation for 6 h, the culture medium was replaced with complete DMEM medium containing 5% allantoic fluid. After 5 days, the supernatant was harvested, passed through a 0.22 μm filter, and inoculated into the allantoic cavities of 10-day-old embryonated SPF eggs. Allantoic fluid was harvested after 96 h and hemagglutination assays (HA) were performed according to the World Organization for Animal Health’s standard method . Viral RNA was extracted from allantoic fluid and the F segment was amplified, cloned, and digested with RsrII (as described above) for identification.
To generate a growth curve, replication of the parental and rescued viruses was examined in DF-1 cells infected at a multiplicity of infection of 10 under one-step growth conditions. After adsorption for 1 h, the DF-1 cell monolayer was incubated in DMEM supplemented with 5% fetal calf serum at 37°C and 5% CO2. At 8-h intervals, the supernatant was collected to determine viral titers in a 50% tissue culture infection dose (TCID50) per milliliter by end-point titration. All titrations were performed in triplicate.
Genetic stability and virulence
The rescued virus was passaged in 9-day-old SPF eggs, obtained from the Harbin Veterinary Research Institute, to assess the genetic stability and pathogenicity of the rescued virus in vivo. After 3 days, the allantoic fluids were examined by HA assay  to confirm effective viral replication. The virus, at a 29 HA titer, was diluted by 0.001% every passage and the inoculated into the allantoic cavity of three 10-day-old embryonated SPF eggs at 37°C. Each virus generation was identified using reverse-transcription polymerase chain reaction (RT-PCR), and the 56–1789 nucleotide region of the NP gene was sequenced after each passage of NA-1. After 10 passages, the mean time until death (MDT) and intracerebral pathogenicity index (ICPI) were determined in embryonated SPF eggs or in SPF chickens .
Quantitative real-time PCR and western blot
The total RNA of each of the virus strains was extracted using TRIZOL reagent (Invitrogen), and cDNA was subsequently synthesized using a random hexamers primer, and then standardized to 1000 ng of total RNA per sample. To assess the transcription level of the NP gene, gene-specific primers were designed. The primers used for these mRNA analyses were: 5′-TGTACCACCAACTGCTTG-3′ and 5′-ATCACGCCACAGTTTCC-3′ for GAPDH, 5′-GGCAGCATGGACATTCCT-3′ and 5′- CCTTGCGGCTTGTTTGAT-3′ for NP. Each real-time PCR reaction (20 μL volume) contained 10 μL SYBR Green real-time PCR master mix (TAKARA), 0.25 μM gene-specific primers, and 1 μL standardized template cDNA. Amplification conditions were as follows: 95°C for 5 min, followed by 35 cycles of 95°C for 30 s, 56°C for 30 s, and 72°C for 30 s. After a final extension at 72°C for 10 min, a melting-curve analysis was performed to ensure specificity of the products. NP gene transcription levels were normalized to that of GAPDH, and transcription relative to the sample with the lowest expression was calculated using the 2-ΔΔCT method .
To assess the effects of deletions in the UTR of the NP gene on translation, NP protein expression by the NDV mutant was analyzed by western blotting . DF-1 cells were infected at 5 PFU/cell, and total proteins were collected at 24 h post-infection. Cells were washed with PBS, scraped, collected by low-speed centrifugation, and lysed with lysis buffer (6.25 mM Tris [pH 6.8], 1% SDS, 10% glycerol, 6 M urea, 0.01% bromophenol blue, 0.01% phenol red) for 30 min on ice. The lysates were clarified by centrifugation at 7,000 × g for 10 min and used for western blot analysis. Samples of the cell lysate were mixed with Laemmli sample buffer and boiled for 3 min prior to electrophoresis. The boiled samples were separated by SDS-polyacrylamide gel electrophoresis on a 12% gel, and the resolved proteins were transferred to nitrocellulose membrane and incubated using a primary monoclonal antibody against the NDV NP protein, or GADPH protein as a loading control. The membranes were washed with PBS-T, and then incubated with the second horseradish peroxidase-conjugated goat anti-rabbit IgG antibody. Bands were visualized using SuperSignal West Pico Substrate (Thermo Scientific). This experiment was performed in triplicate.
Experiment in geese
One-day-old Huoyan geese, obtained from in Jilin Province, Northeast China, were negative for both NDV antigen and maternal antibody. Thirty 1-day-old geese were divided into three equal groups. Purified NA-1 virus, with a hemagglutination assay (HA) titer of 29, was diluted 1/10 in sterile isotonic saline, and 0.05 mL of the diluted virus was injected intracerebrally into each of 10 geese in group 1. Groups 2 and 3 were infected with rNA-1 (HA titer 29) and rNA-1(−) (HA titer 29), respectively, using the same method. The birds were observed for clinical symptoms and mortality daily for 8 days. At each time point the birds were scored as follows: 0, no clinical symptoms; 1, showing signs of illness, and 2, deceased. The intracerebral pathogenicity index (ICPI) results represent the mean score per bird per observation over the 8-day period.
All values are the mean ± SEM of at least three individual samples for the in vitro experiments. We used analysis of variance, as calculated by GraphPad Prism 5, to establish the statistical significance of differences between the experimental groups. Individual inter-group comparisons were made using a two-tailed unpaired t-test with Welch’s correction. Differences were considered significant at P < 0.05 (*).
Animals were treated humanely and with regard for alleviation of suffering. All animal experiments were carried out according to the experimental practices and standards approved by the Animal Welfare and Research Ethics Committee at Jilin University (Approval ID: 20130713–1).
Identification of full-length cDNA
To determine whether the cDNA segments combined in our modified pBR322 vector (pAF from NA-1 and pAF(−) from NA-1 minus the additional nucleotides in the NP 5′-UTR) were full-length copies of the viral genome, six segments of pAF and pAF(−) were amplified and sequenced. Sequencing showed that pAF and pAF(−) were full-length cDNA copies of the parental virus. pAF(−) was identical to pAF except for the 6-bp deletion in the 5′-UTR of the NP gene at nucleotides 1648–1653.
Characterization of the rescued virus
To generate infectious NDV from the cDNA of pAF and pAF(−), plasmids were cotransfected with BHK-21 cells and incubated for 5 days. On day 6, the supernatants from the experimental and control groups were inoculated into the allantoic cavities of SPF 10-day-old embryonated eggs. Allantoic fluid was harvested 96 h later and virus identified via HA assay. The viruses rescued from pAF and pAF(−) were named rNA-1 and rNA-1(−), respectively, and had identical genome sequences apart from the 6-bp deletion in the 5′-UTR of the NP gene. The rescued rNA-1 and rNA-1(−) viruses had HA titers of 24 and 23, respectively.
F segments of rNA-1 and NA-1 that were amplified by 6FP1 and 6FP2, respectively, and digested with RsrII showed that rNA-1 can be differentiated from parental virus NA-1 by restriction fragment length polymorphism (Figure 2).
Transcription and translation of NP
Quantitative real-time PCR analysis showed that expression of NP was significantly lower in rNA-1(−) than in rNA-1 and NA-1 (P < 0.05), and that there was no significant difference in expression levels between NA-1 and rNA-1 (P > 0.05) (Figure 3). Western blot analysis confirmed that there was no significant difference in the ratio of NP to GAPDH between NA-1 and rNA-1 (P > 0.05), but that the ratio in rNA-1(−) was significantly lower than that of rNA-1 or NA-1 (P < 0.05) (Figure 4). This result indicates that the 6-bp deletion in the 5′-UTR of the NP gene at nucleotides 1648–1653 decreases the expression of the NP gene.
Growth and virulence analysis
Growth curves showed that rNA-1(−) replicated more slowly than rNA-1, but that the two strains have similar growth characteristics (Figure 5). Sequencing and alignment of rNA-1 and rNA-1(−) during serial passages in 10-day-old embryonated SPF eggs showed no genetic drift during the 10 separate passages. After the final passage, two standard pathogenicity tests, MDT and ICPI, were performed to compare the virulence of rNA-1, rNA-1(−), and NA-1. The MDT and ICPI of wild-type NA-1 have previously been determined to be 59.6 h and 1.65, respectively. It is generally agreed that NDVs with ICPI values greater than 1.5 are categorized as velogenic strains. The MDT and ICPI values for rNA-1 from the 10th passage were 59 h and 1.62, respectively. The MDT and ICPI values for rNA-1(−) from the 10th passage were 62 h and 1.52, respectively (Table 1). These data demonstrated that strains rNA-1 and rNA-1(−) remained velogenic and genetically stable during serial passages in vivo. Although the virulence of rNA-1(−) did not change significantly in chickens, it was decreased compared with NA-1 and rNA-1.
Experiment in geese
The pathogenicity of the viruses in geese was examined, and the results were calculated by the ICPI method. At 10 days post-infection, all animals from the three groups were dead. The ICPI results for NA-1 and rNA-1 were 1.60, while the value for rNA-1(−) was 1.35. This indicated that rNA-1(−) was less pathogenic than rNA-1 in geese.
Highly pathogenic NDV causes a rapidly spreading illness in poultry, with four worldwide pandemics since its emergence in the early 20th century. The most recent pandemic was caused by genotypes VII and VIII in the late 1980s and 1990s ,. For decades, this virus was not considered to be infectious in mammals and waterfowl. However, cases of waterfowl infection have been frequent in China since the late 1990s. According to recent reports, virulent strains isolated from geese were responsible for outbreaks of ND in waterfowl and belonged to genotype VII ,. It is still unclear whether NDVs circulated amongst chickens have the potential to become highly pathogenic in waterfowl, or whether the specificity of receptors could influence the host range of the virus and viral virulence.
Virulent NDV strains are still isolated frequently from vaccinated birds, demonstrating that NDV remains a sustained threat to commercial flocks. Nowadays, the prevailing NDV strain of geese in China is genotype VII, while vaccine strains such as B1, Clone30, La Sota, and V4, which are used to protect geese from NDV infection, belong to different genotypes. The antigenic differences between the prevailing and vaccine strains might explain current ND outbreaks in vaccinated poultry flocks ,,.
Methods for controlling and preventing ND may be the same as those used for avian influenza virus, with a significant need for epidemiological surveillance. The best control method is the development of homologous vaccines with high efficiency and safety. However, no current commercial vaccine is available for geese for protection against virulent genotype VII NDV. Therefore, we set out to establish a reverse genetics system that allows the genetic manipulation of a virulent strain, which could be used as a vaccine candidate or a viral vector to express foreign genes in our ongoing studies.
The amino acid sequence at the F protein cleavage site is a major determinant of NDV virulence ,, although other genes also have an effect on virulence . Transcriptional and translational control signals may also modulate virulence by controlling protein expression, as has been described for canine distemper virus and NDV ,,. A Typically ICPI is only performed in one day old SPF chickens. However, due to tropism differences that has been observed with different NDV strains in different poultry species, that an ICPI was warranted in geese. That the geese were not SPF due to lack of availability, but that they were NDV antibody negative. The rNA-1 and rNA-1(−) viruses had HA titers of 24 and 23, respectively, in the current study. A correlation between virulence and the efficiency of viral replication has been observed for many viruses . The lower titer of the rNA-1(−) strain may be the result of the decreased levels of transcription and expression of the NP gene following deletion of the six nucleotides from the 5′-UTR. This might also account for the low ICPI value in geese. We must note that the contribution of genes to virulence may be dependent on the particular virus strain used. Additionally, because of the selective pressure from the immune system, the virus can evolve by itself. Therefore, it is not the best way to modify the HN and F genes of velogenic NDV for obtaining live-attenuated vaccine; the greater the pressure, the greater the resistance.
The results of this study demonstrate that the virulence of the rNA-1(−) strain is still at the edge of the velogenic pathotype, and is not significantly different from that of the rNA-1 strain. However, the pathogenicity of rNA-1(−) towards geese is lower than that of rNA-1. Therefore, the deletion of the six nucleotides from the 5′-UTR of the NP gene could contribute to the decrease in virulence of the mutant NA-1 strain.
Fauquet CM, Fargette D: International committee on taxonomy of viruses and the 3,142 unassigned species. Virol J. 2005, 2: 64-10.1186/1743-422X-2-64.
Mayo MA: Virus taxonomy - Houston 2002. Arch Virol. 2002, 147 (5): 1071-1076. 10.1007/s007050200036.
Miller PJ, Decanini EL, Afonso CL: Newcastle disease: evolution of genotypes and the related diagnostic challenges. Infect Genet Evol. 2010, 10 (1): 26-35. 10.1016/j.meegid.2009.09.012.
Chambers P, Millar NS, Bingham RW, Emmerson PT: Molecular cloning of complementary DNA to Newcastle disease virus, and nucleotide sequence analysis of the junction between the genes encoding the haemagglutinin-neuraminidase and the large protein. J Gen Virol. 1986, 67 (Pt 3): 475-486. 10.1099/0022-1317-67-3-475.
Steward M, Vipond IB, Millar NS, Emmerson PT: RNA editing in Newcastle disease virus. J Gen Virol. 1993, 74 (Pt 12): 2539-2547. 10.1099/0022-1317-74-12-2539.
Yusoff K, Tan WS: Newcastle disease virus: macromolecules and opportunities. Avian Pathol. 2001, 30 (5): 439-455. 10.1080/03079450120078626.
Pedersen JC, Senne DA, Woolcock PR, Kinde H, King DJ, Wise MG, Panigrahy B, Seal BS: Phylogenetic relationships among virulent Newcastle disease virus isolates from the 2002–2003 outbreak in California and other recent outbreaks in North America. J Clin Microbiol. 2004, 42 (5): 2329-2334. 10.1128/JCM.42.5.2329-2334.2004.
Kim SH, Nayak S, Paldurai A, Nayak B, Samuel A, Aplogan GL, Awoume KA, Webby RJ, Ducatez MF, Collins PL, Samal SK: Complete genome sequence of a novel Newcastle disease virus strain isolated from a chicken in West Africa. J Virol. 2012, 86 (20): 11394-11395. 10.1128/JVI.01922-12.
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 (3): 239-256. 10.1080/030794503100009783.
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 (7): 1387-1403.
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 (1–2): 206-211. 10.1016/j.jviromet.2006.10.012.
Wang X, Gong Z, Zhao L, Wang J, Sun G, Liu Y, Tao P, Zhang H, Li S, Jiang F, Jiang F, Hu Y, Zhang X: Complete genome sequences of newcastle disease virus strains isolated from three different poultry species in china. Genome Announc 2013, 1(4).
Rauw F, Gardin Y, Palya V, Anbari S, Lemaire S, Boschmans M, van den Berg T, Lambrecht B: Improved vaccination against Newcastle disease by an in ovo recombinant HVT-ND combined with an adjuvanted live vaccine at day-old. Vaccine. 2010, 28 (3): 823-833. 10.1016/j.vaccine.2009.10.049.
Kapczynski DR, King DJ: Protection of chickens against overt clinical disease and determination of viral shedding following vaccination with commercially available Newcastle disease virus vaccines upon challenge with highly virulent virus from the California 2002 exotic Newcastle disease outbreak. Vaccine. 2005, 23 (26): 3424-3433. 10.1016/j.vaccine.2005.01.140.
Miller PJ, King DJ, Afonso CL, Suarez DL: Antigenic differences among Newcastle disease virus strains of different genotypes used in vaccine formulation affect viral shedding after a virulent challenge. Vaccine. 2007, 25 (41): 7238-7246. 10.1016/j.vaccine.2007.07.017.
Anderson DE, von Messling V: Region between the canine distemper virus M and F genes modulates virulence by controlling fusion protein expression. J Virol. 2008, 82 (21): 10510-10518. 10.1128/JVI.01419-08.
Takeda M, Ohno S, Seki F, Nakatsu Y, Tahara M, Yanagi Y: Long untranslated regions of the measles virus M and F genes control virus replication and cytopathogenicity. J Virol. 2005, 79 (22): 14346-14354. 10.1128/JVI.79.22.14346-14354.2005.
Peeters BP, Gruijthuijsen YK, de Leeuw OS, Gielkens AL: Genome replication of Newcastle disease virus: involvement of the rule-of-six. Arch Virol. 2000, 145 (9): 1829-1845. 10.1007/s007050070059.
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 (Pt 1): 131-136.
Krishnamurthy S, Samal SK: Nucleotide sequences of the trailer, nucleocapsid protein gene and intergenic regions of Newcastle disease virus strain Beaudette C and completion of the entire genome sequence. J Gen Virol. 1998, 79 (Pt 10): 2419-2424.
Xu M, Chang S, Ding Z, Gao HW, Wan JY, Liu WS, Liu LN, Gao Y, Xu J: Genomic analysis of Newcastle disease virus strain NA-1 isolated from geese in China. Arch Virol. 2008, 153 (7): 1281-1289. 10.1007/s00705-008-0115-9.
Buchholz UJ, Finke S, Conzelmann KK: Generation of bovine respiratory syncytial virus (BRSV) from cDNA: BRSV NS2 is not essential for virus replication in tissue culture, and the human RSV leader region acts as a functional BRSV genome promoter. J Virol. 1999, 73 (1): 251-259.
Hanson RP, Brandly CA: Identification of vaccine strains of Newcastle disease virus. Science. 1955, 122 (3160): 156-157.
Hierholzer JC, Suggs MT, Hall EC: Standardized viral hemagglutination and hemagglutination-inhibition tests. II. Description and statistical evaluation. Appl Microbiol. 1969, 18 (5): 824-833.
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 (8): 4176-4184. 10.1128/JVI.78.8.4176-4184.2004.
Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.
Huang Z, Elankumaran S, Yunus AS, Samal SK: A recombinant Newcastle disease virus (NDV) expressing VP2 protein of infectious bursal disease virus (IBDV) protects against NDV and IBDV. J Virol. 2004, 78 (18): 10054-10063. 10.1128/JVI.78.18.10054-10063.2004.
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 (1): 49-64. 10.1007/s007050050267.
Tsai HJ, 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 (1–2): 19-30. 10.1016/j.vetmic.2004.09.005.
Ke GM, Yu SW, Ho CH, Chu PY, Ke LY, Lin KH, Tsai YC, Liu HJ, Lin MY: Characterization of newly emerging Newcastle disease viruses isolated during 2002–2008 in Taiwan. Virus Res. 2010, 147 (2): 247-257. 10.1016/j.virusres.2009.11.008.
Mase M, Inoue T, Imada T: Genotyping of Newcastle disease viruses isolated from 2001 to 2007 in Japan. J Vet Med Sci. 2009, 71 (8): 1101-1104. 10.1292/jvms.71.1101.
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 (3): 193-203. 10.1016/S0378-1135(02)00050-0.
Lin MY, Liu HJ, Ke GM: Genetic and antigenic analysis of Newcastle disease viruses from recent outbreaks in Taiwan. Avian Pathol. 2003, 32 (4): 345-350. 10.1080/0307945031000121086.
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.
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.
Dortmans JC, Koch G, Rottier PJ, Peeters BP: Virulence of Newcastle disease virus: what is known so far?. Vet Res. 2011, 42: 122-10.1186/1297-9716-42-122.
Kim SH, Samal SK: Role of untranslated regions in regulation of gene expression, replication, and pathogenicity of Newcastle disease virus expressing green fluorescent protein. J Virol. 2010, 84 (5): 2629-2634. 10.1128/JVI.02049-09.
Yan Y, Rout SN, Kim SH, Samal SK: Role of untranslated regions of the hemagglutinin-neuraminidase gene in replication and pathogenicity of newcastle disease virus. J Virol. 2009, 83 (11): 5943-5946. 10.1128/JVI.00188-09.
The research was supported by the National Natural Scientific Fund (31272561) and the Chinese Special Fund for Agro-scientific Research in the Public Interest (201303033).
The authors declare that they have no competing interests.
JL was involved in experimental design, sampling, sample processing, data analysis, and drafting of manuscript. YC, RY, and YQ were involved in supervision, research conception and design, data analysis, interpretation, and manuscript revision. CD and XL were involved in experimental design, data interpretation, and manuscript revision. ZD and CW were involved in supervision of the experiments and manuscript revision. All authors read and approved the final manuscript.
Jiaxu Liu, Yanlong Cong contributed equally to this work.