key: cord-0767161-mvdbr0eq
authors: Abd El-Hamid, Hatem S.; Shafi, Manal E.; Albaqami, Najah M.; Ellakany, Hany F.; Abdelaziz, Naglaa M.; Abdelaziz, Mohamed N.; Abd El-Hack, Mohamed E.; Taha, Ayman E.; Alanazi, Khalid M.; Elbestawy, Ahmed R.
title: Sequence analysis and pathogenicity of Avian Orthoavulavirus 1 strains isolated from poultry flocks during 2015–2019
date: 2020-07-22
journal: BMC Vet Res
DOI: 10.1186/s12917-020-02470-9
sha: 8022bb0ad4af462bf0fc71c1642d0bb501943b63
doc_id: 767161
cord_uid: mvdbr0eq

BACKGROUND: Newcastle disease (ND) causes severe economic losses in poultry industry worldwide. Egyptian poultry industry suffered from severe economic losses since the isolation of Velogenic Newcastle disease virus (vNDV) genotype VIId in 2011 and up till now despite the use of different vaccination programs. So, this study aimed to isolate and characterize the vNDV from a total of 120 poultry flocks from ten provinces in the Egyptian Delta region with a history of respiratory manifestation, high mortalities or a decrease in egg production between 2015 and 2019. Seventy-three samples’ allantoic fluid (73/120, 60.8%) were positive for hemagglutination with chicken RBCs. These samples were submitted to molecular examination using qRT-PCR specific primers for AOAV-1, highly pathogenic avian influenza (HPAI-H5), low pathogenic avian influenza (LPAI-H9) and infectious bronchitis virus (IBV). RESULTS: Fifty samples (50/120: 41.6%) were confirmed positive for AOAV-1, based on genetic analysis of matrix and fusion protein. The co-infection rate of other respiratory viral diseases examined was 1.6, 14.1, and 4.1%, for HPAI-H5, LPAI-H9, and IBV, respectively. Biologically, the intracerebral pathogenicity index of ten selected AOAV-1 isolates ranged from 1.70 to 1.98, which indicated the velogenic nature of these isolates. All the sixteen sequenced isolates were AOAV-1 genotype VII.1.1. The full F gene sequence of six examined AOAV-1 VII.1.1 isolates contained the seven neutralizing epitopes, and the glycosylation motif of six-potential sites for N linked glycosylation at residues 85, 191, 366, 447, 471, and 541. CONCLUSION: It could be concluded that the high prevalence of AOAV-1 genotype VII.1.1 in the Egyptian chicken flocks despite the intensive vaccination with live and killed ND vaccines, as all the 16 isolates tested were belonged to this genotype. Homologous vaccination is badly needed to control and reduce the spread of AOAV-1 genotype VII.1.1infection in Egyptian poultry flocks.

Newcastle disease (ND) encompasses a critical impact amongst the most economically viral poultry diseases leading to high mortality and morbidity rates in the susceptible poultry. It still causing a drop in egg production even in vaccinated layers [1] . Besides the economic impact through the loss of productive assets, trade restrictions are of significant concern for exporting countries. Within the last 5 years, 109 of 200 member countries have reported the disease to the OIE [2] .

Newcastle disease virus (NDV) or Avian paramyxoviruses 1 (APMV-1) was recently classified to Genus Avian Orthoavulavirus, family Paramyxoviridae and common NDV is recently known as Avian Orthoavulavirus-1 (AOAV-1) [3, 4] . Full fusion (F) gene nucleotide sequencing classified all AOAV-1 into two classes; class I and class II. AOAV-1 class II contains virulent, and nonvirulent viruses, with 21 identified genotypes (I to XXI Genotype). Updated classification criteria of AOAV-1 genotype VII viruses merged all the subgenotypes into only 3, VII.1.1 including previous subgenotypes (b, d, e, j, l), VII.1.2 included previous subgenotypes (f), and VII.2 included the previous subgenotypes (h, i, k) [5] .

