A retrospective study of Newcastle disease in Kenya.

Newcastle disease (ND) is a major constraint to Kenya's poultry production, which is comprised of approximately 80% indigenous chickens (ICs; caged and free-range system) and 20% exotic chickens (intensive system). This study analyzed cases reported as suspected ND in Kenya between 2005 and 2015. Of the suspected 332 ND reported cases from the three production systems in 27 locations within six Kenyan Agro-Ecological Zones (AEZs), 140 diagnosed as infected with avian orthoavulavirus 1 (AOaV-1; formerly Newcastle disease virus) were present in every year in all AEZs. The numbers of AOaV-1-positive cases differed significantly (p < 0.05) between the production systems across the years depending on the season, climate, and location. In the free-range system, both ambient temperatures and season associated significantly (p = 0.001 and 0.02, respectively) with the number of cases, while in the intensive and caged systems, the positive cases correlated significantly with season and relative humidity, respectively (p = 0.05). Regardless of the production systems, the numbers of clinically sick birds positively correlated with the ambient temperatures (r = 0.6; p < 0.05). Failure to detect AOaV-1 in 58% of the ND cases reported, and mortalities exceeding the observed numbers of clinically sick birds suggest deficiencies in the current ND reporting and diagnostic system. Intensive farmers were the slowest in reporting the cases and diagnostic deficiencies were most evident by failure to test the exposure of ICs to natural infection with AOaV-1 and for the AOaV-1-negative cases lack of testing for other pathogens and/or AOaV-1 variants. This study indicates a need for improved surveillance and diagnostics in Kenyan domestic poultry.

ICTV Virus Taxonomy Profile: Hytrosaviridae.

Hytrosaviridae is a family of large, rod-shaped, enveloped entomopathogenic viruses with dsDNA genomes of 120-190 kbp. Hytrosaviruses (also known as salivary gland hypertrophy viruses) primarily replicate in the salivary glands of adult dipteran flies. Hytrosaviruses infecting the haematophagous tsetse fly and the filth-feeding housefly are assigned to two genera, Glossinavirus and Muscavirus, respectively. Whereas muscavirus infections are only overt, glossinavirus infections can be either covert or overt. Overt infections are characterized by diagnostic salivary gland hypertrophy and cause either partial or complete infertility. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the family Hytrosaviridae, which is available at ictv.global/report/hytrosaviridae.

Hytrosavirus genetic diversity and eco-regional spread in Glossina species.

BACKGROUND: The management of the tsetse species Glossina pallidipes (Diptera; Glossinidae) in Africa by the sterile insect technique (SIT) has been hindered by infections of G. pallidipes production colonies with Glossina pallidipes salivary gland hypertrophy virus (GpSGHV; Hytrosaviridae family). This virus can significantly decrease productivity of the G. pallidipes colonies. Here, we used three highly diverged genes and two variable number tandem repeat regions (VNTRs) of the GpSGHV genome to identify the viral haplotypes in seven Glossina species obtained from 29 African locations and determine their phylogenetic relatedness. RESULTS: GpSGHV was detected in all analysed Glossina species using PCR. The highest GpSGHV prevalence was found in G. pallidipes colonized at FAO/IAEA Insect Pest Control Laboratory (IPCL) that originated from Uganda (100%) and Tanzania (88%), and a lower prevalence in G. morsitans morsitans from Tanzania (58%) and Zimbabwe (20%). Whereas GpSGHV was detected in 25-40% of G. fuscipes fuscipes in eastern Uganda, the virus was not detected in specimens of neighboring western Kenya. Most of the identified 15 haplotypes were restricted to specific Glossina species in distinct locations. Seven haplotypes were found exclusively in G. pallidipes. The reference haplotype H1 (GpSGHV-Uga; Ugandan strain) was the most widely distributed, but was not found in G. swynnertoni GpSGHV. The 15 haplotypes clustered into three distinct phylogenetic clades, the largest contained seven haplotypes, which were detected in six Glossina species. The G. pallidipes-infecting haplotypes H10, H11 and H12 (from Kenya) clustered with H7 (from Ethiopia), which presumably corresponds to the recently sequenced GpSGHV-Eth (Ethiopian) strain. These four haplotypes diverged the most from the reference H1 (GpSGHV-Uga). Haplotypes H1, H5 and H14 formed three main genealogy hubs, potentially representing the ancestors of the 15 haplotypes. CONCLUSION: These data identify G. pallidipes as a significant driver for the generation and diversity of GpSGHV variants. This information may provide control guidance when new tsetse colonies are established and hence, for improved management of the virus in tsetse rearing facilities that maintain multiple Glossina species.

