Next generation sequencing of viral RNA genomes

Low-Input, High-Resolution 5' Terminal Filovirus RNA Sequencing with ViBE-Seq

Although next-generation sequencing (NGS) has been instrumental in determining the genomic sequences of emerging RNA viruses, de novo sequence determination often lacks sufficient coverage of the 5' and 3' ends of the viral genomes. Since the genome ends of RNA viruses contain the transcription and genome replication promoters that are essential for viral propagation, a lack of terminal sequence information hinders the efforts to study the replication and transcription mechanisms of emerging and re-emerging viruses. To circumvent this, we have developed a novel method termed ViBE-Seq (Viral Bona Fide End Sequencing) for the high-resolution sequencing of filoviral genome ends using a simple yet robust protocol with high fidelity. This technique allows for sequence determination of the 5' end of viral RNA genomes and mRNAs with as little as 50 ng of total RNA. Using the Ebola virus and Marburg virus as prototypes for highly pathogenic, re-emerging viruses, we show that ViBE-Seq is a reliable technique for rapid and accurate 5' end sequencing of filovirus RNA sourced from virions, infected cells, and tissue obtained from infected animals. We also show that ViBE-Seq can be used to determine whether distinct reverse transcriptases have terminal deoxynucleotidyl transferase activity. Overall, ViBE-Seq will facilitate the access to complete sequences of emerging viruses.

Recent updates on laboratory diagnosis of rabies

Rabies is a lethal viral disease transmitted through the bite of rabid animals. India has a high burden of rabies, contributing to a significant proportion of the global deaths. However, under-reporting of the disease is prevalent due to lack of laboratory confirmation. Laboratory diagnosis of rabies plays a crucial role in differentiating the disease from clinical mimics, initiation of appropriate care, implementing infection control measures and informing disease surveillance. This review provides an overview of the recent advancements in laboratory diagnosis of rabies, aimed at updating physicians involved in diagnosis and management of rabies cases in India.

Phylogenetic Characterization of Rabies Virus Field Isolates Collected from Animals in European Russian Regions in 2009-2022

Rabies is a fatal disease of mammals that poses a high zoonotic risk to humans as well. The distribution of rabies is mainly driven by host animal migration and human-mediated dispersion. To contribute to the global understanding of the rabies virus (RABV) molecular epidemiology, 94 RABV field isolates collected from animals in 13 European Russian regions were phylogenetically characterized using the nearly full-size N gene nucleotide sequences. According to phylogenetic inferences, all isolates belonged to one of the two established phylogenetic groups, either group C (n = 54) or group D (n = 40), which are part of the clade Cosmopolitan of RABVs. Some representatives of group C collected from regions located far apart from each other had a remarkably high level of nucleotide identity. The possibility of the contribution of local bat species to the distribution of RABVs was discussed. Interestingly, over the years, the fraction of group D isolates has been constantly decreasing compared with that of group C isolates. The phylogenetic insights generated herein might have an important contribution to the control and surveillance of animal rabies epidemiology in the region.

Uncovering the endemic circulation of rabies in Cambodia

In epidemiology, endemicity characterizes sustained pathogen circulation in a geographical area, which involves a circulation that is not being maintained by external introductions. Because it could potentially shape the design of public health interventions, there is an interest in fully uncovering the endemic pattern of a disease. Here, we use a phylogeographic approach to investigate the endemic signature of rabies virus (RABV) circulation in Cambodia. Cambodia is located in one of the most affected regions by rabies in the world, but RABV circulation between and within Southeast Asian countries remains understudied. Our analyses are based on a new comprehensive data set of 199 RABV genomes collected between 2014 and 2017 as well as previously published Southeast Asian RABV sequences. We show that most Cambodian sequences belong to a distinct clade that has been circulating almost exclusively in Cambodia. Our results thus point towards rabies circulation in Cambodia that does not rely on external introductions. We further characterize within-Cambodia RABV circulation by estimating lineage dispersal metrics that appear to be similar to other settings, and by performing landscape phylogeographic analyses to investigate environmental factors impacting the dispersal dynamic of viral lineages. The latter analyses do not lead to the identification of environmental variables that would be associated with the heterogeneity of viral lineage dispersal velocities, which calls for a better understanding of local dog ecology and further investigations of the potential drivers of RABV spread in the region. Overall, our study illustrates how phylogeographic investigations can be performed to assess and characterize viral endemicity in a context of relatively limited data.

Large-scale phylogenetic analysis reveals genetic diversity and geographic distribution of rabies virus in South-East and South Asia

South-East Asia (SEA) and South Asia (SA) are two important geographic regions with the most severe enzootic rabies in the world. In these regions, phylogenetic analysis of rabies virus (RABV) has been conducted only at a country level; the results obtained from different countries are scattered and unequal, with a non-uniform system to name RABV genotypes. Therefore, it is difficult to undertake origin-tracking and compare inter-country RABV evolution and transmission. To avoid the confusion in understanding and to generate a panoramic picture of RABV genetic diversity, distribution, and transmission in SEA and SA, the present study conducted a systematic phylogenetic analysis by combining all sequences representing 2368 RABV strains submitted to GenBank by 14 rabies endemic SEA and SA countries. The results showed that RABVs circulating in two regions were classified into four major clades and many subclades: the Asia clade is circulating only in SEA, the Indian subcontinent, and Arctic-like clades only in SA, while the Cosmopolitan clade has been detected in both regions. The results also showed a wide range of hosts were infected by divergent RABV subclades, with dogs being the major transmission source. However, wildlife rabies was also found to be an important issue with 6 wild carnivore species identified as potential sources of spillover risk for sylvatic rabies to humans, domestic animals, and other wild animals. Current findings indicate that the two regions have separate virus clades circulating thus indicating the absence of cross-transmission between the regions. The study emphasizes the importance of phylogenetic analysis in the regions using uniform genotyping and naming systems for rabies surveillance, to coordinate actions of member countries to eliminate dog-mediated human rabies by 2030.

Comparative Neutralization Activity of Commercial Rabies Immunoglobulin against Diverse Lyssaviruses

Novel lyssaviruses, the causative agents of rabies, continue to be described mostly due to increased surveillance in bat hosts. Biologicals for the prevention of rabies in humans have, however, remained largely unchanged for decades. This study aimed to determine if commercial rabies immunoglobulin (RIG) could neutralize diverse lyssaviruses. Two commercial preparations, of human or equine origin, were evaluated against a panel consisting of 13 lyssavirus species. Reduced neutralization was observed for the majority of lyssaviruses compared to rabies virus and was more evident for lyssaviruses outside of phylogroup I. Neutralization of more diverse lyssaviruses only occurred at very high doses, except for Ikoma lyssavirus, which could not be neutralized by the RIG evaluated in this study. The use of RIG is a crucial component of rabies post-exposure prophylaxis and the data generated here indicate that RIG, in its current form, will not protect against all lyssaviruses. In addition, higher doses of RIG may be required for neutralization as the genetic distance from vaccine strains increases. Given the limitations of current RIG preparations, alternative passive immunization options should be investigated.

Whole genome sequencing and phylogenetic characterisation of rabies virus strains from Moldova and north-eastern Romania

BACKGROUND

Rabies is the oldest fatal zoonotic disease recognised as a neglected tropical disease and is caused by an RNA virus belonging to the genus Lyssavirus, family Rhabdoviridae.

METHODOLOGY/PRINCIPAL FINDINGS

