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.

Complete Genome Sequence of an Avian Orthoavulavirus 13 Strain Detected in Ukraine.

We report the complete genome sequence of an avian orthoavulavirus 13 strain, isolated from a white-fronted goose in the Odesa region of Ukraine in 2013. The detection of avian orthoavulavirus 13 in Ukraine confirms that the geographic distribution of this virus extends beyond Asia.

Complete Genome Sequences of Avian Metapneumovirus Subtype B Vaccine Strains from Brazil.

Avian metapneumovirus (aMPV) causes a highly contagious upper respiratory and reproductive disease in chickens, turkeys, and ducks. Here, complete genome sequences of aMPV-B vaccine strains BR/1890/E1/19 (PL21, Nemovac; Boehringer Ingelheim Animal Health, Brazil) and BR/1891/E2/19 (1062; Hipraviar, France) were sequenced and compared with the pathogenic field strain VCO3/60616.

Unique Variants of Avian Coronaviruses from Indigenous Chickens in Kenya.

The avian gamma-coronavirus infectious bronchitis virus (AvCoV, IBV; Coronaviridae family) causes upper respiratory disease associated with severe economic losses in the poultry industry worldwide. Here, we report for the first time in Kenya and the Eastern African region two novel AvCoVs, designated IBV/ck/KE/1920/A374/2017 (A374/17) and AvCoV/ck/KE/1922/A376/2017 (A376/17), inadvertently discovered using random nontargeted next-generation sequencing (NGS) of cloacal swabs collected from indigenous chickens. Despite having genome organization (5'UTR-[Rep1a/1ab-S-3a-3b-E-M-4b-4c-5a-5b-N-6b]-3'UTR), canonical conservation of essential genes and size (~27.6 kb) typical of IBVs, the Kenyan isolates do not phylogenetically cluster with any genotypes of the 37 IBV lineages and 26 unique variants (UVs). Excluding the spike gene, genome sequences of A374/17 and A376/17 are only 93.1% similar to each other and 86.7-91.4% identical to genomes of other AvCoVs. All five non-spike genes of the two isolates phylogenetically cluster together and distinctly from other IBVs and turkey coronaviruses (TCoVs), including the indigenous African GI-26 viruses, suggesting a common origin of the genome backbone of the Kenyan isolates. However, isolate A376/17 contains a TCoV-like spike (S) protein coding sequence and is most similar to Asian TCoVs (84.5-85.1%) compared to other TCoVs (75.6-78.5%), whereas isolate A374/17 contains an S1 gene sequence most similar to the globally distributed lineage GI-16 (78.4-79.5%) and the Middle Eastern lineage GI-23 (79.8-80.2%) viruses. Unanswered questions include the actual origin of the Kenyan AvCoVs, the potential pathobiological significance of their genetic variations, whether they have indeed established themselves as independent variants and subsequently spread within Kenya and to the neighboring east/central African countries that have porous live poultry trade borders, and whether the live-attenuated Mass-type (lineage GI-1)-based vaccines currently used in Kenya and most of the African countries provide protection against these genetically divergent field variants.

Genetic diversity of Newcastle disease viruses circulating in wild and synanthropic birds in Ukraine between 2006 and 2015.

Newcastle disease virus (NDV) infects a wide range of bird species worldwide and is of importance to the poultry industry. Although certain virus genotypes are clearly associated with wild bird species, the role of those species in the movement of viruses and the migratory routes they follow is still unclear. In this study, we performed a phylogenetic analysis of nineteen NDV sequences that were identified among 21,924 samples collected from wild and synanthropic birds from different regions of Ukraine from 2006 to 2015 and compared them with isolates from other continents. In synanthropic birds, NDV strains of genotype II, VI, VII, and XXI of class II were detected. The fusion gene sequences of these strains were similar to strains detected in birds from different geographical regions of Europe and Asia. However, it is noteworthy to mention the isolation of vaccine viruses from synanthropic birds, suggesting the possibility of their role in viral transmission from vaccinated poultry to wild birds, which may lead to the further spreading of vaccine viruses into other regions during wild bird migration. Moreover, here we present the first publicly available complete NDV F gene from a crow (genus Corvus). Additionally, our phylogenetic results indicated a possible connection of Ukrainian NDV isolates with genotype XXI strains circulating in Kazakhstan. Among strains from wild birds, NDVs of genotype 1 of class I and genotype I of class II were detected. The phylogenetic analysis highlighted the possible exchange of these NDV strains between wild waterfowl from the Azov-Black Sea region of Ukraine and waterfowl from different continents, including Europe, Asia, and Africa.

Comparable outcomes from long and short read random sequencing of total RNA for detection of pathogens in chicken respiratory samples.

