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

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

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.

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.

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.

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.

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.

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.

Surveillance and Assessment of Risk Factors for Newcastle Disease Virus from Live Bird Retail Stalls in Lahore District of Pakistan

Live bird markets (LBMs) in Asian countries are considered hubs for the spread of several poultry viruses. In Pakistan, there is a lack of uniformity in practices used in LBMs, which leads to the spread of poultry diseases. A cross-sectional survey was conducted in June-October 2017 to determine the circulation of Newcastle disease virus (NDV) in chickens being sold in live bird retail stalls (LBRSs) and to identify potential risk factors associated with estimated prevalence. A total of 189 stalls (n = 1134 birds) distributed in eight administrative towns of Lahore were visited. A pool of six oropharyngeal swabs was collected from each stall and tested by real-time reverse transcriptase PCR for the presence of NDV. Forty-two out of 189 swabs were found positive with an overall prevalence of 22.22% (95% confidence interval [Cl]: 16.88%-28.67%). Data for 11 potential risk factors acquired through questionnaires were analyzed by survey-weighted logistic regression and prevalence odds ratios (ORs) for associated risk factors were calculated. A final multivariable model identified three risk factors for NDV prevalence in LBRSs, including trading other poultry breeds alongside broilers (OR = 2.41; 95% confidence interval [CI] = 1.5-6.1), purchasing birds from mixed sources (OR = 3.12; 95% CI = 1.4-11.9), and number of birds sold per day (OR = 6.32; 95% CI = 1.9-23.5). Additionally, 24 selected samples were sequenced and phylogenetic analysis of the complete fusion gene (1662 bp) revealed that all isolates belonged to Subgenotype VII.2. This study provides important information on the epidemiology of NDV in Pakistan and highlights the importance of implementing surveillance and biosecurity practices in LBRSs.

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 10^4.5 or 10^3.5 50% embryo infectious dose (EID₅₀/egg) and chickens were challenged at 21 days post-hatch (DPH). In the second experiment, 10^3.5 EID₅₀/egg of each vaccine was administered at 19 DOE, and chickens were challenged at 14 DPH. Chickens vaccinated with 10^3.5 EID₅₀/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.

Zoonotic and Reverse Zoonotic Transmissibility of SARS-CoV-2

The Coronavirus Disease 2019 (COVID-19) is the first known pandemic caused by a coronavirus. Its causative agent, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), appears to be capable of infecting different mammalian species. Recent detections of this virus in pet, zoo, wild, and farm animals have compelled inquiry regarding the zoonotic (animal-to-human) and reverse zoonotic (human-to-animal) transmissibility of SARS-CoV-2 with the potential of COVID-19 pandemic evolving into a panzootic.

It is important to monitor the global spread of disease and to assess the significance of genomic changes to support prevention and control efforts during a pandemic. An understanding of the SARS-CoV-2 epidemiology provides opportunities to prevent the risk of repeated re-infection of humans and requires a robust One Health-based investigation. This review paper describes the known properties and the existing gaps in scientific knowledge about the zoonotic and reverse zoonotic transmissibility of the novel virus SARS-CoV-2 and the COVID-19 disease it causes.

Experimental Infection and Transmission of Newcastle Disease Vaccine Virus in Four Wild Passerines

Our prior work has shown that live poultry vaccines have been intermittently isolated from wild birds sampled during field surveillance studies for Newcastle disease virus (NDV). Thus, we experimentally investigated the susceptibility of four native agriculturally associated wild bird species to the NDV LaSota vaccine and evaluated the shedding dynamics, potential transmission from chickens, and humoral antibody responses. To test susceptibility, we inoculated wild-caught, immunologically NDV-naïve house finches (Haemorhous mexicanus; n = 16), brown-headed cowbirds (Molothrus ater; n = 9), northern cardinals (Cardinalis cardinalis; n = 6), and American goldfinches (Spinus tristis; n = 12) with 0.1 ml (10^6.7 mean embryo infectious doses [EID_50/ml]) of NDV LaSota vaccine via the oculo-nasal route. To test transmission between chickens and wild birds, adult specific-pathogen-free white leghorn chickens were inoculated similarly and cohoused in separate isolators with two to five wild birds of the species listed above. This design resulted in three treatments: wild bird direct inoculation (five groups) and wild bird exposure to one (two groups) or two inoculated chickens (six groups), respectively. Blood and oropharyngeal and cloacal swabs were collected before and after infection with the live vaccine. All wild birds that were directly inoculated with the LaSota vaccine shed virus as demonstrated by virus isolation (VI). Cardinals were the most susceptible species based on shedding viruses from 1 to 11 days postinoculation (dpi) with titers up to 10^4.9 EID₅₀/ml. Although LaSota viruses were shed by all inoculated chickens and were present in the drinking water, most noninoculated wild birds cohoused with these chickens remained uninfected for 14 days as evidenced by VI. However, one American goldfinch tested positive for vaccine transmission by VI at 7 dpi and one house finch tested positive for vaccine transmission by real-time reverse-transcription PCR at 13 dpi. Only one directly inoculated cowbird (out of three) and two cardinals (out of two) developed NDV-specific hemagglutination inhibition antibody titers of 16, 16, and 128, respectively. No clinical signs were detected in the chickens or the wild birds postinoculation.

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.