A deep molecular analysis was conducted on full-length nucleoprotein (N) gene and whole genome sequences of rabies virus from 37 animal brain samples collected between 2012 and 2017 to study the circulation of rabies virus (RABV) variants. The overall aim was to better understand their distribution in Moldova and north-eastern Romania. Both Sanger and high throughput sequencing on Ion Torrent and Illumina platforms were performed. Phylogenetic analysis of the RABV sequences from both Moldova and Romania revealed that all the samples (irrespective of the year of isolation and the species) belonged to a single phylogenetic group: north-eastern Europe (NEE), clustering into three assigned lineages: RO#5, RO#6 and RO#7.

CONCLUSIONS/SIGNIFICANCE

High throughput sequencing of RABV samples from domestic and wild animals was performed for the first time for both countries, providing new insights into virus evolution and epidemiology in this less studied region, expanding our understanding of the disease.

Impact and mitigation of sampling bias to determine viral spread: Evaluating discrete phylogeography through CTMC modeling and structured coalescent model approximations

Bayesian phylogeographic inference is a powerful tool in molecular epidemiological studies, which enables reconstruction of the origin and subsequent geographic spread of pathogens. Such inference is, however, potentially affected by geographic sampling bias. Here, we investigated the impact of sampling bias on the spatiotemporal reconstruction of viral epidemics using Bayesian discrete phylogeographic models and explored different operational strategies to mitigate this impact. We considered the continuous-time Markov chain (CTMC) model and two structured coalescent approximations (Bayesian structured coalescent approximation [BASTA] and marginal approximation of the structured coalescent [MASCOT]). For each approach, we compared the estimated and simulated spatiotemporal histories in biased and unbiased conditions based on the simulated epidemics of rabies virus (RABV) in dogs in Morocco. While the reconstructed spatiotemporal histories were impacted by sampling bias for the three approaches, BASTA and MASCOT reconstructions were also biased when employing unbiased samples. Increasing the number of analyzed genomes led to more robust estimates at low sampling bias for the CTMC model. Alternative sampling strategies that maximize the spatiotemporal coverage greatly improved the inference at intermediate sampling bias for the CTMC model, and to a lesser extent, for BASTA and MASCOT. In contrast, allowing for time-varying population sizes in MASCOT resulted in robust inference. We further applied these approaches to two empirical datasets: a RABV dataset from the Philippines and a SARS-CoV-2 dataset describing its early spread across the world. In conclusion, sampling biases are ubiquitous in phylogeographic analyses but may be accommodated by increasing the sample size, balancing spatial and temporal composition in the samples, and informing structured coalescent models with reliable case count data.

COVID-19 surveillance in wastewater: An epidemiological tool for the monitoring of SARS-CoV-2

The coronavirus disease 2019 (COVID-19) pandemic has prompted a lot of questions globally regarding the range of information about the virus's possible routes of transmission, diagnostics, and therapeutic tools. Worldwide studies have pointed out the importance of monitoring and early surveillance techniques based on the identification of viral RNA in wastewater. These studies indicated the presence of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA in human feces, which is shed via excreta including mucus, feces, saliva, and sputum. Subsequently, they get dumped into wastewater, and their presence in wastewater provides a possibility of using it as a tool to help prevent and eradicate the virus. Its monitoring is still done in many regions worldwide and serves as an early "warning signal"; however, a lot of limitations of wastewater surveillance have also been identified.

Taiwan Bat Lyssavirus: In Vitro and In Vivo Assessment of the Ability of Rabies Vaccine-Derived Antibodies to Neutralise a Novel Lyssavirus

Rabies is a neglected tropical disease. The prototype virus, the rabies virus, still causes tens of thousands of human fatalities annually. Rabies is one member of the genus Lyssavirus. The burden of other lyssaviruses is unclear. The continued emergence of novel lyssaviruses means that assessment of vaccine efficacy against these viruses is critical, as standard rabies vaccines are not efficacious against all lyssaviruses. Taiwan bat lyssavirus (TWBLV) was first reported in 2018 following isolation from Japanese house bats. Since the initial detection and genetic characterisation, no attempts have been made to antigenically define this virus. Due to the inaccessibility of the wildtype isolate, the successful generation of a live recombinant virus, cSN-TWBLV, is described, where the full-length genome clone of the RABV vaccine strain, SAD-B19, was constructed with the glycoprotein of TWBLV. In vitro and in vivo characterization of cSN-TWBLV was undertaken and demonstrated evidence for cross-neutralisation of cSN-TWBLV with phylogroup I -specific sera and rabies virus standard sera. For neutralisation equivalent to 0.5 IU/mL of WHO and World Organisation of Animal Health (WOAH) sera against CVS, 0.5 IU/mL of WOAH sera and 2.5 IU/mL of WHO sera were required to neutralise cSN-TWBLV. In addition, specific sera for ARAV and EBLV-1 exhibited the highest neutralising antibody titres against cSN-TWBLV, compared to other phylogroup I-specific sera.

Unraveling the viral dark matter through viral metagenomics

Viruses are part of the microbiome and have essential roles in immunology, evolution, biogeochemical cycles, health, and disease progression. Viruses influence a wide variety of systems and processes, and the continued discovery of novel viruses is anticipated to reveal new mechanisms influencing the biology of diverse environments. While the identity and roles of viruses continue to be discovered and understood through viral metagenomics, most of the sequences in virome datasets cannot be attributed to known viruses or may be only distantly related to species already described in public sequence databases, at best. Such viruses are known as the viral dark matter. Ongoing discoveries from the viral dark matter have provided insights into novel viruses from a variety of environments, as well as their potential in immunological processes, virus evolution, health, disease, therapeutics, and surveillance. Increased understanding of the viral dark matter will continue with a combination of cultivation, microscopy, sequencing, and bioinformatic efforts, which are discussed in the present review.

Targeted Virome Sequencing Enhances Unbiased Detection and Genome Assembly of Known and Emerging Viruses-The Example of SARS-CoV-2

Targeted virome enrichment and sequencing (VirCapSeq-VERT) utilizes a pool of oligos (baits) to enrich all known—up to 2015—vertebrate-infecting viruses, increasing their detection sensitivity. The hybridisation of the baits to the target sequences can be partial, thus enabling the detection and genomic reconstruction of novel pathogens with <40% genetic diversity compared to the strains used for the baits’ design. In this study, we deploy this method in multiplexed mixes of viral extracts, and we assess its performance in the unbiased detection of DNA and RNA viruses after cDNA synthesis. We further assess its efficiency in depleting various background genomic material. Finally, as a proof-of-concept, we explore the potential usage of the method for the characterization of unknown, emerging human viruses, such as SARS-CoV-2, which may not be included in the baits’ panel. We mixed positive samples of equimolar DNA/RNA viral extracts from SARS-CoV-2, coronavirus OC43, cytomegalovirus, influenza A virus H3N2, parvovirus B19, respiratory syncytial virus, adenovirus C and coxsackievirus A16. Targeted virome enrichment was performed on a dsDNA mix, followed by sequencing on the NextSeq500 (Illumina) and the portable MinION sequencer, to evaluate its usability as a point-of-care (PoC) application. Genome mapping assembly was performed using viral reference sequences. The untargeted libraries contained less than 1% of total reads mapped on most viral genomes, while RNA viruses remained undetected. In the targeted libraries, the percentage of viral-mapped reads were substantially increased, allowing full genome assembly in most cases. Targeted virome sequencing can enrich a broad range of viruses, potentially enabling the discovery of emerging viruses.

The Comparison of Full G and N Gene Sequences From Turkish Rabies Virus Field Strains

