|
| |
 |
|
| ORIGINAL ARTICLE |
|
| Year : 2021 | Volume
: 14
| Issue : 2 | Page : 64-72 |
|
Genomic characterization of velogenic avian orthoavulavirus 1 isolates from poultry workers: Implications to emergence and its zoonotic potential towards public health
Muhammad Zubair Shabbir1, Ruth Helmus Nissly2, Abdul Ahad3, Masood Rabbani1, Levina Lim2, Shubhada K Chothe2, Murugan Subbiah2, Aswathy Sebastian4, Istvan Albert4, Aziz Ul-Rahman5, Bhushan M Jayarao2, Suresh V Kuchipudi2
1 University of Veterinary and Animal Sciences, Lahore, Pakistan
2 Animal Diagnostic Laboratory, Department of Veterinary and Biomedical Sciences, The Pennsylvania State University, University Park, Pennsylvania, USA
3 Chattogram Veterinary and Animal Sciences University, Chattogram, Bangladesh
4 Bioinformatics Consulting Center, The Pennsylvania State University, University Park, Pennsylvania, USA
5 Department of Veterinary and Animal Sciences, Muhammad Nawaz Shareef University of Agriculture, Multan, Pakistan
| Date of Submission | 15-Jan-2020 |
| Date of Decision | 16-Dec-2020 |
| Date of Acceptance | 18-Dec-2020 |
| Date of Web Publication | 20-Jan-2021 |
Correspondence Address:
Suresh V Kuchipudi
Animal Diagnostic Laboratory, Department of Veterinary and Biomedical Sciences, The Pennsylvania State University, University Park, Pennsylvania
USA
 Source of Support: None, Conflict of Interest: None  | 5 |