AOAV-1 (genotype VII) causes fatal infections in poultry and other susceptible birds as they are responsible for the fourth major panzootic of ND worldwide which caused by viruses of genotype VII.1.1 [6] . Similarly, in the Egyptian poultry industry, AOAV-1 is specially subgenotypes VIId and VIIb are incriminated in significant economic losses. The partial sequence of the F gene and phylogenetic analysis of AOAV-1 isolated from different outbreaks in Egypt from 2011 to 2017 belonged to genotype VIId [7] [8] [9] [10] .

Among genotype VII viruses, sub-genotype VIIi has demonstrated an intercontinental spread and, therefore, has a global significance in the perspective of potential fifth panzootic [6] . Most of the AOAV-1 viruses that are pathogenic for chickens have the sequence 112R/K-R-Q/ K/R-K/R-R116 at the C-terminus of the F2 protein and F (phenylalanine) at residue 117 at the N-terminus of the F1 protein. In contrast, the viruses of low virulence have sequences in the same region of 112G/E-K/R-Q-G/ E-R116 and L (leucine) at residue 117 [11, 12] . Thus, there is a need for at least one pair of basic amino acids at residues 116 and 115 plus phenylalanine at residue 117 and a basic amino acid (R) at 113 to enable the virus to show virulence for chickens.

Study of F protein full length during outbreaks of 2016 in poultry flocks in Egyptian provinces revealed that genotype VIId viruses of class II are the most predominant AOAV-1 stains [13] . The phylogenic analysis of AOAV-1 strains circulated in small-scale poultry holdings in Egyptian villages revealed them belonging to genotype VIIb. Also, the viral transmission occurred among neighboring farms over long distances, and coinfections with multiple pathogens were identified [14, 15] . Circulating AOAV-1 in wild birds in borderline or frontier Egyptian provinces containing wetlands was recorded [16] . Consequently, this study aimed to investigate the molecular and biological characterization of the AOAV-1 strains from 2015 through 2019 in the Egyptian Delta region provinces.

Isolation of AOAV-1 in SPF-ECE and slide HA testing

