A robust and cost-effective approach to sequence and analyze complete genomes of small RNA viruses

The first report and biological characterization of Avian Orthoavulavirus 16 in wild migratory waterfowl and domestic poultry in China reveal a potential threat to birds

AbstractThe Avulavirus within the family Paramyxoviridae includes at least 22 different species, and is known to cause different types of infections and even be fatal in multiple avian species. There is limited knowledge of the genetic and biological information of Avulavirus species -2 to 22 in domestic and wild birds and the disease significance of these viruses in birds is not fully determined, although as many as 10 new distinct species have been identified from wild birds and domestic poultry around the world in the last decade. This study aimed to use PCR, virus isolation, and sequencing to genetically and biologically characterize Avian Orthoavulavirus 16 (AOAV-16) in wild birds and domestic poultry collected from different locations in China between 2014 and 2022. Of five isolated AOAV-16 strains (Y1 to Y5), only the Y4 strain had a hemagglutination (HA)-negative result. All of these isolates were low virulent viruses for chickens, except Y3 which was detected simultaneously with avian influenza virus (AIV) of H9N2 subtype. Furthermore, at least four different types of intergenic sequences (IGS) between the HN and L genes junction, and the recombination event as well as interspecific transmission by wild migratory birds, existed within the species AOAV-16. These findings and results of other reported AOAV-16 strains recommend strict control measures to limit contact between wild migratory birds and domestic poultry and imply potential threats to commercial poultry and even public health challenges worldwide.

The 28S rRNA RT-qPCR assay for host depletion evaluation to enhance avian virus detection in Illumina and Nanopore sequencing

Abundant host and bacterial sequences can obscure the detection of less prevalent viruses in untargeted next-generation sequencing (NGS). Efficient removal of these non-targeted sequences is vital for accurate viral detection. This study presents a novel 28S ribosomal RNA (rRNA) RT-qPCR assay designed to assess the efficiency of avian rRNA depletion before conducting costly NGS for the detection of avian RNA viruses. The comprehensive evaluation of this 28S-test focuses on substituting DNase I with alternative DNases in our established depletion protocols and finetuning essential parameters for reliable host rRNA depletion. To validate the effectiveness of the 28S-test, we compared its performance with NGS results obtained from both Illumina and Nanopore sequencing platforms. This evaluation utilized swab samples from chickens infected with highly pathogenic avian influenza virus, subjected to established and modified depletion protocols. Both methods significantly reduced host rRNA levels, but using the alternative DNase had superior performance. Additionally, utilizing the 28S-test, we explored cost- and time-effective strategies, such as reduced probe concentrations and other alternative DNase usage, assessed the impact of filtration pre-treatment, and evaluated various experimental parameters to further optimize the depletion protocol. Our findings underscore the value of the 28S-test in optimizing depletion methods for advancing improvements in avian disease research through NGS.

Next-Generation Sequencing for the Detection of Microbial Agents in Avian Clinical Samples

Direct-targeted next-generation sequencing (tNGS), with its undoubtedly superior diagnostic capacity over real-time PCR (RT-PCR), and direct-non-targeted NGS (ntNGS), with its higher capacity to identify and characterize multiple agents, are both likely to become diagnostic methods of choice in the future. tNGS is a rapid and sensitive method for precise characterization of suspected agents. ntNGS, also known as agnostic diagnosis, does not require a hypothesis and has been used to identify unsuspected infections in clinical samples. Implemented in the form of multiplexed total DNA metagenomics or as total RNA sequencing, the approach produces comprehensive and actionable reports that allow semi-quantitative identification of most of the agents present in respiratory, cloacal, and tissue samples. The diagnostic benefits of the use of direct tNGS and ntNGS are high specificity, compatibility with different types of clinical samples (fresh, frozen, FTA cards, and paraffin-embedded), production of nearly complete infection profiles (viruses, bacteria, fungus, and parasites), production of "semi-quantitative" information, direct agent genotyping, and infectious agent mutational information. The achievements of NGS in terms of diagnosing poultry problems are described here, along with future applications. Multiplexing, development of standard operating procedures, robotics, sequencing kits, automated bioinformatics, cloud computing, and artificial intelligence (AI) are disciplines converging toward the use of this technology for active surveillance in poultry farms. Other advances in human and veterinary NGS sequencing are likely to be adaptable to avian species in the future.

Complete genome sequence of seven virulent Newcastle disease virus isolates of sub-genotype XIII.1.1 from Tanzania

We report the complete genome sequences of seven virulent Newcastle disease viruses (NDVs) that were isolated from chickens from live bird markets in the Arusha, Iringa, Mbeya, and Tanga regions of Tanzania in 2012. Phylogenetic analysis revealed that all isolates belong to sub-genotype XIII.1.1.

RNA sequencing of avian paramyxovirus (Paramyxoviridae, Avulavirinae) isolates from wild mallards in Belgium, 2021: complete genomes and coinfections

We used untargeted RNA sequencing to characterize three Avulavirinae isolates from pooled samples obtained from wild mallards in Belgium in 2021. The complete genome sequences of two avian Orthoavulavirus-1 (AOAV-1) strains and one avian Paraavulavirus-4 (APMV-4) strain were determined confirming hemagglutination inhibition testing of the virus isolates. In addition, the applied sequencing strategy identified an avian influenza virus (AIV) coinfection in all three virus isolates, confirming weak-positive AIV realtime RT-PCR results from the original sample material. In one AOAV-1 isolate, partial sequences covering all genome segments of an AIV of subtype H11N9 could be de novo assembled from the sequencing data. Besides an AIV coinfection, RNA metagenomic data from the APMV-4 isolate also showed evidence of Alpharetrovirus and Megrivirus coinfection. In total, two AOAV-1 of Class II, genotype I.2 and one APMV-4 complete genome sequences were assembled and compared to publicly available sequences, highlighting the importance of surveillance for poultry pathogens in wild birds. Beyond the insights from full genome characterization of virus isolates, untargeted RNA sequencing strategies provide additional insights in the RNA virome of clinical samples as well as their derived virus isolates that are particularly useful when targeting wild avifauna reservoirs of poultry pathogens.

Phylogenetic and antigenic analysis of infectious bronchitis virus isolated from commercial and backyard chickens in Pakistan, 2015-2018

Infectious bronchitis virus (IBV) is a rapidly evolving virus affecting both vaccinated and unvaccinated poultry flocks and is responsible for significant economic losses globally; hence, it is imperative to obtain a deeper understanding of this pathogen. In this study, seven IBV strains were isolated from commercial and backyard poultry flocks during 2015-2018. We obtained full-length IBV genomes of two viruses using the Illumina sequencing method, while five additional viruses were genetically characterized through full-length spike (S1) gene sequencing. Phylogenetic and distance analysis based on complete S1 gene and full-length genome sequences revealed that one IBV isolate belonged to genotype GI-1 and six viruses were clustered within genotype GI-13. Deduced amino acid sequences of GI-13 strains exhibited 31.8-37.2 % divergence with the commonly used classic vaccine strains (M41) and 2.7-12.6 % with variant vaccine strains (4/91) in Pakistan. High evolutionary distances suggest that the IBV viruses circulating in Pakistan are under continuous evolutionary pressure. Moreover, ch/IBV/Pak/AW-2/2017 was found to have originated from an intra-genotypic recombination event between the variant group (GI-23 lineage as a major parent) and variant vaccine strain (4/91-like as a minor parent) and is the first example of recombination within genotype GI-13 in Pakistan. Together, these findings provide genetic and evolutionary insights into the currently circulating IBV genotypes in Pakistan, which could help to better understand the origin, spread and evolution of IBVs, and to ascertain the importance of disease monitoring as well as re-evaluation forof currently used vaccines and vaccination programmes.

Altered human gut virome in patients undergoing antibiotics therapy for Helicobacter pylori

