Evaluation of convenient pretreatment protocols for RNA virus metagenomics in serum and tissue samples

An interlaboratory proficiency test using metagenomic sequencing as a diagnostic tool for the detection of RNA viruses in swine fecal material

Metagenomic shotgun sequencing (mNGS) can serve as a generic molecular diagnostic tool. An mNGS proficiency test (PT) was performed in six European veterinary and public health laboratories to detect porcine astroviruses in fecal material and the extracted RNA. While different mNGS workflows for the generation of mNGS data were used in the different laboratories, the bioinformatic analysis was standardized using a metagenomic read classifier as well as read mapping to selected astroviral reference genomes to assess the semiquantitative representation of astrovirus species mixtures. All participants successfully identified and classified most of the viral reads to the two dominant species. The normalized read counts obtained by aligning reads to astrovirus reference genomes by Bowtie2 were in line with Kraken read classification counts. Moreover, participants performed well in terms of repeatability when the fecal sample was tested in duplicate. However, the normalized read counts per detected astrovirus species differed substantially between participants, which was related to the different laboratory methods used for data generation. Further modeling of the mNGS data indicated the importance of selecting appropriate reference data for mNGS read classification. As virus- or sample-specific biases may apply, caution is needed when extrapolating this swine feces-based PT for the detection of other RNA viruses or using different sample types. The suitability of experimental design to a given pathogen/sample matrix combination, quality assurance, interpretation, and follow-up investigation remain critical factors for the diagnostic interpretation of mNGS results.

IMPORTANCE

Metagenomic shotgun sequencing (mNGS) is a generic molecular diagnostic method, involving laboratory preparation of samples, sequencing, bioinformatic analysis of millions of short sequences, and interpretation of the results. In this paper, we investigated the performance of mNGS on the detection of porcine astroviruses, a model for RNA viruses in a pig fecal material, among six European veterinary and public health laboratories. We showed that different methods for data generation affect mNGS performance among participants and that the selection of reference genomes is crucial for read classification. Follow-up investigation remains a critical factor for the diagnostic interpretation of mNGS results. The paper contributes to potential improvements of mNGS as a diagnostic tool in clinical settings.

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.

Tips and tools to obtain and assess mosquito viromes

Due to their vectorial capacity, mosquitoes (Diptera: Culicidae) receive special attention from health authorities and entomologists. These cosmopolitan insects are responsible for the transmission of many viral diseases, such as dengue and yellow fever, causing huge impacts on human health and justifying the intensification of research focused on mosquito-borne diseases. In this context, the study of the virome of mosquitoes can contribute to anticipate the emergence and/or the reemergence of infectious diseases. The assessment of mosquito viromes also contributes to the surveillance of a wide variety of viruses found in these insects, allowing the early detection of pathogens with public health importance. However, the study of mosquito viromes can be challenging due to the number and complexities of steps involved in this type of research. Therefore, this article aims to describe, in a straightforward and simplified way, the steps necessary for obtention and assessment of mosquito viromes. In brief, this article explores: the capture and preservation of specimens; sampling strategies; treatment of samples before DNA/RNA extraction; extraction methodologies; enrichment and purification processes; sequencing choices; and bioinformatics analysis.

Next-Generation Sequencing for the Detection of Microbial Agents in Avian Clinical Samples

Direct-targeted next-generation sequencing (tNGS), with its undoubtedly superior diagnostic capacity over real-time PCR (RT-PCR), and direct-non-targeted NGS (ntNGS), with its higher capacity to identify and characterize multiple agents, are both likely to become diagnostic methods of choice in the future. tNGS is a rapid and sensitive method for precise characterization of suspected agents. ntNGS, also known as agnostic diagnosis, does not require a hypothesis and has been used to identify unsuspected infections in clinical samples. Implemented in the form of multiplexed total DNA metagenomics or as total RNA sequencing, the approach produces comprehensive and actionable reports that allow semi-quantitative identification of most of the agents present in respiratory, cloacal, and tissue samples. The diagnostic benefits of the use of direct tNGS and ntNGS are high specificity, compatibility with different types of clinical samples (fresh, frozen, FTA cards, and paraffin-embedded), production of nearly complete infection profiles (viruses, bacteria, fungus, and parasites), production of "semi-quantitative" information, direct agent genotyping, and infectious agent mutational information. The achievements of NGS in terms of diagnosing poultry problems are described here, along with future applications. Multiplexing, development of standard operating procedures, robotics, sequencing kits, automated bioinformatics, cloud computing, and artificial intelligence (AI) are disciplines converging toward the use of this technology for active surveillance in poultry farms. Other advances in human and veterinary NGS sequencing are likely to be adaptable to avian species in the future.

Evaluation of extraction and enrichment methods for recovery of respiratory RNA viruses in a metagenomics approach

Viral metagenomics is increasingly applied in viral detection and virome characterization. Different extraction and enrichment techniques may be adopted, however, reports on their effective influence on viral recovery is often conflicting. Using a three step enrichment steps, the effect of three extraction kits and the influence of DNase treatment with or without rRNA removal for respiratory RNA virus recovery from nasopharyngeal swab samples was evaluated. The viral cocktail containing six different RNA viruses pooled in equal volume were subjected to the different extraction and enrichment methods, sequenced using the Illumina MiSeq, and analysed using Genome Detective. The PureLink® Viral RNA/DNA Mini Kit (PureLink) was highly efficient with better recovery of all the viral agents in the cocktail. The use of rRNA treatment resulted in increased viral recovery with PureLink and QIAamp® Viral RNA Mini kit, while having comparable recovery rate as DNase only with the QIAamp® MinElute Virus Spin Kit. The observed low reads and genome coverage of some of the viruses could be attributed to their low abundance. Depending on sample matrix, extraction choice and enrichment strategy may influence recovery of respiratory RNA virus in metagenomics studies, therefore individual evaluation and adoption may be necessary for a robust result.

A novel magnetic fluid for ultra-fast and highly efficient extraction of N1-methylnicotinamide in urine samples

Among the three pre-treatment materials, Fe3O4@HPMC@DMSA NPs were selected to be the best material and were used to perform MSPE-HPLC-UV.

Next Generation Sequencing Approaches to Characterize the Respiratory Tract Virome

The COVID-19 pandemic and heightened perception of the risk of emerging viral infections have boosted the efforts to better understand the virome or complete repertoire of viruses in health and disease, with a focus on infectious respiratory diseases. Next-generation sequencing (NGS) is widely used to study microorganisms, allowing the elucidation of bacteria and viruses inhabiting different body systems and identifying new pathogens. However, NGS studies suffer from a lack of standardization, in particular, due to various methodological approaches and no single format for processing the results. Here, we review the main methodological approaches and key stages for studies of the human virome, with an emphasis on virome changes during acute respiratory viral infection, with applications for clinical diagnostics and epidemiologic analyses.