The embryos of inoculated eggs were died between the 2nd and the 4th days post-inoculation in many samples. Samples with no mortality after the 3rd passage that had no HA activity were confirmed as negative for haemagglutinating viruses. A total of 73 samples were positive with slide HA test constituting 60.8% of total AF samples ( The prevalence of AOAV-I and other viruses in provinces of the Egyptian Delta region Fifty AOAV-1 isolates (50/120: 41.6%) were confirmed, 2/120 isolates (1.6%) were positive for HPAI-H5N1, 17/120 (14.1%) isolates were positive LPAI-H9N2 and 5/120 (4.1%) were IBV positive. The mixed infection of AOAV-1 with H9N2 appeared in two samples. One positive AOAV-1 sample was mixed with H5N1 infection. Two samples of AOAV-1 were mixed with IBV. One AOAV-1 sample was mixed with both H9N2 and IBV. From all the 73 HA positive samples, only seven samples were negative for all of the four viruses under study. The prevalence of virulent AOAV-1 remains the highest higher in the broiler farms in the Delta region representing 41.1% (37/90 flocks) of AOAV-1 isolates, followed by layer flocks 5/13 (38.4%). Despite the lower examined flocks of Balady, Sasso chickens, and pigeons; the confirmed AOAV-1 was high as 5/7 (71.4%) in Balady chickens; 1/3 (33.3%) Sasso chickens and 2/4 (50%) in pigeons. According to the geographical order: incidence of AOAV-1 in the ten Egyptian province farms was as follows: El Beheira 41.2%, Alexandria 30.7%, El Gharbia 53.3%, El Qaliobia 22.2%, El Dakahlia 42.8%, Kafer El Shiekh 66.6%, El Menofia 66.6% and El Giza 0% Marsa Mattrouh 50% and El Sharkia 50% (Table 1) .

Testing only ten AOAV-1 isolates (representing six Egyptian provinces) by ICPI gave values of 1.70-1.98 for all the ten isolates, which indicates the velogenic nature of these isolates (Table 2) .

Sequence analysis of fusion protein, phylogenetic tree and molecular pathotyping of AOAV-1 isolates

The obtained full-length F protein nucleotides sequences of six AOAV-1 isolates submitted for BLASTN analysis revealed 99.2% identity with NDV-Egy-Beh-ck-VII-b-2016-NR726 and NDV-EG-35-2014 AOAV-1 isolates' The comparative alignment of the deduced amino acid of the first six AOAV-1 isolates confirmed the presence of multiple basic amino acids at the positions of 112 to 116 and F phenylalanine aa at the position of 117 (the cleavage site motif of virulent strains RRQKRF) along with the conserved amino acids of K101 and V121. Both are characters of more virulent AOAV-1 viruses confirming the velogenic pathotyping of these isolates (Table 3 and Figs. 1, 2 and 3).

All these six AOAV-1 isolates from 2015 to 2019 submitted for full F genomic sequencing were recognized as AOAV-1 genotype VII.1.1. (Fig. 4 ) based on the complete F protein and according to the classification of Dimitrov et al. [5] . Those six AOAV-1 isolates were clustered with subgenotype VIIj except of the isolate MH445410 which was clustered in the subgenotype VIIb ( Fig. 5) . Also, the other ten AOAV-1 isolates partially sequenced for the F gene were also clustered with the subgenotype VIIj according to Diel et al. [17, 18] .

The nucleotides alignment of the complete F protein sequence (1-1659) and the phylogenetic tree of the six isolates confirmed the high similarity and close relatedness by 98.1-99.3% to the Egyptian isolates from 2011 to 2015 and by 97.4-99.3% identity to the Egyptian isolates from 2016 to 2019. Also, these six isolates showed a high similarity (98.3-99.1%) with the other Middle East strains of genotype VII as well as with the Chinese strains, and with a lower identity of 94.1-95.1% with the Korean strains.

Regarding the nucleotide identity between those six AOAV-1 genotypes, VII.1.1 isolates and the vaccine strains used in the prevention and control of the AOAV-1 infection in Egypt was so far, and this nucleotide identity varied from one vaccine strain to another. For class II, genotype II vaccine strains, an 82.8-83.2% identity was recorded with LaSota, 82.9-83.4% with both Clone 30 and Hitchner B1, and 85.2-86% with VG/GA (Avinew®). While the relatedness between our six isolates and the genotype I vaccine PhyLM42 strains were 84.6-85.2% with V4 which was 85.2-86%, and with D26/76 

which was 85.1-86%. Whilst, the relatedness of these six isolates with genotype VII vaccine (KBNPC415R2L) F gene had higher nucleotide identity that ranged from 93.6 to 94.5%, indicating a closer genetic relationship (Fig. 6 ).

Antigenic analysis of the AOAV-1 isolates 

Newcastle Disease (ND) or AOAV-1 is a highly contagious viral disease that affects more than 240 domestic and wild bird species, with a higher incidence in domestic poultry, resulting in severe outbreaks worldwide. It represents a major depletion in the poultry economy than any other viral disease [5, 19] . To date, the disease occurs on at least six of the seven continents [6] , being an enzootic in many countries, with variable outcomes and subsequently variable economic losses.

Class II of AOAV-1 is more diverse, containing virulent and non-virulent viruses, and the complete analyses identified at least 21 distinct genotypes [5] . Genotype VII of the class II viruses had been associated with the most recent outbreaks in Asia, Europe, Africa, Middle East, and South America [20] [21] [22] [23] [24] . Genotype VII can be further divided into 12 (VIIa-l) subgenotypes according to amino acid substitutions [5, 25] . The previous studies regarding the partial sequence of F protein and phylogenetic analysis of AOAV-1 in the Egyptian poultry flocks revealed that the predominant strain circulating in the poultry field was viruses of subgenotypes b and d [13-16, 26, 27] . The clinical manifestations and PM lesions of ND caused by the confirmed AOAV-1 in the examined poultry flocks as respiratory manifestations, enteric symptoms (greenish diarrhea) and nervous manifestation with high mortalities and decreased egg production and quality [amorphous (oddly shaped), smaller eggs with rough, thin, and whiter bleached shells and poor-quality contents whites and yolks] were similar to those findings of Miller and Koch [5] , Terregino and Capua [28] , Susta et al. [29] , Bertran et al. [30] and Igwe et al. [31] .

Herein, a total of 50/120 (41.6%) samples were positive for AOAV-I, indicating the higher incidence rate of this virus in the Egyptian poultry flocks. The higher rate of isolation of AOAV-1 from infected flocks collected during the winter season in the Delta region may be attributed to increasing the environmental load of AOAV-1 through the high viral stability during the cold and wet Wang et al. [35] studied the risk factors of infectious diseases for poultry in backyard household flocks in rural China. Authors found that they were significantly affected by neighboring commercial poultry and close contact with wild birds. The role of both free-living pigeons in the villages nearby the chicken farms as well as ducks reared in house backyards can't be neglected as a risk factor for chickens in the commercial sector for the NDV infection as most of the sampled farms have open house ventilation directly exposed to the outer environmental air draughts and dust. Artois et al. [36] investigated the risk factors of Avian Influenza in the Northern Egyptian Delta Governorates. They documented that the probability of HPAI (H5N1) presence was associated with a risk factor of the existence of duck farms with significant mortality or drop in egg production in the village near the chicken commercial farms.

Two previously published papers showed the ability of pigeons and ducks in transmitting this virus to chickens. The first study by Ellakany et al. [37] confirmed that incontact chickens to infected pigeons developed severe respiratory, digestive and nervous signs and suffered from the mortality of 2/5 and 3/5 in the in-contact chickens to pigeons infected IM (intra-muscularly) and IN (intra-nasally) respectively. In-contact chickens to infected pigeons, either IM or IN showed high viral shedding titers in both oropharynx and cloaca up to 16 days post-exposure to the infected pigeons. From these results, it could be concluded that free-range pigeons can be considered as an efficient reservoir for NDV-VIId to opened-house raised broilers. A similar experiment by Elbestawy et al. [38] proved that Muscovy ducks infected with the velogenic (AOAV-1) NDV genotype VIId, by the intra-nasal route, shed the virus in their tracheas for up to 16 days PI and their in-contact infected chickens showed clear characteristic signs for NDV and 20% mortality. These in-contact chickens were still shedding the virus in their oropharyngeal and cloacal excreta for up to 11 days post-exposure to the infected ducks. This proved that Muscovy ducks are efficient carriers and considered as a risk factor for commercial chickens for infection with velogenic AOAV-1 -genotype VIId under the conditions of raising poultry in Egyptian villages. Otim et al. [39] evaluated the risk factors for Newcastle disease in Uganda villages. They found that purchasing restocking chickens from the market and neighborhood (hazard ratio [HR] = 1.79), the presence of migratory wild birds (HR = 1.70), and being in agro-ecological zone 1 (HR = 1.48) showed a positive but non-significant association with the risk for ND.

According to the OIE [2] , NDV virulent strains (mesogenic and velogenic) possess 112R/K-R-Q-R/K-R-F117 motif with two pairs of basic amino acids and a Phenylalanine in the N terminus of the F1 (residue 117). The second OIE accepted criteria for virulent NDVs is the ICPI, in one-day-old chicks [5, 40] . AOAV-1 isolates were classified as virulent (mesogenic or velogenic) pathotype when the ICPI score was equal or greater than 0.7 [2] . In this study, the ICPI value for ten confirmed AOAV-1 isolates gave an ICPI score between 1.70 and 1.98 for all the isolates indicating their virulent pathotype nature. Previously, the ICPI values of most Egyptian AOAV-1 genotype VII viruses were recorded to be between 1.66 and 1.97 [9, 16, 27, 41, 42] . This increase in the ICPI score for the Egyptian isolates of AOAV-1 over the last 4 years indicates an increase in their virulence.

Regarding the molecular pathotyping of the other six selected isolates of AOAV-1 chosen for complete F gene sequencing, the comparative alignment of deduced amino acid revealed the presence of two pairs of multiple basic amino acids R/K at positions of 112 to 116 and F phenylalanine aa at position 117 which is the cleavage site motif of virulent strain RRQKRF. Also, the comparative alignments of all cleavage site motifs revealed that these six virulent strains contained the conserved amino acids at positions K101 and V121 of the fusion (F) protein which is the unique features of genotype VII) which agreed with the previous study of Lien et al. [43] . The presence of Q in the motif of virulent strain RRQKRF enhances and increases the virulence [44] . The full-length F gene sequencing of six and the partial sequencing of the other ten AOAV-1 isolates, respectively, showed 97.4-99.3% identity with the previously published Egyptian isolates between 2011 and 2019. This indicates the high prevalence of the AOAV-1 genotype VII.1.1 in the Egyptian poultry flocks.