The rabies infection is a zoonotic viral disease in humans and is spread by both wild and domestic carnivores. This study aimed to molecularly characterize the field strains of the rabies virus circulating in Turkey between 2013 and 2020. Brain samples obtained from 16 infected animals (8 cattle, one donkey, three foxes, three dogs, and one marten) were tested. Full nucleoprotein (N) and glycoprotein (G) gene sequences were used to determine the genetic and antigenic characteristics of the rabies virus field strains. The phylogenetic analyses revealed that the 16 field strains identified in Turkey belonged to the Cosmopolitan lineage.

COVID-19: A systematic review and update on prevention, diagnosis, and treatment

Since the rapid onset of the COVID-19 or SARS-CoV-2 pandemic in the world in 2019, extensive studies have been conducted to unveil the behavior and emission pattern of the virus in order to determine the best ways to diagnosis of virus and thereof formulate effective drugs or vaccines to combat the disease. The emergence of novel diagnostic and therapeutic techniques considering the multiplicity of reports from one side and contradictions in assessments from the other side necessitates instantaneous updates on the progress of clinical investigations. There is also growing public anxiety from time to time mutation of COVID-19, as reflected in considerable mortality and transmission, respectively, from delta and Omicron variants. We comprehensively review and summarize different aspects of prevention, diagnosis, and treatment of COVID-19. First, biological characteristics of COVID-19 were explained from diagnosis standpoint. Thereafter, the preclinical animal models of COVID-19 were discussed to frame the symptoms and clinical effects of COVID-19 from patient to patient with treatment strategies and in-silico/computational biology. Finally, the opportunities and challenges of nanoscience/nanotechnology in identification, diagnosis, and treatment of COVID-19 were discussed. This review covers almost all SARS-CoV-2-related topics extensively to deepen the understanding of the latest achievements (last updated on January 11, 2022).

Identification of full-length circular nucleic acids using long-read sequencing technologies

Unlike the traditional perception in genomic DNA or linear RNA, circular nucleic acids are a class of functional biomolecules with a circular configuration and are often observed in nature. These circular molecules encompass the full spectrum of size and play an important role in organisms, making circular nucleic acids research worthy. Due to the low abundance of most types of circular nucleic acids and the disadvantages of short-read sequencing platforms, accurate and full-length circular nucleic acid sequencing and identification is difficult. In this review, we have provided insights into full-length circular nucleic acid detection methods using long-read sequencing technologies, with a focus on the experimental and bioinformatics strategies to obtain accurate sequences.

Potential Virus Involvement in Alzheimer's Disease: Results from a Phase IIa Trial Evaluating Apovir, an Antiviral Drug Combination

Background:

Accumulating data suggest infectious agents are involved in Alzheimer’s disease (AD). The two primary aims of this trial were to assess safety and efficacy of an antiviral drug combination on AD progression.

Objective:

The trial evaluated whether Apovir, a combination of two antiviral agents, pleconaril (active on enteroviruses) and ribavirin (active on several viruses), could slow AD progression.

Methods:

Sixty-nine patients 60–85 years were treated with Apovir or placebo for 9 months and followed until 12 months after end of treatment. Cognitive tests, safety, biomarkers, drug plasma, and cerebrospinal fluid concentrations were assessed.

Results:

The tolerability of Apovir was compromised as demonstrated by the large drop-out rate and increased frequency and severity of adverse events. The primary endpoint, demonstrating a difference in change from baseline to 9 months between groups in ADAS-cog total score, was not met (p = 0.1809). However, there were observations indicating potential effects on both ADAS-cog and CDR-SB but these effects need to be verified. Also, there was a decrease in cerebrospinal fluid amyloid-β in Apovir at 9 months (p = 0.0330) but no change in placebo.

Conclusion:

This was the first randomized, placebo controlled clinical trial exploring antiviral treatment on AD progression. The trial is considered inconclusive due to the large drop-out rate. New trials are needed to verify if the indications of effect observed can be confirmed and which component(s) in Apovir contributed to such effects. Pleconaril alone may be studied to improve the tolerability and to verify if enterovirus is involved in the disease process.

Assessing Rabies Vaccine Protection against a Novel Lyssavirus, Kotalahti Bat Lyssavirus

Rabies is a fatal encephalitis caused by an important group of viruses within the Lyssavirus genus. The prototype virus, rabies virus, is still the most commonly reported lyssavirus and causes approximately 59,000 human fatalities annually. The human and animal burden of the other lyssavirus species is undefined. The original reports for the novel lyssavirus, Kotalahti bat lyssavirus (KBLV), were based on the detection of viral RNA alone. In this report we describe the successful generation of a live recombinant virus, cSN-KBLV; where the full-length genome clone of RABV vaccine strain, SAD-B19, was constructed with the glycoprotein of KBLV. Subsequent in vitro characterisation of cSN-KBLV is described here. In addition, the ability of a human rabies vaccine to confer protective immunity in vivo following challenge with this recombinant virus was assessed. Naïve or vaccinated mice were infected intracerebrally with a dose of 100 focus-forming units/30 µL of cSN-KBLV; all naïve mice and 8% (n = 1/12) of the vaccinated mice succumbed to the challenge, whilst 92% (n = 11/12) of the vaccinated mice survived to the end of the experiment. This report provides strong evidence for cross-neutralisation and cross-protection of cSN-KBLV using purified Vero cell rabies vaccine.

Full-Genome Sequences and Phylogenetic Analysis of Archived Danish European Bat Lyssavirus 1 (EBLV-1) Emphasize a Higher Genetic Resolution and Spatial Segregation for Sublineage 1a

European bat lyssavirus type 1 (EBLV-1) is the causative agent for almost all reported rabies cases found in European bats. In recent years, increasing numbers of available EBLV-1 full genomes and their phylogenetic analyses helped to further elucidate the distribution and genetic characteristics of EBLV-1 and its two subtypes, namely EBLV-1a and EBLV-1b. Nonetheless, the absence of full-genome sequences from regions with known detections of EBLV-1 still limit the understanding of the phylogeographic relations between viruses from different European regions. In this study, a set of 21 archived Danish EBLV-1 samples from the years 1985 to 2009 was processed for the acquisition of full-genome sequences using a high-throughput sequencing approach. Subsequent phylogenetic analysis encompassing all available EBLV-1 full genomes from databases revealed the Danish sequences belong to the EBLV-1a subtype and further highlighted the distinct, close phylogenetic relationship of Danish, Dutch and German isolates in this region. In addition, the formation of five putative groups nearly exclusively formed by Danish isolates and the overall increased resolution of the EBLV-1a branch indicate a higher genetic diversity and spatial segregation for this sublineage than was previously known. These results emphasize the importance of phylogenetic analyses of full-genome sequences of lyssaviruses for genetic geography.

High Throughput Sequencing for the Detection and Characterization of RNA Viruses

This review aims to assess and recommend approaches for targeted and agnostic High Throughput Sequencing of RNA viruses in a variety of sample matrices. HTS also referred to as deep sequencing, next generation sequencing and third generation sequencing; has much to offer to the field of environmental virology as its increased sequencing depth circumvents issues with cloning environmental isolates for Sanger sequencing. That said however, it is important to consider the challenges and biases that method choice can impart to sequencing results. Here, methodology choices from RNA extraction, reverse transcription to library preparation are compared based on their impact on the detection or characterization of RNA viruses.

Classification and specific primer design for accurate detection of SARS-CoV-2 using deep learning

