Rapid Sequencing of Multiple RNA Viruses in Their Native Form

The Applications of Nanopore Sequencing Technology in Animal and Human Virus Research

In recent years, an increasing number of viruses have triggered outbreaks that pose a severe threat to both human and animal life, as well as caused substantial economic losses. It is crucial to understand the genomic structure and epidemiology of these viruses to guide effective clinical prevention and treatment strategies. Nanopore sequencing, a third-generation sequencing technology, has been widely used in genomic research since 2014. This technology offers several advantages over traditional methods and next-generation sequencing (NGS), such as the ability to generate ultra-long reads, high efficiency, real-time monitoring and analysis, portability, and the ability to directly sequence RNA or DNA molecules. As a result, it exhibits excellent applicability and flexibility in virus research, including viral detection and surveillance, genome assembly, the discovery of new variants and novel viruses, and the identification of chemical modifications. In this paper, we provide a comprehensive review of the development, principles, advantages, and applications of nanopore sequencing technology in animal and human virus research, aiming to offer fresh perspectives for future studies in this field.

Direct RNA sequencing coupled with adaptive sampling enriches RNAs of interest in the transcriptome

Abundant cellular transcripts occupy most of the sequencing reads in the transcriptome, making it challenging to assay for low-abundant transcripts. Here, we utilize the adaptive sampling function of Oxford Nanopore sequencing to selectively deplete and enrich RNAs of interest without biochemical manipulation before sequencing. Adaptive sampling performed on a pool of in vitro transcribed RNAs resulted in a net increase of 22-30% in the proportion of transcripts of interest in the population. Enriching and depleting different proportions of the Candida albicans transcriptome also resulted in a 11-13.5% increase in the number of reads on target transcripts, with longer and more abundant transcripts being more efficiently depleted. Depleting all currently annotated Candida albicans transcripts did not result in an absolute enrichment of remaining transcripts, although we identified 26 previously unknown transcripts and isoforms, 17 of which are antisense to existing transcripts. Further improvements in the adaptive sampling of RNAs will allow the technology to be widely applied to study RNAs of interest in diverse transcriptomes.

The genome of a bunyavirus cannot be defined at the level of the viral particle but only at the scale of the viral population

Bunyaviruses are enveloped negative or ambisense single-stranded RNA viruses with a genome divided into several segments. The canonical view depicts each viral particle packaging one copy of each genomic segment in one polarity named the viral strand. Several opposing observations revealed nonequal ratios of the segments, uneven number of segments per virion, and even packaging of viral complementary strands. Unfortunately, these observations result from studies often addressing other questions, on distinct viral species, and not using accurate quantitative methods. Hence, what RNA segments and strands are packaged as the genome of any bunyavirus remains largely ambiguous. We addressed this issue by first investigating the virion size distribution and RNA content in populations of the tomato spotted wilt virus (TSWV) using microscopy and tomography. These revealed heterogeneity in viral particle volume and amount of RNA content, with a surprising lack of correlation between the two. Then, the ratios of all genomic segments and strands were established using RNA sequencing and qRT-PCR. Within virions, both plus and minus strands (but no mRNA) are packaged for each of the three L, M, and S segments, in reproducible nonequimolar proportions determined by those in total cell extracts. These results show that virions differ in their genomic content but together build up a highly reproducible genetic composition of the viral population. This resembles the genome formula described for multipartite viruses, with which some species of the order Bunyavirales may share some aspects of the way of life, particularly emerging properties at a supravirion scale.

Evaluation of sequencing and PCR-based methods for the quantification of the viral genome formula

Viruses show great diversity in their genome organization. Multipartite viruses package their genome segments into separate particles, most or all of which are required to initiate infection in the host cell. The benefits of such seemingly inefficient genome organization are not well understood. One hypothesised benefit of multipartition is that it allows for flexible changes in gene expression by altering the frequency of each genome segment in different environments, such as encountering different host species. The ratio of the frequency of segments is termed the genome formula (GF). Thus far, formal studies quantifying the GF have been performed for well-characterised virus-host systems in experimental settings using RT-qPCR. However, to understand GF variation in natural populations or novel virus-host systems, a comparison of several methods for GF estimation including high-throughput sequencing (HTS) based methods is needed. Currently, it is unclear how HTS-methods compare a golden standard, such as RT-qPCR. Here we show a comparison of multiple GF quantification methods (RT-qPCR, RT-digital PCR, Illumina RNAseq and Nanopore direct RNA sequencing) using three host plants (Nicotiana tabacum, Nicotiana benthamiana, and Chenopodium quinoa) infected with cucumber mosaic virus (CMV), a tripartite RNA virus. Our results show that all methods give roughly similar results, though there is a significant method effect on genome formula estimates. While the RT-qPCR and RT-dPCR GF estimates are congruent, the GF estimates from HTS methods deviate from those found with PCR. Our findings emphasize the need to tailor the GF quantification method to the experimental aim, and highlight that it may not be possible to compare HTS and PCR-based methods directly. The difference in results between PCR-based methods and HTS highlights that the choice of quantification technique is not trivial.

Rapid whole genome sequencing methods for RNA viruses