Improvements in metagenomic virus detection by simple pretreatment methods

Early detection of pathogens at the point of care helps reduce the threats to human and animal health from emerging pathogens. Initially, the disease-causing agent will be unknown and needs to be identified; this often requires specific laboratory facilities. Here we describe the development of an unbiased detection assay for RNA and DNA viruses using metagenomic Nanopore sequencing and simple methods that can be transferred into a field setting. Human clinical samples containing the RNA virus SARS-CoV-2 or the DNA viruses human papillomavirus (HPV) and molluscum contagiosum virus (MCV) were used as a test of concept. Firstly, the virus detection potential was optimized by investigating different pretreatments for reducing non-viral nucleic acid components. DNase I pretreatment followed by filtration increased the proportion of SARS-CoV-2 sequenced reads > 500-fold compared with no pretreatments. This was sufficient to achieve virus detection with high confidence and allowed variant identification. Next, we tested individual SARS-CoV-2 samples with various viral loads (measured as CT-values determined by RT-qPCR). Lastly, we tested the assay on clinical samples containing the DNA virus HPV and co-infection with MCV to show the assay's detection potential for DNA viruses. This protocol is fast (same day results). We hope to apply this method in other settings for point of care detection of virus pathogens, thus eliminating the need for transport of infectious samples, cold storage and a specialized laboratory.

Comprehensive Evaluation of RNA and DNA Viromic Methods Based on Species Richness and Abundance Analyses Using Marmot Rectal Samples

Viral metagenomics is the most powerful tool to profile viromic composition for a given sample. Different viromic methods, including amplification-free ones, have been developed, but choosing them for different purposes requires comprehensive benchmarks. Here, we assessed the performance of four routinely used methods, i.e., multiple displacement amplification (MDA), direct metagenomic sequencing (MTG), sequence-independent single-primer amplification (SIA), and metatranscriptomic sequencing (MTT), using marmot rectal samples as the templates spiked with five known viruses of different genome types. The obtained clean data were differently contaminated by host and bacterial genomes, resulting in MDA having the most, with ~72.1%, but MTT had only ~7.5% data, useful for follow-up viromic analysis. MDA showed a broader spectrum with higher efficiency to profile the DNA virome, and MTT captured almost all RNA viruses with extraordinary sensitivity; hence, they are advisable in richness-based viromic studies. MTG was weak in capturing single-stranded DNA viruses, and SIA could detect both RNA and DNA viruses but with high randomness. Due to biases to certain types of viruses, the four methods caused different alterations to species abundance compared to the initial virus composition. SIA and MDA introduced greater stochastic errors to relative abundances of species, genus, and family taxa, whereas the two amplification-free methods were more tolerant toward such errors and thus are recommendable in abundance-based analyses. In addition, genus taxon is a compromising analytic level that ensures technically supported and biologically and/or ecologically meaningful viromic conclusions. IMPORTANCE Viral metagenomics can be roughly divided into species richness-based studies and species abundance-based analyses. Viromic methods with different principles have been developed, but rational selection of these techniques according to different purposes requires comprehensive understanding of their properties. By assessing the four most widely used methods using template samples, we found that multiple displacement amplification (MDA) and metatranscriptomic sequencing (MTT) are advisable for species richness-based viromic studies, as they show excellent efficiency to detect DNA and RNA viruses. Meanwhile, metagenomic sequencing (MTG) and MTT are more compatible with stochastic errors of methods introduced into relative abundance of viromic taxa and hence are rational choices in species abundance-based analyses. This study also highlights that MTG needs to tackle host genome contamination and ameliorate the capacity to detect single-stranded DNA viruses in the future, and the MTT method requires an improvement in bacterial rRNA depletion prior to library preparation.

RNA Viruses in Aquatic Ecosystems through the Lens of Ecological Genomics and Transcriptomics

Massive amounts of data from nucleic acid sequencing have changed our perspective about diversity and dynamics of marine viral communities. Here, we summarize recent metatranscriptomic and metaviromic studies targeting predominantly RNA viral communities. The analysis of RNA viromes reaffirms the abundance of lytic (+) ssRNA viruses of the order Picornavirales, but also reveals other (+) ssRNA viruses, including RNA bacteriophages, as important constituents of extracellular RNA viral communities. Sequencing of dsRNA suggests unknown diversity of dsRNA viruses. Environmental metatranscriptomes capture the dynamics of ssDNA, dsDNA, ssRNA, and dsRNA viruses simultaneously, unravelling the full complexity of viral dynamics in the marine environment. RNA viruses are prevalent in large size fractions of environmental metatranscriptomes, actively infect marine unicellular eukaryotes larger than 3 µm, and can outnumber bacteriophages during phytoplankton blooms. DNA and RNA viruses change abundance on hourly timescales, implying viral control on a daily temporal basis. Metatranscriptomes of cultured protists host a diverse community of ssRNA and dsRNA viruses, often with multipartite genomes and possibly persistent intracellular lifestyles. We posit that RNA viral communities might be more diverse and complex than formerly anticipated and that the influence they exert on community composition and global carbon flows in aquatic ecosystems may be underestimated.

Direct sequencing of measles virus complete genomes in the midst of a large-scale outbreak

Measles outbreaks escalated globally despite worldwide elimination efforts. Molecular epidemiological investigations utilizing partial measles virus (MeV) genomes are challenged by reduction in global genotypes and low evolutionary rates. Greater resolution was reached using MeV complete genomes, however time and costs limit the application to numerous samples. We developed an approach to unbiasedly sequence complete MeV genomes directly from patient urine samples. Samples were enriched for MeV using filtration or nucleases and the minimal number of sequence reads to allocate per sample based on its MeV content was assessed using in-silico reduction of sequencing depth. Application of limited-resource sequencing to treated MeV-positive samples demonstrated that 1–5 million sequences for samples with high/medium MeV quantities and 10–15 million sequences for samples with lower MeV quantities are sufficient to obtain >98% MeV genome coverage and over X50 average depth. This approach enables real-time high-resolution molecular epidemiological investigations of large-scale MeV outbreaks.

Microbial metagenomic approach uncovers the first rabbit haemorrhagic disease virus genome in Sub-Saharan Africa