In this paper, deep learning is coupled with explainable artificial intelligence techniques for the discovery of representative genomic sequences in SARS-CoV-2. A convolutional neural network classifier is first trained on 553 sequences from the National Genomics Data Center repository, separating the genome of different virus strains from the Coronavirus family with 98.73% accuracy. The network's behavior is then analyzed, to discover sequences used by the model to identify SARS-CoV-2, ultimately uncovering sequences exclusive to it. The discovered sequences are validated on samples from the National Center for Biotechnology Information and Global Initiative on Sharing All Influenza Data repositories, and are proven to be able to separate SARS-CoV-2 from different virus strains with near-perfect accuracy. Next, one of the sequences is selected to generate a primer set, and tested against other state-of-the-art primer sets, obtaining competitive results. Finally, the primer is synthesized and tested on patient samples (n = 6 previously tested positive), delivering a sensitivity similar to routine diagnostic methods, and 100% specificity. The proposed methodology has a substantial added value over existing methods, as it is able to both automatically identify promising primer sets for a virus from a limited amount of data, and deliver effective results in a minimal amount of time. Considering the possibility of future pandemics, these characteristics are invaluable to promptly create specific detection methods for diagnostics.

Virxicon: A Lexicon Of Viral Sequences

Motivation

Viruses are the most abundant biological entities and constitute a large reservoir of genetic diversity. In recent years, knowledge about them has increased significantly as a result of dynamic development in life sciences and rapid technological progress. This knowledge is scattered across various data repositories, making a comprehensive analysis of viral data difficult.

Results

In response to the need for gathering a comprehensive knowledge of viruses and viral sequences, we developed Virxicon, a lexicon of all experimentally acquired sequences for RNA and DNA viruses. The ability to quickly obtain data for entire viral groups, searching sequences by levels of taxonomic hierarchy—according to the Baltimore classification and ICTV taxonomy—and tracking the distribution of viral data and its growth over time are unique features of our database compared to the other tools.

Availabilityand implementation

Virxicon is a publicly available resource, updated weekly. It has an intuitive web interface and can be freely accessed at http://virxicon.cs.put.poznan.pl/.

Randomly primed, strand-switching, MinION-based sequencing for the detection and characterization of cultured RNA viruses

RNA viruses rapidly mutate, which can result in increased virulence, increased escape from vaccine protection, and false-negative detection results. Targeted detection methods have a limited ability to detect unknown viruses and often provide insufficient data to detect coinfections or identify antigenic variants. Random, deep sequencing is a method that can more fully detect and characterize RNA viruses and is often coupled with molecular techniques or culture methods for viral enrichment. We tested viral culture coupled with third-generation sequencing for the ability to detect and characterize RNA viruses. Cultures of bovine viral diarrhea virus, canine distemper virus (CDV), epizootic hemorrhagic disease virus, infectious bronchitis virus, 2 influenza A viruses, and porcine respiratory and reproductive syndrome virus were sequenced on the MinION platform using a random, reverse primer in a strand-switching reaction, coupled with PCR-based barcoding. Reads were taxonomically classified and used for reference-based sequence building using a stock personal computer. This method accurately detected and identified complete coding sequence genomes with a minimum of 20× coverage depth for all 7 viruses, including a sample containing 2 viruses. Each lineage-typing region had at least 26× coverage depth for all viruses. Furthermore, analyzing the CDV sample through a pipeline devoid of CDV reference sequences modeled the ability of this protocol to detect unknown viruses. Our results show the ability of this technique to detect and characterize dsRNA, negative- and positive-sense ssRNA, and nonsegmented and segmented RNA viruses.

In vitro production of synthetic viral RNAs and their delivery into mammalian cells and the application of viral RNAs in the study of innate interferon responses

Mammalian cells express different types of RNA molecules that can be classified as protein coding RNAs (mRNA) and non-coding RNAs (ncRNAs) the latter of which have housekeeping and regulatory functions in cells. Cellular RNAs are not recognized by cellular pattern recognition receptors (PRRs) and innate immunity is not activated. RNA viruses encode and express RNA molecules that usually differ from cell-specific RNAs and they include for instance 5'capped and 5'mono- and triphosphorylated RNAs, small viral RNAs and viral RNA-protein complexes called vRNPs. These molecules are recognized by certain members of Toll-like receptor (TLR) and RIG-I-like receptor (RLR) families leading to activation of innate immune responses and the production of antiviral cytokines, such as type I and type III interferons (IFNs). Virus-specific ssRNA and dsRNA molecules that mimic the viral genomic RNAs or their replication intermediates can efficiently be produced by bacteriophage T7 DNA-dependent RNA polymerase and bacteriophage phi6 RNA-dependent RNA polymerase, respectively. These molecules can then be delivered into mammalian cells and the mechanisms of activation of innate immune responses can be studied. In addition, synthetic viral dsRNAs can be processed to small interfering RNAs (siRNAs) by a Dicer enzyme to produce a swarm of antiviral siRNAs. Here we describe the biology of RNAs, their in vitro production and delivery into mammalian cells as well as how these molecules can be used to inhibit virus replication and to study the mechanisms of activation of the innate immune system.

Profile of the Spatial Distribution Patterns of the Human and Bacteriophage Virome in a Wastewater Treatment Plant Located in the South of Spain

In wastewater treatment plants, most microbial characterization has focused on bacterial, archaeal, and fungal populations. Due to the difficult isolation, quantification, and identification of viruses, only a limited number of virome studies associated with wastewater treatment plants have been carried out. However, the virus populations play an important role in the microbial dynamics in wastewater treatment systems and the biosafety of effluents. In this work, the viral members present in influent wastewater, mixed liquor (aerobic bioreactor), excess sludge, and effluent water of a conventional activated sludge system for the treatment of urban wastewater were identified. Viral members were observed by transmission electron microscopy and studied through next-generation sequencing studies. The results showed the dominance of bacteriophages in the viral community in all samples, with the dominant viral phylotype classified as Escherichia coli O157 typing phage 7. Moreover, different human viruses, such as Cynomolgus cytomegalovirus and Gammaherpesvirus, were also detected.

Applying next-generation sequencing to unravel the mutational landscape in viral quasispecies

Next-generation sequencing (NGS) has revolutionized the scale and depth of biomedical sciences. Because of its unique ability for the detection of sub-clonal variants within genetically diverse populations, NGS has been successfully applied to analyze and quantify the exceptionally-high diversity within viral quasispecies, and many low-frequency drug- or vaccine-resistant mutations of therapeutic importance have been discovered. Although many works have intensively discussed the latest NGS approaches and applications in general, none of them has focused on applying NGS in viral quasispecies studies, mostly due to the limited ability of current NGS technologies to accurately detect and quantify rare viral variants. Here, we summarize several error-correction strategies that have been developed to enhance the detection accuracy of minority variants. We also discuss critical considerations for preparing a sequencing library from viral RNAs and for analyzing NGS data to unravel the mutational landscape.

Rapid in-country sequencing of whole virus genomes to inform rabies elimination programmes

Genomic surveillance is an important aspect of contemporary disease management but has yet to be used routinely to monitor endemic disease transmission and control in low- and middle-income countries. Rabies is an almost invariably fatal viral disease that causes a large public health and economic burden in Asia and Africa, despite being entirely vaccine preventable. With policy efforts now directed towards achieving a global goal of zero dog-mediated human rabies deaths by 2030, establishing effective surveillance tools is critical. Genomic data can provide important and unique insights into rabies spread and persistence that can direct control efforts. However, capacity for genomic research in low- and middle-income countries is held back by limited laboratory infrastructure, cost, supply chains and other logistical challenges. Here we present and validate an end-to-end workflow to facilitate affordable whole genome sequencing for rabies surveillance utilising nanopore technology. We used this workflow in Kenya, Tanzania and the Philippines to generate rabies virus genomes in two to three days, reducing costs to approximately £60 per genome. This is over half the cost of metagenomic sequencing previously conducted for Tanzanian samples, which involved exporting samples to the UK and a three- to six-month lag time. Ongoing optimization of workflows are likely to reduce these costs further. We also present tools to support routine whole genome sequencing and interpretation for genomic surveillance. Moreover, combined with training workshops to empower scientists in-country, we show that local sequencing capacity can be readily established and sustainable, negating the common misperception that cutting-edge genomic research can only be conducted in high resource laboratories. More generally, we argue that the capacity to harness genomic data is a game-changer for endemic disease surveillance and should precipitate a new wave of researchers from low- and middle-income countries.