RNA viruses are the etiological agents of many infectious diseases. Since RNA viruses are error-prone during genome replication, rapid, accurate and economical whole RNA viral genome sequence determination is highly demanded. Next-generation sequencing (NGS) techniques perform whole viral genome sequencing due to their high-throughput sequencing capacity. However, the NGS techniques involve a significant burden for sample preparation. Since to generate complete viral genome coverage, genomic nucleic acid enrichment is required by reverse transcription PCR using virus-specific primers or by viral particle concentration. Furthermore, conventional NGS techniques cannot determine the 5' and 3' terminal sequences of the RNA viral genome. Therefore, the terminal sequences are determined one by one using rapid amplification of cDNA ends (RACE). However, since some RNA viruses have segmented genomes, the burden of the determination using RACE is proportional to the number of segments. To date, there is only one study attempting whole genome sequencing of multiple RNA viruses without using above mentioned methods, but the generated sequences' accuracy compared to the reference sequences was up to 97% and did not reach 100% due to the low read depth. Hence, we established novel methods, named PCR-NGS and RCA-NGS, that were optimized for an NGS machine, MinION. These methods do not require nucleic acid amplification with virus-specific PCR primers, physical viral particle enrichment, and RACE. These methods enable whole RNA viral genome sequencing by combining the following techniques: (1) removal of unwanted DNA and RNA other than the RNA viral genome by nuclease treatment; (2) the terminal of viral genome sequence determination by barcoded linkers ligation; (3) amplification of the viral genomic cDNA using ligated linker sequences-specific PCR or an isothermal DNA amplification technique, such as rolling circle amplification (RCA). The established method was evaluated using isolated RNA viruses with single-stranded, double-stranded, positive-stranded, negative-stranded, non-segmented or multi-segmented genomes. As a result, all the viral genome sequences could be determined with 100% accuracy, and these mean read depths were greater than 2,500×, at least using either of the methods. This method should allow for easy and economical determination of accurate RNA viral genomes.

Molecular Characteristics of Subgenomic RNAs and the Cap-Dependent Translational Advantage Relative to Corresponding Genomic RNAs of Tomato spotted wilt virus

Tomato spotted wilt virus (TSWV) causes severe viral diseases on many economically important plants of Solanaceae. During the infection process of TSWV, a series of 3'-truncated subgenomic RNAs (sgRNAs) relative to corresponding genomic RNAs were synthesized, which were responsible for the expression of some viral proteins. However, corresponding genomic RNAs (gRNAs) seem to possess the basic elements for expression of these viral proteins. In this study, molecular characteristics of sgRNAs superior to genomic RNAs in viral protein expression were identified. The 3' ends of sgRNAs do not cover the entire intergenic region (IGR) of TSWV genomic RNAs and contain the remarkable A-rich characteristics. In addition, the 3' terminal nucleotides of sgRNAs are conserved among different TSWV isolates. Based on the eIF4E recruitment assay and subsequent northern blot, it is suggested that the TSWV sgRNA, but not gRNA, is capped in vivo; this is why sgRNA is competent for protein expression relative to gRNA. In addition, the 5' and 3' untranslated region (UTR) of sgRNA-Ns can synergistically enhance cap-dependent translation. This study further enriched the understanding of sgRNAs of ambisense RNA viruses.

Transcriptome dataset of six human pathogen RNA viruses generated by nanopore sequencing

Long-read sequencing (LRS) approaches shed new light on the complexity of viral (Kakuk et al., 2021 [1]; Boldogkői et al., 2019 [2]; Depledge et a., 2019 [3]), bacterial (Yan et al., 2018 [4]) and eukaryotic (Tilgner et al., 2014 [5]) transcriptomes. Emerging RNA viruses are zoonotic (Woolhouse et al., 2016 [6]) and create public health problems, e.g. influenza pandemic caused by H1N1 virus in (Fraser et al., 2009 [7]), as well as the current SARS-CoV-2 pandemic (Kim et al., 2020 [8]). In this study, we carried out nanopore sequencing for generating transcriptomic data valuable for structural and kinetic profiling of six important human pathogen RNA viruses, the H1N1 subtype of Influenza A virus (IVA), the Zika virus (ZIKV), the West Nile virus (WNV), the Crimean-Congo hemorrhagic fever virus (CCHFV), the Coxsackievirus [group B serotype 5 (CVB5)] and the Vesicular stomatitis Indiana virus (VSIV), and the response of host cells upon viral infection. The raw sequencing data were filtered during basecalling and only high quality reads (Qscore ≥ 7) were mapped to the appropriate viral and host genomes. Length distribution of sequencing reads were assessed and statistics of data were plotted by the ReadStat.4 python script. The datasets can be used to profile the transcriptomic landscape of RNA viruses, provide information for novel gene annotations, can serve as resource for studying the virus-host interactions, and for the analysis of RNA base modifications. These datasets can be used to compare the different sequencing techniques, library preparation approaches, bioinformatics pipelines, and to analyze the RNA profiles of viruses with small RNA genomes.

Systematic Comparison of Nanopore and Illumina Sequencing for the Detection of Plant Viruses and Viroids Using Total RNA Sequencing Approach