RNA interference-based antiviral immune response against the salivary gland hypertrophy virus in Glossina pallidipes.

BACKGROUND: Glossina pallidipes salivary gland hypertrophy virus (GpSGHV; Hytrosaviridae) is a non-occluded dsDNA virus that specifically infects the adult stages of the hematophagous tsetse flies (Glossina species, Diptera: Glossinidae). GpSGHV infections are usually asymptomatic, but unknown factors can result to a switch to acute symptomatic infection, which is characterized by the salivary gland hypertrophy (SGH) syndrome associated with decreased fecundity that can ultimately lead to a colony collapse. It is uncertain how GpSGHV is maintained amongst Glossina spp. populations but RNA interference (RNAi) machinery, a conserved antiviral defense in insects, is hypothesized to be amongst the host's mechanisms to maintain the GpSGHV in asymptomatic (persistent or latent) infection state. Here, we investigated the involvement of RNAi during GpSGHV infections by comparing the expression of three key RNAi machinery genes, Dicer (DCR), Argonaute (AGO) and Drosha, in artificially virus injected, asymptomatic and symptomatic infected G. pallidipes flies compared to PBS injected (controls) individuals. We further assessed the impact of AGO2 knockdown on virus infection by RT-qPCR quantification of four selected GpSGHV genes, i.e. odv-e66, dnapol, maltodextrin glycosyltransferase (a tegument gene) and SGHV091 (a capsid gene). RESULTS: We show that in response to hemocoelic injections of GpSGHV into G. pallidipes flies, increased virus replication was accompanied by significant upregulation of the expression of three RNAi key genes; AGO1, AGO2 and DCR2, and a moderate increase in the expression of Drosha post injection compared to the PBS-injected controls. Furthermore, compared to asymptomatically infected individuals, symptomatic flies showed significant downregulation of AGO1, AGO2 and Drosha, but a moderate increase in the expression of DCR2. Compared to the controls, knockdown of AGO2 did not have a significant impact on virus infection in the flies as evidenced by unaltered transcript levels of the selected GpSGHV genes. CONCLUSION: The upregulation of the expression of the RNAi genes implicate involvement of this machinery in controlling GpSGHV infections and the establishment of symptomatic GpSGHV infections in Glossina. These findings provide a strategic foundation to understand GpSGHV infections and to control latent (asymptomatic) infections in Glossina spp. and thereby control SGHVs in insect production facilities.

Enhancing vector refractoriness to trypanosome infection: achievements, challenges and perspectives.

With the absence of effective prophylactic vaccines and drugs against African trypanosomosis, control of this group of zoonotic neglected tropical diseases depends the control of the tsetse fly vector. When applied in an area-wide insect pest management approach, the sterile insect technique (SIT) is effective in eliminating single tsetse species from isolated populations. The need to enhance the effectiveness of SIT led to the concept of investigating tsetse-trypanosome interactions by a consortium of researchers in a five-year (2013-2018) Coordinated Research Project (CRP) organized by the Joint Division of FAO/IAEA. The goal of this CRP was to elucidate tsetse-symbiome-pathogen molecular interactions to improve SIT and SIT-compatible interventions for trypanosomoses control by enhancing vector refractoriness. This would allow extension of SIT into areas with potential disease transmission. This paper highlights the CRP's major achievements and discusses the science-based perspectives for successful mitigation or eradication of African trypanosomosis.