Transient gut microbiota alterations have been reported after antibiotic therapy for Helicobacter pylori. However, alteration in the gut virome after H. pylori eradication remains uncertain. Here, we apply metagenomic sequencing to fecal samples of 44 H. pylori-infected patients at baseline, 6-week (N = 44), and 6-month (N = 33) after treatment. Following H. pylori eradication, we discover contraction of the gut virome diversity, separation of virome community with increased community difference, and shifting towards a higher proportion of core virus. While the gut microbiota is altered at 6-week and restored at 6-month, the virome community shows contraction till 6-month after the treatment with enhanced phage-bacteria interactions at 6-week. Multiple courses of antibiotic treatments further lead to lower virus community diversity when compared with treatment naive patients. Our results demonstrate that H. pylori eradication therapies not only result in transient alteration in gut microbiota but also significantly alter the previously less known gut virome community.

Genome Sequencing and Characterization of an Avian Orthoavulavirus 1 VG/GA-like Isolate with a Unique Fusion Cleavage Site Motif

A complete genome sequence of a VG/GA -like strain of avian orthoavulavirus 1 (AOAV-1) was identified by nontargeted next-generation sequencing of an oropharyngeal swab sample collected from a carcass of a 12-mo-old backyard chicken. The isolate has a fusion (F) protein cleavage site motif consistent with a low virulent AOAV-1, but it has a unique motif with phenylalanine at position 117 (¹¹²G-R-Q-G-R↓F¹¹⁷), which is typical for virulent AOAV-1 strains. The one nucleotide difference at the cleavage site compared to other low-virulence viruses made the isolate detectable by F-gene-specific real-time reverse transcription-PCR (rRT-PCR) developed as a diagnostic test to specifically detect virulent strains. The mean death time determined in eggs and intracerebral pathogenicity index determined in chickens classified the isolate as lentogenic. This is the first report of a lentogenic VG/GA-like virus with a phenylalanine residue at position 117 of the F protein cleavage site in the United States. In addition to concern for potential pathogenic shift of the virus through additional changes at the cleavage site, our finding warrants increased awareness of diagnosticians of potential false positive F-gene rRT-PCR tests.

From Clinical Specimen to Whole Genome Sequencing of A(H3N2) Influenza Viruses: A Fast and Reliable High-Throughput Protocol

(1) Background: Over the last few years, there has been growing interest in the whole genome sequencing (WGS) of rapidly mutating pathogens, such as influenza viruses (IVs), which has led us to carry out in-depth studies on viral evolution in both research and diagnostic settings. We aimed at describing and determining the validity of a WGS protocol that can obtain the complete genome sequence of A(H3N2) IVs directly from clinical specimens. (2) Methods: RNA was extracted from 80 A(H3N2)-positive respiratory specimens. A one-step RT-PCR assay, based on the use of a single set of specific primers, was used to retro-transcribe and amplify the entire IV type A genome in a single reaction, thus avoiding additional enrichment approaches and host genome removal treatments. Purified DNA was quantified; genomic libraries were prepared and sequenced by using Illumina MiSeq platform. The obtained reads were evaluated for sequence quality and read-pair length. (3) Results: All of the study specimens were successfully amplified, and the purified DNA concentration proved to be suitable for NGS (at least 0.2 ng/µL). An acceptable coverage depth for all eight genes of influenza A(H3N2) virus was obtained for 90% (72/80) of the clinical samples with viral loads >105 genome copies/mL. The mean depth of sequencing ranged from 105 to 200 reads per position, with the majority of the mean depth values being above 103 reads per position. The total turnaround time per set of 20 samples was four working days, including sequence analysis. (4) Conclusions: This fast and reliable high-throughput sequencing protocol should be used for influenza surveillance and outbreak investigation.

Genome Sequence Variations of Infectious Bronchitis Virus Serotypes From Commercial Chickens in Mexico

New variants of infectious bronchitis viruses (IBVs; Coronaviridae) continuously emerge despite routine vaccinations. Here, we report genome sequence variations of IBVs identified by random non-targeted next generation sequencing (NGS) of vaccine and field samples collected on FTA cards from commercial flocks in Mexico in 2019–2021. Paired-ended sequencing libraries prepared from rRNA-depleted RNAs were sequenced using Illumina MiSeq. IBV RNA was detected in 60.07% (n = 167) of the analyzed samples, from which 33 complete genome sequences were de novo assembled. The genomes are organized as 5'UTR-[Rep1a-Rep1b-S-3a-3b-E-M-4b-4c-5a-5b-N-6b]-3'UTR, except in eight sequences lacking non-structural protein genes (accessory genes) 4b, 4c, and 6b. Seventeen sequences have auxiliary S2' cleavage site located 153 residues downstream the canonically conserved primary furin-specific S1/S2 cleavage site. The sequences distinctly cluster into lineages GI-1 (Mass-type; n = 8), GI-3 (Holte/Iowa-97; n = 2), GI-9 (Arkansas-like; n = 8), GI-13 (793B; n = 14), and GI-17 (California variant; CAV; n = 1), with regional distribution in Mexico; this is the first report of the presence of 793B- and CAV-like strains in the country. Various point mutations, substitutions, insertions and deletions are present in the S1 hypervariable regions (HVRs I-III) across all 5 lineages, including in residues 38, 43, 56, 63, 66, and 69 that are critical in viral attachment to respiratory tract tissues. Nine intra-/inter-lineage recombination events are present in the S proteins of three Mass-type sequences, two each of Holte/Iowa-97 and Ark-like sequence, and one each of 793B-like and CAV-like sequences. This study demonstrates the feasibility of FTA cards as an attractive, adoptable low-cost sampling option for untargeted discovery of avian viral agents in field-collected clinical samples. Collectively, our data points to co-circulation of multiple distinct IBVs in Mexican commercial flocks, underscoring the need for active surveillance and a review of IBV vaccines currently used in Mexico and the larger Latin America region.

Non-target RNA depletion strategy to improve sensitivity of next-generation sequencing for the detection of RNA viruses in poultry

PCR-based assays have become the benchmark for detecting pathogens of poultry and other livestock; however, these techniques are limited in their ability to detect multiple infecting agents, provide limited genetic information on the pathogen, and, for RNA viruses, must be reviewed frequently to assure high sensitivity and specificity. In contrast, untargeted, high-throughput sequencing can rapidly detect all infecting agents in a sample while providing genomic sequence information to allow more in-depth characterization of viruses. Although next-generation sequencing (NGS) offers many advantages, one of its primary limitations is low sensitivity to pathogens given the abundance of host and other non-target sequences in sequencing libraries. We explored methods for improving the sensitivity of NGS to detect respiratory and enteric viruses in poultry from RNA extracts of swab samples. We employed commercial and custom-designed negative enrichment strategies to selectively deplete the most abundant rRNA reads from the host and non-target bacteria; host RNA was diminished from up to 40% of total reads to as low as 3%, and the total number of reads assigned to abundant bacterial classes were reduced greatly. Our treatment resulted in up to a 700-fold increase in the number of viral reads, detection of a greater number of viral agents, and higher average genome coverage for pathogens. Depletion assays added only 2 h to the NGS library preparation workflow. Custom depletion probe design offered significant cost savings (US$7-12 per sample) compared to commercial kits (US$30-50 per sample).

Low Pathogenicity H7N3 Avian Influenza Viruses Have Higher Within-Host Genetic Diversity Than a Closely Related High Pathogenicity H7N3 Virus in Infected Turkeys and Chickens