Co-infections of avian species with different RNA viruses and pathogenic bacteria are often misdiagnosed or incompletely characterized using targeted diagnostic methods, which could affect the accurate management of clinical disease. A non-targeted sequencing approach with rapid and precise characterization of pathogens should help respiratory disease management by providing a comprehensive view of the causes of disease. Long-read portable sequencers have significant potential advantages over established short-read sequencers due to portability, speed, and lower cost. The applicability of short reads random sequencing for direct detection of pathogens in clinical poultry samples has been previously demonstrated. Here we demonstrate the feasibility of long read random sequencing approaches to identify disease agents in clinical samples. Experimental oropharyngeal swab samples (n = 12) from chickens infected with infectious bronchitis virus (IBV), avian influenza virus (AIV) and Mycoplasma synoviae (MS) and field-collected clinical oropharyngeal swab samples (n = 11) from Kenyan live bird markets previously testing positive for Newcastle disease virus (NDV) were randomly sequenced on the MinION platform and results validated by comparing to real time PCR and short read random sequencing in the Illumina MiSeq platform. In the swabs from experimental infections, each of three agents in every RT-qPCR-positive sample (Ct range 19-34) was detectable within 1 h on the MinION platform, except for AIV one agent in one sample (Ct = 36.21). Nine of 12 IBV-positive samples were assigned genotypes within 1 h, as were five of 11 AIV-positive samples. MinION relative abundances of the test agent (AIV, IBV and MS) were highly correlated with RT-qPCR Ct values (R range-0.82 to-0.98). In field-collected clinical swab samples, NDV (Ct range 12-37) was detected in all eleven samples within 1 h of MinION sequencing, with 10 of 11 samples accurately genotyped within 1 h. All NDV-positive field samples were found to be co-infected with one or more additional respiratory agents. These results demonstrate that MinION sequencing can provide rapid, and sensitive non-targeted detection and genetic characterization of co-existing respiratory pathogens in clinical samples with similar performance to the Illumina MiSeq.

Detection and Genome Sequence Analysis of Avian Metapneumovirus Subtype A Viruses Circulating in Commercial Chicken Flocks in Mexico.

Avian metapneumoviruses (aMPV subtypes A-D) are respiratory and reproductive pathogens of poultry. Since aMPV-A was initially reported in Mexico in 2014, there have been no additional reports of its detection in the country. Using nontargeted next-generation sequencing (NGS) of FTA card-spotted respiratory samples from commercial chickens in Mexico, seven full genome sequences of aMPV-A (lengths of 13,288-13,381 nucleotides) were de novo assembled. Additionally, complete coding sequences of genes N (n = 2), P and M (n = 7 each), F and L (n = 1 each), M2 (n = 6), SH (n = 5) and G (n = 2) were reference-based assembled from another seven samples. The Mexican isolates phylogenetically group with, but in a distinct clade separate from, other aMPV-A strains. The genome and G-gene nt sequences of the Mexican aMPVs are closest to strain UK/8544/06 (97.22-97.47% and 95.07-95.83%, respectively). Various amino acid variations distinguish the Mexican isolates from each other, and other aMPV-A strains, most of which are in the G (n = 38), F (n = 12), and L (n = 19) proteins. Using our sequence data and publicly available aMPV-A data, we revised a previously published rRT-PCR test, which resulted in different cycling and amplification conditions for aMPV-A to make it more compatible with other commonly used rRT-PCR diagnostic cycling conditions. This is the first comprehensive sequence analysis of aMPVs in Mexico and demonstrates the value of nontargeted NGS to identify pathogens where targeted virus surveillance is likely not routinely performed.

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.

Novel Recombinant Newcastle Disease Virus-Based In Ovo Vaccines Bypass Maternal Immunity to Provide Full Protection from Early Virulent Challenge.