Another test for virulence detection is the F protein sequence analysis of AOAV-1 which is a class I integral membrane protein present as a trimer in the virion. Initially, F0 is synthesized as an inactive precursor and subsequently cleaved by the cellular proteases into a disulfide-linked F1-F2 complex. The deduced amino acid of the functional domain of F protein of AOAV-1 contains signal peptide from position (1-31 aa), F2 subunit from position (32-116aa) contains the cleavage site (112aa to 116aa). The F1 subunit contains the fusion peptide FP from position (117 aa to 136aa), three hydrophobic heptad repeat (HR) domains (HRa 143aa-185aa, HRb 268aa − 299 aa and HRc 471aa − 500 aa), the transmembrane (TM) domain from (501 aa to 522aa) and cytoplasmic tail (523 to 553 aa) are important for viral infectivity and pathogenicity [45] . The amino acid substitution in neutralizing sites of F protein, especially substitutions in fusion peptides and hydrophobic regions and transmembrane regions, had alteration effect on the fusion activity. It helped the virulent virus to escape neutralizing antibodies as a variant virus [46] [47] [48] . Also, cytoplasmic tails had an importance for virus replication and vaccine competence [49] .

All glycosylation sites were conserved in all tested isolates, which have an important role in the infectivity and antigenicity [50] . A signal peptide region substitution occurred in three isolates MK984236, MK984237, and MK984238 in the position of S30N sharing the same amino acid N30 of genotype II of vaccinal strains. This was previously studied in the Egyptian strain of genotype VIId [13] , suggesting that the use of genotype II vaccine type in control of AOAV1 genotype VII.1.1 in Egypt may be implicated in the emergence of new variants rather than providing benefits against NDV infections. The HRa region in MK984236, and MK984238 had a substitution in the position of D170N that differ from conserved amino acid of neutralizing epitope strain MK984239. The same results were obtained by Selim et al. [27] who examined seven Egyptian strains of genotype VIId and found alterations at different sites of the F gene interest as N-glycosylation sites, epitopes binding sites cysteine residues that may affect virus pathogenicity and may interfere with the classical vaccine's protection. Isolate no. MH445410 had the highest amino acid substitutions rate.