Within-host viral diversity offers a view into the early stages of viral evolution occurring after a virus infects a host. In recent years, advances in deep sequencing have allowed for routine identification of low-frequency variants, which are important sources of viral genetic diversity and can potentially emerge as a major virus population under certain conditions. We examined within-host viral diversity in turkeys and chickens experimentally infected with closely related H7N3 avian influenza viruses (AIVs), specifically one high pathogenicity AIV (HPAIV) and two low pathogenicity AIV (LPAIVs) with different neuraminidase protein stalk lengths. Consistent with the high mutation rates of AIVs, an abundance of intra-host single nucleotide variants (iSNVs) at low frequencies of 2–10% was observed in all samples collected. Furthermore, a small number of common iSNVs were observed between turkeys and chickens, and between directly inoculated and contact-exposed birds. Notably, the LPAIVs have significantly higher iSNV diversities and frequencies of nonsynonymous changes than the HPAIV in both turkeys and chickens. These findings highlight the dynamics of AIV populations within hosts and the potential impact of genetic changes, including mutations in the hemagglutinin gene that confers the high pathogenicity pathotype, on AIV virus populations and evolution.

Novel avian orthoavulavirus 13 in wild migratory waterfowl: biological and genetic considerations

Avian orthoavulavirus 13 (AOAV-13), formerly known as Avian paramyxovirus 13 (APMV-13), is found scatteredly in wild birds around the world. Although four complete genome sequences of AOAV-13 had been identified since the first discovery in Japan in 2003, the information available on the genetic variation and biological characteristics of AOAV-13 is still limited. In the present study, we isolated six AOAV-13 strains from fecal samples of wild migratory waterfowls during annual (2014-2018) viral surveillance of wild bird populations from wetland and domestic poultry of live bird markets (LBMs) in China. The phylogenetic analyses based on the HN and F genes showed that they had very close relationship and the molecular clock estimations showed a low evolutionary rate of AOAV-13. However, Bean goose/Hubei/V97-1/2015 is 1953 nt in size (ORF, 1, 776 nt), which is a unique size and longer than other reported AOAV-13 strains. Additionally, four repeats of conserved sequences "AAAAAT" was presented in the 5'-end trailer region of Swan goose/Hubei/VI49-1/2016, which is unprecedented in the AOAV-13. These findings highlight the importance of continuous monitoring the specific species of APMVs.

The Pathobiology of H7N3 Low and High Pathogenicity Avian Influenza Viruses from the United States Outbreak in 2020 Differs between Turkeys and Chickens

An outbreak caused by H7N3 low pathogenicity avian influenza virus (LPAIV) occurred in commercial turkey farms in the states of North Carolina (NC) and South Carolina (SC), United States in March of 2020. Subsequently, H7N3 high pathogenicity avian influenza virus (HPAIV) was detected on a turkey farm in SC. The infectivity, transmissibility, and pathogenicity of the H7N3 HPAIV and two LPAIV isolates, including one with a deletion in the neuraminidase (NA) protein stalk, were studied in turkeys and chickens. High infectivity [<2 log₁₀ 50% bird infectious dose (BID₅₀)] and transmission to birds exposed by direct contact were observed with the HPAIV in turkeys. In contrast, the HPAIV dose to infect chickens was higher than for turkeys (3.7 log₁₀ BID₅₀), and no transmission was observed. Similarly, higher infectivity (<2–2.5 log₁₀ BID₅₀) and transmissibility were observed with the H7N3 LPAIVs in turkeys compared to chickens, which required higher virus doses to become infected (5.4–5.7 log₁₀ BID₅₀). The LPAIV with the NA stalk deletion was more infectious in turkeys but did not have enhanced infectivity in chickens. These results show clear differences in the pathobiology of AIVs in turkeys and chickens and corroborate the high susceptibility of turkeys to both LPAIV and HPAIV infections.

Avian Paramyxovirus Type 1 in Egypt: Epidemiology, Evolutionary Perspective, and Vaccine Approach

Avian orthoavulavirus 1, formerly known as avian paramyxovirus type-1 (APMV-1), infects more than 250 different species of birds. It causes a broad range of clinical diseases and results in devastating economic impact due to high morbidity and mortality in addition to trade restrictions. The ease of spread has allowed the virus to disseminate worldwide with subjective virulence, which depends on the virus strain and host species. The emergence of new virulent genotypes among global epizootics, including those from Egypt, illustrates the time-to-time genomic alterations that lead to simultaneous evolution of distinct APMV-1 genotypes at different geographic locations across the world. In Egypt, the Newcastle disease was firstly reported in 1947 and continued to occur, despite rigorous prophylactic vaccination, and remained a potential threat to commercial and backyard poultry production. Since 2005, many researchers have investigated the nature of APMV-1 in different outbreaks, as they found several APMV-1 genotypes circulating among various species. The unique intermingling of migratory, free-living, and domesticated birds besides the availability of frequently mobile wild birds in Egypt may facilitate the evolution power of APMV-1 in Egypt. Pigeons and waterfowls are of interest due to their inclusion in Egyptian poultry industry and their ability to spread the infection to other birds either by presence of different genotypes (as in pigeons) or by harboring a clinically silent disease (as in waterfowl). This review details (i) the genetic and pathobiologic features of APMV-1 infections in Egypt, (ii) the epidemiologic and evolutionary events in different avian species, and (iii) the vaccine applications and challenges in Egypt.

Virus Adaptation Following Experimental Infection of Chickens with a Domestic Duck Low Pathogenic Avian Influenza Isolate from the 2017 USA H7N9 Outbreak Identifies Polymorphic Mutations in Multiple Gene Segments

In March 2017, highly pathogenic (HP) and low pathogenic (LP) avian influenza virus (AIV) subtype H7N9 were detected from poultry farms and backyard birds in several states in the southeast United States. Because interspecies transmission is a known mechanism for evolution of AIVs, we sought to characterize infection and transmission of a domestic duck-origin H7N9 LPAIV in chickens and genetically compare the viruses replicating in the chickens to the original H7N9 clinical field samples used as inoculum. The results of the experimental infection demonstrated virus replication and transmission in chickens, with overt clinical signs of disease and shedding through both oral and cloacal routes. Unexpectedly, higher levels of virus shedding were observed in some cloacal swabs. Next generation sequencing (NGS) analysis identified numerous non-synonymous mutations at the consensus level in the polymerase genes (i.e., PA, PB1, and PB2) and the hemagglutinin (HA) receptor binding site in viruses recovered from chickens, indicating possible virus adaptation in the new host. For comparison, NGS analysis of clinical samples obtained from duck specimen collected during the outbreak indicated three polymorphic sides in the M1 segment and a minor population of viruses carrying the D139N (21.4%) substitution in the NS1 segment. Interestingly, at consensus level, A/duck/Alabama (H7N9) had isoleucine at position 105 in NP protein, similar to HPAIV (H7N9) but not to LPAIV (H7N9) isolated from the same 2017 influenza outbreak in the US. Taken together, this work demonstrates that the H7N9 viruses could readily jump between avian species, which may have contributed to the evolution of the virus and its spread in the region.

Testing the genomic stability of the Brazilian yellow fever vaccine strain using next-generation sequencing data

The live attenuated yellow fever (YF) vaccine was developed in the 1930s. Currently, the 17D and 17DD attenuated substrains are used for vaccine production. The 17D strain is used for vaccine production by several countries, while the 17DD strain is used exclusively in Brazil. The cell passages carried out through the seed-lot system of vaccine production influence the presence of quasispecies causing changes in the stability and immunogenicity of attenuated genotypes by increasing attenuation or virulence. Using next-generation sequencing, we carried out genomic characterization and genetic diversity analysis between vaccine lots of the Brazilian YF vaccine, produced by BioManguinhos–Fiocruz, and used during 11 years of vaccination in Brazil. We present 20 assembled and annotated genomes from the Brazilian 17DD vaccine strain, eight single nucleotide polymorphisms and the quasispecies spectrum reconstruction for the 17DD vaccine, through a pipeline here introduced. The V2IDA pipeline provided a relationship between low genetic diversity, maintained through the seed lot system, and the confirmation of genetic stability of lots of the Brazilian vaccine against YF. Our study sets precedents for use of V2IDA in genetic diversity analysis and in silico stability investigation of attenuated viral vaccines, facilitating genetic surveillance during the vaccine production process.