Coevolution of hytrosaviruses and host immune responses.

BACKGROUND: Hytrosaviruses (SGHVs; Hytrosaviridae family) are double-stranded DNA (dsDNA) viruses that cause salivary gland hypertrophy (SGH) syndrome in flies. Two structurally and functionally distinct SGHVs are recognized; Glossina pallidipes SGHV (GpSGHV) and Musca domestica SGHV (MdSGHV), that infect the hematophagous tsetse fly and the filth-feeding housefly, respectively. Genome sizes and gene contents of GpSGHV (~ 190 kb; 160-174 genes) and MdSGHV (~ 124 kb; 108 genes) may reflect an evolution with the SGHV-hosts resulting in differences in pathobiology. Whereas GpSGHV can switch from asymptomatic to symptomatic infections in response to certain unknown cues, MdSGHV solely infects symptomatically. Overt SGH characterizes the symptomatic infections of SGHVs, but whereas MdSGHV induces both nuclear and cellular hypertrophy (enlarged non-replicative cells), GpSGHV induces cellular hyperplasia (enlarged replicative cells). Compared to GpSGHV's specificity to Glossina species, MdSGHV infects other sympatric muscids. The MdSGHV-induced total shutdown of oogenesis inhibits its vertical transmission, while the GpSGHV's asymptomatic and symptomatic infections promote vertical and horizontal transmission, respectively. This paper reviews the coevolution of the SGHVs and their hosts (housefly and tsetse fly) based on phylogenetic relatedness of immune gene orthologs/paralogs and compares this with other virus-insect models. RESULTS: Whereas MdSGHV is not vertically transmitted, GpSGHV is both vertically and horizontally transmitted, and the balance between the two transmission modes may significantly influence the pathogenesis of tsetse virus. The presence and absence of bacterial symbionts (Wigglesworthia and Sodalis) in tsetse and Wolbachia in the housefly, respectively, potentially contributes to the development of SGH symptoms. Unlike MdSGHV, GpSGHV contains not only host-derived proteins, but also appears to have evolutionarily recruited cellular genes from ancestral host(s) into its genome, which, although may be nonessential for viral replication, potentially contribute to the evasion of host's immune responses. Whereas MdSGHV has evolved strategies to counteract both the housefly's RNAi and apoptotic responses, the housefly has expanded its repertoire of immune effector, modulator and melanization genes compared to the tsetse fly. CONCLUSIONS: The ecologies and life-histories of the housefly and tsetse fly may significantly influence coevolution of MdSGHV and GpSGHV with their hosts. Although there are still many unanswered questions regarding the pathogenesis of SGHVs, and the extent to which microbiota influence expression of overt SGH symptoms, SGHVs are attractive 'explorers' to elucidate the immune responses of their hosts, and the transmission modes of other large DNA viruses.

Comprehensive annotation of Glossina pallidipes salivary gland hypertrophy virus from Ethiopian tsetse flies: a proteogenomics approach.

Glossina pallidipes salivary gland hypertrophy virus (GpSGHV; family Hytrosaviridae) can establish asymptomatic and symptomatic infection in its tsetse fly host. Here, we present a comprehensive annotation of the genome of an Ethiopian GpSGHV isolate (GpSGHV-Eth) compared with the reference Ugandan GpSGHV isolate (GpSGHV-Uga; GenBank accession number EF568108). GpSGHV-Eth has higher salivary gland hypertrophy syndrome prevalence than GpSGHV-Uga. We show that the GpSGHV-Eth genome has 190 291 nt, a low G+C content (27.9 %) and encodes 174 putative ORFs. Using proteogenomic and transcriptome mapping, 141 and 86 ORFs were mapped by transcripts and peptides, respectively. Furthermore, of the 174 ORFs, 132 had putative transcriptional signals [TATA-like box and poly(A) signals]. Sixty ORFs had both TATA-like box promoter and poly(A) signals, and mapped by both transcripts and peptides, implying that these ORFs encode functional proteins. Of the 60 ORFs, 10 ORFs are homologues to baculovirus and nudivirus core genes, including three per os infectivity factors and four RNA polymerase subunits (LEF4, 5, 8 and 9). Whereas GpSGHV-Eth and GpSGHV-Uga are 98.1 % similar at the nucleotide level, 37 ORFs in the GpSGHV-Eth genome had nucleotide insertions (n = 17) and deletions (n = 20) compared with their homologues in GpSGHV-Uga. Furthermore, compared with the GpSGHV-Uga genome, 11 and 24 GpSGHV ORFs were deleted and novel, respectively. Further, 13 GpSGHV-Eth ORFs were non-canonical; they had either CTG or TTG start codons instead of ATG. Taken together, these data suggest that GpSGHV-Eth and GpSGHV-Uga represent two different lineages of the same virus. Genetic differences combined with host and environmental factors possibly explain the differential GpSGHV pathogenesis observed in different G. pallidipes colonies.

