What's in a strain? Viral metagenomics identifies genetic variation and contaminating circoviruses in laboratory isolates of pigeon paramyxovirus type 1

Non-Targeted RNA Sequencing: Towards the Development of Universal Clinical Diagnosis Methods for Human and Veterinary Infectious Diseases

Metagenomics offers the potential to replace and simplify classical methods used in the clinical diagnosis of human and veterinary infectious diseases. Metagenomics boasts a high pathogen discovery rate and high specificity, advantages absent in most classical approaches. However, its widespread adoption in clinical settings is still pending, with a slow transition from research to routine use. While longer turnaround times and higher costs were once concerns, these issues are currently being addressed by automation, better chemistries, improved sequencing platforms, better databases, and automated bioinformatics analysis. However, many technical options and steps, each producing highly variable outcomes, have reduced the technology's operational value, discouraging its implementation in diagnostic labs. We present a case for utilizing non-targeted RNA sequencing (NT-RNA-seq) as an ideal metagenomics method for the detection of infectious disease-causing agents in humans and animals. Additionally, to create operational value, we propose to identify best practices for the "core" of steps that are invariably shared among many human and veterinary protocols. Reference materials, sequencing procedures, and bioinformatics standards should accelerate the validation processes necessary for the widespread adoption of this technology. Best practices could be determined through "implementation research" by a consortium of interested institutions working on common samples.

RNA sequencing of avian paramyxovirus (Paramyxoviridae, Avulavirinae) isolates from wild mallards in Belgium, 2021: complete genomes and coinfections

We used untargeted RNA sequencing to characterize three Avulavirinae isolates from pooled samples obtained from wild mallards in Belgium in 2021. The complete genome sequences of two avian Orthoavulavirus-1 (AOAV-1) strains and one avian Paraavulavirus-4 (APMV-4) strain were determined confirming hemagglutination inhibition testing of the virus isolates. In addition, the applied sequencing strategy identified an avian influenza virus (AIV) coinfection in all three virus isolates, confirming weak-positive AIV realtime RT-PCR results from the original sample material. In one AOAV-1 isolate, partial sequences covering all genome segments of an AIV of subtype H11N9 could be de novo assembled from the sequencing data. Besides an AIV coinfection, RNA metagenomic data from the APMV-4 isolate also showed evidence of Alpharetrovirus and Megrivirus coinfection. In total, two AOAV-1 of Class II, genotype I.2 and one APMV-4 complete genome sequences were assembled and compared to publicly available sequences, highlighting the importance of surveillance for poultry pathogens in wild birds. Beyond the insights from full genome characterization of virus isolates, untargeted RNA sequencing strategies provide additional insights in the RNA virome of clinical samples as well as their derived virus isolates that are particularly useful when targeting wild avifauna reservoirs of poultry pathogens.

An Outbreak of Newcastle Disease Virus in the Moscow Region in the Summer of 2022

In August 2022 on a backyard farm in the Moscow region of Russia, mortality was observed among chickens, and all 45 birds of a particular farm died or were slaughtered after the onset of symptoms within a few days. Paramyxovirus was isolated from the diseased birds. Based on the nucleotide sequences of the F and NP gene fragments, it was determined that the virus belonged to subgenotype VII.1 AAvV-1 class II. The cleavage site of the F gene ₁₀₉SGGRRQKRFIG₁₁₉ and T in 546 and 555 position of the NP gene were typical for the velogenic type. The genetically closest NDV isolates were found in Iran. The mean time of death of 10-day-old chicken embryos upon infection with the minimal infectious dose was 52 h, which is typical for the velogenic pathotype. The virus caused 100% death of six-week-old chickens during oral infection as well as 100% mortality of all contact chickens, including those located in remote cages, which proves the ability of the virus to spread not only by the fecal-oral route but also by the aerosol route. That demonstrates a high level of pathogenicity and contagiousness of the isolated strain for chicken. However, mice intranasally infected with high doses of the virus did not die.

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.

Pigeon Circovirus over Three Decades of Research: Bibliometrics, Scoping Review, and Perspectives