Presence of Newcastle disease viruses of sub-genotypes Vc and VIn in backyard chickens and in apparently healthy wild birds from Mexico in 2017

Virulent Newcastle disease viruses (NDV) have been present in Mexico since 1946, and recently, multiple outbreaks have been reported in the country. Here, we characterized eleven NDV isolated from apparently healthy wild birds and backyard chickens in three different locations of Jalisco, Mexico in 2017. Total RNA from NDV was reverse-transcribed, and 1285 nucleotides, which includes 3/4 of the fusion gene, was amplified and sequenced using a long-read MinION sequencing method. The sequences were 99.99-100% identical to the corresponding region obtained using the Illumina MiSeq. Phylogenetic analysis using MinION sequences demonstrated that nine virulent NDV from wild birds belonged to sub-genotypes Vc and VIn, and two backyard chicken isolates were of sub-genotype Vc. The sub-genotype Vc viruses had nucleotide sequence identity that ranged from 97.7 to 98% to a virus of the same sub-genotype isolated from a chicken in Mexico in 2010. Three viruses from pigeons had 96.3-98.7% nucleotide identity to sub-genotype VIn pigeon viruses, commonly referred to as pigeon paramyxovirus, isolated in the USA during 2000-2016. This study demonstrates that viruses of sub-genotype Vc are still present in Mexico, and the detection of this sub-genotype in both chickens and wild birds suggests that transmission among these species may represent a biosecurity risk. This is the first detection and complete genome sequencing of genotype VI NDV from Mexico. In addition, the utilization of an optimized long-read sequencing method for rapid virulence and genotype identification using the Oxford nanopore MinION system is demonstrated.

Whole-genome sequencing of genotype VI Newcastle disease viruses from formalin-fixed paraffin-embedded tissues from wild pigeons reveals continuous evolution and previously unrecognized genetic diversity in the U.S

Background

Newcastle disease viruses (NDV) are highly contagious and cause disease in both wild birds and poultry. A pigeon-adapted variant of genotype VI NDV, often termed pigeon paramyxovirus 1, is commonly isolated from columbids in the United States and worldwide. Complete genomic characterization of these genotype VI viruses circulating in wild columbids in the United States is limited, and due to the genetic variability of the virus, failure of rapid diagnostic detection has been reported. Therefore, in this study, formalin-fixed paraffin-embedded (FFPE) samples were subjected to next-generation sequencing (NGS) to identify and characterize these circulating viruses, providing valuable genetic information. NGS enables multiple samples to be deep-sequenced in parallel. When used on FFPE samples, this methodology allows for retrospective studies of infectious organisms.

Methods

FFPE wild pigeon tissue samples (kidney, liver and spleen) from 10 mortality events in the U.S. between 2010 and 2016 were analyzed using NGS to detect and sequence NDV genomes from randomly amplified total RNA. Results were compared to the previously published immunohistochemistry (IHC) results conducted on the same samples. Additionally, phylogenetic analyses were conducted on the complete and partial fusion gene and complete genome coding sequences.

Results

Twenty-three out of 29 IHC-positive FFPE pigeon samples were identified as positive for NDV by NGS. Positive samples produced an average genome coverage of 99.6% and an average median depth of 199. A previously described sub-genotype (VIa) and a novel sub-genotype (VIn) of NDV were identified as the causative agent of 10 pigeon mortality events in the U.S. from 2010 to 2016. The distribution of these viruses from the North American lineages match the distribution of the Eurasian collared-doves and rock pigeons in the U.S.

Conclusions

This work reports the first successful evolutionary study using deep sequencing of complete NDV genomes from FFPE samples of wild bird origin. There are at least two distinct U.S. lineages of genotype VI NDV maintained in wild pigeons that are continuously evolving independently from each other and have no evident epidemiological connections to viruses circulating abroad. These findings support the hypothesis that columbids are serving as reservoirs of virulent NDV in the U.S.

Electronic supplementary material

The online version of this article (10.1186/s12985-017-0914-2) contains supplementary material, which is available to authorized users.

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

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.

Electronic supplementary material

The online version of this article (10.1186/s12917-017-1211-4) contains supplementary material, which is available to authorized users.

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

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12985-017-0741-5) contains supplementary material, which is available to authorized users.

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