All signal peptides of the six tested AOAV-1 genotype VII.1.1 isolates contained the amino acid at positions of I9, A11, R18, C25 and D170 except for the two strains MK984236 and MK984238, which had substitution at position D170N sharing these amino acids characteristic for subgenotypes (VII b, d, e, g), that represent viruses of 4th panzootic, while the 5th panzootic of AOAV-1 infections were belonged to viruses of subgenotype VII (h and i) [6, 51] .

The higher prevalence of AOAV-1 genotype VII.1.1 (with highly identical F protein sequence) indicated that ND remains enzootic despite the extensive use of live and inactivated vaccines to control the disease. Comparing to the NDV vaccine strains used in Egypt, the nucleotide identity of the six isolates of AOAV-1 genotype VII.1.1 varied from 82.8 to 86% with the [VG/GA strains (Avinew) vaccine and vaccines of genotype I were slightly more identical than those of genotype II, LaSota, clone 30 and Hitchner B1]. Roohani et al. [52] reported that NDV genotype VII strains were higher divergent from the genotypes I or II based vaccines. Also, Xue et al. [53] recorded that the phylogenetic analyses of 11 isolates of genotype VII NDV showed low amino acid similarity with the genotype I and II vaccine strains emphasizing the necessity of development of a genotypematched vaccine. Interestingly, the nucleotide identity with the F protein sequence of genotype VII vaccine (KBNPC415R2L) used in Egypt ranged from 93.6 to 94.5%, indicating the probability of higher protection of this type of vaccine, having the advantage of decreasing the virus multiplication and shedding into the environment to reduce the spread of this devastating viral infection among poultry flock.