Rabbit Haemorrhagic Disease (RHD) causes high morbidity and mortality in rabbits and hares. Here, we report the first genomic characterization of lagovirus GI.2 virus in domestic rabbits from sub-Saharan Africa. We used an unbiased microbial metagenomic Next Generation Sequencing (mNGS) approach to diagnose the pathogen causing the suspected outbreak of RHD in Ibadan, Nigeria. The liver, spleen, and lung samples of five rabbits from an outbreak in 2 farms were analyzed. The mNGS revealed one full and two partial RHDV2 genomes on both farms. Phylogenetic analysis showed close clustering with RHDV2 lineages from Europe (98.6% similarity with RHDV2 in the Netherlands, and 99.1 to 100% identity with RHDV2 in Germany), suggesting potential importation. Subsequently, all the samples were confirmed by RHDV virus-specific RT-PCR targeting the VP60 gene with the expected band size of 398 bp for the five rabbits sampled. Our findings highlight the need for increased genomic surveillance of RHDV2 to track its origin, understand its diversity and to inform public health policy in Nigeria, and Sub-Saharan Africa.

Enrichment approach for unbiased sequencing of respiratory syncytial virus directly from clinical samples

Background: Nasopharyngeal samples contain higher quantities of bacterial and host nucleic acids relative to viruses; presenting challenges during virus metagenomics sequencing, which underpins agnostic sequencing protocols. We aimed to develop a viral enrichment protocol for unbiased whole-genome sequencing of respiratory syncytial virus (RSV) from nasopharyngeal samples using the Oxford Nanopore Technology (ONT) MinION platform. Methods: We assessed two protocols using RSV positive samples. Protocol 1 involved physical pre-treatment of samples by centrifugal processing before RNA extraction, while Protocol 2 entailed direct RNA extraction without prior enrichment. Concentrates from Protocol 1 and RNA extracts from Protocol 2 were each divided into two fractions; one was DNase treated while the other was not. RNA was then extracted from both concentrate fractions per sample and RNA from both protocols converted to cDNA, which was then amplified using the tagged Endoh primers through Sequence-Independent Single-Primer Amplification (SISPA) approach, a library prepared, and sequencing done. Statistical significance during analysis was tested using the Wilcoxon signed-rank test. Results: DNase-treated fractions from both protocols recorded significantly reduced host and bacterial contamination unlike the untreated fractions (in each protocol p<0.01). Additionally, DNase treatment after RNA extraction (Protocol 2) enhanced host and bacterial read reduction compared to when done before (Protocol 1). However, neither protocol yielded whole RSV genomes. Sequenced reads mapped to parts of the nucleoprotein (N gene) and polymerase complex (L gene) from Protocol 1 and 2, respectively. Conclusions: DNase treatment was most effective in reducing host and bacterial contamination, but its effectiveness improved if done after RNA extraction than before. We attribute the incomplete genome segments to amplification biases resulting from the use of short length random sequence (6 bases) in tagged Endoh primers. Increasing the length of the random nucleotides from six hexamers to nine or 12 in future studies may reduce the coverage biases.

A Versatile Processing Workflow to Enable Pathogen Detection in Clinical Samples from Organs Using VIDISCA

In recent years, refined molecular methods coupled with powerful high throughput sequencing technologies have increased the potential of virus discovery in clinical samples. However, host genetic material remains a complicating factor that interferes with discovery of novel viruses in solid tissue samples as the relative abundance of the virus material is low. Physical enrichment processing methods, although usually complicated, labor-intensive, and costly, have proven to be successful for improving sensitivity of virus detection in complex samples. In order to further increase detectability, we studied the application of fast and simple high-throughput virus enrichment methods on tissue homogenates. Probe sonication in high EDTA concentrations, organic extraction with Vertrel™ XF, or a combination of both, were applied prior to chromatography-like enrichment using Capto™ Core 700 resin, after which effects on virus detection sensitivity by the VIDISCA method were determined. Sonication in the presence of high concentrations of EDTA showed the best performance with an increased proportion of viral reads, up to 9.4 times, yet minimal effect on the host background signal. When this sonication procedure in high EDTA concentrations was followed by organic extraction with Vertrel™ XF and two rounds of core bead chromatography enrichment, an increase up to 10.5 times in the proportion of viral reads in the processed samples was achieved, with reduction of host background sequencing. We present a simple and semi-high-throughput method that can be used to enrich homogenized tissue samples for viral reads.

A library preparation optimized for metagenomics of RNA viruses

Our understanding of the viral communities associated to animals has not yet reached the level attained on the bacteriome. This situation is due to, among others, technical challenges in adapting metagenomics using high-throughput sequencing to the study of RNA viromes in animals. Although important developments have been achieved in most steps of viral metagenomics, there is yet a key step that has received little attention: the library preparation. This situation differs from bacteriome studies in which developments in library preparation have largely contributed to the democratisation of metagenomics. Here, we present a library preparation optimized for metagenomics of RNA viruses from insect vectors of viral diseases. The library design allows a simple PCR-based preparation, such as those routinely used in bacterial metabarcoding, that is adapted to shotgun sequencing as required in viral metagenomics. We first optimized our library preparation using mock viral communities and then validated a full metagenomic approach incorporating our preparation in two pilot studies with field-caught insect vectors; one including a comparison with a published metagenomic protocol. Our approach provided a fold increase in virus-like sequences compared to other studies, and nearly-full genomes from new virus species. Moreover, our results suggested conserved trends in virome composition within a population of a mosquito species. Finally, the sensitivity of our approach was compared to a commercial diagnostic PCR for the detection of an arbovirus in field-caught insect vectors. Our approach could facilitate studies on viral communities from animals and the democratization of metagenomics in community ecology of viruses.

Deep viral blood metagenomics reveals extensive anellovirus diversity in healthy humans

Human blood metagenomics has revealed the presence of different types of viruses in apparently healthy subjects. By far, anelloviruses constitute the viral family that is more frequently found in human blood, although amplification biases and contaminations pose a major challenge in this field. To investigate this further, we subjected pooled plasma samples from 120 healthy donors in Spain to high-speed centrifugation, RNA and DNA extraction, random amplification, and massive parallel sequencing. Our results confirm the extensive presence of anelloviruses in such samples, which represented nearly 97% of the total viral sequence reads obtained. We assembled 114 different viral genomes belonging to this family, revealing remarkable diversity. Phylogenetic analysis of ORF1 suggested 28 potentially novel anellovirus species, 24 of which were validated by Sanger sequencing to discard artifacts. These findings underscore the importance of implementing more efficient purification procedures that enrich the viral fraction as an essential step in virome studies and question the suggested pathological role of anelloviruses.

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.

Comparative Metagenomics of Palearctic and Neotropical Avian Cloacal Viromes Reveal Geographic Bias in Virus Discovery