The pigeon circovirus (PiCV), first described in the literature in the early 1990s, is considered one of the most important infectious agents affecting pigeon health. Thirty years after its discovery, the current review has employed bibliometric strategies to map the entire accessible PiCV-related research corpus with the aim of understanding its present research landscape, particularly in consideration of its historical context. Subsequently, developments, current knowledge, and important updates were provided. Additionally, this review also provides a textual analysis examining the relationship between PiCV and the young pigeon disease syndrome (YPDS), as described and propagated in the literature. Our examination revealed that usages of the term ‘YPDS’ in the literature are characterizations that are diverse in range, and neither standard nor equivalent. Guided by our understanding of the PiCV research corpus, a conceptualization of PiCV diseases was also presented in this review. Proposed definitions and diagnostic criteria for PiCV subclinical infection (PiCV-SI) and PiCV systemic disease (PiCV-SD) were also provided. Lastly, knowledge gaps and open research questions relevant to future PiCV-related studies were identified and discussed.

Evaluating methods for Avian avulavirus-1 whole genome sequencing

Background

Avian avulavirus-1 (AAvV-1, previously Newcastle Disease Virus) is responsible for poultry and wild birds' disease outbreaks. Numerous whole genome sequencing methods were reported for this virus. These methods included cloning, specific primers amplification, shotgun PCR approaches, Sequence Independent Single Primer Amplification and next generation sequencing platform kits.

Methods

Three methods were used to sequence 173 Israeli Avian avulavirus-1 field isolates and one vaccine strain (VH). The sequencing was performed on Proton and Ion Torrent Personal Genome Machine and to a lesser extent, Illumina MiSeq and NextSeq sequencers. Target specific primers (SP) and Sequence Independent Single Primer Amplification (SISPA) products sequenced via the Ion torrent sequencer had a high error rate and truncated genomes. All the next generation sequencing platform sequencing kits generated high sequence accuracy and near-complete genomic size.

Results

A high level of mutations was observed in the intergenic regions between the avian avulavirus-1 genes. Within genes, multiple regions are more mutated than the Fusion region currently used for typing.

Conclusions

Our findings suggest that the whole genome sequencing by the Ion torrent sequencing kit is sufficient. However, when higher fidelity is desired, the Illumina NextSeq and Proton torrent sequencing kits were found to be preferable.

Phylogenetic, antigenic and biological characterization of pigeon paramyxovirus type 1 circulating in China

Background

For many years, ND has been one of the most important infectious pigeon diseases in China. In recent years, a high mortality has been observed in ND-infected pigeons in China. Mortality is from 40% to 80% or 100% in some cases.

Methods

The full-length genomes of four pigeon paramyxovirus type 1 (PPMV-1) strains, which were isolated from infected pigeons in China in 2012 and 2013, were sequenced and analyzed to determine the phylogenetic characteristics of PPMV-1 circulating in pigeons of China in recent years. Furthermore, cross hemagglutination inhibition and cross virus neutralization assays, as well as animal experiments were conducted to determine the antigenicity and pathogenicity of those viruses. Proteolytic cleavage sites (residues 112–117) of the F proteins were identified as the typical virulence motif, ¹¹²RRQKR↓F¹¹⁷ for all four PPMV-1 strains investigated.

Results

Phylogenetic analysis based on sequences of complete genomes and F gene revealed that the four PPMV-1 isolates and most of recent isolates in China were highly homologous to European isolates from 1998 to 2011. All those isolates were clustered in one clade of genotype VI NDV, termed as subgroup 4bii f. The R value was calculated based on cross hemagglutination inhibition and cross virus neutralization results, and confirmed antigenic difference of the PPMV-1 strains isolated in 2013 from the LaSota vaccine strain. Several mutations were identified in the surface glycoproteins F and HN, which probably gave rise to those antigenic differences.

Conclusion

Our result suggested that the PPMV-1 strain prevailing in China in the last decade diverged from a common ancestor and was presumably transmitted from Europe. PPMV-1 isolates displayed obvious antigenic differences from vaccine strain LaSota. Even though PPMV-1 did not cause high mortality in experimental pigeons, the infected pigeons were exhibiting viral shedding for 3 weeks after infection, suggesting PPMV-1 is a potential threat to NDV control worldwide.

Electronic supplementary material

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

Characterization and Epidemiology of Pigeon Paramyxovirus Type-1 Viruses (PPMV-1) Isolated in Macedonia