High-throughput sequencing (HTS) has become an important tool for plant virus detection and discovery. Nanopore sequencing has been rapidly developing in the recent years and offers new possibilities for fast diagnostic applications of HTS. With this in mind, a study was completed, comparing the most established HTS platform (MiSeq benchtop sequencer-Illumina), with the MinION sequencer (Oxford Nanopore Technologies) for the detection of plant viruses and viroids. Method comparisons were performed on five selected samples, containing two viroids, which were sequenced using nanopore technology for the first time and 11 plant viruses with different genome organizations. For all samples, sequencing libraries for the MiSeq were prepared from ribosomal RNA-depleted total RNA (rRNA-depleted totRNA) and for MinION sequencing, direct RNA sequencing of totRNA was used. Moreover, for one of the samples, which contained five different plant viruses and a viroid, three additional variations of sample preparation for MinION sequencing were also used: direct RNA sequencing of rRNA-depleted totRNA, cDNA-PCR sequencing of totRNA, and cDNA-PCR sequencing of rRNA-depleted totRNA. Whilst direct RNA sequencing of total RNA was the quickest of the tested approaches, it was also the least sensitive: using this approach, we failed to detect only one virus that was present in a sample at an extremely low titer. All other MinION sequencing approaches showed improved performance with outcomes similar to Illumina sequencing, with cDNA-PCR sequencing of rRNA-depleted totRNA showing the best performance amongst tested nanopore MinION sequencing approaches. Moreover, when enough sequencing data were generated, high-quality consensus viral genome sequences could be reconstructed from MinION sequencing data, with high identity to the ones generated from Illumina data. The results of this study implicate that, when an appropriate sample and library preparation are selected, nanopore MinION sequencing could be used for the detection of plant viruses and viroids with similar performance as Illumina sequencing. Taken as a balance of practicality and performance, this suggests that MinION sequencing may be an ideal tool for fast and affordable virus diagnostics.

Native RNA or cDNA Sequencing for Transcriptomic Analysis: A Case Study on Saccharomyces cerevisiae

Direct sequencing of single molecules through nanopores allows for accurate quantification and full-length characterization of native RNA or complementary DNA (cDNA) without amplification. Both nanopore-based native RNA and cDNA approaches involve complex transcriptome procedures at a lower cost. However, there are several differences between the two approaches. In this study, we perform matched native RNA sequencing and cDNA sequencing to enable relevant comparisons and evaluation. Using Saccharomyces cerevisiae, a eukaryotic model organism widely used in industrial biotechnology, two different growing conditions are considered for comparison, including the poly-A messenger RNA isolated from yeast cells grown in minimum media under respirofermentative conditions supplemented with glucose (glucose growth conditions) and from cells that had shifted to ethanol as a carbon source (ethanol growth conditions). Library preparation for direct RNA sequencing is shorter than that for direct cDNA sequencing. The sequence characteristics of the two methods were different, such as sequence yields, quality score of reads, read length distribution, and mapped on reference ability of reads. However, differential gene expression analyses derived from the two approaches are comparable. The unique feature of direct RNA sequencing is RNA modification; we found that the RNA modification at the 5' end of a transcript was underestimated due to the 3' bias behavior of the direct RNA sequencing. Our comprehensive evaluation from this work could help researchers make informed choices when selecting an appropriate long-read sequencing method for understanding gene functions, pathways, and detailed functional characterization.

The Microbiome in Bronchial Biopsies from Smokers and Ex-Smokers with Stable COPD - A Metatranscriptomic Approach

Current knowledge about the respiratory microbiome is mainly based on 16S ribosomal RNA gene sequencing. Newer sequencing approaches, such as metatranscriptomics, offer the technical ability to measure the viable microbiome response to environmental conditions such as smoking as well as to explore its functional role by investigating host-microbiome interactions. However, knowledge about its feasibility in respiratory microbiome research, especially in lung biopsies, is still very limited. RNA sequencing was performed in bronchial biopsies from clinically stable smokers (n = 5) and ex-smokers (n = 6) with COPD not using (inhaled) steroids. The Trinity assembler was used to assemble non-human reads in order to allow unbiased taxonomical and microbial transcriptional analyses. Subsequently, host-microbiome interactions were analyzed based on associations with host transcriptomic data. Ultra-low levels of microbial mass (0.009%) were identified in the RNA-seq data. Overall, no differences were identified in microbiome diversity or transcriptional profiles of microbial communities or individual microbes between COPD smokers and ex-smokers in the initial test dataset as well as a larger replication dataset. We identified an upregulated host gene set, related to the simultaneous presence of Bradyrhizobium, Roseomonas, Brevibacterium.spp., which were related to PERK-mediated unfolded protein response (UPR) and expression of the microRNA-155-5p. Our results show that metatranscriptomic profiling in bronchial biopsy samples from stable COPD patients yields ultra-low levels of microbial mass. Further, this study illustrates the potential of using transcriptional profiling of the host and microbiome to gain more insight into their interaction in the airways.

Current and Future Methodology for Quantitation and Site-Specific Mapping the Location of DNA Adducts