Our understanding about viruses carried by wild animals is still scarce. The viral diversity of wildlife may be best described with discovery-driven approaches to the study of viral diversity that broaden research efforts towards non-canonical hosts and remote geographic regions. Birds have been key organisms in the transmission of viruses causing important diseases, and wild birds are threatened by viral spillovers associated with human activities. However, our knowledge of the avian virome may be biased towards poultry and highly pathogenic diseases. We describe and compare the fecal virome of two passerine-dominated bird assemblages sampled in a remote Neotropical rainforest in French Guiana (Nouragues Natural Reserve) and a Mediterranean forest in central Spain (La Herrería). We used metagenomic data to quantify the degree of functional and genetic novelty of viruses recovered by examining if the similarity of the contigs we obtained to reference sequences differed between both locations. In general, contigs from Nouragues were significantly less similar to viruses in databases than contigs from La Herrería using Blastn but not for Blastx, suggesting that pristine regions harbor a yet unknown viral diversity with genetically more singular viruses than more studied areas. Additionally, we describe putative novel viruses of the families Picornaviridae, Reoviridae and Hepeviridae. These results highlight the importance of wild animals and remote regions as sources of novel viruses that substantially broaden the current knowledge of the global diversity of viruses.

Phages and their potential to modulate the microbiome and immunity

Bacteriophages (hence termed phages) are viruses that target bacteria and have long been considered as potential future treatments against antibiotic-resistant bacterial infection. However, the molecular nature of phage interactions with bacteria and the human host has remained elusive for decades, limiting their therapeutic application. While many phages and their functional repertoires remain unknown, the advent of next-generation sequencing has increasingly enabled researchers to decode new lytic and lysogenic mechanisms by which they attack and destroy bacteria. Furthermore, the last decade has witnessed a renewed interest in the utilization of phages as therapeutic vectors and as a means of targeting pathogenic or commensal bacteria or inducing immunomodulation. Importantly, the narrow host range, immense antibacterial repertoire, and ease of manipulating phages may potentially allow for their use as targeted modulators of pathogenic, commensal and pathobiont members of the microbiome, thereby impacting mammalian physiology and immunity along mucosal surfaces in health and in microbiome-associated diseases. In this review, we aim to highlight recent advances in phage biology and how a mechanistic understanding of phage-bacteria-host interactions may facilitate the development of novel phage-based therapeutics. We provide an overview of the challenges of the therapeutic use of phages and how these could be addressed for future use of phages as specific modulators of the human microbiome in a variety of infectious and noncommunicable human diseases.

Evaluation of a commercial exogenous internal process control for diagnostic RNA virus metagenomics from different animal clinical samples

Metagenomic next generation sequencing (mNGS) is increasingly recognized as an important complementary tool to targeted human and animal infectious disease diagnostics. It is, however, sensitive to biases and errors that are currently not systematically evaluated by the implementation of quality controls (QC) for the diagnostic use of mNGS. We evaluated a commercial reagent (Mengovirus extraction control kit, CeraamTools, bioMérieux) as an exogenous internal control for mNGS. It validates the integrity of reagents and workflow, the efficient isolation of viral nucleic acids and the absence of inhibitors in individual samples (verified using a specific qRT-PCR). Moreover, it validates the efficient generation of viral sequence data in individual samples (verified by normalized mengoviral read counts in the metagenomic analysis). We show that when using a completely random metagenomics workflow: (1) Mengovirus RNA can be reproducibly detected in different animal sample types (swine feces and sera, wild bird cloacal swabs), except for tissue samples (swine lung); (2) the Mengovirus control kit does not contain other contaminating viruses that may affect metagenomic experiments (using a cutoff of minimum 1 Kraken classified read per million (RPM)); (3) the addition of 2.17 × 10⁶Mengovirus copies/mL of sample does not affect the virome composition of pig fecal samples or wild bird cloacal swab samples; (4) Mengovirus Cq values (using as cutoff the upper limit of the 99 % confidence interval of Cq values for a given sample matrix) allow the identification of samples with poor viral RNA extraction or high inhibitor load; (5) Mengovirus normalized read counts (cutoff RPM > 1) allow the identification of samples where the viral sequences are outcompeted by host or bacterial target sequences in the random metagenomic workflow. The implementation of two QC testing points, a first one after RNA extraction (Mengoviral qRT-PCR) and a second one after metagenomic data analysis provide valuable information for the validation of individual samples and results. Their implementation in addition to external controls validating runs or experiments should be carefully considered for a given sample type and workflow.

Sensitivity and specificity of metatranscriptomics as an arbovirus surveillance tool

The ability to identify all the viruses within a sample makes metatranscriptomic sequencing an attractive tool to screen mosquitoes for arboviruses. Practical application of this technique, however, requires a clear understanding of its analytical sensitivity and specificity. To assess this, five dilutions (1:1, 1:20, 1:400, 1:8,000 and 1:160,000) of Ross River virus (RRV) and Umatilla virus (UMAV) isolates were spiked into subsamples of a pool of 100 Culex australicus mosquitoes. The 1:1 dilution represented the viral load of one RRV-infected mosquito in a pool of 100 mosquitoes. The subsamples underwent nucleic acid extraction, mosquito-specific ribosomal RNA depletion, and Illumina HiSeq sequencing. The viral load of the subsamples was also measured using reverse transcription droplet digital PCR (RT-ddPCR) and quantitative PCR (RT-qPCR). Metatranscriptomic sequencing detected both RRV and UMAV in the 1:1, 1:20 and 1:400 subsamples. A high specificity was achieved, with 100% of RRV and 99.6% of UMAV assembled contigs correctly identified. Metatranscriptomic sequencing was not as sensitive as RT-qPCR or RT-ddPCR; however, it recovered whole genome information and detected 19 other viruses, including four first detections for Australia. These findings will assist arbovirus surveillance programs in utilising metatranscriptomics in routine surveillance activities to enhance arbovirus detection.

dsRNA-Seq: Identification of Viral Infection by Purifying and Sequencing dsRNA

RNA viruses are a major source of emerging and re-emerging infectious diseases around the world. We developed a method to identify RNA viruses that is based on the fact that RNA viruses produce double-stranded RNA (dsRNA) while replicating. Purifying and sequencing dsRNA from the total RNA isolated from infected tissue allowed us to recover dsRNA virus sequences and replicated sequences from single-stranded RNA (ssRNA) viruses. We refer to this approach as dsRNA-Seq. By assembling dsRNA sequences into contigs we identified full length or partial RNA viral genomes of varying genome types infecting mammalian culture samples, identified a known viral disease agent in laboratory infected mice, and successfully detected naturally occurring RNA viral infections in reptiles. Here, we show that dsRNA-Seq is a preferable method for identifying viruses in organisms that don’t have sequenced genomes and/or commercially available rRNA depletion reagents. In addition, a significant advantage of this method is the ability to identify replicated viral sequences of ssRNA viruses, which is useful for distinguishing infectious viral agents from potential noninfectious viral particles or contaminants.

Fecal Viral Diversity of Captive and Wild Tasmanian Devils Characterized Using Virion-Enriched Metagenomics and Metatranscriptomics