Near-Complete Genome Sequences of Five Siciniviruses from North America

Here, we report near-complete genome sequences of Sicinivirus from U.S. poultry flocks in 2003 to 2005 and Mexico in 2019. They show highest nucleotide identity (84.5 to 85.5%) with other members of the Sicinivirus genus. These sequences update knowledge on diversity and contribute to a better understanding of the molecular epidemiology of Sicinivirus.

Here, we report near-complete genome sequences of sicinivirus from U.S. poultry flocks in 2003 to 2005 and Mexico in 2019. They show highest nucleotide identity (84.5 to 85.5%) with other members of the Sicinivirus genus. These sequences update knowledge on diversity and contribute to a better understanding of the molecular epidemiology of sicinivirus.

Whole-Genome Sequence of Avian coronavirus from a 15-Year-Old Sample Confirms Evidence of GA08-like Strain Circulation 4 Years Prior to Its First Reported Outbreak

Here, we report the complete genome sequence of an Avian coronavirus strain GA08-like isolate from a fecal sample from a broiler chicken collected in Georgia, USA, in 2004. The viral genome in this 15-year-old sample provides evidence for the circulation of the GA08-like strain at least 4 years before its first report in 2008.

Here, we report the complete genome sequence of an Avian coronavirus strain GA08-like isolate from a fecal sample from a broiler chicken collected in Georgia in 2004. The viral genome in this 15-year-old sample provides evidence for the circulation of the GA08-like strain at least 4 years before its first report in 2008.

Protection against Different Genotypes of Newcastle Disease Viruses (NDV) Afforded by an Adenovirus-Vectored Fusion Protein and Live NDV Vaccines in Chickens

The efficacy of an adenovirus-vectored Newcastle disease virus (NDV) vaccine expressing the fusion (F) NDV protein (adeno-F) was evaluated against challenges with virulent heterologous and homologous NDV strains to the F protein. In a preliminary study, two different doses (low and high) of adeno-F were tested against a virulent NDV strain containing the homologous NDV F protein, CA02. In a second study, at three weeks post-vaccination, the efficacy of the high dose of adeno-F was compared to a live attenuated NDV vaccine strain (LaSota) against three antigenically distinct virulent NDV challenge strains, one homologous (CA02) and two heterologous (TZ12, EG14) to F in the vectored vaccine. In both experiments, clinical signs, mortality, virus shedding, and humoral response were evaluated. In the first experiment, the survival rates from birds vaccinated with adeno-F at a high and low dose were 100% and 25%, respectively. In the second experiment, birds vaccinated with the high dose of adeno-F had a survival rate of 80%, 75%, and 65% after challenge with the CA02, TZ12, and EG14 viruses, respectively. All of the LaSota-vaccinated birds survived post-challenge no matter the NDV challenge strain. High antibody titers were detected after vaccination with LaSota by HI and ELISA tests. The majority of adeno-F-vaccinated birds had detectable antibodies using the ELISA test, but not using the HI test, before the challenge. The data show that both the similarity of the F protein of the adeno-F vaccine to the challenge virus and the adeno-F vaccination dose affect the efficacy of an adenovirus-vectored NDV vaccine against a virulent NDV challenge.

Surveillance and Genetic Characterization of Virulent Newcastle Disease Virus Subgenotype V.3 in Indigenous Chickens from Backyard Poultry Farms and Live Bird Markets in Kenya

Kenyan poultry consists of ~80% free-range indigenous chickens kept in small flocks (~30 birds) on backyard poultry farms (BPFs) and they are traded via live bird markets (LBMs). Newcastle disease virus (NDV) was detected in samples collected from chickens, wild farm birds, and other domestic poultry species during a 2017–2018 survey conducted at 66 BPFs and 21 LBMs in nine Kenyan counties. NDV nucleic acids were detected by rRT-PCR L-test in 39.5% (641/1621) of 1621 analyzed samples, of which 9.67% (62/641) were NDV-positive by both the L-test and a fusion-test designed to identify the virulent virus, with a majority being at LBMs (64.5%; 40/62) compared to BPFs (25.5%; 22/62). Virus isolation and next-generation sequencing (NGS) on a subset of samples resulted in 32 complete NDV genome sequences with 95.8–100% nucleotide identities amongst themselves and 95.7-98.2% identity with other east African isolates from 2010-2016. These isolates were classified as a new sub-genotype, V.3, and shared 86.5–88.9% and 88.5–91.8% nucleotide identities with subgenotypes V.1 and V.2 viruses, respectively. The putative fusion protein cleavage site (¹¹³R-Q-K-R↓F ¹¹⁷) in all 32 isolates, and a 1.86 ICPI score of an isolate from a BPF chicken that had clinical signs consistent with Newcastle disease, confirmed the high virulence of the NDVs. Compared to genotypes V and VI viruses, the attachment (HN) protein of 18 of the 32 vNDVs had amino acid substitutions in the antigenic sites. A time-scaled phylogeographic analysis suggests a west-to-east dispersal of the NDVs via the live chicken trade, but the virus origins remain unconfirmed due to scarcity of continuous and systematic surveillance data. This study reveals the widespread prevalence of vNDVs in Kenyan backyard poultry, the central role of LBMs in the dispersal and possibly generation of new virus variants, and the need for robust molecular epidemiological surveillance in poultry and non-poultry avian species.

A 24-Year-Old Sample Contributes the Complete Genome Sequence of Fowl Aviadenovirus D from the United States

Here, we report the complete genome sequence of fowl aviadenovirus D (FAdV-D) isolated from a preserved 24-year-old pancreas sample of a broiler chicken embryo. The results of the sequence showed that the viral genome is 44,079 bp long.

Runting and Stunting Syndrome in Broiler Chickens: Histopathology and Association With a Novel Picornavirus

Runting stunting syndrome (RSS) in commercial chickens has been reported worldwide, and although several studies have attempted to clarify the cause and describe the lesions, there are gaps in knowledge of the epidemiology, pathogenesis, and etiology. The study objective was to use commercial chicks naturally affected by RSS to describe the histologic changes of RSS in all segments of the small intestine in chicks of different ages and to identify viral gene sequences in affected chicks and their association with histologic lesions. Chicks lacking clinical signs but from the same houses and from unaffected houses were used as controls. The average weight of affected chicks was significantly lower than expected for their flocks. Macroscopically, the small intestines had paler serosa, with watery, mucoid, or foamy contents and poorly digested food. Histologic lesions were characterized by necrotic crypts, crypt dilation, and flattening of the crypt epithelium. Histomorphometry of the intestines revealed villous atrophy especially in the jejunum and ileum. Histologic changes in other organs were not observed. Random next-generation sequencing of total RNA extracted from formalin-fixed paraffin-embedded tissues detected avian nephritis virus, avian rotavirus, and picornavirus in jejunal segments from 7-day-old chicks. No viruses were detected in the jejunum of 1-day-old chicks. Detection of picornaviral reads was significantly associated (P < .05) with histologic lesions of RSS. Sequence analysis of the picornavirus revealed genetic similarity with the genus Gallivirus. Using in situ hybridization for galliviral nucleic acid sequences, the signal was associated with crypt lesion severity, although signal was detected both in chicks with and without RSS.

Evolution of H9N2 avian influenza viruses in Iran, 2017-2019