Formation of DNA adducts is a key event for a genotoxic mode of action, and their presence is often used as a surrogate for mutation and increased cancer risk. Interest in DNA adducts are twofold: first, to demonstrate exposure, and second, to link DNA adduct location to subsequent mutations or altered gene regulation. Methods have been established to quantitate DNA adducts with high chemical specificity and to visualize the location of DNA adducts, and elegant bio-analytical methods have been devised utilizing enzymes, various chemistries, and molecular biology methods. Traditionally, these highly specific methods cannot be combined, and the results are incomparable. Initially developed for single-molecule DNA sequencing, nanopore-type technologies are expected to enable simultaneous quantitation and location of DNA adducts across the genome. Herein, we briefly summarize the current methodologies for state-of-the-art quantitation of DNA adduct levels and mapping of DNA adducts and describe novel single-molecule DNA sequencing technologies to achieve both measures. Emerging technologies are expected to soon provide a comprehensive picture of the exposome and identify gene regions susceptible to DNA adduct formation.

Direct Sequencing of RNA and RNA Modification Identification Using Nanopore

Direct RNA sequencing (dRNA-seq) simultaneously enables the detection of RNA modifications and characterization of full-length transcripts. In principle, full-length native RNA molecule is translocated through the nanopore by a motor protein while a sensor measures ionic current shifts. Then, the current shifts are interpreted by an algorithm that turn out to RNA sequence. Currently, the standard protocol of dRNA-seq provided by Oxford Nanopore Technologies (ONT) allows to directly ligate and sequence only polyadenylated RNA (poly(A) RNA). Here, we describe a method of dRNA-seq that can be applied for both poly(A) RNA and non-poly(A) tailed-RNA.

Nanopore Whole Transcriptome Analysis and Pathogen Surveillance by a Novel Solid-Phase Catalysis Approach

The requirement of a large input amount (500 ng) for Nanopore direct RNA-seq presents a major challenge for low input transcriptomic analysis and early pathogen surveillance. The high RNA input requirement is attributed to significant sample loss associated with library preparation using solid-phase reversible immobilization (SPRI) beads. A novel solid-phase catalysis strategy for RNA library preparation to circumvent the need for SPRI bead purification to remove enzymes is reported here. This new approach leverages concurrent processing of non-polyadenylated transcripts with immobilized poly(A) polymerase and T4 DNA ligase, followed by directly loading the prepared library onto a flow cell. Whole transcriptome sequencing, using a human pathogen Listeria monocytogenes as a model, demonstrates this new method displays little sample loss, takes much less time, and generates higher sequencing throughput correlated with reduced nanopore fouling compared to the current library preparation for 500 ng input. Consequently, this approach enables Nanopore low-input direct RNA-seq, improving pathogen detection and transcript identification in a microbial community standard with spike-in transcript controls. Besides, as evident in the bioinformatic analysis, the new method provides accurate RNA consensus with high fidelity and identifies higher numbers of expressed genes for both high and low input RNA amounts.

Synthesis of modified nucleotide polymers by the poly(U) polymerase Cid1: application to direct RNA sequencing on nanopores

Understanding transcriptomes requires documenting the structures, modifications, and abundances of RNAs as well as their proximity to other molecules. The methods that make this possible depend critically on enzymes (including mutant derivatives) that act on nucleic acids for capturing and sequencing RNA. We tested two 3' nucleotidyl transferases, Saccharomyces cerevisiae poly(A) polymerase and Schizosaccharomyces pombe Cid1, for the ability to add base and sugar modified rNTPs to free RNA 3' ends, eventually focusing on Cid1. Although unable to polymerize ΨTP or 1meΨTP, Cid1 can use 5meUTP and 4thioUTP. Surprisingly, Cid1 can use inosine triphosphate to add poly(I) to the 3' ends of a wide variety of RNA molecules. Most poly(A) mRNAs efficiently acquire a uniform tract of about 50 inosine residues from Cid1, whereas non-poly(A) RNAs acquire longer, more heterogeneous tails. Here we test these activities for use in direct RNA sequencing on nanopores, and find that Cid1-mediated poly(I)-tailing permits detection and quantification of both mRNAs and non-poly(A) RNAs simultaneously, as well as enabling the analysis of nascent RNAs associated with RNA polymerase II. Poly(I) produces a different current trace than poly(A), enabling recognition of native RNA 3' end sequence lost by in vitro poly(A) addition. Addition of poly(I) by Cid1 offers a broadly useful alternative to poly(A) capture for direct RNA sequencing on nanopores.

Direct RNA nanopore sequencing of Pseudomonas aeruginosa clone C transcriptomes