Whole-genome next-generation sequencing and phylogenetic characterization of viral haemorrhagic septicaemia virus in Korea

Whole-genome next-generation sequencing was used to investigate the local evolution of viral haemorrhagic septicaemia virus, a serious pathogen affecting economically important fish such as rainbow trout and turbot in Europe and olive flounder in Asia. Sequence analysis showed that all isolates were genotype IVa, but could be classified further into four subgroups (K1-K4). In addition, genomic regions encompassing the nucleoprotein, phosphoprotein, matrix protein and non-virion protein genes, as well as the seven non-coding regions, were relatively conserved, whereas glycoprotein and RNA-dependent RNA polymerase genes were variable in the coding region. Taken together, the data demonstrate that whole-genome next-generation sequencing may be useful for future surveillance, prevention and control strategies against viral haemorrhagic septicaemia.

The population structure of Rose rosette virus in the USA

Rose rosette virus (RRV) (genus Emaravirus) is the causal agent of the homonymous disease, the most destructive malady of roses in the USA. Although the importance of the disease is recognized, little sequence information and no full genomes are available for RRV, a multi-segmented RNA virus. To better understand the population structure of the virus we implemented a Hi-Plex PCR amplicon high-throughput sequencing approach to sequence all 7 segments and to quantify polymorphisms in 91 RRV isolates collected from 16 states in the USA. Analysis revealed insertion/deletion (indel) polymorphisms primarily in the 5' and 3' non-coding, but also within coding regions, including some resulting in changes of protein length. Phylogenetic analysis showed little geographical structuring, suggesting that topography does not have a strong influence on virus evolution. Overall, the virus populations were homogeneous, possibly because of regular movement of plants, the recent emergence of RRV and/or because the virus is under strong purification selection to preserve its integrity and biological functions.

Extension of the known distribution of a novel clade C betacoronavirus in a wildlife host

Disease surveillance in wildlife populations presents a logistical challenge, yet is critical in gaining a deeper understanding of the presence and impact of wildlife pathogens. Erinaceus coronavirus (EriCoV), a clade C Betacoronavirus, was first described in Western European hedgehogs (Erinaceus europaeus) in Germany. Here, our objective was to determine whether EriCoV is present, and if it is associated with disease, in Great Britain (GB). An EriCoV-specific BRYT-Green® real-time reverse transcription PCR assay was used to test 351 samples of faeces or distal large intestinal tract contents collected from casualty or dead hedgehogs from a wide area across GB. Viral RNA was detected in 10.8% (38) samples; however, the virus was not detected in any of the 61 samples tested from Scotland. The full genome sequence of the British EriCoV strain was determined using next generation sequencing; it shared 94% identity with a German EriCoV sequence. Multivariate statistical models using hedgehog case history data, faecal specimen descriptions and post-mortem examination findings found no significant associations indicative of disease associated with EriCoV in hedgehogs. These findings indicate that the Western European hedgehog is a reservoir host of EriCoV in the absence of apparent disease.

Pathogenic viruses: Molecular detection and characterization

Pathogenic viruses are viruses that can infect and replicate within human cells and cause diseases. The continuous emergence and re-emergence of pathogenic viruses has become a major threat to public health. Whenever pathogenic viruses emerge, their rapid detection is critical to enable implementation of specific control measures and the limitation of virus spread. Further molecular characterization to better understand these viruses is required for the development of diagnostic tests and countermeasures. Advances in molecular biology techniques have revolutionized the procedures for detection and characterization of pathogenic viruses. The development of PCR-based techniques together with DNA sequencing technology, have provided highly sensitive and specific methods to determine virus circulation. Pathogenic viruses potentially having global catastrophic consequences may emerge in regions where capacity for their detection and characterization is limited. Development of a local capacity to rapidly identify new viruses is therefore critical. This article reviews the molecular biology of pathogenic viruses and the basic principles of molecular techniques commonly used for their detection and characterization. The principles of good laboratory practices for handling pathogenic viruses are also discussed. This review aims at providing researchers and laboratory personnel with an overview of the molecular biology of pathogenic viruses and the principles of molecular techniques and good laboratory practices commonly implemented for their detection and characterization.

•

The continous emergence and re-emergence of pathogenic viruses has become a major threat to public health.

•

PCR-based techniques together with DNA sequencing technology have provided highly sensitive and specific methods to determine virus circulation.

•

Southeast Asia is considered to be vulnerable to potential outbreaks of pathogenic viruses. A number of pathogenic viruses have been reported to circulate in this region. The 2019 novel coronavirus has also been identified in Southeast Asia.

•

Development of local capacity to rapidly identify new viruses is very important.

Transcriptomic Analysis of Testicular Gene Expression in Normal and Cryptorchid Horses

Simple Summary

Cryptorchidism is a common congenital malformation that results in impaired fertility in horses. The high abdominal temperature and the effects of this disease lead to differences in gene expression between retained testes and descended testes (DTs). Here, we focus on the genetic effects of cryptorchidism. All the differentially expressed genes (DEGs) between undescended testes (UDTs) and DTs were analyzed in this study. A total of 84 DEGs were associated with functions related to sperm development and male reproductive performance. Our study has provided fundamental transcriptomic data for future studies on equine testes and cryptorchidism.

Abstract

Testes produce sperm, and investigations into gene expression in the testes will enhance the understanding of the roles of testicular genes in male reproduction. Cryptorchidism, the failure of one or both testes to descend into the scrotal sac, is a common congenital malformation in horses. The major clinical consequence of this abnormality is impaired fertility. The aim of this study was to analyze the expression patterns of testicular genes and to identify the differentially expressed genes (DEGs) in testes between cryptorchid and normal horses. In this study, the gene expression patterns in equine testes and the DEGs between mature descended testes (DTs) and undescended testes (UDTs) were identified by RNA-seq and validated by real-time qPCR. Our results provide comprehensive transcriptomic data on equine testes. The transcriptomic analysis revealed 11 affected genes that were downregulated in UDTs, possibly as a result of the higher temperature in the abdomen than in the scrotal sac. These 11 genes have previously been associated with male reproduction, and their downregulation might explain the impaired fertility of cryptorchid horses. Two homozygous missense mutations detected in horses with cryptorchidism were absent in normal horses and were listed as potential pathogenic mutations; these mutations should be verified in the future.

Rapid in-country sequencing of whole virus genomes to inform rabies elimination programmes