Correlation between structure, protein composition, morphogenesis and cytopathology of Glossina pallidipes salivary gland hypertrophy virus.

The Glossina pallidipes salivary gland hypertrophy virus (GpSGHV) is a dsDNA virus with rod-shaped, enveloped virions. Its 190 kb genome contains 160 putative protein-coding ORFs. Here, the structural components, protein composition and associated aspects of GpSGHV morphogenesis and cytopathology were investigated. Four morphologically distinct structures: the nucleocapsid, tegument, envelope and helical surface projections, were observed in purified GpSGHV virions by electron microscopy. Nucleocapsids were present in virogenic stroma within the nuclei of infected salivary gland cells, whereas enveloped virions were located in the cytoplasm. The cytoplasm of infected cells appeared disordered and the plasma membranes disintegrated. Treatment of virions with 1 % NP-40 efficiently partitioned the virions into envelope and nucleocapsid fractions. The fractions were separated by SDS-PAGE followed by in-gel trypsin digestion and analysis of the tryptic peptides by liquid chromatography coupled to electrospray and tandem mass spectrometry. Using the MaxQuant program with Andromeda as a database search engine, a total of 45 viral proteins were identified. Of these, ten and 15 were associated with the envelope and the nucleocapsid fractions, respectively, whilst 20 were detected in both fractions, most likely representing tegument proteins. In addition, 51 host-derived proteins were identified in the proteome of the virus particle, 13 of which were verified to be incorporated into the mature virion using a proteinase K protection assay. This study provides important information about GpSGHV biology and suggests options for the development of future anti-GpSGHV strategies by interfering with virus-host interactions.

Proteomic footprints of a member of Glossinavirus (Hytrosaviridae): an expeditious approach to virus control strategies in tsetse factories.

The Glossinavirus (Glossina pallidipes salivary gland hypertrophy virus (GpSGHV)) is a rod-shaped enveloped insect virus containing a 190,032 bp-long, circular dsDNA genome. The virus is pathogenic for the tsetse fly Glossina pallidipes and has been associated with the collapse of selected mass-reared colonies. Maintenance of productive fly colonies is critical to tsetse and trypanosomiasis eradication in sub-Saharan Africa using the Sterile Insect Technique. Proteomics, an approach to define the expressed protein complement of a genome, was used to further our understanding of the protein composition, morphology, morphogenesis and pathology of GpSGHV. Additionally, this approach provides potential targets for novel and sustainable molecular-based antiviral strategies to control viral infections in tsetse colonies. To achieve this goal, identification of key protein partners involved in virus transmission is required. In this review, we integrate the available data on GpSGHV proteomics to assess the impact of viral infections on host metabolism and to understand the contributions of such perturbations to viral pathogenesis. The relevance of the proteome findings to tsetse and trypanosomiasis management in sub-Sahara Africa is also considered.

Alternative probe hybridization buffers for target RNA depletion and viral sequence recovery in NGS for poultry samples.