Since its first detection in 1998, avian influenza virus (AIV) subtype H9N2 has been enzootic in Iran. To better understand the evolutionary history of H9N2 viruses in Iran, we sequenced 15 currently circulating H9N2 viruses from domestic poultry during 2017-2019 and performed phylogenetic analysis of complete genome sequences. Phylogenetic analyses indicated that the Iranian H9N2 viruses formed multiple well-supported monophyletic groups within the G1-lineage of H9N2 virus. Our analysis of viral population dynamics revealed an increase in genetic diversity until 2007, corresponding to the multiple introductions and diversification of H9N2 viruses into multiple genetic groups (named Iran 1-4 subgroups), followed by a sudden decrease after 2008. Only the Iran 4 subgroup has survived, expanded, and currently circulates in Iran. The H9N2 viruses possessed many molecular markers associated with mammalian adaption in all gene segments, except neuraminidase gene. Considering the presence of mammalian host-specific markers, the public health threat of H9N2 viruses continues. Molecular analysis showed that Iranian H9N2 strains have continued to evolve and recent strains have multiple amino acid changes and addition of potential N-glycosylation on the antigenic sites of haemagglutinin. Continued antigenic and molecular surveillance of H9N2 viruses in poultry and mammals would be required to monitor further increments in viral evolution and their potential threat to public health.

Complete genome sequence of a subgenotype XXI.1.1 pigeon paramyxovirus type 1 virus (PPMV‑1) isolated from Iran in 2018 and phylogenetic analysis of a possible novel, but unassigned, PPMV-1 group isolated in 2014

Newcastle disease (ND) is one of the most serious infectious and contagious viral diseases in avian species. Recently, several ND outbreaks in pigeon caused by pigeon paramyxovirus serotype-1 (PPMV-1) have been reported from Iran, but unfortunately, phylogenetic studies have been mostly conducted on partial sequence of NDV fusion (F) gene. In addition, a complete genome data of Iranian PPMV-1 strains are not available. In the present study, a PPMV-1, named Avian avulavirus 1/pigeon/Iran/UT-EGV/2018, isolated from an infected pigeon, was subjected to whole-genome sequencing. The isolate showed an MDT of 74 h, thus categorizing it as mesogenic. The phylogenetic analysis based on the F gene sequence revealed the isolate belongs to XXI.1.1 subgenotype (min 0.9 % and max 3 %). To our knowledge, our study is the first study to publish the complete genome of a PPMV-1 from Iran. According to BLAST results, the whole genome of UT-EGV had high homology with some Russian, Egyptian and Ukrainian strains (the highest was 96.55 %). Additionally, we conducted a phylogenetic analysis on five PPMV-1 that we isolated in 2014 to find that they may belong to a completely unreported subgenotype (6 % distance when compared as a group). The information obtained from this study can be useful in preventive measures, including constructing an effective vaccine against PPMV-1 in Iran.

Selection and antigenic characterization of immune-escape mutants of H7N2 low pathogenic avian influenza virus using homologous polyclonal sera

Understanding the dynamics of the selection of influenza A immune escape variants by serum antibody is critical for designing effective vaccination programs for animals, especially poultry where large populations have a short generation time and may be vaccinated with high frequency. In this report, immune-escape mutants of A/turkey/New York/4450/1994 H7N2 low pathogenic avian influenza virus, were selected by serially passaging the virus in the presence of continuously increasing concentrations of homologous chicken polyclonal sera. Amino acid mutations were identified by sequencing the parental hemagglutinin (HA) gene and every 10 passages by both Sanger and deep sequencing, and the antigenic distance of the mutants to the parent strain was determined. Progressively, a total of five amino acid mutations were observed over the course of 30 passages. Based on their absence from the parental virus with deep sequencing, the mutations appear to have developed de novo. The antigenic distance between the selected mutants and the parent strain increased as the number of amino acid mutations accumulated and the concentration of antibodies had to be periodically increased to maintain the same reduction in virus titer during selection. This selection system demonstrates how H7 avian influenza viruses behave under selection with homologous sera, and provides a glimpse of their evolutionary dynamics, which can be applied to developing vaccination programs that maximize the effectiveness of a vaccine over time.

Complete Coding Sequences of Three Chicken Parvovirus Isolates from the United States

Parvoviruses are commonly found in U.S. poultry and are associated with clinical disease. Here, we report the complete coding sequences of three chicken parvoviruses from broiler chickens from commercial farms in the state of Georgia.

Parvoviruses are commonly found in U.S. poultry and are associated with clinical disease. Here, we report the complete coding sequences of three chicken parvoviruses from broiler chickens from commercial farms in the state of Georgia.

Molecular Characterization of Newcastle Disease Viruses Isolated from Chickens in Tanzania and Ghana

Newcastle disease (ND) is one of the most challenging infectious diseases affecting poultry production in Africa, causing major economic losses. To date, Newcastle disease virus isolates from several African countries have been grouped into class II NDV genotypes I, IV, V, VI, VII, XI, XIII, XIV, XVII, XVIII and XXI. Although ND is endemic in many African countries, information on circulating genotypes is still scarce. In Tanzania, outbreaks with genotypes V and XIII have been reported. In West and Central Africa, genotypes XIV, XVII, and XVIII are the most predominant. To investigate other genotypes circulating in Tanzania and Ghana, we performed molecular genotyping on isolates from Tanzania and Ghana using the MinION, a third-generation portable sequencing device from Oxford Nanopore Technologies. Using the MinION, we successfully sequenced the NDV F gene hypervariable region of 24 isolates from Tanzania and four samples from Ghana. In Tanzania, genotypes V, VII and XIII were detected. All isolates from Ghana belonged to genotype XVIII. The data obtained in this study reflect the genetic diversity of NDV in Africa and highlight the importance of surveillance for monitoring the distribution of NDV genotypes and viral evolution.

Infectious Bronchitis Virus Evolution, Diagnosis and Control

RNA viruses are characterized by high mutation and recombination rates, which allow a rapid adaptation to new environments. Most of the emerging diseases and host jumps are therefore sustained by these viruses. Rapid evolution may also hinder the understanding of molecular epidemiology, affect the sensitivity of diagnostic assays, limit the vaccine efficacy and favor episodes of immune escape, thus significantly complicating the control of even well-known pathogens. The history of infectious bronchitis virus (IBV) fits well with the above-mentioned scenario. Despite being known since the 1930s, it still represents one of the main causes of disease and economic losses for the poultry industry. A plethora of strategies have been developed and applied over time, with variable success, to limit its impact. However, they have rarely been evaluated objectively and on an adequate scale. Therefore, the actual advantages and disadvantages of IBV detection and control strategies, as well as their implementation, still largely depend on individual sensibility. The present manuscript aims to review the main features of IBV biology and evolution, focusing on their relevance and potential applications in terms of diagnosis and control.

Complete Genome Sequence of an Avian Metapneumovirus Subtype B Strain from Hungary

Avian metapneumoviruses (aMPVs), which have been reported in many countries, cause acute upper respiratory tract disease in chickens and turkeys. Using next-generation sequencing, we report here the complete genome sequence of an aMPV subtype B strain that was isolated from a turkey in Hungary in 1989.

Avian metapneumoviruses (aMPVs), which have been reported in many countries, cause acute upper respiratory tract disease in chickens and turkeys. Using next-generation sequencing, we report here the complete genome sequence of an aMPV subtype B strain that was isolated from a turkey in Hungary in 1989.

Complete Genome Sequencing, Molecular Epidemiological, and Pathogenicity Analysis of Pigeon Paramyxoviruses Type 1 Isolated in Guangxi, China during 2012-2018

Newcastle disease is an important poultry disease that also affects Columbiform birds. The viruses adapted to pigeons and doves are referred to as pigeon paramyxoviruses 1 (PPMV-1). PPMV-1 are frequently isolated from pigeons worldwide and have the potential to cause disease in chickens. The complete genomes of 18 PPMV-1 isolated in China during 2012–2018 were sequenced by next-generation sequencing (NGS). Comprehensive phylogenetic analyses showed that five of the viruses belong to sub-genotype VI1.2.1.1.2.1 and 13 isolates belong to sub-genotype VI.2.1.1.2.2. The results demonstrate that these sub-genotypes have been predominant in China during the last decade. The viruses of these sub-genotypes have been independently maintained and continuously evolved for over 20 years, and differ significantly from those causing outbreaks worldwide during the 1980s to 2010s. The viral reservoir remains unknown and possibilities of the viruses being maintained in both pigeon farms and wild bird populations are viable. In vivo characterization of the isolates’ pathogenicity estimated mean death times between 62 and 114 h and intracerebral pathogenicity indices between 0.00 and 0.63. Cross-reactivity testing showed minor antigenic differences between the studied viruses and the genotype II LaSota vaccine. These data will facilitate PPMV-1 epidemiology studies, vaccine development, and control of Newcastle disease in pigeons and poultry.

