Molecular Characterization of Complete Genome Sequence of an Avian Coronavirus Identified in a Backyard Chicken from Tanzania

A complete genome sequence of an avian coronavirus (AvCoV; 27,663 bp excluding 3′ poly(A) tail) was determined using nontargeted next-generation sequencing (NGS) of an oropharyngeal swab from a backyard chicken in a live bird market in Arusha, Tanzania. The open reading frames (ORFs) of the Tanzanian strain TZ/CA127/19 are organized as typical of gammaCoVs (Coronaviridae family): 5′UTR-[ORFs 1a/1b encoding replicase complex (Rep1ab) non-structural peptides nsp2-16]-[spike (S) protein]-[ORFs 3a/3b]-[small envelop (E) protein]-[membrane (M) protein]-[ORFs 4a/4c]-[ORFs 5a/5b]-[nucleocapsid (N) protein]-[ORF6b]-3′UTR. The structural (S, E, M and N) and Rep1ab proteins of TZ/CA127/19 contain features typically conserved in AvCoVs, including the cleavage sites and functional motifs in Rep1ab and S. Its genome backbone (non-spike region) is closest to Asian GI-7 and GI-19 infectious bronchitis viruses (IBVs) with 87.2–89.7% nucleotide (nt) identities, but it has a S gene closest (98.9% nt identity) to the recombinant strain ck/CN/ahysx-1/16. Its 3a, 3b E and 4c sequences are closest to the duck CoV strain DK/GD/27/14 at 99.43%, 100%, 99.65% and 99.38% nt identities, respectively. Whereas its S gene phylogenetically cluster with North American TCoVs and French guineafowl COVs, all other viral genes group monophyletically with Eurasian GI-7/GI-19 IBVs and Chinese recombinant AvCoVs. Detection of a 4445 nt-long recombinant fragment with breakpoints at positions 19,961 and 24,405 (C- and N-terminus of nsp16 and E, respectively) strongly suggested that TZ/CA127/19 acquired its genome backbone from an LX4-type (GI-19) field strain via recombination with an unknown AvCoV. This is the first report of AvCoV in Tanzania and leaves unanswered the questions of its emergence and the biological significance

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

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

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

Abstract Background Next-generation sequencing (NGS) allows ultra-deep sequencing of nucleic acids. The use of sequence-independent amplification of viral nucleic acids without utilization of target-specific primers provides advantages over traditional sequencing methods and allows detection of unsuspected variants and co-infecting agents. However, NGS is not widely used for small RNA viruses because of incorrectly perceived cost estimates and inefficient utilization of freely available bioinformatics tools. Methods In this study, we have utilized NGS-based random sequencing of total RNA combined with barcode multiplexing of libraries to quickly, effectively and simultaneously characterize the genomic sequences of multiple avian paramyxoviruses. Thirty libraries were prepared from diagnostic samples amplified in allantoic fluids and their total RNAs were sequenced in a single flow cell on an Illumina MiSeq instrument. After digital normalization, data were assembled using the MIRA assembler within a customized workflow on the Galaxy platform. Results Twenty-eight avian paramyxovirus 1 (APMV-1), one APMV-13, four avian influenza and two infectious bronchitis virus complete or nearly complete genome sequences were obtained from the single run. The 29 avian paramyxovirus genomes displayed 99.6% mean coverage based on bases with Phred quality scores of 30 or more. The lower and upper quartiles of sample median depth per position for those 29 samples were 2984 and 6894, respectively, indicating coverage across samples sufficient for deep variant analysis. Sample processing and library preparation took approximately 25–30 h, the sequencing run took 39 h, and processing through the Galaxy workflow took approximately 2–3 h. The cost of all steps, excluding labor, was estimated to be 106 USD per sample. Conclusions This work describes an efficient multiplexing NGS approach, a detailed analysis workflow, and customized tools for the characterization of the genomes of RNA viruses. The combination of multiplexing NGS technology with the Galaxy workflow platform resulted in a fast, user-friendly, and cost-efficient protocol for the simultaneous characterization of multiple full-length viral genomes. Twenty-nine full-length or near-full-length APMV genomes with a high median depth were successfully sequenced out of 30 samples. The applied de novo assembly approach also allowed identification of mixed viral populations in some of the samples

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