Non-targeted next generation sequencing (NGS) is widely applied to identify the diversity of pathogens in field samples. However, abundance of host RNA (especially rRNA) and other environmental nucleic acids can reduce the abundance of pathogen specific reads of interest, reduce depth of coverage and increase surveillance costs. We presently deplete chicken- and selected bacterial-specific rRNAs in poultry field RNA samples with complementary DNA probes in a commercially available probe hybridization buffer followed by digestion of the RNA:DNA hybrids with RNase H. Because the current buffer is an expensive special order reagent of proprietary composition, we tested in-house and other commercially available buffers and identified a viable alternative that yields equivalent host rRNA depletion and viral-specific reads in poultry samples as the current special order reagent but at a reduced cost.

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.

Virulent Newcastle disease virus genotypes V.3, VII.2, and XIII.1.1 and their coinfections with infectious bronchitis viruses and other avian pathogens in backyard chickens in Tanzania.

Oropharyngeal (OP) and cloacal (CL) swabs from 2049 adult backyard chickens collected at 12 live bird markets, two each in Arusha, Dar es Salaam, Iringa, Mbeya, Morogoro and Tanga regions of Tanzania were screened for Newcastle disease virus (NDV) using reverse transcription real-time PCR (rRT-PCR). The virus was confirmed in 25.23% of the birds (n = 517; rRT-PCR CT ≤ 30), with the highest positivity rates observed in birds from Dar es Salaam region with higher prevalence during the dry season (September-November 2018) compared to the rainy season (January and April-May 2019). Next-generation sequencing of OP/CL samples of 20 out of 32 birds that had high amounts of viral RNAs (CT ≤ 25) resulted in the assembly of 18 complete and two partial genome sequences (15,192 bp and 15,045-15,190 bp in length, respectively) of NDV sub-genotypes V.3, VII.2 and XIII.1.1 (n = 1, 13 and 4 strains, respectively). Two birds had mixed NDV infections (V.3/VII.2 and VII.2/XIII.1.1), and nine were coinfected with viruses of families Astroviridae, Coronaviridae, Orthomyxoviridae, Picornaviridae, Pneumoviridae, and Reoviridae. Of the coinfecting viruses, complete genome sequences of two avastroviruses (a recombinant chicken astrovirus antigenic group-Aii and avian nephritis virus genogroup-5) and two infectious bronchitis viruses (a turkey coronavirus-like recombinant and a GI-19 virus) were determined. The fusion (F) protein F1/F2 cleavage sites of the Tanzanian NDVs have the consensus motifs 112 RRRKR↓F 117 (VII.2 strains) and 112 RRQKR↓F 117 (V.3 and XIII.1.1 strains) consistent with virulent virus; virulence was confirmed by intracerebral pathogenicity index scores of 1.66-1.88 in 1-day-old chicks using nine of the 20 isolates. Phylogenetically, the complete F-gene and full genome sequences regionally cluster the Tanzanian NDVs with, but distinctly from, other strains previously reported in eastern and southern African countries. These data contribute to the understanding of NDV epidemiology in Tanzania and the region.

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 (112G-R-Q-G-R↓F117), 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.

Secuenciación y caracterización del genoma de un aislado similar a VG/GA del ortoavulavirus aviar 1 con un motivo único en el sitio de disociación del gene de fusión. Se identificó una secuencia genómica completa de una cepa similar a la cepa Villegas-Glisson/Universidad de Georgia (VG/GA) del ortoavulavirus aviar 1 (AOAV-1) mediante secuenciación no dirigida de nueva generación de una muestra de hisopo orofaríngeo recolectada de una gallina muerta de traspatio de 12 meses. El aislado tiene un motivo en el sitio de disociación de la proteína de fusión (F) consistente con un ortoavulavirus aviar de baja virulencia, pero tiene un motivo único con fenilalanina en la posición 117 (112G-R-Q-G-R↓F117), que es típico para cepas virulentas del AOAV-1. La diferencia de un nucleótido en el sitio de escisión en comparación con otros virus de baja virulencia hizo que el aislado fuera detectable mediante transcripción reversa y PCR en tiempo real en tiempo real específica del gene F (rtRT-PCR) desarrollada como una prueba de diagnóstico para detectar específicamente a las cepas virulentas. El tiempo medio de muerte determinado en huevos y el índice de patogenicidad intracerebral determinado en pollos clasificaron al aislado como lentogénico. Este es el primer informe en los Estados Unidos de un virus lentogénico similar a VG/GA con un residuo de fenilalanina en la posición 117 del sitio de disociación de la proteína F. Además de la preocupación por el posible cambio patogénico del virus a través de cambios adicionales en el sitio de disociación, nuestro contribuye con un mayor conocimiento por parte del personal de diagnóstico acerca de posibles falsos positivos en las pruebas rtRT-PCR del gene F.

