Genomic and Pathogenic Characteristics of Virulent Newcastle Disease Virus Isolated from Chicken in Live Bird Markets and Backyard Flocks in Kenya
Newcastle disease (ND) causes significant economic losses in the poultry industry in developing countries. In Kenya, despite rampant annual ND outbreaks, implementation of control strategies is hampered by a lack of adequate knowledge on the circulating and outbreak causing-NDV strains. This study reports the first complete genome sequences of NDV from backyard chicken in Kenya. The results showed that all three isolates are virulent, as assessed by the mean death time (MDT) and intracerebral pathogenicity index (ICPI) in specific antibody negative (SAN) embryonated eggs and 10-day-old chickens, respectively. Also, the polybasic amino acid sequence at the fusion-protein cleavage site had the motif 112RRQKRFV118. Histopathological findings in four-week-old SPF chicken challenged with the NDV isolates KE001, KE0811, and KE0698 showed multiple organ involvement at five days after infection with severe effects seen in lymphoid tissues and blood vessels. Analysis of genome sequences obtained from the three isolates showed that they were 15192 base pair (bp) in length and had genomic features consistent with other NDV strains, the functional sites within the coding sequence being highly conserved in the sequence of the three isolates. Amino acid residues and substitutions in the structural proteins of the three isolates were similar to the newly isolated Tanzanian NDV strain (Mbeya/MT15). A similarity matrix showed a high similarity of the isolates to NDV strains of class II genotype V (89–90%) and subgenotype Vd (95–97%). Phylogenetic analysis confirmed that the three isolates are closely related to NDV genotype V strains but form a distinct cluster together with NDV strains from the East African countries of Uganda and Tanzania to form the newly characterized subgenotype Vd. Our study provides the first description of the genomic and pathological characteristics of NDV of subgenotype Vd and lays a baseline in understanding the evolutionary dynamics of NDV and, in particular, Genotype V. This information will be useful in the development of specific markers for detection of viruses of genotype V and generation of genotype matched vaccines.
Newcastle disease (ND) is a highly contagious disease of poultry that is notifiable to the World Organization for Animal Health (OIE) . In developing countries, ND has a significant economic impact on the poultry sector due to the high morbidity and mortality . Newcastle disease is caused by Newcastle disease virus (NDV), a single-stranded negative-sense RNA virus classified in the genus Avulavirus within the family Paramyxoviridae in the Order Mononegavirales . Newcastle disease virus consists of a 15 kb genome with six open reading frames (ORF) encoding six major structural proteins, namely, nucleoprotein (NP), phosphoprotein (P), matrix protein (M), fusion-protein (F), hemagglutinin neuraminidase (HN), and the RNA-dependent RNA polymerase (L) in the order, 3′-NP-P-M-F-HN-L-5′ [4, 5]. Each of these genes is flanked by short extragenic “leader” and “trailer” sequences. Besides, each gene starts with a conserved gene start (GS) and ends with a conserved sequence: gene end (GE) .
Newcastle disease virus, although belonging to one serotype, is highly diverse genetically and antigenically. NDV isolates from the world are grouped into two distinct classes (class I or II), based on genome lengths and nucleotide sequences. Class I NDV is less divergent and belongs to a single genotype, whereas class II isolates are more diverse . Based on the analysis of the complete coding region of the fusion gene, class II NDV has been divided into 18 genotypes [7–10]. Newcastle disease virus belonging to class II, which was isolated between the 1930s and 1960s, includes genotypes I, II, III, and IV. These are considered “ancient” and have genome sizes of 15,186 nucleotides (nt) . The currently circulating or “recent” class II genotypes were isolated after the 1960s (V, VI, VII, VIII, and X–XVIII). The latter genotypes have an insertion of six nucleotides into the 5′ noncoding region of the nucleoprotein gene. They, therefore, possess 15192 nt in their genomes . More genetic distance has been observed between “ancient” and the “recent” genotypes . Some of the genotypes are further divided into subgenotypes according to set criteria .
The most recent genotypes (XIV–XVIII) have mostly been isolated on the African continent [8, 10, 12]. A study on isolates from Southern Africa as well as from Namibia, Zambia, Zimbabwe, and Mozambique revealed the presence of genotypes VIII and VII [13, 14]. In Zambia, genotype XIII has also been isolated from poultry . A comparative analysis of APMV-1 sequences from West Africa indicated that the viruses correspond to genotypes I and VI and novel genotype XVIII . Novel genotypes XIV and XVII have also been isolated from Western and Central African countries [17, 18]. Isolates from Ethiopia were assigned to genotypes VII and VI [19–21]. Isolates from Tanzania were assigned to genotypes I, V, and Via, respectively . Previous studies reported that isolates from Sudan belonged to genotypes VII and VI [23, 24]. A study in Eastern Uganda in 2001 revealed that they belonged to genotype VIa , 2004) . A more recent study in Uganda isolated strains of genotype Vd from poultry in LBMs . It is suggested that increased surveillance on the African Continent may lead to the discovery of more genotypes and subgenotypes . This is thought to be a result of the predominance of the backyard poultry management system which encourages the presence and spread of varied strains of velogenic NDV which threatens commercial poultry enterprises .
In Kenya, like other developing countries, the high population of poultry is kept on smallholder backyard farms. Poultry production is important for the livelihoods of many rural households as a source of protein and income and fulfills the sociocultural role. Although Newcastle disease is endemic in Kenya, potential reservoirs of NDV virus exist which are responsible for the annual epidemic outbreaks of ND causing annual epidemic outbreaks that result in high bird mortality, economic losses, and negative impact on the livelihoods of rural populations . The introduction of birds from informal live bird markets (LBMs) poses an important source of introduction of NDV into a flock . Live bird markets may enhance the movement of infected birds and the spread of NDV is rampant in Kenya and may be important in the epidemiology of ND . However, until now, only a few studies have reported the genetic and pathological characteristics of NDV circulating in Kenya. This study for the first time characterized and compared the genomic and pathological characteristics of NDV from live bird markets and backyard poultry farms in Kenya to establish the epidemiological link and role of live bird markets and backyard poultry in the spread of NDV in Kenya.
2. Materials and Methods
Cloacal swabs and tissue samples for the study came from chicken showing clinical symptoms of Newcastle disease. These samples were obtained from live bird markets and backyard poultry farms in Kenya during a survey of Newcastle disease virus described by . We inoculated samples from individual birds in three 10-day-old specific antibody negative (SAN) embryonated chicken eggs through the allantoic cavity. Three days after inoculation, the allantoic fluid was harvested and clarified by centrifugation at 4000 × g for 20 minutes at 4°C. We used the supernatant to run the hemagglutination inhibition test using specific antisera to the reference strains of NDV. Three NDV isolates designated KE001/2015 (Makueni), KE0811/2016 (Mombasa), and KE0698/2016 (Malaba) showed high HA titer (>24). The allantoic fluid of the three isolates was aliquoted for use in pathogenicity tests and genome sequencing and characterization.
2.2. Study Clearance and Permits of Compliance
The Kenya Agricultural and Livestock Research Organization (KALRO), Veterinary Science Research Institute- (VSRI-) Institutional Animal Care and Use Committee (IACUC), approved poultry handling procedures used in the study (protocol number: KALRO/VSRI/IACUC/1/23082017) according to the appropriate animal guidelines . We carried out the study in the biosafety level II isolation facility at KALRO-VSRI.
2.3. Biological Characterization of the NDV Isolates
The pathogenic potential of the three NDV isolates was evaluated using a standard assay to determine mean death time (MDT) in 10-day-old specific antibody negative (SAN) embryonated eggs and the intracerebral pathogenicity index (ICPI) in 1-day-old chicks . Briefly, ICPI was performed in 1-day-old chicks by inoculation of 50 μL of allantoic fluid diluted 10-fold in PBS without antibiotics as recommended . We calculated the intracerebral pathogenicity index as the mean of each observation for each bird. We observed the birds every 24 hours for 1 week. At each observation, the birds were scored 0 if normal, 1 if sick, and 2 if dead. We humanely sacrificed birds that were too ill and unable to eat by cervical dislocation and scored them as dead in the next observation. We gave dead birds a score of two at each of the remaining observations.
We considered Newcastle disease virus strains with ICPI score of 0.7 lentogenic, while those with a score of more than 1.5 as velogenic. On the other hand, those strains with ICPI scores between 0.7 and 1.5 were mesogenic. For MDT, allantoic fluid with virus was diluted in a series of 10-fold dilutions with phosphate-buffered saline (PBS) and each dilution inoculated into the allantoic cavity of specific pathogen-free 10-day-old embryonated eggs and incubated at 37°C. We observed the eggs every 8 hours for 7 days and the time of embryo death was recorded. The MDT was recorded as the mean time (hours) for the minimum lethal dose of the virus to kill all the embryos . We considered isolates with MDT ≤60 hours velogenic: MDT between 61 and 90 hours was mesogenic, while those with MDT >90 hours were lentogenic.
2.4. Pathological Characteristics of the NDV Isolates
To determine the pathological characteristics of the NDV isolates, we randomly allocated forty (40) SAN chicken aged four (4) weeks into four experimental groups consisting of three NDV infected groups (group 1-KE001/2015, n = 10; group 2-KE0698/2016, n = 10; and group 3-KE0811/2016, n = 10) and one control group (n = 10). We challenged birds in the NDV infected groups via intraocular inoculation with respective virus suspension (0.1 ml) containing 106 50% lethal dose (LD50). We inoculated the control group with 0.1 ml PBS and observed all birds daily for clinical signs of disease. We euthanized and necropsied two birds from each group at day 5 after inoculation and all birds with severe clinical symptoms. We monitored the birds daily until the end of the experiment, two weeks after challenge. Thereafter, we euthanized and necropsied all remaining birds. At necropsy, we collected tissues of trachea, spleen, brain, lung, liver, bursa of Fabricius, proventriculus, and caecal tonsils and fixed them in 10% buffered formalin. The samples were embedded into the paraffin and sectioned at 3 μm and routinely deparaffinized and stained with hematoxylin and eosin for histopathological examination.
2.5. RNA Extraction and Whole-Genome Sequencing
We extracted viral genomic RNA from the allantoic fluid using Trizol reagent (Life Technologies, USA) according to the manufacturer’s instructions. The concentration and purity of extracted RNA were measured using a spectrophotometer (ND-1000, Nanodrop Technologies, USA) and Qubit Fluorimeter (Invitrogen, Carlsbad, CA, USA), and the integrity of RNA was visualized by electrophoresis in a 1.2% formaldehyde agarose gel stained with GelRed and also using the Bioanalyzer. RNA extracts were selected for library preparation when the 260/280 purity index was ≥2.0 and the integral RNA in electrophoresis and Bioanalyzer measurements was greater than eight (RIN > 8). The concentration and purity of the extracted RNA were determined using a Qubit Fluorimeter (Thermo Scientific, Walton, Massachusetts, USA) before storage at −80°C until use.
The nucleotide sequence of the full viral genome was determined using next-generation sequencing (NGS). The cDNA libraries were prepared using a TruSeq RNA sample preparation kit (Illumina, CA, USA). Briefly, mRNA was prepared by purification and fragmentation of 2 μg of total RNA using oligo (dT) magnetic beads and was used as a template for cDNA synthesis by random hexamer priming. Paired-end (2 × 100 reads) sequencing was done using Miseq (Illumina). The reads of the full-length viral genome were assembled de novo with SPAdes assembler version 3.10.1 . We confirmed the NGS results with conventional Reverse-Transcription polymerase chain reaction (RT-PCR) using virus-specific primers. This involved purification of amplified products using Qiagen PCR Purification Kit and sequencing purified products by the Sanger method at Macrogen, Korea. The sequences obtained by the Sanger method were 100% identical to the NGS sequences.
2.6. Phylogenetic Analysis
We obtained complete sequences of reference strains for each NDV genotype from GenBank including the recent Tanzania isolate MK583011 (Supplementary Table S1). We also compared the study isolates with other NDV isolates from Africa, which are present in the GenBank. These sequences from GenBank and the complete coding sequences of the four NDV isolates were aligned using MUSCLE v. 3.8.31 [33, 34]. Phylogenetic and molecular evolutionary analyses were conducted using MEGA (Molecular Evolutionary Genetics Analysis) version 6.0 . We constructed phylogenetic trees using the Maximum Likelihood (ML) method and estimated the tree using best-fit general time-reversible (GTR) model of nucleotide substitution with gamma-distributed rate variation among sites. We employed a bootstrap resampling process (1000 replications) to assess the robustness of individual nodes of phylogeny. Evolutionary distances were computed by the Pairwise Distance method using the Maximum Composite Likelihood method .
3. Results and Discussion
3.1. Biological Characterization
The three isolated viruses were virulent by the OIE standard criteria. Table 1 shows the MDT in embryonated eggs and ICPI values for the three NDV isolates as well as the amino acid sequence at the cleavage site. The ICPI for the three isolates was >1.5, while the MDT was <60 hours indicative of velogenic strains. Further, the amino acid sequence motif at the cleavage site of the fusion gene concurred with the ICPI and MDT scores of the three isolates that indicated their velogenic nature.
3.2. Pathological Lesions
All chickens inoculated with NDV study isolates were successfully infected and showed signs of depression, mucoid diarrhea, leg paralysis, comb cyanosis, and head edema and facial swelling with mucoid exudate from the nostrils. Gross lesions consisted of enlargement of organs such as the liver, spleen, and kidneys in the majority of the birds. Gross lesions were also evident in the mucosa of the proventriculus, spleen, liver, and small intestines (Figure 1).
These lesions included hemorrhagic and necrotic patches and blood vessel congestion. The spleen had numerous multifocal whitish spots that were evident on the surface. The trachea also showed mucosal hemorrhages and necrosis with a few birds showing pneumonic consolidation in the lungs with the discharge of caseous exudates from the air sacs. The intestines showed darkened areas of necrosis that were evident on the serosa surface. These necrotic areas in the lumen of the intestines were covered with mucoid content underneath of which were ulcerative foci. The bursa of Fabricius showed marked edema (Figure 1).
All groups infected with the NDV study isolates presented similar histological changes in the multiple organs and tissues at 5 days after challenge. The most extensive damage involved blood vessels and lymphoid tissues. We observed vascular changes including hemorrhages, congestion, and edema in majority of the organs. Severe congestion of blood vessels and diffuse hemorrhages were present in the liver, small intestines, and kidneys.
We observed depletion of the lymphoid tissue in lymphoid organs such as the spleen, bursa, and payers patches. On the other hand, organs such as lungs, small intestines, and liver showed lymphocyte infiltration (Figure 2). The intestinal mucosal epithelium was denuded with loss of villi structure. The findings resemble pathological presentation in infection by highly virulent viruses .
3.3. Genome Similarity Analysis
Similarity search of the three isolates on the nucleotide databases through the BLAST algorithm yielded the highest similarity (97%) to the genome sequence of a previous Kenyan Newcastle disease virus (NDV) isolate, A48 (Accession JQ217420), and Tanzanian isolate, Tanzania/Mbeya/MT15/2012. The isolates also showed high similarity (92%) to the complete genome sequence Newcastle disease virus strain, Largo/71 from the USA classified as genotype V clade II.
The percentage similarity for nucleotide sequences encoding the different NDV genes was compared between isolate KE001/2015 (Makueni) and other NDV strains including the study strains: KE0811/2016 (Mombasa) and KE0698/2016 (Malaba) and others from the GenBank (NCBI). Study isolates were highly similar to NDV strains grouped under genotype V ranging from 89 to 97 (Table 2). On the other hand, the similarity of the study isolates to NDV strains of genotypes I and II ranged between 82 and 83 percent. The similarity of the coding regions of KE001/2015 (Makueni) to NDV vaccine strains (LaSota, BHG, and I-2) ranged from 77 to 85%, while the similarity to gene sequences of genotype V strains was the highest ranging from 87 to 98%. We observed the lowest similarity between the study and vaccine strains in the coding sequence of the V and P gene.
3.4. Genomic Features of Coding and Noncoding Regions of Isolates
The genomes of the three isolates had a length of 15192 nucleotides with the insertion of six nt (cytosine: “CCCCC (T) C”) in the 5′ UTR noncoding region of the NP gene between position 1643 nt and 1650 nt. The genome of the three isolates was organized as 3′-NP-P-M-F-HN-L-5′. The genomes started with a leader and trailer sequences at the 3′ and 5′ end with lengths of 330 55 nt and 114 nt, respectively (Table 3). Similarly, the isolates had a considerably conserved 10 nt sequence at the start of every gene (GS) and another one at the end (GE) of each gene. The open reading frame of each gene was flanked with 3′ and 5′ end untranslated regions (UTRs) of varying lengths on each end (Table 3). Between GE of one gene and GS of the next gene was a conserved sequence, intergenic sequence (IGS). The length of the nucleotide and amino acid sequence of each gene was varied with the RNA-dependent RNA polymerase (L) gene presenting the longest gene.
3.5. Genomic Features of Protein-Coding Regions of Isolates
In multiple sequence alignment, the NP protein of the study isolates showed the highest similarity to NDV strain A48/Kenya and Mt/Tanzania. The amino acid sequence of the three isolates had three conserved domains on the NP protein: region 1 spanning between aa 171 and 181 with the motif: 171QVWVTVAKAMT181; region 2 spanning between aa 267 and277 with the motif: 267FFLTLKYGINT277; and the third conserved region between aa 322 and 336 with the motif: 322FAPAEYAQLYSFAMG336. The latter region found in the NP ORF (Open reading frame) of all paramyxoviruses of the genera Avulavirus is a 15 amino acid region and was similarly conserved in all genotype representative strains of NDV used in this study. These regions are important for the formation of the NP-P complexes and regulation of viral replication [37–39].
The phosphoprotein (P) gene had the most variable sequence high variability between the three isolates and other NDV strains but had a conserved site region at nucleotide positions 394–401 (5′-AAAAAGGG-3′) corresponding to RNA editing site in the P gene . At this site, insertion of one G gives rise to the V protein, while the addition of two Gs gives rise to W protein. The V protein sequence of the study isolates aligned with that of reference genotypes of NDV. The length of the V protein is 239 aa for the NDV strains (Figure 3). The V protein forms into unique folds through seven cysteine amino acids at the carboxyl- (C-) terminus. For the three isolates, the cysteines are at positions 196, 200, 212, 214, 217, 221, and 224. The V protein of the study isolates was 96.5% similar to that of Tanzanian isolate: Mbeya/MT15/2012 (MK583011) and 90.8% similar to strain Largo/71.
The matrix protein of the three isolates had conserved sites among them: 40QYRIQRLDSWIDSKE54, 171YKVNFVSLTVVPRKD185, 213EVDPRSPLVKSLSRS227, 247KKGKKVTFDKLERKIRR263, and 349KIEKRHTIAKYNPFK363, which were similar in other NDV strains except for a few substitutions. These sites correspond to immunodominant epitopes of the matrix protein . For the M protein to be localized in the nucleus, it utilizes nuclear localization signals (NLS). For the study isolates, the NLS was 246DKKGKKVTFDKLERKIRR263.
In the fusion protein sequence of the three isolates, the 5 potential asparagine (N) linked glycosylation sites, 85, 191, 366, 447, and 471, as well as cysteine residue positions, 76,199, 338, 347, 362, 370, 394, 399, 401, and 424 and individual residues, D72, E74, A75, K78, A79, and L343, were highly conserved in the study isolates. Besides, the study strains and other genotypes shared similar neutralizing antibody epitope amino acid residues of the fusion gene at positions 151 to 171, 151ILRLKESIAATNEAVHEVTDG171. The conserved sites have been reported as important functional sites of the F gene [42, 43]. The key function of the F protein is to initiate the fusion of viral surface to the host cell membrane. The amino acid motif at the F gene cleavage site is a major determinant of virulence. For the study isolates, the cleavage site had the motif 112R-R-Q-K-R-F117 corresponding to the cleavage site of velogenic viruses.
The HN protein plays a critical role in virus infectivity by binding sialic acid-containing receptors. We analyzed the HN protein of study isolates to identify the amino acid residues essential for receptor recognition. Fourteen residues essential for receptor recognition, R174, I175, D198, K236, E258, Y299, Y317, E401, R416, R498, R516, Y526, and E547, in the HN protein of the study isolates were similar to that of other NDV reference strains. To further understand the structure of the HN protein of the study isolates, the seven antigenic sites within the HN protein involved in the formation of the three-dimensional HN molecule were compared across the NDV reference strains (Table 4).
Compared to the strains I-2, we observed substitutions at different antigenic sites in the HN protein of the study isolates including D569N at antigenic site 2, N263R at site 3, R321K and K333R at site 4, and D494N at site 12. For the study isolates, the sequence at the antigenic determinant linear epitope of HN at position 345–353 had the motif 345PDEQDYQVR353 with an I/V substitution at aa 352. However, there was a substitution E347K for one study isolate, KE0698 (MN685356). Like other NDV strains, the HN protein of study isolates had a conserved amino acid sequence 79DVIDRVYKQVALESPLALLSTESIIMNAITSLSYQIN115 within the dominant linear antigenic domain of the HN gene which spans aa 53 to 192 . This was similar in NDV strains except three substitutions at residues V81I, N98S, and T102I observed in the study isolates. The sialic acid binding site of the HN protein  of the study isolates, like other NDV strains, was located at aa positions 234–239. For the study isolates, the binding site had the motif 234NRKSCS239. Other amino acid residues that form HN antigen epitopes were similar in the study isolates and other NDV strains including K138, P244, F277, T380, T409, and H482. These sites correspond to the neuraminidase active site that binds the fusion protein .
The RNA-dependent RNA polymerase (L) gene of the study isolates and other NDV strains had a highly conserved region between aa 740 and 753 (740SHCRVACMVQGDNQ753) which corresponds to the putative active site of the L gene essential to polymerase function . Other conserved domains of the L gene were observed between aa 503 and 607 and 418 between aa 634 and 854. The L gene protein was highly conserved among the NDV structural 419 proteins in the study isolates.
3.6. Phylogenetic Analysis
Eighteen genotypes (I–XVIII) of NDV class II have been described [5, 17]. Phylogenetic analysis was based on the complete genome sequences including genome sequences of known strains isolated from African countries including South Africa, West Africa, and East Africa (Figure 4). The phylogenetic analysis of complete genome sequences revealed that the study isolates were clustered separately from the sequences from South Africa, which were classified as genotypes VII and XII and those of West Africa that were classified in genotypes XIV and XVIII. However, the study isolates clustered together with others from East Africa (Tanzanian and previous Kenyan strains) as well as other strains from the United States of America grouped in genotype V (Figure 4).
Within genotype V, the East African strains including study strains formed a closer cluster that has been referred to as new subgenotype Vd . The genotype V of NDV is separated into four subgenotypes (Va–d). With the earlier classified subgenotypes, Va, Vb, and Vc hav been isolated mainly from Europe, Central America, and South America [9, 48, 49]. To confirm the results obtained by genome analysis, we performed a phylogenetic analysis based on the complete F gene coding sequence, which incorporated more sequences representing the known subgenotypes of genotype V (Figure 4(a)).
Complete genome, protein level, and biological analysis of three NDV isolates from backyard chicken in Kenya demonstrated that the isolates have the biological and pathological characteristics of virulent viruses. Further, the study demonstrated that the isolates belong to genotype V and subgenotype Vd. This is, therefore, the first genomic and pathological characterization of NDV of subgenotype Vd. Subgenotype Vd seems to be responsible for current outbreaks in Kenya and parts of the East African region. NDV isolates, sequences, and findings obtained from the study will, therefore, be valuable in analyzing the nature and diversity of NDV in the East African region. The study findings contribute to the currently available genome data on NDV in Africa and worldwide. Through this work, we also establish local and regional reference viruses for future studies in the development of improved control and diagnostic strategies.
The full-genome sequence of the NDV isolates analyzed in this study has been submitted to GenBank under Accession nos MN685354, MN685355, and MN685356.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
This research work was supported by the International Foundation for Science (IFS), Stockholm Sweden through Research Grant 13_B/32827. The authors acknowledge the International Foundation for Science (IFS), Stockholm, Sweden, for financial support. The authors would also like to thank Moses Ogugo, Samuel Imujaro, and John Odera for their support in the field sampling and laboratory analysis. The authors also acknowledge the Country Livestock and Veterinary field staff for their technical support.
Table S1 describes the accession numbers, strain name, year, and country of isolation of reference sequences used for comparison with study sequences. The references are of known genotypes of Newcastle disease virus and are available in the GenBank including those isolated from Africa. Figure S2 shows the whole phylogenetic tree representing full coding sequences of the fusion gene of 306 strains of NDV available in the GenBank including sequences obtained from this study and all known isolates from Africa. The fusion gene is useful in the characterization of NDV. Therefore, this Maximum Likelihood tree is useful for comparison with the complete genome phylogenetic tree that is presented in Figure 4. The two trees are similar in their characterization of study isolates. (Supplementary Materials)
OIE, “Newcastle disease: general infectious disease sheets,” 2015, https://www.oie.int/fileadmin/Home/eng/Media_Center/docs/pdf/Disease_cards/NEWCAS-EN.pdf.View at: Google Scholar
D. G. Diel, L. H. A. da Silva, H. Liu, Z. Wang, P. J. Miller, and C. L. Afonso, “Genetic diversity of avian paramyxovirus type 1: proposal for a unified nomenclature and classification system of Newcastle disease virus genotypes,” Infection, Genetics and Evolution, vol. 12, no. 8, pp. 1770–1779, 2012.View at: Publisher Site | Google Scholar
S. C. Courtney, L. Susta, D. Gomez et al., “Highly divergent virulent isolates of Newcastle disease virus from the Dominican Republic are members of a new genotype that may have evolved unnoticed for over 2 decades,” Journal of Clinical Microbiology, vol. 51, no. 2, pp. 508–517, 2013.View at: Publisher Site | Google Scholar
C. J. Snoeck, A. A. Owoade, E. Couacy-Hymann et al., “High genetic diversity of Newcastle disease virus in poultry in west and central Africa: cocirculation of genotype XIV and newly defined genotypes XVII and XVIII,” Journal of Clinical Microbiology, vol. 51, no. 7, pp. 2250–2260, 2013.View at: Publisher Site | Google Scholar
M. O. Otim, H. Christensen, P. H. Jorgensen, K. J. Handberg, and M. Bisgaard, “Molecular characterization and phylogenetic study of Newcastle disease virus isolates from recent outbreaks in eastern Uganda,” Journal of Clinical Microbiology, vol. 42, no. 6, pp. 2802–2805, 2004.View at: Publisher Site | Google Scholar
I. N. Ogali, L. W. Wamuyu, J. K. Lichoti, E. O. Mungube, B. Agwanda, and S. C. Ommeh, “Molecular characterization of Newcastle disease virus from backyard poultry farms and live bird markets in Kenya,” International Journal of Microbiology, vol. 2018, Article ID 2368597, 11 pages, 2018.View at: Publisher Site | Google Scholar
Federation of Animal Science Societies, “Guide for the care and use of agricultural animals in research and teaching,” January 2010, https://www.asas.org/docs/default-source/default-document-library/ag_guide_3rded.pdf?sfvrsn=4.View at: Google Scholar
S. Akhtar, M. A. Muneer, K. Muhammad et al., “Molecular characterization and epitope mapping of fusion (F) and hemagglutinin (HN) genes of avian paramyxovirus serotype i from peacocks in Pakistan,” Pakistan Journal of Zoology, vol. 49, no. 2, pp. 755–759, 2017.View at: Publisher Site | Google Scholar
P. Yuan, K. A. Swanson, G. P. Leser, R. G. Paterson, R. A. Lamb, and T. S. Jardetzky, “Structure of the Newcastle disease virus hemagglutinin-neuraminidase (HN) ectodomain reveals a four-helix bundle stalk,” Proceedings of the National Academy of Sciences, vol. 108, no. 36, pp. 14920–14925, 2011.View at: Publisher Site | Google Scholar
D. E. Anderson and L. Wang, The Biology of Paramyxoviruses, Caister Academic Press, Norfolk, UK, 2015.
E. Wehmann, D. Ujvári, H. Mazija et al., “Genetic analysis of Newcastle disease virus strains isolated in Bosnia-Herzegovina, Croatia, Slovenia and Yugoslavia, reveals the presence of only a single genotype, V, between 1979 and 2002,” Veterinary Microbiology, vol. 94, no. 4, pp. 269–281, 2003.View at: Publisher Site | Google Scholar