lunes, 18 de julio de 2016

Use of Unamplified RNA/cDNA–Hybrid Nanopore Sequencing for Rapid Detection and Characterization of RNA Viruses - Volume 22, Number 8—August 2016 - Emerging Infectious Disease journal - CDC

Use of Unamplified RNA/cDNA–Hybrid Nanopore Sequencing for Rapid Detection and Characterization of RNA Viruses - Volume 22, Number 8—August 2016 - Emerging Infectious Disease journal - CDC







Volume 22, Number 8—August 2016

Dispatch

Use of Unamplified RNA/cDNA–Hybrid Nanopore Sequencing for Rapid Detection and Characterization of RNA Viruses

Andy Kilianski, Pierce A. Roth, Alvin T. Liem, Jessica M. Hill, Kristen L. Willis, Rebecca D. Rossmaier, Andrew V. Marinich, Michele N. Maughan, Mark A. Karavis, Jens H. Kuhn, Anna N. Honko, and C. Nicole RosenzweigComments to Author 
Author affiliations: US Army Edgewood Chemical Biological Center, Aberdeen Proving Ground, Maryland, USA (A. Kilianski, P.A. Roth, A.T. Liem, J.M. Hill, K.L. Willis, R.D. Rossmaier, A.V. Marinich, M.N. Maughan, M.A. Karavis, C.N. Rosenzweig)Defense Threat Reduction Agency, Fort Belvoir, Virginia, USA (K.L. Willis)National Institutes of Health, Fort Detrick, Frederick, Maryland, USA (J.H. Kuhn, A.N. Honko)

Abstract

Nanopore sequencing, a novel genomics technology, has potential applications for routine biosurveillance, clinical diagnosis, and outbreak investigation of virus infections. Using rapid sequencing of unamplified RNA/cDNA hybrids, we identified Venezuelan equine encephalitis virus and Ebola virus in 3 hours from sample receipt to data acquisition, demonstrating a fieldable technique for RNA virus characterization.
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. .

The Study

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 (14). 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.
Thumbnail of Use of unamplified RNA/cDNA−-hybrid nanopore sequencing for genomic characterization of Venezuelan equine encephalitis virus (VEEV) TC-83. A) Sample preparation workflow for nanopore sequencing. First, viral RNA from BHK21 cell cultures of VEEV TC-83 was isolated, then single strand complimentary DNA (cDNA) was synthesized. The resulting RNA/cDNA hybrids were then prepared for nanopore sequencing and sequenced with data analysis occurring in real time. B) Genome organization and seq
Figure. Use of unamplified RNA/cDNA−-hybrid nanopore sequencing for genomic characterization of Venezuelan equine encephalitis virus (VEEV) TC-83. A) Sample preparation workflow for nanopore sequencing. First, viral RNA from BHK21 cell cultures of...
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 (911). 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.

Conclusions

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.

Acknowledgments

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.

References

  1. Kilianski AHaas JLCorriveau EJLiem ATWillis KLKadavy DRBacterial and viral identification and differentiation by amplicon sequencing on the MinION nanopore sequencer. Gigascience. 2015;4:12. PubMed
  2. Hoenen TGroseth ARosenke KFischer RJHoenen AJudson SDNanopore sequencing as a rapidly deployable Ebola outbreak tool. Emerg Infect Dis2016;22:3314Epub ahead of printDOIPubMed
  3. Quick JLoman NJDuraffour SSimpson JTSeveri ECowley LReal-time, portable genome sequencing for Ebola surveillance. Nature.2016;530:22832.PubMed
  4. Wang JMoore NEDeng Y-MEccles DAHall RJ. MinION nanopore sequencing of an influenza genome. Front Microbiol. 2015;6:766.
  5. Greninger ALNaccache SNFederman SYu GMbala PBres VRapid metagenomic identification of viral pathogens in clinical samples by real-time nanopore sequencing analysis. Genome Med. 2015;7:99.
  6. Quick JAshton PCalus SChatt CGossain SHawker JRapid draft sequencing and real-time nanopore sequencing in a hospital outbreak ofSalmonella. Genome Biol2015;16:114DOIPubMed
  7. Ashton PMNair SDallman TRubino SRabsch WMwaigwisya SMinION nanopore sequencing identifies the position and structure of a bacterial antibiotic resistance island. Nat Biotechnol2014;33:296300Epub ahead of printDOIPubMed
  8. Kilianski ACarcel PYao SRoth PSchulte JDonarum GBPathosphere.org: pathogen detection and characterization through a web-based, open source informatics platform. BMC Bioinformatics2015;16:416DOIPubMed
  9. Mate SEKugelman JRNyenswah TGLadner JTWiley MRCordier-Lassalle TMolecular evidence of sexual transmission of Ebola virus. N Engl J Med2015;373:244854DOIPubMed
  10. Park DJDudas GWohl SGoba AWhitmer SLMAndersen KGEbola virus epidemiology, transmission, and evolution during seven months in Sierra Leone. Cell2015;161:151626DOIPubMed
  11. Kugelman JRSanchez-Lockhart MAndersen KGGire SPark DJSealfon REvaluation of the potential impact of Ebola virus genomic drift on the efficacy of sequence-based candidate therapeutics. MBiol. 2015;6. PubMed
  12. Gire SKGoba AAndersen KGSealfon RSGPark DJKanneh LGenomic surveillance elucidates Ebola virus origin and transmission during the 2014 outbreak. Science2014;345:136972DOIPubMed
  13. Carroll MWMatthews DAHiscox JAElmore MJPollakis GRambaut ATemporal and spatial analysis of the 2014–2015 Ebola virus outbreak in West Africa. Nature2015;524:97101DOIPubMed
  14. Ip CLCLoose MTyson JRde Cesare MBrown BLJain MMinION analysis and reference consortium: phase 1 data release and analysis.F1000Research. 2015;4:1075. DOIPubMed
  15. 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. PubMed

Figure

Table

Technical Appendix

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


DOI: 10.3201/eid2208.160270

No hay comentarios:

Publicar un comentario