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.


Introduction
Newcastle disease (ND) is a highly contagious disease of poultry that is notifiable to the World Organization for Animal Health (OIE) [1]. In developing countries, ND has a significant economic impact on the poultry sector due to the high morbidity and mortality [2]. 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 [3]. 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) [6].
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 [5]. Based on the analysis of the complete coding region of the fusion gene, class II NDV has been divided into 18 genotypes [7][8][9][10]. Newcastle disease virus belonging to class II, which was isolated between the 1930s and 1960s, includes genotypes I, II, III, and IV. ese are considered "ancient" and have genome sizes of 15,186 nucleotides (nt) [11]. e currently circulating or "recent" class II genotypes were isolated after the 1960s (V, VI, VII, VIII, and X-XVIII). e latter genotypes have an insertion of six nucleotides into the 5′ noncoding region of the nucleoprotein gene. ey, therefore, possess 15192 nt in their genomes [11]. More genetic distance has been observed between "ancient" and the "recent" genotypes [5]. Some of the genotypes are further divided into subgenotypes according to set criteria [7]. e 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 [15]. A comparative analysis of APMV-1 sequences from West Africa indicated that the viruses correspond to genotypes I and VI and novel genotype XVIII [16]. 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][20][21]. Isolates from Tanzania were assigned to genotypes I, V, and Via, respectively [22]. 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) [25]. A more recent study in Uganda isolated strains of genotype Vd from poultry in LBMs [26]. It is suggested that increased surveillance on the African Continent may lead to the discovery of more genotypes and subgenotypes [27]. is 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 [28].
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 [1]. e introduction of birds from informal live bird markets (LBMs) poses an important source of introduction of NDV into a flock [29]. 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 [30]. However, until now, only a few studies have reported the genetic and pathological characteristics of NDV circulating in Kenya. is 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.

Viruses.
Cloacal swabs and tissue samples for the study came from chicken showing clinical symptoms of Newcastle disease. ese samples were obtained from live bird markets and backyard poultry farms in Kenya during a survey of Newcastle disease virus described by [30]. We inoculated samples from individual birds in three 10-day-old specific antibody negative (SAN) embryonated chicken eggs through the allantoic cavity. ree 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. ree NDV isolates designated KE001/2015 (Makueni), KE0811/2016 (Mombasa), and KE0698/2016 (Malaba) showed high HA titer (>2 4 ). e allantoic fluid of the three isolates was aliquoted for use in pathogenicity tests and genome sequencing and characterization.

Study Clearance and Permits of Compliance.
e 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 [31]. We carried out the study in the biosafety level II isolation facility at KALRO-VSRI.

Biological Characterization of the NDV Isolates.
e 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 [1]. Briefly, ICPI was performed in 1day-old chicks by inoculation of 50 μL of allantoic fluid diluted 10-fold in PBS without antibiotics as recommended [1]. 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.
e MDT was recorded as the mean time (hours) for the minimum lethal dose of the virus to kill all the embryos [1]. We considered isolates with MDT ≤60 hours velogenic: MDT between 61 and 90 hours was mesogenic, while those with MDT >90 hours were lentogenic.

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 10 6 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. ereafter, 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. e samples were embedded into the paraffin and sectioned at 3 μm and routinely deparaffinized and stained with hematoxylin and eosin for histopathological examination.

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. e 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). e concentration and purity of the extracted RNA were determined using a Qubit Fluorimeter ( ermo Scientific, Walton, Massachusetts, USA) before storage at −80°C until use. e nucleotide sequence of the full viral genome was determined using next-generation sequencing (NGS). e 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). e reads of the full-length viral genome were assembled de novo with SPAdes assembler version 3.10.1 [32]. We confirmed the NGS results with conventional Reverse-Transcription polymerase chain reaction (RT-PCR) using virus-specific primers. is involved purification of amplified products using Qiagen PCR Purification Kit and sequencing purified products by the Sanger method at Macrogen, Korea. e sequences obtained by the Sanger method were 100% identical to the NGS sequences.

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. ese 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 [35]. We constructed phylogenetic trees using the Maximum Likelihood (ML) method and estimated the tree using bestfit 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 [35].

Biological Characterization.
e 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. e 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.

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 International Journal of Microbiology 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). ese lesions included hemorrhagic and necrotic patches and blood vessel congestion. e spleen had numerous multifocal whitish spots that were evident on the surface. e 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. e intestines showed darkened areas of necrosis that were evident on the serosa surface. ese necrotic areas in the lumen of the intestines were covered with mucoid content underneath of which were ulcerative foci. e 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. e 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). e intestinal mucosal epithelium was denuded with loss of villi structure. e findings resemble pathological presentation in infection by highly virulent viruses [36].

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. e 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. e 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. e 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.