Complete Genome Sequence of Avian Coronavirus Strain GA08 (GI-27 Lineage)

Avian coronavirus, also known as infectious bronchitis virus, is a highly contagious respiratory pathogen of chickens that is responsible for major economic losses to the poultry industry around the globe. Here, we report the complete genome sequence of strain GA08 of the GI-27 lineage, isolated from a fecal sample from a broiler chicken collected in Georgia in 2015.

Avian coronavirus, also known as infectious bronchitis virus, is a highly contagious respiratory pathogen of chickens that is responsible for major economic losses to the poultry industry around the globe. Here, we report the complete genome sequence of strain GA08 of the GI-27 lineage, isolated from a fecal sample from a broiler chicken collected in Georgia in 2015.

A 25-Year-Old Sample Contributes the Complete Genome Sequence of Avian Coronavirus Vaccine Strain ArkDPI, Reisolated from Commercial Broilers in the United States

Here, we report the complete genome sequence of Avian coronavirus strain ArkDPI of the GI-9 lineage, isolated from broiler chickens in North Georgia in 1994. This is the complete genome sequence of this vaccine strain, reisolated from broilers in the United States.

Here, we report the complete genome sequence of Avian coronavirus strain ArkDPI of the GI-9 lineage, isolated from broiler chickens in North Georgia in 1994. This is the complete genome sequence of this vaccine strain, reisolated from broilers in the United States.

Complete Genome Sequences of 11 Newcastle Disease Virus Isolates of Subgenotype VII.2 from Indonesia

We report the complete genome sequences of 11 virulent Newcastle disease viruses. The isolates were obtained from vaccinated broiler and layer chickens in three different provinces of Indonesia in 2013 and 2014. Phylogenetic analysis revealed that all isolates belong to subgenotype VII.2 in the class II cluster.

We report the complete genome sequences of 11 virulent Newcastle disease viruses. The isolates were obtained from vaccinated broiler and layer chickens in three different provinces of Indonesia in 2013 and 2014. Phylogenetic analysis revealed that all isolates belong to subgenotype VII.2 in the class II cluster.

Isolation and genetic characterization of virulent strains of avian paramyxovirus-1 from multiple avian species in Azad Jammu and Kashmir 2017-2018

Despite intensive vaccination, endemicity of Avian paramyxoviruses-1 (APMV-1) is a significant problem in developing countries in Africa, Middle East, and Asia. Given the importance of APMV-1 in poultry and multiple non-poultry avian species, it is important to continue surveillance programs, routine monitoring and characterization of field isolates in the region where viruses are endemic. The purpose of this study was to pathotyped and genetically characterized 21 APMV-1s isolated from multiple avian species reared in different regions of Azad Jammu and Kashmir (AJK). Phylogenetic analysis based on complete fusion (F) gene sequences showed that 17 APMV-1 isolates obtained from commercial poultry and backyard birds belonged to sub-genotype VIIi. Though, one pigeon-origin APMV-1 isolate was clustered in sub-genotype VIg and three in recently designated new sub-genotype VIm of genotype VI. The pigeon-origin isolates had the following two motifs ¹¹³-RKKR↓F-¹¹⁷ and ¹¹³-RQRR↓F-¹¹⁷, while all other isolates had the polybasic amino acid sequence ¹¹³-RQKR↓F-¹¹⁷ at the F-cleavage site, which is characteristic of virulent APMV-1 strains. These results are consistent with the five viruses that had intracerebral pathogenicity indices (ICPIs) of between 1.50 and 1.73, corresponding to a velogenic pathotype. The APMV-1s isolated from commercial poultry and backyard birds in this study showed low nucleotide distance (0.3-0.9%) and genetically closely related (> 97%) to viruses repeatedly isolated (2011-2017) from multiple avian species in other states of Pakistan. Strengthened surveillance programs in both commercial poultry and backyard flocks are needed to better assess the commercial-backyard bird interface and form a basis for evidence-based measures to limit and prevent APMV-1 transmission.

Genetic stability of a Newcastle disease virus vectored infectious laryngotracheitis virus vaccine after serial passages in chicken embryos

Previously, we have demonstrated that the recombinant Newcastle disease virus (NDV) expressing the infectious laryngotracheitis virus (ILTV) glycoprotein D (gD) conferred protection against both virulent NDV and ILTV challenges in chickens. In this study, we evaluated the genetic stability of the recombinant vaccine after eight serial passages in embryonated chicken eggs (ECE). The vaccine master seed virus at the original egg-passage level 3 (EP3) was diluted and passaged in three separate repetitions (A, B and C) in ECE eight times (EP4 to EP11). RT-PCR analysis of the vaccine seed and egg-passaged virus stocks showed that there was no detectable insertion/deletion in the ILTV gD insert region. Next-generation sequencing analysis of the EP3 and EP11 virus stocks confirmed their genome integrity and revealed a total of thirteen single-nucleotide polymorphisms (SNPs). However, none of these SNPs were located in the ILTV gD insert or any of the known critical biological determinant positions. Virological and immunofluorescent assays provided additional evidence that the EP11 virus stocks retained their growth kinetics, low pathogenicity, and robust level of gD expression comparable to that of the vaccine master seed virus. This indicated that the SNPs were non-detrimental sporadic mutations. These results demonstrated that the insertion of ILTV gD gene into the NDV LaSota backbone did not significantly affect the genetic stability of the recombinant virus and that the rLS/ILTV-gD virus is a safe and genetically stable vaccine candidate after at least eight serial passages in ECE.

Genetic characterization and pathogenesis of the first H9N2 low pathogenic avian influenza viruses isolated from chickens in Kenyan live bird markets

Poultry production plays an important role in the economy and livelihoods of rural households in Kenya. As part of a surveillance program, avian influenza virus (AIV)-specific real-time RT-PCR (RRT-PCR) was used to screen 282 oropharyngeal swabs collected from chickens at six live bird markets (LBMs) and 33 backyard poultry farms in Kenya and 8 positive samples were detected. Virus was isolated in eggs from five samples, sequenced, and identified as H9N2 low pathogenic AIV (LPAIV) G1 lineage, with highest nucleotide sequence identity (98.6-99.9%) to a 2017 Ugandan H9N2 isolate. The H9N2 contained molecular markers for mammalian receptor specificity, implying their zoonotic potential. Virus pathogenesis and transmissibility was assessed by inoculating low and medium virus doses of a representative Kenyan H9N2 LPAIV isolate into experimental chickens and exposing them to naïve uninfected chickens at 2 -days post inoculation (dpi). Virus shedding was determined at 2/4/7 dpi and 2/5 days post placement (dpp), and seroconversion determined at 14 dpi/12 dpp. None of the directly-inoculated or contact birds exhibited any mortality or clinical disease signs. All directly-inoculated birds in the low dose group shed virus during the experiment, while only one contact bird shed virus at 2 dpp. Only two directly-inoculated birds that shed high virus titers seroconverted in that group. All birds in the medium dose group shed virus at 4/7 dpi and at 5 dpp, and they all seroconverted at 12/14 dpp. This is the first reported detection of H9N2 LPAIV from Kenya and it was shown to be infectious and transmissible in chickens by direct contact and represents a new disease threat to poultry and potentially to people.