The Tasmanian devil is an endangered carnivorous marsupial threatened by devil facial tumor disease (DFTD). While research on DFTD has been extensive, little is known about viruses in devils and whether any are of potential conservation relevance for this endangered species. Using both metagenomics based on virion enrichment and sequence-independent amplification (virion-enriched metagenomics) and metatranscriptomics based on bulk RNA sequencing, we characterized and compared the fecal viromes of captive and wild devils. A total of 54 fecal samples collected from two captive and four wild populations were processed for virome characterization using both approaches. In total, 24 novel marsupial-related viruses, comprising a sapelovirus, astroviruses, rotaviruses, picobirnaviruses, parvoviruses, papillomaviruses, polyomaviruses, and a gammaherpesvirus, were identified, as well as known mammalian pathogens such as rabbit hemorrhagic disease virus 2. Captive devils showed significantly lower viral diversity than wild devils. Comparison of the two virus discovery approaches revealed substantial differences in the number and types of viruses detected, with metatranscriptomics better suited for RNA viruses and virion-enriched metagenomics largely identifying more DNA viruses. Thus, the viral communities revealed by virion-enriched metagenomics and metatranscriptomics were not interchangeable and neither approach was able to detect all viruses present. An integrated approach using both virion-enriched metagenomics and metatranscriptomics constitutes a powerful tool for obtaining a complete overview of both the taxonomic and functional profiles of viral communities within a sample.IMPORTANCE The Tasmanian devil is an iconic Australian marsupial that has suffered an 80% population decline due to a contagious cancer, devil facial tumor disease, along with other threats. Until now, viral discovery in this species has been confined to one gammaherpesvirus (dasyurid herpesvirus 2 [DaHV-2]), for which captivity was identified as a significant risk factor. Our discovery of 24 novel marsupial-associated RNA and DNA viruses, and that viral diversity is lower in captive than in wild devils, has greatly expanded our knowledge of gut-associated viruses in devils and provides important baseline information that will contribute to the conservation and captive management of this endangered species. Our results also revealed that a combination of virion-enriched metagenomics and metatranscriptomics may be a more comprehensive approach for virome characterization than either method alone. Our results thus provide a springboard for continuous improvements in the way we study complex viral communities.

Fecal Viral Diversity of Captive and Wild Tasmanian Devils Characterized Using Virion-Enriched Metagenomics and Metatranscriptomics

The Tasmanian devil is an iconic Australian marsupial that has suffered an 80% population decline due to a contagious cancer, devil facial tumor disease, along with other threats. Until now, viral discovery in this species has been confined to one gammaherpesvirus (dasyurid herpesvirus 2 [DaHV-2]), for which captivity was identified as a significant risk factor. Our discovery of 24 novel marsupial-associated RNA and DNA viruses, and that viral diversity is lower in captive than in wild devils, has greatly expanded our knowledge of gut-associated viruses in devils and provides important baseline information that will contribute to the conservation and captive management of this endangered species. Our results also revealed that a combination of virion-enriched metagenomics and metatranscriptomics may be a more comprehensive approach for virome characterization than either method alone. Our results thus provide a springboard for continuous improvements in the way we study complex viral communities.

The Tasmanian devil is an endangered carnivorous marsupial threatened by devil facial tumor disease (DFTD). While research on DFTD has been extensive, little is known about viruses in devils and whether any are of potential conservation relevance for this endangered species. Using both metagenomics based on virion enrichment and sequence-independent amplification (virion-enriched metagenomics) and metatranscriptomics based on bulk RNA sequencing, we characterized and compared the fecal viromes of captive and wild devils. A total of 54 fecal samples collected from two captive and four wild populations were processed for virome characterization using both approaches. In total, 24 novel marsupial-related viruses, comprising a sapelovirus, astroviruses, rotaviruses, picobirnaviruses, parvoviruses, papillomaviruses, polyomaviruses, and a gammaherpesvirus, were identified, as well as known mammalian pathogens such as rabbit hemorrhagic disease virus 2. Captive devils showed significantly lower viral diversity than wild devils. Comparison of the two virus discovery approaches revealed substantial differences in the number and types of viruses detected, with metatranscriptomics better suited for RNA viruses and virion-enriched metagenomics largely identifying more DNA viruses. Thus, the viral communities revealed by virion-enriched metagenomics and metatranscriptomics were not interchangeable and neither approach was able to detect all viruses present. An integrated approach using both virion-enriched metagenomics and metatranscriptomics constitutes a powerful tool for obtaining a complete overview of both the taxonomic and functional profiles of viral communities within a sample.

IMPORTANCE The Tasmanian devil is an iconic Australian marsupial that has suffered an 80% population decline due to a contagious cancer, devil facial tumor disease, along with other threats. Until now, viral discovery in this species has been confined to one gammaherpesvirus (dasyurid herpesvirus 2 [DaHV-2]), for which captivity was identified as a significant risk factor. Our discovery of 24 novel marsupial-associated RNA and DNA viruses, and that viral diversity is lower in captive than in wild devils, has greatly expanded our knowledge of gut-associated viruses in devils and provides important baseline information that will contribute to the conservation and captive management of this endangered species. Our results also revealed that a combination of virion-enriched metagenomics and metatranscriptomics may be a more comprehensive approach for virome characterization than either method alone. Our results thus provide a springboard for continuous improvements in the way we study complex viral communities.

Metagenomic analysis of viral nucleic acid extraction methods in respiratory clinical samples

BACKGROUND

Numerous protocols for viral enrichment and genome amplification have been created. However, the direct identification of viral genomes from clinical specimens using next-generation sequencing (NGS) still has its challenges. As a selected viral nucleic acid extraction method may determine the sensitivity and reliability of NGS, it is still valuable to evaluate the extraction efficiency of different extraction kits using clinical specimens directly.

RESULTS

In this study, we performed qRT-PCR and viral metagenomic analysis of the extraction efficiency of four commonly used Qiagen extraction kits: QIAamp Viral RNA Mini Kit (VRMK), QIAamp MinElute Virus Spin Kit (MVSK), RNeasy Mini Kit (RMK), and RNeasy Plus Micro Kit (RPMK), using a mixed respiratory clinical sample without any pre-treatment. This sample contained an adenovirus (ADV), influenza virus A (Flu A), human parainfluenza virus 3 (PIV3), human coronavirus OC43 (OC43), and human metapneumovirus (HMPV). The quantity and quality of the viral extracts were significantly different among these kits. The highest threshold cycle(Ct)values for ADV and OC43 were obtained by using the RPMK. The MVSK had the lowest Ct values for ADV and PIV3. The RMK revealed the lowest detectability for HMPV and PIV3. The most effective rate of NGS data at 67.47% was observed with the RPMK. The other three kits ranged between 12.1-26.79% effectiveness rates for the NGS data. Most importantly, compared to the other three kits the highest proportion of non-host reads was obtained by the RPMK. The MVSK performed best with the lowest Ct value of 20.5 in the extraction of ADV, while the RMK revealed the best extraction efficiency by NGS analysis.