We have characterized in this study 10 PPMV-1 isolated from domestic pigeons and one PPMV-1 isolated from a feral pigeon in the period 2007-2012, using both classical methods (HI test and ICPI test) and molecular methods (RT-qPCR, RT-PCR, and nucleotide sequencing). Using phylogenetic analysis of partial fusion gene sequences, these viruses clustered with recent European PPMV-1 isolates (EU/re) within the genotype VIb/1. All isolates possessed virulent cleavage site motifs with variable morbidity and mortality in pigeons. The intracerebral pathogenecity indices of the five isolates ranged from 0.59 to 1.53. The repetitive isolation of PPMV-1 viruses for several consecutive years led toward establishing enzootic presence of the disease in pigeons. A high nucleotide sequence homology between the Macedonian isolates and EU/re isolates was shown. Co-circulation of different isolates in the same holdings was detected. This is the first study to extensively describe the molecular epidemiology of PPMV-1 isolated in Macedonia.

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.

Genomic Characterizations of Six Pigeon Paramyxovirus Type 1 Viruses Isolated from Live Bird Markets in China during 2011 to 2013

The genomes of six pigeon paramyxovirus type 1 (PPMV-1) isolated from symptomless pigeons in live poultry markets during the national active surveillance from 2011 to 2013 were sequenced and analyzed in this study. The complete genome lengths of all isolates were 15,192 nucleotides with the gene order of 3’-NP-P-M-F-HN-L-5’. All isolates had the same motif of ¹¹²RRQKRF¹¹⁷ at the cleavage site of the fusion protein, which was typical of velogenic Newcastle disease virus (NDV). Several mutations were identified in the functional domains of F and HN proteins, including fusion peptide, heptad repeat region, transmembrane domains and neutralizing epitopes. Phylogenetic analysis based on sequences of complete genomes and six genes revealed that all isolates belonged to genotype VI in class II, but at least 2 sub-genotypes were identified. Most isolates were placed into sub-genotype VIb with the exception of pi/GX/1015/13, which was classified in sub-genotype VIa. The obvious antigenic difference between PPMV-1 isolates and La Sota strain was found based on the R-value calculated by cross hemagglutination inhibition (HI) assay. These results provided the evidence that PPMV-1 could be detected from healthy pigeons, and our study may be useful in designing vaccines used in pigeon, and developing molecular diagnostic tools to monitor and prevent future PPMV-1 outbreaks.

A survey of overlooked viral infections in biological experiment systems

It is commonly accepted that there are many unknown viruses on the planet. For the known viruses, do we know their prevalence, even in our experimental systems? Here we report a virus survey using recently published small (s)RNA sequencing datasets. The sRNA reads were assembled and contigs were screened for virus homologues against the NCBI nucleotide (nt) database using the BLASTn program. To our surprise, approximately 30% (28 out of 94) of publications had highly scored viral sequences in their datasets. Among them, only two publications reported virus infections. Though viral vectors were used in some of the publications, virus sequences without any identifiable source appeared in more than 20 publications. By determining the distributions of viral reads and the antiviral RNA interference (RNAi) pathways using the sRNA profiles, we showed evidence that many of the viruses identified were indeed infecting and generated host RNAi responses. As virus infections affect many aspects of host molecular biology and metabolism, the presence and impact of viruses needs to be actively investigated in experimental systems.

Next-generation sequencing in veterinary medicine: how can the massive amount of information arising from high-throughput technologies improve diagnosis, control, and management of infectious diseases?

The development of high-throughput molecular technologies and associated bioinformatics has dramatically changed the capacities of scientists to produce, handle, and analyze large amounts of genomic, transcriptomic, and proteomic data. A clear example of this step-change is represented by the amount of DNA sequence data that can be now produced using next-generation sequencing (NGS) platforms. Similarly, recent improvements in protein and peptide separation efficiencies and highly accurate mass spectrometry have promoted the identification and quantification of proteins in a given sample. These advancements in biotechnology have increasingly been applied to the study of animal infectious diseases and are beginning to revolutionize the way that biological and evolutionary processes can be studied at the molecular level. Studies have demonstrated the value of NGS technologies for molecular characterization, ranging from metagenomic characterization of unknown pathogens or microbial communities to molecular epidemiology and evolution of viral quasispecies. Moreover, high-throughput technologies now allow detailed studies of host-pathogen interactions at the level of their genomes (genomics), transcriptomes (transcriptomics), or proteomes (proteomics). Ultimately, the interaction between pathogen and host biological networks can be questioned by analytically integrating these levels (integrative OMICS and systems biology). The application of high-throughput biotechnology platforms in these fields and their typical low-cost per information content has revolutionized the resolution with which these processes can now be studied. The aim of this chapter is to provide a current and prospective view on the opportunities and challenges associated with the application of massive parallel sequencing technologies to veterinary medicine, with particular focus on applications that have a potential impact on disease control and management.