Rapid evolution of Mexican H7N3 highly pathogenic avian influenza viruses in poultry

Highly pathogenic avian influenza (HPAI) virus subtype H7N3 has been circulating in poultry in Mexico since 2012 and vaccination has been used to control the disease. In this study, eight Mexican H7N3 HPAI viruses from 2015–2017 were isolated and fully sequenced. No evidence of reassortment was detected with other avian influenza (AI) viruses, but phylogenetic analyses show divergence of all eight gene segments into three genetic clusters by 2015, with 94.94 to 98.78 percent nucleotide homology of the HA genes when compared to the index virus from 2012. The HA protein of viruses from each cluster showed a different number of basic amino acids (n = 5–7) in the cleavage site, and six different patterns at the predicted N-glycosylation sites. Comparison of the sequences of the Mexican lineage H7N3 HPAI viruses and American ancestral wild bird AI viruses to characterize the virus evolutionary dynamics showed that the nucleotide substitution rates in PB2, PB1, PA, HA, NP, and NS genes greatly increased once the virus was introduced into poultry. The global nonsynonymous and synonymous ratios imply strong purifying selection driving the evolution of the virus. Forty-nine positively selected sites out of 171 nonsynonymous mutations were identified in the Mexican H7N3 HPAI viruses, including 7 amino acid changes observed in higher proportion in North American poultry origin AI viruses isolates than in wild bird-origin viruses. Continuous monitoring and molecular characterization of the H7N3 HPAI virus is important for better understanding of the virus evolutionary dynamics and further improving control measures including vaccination.

First Complete Genome Sequence of Currently Circulating Infectious Bronchitis Virus Strain DMV/1639 of the GI-17 Lineage

Avian infectious bronchitis virus is the causative agent of a highly contagious disease that results in severe economic losses to the poultry industry worldwide. Here, we report the first coding-complete genome sequence of strain DMV/1639 of the GI-17 lineage, isolated from broiler chickens in Georgia in 2019.

Avian infectious bronchitis virus is the causative agent of a highly contagious disease that results in severe economic losses to the poultry industry worldwide. Here, we report the first coding-complete genome sequence of strain DMV/1639 of the GI-17 lineage, isolated from broiler chickens in Georgia in 2019.

Pathogenicity and genomic changes of a 2016 European H5N8 highly pathogenic avian influenza virus (clade 2.3.4.4) in experimentally infected mallards and chickens

Highly pathogenic avian influenza H5N8 clade 2.3.4.4 virus caused outbreaks in poultry and unusually high mortality in wild birds in 2016-2017. The pathobiology of one of these viruses was examined in mallards and chickens. High mortality and transmission to direct contacts were observed in mallards inoculated with medium and high doses of the virus. However, in chickens, high mortality occurred only when birds are given the high virus dose and no transmission was observed, indicating that the virus was better adapted to mallards. In comparison with the virus inoculum, viral sequences obtained from the chickens had a higher number of nucleotide changes but lower intra-host genomic diversity than viral sequences obtained from the mallards. These observations are consistent with population bottlenecks occurring when viruses infect and replicate in a host that it is not well adapted to. Whether these observations apply to influenza viruses in general remains to be determined.

The Emergence of Avian Orthoavulavirus 13 in Wild Migratory Waterfowl in China Revealed the Existence of Diversified Trailer Region Sequences and HN Gene Lengths within this Serotype

Avian orthoavulavirus 13 (AOAV-13), also named avian paramyxovirus 13 (APMV-13), has been found sporadically in wild birds around the world ever since the discovery of AOAV-13 (AOAV-13/wild goose/Shimane/67/2000) in a wild goose from Japan in 2000. However, there are no reports of AOAV-13 in China. In the present study, a novel AOAV-13 virus (AOAV-13/wild goose/China/Hubei/V93-1/2015), isolated from a wild migratory waterfowl in a wetland of Hubei province of China, during active surveillance from 2013 to 2018, was biologically and genetically characterized. Phylogenetic analyses demonstrated a very close genetic relationship among all AOAV-13 strains, as revealed by very few genetic variations. Moreover, pathogenicity tests indicated that the V93-1 strain is a low virulent virus for chickens. However, the genome of the V93-1 virus was found to be 16,158 nucleotides (nt) in length, which is 12 nt or 162 nt longer than the other AOAV-13 strains that have been reported to date. The length difference of 12 nt in strain V93-1 is due to the existence of three repeats of the conserved sequence, “AAAAAT”, in the 5′-end trailer of the genome. Moreover, the HN gene of the V93-1 virus is 2070 nt in size, encoding 610 aa, which is the same size as the AOAV-13 strain from Japan, whereas that of two strains from Ukraine and Kazakhstan are 2080 nt in length, encoding 579 aa. We describe a novel AOAV-13 in migratory waterfowl in China, which suggests that diversified trailer region sequences and HN gene lengths exist within serotype AOAV-13, and highlight the need for its constant surveillance in poultry from live animal markets, and especially migratory birds.

First Complete Genome Sequence of a Subgenotype Vd Newcastle Disease Virus Isolate

A Newcastle disease virus was isolated from a chicken from a live bird market in the Mbeya region of Tanzania. Complete genome characterization of the isolate identified it as a member of subgenotype Vd. This is the first complete genome sequence of this subgenotype.

A Newcastle disease virus was isolated from a chicken from a live bird market in the Mbeya region of Tanzania. Complete genome characterization of the isolate identified it as a member of subgenotype Vd. This is the first complete genome sequence of this subgenotype.

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

Several Avian paramyxoviruses 1 (synonymous with Newcastle disease virus or NDV, used hereafter) classification systems have been proposed for strain identification and differentiation. These systems pioneered classification efforts; however, they were based on different approaches and lacked objective criteria for the differentiation of isolates. These differences have created discrepancies among systems, rendering discussions and comparisons across studies difficult. Although a system that used objective classification criteria was proposed by Diel and co-workers in 2012, the ample worldwide circulation and constant evolution of NDV, and utilization of only some of the criteria, led to identical naming and/or incorrect assigning of new sub/genotypes. To address these issues, an international consortium of experts was convened to undertake in-depth analyses of NDV genetic diversity. This consortium generated curated, up-to-date, complete fusion gene class I and class II datasets of all known NDV for public use, performed comprehensive phylogenetic neighbor-Joining, maximum-likelihood, Bayesian and nucleotide distance analyses, and compared these inference methods. An updated NDV classification and nomenclature system that incorporates phylogenetic topology, genetic distances, branch support, and epidemiological independence was developed. This new consensus system maintains two NDV classes and existing genotypes, identifies three new class II genotypes, and reduces the number of sub-genotypes. In order to track the ancestry of viruses, a dichotomous naming system for designating sub-genotypes was introduced. In addition, a pilot dataset and sub-trees rooting guidelines for rapid preliminary genotype identification of new isolates are provided. Guidelines for sequence dataset curation and phylogenetic inference, and a detailed comparison between the updated and previous systems are included. To increase the speed of phylogenetic inference and ensure consistency between laboratories, detailed guidelines for the use of a supercomputer are also provided. The proposed unified classification system will facilitate future studies of NDV evolution and epidemiology, and comparison of results obtained across the world.

Virulent Newcastle disease viruses from chicken origin are more pathogenic and transmissible to chickens than viruses normally maintained in wild birds