It had been established that homologous genotype vaccines are more effective in preventing virus replication and shedding after challenge than heterologous NDV genotype VII vaccine [52, 54, 55] . However, Sultan et al. [56] recommended the combined vaccination regime of inactivated genotype VII-matched and live genotype II vaccines for three times along the rearing period in commercial layers for better protection against clinical disease, mortality and viral shedding confirming the results of Sedeik et al. [57] . The aforementioned authors reported that the use of either inactivated genotype VII or II once alone induced the insufficient protection against clinical disease, mortality, and viral shedding in broiler chickens.

Based on the partial and complete sequence of F gene and according to the updated classification of AOAV-1, we can conclude the high prevalence of AOAV-1 genotype VII.1.1 in the Egyptian chicken flocks despite the intensive vaccination with live and killed ND vaccines, as all the 16 isolates tested belonged to this genotype. Current vaccination strategies in the Egyptian chicken flocks don't appear to control effectively AOAV-1 genotype VII.1.1, given the high prevalence of this genotype in the vaccinated birds surveyed here. Better vaccination targets are needed to control this genotype and reduce the spread of this devastating viral infection in Egyptian poultry flocks.

A total of 120 various poultry flocks or backyard household places (broiler, broiler breeder, commercial layer, native Saso and Balady chickens and pigeons) from 10 Egyptian provinces (El Beheira, Alexandria, El Gharbia, EL Qaliuobia, El Dakahlia, Kafr El Shiekh, El Menofia, El Giza, Marsa Mattrouh, and El Sharkia) were examined during 2015-2019 ( Table 4 ). The locations were identified following consultation with farm owners who had contacted us to discuss their current situation. The household sources capacity started from 50 birds, while the organized commercial flocks' capacity ranged from 3000 to 70.000 birds. The age of commercial broilers ranged from 18 to 36 days; layers: 52-315 days; Balady: 17-56 days; pigeons: 220-300 days; broiler breeders: 55-78 days and Saso: 45-55 days old.

These examined poultry flocks suffered from high mortality ranged from 10 to 35% within the last week before sampling. Typical clinical signs and post-mortem lesions suspecting ND were seen before sampling. All flocks were vaccinated at least with live with or without inactivated vaccine regimens for NDV. Samples were collected from trachea, lung, spleen, kidney, liver, and brain tissues and/or swabs of sick and freshly dead birds. Samples were labeled and transported immediately on ice either to the laboratory of poultry and fish diseases department, Faculty of Veterinary Medicine Damanhour University, or Reference Laboratory for Veterinary Quality Control on Poultry production (RLQP) Biotechnology unit, Animal Health Research Institute. A single sample pool for each examined flock was processed except the brain tissues were processed separately according to the protocol of OIE [2] .

Tissue homogenate and/or swabs from 3 to 5 birds/ each flock (as a pooled sample) was inoculated via allantoic sac of 3 Specific Pathogen Free Embryonated Chicken Eggs (SPF-ECE) 10 days old (Kom Oshem, SPF Farm, El Fayoum) according to OIE [2] . The harvested allantoic fluid of each sample was tested for haemagglutination activity by slide HA test after each passage [28] .

Total viral RNA was extracted from all positive HA allantoic fluids using QIAamp Viral RNA Mini Kit (QIAG EN, Germany) according to manufacturer instructions.

Detection of virulent AOAV-1 and co-infection with H5, H9, and IBV by quantitative or real-time reverse transcriptase polymerase chain reaction (qRT-PCR)

Oligonucleotide primers and probes used were supplied from Metabion; PCR reaction was conducted using realtime PCR (Step one Applied biosystems) according to manufacturer instructions. The oligonucleotide primers and probes used according to Wise et al. [58] and Creelan et al. [59] for NDV matrix and velogenic F gene, respectively; L ndt et al. [60] for highly pathogenic avian influenza (HPAI), H5; Ben Shabat et al. [61] for low pathogenic avian influenza (LPAI) H9 and Callison et al. [62] for infectious bronchitis virus (IBV).

Partial F gene sequencing for 10 AOAV-1 suspected isolates using primer sets designed by Selim et al. [27] was performed. Also, a full gene sequencing was carried on for the other six suspected AOAV-1 isolates rather than the above mentioned ten strains. Purification through QIAquick Gel Extraction Kit (Qiagen, Germany) and the one-step rRT-PCR reaction was used according to the manufacturer instructions. Sequence analysis using Big Dye Terminator V3.1 cycle sequencing kit (Perkinelmer, Foster city, CA) and Applied Biosystems 3130 genetic analyzer (ABI, USA) were done.