CONCLUSIONS

The evaluation of viral nucleic acid extraction efficiency is different between NGS and qRT-PCR analysis. The RPMK was most applicable for the metagenomic analysis of viral RNA and enabled more sensitive identification of the RNA virus genome in respiratory clinical samples. In addition, viral RNA extraction kits were also applicable for metagenomic analysis of the DNA virus. Our results highlighted the importance of nucleic acid extraction kit selection, which has a major impact on the yield and number of viral reads by NGS analysis. Therefore, the choice of extraction method for a given viral pathogen needs to be carefully considered.

A Versatile Sample Processing Workflow for Metagenomic Pathogen Detection

Metagenomics is currently the only generic method for pathogen detection. Starting from RNA allows the assessment of the whole sample community including RNA viruses. Here we present our modular concerted protocol for sample processing for diagnostic metagenomics analysis of human, animal, and food samples. The workflow does not rely on dedicated amplification steps at any stage in the process and, in contrast to published methods, libraries prepared accordingly will yield only minute amounts of unclassifiable reads. We confirmed the performance of the approach using a spectrum of pathogen/matrix-combinations showing it has the potential to become a commonly usable analytical framework.

Simple and label-free pathogen enrichment via homobifunctional imidoesters using a microfluidic (SLIM) system for ultrasensitive pathogen detection in various clinical specimens

Diseases caused by pathogenic microorganisms including bacteria and viruses can cause serious medical issues including death and result in huge economic losses. Despite the myriad of recent advances in the rapid and accurate detection of pathogens, large volume clinical samples with a low concentration of pathogens continue to present challenges for diagnosis and surveillance. We here report a simple and label-free approach via homobifunctional imidoesters (HIs) with a microfluidic platform (SLIM) to efficiently enrich and extract pathogens at low concentrations from clinical samples. The SLIM system consists of an assembled double microfluidic chip for streamlining large volume processing and HIs for capturing pathogens and isolating nucleic acids by both electrostatic and covalent interaction without a chaotropic detergent or bulky instruments. The SLIM system significantly increases the enrichment and extraction rate of pathogens (up to 80% at 10 CFU (colony forming unit) in a 1 mL volume within 50 min). We demonstrated its clinical utility in large sample volumes from 46 clinical specimens including environmental swabs, saliva, and blood plasma. The SLIM system showed higher sensitivity with these samples and could detect pathogens that were below the threshold of detection with other methods. Finally, by combining our SLIM approach with an isothermal optical sensor, pathogens could be detected at a very high sensitivity in blood plasma samples within 80 min via enrichment, extraction and detection steps. Our SLIM system thus provides a simple, reliable, cost-effective and ultrasensitive pathogen diagnosis platform for use with large volume clinical samples and would thus have significant utility for various infectious diseases.

•

SLIM system significantly increases the enrichment and extraction rate of pathogens.

•

Demonstrated its clinical utility in large sample volumes from 46 clinical specimens.

•

A simple, reliable, cost-effective and ultrasensitive pathogen diagnosis platform.

Complete coding sequence of a novel picorna-like virus in a blackbird infected with Usutu virus

Using random high-throughput RNA sequencing, the complete coding sequence of a novel picorna-like virus (a 9,228-nt contig containing 212,202 reads) was determined from a blackbird (Turdus merula) infected with Usutu virus. This sequence shares only 36% amino acid sequence identity with its closest homolog, arivirus 1, (an unclassified member of the order Picornavirales), and shares its dicistronic genome arrangement. The new virus was therefore tentatively named "blackbird arilivirus" (ari-like virus). The nearly complete genome sequence consists of at least 9,228 nt and contains two open reading frames (ORFs) encoding the nonstructural polyprotein (2235 amino acids) and structural polyprotein (769 amino acids). Two TaqMan RT-qPCR assays specific for ORF1 confirmed the presence of high levels of this novel virus in the original sample. Nucleotide composition analysis suggests that blackbird arilivirus is of dietary (plant) origin.

Protocol for Generating Infectious RNA Viromes from Complex Biological Samples

This chapter proposes a simple, standardized protocol for generating RNA viromes from complex host-associated biological samples of various origins. Compared to other existing protocols to generate RNA viromes, this protocol preserves the infectivity of viral particles and allows for downstream applications such as viral characterization and isolation tests.

Identification of sapovirus GV.2, astrovirus VA3 and novel anelloviruses in serum from patients with acute hepatitis of unknown aetiology

Hepatitis is a general term meaning inflammation of the liver, which can be caused by a variety of viruses. However, a substantial number of cases remain with unknown aetiology. We analysed the serum of patients with clinical signs of hepatitis using a metagenomics approach to characterize their viral species composition. Four pools of patients with hepatitis without identified aetiological agents were evaluated. Additionally, one pool of patients with hepatitis E (HEV) and pools of healthy volunteers were included as controls. A high diversity of anelloviruses, including novel sequences, was found in pools from patients with hepatitis of unknown aetiology. Moreover, viruses recently associated with gastroenteritis as sapovirus GV.2 and astrovirus VA3 were also detected only in those pools. Besides, most of the HEV genome was recovered from the HEV pool. Finally, GB virus C and human endogenous retrovirus were found in the HEV and healthy pools. Our study provides an overview of the virome in serum from hepatitis patients suggesting a potential role of these viruses not previously described in cases of hepatitis. However, further epidemiologic studies are necessary to confirm their contribution to the development of hepatitis.

Optimization and validation of sample preparation for metagenomic sequencing of viruses in clinical samples

BACKGROUND

Sequence-specific PCR is the most common approach for virus identification in diagnostic laboratories. However, as specific PCR only detects pre-defined targets, novel virus strains or viruses not included in routine test panels will be missed. Recently, advances in high-throughput sequencing allow for virus-sequence-independent identification of entire virus populations in clinical samples, yet standardized protocols are needed to allow broad application in clinical diagnostics. Here, we describe a comprehensive sample preparation protocol for high-throughput metagenomic virus sequencing using random amplification of total nucleic acids from clinical samples.

RESULTS

In order to optimize metagenomic sequencing for application in virus diagnostics, we tested different enrichment and amplification procedures on plasma samples spiked with RNA and DNA viruses. A protocol including filtration, nuclease digestion, and random amplification of RNA and DNA in separate reactions provided the best results, allowing reliable recovery of viral genomes and a good correlation of the relative number of sequencing reads with the virus input. We further validated our method by sequencing a multiplexed viral pathogen reagent containing a range of human viruses from different virus families. Our method proved successful in detecting the majority of the included viruses with high read numbers and compared well to other protocols in the field validated against the same reference reagent. Our sequencing protocol does work not only with plasma but also with other clinical samples such as urine and throat swabs.