DOI: 10.4103/1995-7645.306762
Objective: To carry out the genetic characterization and evolutionary analysis of three avian orthoavulavirus 1 (AOAV-1) isolates from poultry workers with respiratory symptoms.
Methods: Using Illumina MiSeq, whole-genome sequencing was carried out to assess the evolutionary dynamics of three AOAV-1 isolates. A phylogenetic and comparative analysis of all coding genes was done using bioinformatics tools.
Results: Phylogenetic analysis and genetic distance estimation suggested a close relationship among human- and avian-originated velogenic strains of genotype XIII, sub-genotype XIII.2.1. Several substitutions in the significant structural and biological motifs were exclusively identified in the human-originated strains.
Conclusions: To our knowledge, this is the first report of a velogenic AOAV-1 isolate from natural infection of the human upper respiratory tract. Our findings highlight the evolution and zoonotic potential of velogenic AOAV-1 in a disease endemic setting. Keywords: Avian orthoavulavirus 1; Human originated strain; Zoonotic potential; Evolution; Genotype XIII; Poultry workers
How to cite this article:
Shabbir MZ, Nissly RH, Ahad A, Rabbani M, Lim L, Chothe SK, Subbiah M, Sebastian A, Albert I, Ul-Rahman A, Jayarao BM, Kuchipudi SV. Genomic characterization of velogenic avian orthoavulavirus 1 isolates from poultry workers: Implications to emergence and its zoonotic potential towards public health. Asian Pac J Trop Med 2021;14:64-72 |
How to cite this URL:
Shabbir MZ, Nissly RH, Ahad A, Rabbani M, Lim L, Chothe SK, Subbiah M, Sebastian A, Albert I, Ul-Rahman A, Jayarao BM, Kuchipudi SV. Genomic characterization of velogenic avian orthoavulavirus 1 isolates from poultry workers: Implications to emergence and its zoonotic potential towards public health. Asian Pac J Trop Med [serial online] 2021 [cited 2023 Apr 1];14:64-72. Available from: https://www.apjtm.org/text.asp?2021/14/2/64/306762 |
| 1. Introduction | |  |
The International Committee on Taxonomy of Viruses has created a new subfamily Avulavirinae and classified previously designated avian avulaviruses into three genera: Orthoavulavirus, Metaavulavirus, and Paraavulavirus[1]. Avian orthoavulavirus 1 (AOAV-1) belongs to the genus Orthoavulavirus and causes Newcastle disease in multiple susceptible hosts including several avian species[2]. The disease is characterized by the rapid onset of respiratory and nervous symptoms in susceptible species of birds including chickens, pigeons, and turkeys[2].
AOAV-1 is an enveloped, non-segmented negative-sense single-stranded RNA virus with (14.9-17.4) kb in genome length[1]. Based on the clinical disease produced in chicken and typical residue patterns at the fusion protein cleavage site (112RRKQR ↓ F117), the strains of AOAV-1 are classified into three distinct pathotypes (lentogenic, mesogenic, and velogenic strains). The lentogenic strains produce mild infection with low or no mortality; mesogenic strains yield an intermediate virulence, cause respiratory symptoms and may affect egg quality and production; and the velogenic strains are highly virulent and cause severe nervous and respiratory infection with up to 90% mortality in infected flocks[2].
Viruses of the family Paramyxoviridae have the potential to infect a wide range of susceptible hosts including birds, reptiles, mammals, and aquatic species[3]. Besides avian species, AOAV-1 infection has been observed in unusual hosts including humans, mink, swine, and cattle, highlighting the zoonotic potential of this virus[2],[4],[5],[6],[7]. Indeed, the zoonotic potential of AOAV-1 has been recognized for more than 70 years, which was isolated from human subjects, including poultry workers exhibiting conjunctivitis. Accidental exposure of laboratory workers to AOAV-1 by infected liquid splashed into the eyes has been associated with conjunctivitis, and in rare cases, systemic infection with symptoms including headache and general malaise. Occasionally, upper respiratory tract infection and/or general flulike symptoms have been observed in infected humans[8],[9]. Recently, two cases of AOAV-1 infection have been reported from immuno-compromised human individuals who died from deep respiratory tract infection (pneumonia) where the isolated AOAV-1 was presumed to be the cause of infection[4].
In the current report, three closely related AOAV-1 strains that were isolated from nasal swabs of poultry workers, exhibiting mild upper respiratory symptoms, were sequenced for their complete genome and genetically characterized. Unlike previous cases from immuno-compromised patients, the velogenic AOAV-1 isolates described in this study were associated with upper respiratory tract infection. Phylogenomic characterization revealed a close relationship of these human-originated isolates with strains known to exist in chickens from the same geographical region. The human isolates demonstrated several unique amino acid substitutions in the important structural and biological motifs of fusion and hemagglutinin proteins. The finding of AOAV-1 infection in association with upper respiratory tract disease in humans raises the possibility for this virus to adapt to the human host and may gain the ability to transmit and cause infection in other healthy humans. Therefore, there is a need for continuous monitoring and surveillance of AOAV-1 in humans particularly in those linked with poultry settings.
| 2. Materials and methods | | |
2.1. Ethical approval
All essential procedures were approved by the Ethical Review Committee for the Use of Laboratory Animals of the University of Veterinary and Animal Sciences, Lahore, Pakistan (approval number: ERCULA/DR-1349). All applicable international, national, and institutional guidelines for the care and use of animals were followed.
2.2. Sample collection and processing
We previously found extensive seropositivity to AOAV-1 among poultry workers in Pakistan through modified horse red blood cells haemagglutination inhibition assay[10]. Alongside blood collection from poultry workers in these studies, nasal swab samples were collected aseptically from everyone. Swabs were introduced into both nostrils and rotated 3-4 times, then placed in tubes of 2.0 mL brain heart infusion media containing 200 μg/mL gentamicin 5 μg/mL amphotericin B, stored at -80 °C until processing. The tubes were centrifuged at 1 000 ×g for 5 minutes. The resulting supernatant was treated with penicillin (10 000 IU/mL), streptomycin (10 000 μg/mL), amphotericin B (20 μg/mL), gentamicin (1 000 μg/mL), and kanamycin sulfate (600 μg/mL) for 20 minutes at 25 °C and then filtered through a 0.2 μm pore syringe filter. Nine-day old embryonated chicken eggs were inoculated with 0.2 mL of the treated sample via chorioallantoic sac route. Embryos were monitored daily, and all were dead within 60 hours, which is a typical characteristic of velogenic strains of AOAV-1 and pathogenic avian influenza A viruses[11]. The chorioallantoic fluid (CAF) was acquired and processed by spot agglutination test using 10% washed chicken red blood cells in saline as previously described[11]. Extracted RNA from CAF samples showing haemagglutination was tested by real-time reverse-transcript polymerase chain reaction (RT-PCR) for AOAV-1 matrix gene using M+4100 (5′-AGTGATGTGCTCGGACCTTC-3′) and M-4220 (5′-CCTGAGGAGAGGCATTTGCTA-3′) primers and M+4169 probe[12], and using M+25 and M-124 primers and M+64 probe for influenza A virus matrix gene[13].
2.3. Whole-genome sequencing
Viral RNA from the hemagglutination-positive CAF samples was extracted using a MagMAX AI/ND viral RNA extraction kit (Applied Biosystems, Foster City, CA) according to the manufacturer's instructions, and the extracted RNA was stored at -80 °C until processing for next-generation sequencing. A barcoded library was prepared from each RNA sample using the Illumina TruSeq Stranded mRNA kit (San Diego, CA). The manufacturer's protocol was followed using the “purified mRNA input” recommendation, which skips polyA RNA enrichment and begins with RNA fragmentation and priming. An equimolar pool of libraries was made and sequenced on an Illumina MiSeq instrument using 150 bp× 150 bp paired-end sequencing according to the manufacturer's protocol. Between 29947 and 82808 reads were assembled for each genome. Following FastQC, the paired-end reads were merged using FLASH version 1.2.11 and then run against the NCBI nr-protein database using DIAMOND version 0.8.33 and Kraken version 1.1. Results were visualized with MEGAN 6[14], revealing close relation to AOAV-1.
2.4. Viral genome assembly and analysis
Reads belonging to all viruses were extracted, and the sequence was assembled using SPAdes version 3.9.1[15]. From the assembly, the largest contigs (>500 bases) were extracted, and these were used as a query for BLAST against the nr-genome database. The top hit was goose paramyxovirus SF02 (NC_005036) and the second hit was Newcastle disease virus strain B1 (NC_002617). Contigs were aligned against each of these reference strains separately. Nodes were reversely complemented as required and concatenated based on alignment with the goose paramyxovirus reference strain (NC_005036). Where needed, gaps were filled based on alignment with the reference strain.
2.5. Comparative genomic characterization
Assembled sequences were aligned with AOAV-1 strains representing different genotypes accessed from NCBI (http://www.ncbi.nlm. nih.gov/) using ClustalW methods in BioEdit® version 5.0.6[16] for comparative genomic and phylogenetic analysis and prediction of deduced amino acid substitution sites for all coding genes and genetic distance. Consensus sequences from isolates belonging to the same geography isolated within a specific timeframe were used for the identification of residue substitutions. The conserved domains, functionally and structurally important motifs and unique substitutions in open reading frames were predicted using ORF Finder (http://www.ncbi.nlm.nih.gov/gorf/gorf.html), Conserved Domain Prediction (http://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi), and HMMTOP (http://www.enzim.hu/hmmtop/index.php). The potential N-glycosylation sites were predicted using NetNGlyc 1.0 server (http:// www.cbs.dtu.dk/services/NetNGlyc) if the encoded residue was N-X-T/S, where X denoted any residue except proline and accepted if the G-score was 0.5.
2.6. Evolutionary analysis
Phylogenetic analysis and evolutionary distance estimation were performed using MEGA® version 6.0 software[17]. A phylogenetic tree was constructed using aligned complete genome sequences from AOAVs via the neighbor-joining method with 1 000 bootstrap replications. To determine the nucleotide identity and divergence of isolates compared with the AOAVs belonging to genotype XIII, sequence comparison analysis was performed using whole-genome sequences. To determine sub-genotype categories, complete F gene sequences from group XIII isolates were analyzed using the maximum likelihood method with 1 000 bootstrap replications. To estimate the mean inter-population evolutionary diversity (mean evolutionary distance/genetic distance) between sub-genotypes of XIII, the analysis was performed using complete F gene sequences. This analysis was performed via the maximum composite likelihood method (d: Transitions+Transversions model). The rate and pattern of substitutions among sites were modeled with a gamma distribution (parameter=1 with homogenous lineage pattern)[18].
2.7. Recombination analysis
For detection of an occurrence of putative recombination events within the same genotype and/or with vaccine strains, the complete genome sequences of human-originated study strains, avian-originated AOAV-1 isolates from genotype XIII, and vaccine strains [(genotypes Π (LaSota, GenBank accession number AF077761) and III (Mukteshwar, GenBank accession number EF201805)] were used. Recombination analysis was performed for identification of putative breakpoints using SimPlot, GARD (http://www. datamonkey.org/GARD), DAMBE[19], and RDP4 version 4.95[20]. Initially, all enriched sequences were used to check the occurrence of recombination events. Then the outcomes shown by RDP4 were further investigated to enhance its accuracy, clarity, and reliability. This analysis employed seven different recombination algorithm methods (RDP4, GENECONV, BootScan, MaxChi, Chimaera, SiScan and 3Seq) to reveal putative recombinant and parent isolates at P<0.001. Putative recombination events were assumed to have occurred only when they were consistently identified by at least four of the recombination algorithms at a probability threshold of 0.05.
| 3. Results | | |
3.1. Viral agent isolation
CAF from eggs inoculated with material derived from nasal swabs, which were isolated from poultry workers, demonstrated the haemagglutination reaction to horse red blood cells and death of embryos within 60 hours, suggesting infection of embryos with a viral agent. However, the CAF tested negative for both avian influenza A virus and AOAV-1 by real-time RT-PCR. Therefore, viral RNA was extracted from CAF and analyzed by next-generation sequencing to identify the viral agent.
The individuals from whom these nasal swabs originated demonstrated serum antibody titers against the lentogenic AOAV-1 LaSota strain (GenBank accession number: AF077761). Endpoint titers of AOAV-1-specific antibody of these individuals ranged from 160 to 320 (254.0±53.3). These individuals had worked with poultry from 9 months to 8 years and age ranged from 21 to 41 years. One individual reported mild respiratory symptoms, while none of the individuals presented symptoms associated with human AOAV-1 infection, including conjunctivitis.
3.2. AOAV-1 genome assembly
Total reads of 50 to 151 nt length for the three samples were 534941 (PK1, 92.4% of all reads), 696870 (PK2, 91.70%), and 629717 (PK3, 95.70%). Visualization of DIAMOND blast analysis of merged paired-end reads (262277, 228934, and 192498 for PK1, PK2, and PK3, respectively) surprisingly suggested a close match of the sequences with Avulaviruses. Phylogenetic analysis of complete genome sequences showed that these three sequences clustered with AOAV-1 strains of genotype XIII isolated from avian species from Pakistan [Figure 1]. The M gene sequence of the isolates was aligned with the primers and probe [Table 1], which target the M gene that had been used for screening the CAF samples for AOAV-1. The nucleotide identity of the probe with the three human isolates was 84% (21/25) [Table 2]. The complete genome sequences of these AOAV-1 isolates have been deposited in GenBank under the accession numbers MH019281 (PK1), MH019282 (PK2), and MH019283 (PK3). | Figure 1: Phylogenetic relationship of novel human avian orthoavulavirus 1 (AOAV-1) isolates with previously described AOAV-1 strains. Aligned whole-genome sequences of isolates representative of all AOAV-1 genotypes were used to construct phylogenetic trees via the neighbor-joining method with 1 000 bootstrap replications.
Click here to view |
 | Table 1: Sequence of human Pakistani isolates of AOAV-1 compared with sequence of real-time RT-PCR probe[12].
Click here to view |
 | Table 2: Percent divergence of complete genome of human and avian genotype XIII AOAV-1 isolates.
Click here to view |
3.3. Phylogenetic and evolutionary analysis
In the initial phylogenetic analysis of complete genome sequences, distinct clades were observed within genotype XIII with AOAV-1 strains clustering based on location and host species [Figure 1]. The Indian, European, Iranian, and Pakistani avian AOAV-1 strains were clustered into XIII.1.1, XIII.1.2, XIII.2.1 and XIII.2.1 sub-genotypes, as described previously[21],[22],[23],[24]. The human-originated isolates clustered in a distinct clade closely related to avian originated XIII.2.1 isolates from Pakistan [Figure 2]. | Figure 2: Phylogenetic tree was constructed using complete F gene sequence of understudy human originated avian orthoavulavirus 1 (AOAV-1) isolates with previously described AOAV-1 strains. Tree was constructed using the maximum likelihood method with 1 000 bootstrap replications.
Click here to view |
Among human-originated AOAV-1 isolates, 0%-0.01% nucleotide divergence (99.9%-100.0% nucleotide identity) was observed. In contrast, the human-originated strains were 4.39%-4.52% divergent from Pakistani avian AOAV-1 strains during 2010, 5.80%-5.97% divergent from Indian avian AOAV-1 strains isolated during 2013, and 10.09%-10.46% divergent from Indian avian AOAV-1 strains isolated during 2014-2015 [Table 2].
For genetic distance among AOAV strains from genotype XIII, the mean inter-population evolutionary distance was found to be greatest (0.097 1) between sub-genotypes XIII.1.1 and XIII.2.2, with smaller distances between sub-genotypes XIII.1.1 and XIII.2.1 (0.084 3), sub-genotypes XIII.1.2 and XIII.2.2 (0.079 4), sub-genotypes XIII.2.1 and XIII.2.2 (0.068 4), sub-genotypes XIII.1.2 and XIII.2.1 (0.060 5), and sub-genotypes XIII.1.1 and XIII.1.2 (0.054 9) [Figure 2].
3.4. Comparative genomic and residue characterization
All six coding genes of the human-originated AOAVs shared the same nucleotide length with avian-originated AOAVs from genotype XIII. However, a 7 nt insertion (PK2) and 20 nt deletions (PK3) were observed in the non-coding region between NP and P genes whereas, all three human isolates contained a 3 nt insertion in the non-coding region between HN and L genes [Table 3]. Based on predicted amino acid sequences, the cleavage activation site of the fusion protein from the human isolates was found to be 112RRQKR ↓ F117, consistent with velogenic strains of AOAV-1. The combination of these findings and the time until death in embryonated eggs indicate that these isolates are velogenic strains. | Table 3: Comparative genomic characterization of human-originated AOAV-1 with avian-originated classical AOAV-1 isolates.
Click here to view |
Amino acid residue analysis of all coding genes revealed conservation of most of the functional motifs among classical genotype XIII AOAVs strains originating from avian hosts and the human originated AOAVs isolates [Table 4]. A few amino acid substitutions were found in two motifs, the signal peptide motif in the F protein and the hydrophobic anchor region in the HN protein; however, variation was also observed in these motifs within the avian AOAV-1 strains. The number of glycosylation sites in all coding genes (3 in NP, 2 in P, 3 in M, 6 in F, 4 in HN and 8 in L) was conserved except in the HN coding sequence of a Pakistani strain isolated in 2010 and another Indian strain isolated in 2015, where only three glycosylation sites existed due to substitution at position 508 (N508S). Besides, a negligible substitution (R434K) was identified in one glycosylation site of the HN protein (433NRT435) in all avian-originated AOAVs as compared to human-originated AOAVs. A total of 28 single amino acid substitutions were exclusive to human-originated AOAV-1 strains as compared to avian-originated classical AOAV-1 strains isolated from Pakistan and India. | Table 4: Amino acid characterization of avian and human AOAV-1 isolates from Pakistan and India.
Click here to view |
3.5. Detection of potential recombination events
SimPlot showed similarities in non-coding intergenic regions among selected strains. No potential recombination events for the human isolates were found by detecting putative recombination events or breakpoints integrated shown in RDP4 software (data not shown).
| 4. Discussion | | |
The present report on AOAV-1 infection in humans provides evidence of the zoonotic potential as well as a continued evolution among genotype XIII viruses. The study not only identifies a deficiency in the M-gene-based screening assay used for the initial identification of AOAV-1 but also signifies the importance of advanced sequencing tools (e.g. next-generation sequencing) and associated bioinformatics pipelines for the identification of other unknown or unleashed viruses of public health significance.
Based on evolutionary distance and phylogenetic analyses, the three human AOAV-1 isolates clustered together with avian originated strain from Pakistan. The study isolates had an average nucleotide distance of 7.9% (XIII.1.1), 6.0% (XIII.1.2), 2.8% (XIII.2.1), and 6.1% (XIII.2.2) and therefore was classified within sub-genotype XIII.2.1 as per recently proposed classification criteria[24]. Historically, the most ancestral strain of genotype XIII AOAV-1 was firstly isolated from a cockatoo in India in 1982. Over a period of time, the virus is believed to have an ongoing evolution and this is fairly evidenced by the identification and classification of variants or sub-genotype within genotype XIII from countries including Iran[21], Pakistan[25], Bangladesh[22], and India[23]. It is important to indicate that the previous viruses of genotype XIII have exclusively been reported from multiple avian hosts; nevertheless, to the best of our knowledge, this is the first evidence of isolation, identification, and genomic characterization of genotype XIII (XIII.2.1) viruses among human in close contact with poultry.
The isolates encoded typical cleavage motif of the F protein and exhibited mean time-to-death in embryonated chicken eggs, consistent with the definition of velogenic virulent AOAV-1 strains. Theoretically, virulent viruses need only a short period of time for replication in a new environment and therefore demonstrate accelerated rates of adaptation[26],[27]. This aspect could be correlated with the zoonotic potential of these velogenic strains for causing infection in humans. In addition, some substitutions were exclusively observed in human-originated strains, highlighting the capability of AOAV-1 as a mutable RNA virus for host adaptation and survival of fitness in the new environment. In fact, a continuous accumulation of several point mutations leading to amino acid substitutions combined with host immune pressure drives the evolution of AOAV-1 and may contribute to the emergence of novel virulent strains with a potentially altered host tropism[22],[23],[24]. Notably, several nucleotide substitutions in the functional motifs of the F protein signal peptide and HN protein hydrophobic signal anchor were observed in the human isolates, similar to a fatal human AOAV-1 isolate[4]. This is important because F protein mediates fusion between the viral lipid membrane and the host cell membrane, and the HN protein promotes fusion between the two membranes through interaction with the F protein. Therefore, further investigation of the effect of these substitutions on viral attachment and/or replication of these strains in human and avian cell types is warranted.
Previous studies have demonstrated the ability of AOAV-1 to replicate in hamsters, rats, guinea pigs, cattle, swine, sheep, and non-human primates[28],[29]. When passaged intracerebrally in hamsters, a virulent strain of AOAV-1 was able to infect hamsters by multiple routes and cause disease in mice, sheep, cattle, and rhesus macaques[30],[31],[32]. Together, these findings in the past decades support the possibility for AOAV-1 to adapt to a mammalian host. Indeed, the use of AOAV-1 as a vaccine vector and oncolytic therapy in humans and other mammals shows the reality of this potential[33]. Until the last decade, only one natural infection of a non-avian livestock animal with AOAV-1 had been described in a bovine calf[30]. Since 2009, increasing numbers of natural infection of mammals with AOAV-1 have been described, including healthy sheep and swine demonstrating deadly respiratory symptoms, and mink with encephalitis[5],[31],[34]. Notably, the more recent isolations from mammals have been associated with the disease, and viruses isolated from mink and cattle were velogenic strains.
Human infection with AOAV-1 was first described in 1943 by Burnet and others following an accidental ocular exposure of a laboratory worker to a vaccine strain of AOAV-1. Since then, additional laboratory exposures, as well as occupational exposure in commercial poultry workers have produced clinical symptoms mainly of conjunctivitis. Two cases of deep respiratory infection with AOAV-1 have been implicated as the cause of death in immuno-compromised humans[4], highlighting the zoonotic potential of AOAV-1 in susceptible individuals. Evidence of AOAV-1 causing upper respiratory symptoms in humans is little[35],[36],[37] and has not been described in detail till now. In the present study, we sequenced and genetically characterized three virulent AOAV-1 isolates from one individual with mild upper respiratory symptoms and from two individuals demonstrating no symptoms and reporting no recent respiratory symptom. To our knowledge, this is the first report of a velogenic AOAV-1 isolate from natural infection of the human upper respiratory tract.
Evidence for the ability of AOAV-1 to contribute to serious illness and death in humans is rare, but it has repeatedly been suggested that this virus should be included as a suspect in the inexplicable respiratory disease of non-avian species[4],[5],[31]. The new evidence presented in this study further highlights the need for further investigation of AOAV-1 evolution in mammals to investigate its zoonotic potential. Monitoring the incidence of AOAV-1 and its genetic diversity in humans may prove to be an important tool in predicting or preventing the adaptation of AOAV-1 to not only human lungs but also in other mammals as suggested in previous studies[31]. The ability to adapt to these diverse hosts raises concerns that AOAV-1 replicating in mammals, including humans, may acquire the ability to spread between mammals of the same species and/or cause more severe morbidity than historically observed. Indeed, the AOAV-1 isolated from a patient with fatal respiratory failure led to severe pneumonia when experimentally inoculated to cynomolgus macaques[32]. Continued surveillance of AOAV-1 in mammals, especially humans who are likely to contract the virus, is advisable to be prepared for possible epizootics. In addition, studies aimed at understanding the pathogenesis of AOAV-1 in mammalian hosts, including non-human primates[32], should be further pursued to better understand the potential for AOAVs to cause mild or severe infection in healthy individuals and immuno-compromised patients, respectively.
We characterized three human-originated velogenic AOAV-1 strains and found them closely related to those strains isolated from avian species within the same geography. The highest genetic homology of viruses from different origins (human and avian species) raised the concern towards the adaptation of novel host with implications of its zoonotic potential. It may suggest that the evolutionary dynamics of the virus occur naturally to adapt to the environment or a new host or perhaps confer virulent viruses with additional evolutionary advantages to allow persistence in nature.
Conflict of interest statement
The authors declare that they have no conflict of interest.
Acknowledgments
This work was supported by the startup research grant of the Pennsylvania State University (S.V.K.). The authors wish to thank Craig Praul and the Penn State Genomics Core Facility, University Park, PA for useful advice and guidance for planning the whole genome sequencing.
Authors’ contributions
SVK and MZS designed the study. MZS, MR, AR, AA performed virus isolation and genome sequencing. RHN, LL, SKC, AS, IA, BMG analyzed the data. SVK, RHN and MZS wrote the initial draft manuscript all authors revised the manuscript.
| References | | |
| 1. | Amarasinghe GK, Ayllón MA, Bào Y, Basler CF, Bavari S, Blasdell KR, et al. Taxonomy of the order Mononegavirales: Update 2019. Arch Virol 2019; 164(7): 1967-1980.  |
| 2. | Kaleta EF, Baldauf C. Newcastle disease in free-living and pet birds. In: Newcastle disease. Boston: Springer; 1988, p. 197-246. |
| 3. | Samal SK. The biology of paramyxoviruses. Norfolk: Caister Academic Press; 2011, p. 469. |
| 4. | Goebel SJ, Taylor J, Barr BC, Kiehn TE, Castro-Malaspina HR, Hedvat CV, et al. Isolation of avian paramyxovirus 1 from a patient with a lethal case of pneumonia. J Virol 2017; 81(22): 12709-12714. |
| 5. | Zhao P, Sun L, Sun X, Li S, Zhang W, Pulscher LA, et al. Newcastle disease virus from domestic mink, China, 2014. Vet Microbiol 2017; 198: 104-107. |
| 6. | Shabbir MZ, Nissly RH, Ahad A, Rabbani M, Chothe SK, Sebastian A, et al. Complete genome sequences of three related avian avulavirus 1 isolates from poultry farmers in Pakistan. Genome Announc 2018; 6(18): e00361-18. doi: 10.1128/genomeA.00361-18. |
| 7. | Ul-Rahman A, Shabbir MZ. A comparative phylogenomic analysis of avian avulavirus 1 isolated from non-avian hosts: Conquering new frontiers of zoonotic potential among species. Arch Virol 2019; 164(7): 1771-1780. |
| 8. | Capua I, Alexander DJ. Human health implications of avian influenza viruses and paramyxoviruses. Europ J Clin Microb Infect Dis 2004; 23(1): 1-6. |
| 9. | Chang PW. Newcastle disease. In: Beran GW, ed. CRC handbook series in zoonoses. Boca Raton: CRC Press; 1981, p. 261-274. |
| 10. | Ahad A, Rabbani M, Yaqub T, Mahmood A, Kuthu ZH, Shabbir MZ, et al. Detection of antibody to Newcastle disease virus in human sera in Pakistan. J Anim Plant Sci 2012; 23(4): 990-994. |
| 11. | Office International des Epizooties World Organization for Animal Health. Manual of diagnostic tests and vaccines for terrestrial animals, Paris, 2012. [Online]. Available from: https://www.oie.int/doc/ged/D12009.PDF. [Accessed on 20 January 2020]. |
| 12. | Wise MG, Suarez DL, Seal BS, Pedersen JC, Senne DA, King DJ, et al. Development of a real-time reverse-transcription PCR for detection of Newcastle disease virus RNA in clinical samples. J Clin Microbiol 2004; 42(1): 329-338. |
| 13. | Spackman E, Senne DA, Myers TJ, Bulaga LL, Garber LP, Perdue ML, et al. 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(9): 3256-3260. |
| 14. | Huson DH, Beier S, Flade I, Górska A, El-Hadidi M, Mitra S, et al. MEGAN community edition-interactive exploration and analysis of large-scale microbiome sequencing data. PLoS Comp Biol 2016; 12(6):e1004957. |
| 15. | Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J Comp Biol 2012; 19(5): 455-477. |
| 16. | Hall TA. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp 1999; 41: 95-98. |
| 17. | Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol Biol Evol 2013; 30(12): 2725-2729. |
| 18. | Tamura K, Kumar S. Evolutionary distance estimation under heterogeneous substitution pattern among lineages. Mol Biol Evol 2002; 19(10): 1727-1736. |
| 19. | Xia X, Xie Z. DAMBE: Software package for data analysis in molecular biology and evolution. J Heredity 2001; 92(4): 371-373. |
| 20. | Martin DP, Murrell B, Golden M, Khoosal A, Muhire B. RDP4: Detection and analysis of recombination patterns in virus genomes. Virus Evol 2015; 1(1): vev003. doi: 10.1093/ve/vev003. |
| 21. | Ebrahimi MM, Shahsavandi S, Moazenijula G, Shamsara M. Phylogeny and evolution of Newcastle disease virus genotypes isolated in Asia during 2008-2011. Virus Genes 2012; 45(1): 63-68. |
| 22. | Barman LR, Nooruzzaman M, Sarker RD, Rahman MT, Saife MR, Giasuddin M, et al. Phylogenetic analysis of Newcastle disease viruses from Bangladesh suggests continuing evolution of genotype XIII. Arch Virol 2017; 162(10): 3177-31782. |
| 23. | Nath B, Kumar S. Emerging variant of genotype XIII Newcastle disease virus from Northeast India. Acta Trop 2017; 172: 64-69. |
| 24. | Dimitrov KM, Abolnik C, Afonso CL, Albina E, Bahl J, Berg M, et al. Updated unified phylogenetic classification system and revised nomenclature for Newcastle disease virus. Infect Genet Evol 2019; 103917. doi:10.1016/j.meegid.2019.103917. |
| 25. | 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. |
| 26. | Perozo F, Merino R, Afonso CL, Villegas P, Calderon N. Biological and phylogenetic characterization of virulent Newcastle disease virus circulating in Mexico. Avian Dis 2008; 52(3): 472-479. |
| 27. | Qin ZM, Tan LT, Xu HY, Ma BC, Wang YL, Yuan XY, et al. Pathotypical characterization and molecular epidemiology of Newcastle disease virus isolates from different hosts in China from 1996 to 2005. J Clin Microbiol 2008; 46(2): 601-611. |
| 28. | Hofstad MS. Experimental inoculation of swine and sheep with Newcastle disease virus. Cornell Vet 1950; 40: 190-197. |
| 29. | Samuel AS, Subbiah M, Shive H, Collins PL, Samal SK. Experimental infection of hamsters with avian paramyxovirus serotypes 1 to 9. Vet Res 2011; 42(1): 38. |
| 30. | Yates VJ, Fry DE, Henderson Jr. BW. Isolation of Newcastle disease virus from a calf. J Am Vet Med Assoc 1952; 120: 149-150. |
| 31. | Sharma B, Pokhriyal M, Rai GK, Saxena M, Ratta B, Chaurasia M, et al. Isolation of Newcastle disease virus from a non-avian host (sheep) and its implications. Arch Virol 2012; 157(8): 1565-1567. |
| 32. | Kuiken T, Buijs P, Van Run P, Van Amerongen G, Koopmans M, Van den Hoogen B. Pigeon paramyxovirus type 1 from a fatal human case induces pneumonia in experimentally infected cynomolgus macaques ( Macaca fascicularis). Vet Res 2017; 48(1): 80. |
| 33. | Zamarin D, Palese P. Oncolytic Newcastle disease virus for cancer therapy: Old challenges and new directions. Future Microbiol 2012; 7(3): 347-367. |
| 34. | Chen S, Hao H, Wang X, Du E, Liu H, Yang T, et al. Genomic characterisation of a lentogenic Newcastle disease virus strain HX01 isolated from sick pigs in China. Virus Genes 2013; 46(2): 264-270. |
| 35. | Mitchell CA, Walker RV. Note on the infection of a person with Newcastle disease virus. Canadian J Comp Med Vet Sci 1951; 15(9): 226. |
| 36. | Quinn RW, Hanson RP, Brown JW, Brandly CA. Newcastle disease virus in man: Results of studies in five cases. J Lab Clin Med 1952; 40(5): 736-743. |
| 37. | Dardiri AH, Yates VJ, Flanagan TD. The reaction to infection with the B1 strain of Newcastle disease virus in man. American J Vet Res 1962; 23: 918. |
[Figure 1], [Figure 2]
[Table 1], [Table 2], [Table 3], [Table 4]
|
Search Pubmed for