Genomic surveillance is an important aspect of contemporary disease management but has yet to be used routinely to monitor endemic disease transmission and control in low- and middle-income countries. Rabies is an almost invariably fatal viral disease that causes a large public health and economic burden in Asia and Africa, despite being entirely vaccine preventable. With policy efforts now directed towards achieving a global goal of zero dog-mediated human rabies deaths by 2030, establishing effective surveillance tools is critical. Genomic data can provide important and unique insights into rabies spread and persistence that can direct control efforts. However, capacity for genomic research in low- and middle-income countries is held back by limited laboratory infrastructure, cost, supply chains and other logistical challenges. Here we present and validate an end-to-end workflow to facilitate affordable whole genome sequencing for rabies surveillance utilising nanopore technology. We used this workflow in Kenya, Tanzania and the Philippines to generate rabies virus genomes in two to three days, reducing costs to approximately £60 per genome. This is over half the cost of metagenomic sequencing previously conducted for Tanzanian samples, which involved exporting samples to the UK and a three- to six-month lag time. Ongoing optimization of workflows are likely to reduce these costs further. We also present tools to support routine whole genome sequencing and interpretation for genomic surveillance. Moreover, combined with training workshops to empower scientists in-country, we show that local sequencing capacity can be readily established and sustainable, negating the common misperception that cutting-edge genomic research can only be conducted in high resource laboratories. More generally, we argue that the capacity to harness genomic data is a game-changer for endemic disease surveillance and should precipitate a new wave of researchers from low- and middle-income countries.

Whole genome sequencing and phylogenetic analysis of human metapneumovirus strains from Kenya and Zambia

BACKGROUND

Human metapneumovirus (HMPV) is an important cause of acute respiratory illness in young children. Whole genome sequencing enables better identification of transmission events and outbreaks, which is not always possible with sub-genomic sequences.

RESULTS

We report a 2-reaction amplicon-based next generation sequencing method to determine the complete genome sequences of five HMPV strains, representing three subgroups (A2, B1 and B2), directly from clinical samples. In addition to reporting five novel HMPV genomes from Africa we examined genetic diversity and sequence patterns of publicly available HMPV genomes. We found that the overall nucleotide sequence identity was 71.3 and 80% for HMPV group A and B, respectively, the diversity between HMPV groups was greater at amino acid level for SH and G surface protein genes, and multiple subgroups co-circulated in various countries. Comparison of sequences between HMPV groups revealed variability in G protein length (219 to 241 amino acids) due to changes in the stop codon position. Genome-wide phylogenetic analysis showed congruence with the individual gene sequence sets except for F and M2 genes.

CONCLUSION

This is the first genomic characterization of HMPV genomes from African patients.

Evaluation of Direct Metagenomics and Target Enriched Approaches for High-throughput Sequencing of Field Rabies Viruses

INTRODUCTION

High-throughput sequencing (HTS) identifies random viral fragments in environmental samples metagenomically. High reliability gains it broad application in virus evolution, host-virus interaction, and pathogenicity studies. Deep sequencing of field samples with content of host genetic material and bacteria often produces insufficient data for metagenomics and must be preceded by target enrichment. The main goal of the study was the evaluation of HTS for complete genome sequencing of field-case rabies viruses (RABVs).

MATERIAL AND METHODS

The material was 23 RABVs isolated mainly from red foxes and one European bat lyssavirus-1 isolate propagated in neuroblastoma cells. Three methods of RNA isolation were tested for the direct metagenomics and RABV-enriched approaches. Deep sequencing was performed with a MiSeq sequencer (Illumina) and reagent v3 kit. Bioinformatics data were evaluated by Kraken and Centrifuge software and de novo assembly was done with metaSPAdes.

RESULTS

Testing RNA extraction procedures revealed the deep sequencing scope superiority of the combined TRIzol/column method. This HTS methodology made it possible to obtain complete genomes of all the RABV isolates collected in the field. Significantly greater rates of RABV genome coverages (over 5,900) were obtained with RABV enrichment. Direct metagenomic studies sequenced the full length of 6 out of 16 RABV isolates with a medium coverage between 1 and 71.

CONCLUSION

Direct metagenomics gives the most realistic illustration of the field sample microbiome, but with low coverage. For deep characterisation of viruses, e.g. for spatial and temporal phylogeography during outbreaks, target enrichment is recommended as it covers sequences much more completely.

PAIPline: pathogen identification in metagenomic and clinical next generation sequencing samples

Motivation

Next generation sequencing (NGS) has provided researchers with a powerful tool to characterize metagenomic and clinical samples in research and diagnostic settings. NGS allows an open view into samples useful for pathogen detection in an unbiased fashion and without prior hypothesis about possible causative agents. However, NGS datasets for pathogen detection come with different obstacles, such as a very unfavorable ratio of pathogen to host reads. Alongside often appearing false positives and irrelevant organisms, such as contaminants, tools are often challenged by samples with low pathogen loads and might not report organisms present below a certain threshold. Furthermore, some metagenomic profiling tools are only focused on one particular set of pathogens, for example bacteria.

Results

We present PAIPline, a bioinformatics pipeline specifically designed to address problems associated with detecting pathogens in diagnostic samples. PAIPline particularly focuses on userfriendliness and encapsulates all necessary steps from preprocessing to resolution of ambiguous reads and filtering up to visualization in a single tool. In contrast to existing tools, PAIPline is more specific while maintaining sensitivity. This is shown in a comparative evaluation where PAIPline was benchmarked along other well-known metagenomic profiling tools on previously published well-characterized datasets. Additionally, as part of an international cooperation project, PAIPline was applied to an outbreak sample of hemorrhagic fevers of then unknown etiology. The presented results show that PAIPline can serve as a robust, reliable, user-friendly, adaptable and generalizable stand-alone software for diagnostics from NGS samples and as a stepping stone for further downstream analyses.

Availability and implementation

PAIPline is freely available under https://gitlab.com/rki_bioinformatics/paipline.

Isolation and characterization of a variant duck orthoreovirus causing spleen necrosis in Peking ducks, China

Since December 2017, an infectious disease has caused economic hardship for duck farms and breeding ducks in many regions of China. This disease characterized by spleen necrosis and swelling, is due to a variant strain of duck orthoreovirus (DRV) (Duck/N-DRV-XT18/China/2018), which we isolated from the spleen of diseased ducks. After isolating the virus, we used next-generation sequencing technology to determine the entire genomic of the virus. Our phylogenetic analysis of 10 genomic segments showed that the N-DRV-XT18 strain is closely related to orthoreovirus isolates derived from ducks and geese, with nucleotide sequence identities for 10 genomic fragments ranging between 49.8% and 99.3%. In contract, the nucleotide sequence of N-DRV-XT18 genomic fragments are only 38.6% to 78.8% similar to the chicken orthoreovirus isolate. Therefore, we determined that this pathogen, causing duck spleen necrosis, is a new variant of a duck orthoreovirus that is significantly different from any previously reported waterfowl-derived othoreovirus.

Comparison of intra- and inter-host genetic diversity in rabies virus during experimental cross-species transmission

The development of high-throughput genome sequencing enables accurate measurements of levels of sub-consensus intra-host virus genetic diversity and analysis of the role played by natural selection during cross-species transmission. We analysed the natural and experimental evolution of rabies virus (RABV), an important example of a virus that is able to make multiple host jumps. In particular, we (i) analyzed RABV evolution during experimental host switching with the goal of identifying possible genetic markers of host adaptation, (ii) compared the mutational changes observed during passage with those observed in natura, and (iii) determined whether the colonization of new hosts or tissues requires adaptive evolution in the virus. To address these aims, animal infection models (dog and fox) and primary cell culture models (embryo brain cells of dog and fox) were developed and viral variation was studied in detail through deep genome sequencing. Our analysis revealed a strong unidirectional host evolutionary effect, as dog-adapted rabies virus was able to replicate in fox and fox cells relatively easily, while dogs or neuronal dog cells were not easily susceptible to fox adapted-RABV. This suggests that dog RABV may be able to adapt to some hosts more easily than other host variants, or that when RABV switched from dogs to red foxes it lost its ability to adapt easily to other species. Although no difference in patterns of mutation variation between different host organs was observed, mutations were common following both in vitro and in vivo passage. However, only a small number of these mutations also appeared in natura, suggesting that adaptation during successful cross-species virus transmission is a complex, multifactorial evolutionary process.