Bioinformatic and phylogenetic analysis of F protein of AOAV-1

Relationships of the full F gene were compared with previously published NDV reference strains available in the public database (BLASTn, NCBI, USA) (http://www. ncbi.nlm.nih.gov/BLAST), using Bioedit software version 7.2.4. Nucleotides similarity and divergence were conducted in the MegAlign program of the Laser gene package (DNASTAR Inc., Madison, WI, USA). Molecular Evolutionary Genetics Analysis (MEGA) X version 2018 across computing platforms Software [63] was used in the construction of phylogenetic trees for the F gene nucleotide and amino acid sequences of the isolates. Algorithms to a matrix of pairwise distances estimated using the Maximum Composite Likelihood (MCL) approach and Tamura-Nei model [64] .

Intracerebral Pathogenicity Index (ICPI) was carried out for evaluation of the pathogenicity of 10 isolates of AOAV-1 from different provinces according to procedures of OIE [2] . 

This project was funded by the Deanship of Scientific Research (DSR) at King Abdulaziz University, Jeddah; the authors, therefore, acknowledge with thanks DSR for technical and financial support.

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request. All the sequences are registered on the NCBI database, and these are their links: AOAV1-Eg-Ch- 

Newcastle disease. Diseases of poultry. 14th. editor-in-chief, David E. Swayne; associate editors, Martine Boulianne

Chapter 3. 3.14. In: Manual of Diagnostic Tests and Vaccines for Terrestrial Animals

Taxonomy of the order Mononegavirales: update

ICTV: International committee on taxonomy of viruses. Virus Taxonomy: (2019) 2018b Release, Available at

Updated unified phylogenetic classification system and revised nomenclature for Newcastle disease virus

Identification of new subgenotypes of virulent Newcastle disease virus with potential panzootic features

Isolation, molecular characterization and pathotyping of Newcastle disease viruses from field outbreaks among broiler flocks in Egypt from 2014-2015

Detection and characterization of Newcastle disease virus in formalin-fixed, paraffin-embedded tissues from commercial broilers in Egypt

Genotypic and pathogenic characterization of genotype VII Newcastle disease viruses isolated from commercial farms in Egypt and evaluation of heterologous antibody responses

Avian influenza virus and Newcastle virus surveillance and characterization in broiler and layer chicken flocks in Egypt

Phylogenetic diversity among low virulence Newcastle disease viruses from waterfowl and shorebirds and comparison of genotype distributions to poultry-origin isolates

Antigenic and immunogenic investigation of the virulence motif of the Newcastle disease virus fusion protein

Evolutionary insights into the fusion protein of Newcastle disease virus isolated from vaccinated chickens in 2016 in Egypt

Isolation and identification of H9N2 avian influenza and Newcastle disease viruses co-infections in chicken

Investigation of suspected Newcastle disease (ND) outbreaks in Egypt uncovers a high virus velogenic ND virus burden in small-scale holdings and the presence of multiple pathogens

Biological characterization of wild-bird-origin avian avulavirus 1 and efficacy of currently applied vaccines against potential infection in commercial poultry

Genetic diversity of avian paramyxovirus type 1: proposal for a unified nomenclature and classification system of Newcastle disease virus genotypes

Characterization of Newcastle disease viruses isolated from cormorant and gull species in the United States in 2010

Molecular evolution of the Newcastle disease virus matrix protein gene and phylogenetic relationships among the paramyxoviridae

Phylogenetic and biological characterization of Newcastle disease virus isolates from Pakistan

Newcastle disease outbreaks in recent years in Western Europe were caused by an old (VI) and a novel genotype (VII)

Evolutionary dynamics of Newcastle disease virus

Biological and phylogenetic characterization of a genotype VII Newcastle disease virus from Venezuela: efficacy of field vaccination

Phylogenetic and pathotypical analysis of two virulent Newcastle disease viruses isolated from domestic ducks in China

NDV subgenotype VII (L) is currently circulating in commercial broiler farms of Iran

Isolation and molecular characterization of Newcastle disease virus genotypes II and VIId in Egypt between

Molecular characterization of full fusion protein (F) of Newcastle disease virus genotype VIId isolated from Egypt during 2012-2016

Conventional diagnosis of Newcastle disease virus infection. Pages 123-125 in avian influenza and Newcastle disease. A field and laboratory

Clinicopathological characterization in poultry of three strains of Newcastle disease virus isolated from recent outbreaks

Avian influenza virus and Newcastle disease virus

Pathology and distribution of Velogenic Viscerotropic Newcastle disease virus in the reproductive system of vaccinated and unvaccinated laying hens (Gallus gallus domesticus) by Immunohistochemical Labelling

Natural co-infection caused by avian influenza H9 subtype and infectious bronchitis viruses in broiler chicken farms. Veterinarski Archive

Impact of treatments for recycled broiler litter on the viability and infectivity of microoganisms

Repeated isolation of virulent Newcastle disease viruses in poultry and captive nonpoultry avian species in Pakistan from

Risk factors for infectious diseases in backyard poultry farms in the Poyang Lake area

Avian influenza A (H5N1) outbreaks in different poultry farm types in Egypt: the effect of vaccination, closing status and farm size

Role of pigeons in the transmission of avian Avulavirus (Newcastle disease-genotype VIId) to chickens

Muscovy ducks infected with velogenic Newcastle disease virus (genotype VIId) act as carriers to infect in-contact chickens

A study of risk factors associated with Newcastle disease epidemics in village freerange chickens in Uganda. Trop Anim Health Prod

Rapid pathotyping of Newcastle disease virus by pyrosequencing

Emergence of a new genetic lineage of Newcastle disease virus in west and Central Africa-implications for diagnosis and control

Pathogenesis of new strains of Newcastle disease virus from Israel and Pakistan

Phylogenetic characterization of Newcastle disease viruses isolated in Taiwan during

Comprehensive analysis of amino acid sequence diversity at the F protein cleavage site of Newcastle disease virus in fusogenic activity

Mutations in the fusion peptide and heptad repeat regions of the Newcastle disease virus fusion protein block fusion

Characterization of a recombinant Newcastle disease vaccine strain

A vaccine candidate of attenuated genotype VII Newcastle disease virus generated by reverse genetics

Characterization of complete genome sequence of genotype VI and VII velogenic Newcastle disease virus from Japan

Genotype matched Newcastle disease virus vaccine confers improved protection against genotype XII challenge: The importance of cytoplasmic tails in viral replication and vaccine design

N-glycans on Nipah virus fusion protein protect against neutralization but reduce membrane fusion and viral entry

Molecular characterization of new emerging sub-genotype VIIh Newcastle disease viruses in China

Characterization of genotype VII Newcastle disease virus (NDV) isolated from NDV vaccinated chickens, and the efficacy of LaSota and recombinant genotype VII vaccines against challenge with velogenic NDV

Genetic diversity of the genotype VII Newcastle disease virus: identification of a novel VIIj sub-genotype

Antigenic differences among Newcastle disease virus strains of different genotypes used in vaccine formulation affect viral shedding after a virulent challenge

A novel genotype VII Newcastle disease virus vaccine candidate generated by mutation in the L and F genes confers improved protection in chickens

Protective efficacy of the Newcastle disease virus genotype VII-matched vaccine in commercial layers

Comparative efficacy of commercial inactivated Newcastle disease virus vaccines against Newcastle disease virus genotype VII in broiler chickens

Development of a real-time reverse-transcription PCR for detection of Newcastle disease virus RNA in clinical samples

Detection and differentiation of pathogenicity of avian paramyxovirus serotype 1 from field cases using one-step reverse transcriptase-polymerase chain reaction

Pathogenesis of highly pathogenic avian influenza A/Turkey/Turkey/1/2005 H5N1 in Pekin ducks (Anas platyrhynchos) infected experimentally

Development of a real-time TaqMan RT-PCR assay for the detection of H9N2 avian influenza viruses

Molecular characterization of infectious bronchitis virus isolates foreign to the United States and comparison with United States isolates

MEGA X: molecular evolutionary genetics analysis across computing platforms

MEGA6: molecular evolutionary genetics analysis version 6.0

Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations