Volume 22, Number 8—August 2016
Use of Unamplified RNA/cDNA–Hybrid Nanopore Sequencing for Rapid Detection and Characterization of RNA Viruses
Portable and reliable molecular epidemiology techniques and field approaches for assessing virus genomes are desired to inform clinical diagnostics and public health operations. Need for such methods has been highlighted by the recent Middle East respiratory syndrome and Ebola virus disease (EVD) epidemics, during which it became necessary to characterize novel viruses and to evaluate genetic drift, transmission chains, and zoonotic introductions.
To determine if nanopore sequencing can be used as an accelerated viral genome sequencing tool, we utilized a rapid cDNA/RNA–hybrid library preparation procedure to sequence cell cultures of Venezuelan equine encephalitis virus vaccine (VEEV) strain TC-83 or Ebola virus (EBOV) isolate Makona-C05 stock IRF0137. To evaluate nanopore sequencing for rapid, field-deployable pathogen characterization, we collected raw read data and statistics for VEEV and EBOV sequence runs on the MinION sequencing device (Oxford Nanopore Technologies, Oxford, UK). To determine the level of identification and accuracy of genome characterization over sequencing runtime, these reads were then mapped to VEEV and EBOV genomes and to reference databases (RefSeq [www.ncbi.nlm.nih.gov/RefSeq/]). From the results of these analyses, we determined that the current and future versions of nanopore sequencing technology can be used to rapidly identify and characterize pathogens. .
This approach for pathogen identification and characterization differs from the previously used methods on the MinION platform. Biased techniques, such as amplicon sequencing, have proven to be effective in complex sample backgrounds in which titers of the target pathogen might be low, but such approaches limit characterization to known pathogens and require additional viral genome amplification (1–4). Unbiased techniques that require viral genome amplification (5) or that have been optimized for bacterial genomes (6,7) require longer sample and library preparation times, but can detect low pathogen titers or create highly accurate genomic data. We sequenced unamplified poly(A)-tailed viral RNA using rapid cDNA library preparation coupled with real-time data analysis to determine its potential application for pathogen genomic characterization.
VEEV has a single-stranded, linear, poly(A)-tailed RNA genome. Thus, poly-dT primers can be used for cDNA production without further genomic RNA manipulation. The workflow to isolate the RNA and prepare it for sequencing (Technical Appendix[PDF - 160 KB - 10 pages]) took ≈3 hours from the initiation of sample processing to data acquisition on MinION (Figure, panel A). The sequencing of VEEV genomic RNA/cDNA hybrids attained in hours by using MinION revealed reads that mapped to the VEEV TC-83 genome within minutes by using the LAST (Computational Biology Research Consortium, Tokyo, Japan) multiple sequence alignment program (Technical Appendix[PDF - 160 KB - 10 pages]; Figure, panel B [2,4,6]). The coverage increased from 15–60 min from the 3′ end of the VEEV genome with reads aligning directionally from the 3′ to 5′ end of the VEEV genome (Figure, panel B). These alignment characteristics are indicative of the poly-dT priming strategy for poly(A)-tailed RNA.
To determine if the reads generated from VEEV TC-83 would align to the correct viral genome within a set of reference sequences, we used the viral genome reference sequences (http://www.ncbi.nlm.nih.gov/genomes/GenomesGroup.cgi?taxid=10239&opt=Virus), plus the VEEV TC-83 genome database (GenBank accession no. L01443). We then used LAST to align nanopore reads against this set of references (Technical Appendix[PDF - 160 KB - 10 pages]). These alignments were used to generate a top hit table and the associated read alignment statistics against each hit. VEEV TC-83 was the top hit based on LAST alignment versus virus RefSeq genome sequences; wild-type VEEV placed second (Table). VEEV TC-83 was also identified as the top hit when the 15- and 60–min sequencing datasets were compared with alphavirus genome sequences (Technical Appendix[PDF - 160 KB - 10 pages]) (Table), demonstrating accuracy and depth achieved in a short time. We also analyzed the VEEV TC-83 dataset using the cloud-based metagenomic detection platforms Pathosphere (8) and One Codex (www.onecodex.com), and found that the sample contained VEEV (Technical Appendix[PDF - 160 KB - 10 pages]).
Molecular epidemiology, including use of viral genomics, played a major role during the 2013–2016 EVD response, informing contact tracing, diagnostic operability, and public health measures (9–11). To determine if EBOV is amenable to the same rapid sequencing methodology that was used for VEEV, unamplified negative-stranded RNA isolated from EBOV in Trizol (Thermofisher Scientific, http://www.thermofisher.com/us/en/home/brands/product-brand/trizol.html) was poly(A)-tailed, a single complementary strand of cDNA synthesized, and RNA/cDNA hybrids sequenced. The EBOV samples sequenced on MinION rapidly provided usable, accurate data, despite less raw data than the VEEV TC-83 dataset (137kbp for EBOV versus 2.4Mbp for VEEV at 60 min). Using 15-and 60-min time points and an identical alignment strategy to VEEV TC-83 above, we detected EBOV as the top hit within the sequencing dataset when compared to all virus RefSeq sequences (Table).
Despite success against the RefSeq database, the lack of depth within the dataset did not enable differentiation between the EBOV isolate sequenced here and the >1,500 EBOV draft genomes sequenced during the 2013–2016 outbreak (10,12,13), which indicates a limitation in this sequencing approach for negative-stranded RNA viruses. The poly(A)-tailing method was chosen because the reverse transcription primer adapters designed by Oxford Nanopore were developed to interact directly with the motor protein necessary for guiding DNA through the nanopores. This method greatly reduced preparation time and eliminated need for adaptor ligation reagents. This approach can be revisited for sequencing negative-strand RNA viruses (2,5). Despite this limitation, the RefSeq alignments and nearest neighbor calls were possible with limited data, demonstrating the potential power of long-read rapid sequencing on nanopore platforms.
The current Middle East respiratory syndrome, EVD, and Zika virus disease outbreaks illustrate the necessity for rapid characterization of pathogens for environmental detection, clinical evaluation, and epidemiologic investigation. To determine whether nanopore sequencing can fill this role in a fieldable platform, we tested an RNA/cDNA–hybrid sequencing approach on VEEV TC-83 (a positive-stranded RNA virus) and EBOV (a negative-stranded RNA virus) prepared from cell-culture supernatants. This method definitively identified VEEV TC-83 and differentiated it from wild-type VEEV in ≈3 hours, including only 15 min of data acquisition on MinION. VEEV TC-83 was also differentiated from other alphavirus genomes, facilitating strain-level identification of TC-83. EBOV was also identified rapidly by this approach, differentiating the virus in the sample analyzed here from available virus reference genomes. However, variant/isolate level characterization was not possible due to limited data generated from the RNA/cDNA–hybrid approach.
The method applied here is greatly accelerated compared to traditional next-generation sequencing library preparation, and was used with reagents and equipment suitable for austere conditions (e.g., little need for cold chain, steps not requiring PCR). This study confirmed the possibility of accurate RNA virus genome characterization from RNA/cDNA hybrids by using limited sample manipulation, albeit from relatively pure samples. If samples derived directly from clinical matrices (e.g., blood, saliva) were used, this method would probably not support the necessary depth to characterize virus genomes unless the pathogen titer within these samples was high. As the depth of sequence data obtained from nanopore sequencing approaches continues to improve (14) and other pore types (such as RNA-specific sequencing pores) are integrated into commercial products, these unamplified techniques can transition from the laboratory to the field for more complex analysis.
Utilization of nanopore sequencing in Western Africa (2,3) has demonstrated potential for its use, and newly developed methods like this RNA/cDNA–hybrid approach can be integrated into fieldable protocols. For the emerging Zika virus, insufficiently high virus titers in clinical samples usually necessitates virus culture before genomic sequencing (15). Genomic Zika virus isolate characterization efforts would greatly benefit from the approaches outlined here, especially regarding materials needed for genomic library preparation and the time reduction for strain-level identification (15). By preparing and sequencing RNA/cDNA hybrids, the sample-to-answer time for RNA sequencing is greatly reduced, providing pathogen identification and characterization rapidly to inform future decision making.
Dr. Kilianski is a National Research Council fellow in the BioSciences Division at Edgewood Chemical Biological Center. His research focuses on biosurveillance, emerging viral pathogens, and the identification and characterization of novel agents that threaten today’s warfighter.
The authors acknowledge the efforts of the Sabeti Laboratory and Broad Institute of MIT and Harvard for providing the consensus sequence for EBOV/Mak-C05
This project was funded through the United States Army and Edgewood Chemical Biological Center 219 Innovative Project initiative to A.K. A.K. is also supported by the US National Academy of Science and the US Defense Threat Reduction Agency as a National Research Council fellow. A.K. is part of the Oxford Nanopore MinION Access Program, and has received free materials for this research. This work was funded in part through the prime contract of Battelle Memorial Institute, Columbus, Ohio, with the US National Institute of Allergy and Infectious Diseases under Contract no. HHSN272200700016I. A subcontractor to Battelle Memorial Institute who performed this work is: J.H.K., an employee of Tunnell Government Services, Inc, Bethesda, MD. Information in this report is cleared for public release and distribution is unlimited.
- Kilianski A, Haas JL, Corriveau EJ, Liem AT, Willis KL, Kadavy DR, Bacterial and viral identification and differentiation by amplicon sequencing on the MinION nanopore sequencer. Gigascience. 2015;4:12.
- Hoenen T, Groseth A, Rosenke K, Fischer RJ, Hoenen A, Judson SD, Nanopore sequencing as a rapidly deployable Ebola outbreak tool. Emerg Infect Dis. 2016;22:331–4; Epub ahead of print.
- Quick J, Loman NJ, Duraffour S, Simpson JT, Severi E, Cowley L, Real-time, portable genome sequencing for Ebola surveillance. Nature.2016;530:228–32.
- Wang J, Moore NE, Deng Y-M, Eccles DA, Hall RJ. MinION nanopore sequencing of an influenza genome. Front Microbiol. 2015;6:766.
- Greninger AL, Naccache SN, Federman S, Yu G, Mbala P, Bres V, Rapid metagenomic identification of viral pathogens in clinical samples by real-time nanopore sequencing analysis. Genome Med. 2015;7:99.
- Quick J, Ashton P, Calus S, Chatt C, Gossain S, Hawker J, Rapid draft sequencing and real-time nanopore sequencing in a hospital outbreak ofSalmonella. Genome Biol. 2015;16:114.
- Ashton PM, Nair S, Dallman T, Rubino S, Rabsch W, Mwaigwisya S, MinION nanopore sequencing identifies the position and structure of a bacterial antibiotic resistance island. Nat Biotechnol. 2014;33:296–300; Epub ahead of print.
- Kilianski A, Carcel P, Yao S, Roth P, Schulte J, Donarum GB, Pathosphere.org: pathogen detection and characterization through a web-based, open source informatics platform. BMC Bioinformatics. 2015;16:416.
- Mate SE, Kugelman JR, Nyenswah TG, Ladner JT, Wiley MR, Cordier-Lassalle T, Molecular evidence of sexual transmission of Ebola virus. N Engl J Med. 2015;373:2448–54.
- Park DJ, Dudas G, Wohl S, Goba A, Whitmer SLM, Andersen KG, Ebola virus epidemiology, transmission, and evolution during seven months in Sierra Leone. Cell. 2015;161:1516–26.
- Kugelman JR, Sanchez-Lockhart M, Andersen KG, Gire S, Park DJ, Sealfon R, Evaluation of the potential impact of Ebola virus genomic drift on the efficacy of sequence-based candidate therapeutics. MBiol. 2015;6.
- Gire SK, Goba A, Andersen KG, Sealfon RSG, Park DJ, Kanneh L, Genomic surveillance elucidates Ebola virus origin and transmission during the 2014 outbreak. Science. 2014;345:1369–72.
- Carroll MW, Matthews DA, Hiscox JA, Elmore MJ, Pollakis G, Rambaut A, Temporal and spatial analysis of the 2014–2015 Ebola virus outbreak in West Africa. Nature. 2015;524:97–101.
- Ip CLC, Loose M, Tyson JR, de Cesare M, Brown BL, Jain M, MinION analysis and reference consortium: phase 1 data release and analysis.F1000Research. 2015;4:1075.
- Faria NR. Azevedo R do S da S, Kraemer MUG, Souza R, Cunha MS, Hill SC, et al. Zika virus in the Americas: early epidemiological and genetic findings. Science. 2016;352:345–9.
Suggested citation for this article: Kilianski A, Roth PA, Liem AT, Hill JM, Willis KL, Rossmaier RD, et al. Use of unamplified RNA/cDNA–hybrid nanopore sequencing for rapid detection and characterization of RNA virus. Emerg Infect Dis. 2016 Aug [date cited]. http://dx.doi.org/10.3201/eid2208.160270