Presence of vaccine-derived newcastle disease viruses in wild birds

Our study demonstrates the repeated isolation of vaccine-derived Newcastle disease viruses from different species of wild birds across four continents from 1997 through 2014. The data indicate that at least 17 species from ten avian orders occupying different habitats excrete vaccine-derived Newcastle disease viruses. The most frequently reported isolates were detected among individuals in the order Columbiformes (n = 23), followed in frequency by the order Anseriformes (n = 13). Samples were isolated from both free-ranging (n = 47) and wild birds kept in captivity (n = 7). The number of recovered vaccine-derived viruses corresponded with the most widely utilized vaccines, LaSota (n = 28) and Hitchner B1 (n = 19). Other detected vaccine-derived viruses resembled the PHY-LMV2 and V4 vaccines, with five and two cases, respectively. These results and the ubiquitous and synanthropic nature of wild pigeons highlight their potential role as indicator species for the presence of Newcastle disease virus of low virulence in the environment. The reverse spillover of live agents from domestic animals to wildlife as a result of the expansion of livestock industries employing massive amounts of live virus vaccines represent an underappreciated and poorly studied effect of human activity on wildlife119sem informaçã

Presence of Vaccine-Derived Newcastle Disease Viruses in Wild Birds

<div><p>Our study demonstrates the repeated isolation of vaccine-derived Newcastle disease viruses from different species of wild birds across four continents from 1997 through 2014. The data indicate that at least 17 species from ten avian orders occupying different habitats excrete vaccine-derived Newcastle disease viruses. The most frequently reported isolates were detected among individuals in the order <i>Columbiformes</i> (n = 23), followed in frequency by the order <i>Anseriformes</i> (n = 13). Samples were isolated from both free-ranging (n = 47) and wild birds kept in captivity (n = 7). The number of recovered vaccine-derived viruses corresponded with the most widely utilized vaccines, LaSota (n = 28) and Hitchner B1 (n = 19). Other detected vaccine-derived viruses resembled the PHY-LMV2 and V4 vaccines, with five and two cases, respectively. These results and the ubiquitous and synanthropic nature of wild pigeons highlight their potential role as indicator species for the presence of Newcastle disease virus of low virulence in the environment. The reverse spillover of live agents from domestic animals to wildlife as a result of the expansion of livestock industries employing massive amounts of live virus vaccines represent an underappreciated and poorly studied effect of human activity on wildlife.</p></div

Phylogenetic tree of isolates and their relationship to class II NDV viruses.

<p>Phylogenetic analysis based on the complete nucleotide sequence of the fusion gene of isolates representing NDV class II. The evolutionary history was inferred by using the Maximum Likelihood method based on Tamura 3-parameter model with 500 bootstrap replicates [<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0162484#pone.0162484.ref070" target="_blank">70</a>]. The tree with the highest log likelihood (-108983.3717) is shown. A discrete Gamma distribution was used to model evolutionary rate differences among sites (4 categories (+G, parameter = 0.0936). The rate variation model allowed for some sites to be evolutionarily invariable ([+I], 39.7777% sites). The tree is drawn to scale with branch lengths measured in the number of substitutions per site and the percentage of trees in which the associated taxa clustered together are shown below the branches. The analysis involved 81 nucleotide sequences with a total of 1662 positions in the final dataset. Isolates studied in this work are designated in front of the taxa name as follows: USA—●; Ukraine—○; Brazil—□, Bulgaria—■. Evolutionary analyses were conducted in MEGA6 [<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0162484#pone.0162484.ref067" target="_blank">67</a>]. The Roman numerals presented in the taxa names in the phylogenetic trees represent the respective genotype for each isolate, followed by the GenBank identification number, host name (if available), country of isolation, strain designation and country of isolation.</p