Five, class II, virulent Newcastle disease virus (vNDV) isolates of different genotypes from different host species were evaluated for their ability to infect, cause disease, and transmit to naïve chickens. Groups of five birds received a low, medium, or high dose, by the oculonasal route, of one of the following vNDV: three chicken-origin, one cormorant-origin, and one pigeon-origin. Three naïve birds were added to each group at two days post-inoculation (DPI) to evaluate transmission. Virus shedding was quantified from swabs (2/4/7 DPI), and seroconversion was evaluated at 14 DPI. All inoculated and contact birds in the chicken-origin vNDV groups succumbed to infection, displaying clinical signs typical of Newcastle disease and shed virus titers above 6 log₁₀ EID₅₀/ml. Birds receiving a high and medium dose of the cormorant virus showed primarily neurological clinical signs with 80% and 60% mortality, respectively. The chickens showing clinical disease shed virus at titers below 4 log₁₀ EID₅₀/ml, and the remaining bird in the high dose group seroconverted with a high HI titer. For the pigeon-origin virus, no clinical signs were observed in any of the birds, but all 5 chickens in the high challenge dose and one bird in the medium challenge group shed virus at mean titers of 3.1 and 2.2 log₁₀ EID₅₀/ml, respectively. Overall, the chicken-origin viruses infected chickens and efficiently transmitted to naïve birds, while the cormorant- and pigeon-origin viruses infected chickens only at the higher doses and did not transmit to other birds.

Global phylodynamic analysis of avian paramyxovirus-1 provides evidence of inter-host transmission and intercontinental spatial diffusion

BACKGROUND

Avian avulavirus (commonly known as avian paramyxovirus-1 or APMV-1) can cause disease of varying severity in both domestic and wild birds. Understanding how viruses move among hosts and geography would be useful for informing prevention and control efforts. A Bayesian statistical framework was employed to estimate the evolutionary history of 1602 complete fusion gene APMV-1 sequences collected from 1970 to 2016 in order to infer viral transmission between avian host orders and diffusion among geographic regions. Ancestral states were estimated with a non-reversible continuous-time Markov chain model, allowing transition rates between discrete states to be calculated. The evolutionary analyses were stratified by APMV-1 classes I (n = 198) and II (n = 1404), and only those sequences collected between 2006 and 2016 were allowed to contribute host and location information to the viral migration networks.

RESULTS

While the current data was unable to assess impact of host domestication status on APMV-1 diffusion, these analyses supported the sharing of APMV-1 among divergent host taxa. The highest supported transition rate for both classes existed from domestic chickens to Anseriformes (class I:6.18 transitions/year, 95% highest posterior density (HPD) 0.31-20.02, Bayes factor (BF) = 367.2; class II:2.88 transitions/year, 95%HPD 1.9-4.06, BF = 34,582.9). Further, among class II viruses, domestic chickens also acted as a source for Columbiformes (BF = 34,582.9), other Galliformes (BF = 34,582.9), and Psittaciformes (BF = 34,582.9). Columbiformes was also a highly supported source to Anseriformes (BF = 322.0) and domestic chickens (BF = 402.6). Additionally, our results provide support for the diffusion of viruses among continents and regions, but no interhemispheric viral exchange between 2006 and 2016. Among class II viruses, the highest transition rates were estimated from South Asia to the Middle East (1.21 transitions/year; 95%HPD 0.36-2.45; BF = 67,107.8), from Europe to East Asia (1.17 transitions/year; 95%HPD 0.12-2.61; BF = 436.2) and from Europe to Africa (1.06 transitions/year, 95%HPD 0.07-2.51; BF = 169.3).

CONCLUSIONS

While migration appears to occur infrequently, geographic movement may be important in determining viral diversification and population structure. In contrast, inter-order transmission of APMV-1 may occur readily, but most events are transient with few lineages persisting in novel hosts.

Genomic comparison of Newcastle disease viruses isolated in Nigeria between 2002 and 2015 reveals circulation of highly diverse genotypes and spillover into wild birds

Newcastle disease virus (NDV) has a wide avian host range and a high degree of genetic variability, and virulent strains cause Newcastle disease (ND), a worldwide concern for poultry health. Although NDV has been studied in Nigeria, genetic information about the viruses involved in the endemicity of the disease and the transmission that likely occurs at the poultry-wildlife interface is still largely incomplete. Next-generation and Sanger sequencing was performed to provide complete (n = 73) and partial genomic sequence data (n = 38) for NDV isolates collected from domestic and wild birds in Nigeria during 2002-2015, including the first complete genome sequences of genotype IV and subgenotype VIh from the African continent. Phylogenetic analysis revealed that viruses of seven different genotypes circulated in that period, demonstrating high genetic diversity of NDV for a single country. In addition, a high degree of similarity between NDV isolates from domestic and wild birds was observed, suggesting that spillovers had occurred, including to three species that had not previously been shown to be susceptible to NDV infection. Furthermore, the first spillover of a mesogenic Komarov vaccine virus is documented, suggesting a previous spillover and evolution of this virus. The similarities between viruses from poultry and multiple bird species and the lack of evidence for host adaptation in codon usage suggest that transmission of NDV between poultry and non-poultry birds occurred recently. This is especially significant when considering that some viruses were isolated from species of conservation concern. The high diversity of NDV observed in both domestic and wild birds in Nigeria emphasizes the need for active surveillance and epidemiology of NDV in all bird species.

Enhanced phylogenetic resolution of Newcastle disease outbreaks using complete viral genome sequences from formalin-fixed paraffin-embedded tissue samples

Highly virulent Newcastle disease virus (NDV) causes Newcastle disease (ND), which is a threat to poultry production worldwide. Effective disease management requires approaches to accurately determine sources of infection, which involves tracking of closely related viruses. Next-generation sequencing (NGS) has emerged as a research tool for thorough genetic characterization of infectious organisms. Previously formalin-fixed paraffin-embedded (FFPE) tissues have been used to conduct retrospective epidemiological studies of related but genetically distinct viruses. However, this study extends the applicability of NGS for complete genome analysis of viruses from FFPE tissues to track the evolution of closely related viruses. Total RNA was obtained from FFPE spleens, lungs, brains, and small intestines of chickens in 11 poultry flocks during disease outbreaks in Pakistan. The RNA was randomly sequenced on an Illumina MiSeq instrument and the raw data were analyzed using a custom data analysis pipeline that includes de novo assembly. Genomes of virulent NDV were detected in 10/11 birds: eight nearly complete (> 95% coverage of concatenated coding sequence) and two partial genomes. Phylogeny of the NDV complete genome coding sequences was compared to current methods of analysis based on the full and partial fusion genes and determined that the approach provided a better phylogenetic resolution. Two distinct lineages of sub-genotype VIIi NDV were identified to be simultaneously circulating in Pakistani poultry. Non-targeted NGS of total RNA from FFPE tissues coupled with de novo assembly provided a reliable, safe, and affordable method to conduct epidemiological and evolutionary studies to facilitate management of ND in Pakistan.

Intercontinental spread of Asian-origin H7 avian influenza viruses by captive bird trade in 1990's

Wild bird migration and illegal trade of infected poultry, eggs, and poultry products have been associated with the spread of avian influenza viruses (AIV). During 1992-1996, H7N1 and H7N8 low pathogenic AIV (LPAIV) were identified from captive wild birds; such as Pekin robin (Leiothrix lutea), magpie robin (Copsychus saularis), flycatcher sp. (genus Empidonax), a species of softbill and parakeet, sun conure (Aratinga solstitialis), painted conure (Pyrrhura picta), fairy bluebird (Irena puella), and common iora (Aegithina tiphia), kept in aviaries or quarantine stations in England, The Netherlands, Singapore and the United States (U.S.). In this study, we sequenced these H7 viruses isolated from quarantine facilities and aviaries using next-generation sequencing and conducted a comparative phylogenetic analysis of complete genome sequences to elucidate spread patterns. The complete genome sequencing and phylogenetic analysis suggested that H7 viruses originated from a common source, even though they were obtained from birds in distant geographical regions. All H7N1 and H7N8 viruses were LPAIV, except a H7N1 highly pathogenic AIV (HPAIV), A/Pekin robin/California/30412/1994(H7N1) virus. Our results support the continued need for regulation of the captive wild bird trade to reduce the risk of introduction and dissemination of both LPAIV and HPAIV throughout the world.