The transcriptomes of Pseudomonas aeruginosa clone C isolates NN2 and SG17M during the mid-exponential and early stationary phase of planktonic growth were evaluated by direct RNA sequencing on the nanopore platform and compared with established short-read cDNA sequencing on the Illumina platform. Fifty to ninety percent of the sense RNAs turned out to be rRNA molecules followed by similar proportions of mRNA transcripts and non-coding RNAs. Both platforms detected similar proportions of uncharged tRNAs and 29 yet undescribed antisense tRNAs. For example, the rarest arginine codon was paired with the most abundant tRNA^Arg, and the tRNA^Arg gene is missing for the most frequent arginine codon. More than 90% of the antisense RNA molecules were complementary to a coding sequence. The antisense RNAs were evenly distributed in the genomes. Direct RNA sequencing identified more than 4,000 distinct non-overlapping antisense RNAs during exponential and stationary growth. Besides highly expressed small antisense RNAs less than 200 bases in size, a population of longer antisense RNAs was sequenced that covered a broad range of a few hundred to thousands of bases and could be complementary to a contig of several genes. In summary, direct RNA sequencing identified yet undescribed RNA molecules and an unexpected composition of the pools of tRNAs, sense and antisense RNAs. IMPORTANCE Genome-wide gene expression of bacteria is commonly studied by high-throughput sequencing of size-selected cDNA fragment libraries of reverse-transcribed RNA preparations. However, the depletion of ribosomal RNAs, enzymatic reverse transcription and the fragmentation, size selection and amplification during library preparation lead to inevitable losses of information about the initial composition of the RNA pool. We demonstrate that direct RNA sequencing on the nanopore platform can overcome these limitations. Nanopore sequencing of total RNA yielded novel insights into the Pseudomonas aeruginosa transcriptome that - if replicated in other species - will change our view of the bacterial RNA world. The discovery of sense - antisense pairs of tmRNA, tRNAs and mRNAs indicates a further and unknown level of gene regulation in bacteria.

Nanopore sequencing technology, bioinformatics and applications

Rapid advances in nanopore technologies for sequencing single long DNA and RNA molecules have led to substantial improvements in accuracy, read length and throughput. These breakthroughs have required extensive development of experimental and bioinformatics methods to fully exploit nanopore long reads for investigations of genomes, transcriptomes, epigenomes and epitranscriptomes. Nanopore sequencing is being applied in genome assembly, full-length transcript detection and base modification detection and in more specialized areas, such as rapid clinical diagnoses and outbreak surveillance. Many opportunities remain for improving data quality and analytical approaches through the development of new nanopores, base-calling methods and experimental protocols tailored to particular applications.

Metagenomic Studies of Viruses in Weeds and Wild Plants: A Powerful Approach to Characterise Variable Virus Communities

High throughput sequencing (HTS) has revolutionised virus detection and discovery, allowing for the untargeted characterisation of whole viromes. Viral metagenomics studies have demonstrated the ubiquity of virus infection - often in the absence of disease symptoms - and tend to discover many novel viruses, highlighting the small fraction of virus biodiversity described to date. The majority of the studies using high-throughput sequencing to characterise plant viromes have focused on economically important crops, and only a small number of studies have considered weeds and wild plants. Characterising the viromes of wild plants is highly relevant, as these plants can affect disease dynamics in crops, often by acting as viral reservoirs. Moreover, the viruses in unmanaged systems may also have important effects on wild plant populations and communities. Here, we review metagenomic studies on weeds and wild plants to show the benefits and limitations of this approach and identify knowledge gaps. We consider key genomics developments that are likely to benefit the field in the near future. Although only a small number of HTS studies have been performed on weeds and wild plants, these studies have already discovered many novel viruses, demonstrated unexpected trends in virus distributions, and highlighted the potential of metagenomics as an approach.

Nanopore Assay Reveals Cell-Type-Dependent Gene Expression of Vesicular Stomatitis Indiana Virus and Differential Host Cell Response

Vesicular stomatitis Indiana virus (VSIV) of genus Vesiculovirus, species IndianaVesiculovirus (formerly as Vesicular stomatitis virus, VSV) causes a disease in livestock that is very similar to the foot and mouth disease, thereby an outbreak may lead to significant economic loss. Long-read sequencing (LRS) -based approaches already reveal a hidden complexity of the transcriptomes in several viruses. This technique has been utilized for the sequencing of the VSIV genome, but our study is the first for the application of this technique for the profiling of the VSIV transcriptome. Since LRS is able to sequence full-length RNA molecules, it thereby provides more accurate annotation of the transcriptomes than the traditional short-read sequencing methods. The objectives of this study were to assemble the complete transcriptome of using nanopore sequencing, to ascertain cell-type specificity and dynamics of viral gene expression, and to evaluate host gene expression changes induced by the viral infection. We carried out a time-course analysis of VSIV gene expression in human glioblastoma and primate fibroblast cell lines using a nanopore-based LRS approach and applied both amplified and direct cDNA sequencing (as well as cap-selection) for a fraction of samples. Our investigations revealed that, although the VSIV genome is simple, it generates a relatively complex transcriptomic architecture. In this study, we also demonstrated that VSIV transcripts vary in structure and exhibit differential gene expression patterns in the two examined cell types.

Perspectives and Benefits of High-Throughput Long-Read Sequencing in Microbial Ecology

Short-read, high-throughput sequencing (HTS) methods have yielded numerous important insights into microbial ecology and function. Yet, in many instances short-read HTS techniques are suboptimal, for example, by providing insufficient phylogenetic resolution or low integrity of assembled genomes. Single-molecule and synthetic long-read (SLR) HTS methods have successfully ameliorated these limitations. In addition, nanopore sequencing has generated a number of unique analysis opportunities, such as rapid molecular diagnostics and direct RNA sequencing, and both Pacific Biosciences (PacBio) and nanopore sequencing support detection of epigenetic modifications. Although initially suffering from relatively low sequence quality, recent advances have greatly improved the accuracy of long-read sequencing technologies. In spite of great technological progress in recent years, the long-read HTS methods (PacBio and nanopore sequencing) are still relatively costly, require large amounts of high-quality starting material, and commonly need specific solutions in various analysis steps. Despite these challenges, long-read sequencing technologies offer high-quality, cutting-edge alternatives for testing hypotheses about microbiome structure and functioning as well as assembly of eukaryote genomes from complex environmental DNA samples.