CONCLUSIONS

The workflow for virus metagenomic sequencing that we established proved successful in detecting a variety of viruses in different clinical samples. Our protocol supplements existing virus-specific detection strategies providing opportunities to identify atypical and novel viruses commonly not accounted for in routine diagnostic panels.

Efficient and unbiased metagenomic recovery of RNA virus genomes from human plasma samples

RNA viruses cause significant human pathology and are responsible for the majority of emerging zoonoses. Mainstream diagnostic assays are challenged by their intrinsic diversity, leading to false negatives and incomplete characterisation. New sequencing techniques are expanding our ability to agnostically interrogate nucleic acids within diverse sample types, but in the clinical setting are limited by overwhelming host material and ultra-low target frequency. Through selective host RNA depletion and compensatory protocol adjustments for ultra-low RNA inputs, we are able to detect three major blood-borne RNA viruses – HIV, HCV and HEV. We recovered complete genomes and up to 43% of the genome from samples with viral loads of 10⁴ and 10³ IU/ml respectively. Additionally, we demonstrated the utility of this method in detecting and characterising members of diverse RNA virus families within a human plasma background, some present at very low levels. By applying this method to a patient sample series, we have simultaneously determined the full genome of both a novel subtype of HCV genotype 6, and a co-infecting human pegivirus. This method builds upon earlier RNA metagenomic techniques and can play an important role in the surveillance and diagnostics of blood-borne viruses.

Evolution of selective-sequencing approaches for virus discovery and virome analysis

•

Description of virus enrichment techniques for metagenomics based virome analysis.

•

Usefulness of recently developed virome capture sequencing techniques.

•

Perspective on negative and positive selection approaches for virome analysis.

Recent advances in sequencing technologies have transformed the field of virus discovery and virome analysis. Once mostly confined to the traditional Sanger sequencing based individual virus discovery, is now entirely replaced by high throughput sequencing (HTS) based virus metagenomics that can be used to characterize the nature and composition of entire viromes. To better harness the potential of HTS for the study of viromes, sample preparation methodologies use different approaches to exclude amplification of non-viral components that can overshadow low-titer viruses. These virus-sequence enrichment approaches mostly focus on the sample preparation methods, like enzymatic digestion of non-viral nucleic acids and size exclusion of non-viral constituents by column filtration, ultrafiltration or density gradient centrifugation. However, recently a new approach of virus-sequence enrichment called virome-capture sequencing, focused on the amplification or HTS library preparation stage, was developed to increase the ability of virome characterization. This new approach has the potential to further transform the field of virus discovery and virome analysis, but its technical complexity and sequence-dependence warrants further improvements. In this review we discuss the different methods, their applications and evolution, for selective sequencing based virome analysis and also propose refinements needed to harness the full potential of HTS for virome analysis.

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.

Complete Coding Sequence of Usutu Virus Strain Gracula religiosa/U1609393/Belgium/2016 Obtained from the Brain Tissue of an Infected Captive Common Hill Myna (Gracula religiosa)

The complete and annotated coding sequence and partial noncoding sequence of an Usutu virus genome were sequenced from RNA extracted from a clinical brain tissue sample obtained from a common hill myna (Gracula religiosa), demonstrating close homology with Usutu viruses circulating in Europe.

The blood DNA virome in 8,000 humans

The characterization of the blood virome is important for the safety of blood-derived transfusion products, and for the identification of emerging pathogens. We explored non-human sequence data from whole-genome sequencing of blood from 8,240 individuals, none of whom were ascertained for any infectious disease. Viral sequences were extracted from the pool of sequence reads that did not map to the human reference genome. Analyses sifted through close to 1 Petabyte of sequence data and performed 0.5 trillion similarity searches. With a lower bound for identification of 2 viral genomes/100,000 cells, we mapped sequences to 94 different viruses, including sequences from 19 human DNA viruses, proviruses and RNA viruses (herpesviruses, anelloviruses, papillomaviruses, three polyomaviruses, adenovirus, HIV, HTLV, hepatitis B, hepatitis C, parvovirus B19, and influenza virus) in 42% of the study participants. Of possible relevance to transfusion medicine, we identified Merkel cell polyomavirus in 49 individuals, papillomavirus in blood of 13 individuals, parvovirus B19 in 6 individuals, and the presence of herpesvirus 8 in 3 individuals. The presence of DNA sequences from two RNA viruses was unexpected: Hepatitis C virus is revealing of an integration event, while the influenza virus sequence resulted from immunization with a DNA vaccine. Age, sex and ancestry contributed significantly to the prevalence of infection. The remaining 75 viruses mostly reflect extensive contamination of commercial reagents and from the environment. These technical problems represent a major challenge for the identification of novel human pathogens. Increasing availability of human whole-genome sequences will contribute substantial amounts of data on the composition of the normal and pathogenic human blood virome. Distinguishing contaminants from real human viruses is challenging.

Novel sequencing technologies offer insight into the virome in human samples. Here, we identify the viral DNA sequences in blood of over 8,000 individuals undergoing whole genome sequencing. This approach serves to identify 94 viruses; however, many are shown to reflect widespread DNA contamination of commercial reagents or of environmental origin. While this represents a significant limitation to reliably identify novel viruses infecting humans, we could confidently detect sequences and quantify abundance of 19 human viruses in 42% of individuals. Ancestry, sex, and age were important determinants of viral prevalence. This large study calls attention on the challenge of interpreting next generation sequencing data for the identification of novel viruses. However, it serves to categorize the abundance of human DNA viruses using an unbiased technique.

Tunable and label-free virus enrichment for ultrasensitive virus detection using carbon nanotube arrays

Viral infectious diseases can erupt unpredictably, spread rapidly, and ravage mass populations. Although established methods, such as polymerase chain reaction, virus isolation, and next-generation sequencing have been used to detect viruses, field samples with low virus count pose major challenges in virus surveillance and discovery. We report a unique carbon nanotube size-tunable enrichment microdevice (CNT-STEM) that efficiently enriches and concentrates viruses collected from field samples. The channel sidewall in the microdevice was made by growing arrays of vertically aligned nitrogen-doped multiwalled CNTs, where the intertubular distance between CNTs could be engineered in the range of 17 to 325 nm to accurately match the size of different viruses. The CNT-STEM significantly improves detection limits and virus isolation rates by at least 100 times. Using this device, we successfully identified an emerging avian influenza virus strain [A/duck/PA/02099/2012(H11N9)] and a novel virus strain (IBDV/turkey/PA/00924/14). Our unique method demonstrates the early detection of emerging viruses and the discovery of new viruses directly from field samples, thus creating a universal platform for effectively remediating viral infectious diseases.