Newcastle disease (ND) is one of the most economically important poultry diseases. Despite intensive efforts with current vaccination programs, this disease still occurs worldwide, causing significant mortality even in vaccinated flocks. This has been partially attributed to a gap in immunity during the post-hatch period due to the presence of maternal antibodies that negatively impact the replication of the commonly used live vaccines. In ovo vaccines have multiple advantages and present an opportunity to address this problem. Currently employed in ovo ND vaccines are recombinant herpesvirus of turkeys (HVT)-vectored vaccines expressing Newcastle disease virus (NDV) antigens. Although proven efficient, these vaccines have some limitations, such as delayed immunogenicity and the inability to administer a second HVT vaccine post-hatch. The use of live ND vaccines for in ovo vaccination is currently not applicable, as these are associated with high embryo mortality. In this study, recombinant NDV-vectored experimental vaccines containing an antisense sequence of avian interleukin 4 (IL4R) and their backbones were administered in ovo at different doses in 18-day-old commercial eggs possessing high maternal antibodies titers. The hatched birds were challenged with virulent NDV at 2 weeks-of-age. Post-hatch vaccine shedding, post-challenge survival, challenge virus shedding, and humoral immune responses were evaluated at multiple timepoints. Recombinant NDV (rNDV) vaccinated birds had significantly reduced post-hatch mortality compared with the wild-type LaSota vaccine. All rNDV vaccines were able to penetrate maternal immunity and induce a strong early humoral immune response. Further, the rNDV vaccines provided protection from clinical disease and significantly decreased virus shedding after early virulent NDV challenge at two weeks post-hatch. The post-challenge hemagglutination-inhibition antibody titers in the vaccinated groups remained comparable with the pre-challenge titers, suggesting the capacity of the studied vaccines to prevent efficient replication of the challenge virus. Post-hatch survival after vaccination with the rNDV-IL4R vaccines was dose-dependent, with an increase in survival as the dose decreased. This improved survival and the dose-dependency data suggest that novel attenuated in ovo rNDV-based vaccines that are able to penetrate maternal immunity to elicit a strong immune response as early as 14 days post-hatch, resulting in high or full protection from virulent challenge, show promise as a contributor to the control of Newcastle disease.

A Novel Recombinant Newcastle Disease Vaccine Improves Post- In Ovo Vaccination Survival with Sustained Protection against Virulent Challenge.

In ovo vaccination has been employed by the poultry industry for over 20 years to control numerous avian diseases. Unfortunately, in ovo live vaccines against Newcastle disease have significant limitations, including high embryo mortality and the inability to induce full protection during the first two weeks of life. In this study, a recombinant live attenuated Newcastle disease virus vaccine containing the antisense sequence of chicken interleukin 4 (IL-4), rZJ1*L-IL4R, was used. The rZJ1*L-IL4R vaccine was administered in ovo to naïve specific pathogen free embryonated chicken eggs (ECEs) and evaluated against a homologous challenge. Controls included a live attenuated recombinant genotype VII vaccine based on the virus ZJ1 (rZJ1*L) backbone, the LaSota vaccine and diluent alone. In the first of two experiments, ECEs were vaccinated at 18 days of embryonation (DOE) with either 104.5 or 103.5 50% embryo infectious dose (EID50/egg) and chickens were challenged at 21 days post-hatch (DPH). In the second experiment, 103.5 EID50/egg of each vaccine was administered at 19 DOE, and chickens were challenged at 14 DPH. Chickens vaccinated with 103.5 EID50/egg of rZJ1*L-IL4R had hatch rates comparable to the group that received diluent alone, whereas other groups had significantly lower hatch rates. All vaccinated chickens survived challenge without displaying clinical disease, had protective hemagglutination inhibition titers, and shed comparable levels of challenge virus. The recombinant rZJ1*L-IL4R vaccine yielded lower post-vaccination mortality rates compared with the other in ovo NDV live vaccine candidates as well as provided strong protection post-challenge.

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.

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

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 (113R-Q-K-R↓F 117) 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.

Virulence during Newcastle Disease Viruses Cross Species Adaptation.

The hypothesis that host adaptation in virulent Newcastle disease viruses (NDV) has been accompanied by virulence modulation is reviewed here. Historical records, experimental data, and phylogenetic analyses from available GenBank sequences suggest that currently circulating NDVs emerged in the 1920-1940's from low virulence viruses by mutation at the fusion protein cleavage site. These viruses later gave rise to multiple virulent genotypes by modulating virulence in opposite directions. Phylogenetic and pathotyping studies demonstrate that older virulent NDVs further evolved into chicken-adapted genotypes by increasing virulence (velogenic-viscerotropic pathotypes with intracerebral pathogenicity indexes [ICPIs] of 1.6 to 2), or into cormorant-adapted NDVs by moderating virulence (velogenic-neurotropic pathotypes with ICPIs of 1.4 to 1.6), or into pigeon-adapted viruses by further attenuating virulence (mesogenic pathotypes with ICPIs of 0.9 to 1.4). Pathogenesis and transmission experiments on adult chickens demonstrate that chicken-adapted velogenic-viscerotropic viruses are more capable of causing disease than older velogenic-neurotropic viruses. Currently circulating velogenic-viscerotropic viruses are also more capable of replicating and of being transmitted in naïve chickens than viruses from cormorants and pigeons. These evolutionary virulence changes are consistent with theories that predict that virulence may evolve in many directions in order to achieve maximum fitness, as determined by genetic and ecologic constraints.

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.

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.

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.

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.