Application of Oxford Nanopore Technology to Plant Virus Detection

The adoption of Oxford Nanopore Technologies (ONT) sequencing as a tool in plant virology has been relatively slow despite its promise in more recent years to yield large quantities of long nucleotide sequences in real time without the need for prior amplification. The portability of the MinION and Flongle platforms combined with lowering costs and continued improvements in read accuracy make ONT an attractive method for both low- and high-scale virus diagnostics. Here, we provide a detailed step-by-step protocol using the ONT Flongle platform that we have developed for the routine application on a range of symptomatic post-entry quarantine and domestic surveillance plant samples. The aim of this methods paper is to highlight ONT’s feasibility as a valuable component to the diagnostician’s toolkit and to hopefully stimulate other laboratories towards the eventual goal of integrating high-throughput sequencing technologies as validated plant virus diagnostic methods in their own right.

De novo Sequencing of Novel Mycoviruses From Fusarium sambucinum: An Attempt on Direct RNA Sequencing of Viral dsRNAs

An increasing number of viruses are continuously being found in a wide range of organisms, including fungi. Recent studies have revealed a wide viral diversity in microbes and a potential importance of these viruses in the natural environment. Although virus exploration has been accelerated by short-read, high-throughput sequencing (HTS), and viral de novo sequencing is still challenging because of several biological/molecular features such as micro-diversity and secondary structure of RNA genomes. This study conducted de novo sequencing of multiple double-stranded (ds) RNA (dsRNA) elements that were obtained from fungal viruses infecting two Fusarium sambucinum strains, FA1837 and FA2242, using conventional HTS and long-read direct RNA sequencing (DRS). De novo assembly of the read data from both technologies generated near-entire genomic sequence of the viruses, and the sequence homology search and phylogenetic analysis suggested that these represented novel species of the Hypoviridae, Totiviridae, and Mitoviridae families. However, the DRS-based consensus sequences contained numerous indel errors that differed from the HTS consensus sequences, and these errors hampered accurate open reading frame (ORF) prediction. Although with its present performance, the use of DRS is premature to determine viral genome sequences, the DRS-mediated sequencing shows great potential as a user-friendly platform for a one-shot, whole-genome sequencing of RNA viruses due to its long-reading ability and relative structure-tolerant nature.

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.

Decoding the epitranscriptional landscape from native RNA sequences

Traditional epitranscriptomics relies on capturing a single RNA modification by antibody or chemical treatment, combined with short-read sequencing to identify its transcriptomic location. This approach is labor-intensive and may introduce experimental artifacts. Direct sequencing of native RNA using Oxford Nanopore Technologies (ONT) can allow for directly detecting the RNA base modifications, although these modifications might appear as sequencing errors. The percent Error of Specific Bases (%ESB) was higher for native RNA than unmodified RNA, which enabled the detection of ribonucleotide modification sites. Based on the %ESB differences, we developed a bioinformatic tool, epitranscriptional landscape inferring from glitches of ONT signals (ELIGOS), that is based on various types of synthetic modified RNA and applied to rRNA and mRNA. ELIGOS is able to accurately predict known classes of RNA methylation sites (AUC > 0.93) in rRNAs from Escherichiacoli, yeast, and human cells, using either unmodified in vitro transcription RNA or a background error model, which mimics the systematic error of direct RNA sequencing as the reference. The well-known DRACH/RRACH motif was localized and identified, consistent with previous studies, using differential analysis of ELIGOS to study the impact of RNA m6A methyltransferase by comparing wild type and knockouts in yeast and mouse cells. Lastly, the DRACH motif could also be identified in the mRNA of three human cell lines. The mRNA modification identified by ELIGOS is at the level of individual base resolution. In summary, we have developed a bioinformatic software package to uncover native RNA modifications.

Decoding the epitranscriptional landscape from native RNA sequences

Traditional epitranscriptomics relies on capturing a single RNA modification by antibody or chemical treatment, combined with short-read sequencing to identify its transcriptomic location. This approach is labor-intensive and may introduce experimental artifacts. Direct sequencing of native RNA using Oxford Nanopore Technologies (ONT) can allow for directly detecting the RNA base modifications, although these modifications might appear as sequencing errors. The percent Error of Specific Bases (%ESB) was higher for native RNA than unmodified RNA, which enabled the detection of ribonucleotide modification sites. Based on the %ESB differences, we developed a bioinformatic tool, epitranscriptional landscape inferring from glitches of ONT signals (ELIGOS), that is based on various types of synthetic modified RNA and applied to rRNA and mRNA. ELIGOS is able to accurately predict known classes of RNA methylation sites (AUC > 0.93) in rRNAs from Escherichiacoli, yeast, and human cells, using either unmodified in vitro transcription RNA or a background error model, which mimics the systematic error of direct RNA sequencing as the reference. The well-known DRACH/RRACH motif was localized and identified, consistent with previous studies, using differential analysis of ELIGOS to study the impact of RNA m6A methyltransferase by comparing wild type and knockouts in yeast and mouse cells. Lastly, the DRACH motif could also be identified in the mRNA of three human cell lines. The mRNA modification identified by ELIGOS is at the level of individual base resolution. In summary, we have developed a bioinformatic software package to uncover native RNA modifications.