Development and evaluation of a non-ribosomal random PCR and next-generation sequencing based assay for detection and sequencing of hand, foot and mouth disease pathogens

Background

Hand, foot and mouth disease (HFMD) has become a major public health problem across the Asia-Pacific region, and is commonly caused by enterovirus A71 (EV-A71) and coxsackievirus A6 (CV-A6), CV-A10 and CV-A16. Generating pathogen whole-genome sequences is essential for understanding their evolutionary biology. The frequent replacements among EV serotypes and a limited numbers of available whole-genome sequences hinder the development of overlapping PCRs for whole-genome sequencing.

We developed and evaluated a non-ribosomal random PCR (rPCR) and next-generation sequencing based assay for sequence-independent whole-genome amplification and sequencing of HFMD pathogens. A total of 16 EV-A71/CV-A6/CV-A10/CV-A16 PCR positive rectal/throat swabs (Cp values: 20.9–33.3) were used for assay evaluation.

Results

Our assay evidently outperformed the conventional rPCR in terms of the total number of EV-A71 reads and the percentage of EV-A71 reads: 2.6 % (1275/50,000 reads) vs. 0.1 % (31/50,000) and 6 % (3008/50,000) vs. 0.9 % (433/50,000) for two samples with Cp values of 30 and 26, respectively. Additionally the assay could generate genome sequences with the percentages of coverage of 94–100 % of 4 different enterovirus serotypes in 73 % of the tested samples, representing the first whole-genome sequences of CV-A6/10/16 from Vietnam, and could assign correctly serotyping results in 100 % of 24 tested specimens. In all but three the obtained consensuses of two replicates from the same sample were 100 % identical, suggesting that our assay is highly reproducible.

Conclusions

In conclusion, we have successfully developed a non-ribosomal rPCR and next-generation sequencing based assay for sensitive detection and direct whole-genome sequencing of HFMD pathogens from clinical samples.

Electronic supplementary material

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

Complete Genome Sequences of Three African Foot-and-Mouth Disease Viruses from Clinical Samples Isolated in 2009 and 2010

The complete genome sequences of three foot-and-mouth disease viruses (one virus of each serotype SAT1, SAT2 and O) were directly sequenced from RNA extracted from clinical bovine samples, demonstrating the feasibility of full-genome sequencing from strong positive samples taken from symptomatic animals.

Complete Genome Sequence of Bovine Polyomavirus Type 1 from Aborted Cattle, Isolated in Belgium in 2014

The complete and fully annotated genome sequence of a bovine polyomavirus type 1 (BPyV/BEL/1/2014) from aborted cattle was assembled from a metagenomics data set. The 4,697-bp circular dsDNA genome contains 6 protein-coding genes. Bovine polyomavirus is unlikely to be causally related to the abortion cases.

Viral metagenomics and blood safety

The characterization of the human blood-associated viral community (also called blood virome) is essential for epidemiological surveillance and to anticipate new potential threats for blood transfusion safety. Currently, the risk of blood-borne agent transmission of well-known viruses (HBV, HCV, HIV and HTLV) can be considered as under control in high-resource countries. However, other viruses unknown or unsuspected may be transmitted to recipients by blood-derived products. This is particularly relevant considering that a significant proportion of transfused patients are immunocompromised and more frequently subjected to fatal outcomes. Several measures to prevent transfusion transmission of unknown viruses have been implemented including the exclusion of at-risk donors, leukocyte reduction of donor blood, and physicochemical treatment of the different blood components. However, up to now there is no universal method for pathogen inactivation, which would be applicable for all types of blood components and, equally effective for all viral families. In addition, among available inactivation procedures of viral genomes, some of them are recognized to be less effective on non-enveloped viruses, and inadequate to inactivate higher viral titers in plasma pools or derivatives. Given this, there is the need to implement new methodologies for the discovery of unknown viruses that may affect blood transfusion. Viral metagenomics combined with High Throughput Sequencing appears as a promising approach for the identification and global surveillance of new and/or unexpected viruses that could impair blood transfusion safety.

Brain Meta-Transcriptomics from Harbor Seals to Infer the Role of the Microbiome and Virome in a Stranding Event

Marine diseases are becoming more frequent, and tools for identifying pathogens and disease reservoirs are needed to help prevent and mitigate epizootics. Meta-transcriptomics provides insights into disease etiology by cataloguing and comparing sequences from suspected pathogens. This method is a powerful approach to simultaneously evaluate both the viral and bacterial communities, but few studies have applied this technique in marine systems. In 2009 seven harbor seals, Phoca vitulina, stranded along the California coast from a similar brain disease of unknown cause of death (UCD). We evaluated the differences between the virome and microbiome of UCDs and harbor seals with known causes of death. Here we determined that UCD stranded animals had no viruses in their brain tissue. However, in the bacterial community, we identified Burkholderia and Coxiella burnetii as important pathogens associated with this stranding event. Burkholderia were 100% prevalent and ~2.8 log2 fold more abundant in the UCD animals. Further, while C. burnetii was found in only 35.7% of all samples, it was highly abundant (~94% of the total microbial community) in a single individual. In this harbor seal, C. burnetii showed high transcription rates of invading and translation genes, implicating it in the pathogenesis of this animal. Based on these data we propose that Burkholderia taxa and C. burnetii are potentially important opportunistic neurotropic pathogens in UCD stranded harbor seals.

Modular approach to customise sample preparation procedures for viral metagenomics: a reproducible protocol for virome analysis

A major limitation for better understanding the role of the human gut virome in health and disease is the lack of validated methods that allow high throughput virome analysis. To overcome this, we evaluated the quantitative effect of homogenisation, centrifugation, filtration, chloroform treatment and random amplification on a mock-virome (containing nine highly diverse viruses) and a bacterial mock-community (containing four faecal bacterial species) using quantitative PCR and next-generation sequencing. This resulted in an optimised protocol that was able to recover all viruses present in the mock-virome and strongly alters the ratio of viral versus bacterial and 16S rRNA genetic material in favour of viruses (from 43.2% to 96.7% viral reads and from 47.6% to 0.19% bacterial reads). Furthermore, our study indicated that most of the currently used virome protocols, using small filter pores and/or stringent centrifugation conditions may have largely overlooked large viruses present in viromes. We propose NetoVIR (Novel enrichment technique of VIRomes), which allows for a fast, reproducible and high throughput sample preparation for viral metagenomics studies, introducing minimal bias. This procedure is optimised mainly for faecal samples, but with appropriate concentration steps can also be used for other sample types with lower initial viral loads.