Genome Sequences and Characterization of Chicken Astrovirus and Avian Nephritis Virus from Tanzanian Live Bird Markets.

The enteric chicken astrovirus (CAstV) and avian nephritis virus (ANV) are the type species of the genus Avastrovirus (AAstV; Astroviridae family), capable of causing considerable production losses in poultry. Using next-generation sequencing of a cloacal swab from a backyard chicken in Tanzania, we assembled genome sequences of ANV and CAstV (6918 nt and 7318 nt in length, respectively, excluding poly(A) tails, which have a typical AAstV genome architecture (5'-UTR-ORF1a-ORF1b-ORF2-'3-UTR). They are most similar to strains ck/ANV/BR/RS/6R/15 (82.72%) and ck/CAstV/PL/G059/14 (82.23%), respectively. Phylogenetic and sequence analyses of the genomes and the three open reading frames (ORFs) grouped the Tanzanian ANV and CAstV strains with Eurasian ANV-5 and CAstV-Aii viruses, respectively. Compared to other AAstVs, the Tanzanian strains have numerous amino acid variations (substitutions, insertions and deletions) in the spike region of the capsid protein. Furthermore, CAstV-A has a 4018 nt recombinant fragment in the ORF1a/1b genomic region, predicted to be from Eurasian CAstV-Bi and Bvi parental strains. These data should inform future epidemiological studies and options for AAstV diagnostics and vaccines.

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.

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.

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.

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.

Phylogenetic analysis, molecular changes, and adaptation to chickens of Mexican lineage H5N2 low-pathogenic avian influenza viruses from 1994 to 2019.

The Mexican lineage H5N2 low pathogenic avian influenza viruses (LPAIVs) were first detected in 1994 and mutated to highly pathogenic avian influenza viruses (HPAIVs) in 1994-1995 causing widespread outbreaks in poultry. By using vaccination and other control measures, the HPAIVs were eradicated but the LPAIVs continued circulating in Mexico and spread to several other countries. To get better resolution of the phylogenetics of this virus, the full genome sequences of 44 H5N2 LPAIVs isolated from 1994 to 2011, and 6 detected in 2017 and 2019, were analysed. Phylogenetic incongruence demonstrated genetic reassortment between two separate groups of the Mexican lineage H5N2 viruses between 2005 and 2010. Moreover, the recent H5N2 viruses reassorted with previously unidentified avian influenza viruses. Bayesian phylogeographic results suggested that mechanical transmission involving human activity is the most probable cause of the virus spillover to Central American, Caribbean, and East Asian countries. Increased infectivity and transmission of a 2011 H5N2 LPAIV in chickens compared to a 1994 virus demonstrates improved adaptation to chickens, while low virus shedding, and limited contact transmission was observed in mallards with the same 2011 virus. The sporadic increase in basic amino acids in the HA cleavage site, changes in potential N-glycosylation sites in the HA, and truncations of PB1-F2 should be further examined in relation to the increased infectivity and transmission in poultry. The genetic changes that occur as this lineage of H5N2 LPAIVs continues circulating in poultry is concerning not only because of the effect of these changes on vaccination efficacy, but also because of the potential of the viruses to mutate to the highly pathogenic form. Continued vigilance and surveillance efforts, and the pathogenic and genetic characterization of circulating viruses, are required for the effective control of this virus.

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.