Amplicon-Based, Next-Generation Sequencing Approaches to Characterize Single Nucleotide Polymorphisms of Orthohantavirus Species

Whole-genome sequencing (WGS) of viruses from patient or environmental samples can provide tremendous insight into the epidemiology, drug resistance or evolution of a virus. However, we face two common hurdles in obtaining robust sequence information; the low copy number of viral genomes in specimens and the error introduced by WGS techniques. To optimize detection and minimize error in WGS of hantaviruses, we tested four amplification approaches and different amplicon pooling methods for library preparation and examined these preparations using two sequencing platforms, Illumina MiSeq and Oxford Nanopore Technologies MinION. First, we tested and optimized primers used for whole segment PCR or one kilobase amplicon amplification for even coverage using RNA isolated from the supernatant of virus-infected cells. Once optimized we assessed two sources of total RNA, virus-infected cells and supernatant from the virus-infected cells, with four variations of primer pooling for amplicons, and six different amplification approaches. We show that 99–100% genome coverage was obtained using a one-step RT-PCR reaction with one forward and reverse primer. Using a two-step RT-PCR with three distinct tiling approaches for the three genomic segments (vRNAs), we optimized primer pooling approaches for PCR amplification to achieve a greater number of aligned reads, average depth of genome, and genome coverage. The single nucleotide polymorphisms identified from MiSeq and MinION sequencing suggested intrinsic mutation frequencies of ~10⁻⁵-10⁻⁷ per genome and 10⁻⁴-10⁻⁵ per genome, respectively. We noted no difference in the coverage or accuracy when comparing WGS results with amplicons amplified from RNA extracted from infected cells or supernatant of these infected cells. Our results show that high-throughput diagnostics requiring the identification of hantavirus species or strains can be performed using MiSeq or MinION using a one-step approach. However, the two-step MiSeq approach outperformed the MinION in coverage depth and accuracy, and hence would be superior for assessment of genomes for epidemiology or evolutionary questions using the methods developed herein.

Characterization of Emerging Swine Viral Diseases through Oxford Nanopore Sequencing Using Senecavirus A as a Model

Emerging viral infectious diseases present a major threat to the global swine industry. Since 2015, Senecavirus A (SVA) has been identified as a cause of vesicular disease in different countries and is considered an emerging disease. Despite the growing concern about SVA, there is a lack of preventive and diagnostic strategies, which is also a problem for all emerging infectious diseases. Using SVA as a model, we demonstrated that Oxford Nanopore MinION sequencing could be used as a robust tool for the investigation and surveillance of emerging viral diseases. Our results identified that MinION sequencing allowed for rapid, unbiased pathogen detection at the species and strain level for clinical cases. SVA whole genome sequences were generated using both direct RNA sequencing and PCR-cDNA sequencing methods, with an optimized consensus accuracy of 94% and 99%, respectively. The advantages of direct RNA sequencing lie in its shorter turnaround time, higher analytical sensitivity and its quantitative relationship between input RNA and output sequencing reads, while PCR-cDNA sequencing excelled at creating highly accurate sequences. This study developed whole genome sequencing methods to facilitate the control of SVA and provide a reference for the timely detection and prevention of other emerging infectious diseases.

Porcine reproductive and respiratory syndrome virus whole-genome sequencing efficacy with field clinical samples using a poly(A)-tail viral genome purification method

The genomic surveillance of porcine reproductive and respiratory syndrome virus (PRRSV) is based on sequencing of the ORF5 gene of the virus, which covers only 4% of the entire viral genome. It is expected that PRRSV whole-genome sequencing (WGS) will improve PRRSV genomic data and allow better understanding of clinical discrepancies observed in the field when using ORF5 sequencing. Our main objective was to implement an efficient method for WGS of PRRSV from clinical samples. The viral genome was purified using a poly(A)-tail viral genome purification method and sequenced using Illumina technology. We tested 149 PRRSV-positive samples: 80 sera, 33 lungs, 33 pools of tissues, 2 oral fluids, and 1 processing fluid (i.e., castration liquid). Overall, WGS of 67.1% of PRRSV-positive cases was successful. The viral load, in particular for tissues, had a major impact on the PRRSV WGS success rate. Serum was the most efficient type of sample to conduct PRRSV WGS poly(A)-tail assays, with a success rate of 76.3%, and this result can be explained by improved sequencing reads dispersion matching throughout the entire viral genome. WGS was unsuccessful for all pools of tissue and lung samples with Cq values > 26.5, whereas it could still be successful with sera at Cq ≤ 34.1. Evaluation of results of highly qualified personnel confirmed that laboratory skills could affect PRRSV WGS efficiency. Oral fluid samples seem very promising and merit further investigation because, with only 2 samples of low viral load (Cq = 28.8, 32.8), PRRSV WGS was successful.

