Phylogenetic assessment reveals continuous evolution and circulation of pigeon-derived virulent avian avulaviruses 1 in Eastern Europe, Asia, and Africa
BMC Veterinary Research volume 13, Article number: 291 (2017)
The remarkable diversity and mobility of Newcastle disease viruses (NDV) includes virulent viruses of genotype VI. These viruses are often referred to as pigeon paramyxoviruses 1 because they are normally isolated and cause clinical disease in birds from the Columbidae family. Genotype VI viruses occasionally infect, and may also cause clinical disease in poultry. Thus, the evolution, current spread and detection of these viruses are relevant to avian health.
Here, we describe the isolation and genomic characterization of six Egyptian (2015), four Pakistani (2015), and two Ukrainian (2007, 2013) recent pigeon-derived NDV isolates of sub-genotype VIg. These viruses are closely related to isolates from Kazakhstan, Nigeria and Russia. In addition, eight genetically related NDV isolates from Pakistan (2014–2016) that define a new sub-genotype (VIm) are described. All of these viruses, and the ancestral Bulgarian (n = 2) and South Korean (n = 2) viruses described here, have predicted virulent cleavage sites of the fusion protein, and those selected for further characterization have intracerebral pathogenicity index assay values characteristic of NDV of genotype VI (1.31 to 1.48). A validated matrix gene real-time RT-PCR (rRT-PCR) NDV test detect all tested isolates. However, the validated rRT-PCR test that is normally used to identify the virulent fusion gene fails to detect the Egyptian and Ukrainian viruses due to mismatches in primers and probe. A new rapid rRT-PCR test to determine the presence of virulent cleavage sites for viruses from sub-genotypes VIg was developed and evaluated on these and other viruses.
We describe the almost simultaneous circulation and continuous evolution of genotype VI Newcastle disease viruses in distant locations, suggesting epidemiological connections among three continents. As pigeons are not migratory, this study suggests the need to understand the possible role of human activity in the dispersal of these viruses. Complete genomic characterization identified previously unrecognized genetic diversity that contributes to diagnostic failure and will facilitate future evolutionary studies. These results highlight the importance of conducting active surveillance on pigeons worldwide and the need to update existent rapid diagnostic protocols to detect emerging viral variants and help manage the disease in affected regions.
Virulent Newcastle disease viruses (NDV), synonymous with avian avulaviruses 1 (AAvV-1), are the causative agents of Newcastle disease (ND). Newcastle disease is one of the most important infectious diseases of poultry because of its worldwide distribution and devastating economic effects for the poultry industry [1,2,3]. The disease is highly contagious, presents high morbidity and mortality, and is classified as a notifiable disease by the World Organisation for Animal Health (OIE) . AAvV-1 (along with another 12 serotypes, namely AAvV 2–13) belongs to genus Avulavirus of the family Paramyxoviridae [4, 5]. Several additional AAvV, have been recently proposed as potential new serotypes [6,7,8,9]. Avian avulaviruses 1 are enveloped, have a single stranded, non-segmented, negative sense RNA genome with helical capsid symmetry . Three different genomic sizes (15,186 nucleotides [nt], 15,192 nt and 15,198 nt) have been identified so far, and all known AAvV-1 are divided into two major genetic groups, class I and class II . Currently, AAvV-1 isolates are classified into 18 class II and one class I genotypes based on the complete coding sequences of the fusion protein [12, 13]. For the purpose of this work, the taxa name Newcastle disease virus will be used.
Pigeons may be infected with NDV of all genotypes, but are particularly susceptible to the genotype VI genetic variants, also referred to as pigeon paramyxovirus 1 (PPMV-1). These genotype VI viruses are mainly isolated from Rock Pigeons (Columba livia), but also have been isolated from other members of the family Columbidae, e.g. feral Eurasian Collared Doves (Streptopelia decaocto) . Viruses of genotype VI are endemic in the pigeon population throughout the world, and can be distinguished as being variants of NDV by the patterns produced in a hemagglutination inhibition (HI) assay [15, 16] when tested with a panel of monoclonal antibodies. Genotype VI viruses have evolved rapidly since their emergence in the Middle East during the 1960s and are currently divided into at least eleven sub-genotypes (and one more putative for a group of viruses identified in Ethiopia), namely VIa-VIj (respectively VIl) [13, 17].
Commonly, these viruses are an example of viruses that have a cleavage site motif that is generally associated with virulent viruses. However, most of them are considered to be of intermediate or low virulence for chickens, as assessed through the intracerebral pathogenicity index (ICPI) test [18, 19]. Nevertheless, their pathogenicity might be increased after passage in poultry species [20, 21].
Diagnostic testing and rapid detection of NDV are important steps in controlling an ND outbreak. Virus isolation in specific pathogen free (SPF) embryonated chicken eggs (ECE), followed by identification using hemagglutination (HA) and HI assays with a NDV-monospecific antiserum  is considered “gold standard” for NDV diagnostics. This approach is time-consuming and laborious and often requires up to ten days . Molecular diagnostics assays are a viable alternative to classical diagnostic assays and are widely used. Several protocols for the detection of NDV by reverse-transcription PCR (RT-PCR) have been published in the last decade . Real-time RT-PCR (rRT-PCR) is a rapid screening assay allowing for detection and pathotyping of NDV directly from diagnostic specimens and different protocols based on the use of hydrolysis probes, SybrGreen or LUX primers, have been published [24,25,26]. Real-time RT-PCR is highly dependent on genetic similarity between the primers/probe and the target genome, and protocols have to be updated as the genomes of pathogens accumulate mutations over time.
To improve our understanding of the distribution and evolution of NDV of genotype VI, viruses of this genotype circulating in five countries from 1982 to 2016 were isolated and studied. For this purpose, the following studies were done: i) isolation and biological and phylogenetic characterization of genotype VI viruses from different geographical locations; ii) analyses of the complete fusion protein gene coding sequences and complete genome sequences obtained from the studied viruses; iii) evaluation of their epidemiological relation to other circulating NDVs; and iv) optimization of rRT-PCR diagnostic test for detection of virulent variants of sub-genotype VIg as a result of identified failure of the current validated diagnostic protocol.
Sample collection and isolates background data
The samples and viruses studied here were collected from different species (pigeons, chickens and quail) in Egypt, Pakistan, South Korea, Ukraine and Bulgaria, representing three different continents (Africa, Asia and Europe). During April 2015, one hundred sixty-seven (n = 167) oropharyngeal and cloacal/fecal swabs were collected in Egypt from apparently healthy Rock Pigeons. Samples were collected from pigeon lofts in El Fayoum and Qena provinces, and a live bird market (LBM) in Cairo named Souq al-Goma’a (also named Souq Sayeda Aisha or Friday Market). The swabs were placed immediately into tubes with 3 ml of brain-heart-infusion broth (Difco, New Zealand) supplemented with penicillin G (10,000 IU /ml), amphotericin B (20 μg/ml), and gentamycin (1000 μg/ml). After collection, the samples were labelled, stored on ice, transported to the lab where kept frozen at −76 °C, and shipped on dry ice to the Southeast Poultry Research Laboratory (SEPRL) of the United States Department of Agriculture (USDA). Additionally, three more samples were collected from pigeons in Pakistan in 2015 and shipped to SEPRL, USDA. Six additional NDV from repositories in Bulgaria (n = 2), Ukraine (n = 2), and South Korea (n = 2), were also sent to SEPRL for further characterization. The Bulgarian and the Ukrainian viruses were passaged in eggs 3 times and 2 times, respectively. The Korean isolates were passaged 2 times at SEPRL and passage information was missing prior to their receiving. Nine more NDV, isolated between 2014 and 2016 from healthy and diseased pigeons (pet and zoo birds), were studied at the Quality Operations Laboratory (QOL) at the University of Veterinary and Animal Sciences (UVAS), Lahore, Pakistan.
Virus isolation, virus propagation and intracerebral pathogenicity index test
Initial screening of all Egyptian samples employing the NDV and avian influenza matrix gene rRT-PCR assays [24, 27], revealed that out of 167 oropharyngeal and cloacal swabs, 71 samples had cycle threshold (Ct) values ≤35 (40.3%) in the NDV assay and all samples were negative in the avian influenza assay. At SEPRL, thirty-one (selected to achieve representativeness of all locations) of the 71 rRT-PCR positive Egyptian samples and three Pakistani samples from 2015 (designated 21A, 22A and 25A) were selected for further studies and inoculated into 9-to-11-day-old SPF ECE, following standard procedures . The SPF ECE and chickens used to characterize these viruses were from the SEPRL SPF White Leghorn flock. The allantoic fluids from both eggs with embryo mortality and embryos alive at the end of the incubation period were collected and tested by hemagglutination assay. All the hemagglutinating agents were confirmed as NDV using a HI assay with NDV specific antiserum . Additionally, The NDV obtained from repositories in Ukraine, Bulgaria and South Korea were propagated into 9-to-11-day-old SPF ECE following the same procedures . Assessment of the virulence in vivo of three selected viruses was done by the ICPI test using one-day-old SPF chickens following established procedures .
Complete fusion protein gene sequencing
RNA from the two South Korean isolates was extracted from infected allantoic fluid using TRIzol LS Reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer’s instructions at SEPRL. RT-PCR was performed to amplify the complete F gene using the Superscript™ III One-step RT-PCR kit with Platinum Taq DNA polymerase (Invitrogen, Carlsbad, CA, USA), per manufacturer’s instructions using a set of F-gene specific primers (4331F/5090R, MSF1/NDVR2, 4927F/5673R and 5491F/6341R)  (see Additional file 1: Table S1). The PCR amplicons were processed and sequenced as described previously . RNA isolation and nucleotide sequencing of nine Pakistani viruses were performed in Quality operations laboratory at University of Veterinary and Animal Science in Pakistan as follows: RNA extraction and RT-PCR F-gene amplification were performed as described above using a set of F-gene specific primers (see Additional file 1: Table S1). The amplicons were electrophoresed using 1% agarose gel and purified using QIAquick® Gel Extraction Kit (Qiagen, Valencia, CA, USA). The purified products were sequenced using the ABI 3130 automated sequencer (Applied Biosystem Inc., Foster City, CA, USA), as described previously .
Complete genome sequencing using next-generation sequencing (NGS)
Total viral RNA of six Egyptian (positive in virus isolation assay), two Ukrainian, three Pakistani, and two Bulgarian viruses were extracted from the infected allantoic fluids using QIAmp® Viral RNA Mini Kit (Qiagen, USA) according to manufacturer’s instructions. The recovered RNA was quantified using the Qubit® RNA HS Assay Kit (Life Technologies, USA) in the Qubit® Fluorometer instrument (Invitrogen, USA). Newcastle disease virus RNA was captured and enriched using Sera-Mag beads (GE Healthcare Life Sciences, USA), and three biotin-labeled oligonucleotide probes targeting three different positions in the NDV genome: 1) 8-AGA GAA TCT GTG AGG TAC GA/3Bio/, 2) 5905-TTC TCA AGT CAT CGT GAC AG/3Bio/, and 3) 12226-CCC TGC ATC TCT CTA CAG/3Bio/. Reverse transcription was performed using the Moloney Murine Leukemia Virus Reverse Transcriptase (M-MLV RT) kit (Thermo Scientific, USA) according to the manufacturer’s instructions. The cDNA products were recovered and purified using the Agencourt® RNAClean® XP beads (Beckman Coulter, USA) according to the manufacturer’s instructions and quantified using Qubit® ssDNA Assay Kit (Thermo Fisher Scientific, USA) in the Qubit instrument. The purified cDNA products were tagmented and amplified for NGS by using 1 ng/5 μl (0.25 ng/μl in water) of the cDNA product employing the Nextera XT DNA Library Preparation Kit (Illumina, USA) following manufacturer’s protocol. The two Bulgarian samples were processed using the KAPA Stranded RNA-Seq Library Preparation Kit for Illumina platforms (Kapa Biosystems, USA) according to the manufacturer’s instructions.
The distribution size and concentration of DNA in the prepared libraries were checked on a Bioanalyzer 2100 and Qubit instrument using Agilent High Sensitivity DNA Kit (Agilent Technologies, Germany) and Qubit® dsDNA HS Assay Kit (Life Technologies, USA), respectively. Paired-end sequencing (2 × 250 base pairs) of the generated libraries was performed on an Illumina MiSeq instrument using the 500 cycle MiSeq Reagent Kit version 2 (Illumina, USA). Raw sequence data were analyzed and assembled using MIRA version 3.4.1  within a customized workflow on the Galaxy platform  as described previously . Short internal gaps at the 3′ UTR of the nucleoprotein gene were closed using Sanger technology with primers designed using the sequences obtained from NGS (see Additional file 1: Table S1). The 5′ and 3′ ends of the genomes reported here were sequenced and confirmed as described previously .
Collection of sequences and phylogenetic analyses
All available NDV complete genome and complete fusion protein gene coding sequences were downloaded from GenBank and curated, resulting in two large datasets (n = 331 and n = 1406, respectively). Together with the sequences obtained in the current study, all sequences were aligned using Multiple Alignment with Fast Fourier Transformation (MAFFT v7.017)  as implemented in the Geneious software v8.1.4 . The datasets were used for two preliminary maximum likelihood phylogenetic analyses using MEGA6 . Based on evolutionary relatedness to the viruses sequenced here, two smaller datasets (including newly obtained sequences) of complete fusion protein gene coding sequences (n = 82) (see Additional file 1: Table S2) and complete genome sequences (n = 83) (see Additional file 1: Table S3) were parsed from the larger datasets. Sequences of selected representative isolates from other NDV genotypes were also included in each of the datasets. The coding regions of the complete genome and complete fusion gene were used to construct final phylogenetic trees using MEGA6. To select best-fit substitution model, the Bayesian Information Criterion (BIC) and corrected Akaike Information Criterion (AICc) values were estimated using MEGA6. The General Time Reversible (GTR) model as implemented in MEGA6 with a discrete gamma distribution (4 categories [+G, parameter = 0.6519 for the full fusion gene tree and 0.5984 for the complete genome tree]) with 1000 bootstrap replicates was used in all data analysis. The phylogenetic trees were visualized and edited using tree explorer implemented in MEGA6 and branch lengths are proportional to the differences between the isolates. The evolutionary distances were inferred using MEGA6 and showed as the average number of base substitutions per site. Analyses were conducted using the maximum composite likelihood model  with a gamma distribution (shape parameter = 1) of rate variation among sites. For all analyses, the codon positions included were the 1st, 2nd, 3rd, and noncoding and positions containing gaps and missing data were eliminated from the datasets. The classification criteria proposed by Diel et al. for naming sub-genotypes and genotypes were followed in the current study .
Real-time reverse transcription polymerase chain reaction (rRT-PCR)
Allantoic fluids from all of the samples studied at SEPRL were subjected to the USDA-validated matrix gene rRT-PCR (M-gene assay), as well as to the USDA-validated fusion gene rRT-PCR assay (F-gene assay) described previously by Wise et al. . An additional test was performed using pigeon-specific fusion protein gene assay as described previously by Kim et al. . AgPath one-Step rRT-PCR Reagents (Thermo Scientific, USA) was used and the reactions were carried out in the Cepheid Real-Time Thermal Cycler (Life Science, USA). The samples that had a Ct value ≤35 were considered positive in all assays. Due to the failure of the validated F-gene rRT-PCR assay  and pigeon-specific fusion protein gene assay  to detect some of the genotype VI viruses studied here, the probe and the forward primer described by Kim et al.,  were optimized as follows: F-4876 5′–[6-FAM] AAG CGY TTC TGT CTC YTT CCT CCT [BHQ_1]–3′ and (F + 4837) 5′- TGA TTC CAT CCG CAG GAT ACA AG -3′. Additionally, the reverse primer for the pigeon-specific F-gene assay (F-4837) was replaced with a new primer (F-4943) 5′- GCT GCT GTT ATC TGT GCC GA-3′. The optimized probe and newly designed primers were analyzed by OligoAnalyzer 3.1. (Integrated DNA Technologies, USA, https://eu.idtdna.com/calc/analyzer ) tool and checked for self-annealing, hairpin loops and heterodimers. The sequences of the optimized primers and probes along with the tested fifteen NDV and selected representatives from other genotypes were aligned and compared to the M-gene and F-gene primers and probes and the pigeon-specific fusion probe to identify variable sites that determined the different outcomes of the F-gene rRT-PCR assay [24, 39] (see Additional file 2: Figs. S1, S2, and S3).
Virus isolation and biological properties
Twenty-four NDV isolates from pigeons, chickens and quail, isolated in five different countries were studied here (Table 1). Most samples were obtained from pigeons (pet pigeons living in lofts, n = 12; zoo pigeons, n = 2; live bird markets, n = 2; and pigeons of unknown habitats, n = 6). Out of 31 oropharyngeal and cloacal Egyptian swabs from healthy birds inoculated in eggs, six contained live NDV and were fully characterized. The three Pakistani clinical samples were found to be positive for NDV by virus isolation and HI. All viruses obtained from repositories were successfully re-isolated after propagation and confirmed as NDV by HI. The ICPI values of selected viruses (pigeon/Egypt/Giza/11/2015, pigeon/Ukraine/Doneck/3/2007 and pigeon/Pakistan/Lahore/25A/2015) were 1.31, 1.48 and 1.46, respectively, and these indexes are typical for NDV with moderate virulence (mesogenic) in chickens .
Complete fusion protein gene and complete genome sequencing
The complete genome sequences of thirteen of the studied viruses were obtained at SEPRL, and the genome characteristics of these complete genomes are typical for NDV (See Additional file 1: Table S4). Additional eleven complete fusion gene coding sequences were also obtained and utilized for prediction of virulence and phylogenetic analysis. The fusion protein cleavage site of NDV is a major virulence determinant marker , and the deduced amino acid sequences of the fusion protein cleavage site revealed that all 24 isolates had multiple basic amino acids at positions 113–116, and a phenylalanine residue at position 117 (Table 1). The cleavage site motif of the Egyptian and Ukrainian viruses was 113RQKR↓F117 (n = 8). The Pakistani viruses (n = 12) had the following three motifs, 113RKKRF117, 113RQRR↓F117, and 113RQKR↓F117, sharing the last motif with the Bulgarian viruses (n = 2). The South Korean (n = 2) viruses had a 13RRKR↓F117 motif (Table 1). All of these motifs are characteristic of virulent NDV [3, 41].
Distance and phylogenetic analyses
In order to determine the phylogenetic relationship between the studied viruses and other NDV isolated worldwide, the complete fusion gene coding sequences obtained in the current work, along with sequences of highly related viruses, were used to construct a phylogenetic tree (n = 82) (Fig. 1 and Additional file 1: Table S2). The full fusion coding region of all available class II NDV and of genotype VI (n = 1430 and n = 281, respectively, both including the 24 sequences from the current study) was also used to estimate evolutionary distances (Tables 2 and 3). A second phylogenetic analysis was performed using 83 complete genome-concatenated coding sequences of viruses pertaining to class II (see Additional file 1: Table S3 and Additional file 2: Fig. S4). Both, the full fusion gene and the complete genome phylogenetic analyses displayed similar topology confirming the phylogenetic classification of the viruses studied here into different sub-genotypes related to viruses previously isolated in Asia, Africa and Eastern Europe. The phylogenetic analyses (Fig. 1 and Additional file 2: Fig. S4) demonstrated that these 24 NDV studied here the clustered with viruses of genotype VI. The isolates from pigeons grouped with previously described viruses from columbine birds, while the chicken isolates grouped with other viruses from chickens. The topology of the full fusion phylogenetic tree (Fig. 1) indicated the viruses from Egypt (n = 6) (0%–3.1% genetic distance among themselves) and Ukraine (n = 2) (1.9% genetic distance between themselves) isolated from pigeons grouped together and clustered within sub-genotype VIg. These viruses were closely related to viruses isolated from pigeons in Russia, Ukraine, Kazakhstan, and Nigeria during 2005–2014 (Fig. 1) [42,43,44]. These results were consistent with the results of the complete genome analysis (see Fig. S4). Four Pakistani viruses isolated during 2015–2016 clustered together in a separate monophyletic branch within sub-genotype VIg (Fig. 1). The genetic distance of these Pakistani viruses (4.7%) to the rest of the VIg viruses shows higher nucleotide diversity within the sub-genotype.
The remaining eight Pakistani viruses (0%–3.6% genetic distance among themselves, see Additional file 1: Table S5) isolated from pigeons during 2014–2016 did not cluster within any of the previously known sub-genotypes within genotype VI. These Pakistani isolates grouped with a Pakistani virus isolated from pigeon (KU885949/pigeon/Pakistan/MZS-UVAS/2014) creating a separate branch in the phylogenetic tree (Fig. 1) and were 8.5% to 12% distant from the rest of the sub-genotypes in genotype VI (Table 2). This new group of viruses fulfills all classification criteria set by Diel et al.  and was named as a novel, hitherto undescribed, sub-genotype of class II genotype VI, namely sub-genotype VIm. The viruses from the newly designated sub-genotype VIm were also more than 10% (11.2–24.2%) distant from all known NDV class II genotypes (Table 3).
The topology of the phylogenetic tree based on the full fusion gene sequences revealed that the Bulgarian virus (Chicken/Bulgaria/Mokresh/1982) grouped together with isolates from Ethiopia (during 2011–2012)  (Fig.1). However, the Bulgarian isolate was quite divergent from other members in this sub-genotype showing high genetic distance (10–11%). When this Bulgarian isolate was phylogenetically analyzed based on complete genome, it shared common ancestor with sub-genotypes VIg and VIm described above (see Additional file 2: Fig. S4). The other Bulgarian virus (chicken/Bulgaria/Dolno Linevo/1992), as well as both South Korean viruses obtained from a chicken and a quail, grouped within sub-genotype VIc with other chicken isolates from Sweden, Japan and China indicating that they belong to the older genotype VI sub-genotype (Fig. 1).
Rapid diagnostic and development of new primers and probe
While all tested viruses were positive with the M-gene rRT-PCR assay , the F-gene rRT-PCR assay used to identify virulent NDV  failed to detect eight viruses from Egypt (n = 6) and Ukraine (n = 2) (Table 4). The pigeon-specific F-gene assay  failed to detect all fifteen NDV (the viruses sequenced in Pakistan were not tested, designated in Table 1). The genomic sequences of the viruses that failed detection were used to design new primers and probe based on a previous assay . The new pigeon-specific F-gene test successfully detected all the genotype VI NDV (n = 15) that were not detected by the previously available protocol (Table 4). Upon further analysis, the new pigeon-specific fusion gene rRT-PCR assay was evaluated utilizing total RNA from allantoic fluids infected with virulent viruses of different genotypes and also side by side comparison of equal amount of total RNAs were performed with the M-gene and F-gene assays . The results revealed that while the new fusion test recognized most effectively viruses of sub-genotype VIg, although with higher Ct values, it also recognized virulent viruses of seven different genotypes (II, Vb, VIa, VIIi, XII, XIVb, and XVIIa) (Table 4). While three Pakistani viruses (pigeon/Pakistan/Lahore/21A/2015, pigeon/Pakistan/Lahore/22A/2015 and pigeon/Pakistan/Lahore/25A/2015) were detected by the regular F-gene assay, using the new set of primers and probe resulted in four to nine lower Ct values (Table 4).
Here, we present the relationship between NDV of genotype VI isolated from Egypt, Ukraine, Pakistan, Bulgaria and South Korea based on their full fusion and complete genome characterization. Our data confirm the concurrent evolution and mobility of viruses of two sub-genotypes of genotype VI NDV across 3 continents. Viruses of genotype VI were first isolated from pigeons in the Middle East in 1960s and spread rapidly throughout Northern Africa to Europe and the rest of the world [45,46,47]. In Egypt, NDV was first identified in 1947  on the basis of virus isolation into ECE and serologically by HI tests. Clinical signs in pigeons consistent with ND have been seen in Egypt since early 1981, and infection with NDV was serologically confirmed in diseased pigeons in the delta area of the Nile in 1984 [49, 50]. In Bulgaria, ND was first detected in 1943 , however the viruses of genotype VI were first found in Bulgaria in the mid-1970s . To the best of our knowledge, the first identification of genotype VI NDV in Ukraine, Pakistan and South Korea has not been previously documented. The phylogenetic relationship among Eastern European, African and Asian viruses suggest the circulation of related viruses in pigeons across three continents.
Pigeons are not migratory and the circulation of these closely related viruses isolated from pigeons in six distant countries (Fig. 1) within eight years (2007–2015) is epidemiologically important. In many countries, pigeons (Columba livia) live freely as synantropic birds, or are bred for a variety of different purposes, such as a source of meat, pet companion birds, or for laboratory experiments in biology and cognitive science. Viruses of genotype VI have been reported previously to circulate in apparently healthy pigeons [52,53,54], and in the current work it was confirmed that at least sub-genotype VIg is seemingly maintained in healthy pigeons kept in captivity in Egypt (all samples from Egypt were collected from healthy pigeons). However, the mechanism of spread of these viruses at long distances remains unknown. A possible explanation for the spread of genotype VI NDV is the contact between columbid birds during competition flights, exhibitions, or due to the intensive international trade of such birds . Other possibility is the international trade of live birds from other species or avian products between countries, either by legal or illegal routes of importation and exportation. In addition, NDV of genotype VI have also been isolated from birds from non-columbidae species kept in captivity and from wild birds, including partridges, pheasants, swans, falcons, blackbirds, cockatoos, budgerigars, raptors, partridges, crested ibises, waterfowl, starlings, pintails, gannets, and buzzards [44, 56,57,58,59]. However, while other NDV (e.g. genotype VII) are known to infect wild birds and this possibility can’t be excluded, no evidence exists that wild birds, especially migratory, play a role in the spread of genotype VI viruses. The presence of genotype VI viruses in fecal and oral swabs suggests that viral replication, which could result in virus transmission and possible outbreaks in poultry, as seen previously, is occurring in pigeons. Although sporadic, ND outbreaks in poultry caused by viruses of genotype VI have been reported [60,61,62,63], and the potential of the virus to cause clinical disease in poultry must not be under estimated.
The phylogenetic analyses revealed the complexity of NDV genotype VI and the challenges in the classification of its sub-genotypes. As shown in Table 2, some sub-genotypes have more than 10% distance compared to the rest of the sub-genotypes within genotype VI. In addition, the newly designated sub-genotype VIm is more than 10% distant from the rest of the genotypes in NDV class II (Table 3). These viruses topologically fall into a lower-order group (see Fig. 1 and Additional file 2: Fig. S4) compared to the existing genotype VI. Some of these groups of viruses meet, or will eventually meet the classification criteria for consideration as new genotypes. However, as naming new genotypes that originated from existing sub-genotypes may create confusion in the field of NDV classification, here, the most diverse group of viruses in genotype VI is named as sub-genotype VIm. In the future it may be appropriate to utilize rules similar to those put forth the WHO/OIE/FAO H5N1 Evolution Working Group for the nomenclature of highly pathogenic H5N1 avian influenza viruses  as already proposed by Susta et al. . According to those rules, when a group is split into subgroups of higher order, the newly named groups remain part of the existing original lower-order group (e.g., VIa.1 or VIi.1 and so on). We believe that the criteria for naming new genotypes need to be updated; however, this has to be done based on international consensus rather than by individual scientific teams .
The failure of the validated F-gene rRT-PCT assay  designed to specifically detect viruses from the outbreak that occurred in southern California in 2002–2003 to detect the sub-genotype VIg Egyptian and the Ukrainian isolates (Table 4) was attributed to the mismatches in the probe and/or the primers. Here, several new variable sites that resulted in nucleotide mismatches were identified, suggesting continuous variation at the site used for the design of the test (see Additional file 2: Figsure S1, S2, and S3). The pigeon-specific probe, designed to detect dove/Italy/2736/2000 and US pigeon viruses  also had significant mismatches to the viruses tested here. The increasing number of mismatches of different probes and primers indicates that genotype VI is composed of a highly diverse group of viruses that is not covered by a single test. The newly developed pigeon test successfully identified the genotype VI isolates tested here (including the sub-genotype VIg); however, the new fusion test is not genotype specific, and positive detection of viruses of other genotypes is also possible (although with decreased sensitivity). Albeit also dependent on specific primers and not routinely performed in all laboratories, sequencing of the fusion gene will provide data allowing definitive diagnostics and characterization.
In summary, genotype VI NDV continue to evolve in Africa, Asia and Europe suggesting the need for a constant surveillance and characterization of these viruses. The described epidemiological connections among viruses isolated from non-migratory birds on three continents creates uncertainty regarding the mechanism of spread of these viruses. Further studies are needed to elucidate this issue and understand the possible role of human activity in the dispersal of these viruses. Complete genomic characterization identified previously unrecognized genetic diversity that contributes to diagnostic failure and will facilitate future evolutionary studies. These results highlight the importance of updating existent rapid diagnostic protocols to detect emerging viral variants and help manage the disease in affected regions. The obtained rRT-PCR results suggest that the developed pigeon-specific F-gene assay, run in conjunction with the USDA-validated M-gene rRT-PCR assay, should be effective in detecting viruses from sub-genotype VIg, and those from the newly classified sub-genotype VIm.
Aldous EW, Alexander DJ. Detection and differentiation of Newcastle disease virus avian paramyxovirus type 1. Avian Pathol. 2001;30:117–28. doi: 10.1080/03079450120044515
Miller PJ, Decanini EL, Afonso CL. Newcastle disease: Evolution of genotypes and the related diagnostic challenges. Infect Genet Evol. 2010;10:26–35. doi: 10.1016/j.meegid.2009.09.012
OIE. Newcastle disease. In Biological Standards Commission, Manual of diagnostic tests and vaccines for terrestrial animals: mammals, birds and bees, vol. Volume 1, Part 2, Chapter 2.3.14, 7th edition. pp. 555–74. Paris, France: World Organisation for Animal Health; 2012:555–574.
Afonso CL, Amarasinghe GK, Banyai K, Bao Y, Basler CF, Bavari S, Bejerman N, Blasdell KR, Briand FX, Briese T, et al. Taxonomy of the order Mononegavirales: update 2016. Arch Virol. 2016;161:2351–60. doi: 10.1007/s00705-016-2880-1
Goraichuk I, Sharma P, Stegniy B, Muzyka D, Pantin-Jackwood MJ, Gerilovych A, Solodiankin O, Bolotin V, Miller PJ, Dimitrov KM, Afonso CL. Complete Genome Sequence of an Avian Paramyxovirus Representative of Putative New Serotype 13. Genome Announc. 2016;4:e00729–16. doi: 10.1128/genomeA.00729-16
Thampaisarn R, Bui VN, Trinh DQ, Nagai M, Mizutani T, Omatsu T, Katayama Y, Gronsang D, Le DH, Ogawa H, Imai K. Characterization of avian paramyxovirus serotype 14, a novel serotype, isolated from a duck fecal sample in Japan. Virus Res. 2016;228:46–57. doi: 10.1016/j.virusres.2016.11.018
Thomazelli LM, Araujo J, Oliveira DB, Sanfilippo L, Ferreira CS, Brentano L, Pelizari VH, Nakayama C, Duarte R, Hurtado R, et al. Newcastle disease virus in penguins from King George Island on the Antarctic region. Vet Microbiol. 2010; doi: 10.1016/j.vetmic.2010.05.006
Lee HJ, Kim JY, Lee YJ, Lee EK, Song BM, Lee HS, Choi KS. A Novel Avian Paramyxovirus (Putative Serotype 15) Isolated from Wild Birds. Front Microbiol. 2017;8:786. doi: 10.3389/fmicb.2017.00786
Neira V, Tapia R, Verdugo C, Barriga G, Mor S, Ng TFF, Garcia V, Del Rio J, Rodrigues P, Briceno C, et al. Novel Avulaviruses in Penguins. Antarctica Emerg Infectious Diseases. 2017;23:1212–4. doi: 10.3201/eid2307.170054
Seal BS, King DJ, Sellers HS. The Avian Response to Newcastle Disease Virus. Dev Comp Immunol. 2000;24:257–68. doi: 10.1016/S0145-305X(99)00077-4
Czeglédi A, Ujvári D, Somogyi E, Wehmann E, Werner O, Lomniczi B. Third genome size category of avian paramyxovirus serotype 1 (Newcastle disease virus) and evolutionary implications. Virus Res. 2006;120:36–48. doi: 10.1016/j.virusres.2005.11.009
Diel DG, da Silva LH, Liu H, Wang Z, Miller PJ, Afonso CL. Genetic diversity of avian paramyxovirus type 1: proposal for a unified nomenclature and classification system of Newcastle disease virus genotypes. Infect Genet Evol. 2012;12:1770–9. doi: 10.1016/j.meegid.2012.07.012
Dimitrov KM, Ramey AM, Qiu X, Bahl J, Afonso CL. Temporal, geographic, and host distribution of avian paramyxovirus 1 Newcastle disease virus. Infect Genet Evol. 2016;39:22–34. doi: 10.1016/j.meegid.2016.01.008
Terregino C, Cattoli G, Grossele B, Bertoli E, Tisato E, Capua I. Characterization of Newcastle disease virus isolates obtained from Eurasian collared doves (Streptopelia decaocto) in Italy. Avian Pathol. 2003;32:63–8. doi: 10.1080/0307945021000070732
Collins MS, Alexander DJ, Brockman S, Kemp PA, Manvell RJ. Evaluation of mouse monoclonal antibodies raised against an isolate of the variant avian paramyxovirus type 1 responsible for the current panzootic in pigeons. Arch Virol. 1989;104:53–62. doi: 10.1007/BF01313807
Ujvari D, Wehmann E, Kaleta EF, Werner O, Savic V, Nagy E, Czifra G, Lomniczi B. Phylogenetic analysis reveals extensive evolution of avian paramyxovirus type 1 strains of pigeons (Columba livia) and suggests multiple species transmission. Virus Res. 2003;96:63–73. doi: 10.1016/S0168-1702(03)00173-4
Xue C, Xu X, Yin R, Qian J, Sun Y, Wang C, Ding C, Yu S, Hu S, Liu X, et al. Identification and pathotypical analysis of a novel VIk sub-genotype Newcastle disease virus obtained from pigeon in China. Virus Res. 2017;238:1–7. doi: 10.1016/j.virusres.2017.05.011
Dortmans JC, Fuller CM, Aldous EW, Rottier PJ, Peeters BP. Two genetically closely related pigeon paramyxovirus type 1 (PPMV-1) variants with identical velogenic fusion protein cleavage sites but with strongly contrasting virulence. Vet Microbiol. 2009; doi: 10.1016/j.vetmic.2009.11.021.
Meulemans G, van den Berg TP, Decaesstecker M, Boschmans M. Evolution of pigeon Newcastle disease virus strains. Avian Pathol. 2002;31:515–9. doi: 10.1080/0307945021000005897
Collins MS, Strong I, Alexander DJ. Evaluation of the molecular basis of pathogenicity of the variant Newcastle disease viruses termed “pigeon PMV-1 viruses”. Arch Virol. 1994;134:403–11. doi: 10.1007/BF01310577
Fuller CM, Collins MS, Easton AJ, Alexander DJ. Partial characterisation of five cloned viruses differing in pathogenicity, obtained from a single isolate of pigeon paramyxovirus type 1 (PPMV-1) following passage in fowls' eggs. Arch Virol. 2007;152:1575–182. doi: 10.1007/s00705-007-0963-8
Dimitrov KM, Clavijo A, Sneed L. RNA extraction for molecular detection of Newcastle disease virus – comparative study of three methods. Rev Méd Vét. 2014;165:172–5.
Cattoli G, Monne I. Molecular Diagnosis of Newcastle Disease Virus. In: Capua I, Alexander DJ, editors. . Springer Milan: Avian Influenza and Newcastle Disease: A Field and Laboratory Manual; 2009. p. 127–32.
Wise MG, Suarez DL, Seal BS, Pedersen JC, Senne DA, King DJ, Kapczynski DR, Spackman E. Development of a real-time reverse-transcription PCR for detection of newcastle disease virus RNA in clinical samples. J Clin Microbiol. 2004;42:329–38. doi: 10.1128/JCM.42.1.329-338.2004
Pham HM, Konnai S, Usui T, Chang KS, Murata S, Mase M, Ohashi K, Onuma M. Rapid detection and differentiation of Newcastle disease virus by real-time pcr with melting-curve analysis. Arch Virol. 2005;150:2429–38. doi: 10.1007/s00705-005-0603-0
Antal M, Farkas T, German P, Belak S, Kiss I. Real-time reverse transcription-polymerase chain reaction detection of Newcastle disease virus using light upon extension fluorogenic primers. J Vet Diagn Investig. 2007;19:400–4. 17609351.
Spackman E, Senne DA, Myers TJ, Bulaga LL, Garber LP, Perdue ML, Lohman K, Daum LT, Suarez DL. Development of a real-time reverse transcriptase PCR assay for type A influenza virus and the avian H5 and H7 hemagglutinin subtypes. J Clin Microbiol. 2002;40:3256–60. doi: 10.1128/JCM.40.9.3256-3260.2002
Alexander DJ, Swayne DE. Newcastle disease virus and other avian paramyxoviruses. In: Swayne DE, Glisson JR, Jackwood MW, Pearson JE, Reed WM, editors. A Laboratory Manual for the Isolation and Identification of Avian Pathogens. 4th ed. Kennett Square, PA: The American Association of Avian Pathologists; 1998. p. 156–63.
Miller PJ, Dimitrov KM, Williams-Coplin D, Peterson MP, Pantin-Jackwood MJ, Swayne DE, Suarez DL, Afonso CL. International biological engagement programs facilitate Newcastle disease epidemiological studies. Front Public Health. 2015;3:235. doi: 10.3389/fpubh.2015.00235
Miller PJ, Haddas R, Simanov L, Lublin A, Rehmani SF, Wajid A, Bibi T, Khan TA, Yaqub T, Setiyaningsih S, Afonso CL. Identification of new sub-genotypes of virulent Newcastle disease virus with potential panzootic features. Infect Genet Evol. 2015;29:216–29. doi: 10.1016/j.meegid.2014.10.032
Chevreux B, Wetter T, Suhai S. Genome sequence assembly using trace signals and additional sequence information. In: Computer Science and Biology, vol. 99. Hanover, Germany; 1999. p. 45–56.
Afgan E, Baker D, Van den Beek M, Blankenberg D, Bouvier D, Čech M, Chilton J, Clements D, Coraor N, Eberhard C. The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2016 update. Nucleic Acids Res. 2016:gkw343. doi: 10.1093/nar/gkw343
Dimitrov KM, Sharma P, Volkening JD, Goraichuk IV, Wajid A, Rehmani SF, Basharat A, Shittu I, Joannis TM, Miller PJ, Afonso CL. A robust and cost-effective approach to sequence and analyze complete genomes of small RNA viruses. Virol J. 2017;14:72. doi: 10.1186/s12985-017-0741-5
Brown PA, Briand F-X, Guionie O, Lemaitre E, Courtillon C, Henry A, Jestin V, Eterradossi N. An alternative method to determine the 5′ extremities of non-segmented, negative sense RNA viral genomes using positive replication intermediate 3′ tailing: Application to two members of the Paramyxoviridae family. J Virol Methods. 2013;193:121–7. doi: 10.1016/j.jviromet.2013.05.007
Katoh K, Misawa K, Kuma K, Miyata T. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002;30:3059–66. doi: 10.1093/nar/gkf436
Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, Buxton S, Cooper A, Markowitz S, Duran C, et al. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28:1647–9. doi: 10.1093/bioinformatics/bts199
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9. doi: 10.1093/molbev/mst197
Tamura K, Nei M, Kumar S. Prospects for inferring very large phylogenies by using the neighbor-joining method. Proc Natl Acad Sci U S A. 2004;101:11030–5. doi: 10.1073/pnas.0404206101
Kim LM, Afonso CL, Suarez DL. Effect of probe-site mismatches on detection of virulent Newcastle disease viruses using a fusion-gene real-time reverse transcription polymerase chain reaction test. J Vet Diagn Investig. 2006;18:519–28. 17121078.
Glickman RL, Syddall RJ, Iorio RM, Sheehan JP, Bratt MA. Quantitative Basic Residue Requirements in the Cleavage-Activation Site of the Fusion Glycoprotein as a Determinant of Virulence for Newcastle Disease Virus. J Virol. 1988;62:354–6. 3275436.
Miller PJ, Koch G. Newcastle disease. In: Swayne DE, Glisson JR, McDougald LR, Nolan LK, Suarez DL, Nair V, editors. Diseases of Poultry. 13th ed. Hoboken, New Jersey: Wiley-Blackwell; 2013. p. 89–138.
Pchelkina IP, Manin TB, Kolosov SN, Starov SK, Andriyasov AV, Chvala IA, Drygin VV, Yu Q, Miller PJ, Suarez DL. Characteristics of pigeon paramyxovirus serotype-1 isolates (PPMV-1) from the Russian Federation from 2001 to 2009. Avian Dis. 2013;57:2–7. doi: 10.1637/10246-051112-Reg.1
Yurchenko KS, Sivay MV, Glushchenko AV, Alkhovsky SV, Shchetinin AM, Shchelkanov MY, Shestopalov AM. Complete Genome Sequence of a Newcastle Disease Virus Isolated from a Rock Dove (Columba livia) in the Russian Federation. Genome Announc. 2015:3e01514–4. doi: 10.1128/genomeA.01514-14
Benson DA, Clark K, Karsch-Mizrachi I, Lipman DJ, Ostell J, Sayers EW. GenBank. Nucleic Acids Res. 2015;43:D30. doi: 10.1093/nar/gku1216
Czeglédi A, Herczeg J, Hadjiev G, Doumanova L, Wehmann E, Lomniczi B. The occurrence of five major Newcastle disease virus genotypes (II, IV, V, VI and VIIb) in Bulgaria between 1959 and 1996. Epidemiol Infect. 2002;129:679–88. doi: 10.1017/S0950268802007732
Lomniczi B, Wehmann E, Herczeg J, Ballagi-Pordany A, Kaleta EF, Werner O, Meulemans G, Jorgensen PH, Mante AP, Gielkens AL, et al. Newcastle disease outbreaks in recent years in western Europe were caused by an old (VI) and a novel genotype (VII). Arch Virol. 1998;143:49–64. doi: 10.1007/s007050050267
Alexander DJ, Russell PH, Parsons G, Elzein EMEA, Ballouh A, Cernik K, Engstrom B, Fevereiro M, Fleury HJA, Guittet M, et al. Antigenic and biological characterisation of avian paramyxovirus type I isolates from pigeons an international collaborative study. Avian Pathol. 1985;14:365–76. doi: 10.1080/03079458508436238
Daubney R, Mansy W. The occurrence of Newcastle disease in Egypt. J Comp Pathol Ther. 1948;58:189–200. doi: 10.1016/S0368-1742(48)80019-6
Kaleta E, Alexander D, Russell P. The first isolation of the avian PMV-1 virus responsible for the current panzootic in pigeons. Avian Pathol. 1985;14:553–7. doi: 10.1080/03079458508436258
Vindevogel H, Duchatel JP. Panzootic Newcastle disease virus in pigeons. In: Alexander DJ, editor. Newcastle disease. Boston: Kluwer Academic Publishers; 1988. p. 184–96.
Semerdzhiev B. Fowl plague in Bulgaria. Vet Sbirka. 1946;49:105–17.
Wang J, Liu H, Liu W, Zheng D, Zhao Y, Li Y, Wang Y, Ge S, Lv Y, Zuo Y, et al. Genomic Characterizations of Six Pigeon Paramyxovirus Type 1 Viruses Isolated from Live Bird Markets in China during 2011 to 2013. PLoS One. 2015;10:e0124261. doi: 10.1371/journal.pone.0124261
Byarugaba DK, Mugimba KK, Omony JB, Okitwi M, Wanyana A, Otim MO, Kirunda H, Nakavuma JL, Teillaud A, Paul MC, Ducatez MF. High pathogenicity and low genetic evolution of avian paramyxovirus type I (Newcastle disease virus) isolated from live bird markets in Uganda. Virol J. 2014;11:173. doi: 10.1186/1743-422X-11-173
Teske L, Ryll M, Rautenschlein S. Epidemiological investigations on the role of clinically healthy racing pigeons as a reservoir for avian paramyxovirus-1 and avian influenza virus. Avian Pathol. 2013;42:557–65. doi: 10.1080/03079457.2013.852157
Aldous EW, Fuller CM, Ridgeon JH, Irvine RM, Alexander DJ, Brown IH. The evolution of pigeon paramyxovirus type 1 (PPMV-1) in Great Britain: a molecular epidemiological study. Transbound Emerg Dis. 2014;61:134–9. doi: 10.1111/tbed.12006
Alexander DJ. Newcastle Disease and other Avian Paramyxoviridae Infections. In: Calnek BW, Barnes HJ, Beard CW, McDougald LR, Saif YM, editors. Diseases of Poultry. 10th ed. Ames, Iowa: Iowa State University Press; 1997. p. 541–69.
Aldous EW, Fuller CM, Mynn JK, Alexander DJ. A molecular epidemiological investigation of isolates of the variant avian paramyxovirus type 1 virus (PPMV-1) responsible for the 1978 to present panzootic in pigeons. Avian Pathol. 2004;33:258–69. https://doi.org/10.1080/0307945042000195768
Irvine RM, Aldous EW, Manvell RJ, Cox WJ, Ceeraz V, Fuller CM, Wood AM, Milne JC, Wilson M, Hepple RG, et al. Outbreak of Newcastle disease due to pigeon paramyxovirus type 1 in grey partridges (Perdix perdix) in Scotland in October 2006. Vet Rec. 2009;165:531–5. doi: 10.1136/vr.165.18.531
Krapez U, Steyer AF, Slavec B, Barlic-Maganja D, Dovc A, Racnik J, Rojs OZ. Molecular characterization of avian paramyxovirus type 1 (Newcastle disease) viruses isolated from pigeons between 2000 and 2008 in Slovenia. Avian Dis. 2010;54:1075–80. doi: 10.1637/9161-111709-ResNote.1
Alexander DJ, Parsons G, Marshall R. Infection of fowls with Newcastle disease virus by food contaminated with pigeon faeces. Vet Rec. 1984;115:601–2. 6523693.
Alexander DJ. Newcastle disease in the European Union 2000 to 2009. Avian Pathol. 2011;40:547–58. https://doi.org/10.1080/03079457.2011.618823
Alexander DJ, Wilson GW, Russell PH, Lister SA, Parsons G. Newcastle disease outbreaks in fowl in Great Britain during 1984. Vet Rec. 1985;117:429–34. doi: 10.1136/vr.117.17.429
Abolnik C, Gerdes GH, Kitching J, Swanepoel S, Romito M, Bisschop SP. Characterization of pigeon paramyxoviruses (Newcastle disease virus) isolated in South Africa from 2001 to 2006. Onderstepoort J Vet Res. 2008;75:147–52.
WHO. Continued evolution of highly pathogenic avian influenza A (H5N1): updated nomenclature. Influenza Other Respir Viruses. 2012;6:1–5. doi: 10.1111/j.1750-2659.2011.00298.x
Susta L, Dimitrov KM, Miller PJ, Brown CC, Afonso CL. Reply to “May Newly Defined Subgenotypes Va and Vb of Newcastle Disease Virus in Poultry Be Considered Two Different Genotypes?”. J Clin Microbiol. 2016;54:2205–6. doi: 10.1128/jcm.00914-16
Kim LM, King DJ, Suarez DL, Wong CW, Afonso CL. Characterization of class I Newcastle disease virus isolates from Hong Kong live bird markets and detection using real-time reverse transcription-PCR. J Clin Microbiol. 2007;45:1310–4. doi: 10.1128/JCM.02594-06
The authors gratefully acknowledge Tim Olivier for his technical assistance and Dr. Hassan El-Sayed Helal for his help in collecting Egyptian clinical samples. The mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U.S. Department of Agriculture. The USDA is an equal opportunity provider and employer.
This work was supported by the U.S. Department of Agriculture ARS CRIS 6040–32,000-072, the U.S. Department of State BEP/CRDF NDV 31063, and by the Egyptian Cultural and Education Bureau.
Availability of data and materials
The full-length fusion gene coding sequences and complete genome sequences (n = 11 and n = 13, respectively) of NDV obtained in this study were submitted to GenBank and are available under the accession numbers KU862297 to KU862299, KY042125 to KY042143, and KX236100 to KX236101.
CLA, PJM, MS, and KMD conceived this project. MS, AW, AB, SFR, DWC and DVM coordinated field sampling efforts and isolated viruses. MS, PS, IVG, DWC and KMD prepared and sequenced the viral isolates. MS, DWC, KMD, and CLA optimized the diagnostic tests. MS, KMD, PJM and CLA conducted analyses of the data. MS, KMD, PJM, and CLA wrote and edited the manuscript. All authors read and approved the final manuscript.
All experiments were conducted complying with protocols reviewed and approved by the SEPRL institutional biosafety committee and were conducted with appropriate measures to maintain biosecurity and biosafety. General care of chickens was provided in accordance with the procedures reviewed and approved by the SEPRL Institutional Animal Care and Use Committee, as outlined in the Guide for the Care and Use of Agricultural Animals in Agricultural Research and Teaching.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Nucleotide sequences of primers used in PCR amplification, and sequencing of the NDV isolates used in this study. Table S2. List of the NDV used for construction of full fusion phylogenetic tree presented in Fig. 1. Highlighted in bold font are the viruses studied in the current work. Table S3. List of the NDV used for construction of complete genome phylogenetic tree presented in Additional file 2: Fig. S4. Highlighted in bold font are the viruses studied in the current work. Table S4. Characteristics of the thirteen complete genomes of Newcastle disease viruses of genotype VI sequenced in this study. Table S5. Estimated pairwise evolutionary distances among viruses of the new sub-genotype VIm. (DOCX 60 kb)
Mismatches between the tested viruses and: A) previous fusion probe designed by Wise et al. ; B) pigeon-specific fusion probe designed by Kim et al. ; and C) the optimized pigeon specific probe in this study, respectively. Sequences are in order of 5′ to 3′. Figure S2. A and B Mismatches between the tested viruses and: A) previous fusion forward primer designed by Wise et al. ; and B) the optimized fusion forward primer in this study, respectively. Sequences are in order of 5′ to 3′. Figure S3. A and B Mismatches between the tested viruses and: A) previous fusion reverse primers designed by Wise et al. ; B) and the new fusion reverse in this study, respectively. Sequences are in order of 5′ to 3′. Figure S4. Phylogenetic analysis based on the complete genome concatenated coding sequence of viruses representing NDV class II. Only bootstrap values greater or equal to 60% are visualized. There were a total of 13,697 positions in the final dataset. The strains sequenced in this study are highlighted in bold font and have a circle symbol in front the taxa name. Provisional designation of genotypes is indicated on the right. (DOCX 1990 kb)
About this article
Cite this article
Sabra, M., Dimitrov, K.M., Goraichuk, I.V. et al. Phylogenetic assessment reveals continuous evolution and circulation of pigeon-derived virulent avian avulaviruses 1 in Eastern Europe, Asia, and Africa. BMC Vet Res 13, 291 (2017). https://doi.org/10.1186/s12917-017-1211-4