Understanding the mechanisms that underpin the cross-species transmission and host adaptation of rabies virus (RABV) remains an important part of the ongoing goal to reduce and eliminate rabies. We utilized next-generation sequencing to perform a deep comparative analysis of the genomic evolution of RABV subpopulations during host adaptation in culture and in animals, with the aim of determining the molecular mechanisms involved in the host-species or tissue adaptation of rabies virus. In particular, we aimed to determine whether experimental evolution can recapitulate evolution in nature. Our results suggest that a limited number of mutations that appeared following both in vitro and in vivo passage were observed in natura. This study also suggests that dog RABV may be able to adapt to some hosts more easily than other host variants.

A novel method for the capture-based purification of whole viral native RNA genomes

Current technologies for targeted characterization and manipulation of viral RNA primarily involve amplification or ultracentrifugation with isopycnic gradients of viral particles to decrease host RNA background. The former strategy is non-compatible for characterizing properties innate to RNA strands such as secondary structure, RNA-RNA interactions, and also for nanopore direct RNA sequencing involving the sequencing of native RNA strands. The latter strategy, ultracentrifugation, causes loss in genomic information due to its inability to retrieve unassembled viral RNA. To address this, we developed a novel application of current nucleic acid hybridization technologies for direct characterization of RNA. In particular, we modified a current enrichment protocol to capture whole viral native RNA genomes for downstream RNA assays to circumvent the abovementioned problems. This technique involves hybridization of biotinylated baits at 500 nucleotides (nt) intervals, stringent washes and release of free native RNA strands using DNase I treatment, with a turnaround time of about 6 h 15 min. RT-qPCR was used as the primary proof of concept that capture-based purification indeed removes host background. Subsequently, capture-based purification was applied to direct RNA sequencing as proof of concept that capture-based purification can be coupled with downstream RNA assays. We report that this protocol was able to successfully purify viral RNA by 561- to 791-fold. We also report that application of this protocol to direct RNA sequencing yielded a reduction in human host RNA background by 1580-fold, a 99.91% recovery of viral genome with at least 15× coverage, and a mean coverage across the genome of 120×. This report is, to the best of our knowledge, the first description of a capture-based purification method for assays that involve direct manipulation or characterisation of native RNA. This report also describes a successful application of capture-based purification as a direct RNA sequencing strategy that addresses certain limitations of current strategies in sequencing RNA viral genomes.

Diversity and Evolution of Viral Pathogen Community in Cave Nectar Bats (Eonycteris spelaea)

Bats are unique mammals, exhibit distinctive life history traits and have unique immunological approaches to suppression of viral diseases upon infection. High-throughput next-generation sequencing has been used in characterizing the virome of different bat species. The cave nectar bat, Eonycteris spelaea, has a broad geographical range across Southeast Asia, India and southern China, however, little is known about their involvement in virus transmission. Here we investigate the diversity and abundance of viral communities from a colony of Eonycteris spelaea residing in Singapore. Our results detected 47 and 22 different virus families from bat fecal and urine samples, respectively. Among these, we identify a large number of virus families including Adenoviridae, Flaviviridae, Reoviridae, Papillomaviridae, Paramyxoviridae, Parvoviridae, Picornaviridae, and Polyomaviridae. In most cases, viral sequences from Eonycteris spelaea are genetically related to a group of bat viruses from other bat genera (e.g., Eidolon, Miniopterus, Rhinolophus and Rousettus). The results of this study improve our knowledge of the host range, spread and evolution of several important viral pathogens. More significantly, our findings provide a baseline to study the temporal patterns of virus shedding and how they correlate with bat phenological trends.

Identification and genetic characterization of a novel Orthobunyavirus species by a straightforward high-throughput sequencing-based approach

Identification and characterization of novel unknown viruses is of great importance. The introduction of high-throughput sequencing (HTS)-based methods has paved the way for genomics-based detection of pathogens without any prior assumptions about the characteristics of the organisms. However, the use of HTS for the characterization of viral pathogens from clinical samples remains limited. Here, we report the identification of a novel Orthobunyavirus species isolated from horse plasma. The identification was based on a straightforward HTS approach. Following enrichment in cell culture, RNA was extracted from the growth medium and rapid library preparation, HTS and primary bioinformatic analyses were performed in less than 12 hours. Taxonomical profiling of the sequencing reads did not reveal sequence similarities to any known virus. Subsequent application of de novo assembly tools to the sequencing reads produced contigs, of which three showed some similarity to the L, M, and S segments of viruses belonging to the Orthobunyavirus genus. Further refinement of these contigs resulted in high-quality, full-length genomic sequences of the three genomic segments (L, M and S) of a novel Orthobunyavirus. Characterization of the genomic sequence, including the prediction of open reading frames and the inspection of consensus genomic termini and phylogenetic analysis, further confirmed that the novel virus is indeed a new species, which we named Ness Ziona virus.

Rapid Sequencing of Multiple RNA Viruses in Their Native Form

Long-read nanopore sequencing by a MinION device offers the unique possibility to directly sequence native RNA. We combined an enzymatic poly-A tailing reaction with the native RNA sequencing to (i) sequence complex population of single-stranded (ss)RNA viruses in parallel, (ii) detect genome, subgenomic mRNA/mRNA simultaneously, (iii) detect a complex transcriptomic architecture without the need for assembly, (iv) enable real-time detection. Using this protocol, positive-ssRNA, negative-ssRNA, with/without a poly(A)-tail, segmented/non-segmented genomes were mixed and sequenced in parallel. Mapping of the generated sequences on the reference genomes showed 100% length recovery with up to 97% identity. This work provides a proof of principle and the validity of this strategy, opening up a wide range of applications to study RNA viruses.

Bats and Viruses: Emergence of Novel Lyssaviruses and Association of Bats with Viral Zoonoses in the EU

Bats in the EU have been associated with several zoonotic viral pathogens of significance to both human and animal health. Virus discovery continues to expand the existing understating of virus classification, and the increased interest in bats globally as reservoirs or carriers of zoonotic agents has fuelled the continued detection and characterisation of new lyssaviruses and other viral zoonoses. Although the transmission of lyssaviruses from bat species to humans or terrestrial species appears rare, interest in these viruses remains, through their ability to cause the invariably fatal encephalitis—rabies. The association of bats with other viral zoonoses is also of great interest. Much of the EU is free of terrestrial rabies, but several bat species harbor lyssaviruses that remain a risk to human and animal health. Whilst the rabies virus is the main cause of rabies globally, novel related viruses continue to be discovered, predominantly in bat populations, that are of interest purely through their classification within the lyssavirus genus alongside the rabies virus. Although the rabies virus is principally transmitted from the bite of infected dogs, these related lyssaviruses are primarily transmitted to humans and terrestrial carnivores by bats. Even though reports of zoonotic viruses from bats within the EU are rare, to protect human and animal health, it is important characterise novel bat viruses for several reasons, namely: (i) to investigate the mechanisms for the maintenance, potential routes of transmission, and resulting clinical signs, if any, in their natural hosts; (ii) to investigate the ability of existing vaccines, where available, to protect against these viruses; (iii) to evaluate the potential for spill over and onward transmission of viral pathogens in novel terrestrial hosts. This review is an update on the current situation regarding zoonotic virus discovery within bats in the EU, and provides details of potential future mechanisms to control the threat from these deadly pathogens.

