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

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

Phylogenetic assessment reveals continuous evolution and circulation of pigeon-derived virulent avian avulaviruses 1 in Eastern Europe, Asia, and Africa

Abstract Background The remarkable diversity and mobility of Newcastle disease viruses (NDV) includes virulent viruses of genotype VI. These viruses are often referred to as pigeon paramyxoviruses 1 because they are normally isolated and cause clinical disease in birds from the Columbidae family. Genotype VI viruses occasionally infect, and may also cause clinical disease in poultry. Thus, the evolution, current spread and detection of these viruses are relevant to avian health. Results Here, we describe the isolation and genomic characterization of six Egyptian (2015), four Pakistani (2015), and two Ukrainian (2007, 2013) recent pigeon-derived NDV isolates of sub-genotype VIg. These viruses are closely related to isolates from Kazakhstan, Nigeria and Russia. In addition, eight genetically related NDV isolates from Pakistan (2014–2016) that define a new sub-genotype (VIm) are described. All of these viruses, and the ancestral Bulgarian (n = 2) and South Korean (n = 2) viruses described here, have predicted virulent cleavage sites of the fusion protein, and those selected for further characterization have intracerebral pathogenicity index assay values characteristic of NDV of genotype VI (1.31 to 1.48). A validated matrix gene real-time RT-PCR (rRT-PCR) NDV test detect all tested isolates. However, the validated rRT-PCR test that is normally used to identify the virulent fusion gene fails to detect the Egyptian and Ukrainian viruses due to mismatches in primers and probe. A new rapid rRT-PCR test to determine the presence of virulent cleavage sites for viruses from sub-genotypes VIg was developed and evaluated on these and other viruses. Conclusions We describe the almost simultaneous circulation and continuous evolution of genotype VI Newcastle disease viruses in distant locations, suggesting epidemiological connections among three continents. As pigeons are not migratory, this study suggests the need to understand the possible role of human activity in the dispersal of these viruses. Complete genomic characterization identified previously unrecognized genetic diversity that contributes to diagnostic failure and will facilitate future evolutionary studies. These results highlight the importance of conducting active surveillance on pigeons worldwide and the need to update existent rapid diagnostic protocols to detect emerging viral variants and help manage the disease in affected regions

Rapid virulence prediction and identification of Newcastle disease virus genotypes using third-generation sequencing

Abstract Background Newcastle disease (ND) outbreaks are global challenges to the poultry industry. Effective management requires rapid identification and virulence prediction of the circulating Newcastle disease viruses (NDV), the causative agent of ND. However, these diagnostics are hindered by the genetic diversity and rapid evolution of NDVs. Methods An amplicon sequencing (AmpSeq) workflow for virulence and genotype prediction of NDV samples using a third-generation, real-time DNA sequencing platform is described here. 1D MinION sequencing of barcoded NDV amplicons was performed using 33 egg-grown isolates, (15 NDV genotypes), and 15 clinical swab samples collected from field outbreaks. Assembly-based data analysis was performed in a customized, Galaxy-based AmpSeq workflow. MinION-based results were compared to previously published sequences and to sequences obtained using a previously published Illumina MiSeq workflow. Results For all egg-grown isolates, NDV was detected and virulence and genotype were accurately predicted. For clinical samples, NDV was detected in ten of eleven NDV samples. Six of the clinical samples contained two mixed genotypes as determined by MiSeq, of which the MinION method detected both genotypes in four samples. Additionally, testing a dilution series of one NDV isolate resulted in NDV detection in a dilution as low as 101 50% egg infectious dose per milliliter. This was accomplished in as little as 7 min of sequencing time, with a 98.37% sequence identity compared to the expected consensus obtained by MiSeq. Conclusion The depth of sequencing, fast sequencing capabilities, accuracy of the consensus sequences, and the low cost of multiplexing allowed for effective virulence prediction and genotype identification of NDVs currently circulating worldwide. The sensitivity of this protocol was preliminary tested using only one genotype. After more extensive evaluation of the sensitivity and specificity, this protocol will likely be applicable to the detection and characterization of NDV

Additional file 1: Table S1. of Phylogenetic assessment reveals continuous evolution and circulation of pigeon-derived virulent avian avulaviruses 1 in Eastern Europe, Asia, and Africa

Nucleotide sequences of primers used in PCR amplification, and sequencing of the NDV isolates used in this study. Table S2. List of the NDV used for construction of full fusion phylogenetic tree presented in Fig. 1. Highlighted in bold font are the viruses studied in the current work. Table S3. List of the NDV used for construction of complete genome phylogenetic tree presented in Additional file 2: Fig. S4. Highlighted in bold font are the viruses studied in the current work. Table S4. Characteristics of the thirteen complete genomes of Newcastle disease viruses of genotype VI sequenced in this study. Table S5. Estimated pairwise evolutionary distances among viruses of the new sub-genotype VIm. (DOCX 60 kb

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.Online Resource 1 Maximum composite likelihood tree with no collapsed branches, constructed using complete fusion gene coding sequences. Genotypes and subgenotypes of viruses are presented with Roman numerals and lowercase letters in each taxon name.Online Resource 2 Maximum composite likelihood tree constructed using complete genome coding sequences. Genotypes and subgenotypes of viruses are presented with Roman numerals and lowercase letters in each taxon name. Red lettering indicates viruses sequenced for this study.Online Resource 3 List of sequences used for the maximum composite likelihood tree constructed using complete fusion gene coding sequences (Fig. 2, Online Resource 1). Isolates indicated in bold were sequenced for this study.Online Resource 4 List of sequences used for the maximum composite likelihood tree constructed using complete gene coding sequences (Online Resource 2). Isolates indicated in bold were sequenced for this study.Online Resource 5 Complete genome relative synonymous codon usage (RSCU) values among poultry and wild bird viruses.Online Resource 6 Number of codons used in the reference Gallus gallus dataset as implemented in DAMBE.The Agricultural Research Service (ARS) and supported by the USDA Current Research Information System (CRIS) (number 6612-32000-072-00D) and partially funded by The Defense Threat Reduction Agency (DTRA) (FRCALL12-6-2-0015).http://link.springer.com/journal/7052020-08-01hj2019Production Animal Studie

Additional file 2: Figure S1 A, B and C. of Phylogenetic assessment reveals continuous evolution and circulation of pigeon-derived virulent avian avulaviruses 1 in Eastern Europe, Asia, and Africa

Mismatches between the tested viruses and: A) previous fusion probe designed by Wise et al. [24]; B) pigeon-specific fusion probe designed by Kim et al. [66]; and C) the optimized pigeon specific probe in this study, respectively. Sequences are in order of 5′ to 3′. Figure S2. A and B Mismatches between the tested viruses and: A) previous fusion forward primer designed by Wise et al. [24]; and B) the optimized fusion forward primer in this study, respectively. Sequences are in order of 5′ to 3′. Figure S3. A and B Mismatches between the tested viruses and: A) previous fusion reverse primers designed by Wise et al. [24]; B) and the new fusion reverse in this study, respectively. Sequences are in order of 5′ to 3′. Figure S4. Phylogenetic analysis based on the complete genome concatenated coding sequence of viruses representing NDV class II. Only bootstrap values greater or equal to 60% are visualized. There were a total of 13,697 positions in the final dataset. The strains sequenced in this study are highlighted in bold font and have a circle symbol in front the taxa name. Provisional designation of genotypes is indicated on the right. (DOCX 1990 kb

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

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