High-Throughput Next-Generation Sequencing of Polioviruses

ABSTRACT The poliovirus (PV) is currently targeted for worldwide eradication and containment. Sanger-based sequencing of the viral protein 1 (VP1) capsid region is currently the standard method for PV surveillance. However, the whole-genome sequence is sometimes needed for higher resolution global surveillance. In this study, we optimized whole-genome sequencing protocols for poliovirus isolates and FTA cards using next-generation sequencing (NGS), aiming for high sequence coverage, efficiency, and throughput. We found that DNase treatment of poliovirus RNA followed by random reverse transcription (RT), amplification, and the use of the Nextera XT DNA library preparation kit produced significantly better results than other preparations. The average viral reads per total reads, a measurement of efficiency, was as high as 84.2% ± 15.6%. PV genomes covering >99 to 100% of the reference length were obtained and validated with Sanger sequencing. A total of 52 PV genomes were generated, multiplexing as many as 64 samples in a single Illumina MiSeq run. This high-throughput, sequence-independent NGS approach facilitated the detection of a diverse range of PVs, especially for those in vaccine-derived polioviruses (VDPV), circulating VDPV, or immunodeficiency-related VDPV. In contrast to results from previous studies on other viruses, our results showed that filtration and nuclease treatment did not discernibly increase the sequencing efficiency of PV isolates. However, DNase treatment after nucleic acid extraction to remove host DNA significantly improved the sequencing results. This NGS method has been successfully implemented to generate PV genomes for molecular epidemiology of the most recent PV isolates. Additionally, the ability to obtain full PV genomes from FTA cards will aid in facilitating global poliovirus surveillance.

[1]  Alison L. Van Eenennaam,et al.  A Metagenomics and Case-Control Study To Identify Viruses Associated with Bovine Respiratory Disease , 2015, Journal of Virology.

[2]  Terry Ng,et al.  An ensemble strategy that significantly improves de novo assembly of microbial genomes from metagenomic next-generation sequencing data , 2015, Nucleic acids research.

[3]  E. Delwart,et al.  Comparing viral metagenomics methods using a highly multiplexed human viral pathogens reagent , 2014, Journal of Virological Methods.

[4]  M. Pallansch,et al.  Detection of vaccine-derived polioviruses in Mexico using environmental surveillance. , 2014, The Journal of infectious diseases.

[5]  S. Lowther,et al.  Environmental surveillance for polioviruses in the Global Polio Eradication Initiative. , 2014, The Journal of infectious diseases.

[6]  T. Ng,et al.  Preservation of viral genomes in 700-y-old caribou feces from a subarctic ice patch , 2014, Proceedings of the National Academy of Sciences.

[7]  Dmitry Antipov,et al.  Assembling Single-Cell Genomes and Mini-Metagenomes From Chimeric MDA Products , 2013, J. Comput. Biol..

[8]  Rachel L. Marine,et al.  High Variety of Known and New RNA and DNA Viruses of Diverse Origins in Untreated Sewage , 2012, Journal of Virology.

[9]  O. Kew.  Reaching the last one per cent: progress and challenges in global polio eradication. , 2012, Current opinion in virology.

[10]  B. Semler,et al.  Viral subversion of host functions for picornavirus translation and RNA replication. , 2012, Future virology.

[11]  M. Roivainen,et al.  Role of environmental poliovirus surveillance in global polio eradication and beyond , 2011, Epidemiology and Infection.

[12]  Marcel Martin.  Cutadapt removes adapter sequences from high-throughput sequencing reads , 2011 .

[13]  K. Chumakov,et al.  Massively parallel sequencing for monitoring genetic consistency and quality control of live viral vaccines , 2010, Proceedings of the National Academy of Sciences.

[14]  N. Nathanson,et al.  From Emergence to Eradication: The Epidemiology of Poliomyelitis Deconstructed , 2010, American journal of epidemiology.

[15]  Ning Ma,et al.  BLAST+: architecture and applications , 2009, BMC Bioinformatics.

[16]  A. Vincent,et al.  Rapid Group-, Serotype-, and Vaccine Strain-Specific Identification of Poliovirus Isolates by Real-Time Reverse Transcription-PCR Using Degenerate Primers and Probes Containing Deoxyinosine Residues , 2009, Journal of Clinical Microbiology.

[17]  S. West,et al.  The Analysis of Count Data: A Gentle Introduction to Poisson Regression and Its Alternatives , 2009, Journal of personality assessment.

[18]  Cole Trapnell,et al.  Ultrafast and memory-efficient alignment of short DNA sequences to the human genome , 2009, Genome Biology.

[19]  Christopher F. Baum.  Stata Tip 63: Modeling Proportions , 2008 .

[20]  O. Kew,et al.  Calibration of Multiple Poliovirus Molecular Clocks Covering an Extended Evolutionary Range , 2008, Journal of Virology.

[21]  A. Vincent,et al.  Modulation of Poliovirus Replicative Fitness in HeLa Cells by Deoptimization of Synonymous Codon Usage in the Capsid Region , 2006, Journal of Virology.

[22]  A. Bruu,et al.  Enteroviruses: Polioviruses, Coxsackieviruses, Echoviruses and Newer Enteroviruses , 2003 .

[23]  S. Emerson,et al.  A virus discovery method incorporating DNase treatment and its application to the identification of two bovine parvovirus species , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[24]  International Human Genome Sequencing Consortium.  Initial sequencing and analysis of the human genome , 2001, Nature.

[25]  D. Wood,et al.  L20B cells simplify culture of polioviruses from clinical samples , 1999, Journal of medical virology.