Whole-Genome Sequencing of Human Enteroviruses from Clinical Samples by Nanopore Direct RNA Sequencing

Enteroviruses are small RNA viruses that affect millions of people each year by causing an important burden of disease with a broad spectrum of symptoms. In routine diagnostic laboratories, enteroviruses are identified by PCR-based methods, often combined with partial sequencing for genotyping. In this proof-of-principle study, we assessed direct RNA sequencing (DRS) using nanopore sequencing technology for fast whole-genome sequencing of viruses directly from clinical samples. The approach was complemented by sequencing the corresponding viral cDNA via Illumina MiSeq sequencing. DRS of total RNA extracted from three different enterovirus-positive stool samples produced long RNA fragments, covering between 59% and 99.6% of the most similar reference genome sequences. The identification of the enterovirus sequences in the samples was confirmed by short-read cDNA sequencing. Sequence identity between DRS and Illumina MiSeq enterovirus consensus sequences ranged between 94% and 97%. Here, we show that nanopore DRS can be used to correctly identify enterovirus genotypes from patient stool samples with high viral load and that the approach also provides rich metatranscriptomic information on sample composition for all life domains.

Efficient Method for Molecular Characterization of the 5' and 3' Ends of the Dengue Virus Genome

Dengue is a mosquito-borne disease that is of major importance in public health. Although it has been extensively studied at the molecular level, sequencing of the 5′ and 3′ ends of the untranslated regions (UTR) commonly requires specific approaches for completion and corroboration. The present study aimed to characterize the 5′ and 3′ ends of dengue virus types 1 to 4. The 5′ and 3′ ends of twenty-nine dengue virus isolates from acute infections were amplified through a modified protocol of the rapid amplification cDNA ends approach. For the 5′ end cDNA synthesis, specific anti-sense primers for each serotype were used, followed by polyadenylation of the cDNA using a terminal transferase and subsequent PCR amplification with oligo(dT) and internal specific reverse primer. At the 3′ end of the positive-sense viral RNA, an adenine tail was directly synthetized using an Escherichia coli poly(A) polymerase, allowing subsequent hybridization of the oligo(dT) during cDNA synthesis. The incorporation of the poly(A) tail at the 5′ and 3′ ends of the dengue virus cDNA and RNA, respectively, allowed for successful primer hybridization, PCR amplification and direct sequencing. This approach can be used for completing dengue virus genomes obtained through direct and next-generation sequencing methods.

Actinoporins: From the Structure and Function to the Generation of Biotechnological and Therapeutic Tools

Actinoporins (APs) are a family of pore-forming toxins (PFTs) from sea anemones. These biomolecules exhibit the ability to exist as soluble monomers within an aqueous medium or as constitutively open oligomers in biological membranes. Through their conformational plasticity, actinoporins are considered good candidate molecules to be included for the rational design of molecular tools, such as immunotoxins directed against tumor cells and stochastic biosensors based on nanopores to analyze unique DNA or protein molecules. Additionally, the ability of these proteins to bind to sphingomyelin (SM) facilitates their use for the design of molecular probes to identify SM in the cells. The immunomodulatory activity of actinoporins in liposomal formulations for vaccine development has also been evaluated. In this review, we describe the potential of actinoporins for use in the development of molecular tools that could be used for possible medical and biotechnological applications.

Report of the 2019 NIST-FDA workshop on standards for next generation sequencing detection of viral adventitious agents in biologics and biomanufacturing

Adventitious virus testing assures product safety by demonstrating the absence of viruses that could be unintentionally introduced during the manufacturing process. The capabilities of next-generation sequencing (NGS) for broad virus detection in biologics have been demonstrated by the detection of known and novel viruses that were previously missed using the recommended routine assays for adventitious agent testing. A meeting was co-organized by the National Institute of Standards and Technology and the U.S. Food and Drug Administration on September 18-19, 2019 in Gaithersburg, Maryland, USA, to facilitate standardization of NGS technologies for applications of adventitious virus testing in biologics. The goal was to assess the currently used standards for virus detection by NGS and their public availability, and to identify additional needs for different types of reference materials and standards (natural and synthetic). The meeting focused on the NGS processes from sample preparation through sequencing but did not thoroughly cover bioinformatics, since this was considered to be the topic of a separate meeting.

Current challenges and best-practice protocols for microbiome analysis

Analyzing the microbiome of diverse species and environments using next-generation sequencing techniques has significantly enhanced our understanding on metabolic, physiological and ecological roles of environmental microorganisms. However, the analysis of the microbiome is affected by experimental conditions (e.g. sequencing errors and genomic repeats) and computationally intensive and cumbersome downstream analysis (e.g. quality control, assembly, binning and statistical analyses). Moreover, the introduction of new sequencing technologies and protocols led to a flood of new methodologies, which also have an immediate effect on the results of the analyses. The aim of this work is to review the most important workflows for 16S rRNA sequencing and shotgun and long-read metagenomics, as well as to provide best-practice protocols on experimental design, sample processing, sequencing, assembly, binning, annotation and visualization. To simplify and standardize the computational analysis, we provide a set of best-practice workflows for 16S rRNA and metagenomic sequencing data (available at https://github.com/grimmlab/MicrobiomeBestPracticeReview).

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.