Genomic Features of Coding and Noncoding Regions of
Isolates.
e 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. e genome of the three isolates was organized as 3′-NP-P-M-F-HN-L-5′. e 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. e 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 A H H C C Figure 1: Gross pathology lesions observed on organs of specificpathogen-free chicken infected with NDV isolates at 5 days after challenge. First row: (a) proventriculus showed patches of reddened areas on the mucosa (hemorrhagic patches), (b) the spleen was oedematous and showed multifocal white specks (necrotic points-A), and (c) the liver was oedematous and showed focal reddened hemorrhagic areas (H). Second row: ((d-f )) small intestines showed congestion of blood vessels (C), darkened patches of necrosis (arrow), and diffuse hemorrhages (H) on the mucosa. 4 International Journal of Microbiology Control KE0811 KE001 Spleen Small intestine Liver Figure 2: Photomicrographs showing hematoxylin and eosin (HE) staining on sections of the spleen, small intestines, and liver at 5 days after challenge. Magnification for panels of spleen and liver and 100x for small intestines × 400. First row: spleen tissue samples from isolate: KE0811 and KE001 showed lymphocyte depletion in the follicles and inflammatory exudates (star). Second row: small intestine tissue samples from NDV challenged groups showed multifocal mucosal hemorrhages (H), lymphocytic infiltration of the submucosa (star), and sloughing of mucosa (arrow). ird row: liver tissue samples of NDV challenged group showed congestion of hepatic blood vessels (C) and lymphocyte infiltration of periportal areas (arrow).   was a conserved sequence, intergenic sequence (IGS). e 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.

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. e 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: 171 QVWVTVAKAMT 181 ; region 2 spanning between aa 267 and277 with the motif: 267 FFLTLKYGINT 277 ; and the third conserved region between aa 322 and 336 with the motif: 322 FAPAEYAQLYSFAMG 336 . e 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. ese regions are important for the formation of the NP-P complexes and regulation of viral replication [37][38][39]. e 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 [40]. At this site, insertion of one G gives rise to the V protein, while the addition of two Gs gives rise to W protein. e V protein sequence of the study isolates aligned with that of reference genotypes of NDV. e length of the V protein is 239 aa for the NDV strains ( Figure 3). e 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   , which were similar in other NDV strains except for a few substitutions. ese sites correspond to immunodominant epitopes of the matrix protein [41]. For the M protein to be localized in the nucleus, it utilizes nuclear localization signals (NLS). For the study isolates, the NLS was 246 DKKGKKVTFDKLERKIRR 263 .
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, 151 ILRLKESIAATNEAVHEVTDG 171 . e conserved sites have been reported as important functional sites of the F gene [42,43]. e key function of the F protein is to initiate the fusion of viral surface to the host cell membrane. e amino acid motif at the F gene cleavage site is a major determinant of virulence. For the study Table 4: Comparison of the amino acid residues at the antigenic sites of the hemagglutination neuraminidase protein of the Newcastle disease virus.
Amino acid residues at different antigenic sites 1  2  3  4  12  14  23  345 513 514 521 569 263 287 321 333 356 494 516 347 350 353 193 194 201 Virus Strain (Accession) Genotype I2 (AY935499) Highlighted are sequences of isolates from Kenya that were collected in the study. e tree was constructed based on the GTR model with a discrete gamma distribution used to model evolutionary rate differences among sites (2 categories (+G), parameter � 0.33) [35]. (b) A phylogenetic tree constructed using 122 nucleotide sequences including most known African isolates. e tree is based on the General Time Reversible (GTR) model with a discrete gamma distribution used to model evolutionary rate differences among sites (2 categories (+G), parameter � 0.4054). Highlighted are sequences of Kenyan isolates collected in the study.
isolates, the cleavage site had the motif 112 R-R-Q-K-R-F 117 corresponding to the cleavage site of velogenic viruses. e 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 (  [44]. is was similar in NDV strains except three substitutions at residues V81I, N98S, and T102I observed in the study isolates. e sialic acid binding site of the HN protein [45] 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 234 NRKSCS 239 . 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. ese sites correspond to the neuraminidase active site that binds the fusion protein [46]. e RNA-dependent RNA polymerase (L) gene of the study isolates and other NDV strains had a highly conserved region between aa 740 and 753 ( 740 SHCRVACMVQGDNQ 753 ) which corresponds to the putative active site of the L gene essential to polymerase function [47]. Other conserved domains of the L gene were observed between aa 503 and 607 and 418 between aa 634 and 854. e 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). e 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 [5]. e 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)).

Conclusion
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. is 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. e study findings contribute to the currently available genome data on NDV in Africa and worldwide. rough this work, we also establish local and regional reference viruses for future studies in the development of improved control and diagnostic strategies.
Data Availability e 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
e authors declare that there are no conflicts of interest regarding the publication of this paper. Table S1 describes the accession numbers, strain name, year, and country of isolation of reference sequences used for comparison with study sequences. e 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. e fusion gene is useful in the characterization of NDV. erefore, this Maximum Likelihood tree is useful for comparison with the complete genome phylogenetic tree that is presented in Figure 4. e two trees are similar in their characterization of study isolates. (Supplementary Materials)