Viral Enrichment Methods Affect the Detection but Not Sequence Variation of West Nile Virus in Equine Brain Tissue

West Nile virus (WNV), a small, positive sense, single stranded RNA virus continues to encroach into new locales with emergence of new viral variants. Neurological disease in the equine can be moderate to severe in the face of low to undetectable virus loads. Physical methods of virus enrichment may increase sensitivity of virus detection and enhance analysis of viral diversity, especially for deep sequencing studies. However, the use of these techniques is limited mainly to non-neural tissues. We investigated the hypothesis that elimination of equine brain RNA enhances viral detection without limiting viral variation. Eight different WNV viral RNA enrichment and host RNA separation methods were evaluated to determine if elimination of host RNA enhanced detection of WNV and increase the repertoire of virus variants for sequencing. Archived brain tissue from 21 different horses was inoculated with WNV, homogenized, before enrichment and separation. The protocols utilized combinations of low-speed centrifugation, syringe filtration, and nuclease treatment. Viral and host RNA were analyzed using real-time PCR targeting the WNV Envelope (E) protein and equine G3PDH to determine relative sensitivity for WNV and host depletion, respectively. To determine the effect of these methods on viral variation, deep sequencing of the E protein was performed. Our results demonstrate that additional separation and enrichment methods resulted in loss of virus in the face of host RNA depletion. DNA sequencing showed no significant difference in total sequence variation between the RNA enrichment protocols. For equine brain infected with WNV, direct RNA extraction followed by host RNA depletion was most suitable. This study highlights the importance of evaluating viral enrichment and separation methods according to tissue type before embarking on studies where quantification of virus and viral variants is essential to the outcome of the study.

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

BACKGROUND

Newcastle disease (ND) outbreaks are global challenges to the poultry industry. Effective management requires rapid identification and virulence prediction of the circulating Newcastle disease viruses (NDV), the causative agent of ND. However, these diagnostics are hindered by the genetic diversity and rapid evolution of NDVs.

METHODS

An amplicon sequencing (AmpSeq) workflow for virulence and genotype prediction of NDV samples using a third-generation, real-time DNA sequencing platform is described here. 1D MinION sequencing of barcoded NDV amplicons was performed using 33 egg-grown isolates, (15 NDV genotypes), and 15 clinical swab samples collected from field outbreaks. Assembly-based data analysis was performed in a customized, Galaxy-based AmpSeq workflow. MinION-based results were compared to previously published sequences and to sequences obtained using a previously published Illumina MiSeq workflow.

RESULTS

For all egg-grown isolates, NDV was detected and virulence and genotype were accurately predicted. For clinical samples, NDV was detected in ten of eleven NDV samples. Six of the clinical samples contained two mixed genotypes as determined by MiSeq, of which the MinION method detected both genotypes in four samples. Additionally, testing a dilution series of one NDV isolate resulted in NDV detection in a dilution as low as 101 50% egg infectious dose per milliliter. This was accomplished in as little as 7 min of sequencing time, with a 98.37% sequence identity compared to the expected consensus obtained by MiSeq.

CONCLUSION

The depth of sequencing, fast sequencing capabilities, accuracy of the consensus sequences, and the low cost of multiplexing allowed for effective virulence prediction and genotype identification of NDVs currently circulating worldwide. The sensitivity of this protocol was preliminary tested using only one genotype. After more extensive evaluation of the sensitivity and specificity, this protocol will likely be applicable to the detection and characterization of NDV.

An optimized methodology for whole genome sequencing of RNA respiratory viruses from nasopharyngeal aspirates

Over the last decade, the number of viral genome sequences deposited in available databases has grown exponentially. However, sequencing methodology vary widely and many published works have relied on viral enrichment by viral culture or nucleic acid amplification with specific primers rather than through unbiased techniques such as metagenomics. The genome of RNA viruses is highly variable and these enrichment methodologies may be difficult to achieve or may bias the results. In order to obtain genomic sequences of human respiratory syncytial virus (HRSV) from positive nasopharyngeal aspirates diverse methodologies were evaluated and compared. A total of 29 nearly complete and complete viral genomes were obtained. The best performance was achieved with a DNase I treatment to the RNA directly extracted from the nasopharyngeal aspirate (NPA), sequence-independent single-primer amplification (SISPA) and library preparation performed with Nextera XT DNA Library Prep Kit with manual normalization. An average of 633,789 and 1,674,845 filtered reads per library were obtained with MiSeq and NextSeq 500 platforms, respectively. The higher output of NextSeq 500 was accompanied by the increasing of duplicated reads percentage generated during SISPA (from an average of 1.5% duplicated viral reads in MiSeq to an average of 74% in NextSeq 500). HRSV genome recovery was not affected by the presence or absence of duplicated reads but the computational demand during the analysis was increased. Considering that only samples with viral load ≥ E+06 copies/ml NPA were tested, no correlation between sample viral loads and number of total filtered reads was observed, nor with the mapped viral reads. The HRSV genomes showed a mean coverage of 98.46% with the best methodology. In addition, genomes of human metapneumovirus (HMPV), human rhinovirus (HRV) and human parainfluenza virus types 1–3 (HPIV1-3) were also obtained with the selected optimal methodology.

Defining objective clusters for rabies virus sequences using affinity propagation clustering

Rabies is caused by lyssaviruses, and is one of the oldest known zoonoses. In recent years, more than 21,000 nucleotide sequences of rabies viruses (RABV), from the prototype species rabies lyssavirus, have been deposited in public databases. Subsequent phylogenetic analyses in combination with metadata suggest geographic distributions of RABV. However, these analyses somewhat experience technical difficulties in defining verifiable criteria for cluster allocations in phylogenetic trees inviting for a more rational approach. Therefore, we applied a relatively new mathematical clustering algorythm named ‘affinity propagation clustering’ (AP) to propose a standardized sub-species classification utilizing full-genome RABV sequences. Because AP has the advantage that it is computationally fast and works for any meaningful measure of similarity between data samples, it has previously been applied successfully in bioinformatics, for analysis of microarray and gene expression data, however, cluster analysis of sequences is still in its infancy. Existing (516) and original (46) full genome RABV sequences were used to demonstrate the application of AP for RABV clustering. On a global scale, AP proposed four clusters, i.e. New World cluster, Arctic/Arctic-like, Cosmopolitan, and Asian as previously assigned by phylogenetic studies. By combining AP with established phylogenetic analyses, it is possible to resolve phylogenetic relationships between verifiably determined clusters and sequences. This workflow will be useful in confirming cluster distributions in a uniform transparent manner, not only for RABV, but also for other comparative sequence analyses.

Rabies is one of the oldest known zoonoses, caused by lyssaviruses. In recent years, more than 21,000 nucleotide sequences for rabies viruses (RABV) have been deposited in public databases. In this study, a novel mathematical approach called affinity propagation (AP) clustering, a highly powerful tool, to verifiably divide full genome RABV sequences into genetic clusters, was used. A panel of existing and novel RABV full genome sequences was used to demonstrate the application of AP for RABV clustering. Using a combination of AP with established phylogenetic analyses is useful in resolving phylogenetic relationships between more objectively determined clusters and sequences. This workflow will help to substantiate a transparent cluster distribution, not only for RABV, but also for other comparative sequence analyses.

Lyssavirus in Indian Flying Foxes, Sri Lanka

A novel lyssavirus was isolated from brains of Indian flying foxes (Pteropus medius) in Sri Lanka. Phylogenetic analysis of complete virus genome sequences, and geographic location and host species, provides strong evidence that this virus is a putative new lyssavirus species, designated as Gannoruwa bat lyssavirus.