Optimized SMRT-UMI protocol produces highly accurate sequence datasets from diverse populations - application to HIV-1 quasispecies

Abstract

Pathogen diversity resulting in quasispecies can enable persistence and adaptation to host defenses and therapies. However, accurate quasispecies characterization can be impeded by errors introduced during sample handling and sequencing which can require extensive optimizations to overcome. We present complete laboratory and bioinformatics workflows to overcome many of these hurdles. The Pacific Biosciences single molecule real-time platform was used to sequence PCR amplicons derived from cDNA templates tagged with universal molecular identifiers (SMRT-UMI). Optimized laboratory protocols were developed through extensive testing of different sample preparation conditions to minimize between-template recombination during PCR and the use of UMI allowed accurate template quantitation as well as removal of point mutations introduced during PCR and sequencing to produce a highly accurate consensus sequence from each template. Handling of the large datasets produced from SMRT-UMI sequencing was facilitated by a novel bioinformatic pipeline, Probabilistic Offspring Resolver for Primer IDs (PORPIDpipeline), that automatically filters and parses reads by sample, identifies and discards reads with UMIs likely created from PCR and sequencing errors, generates consensus sequences, checks for contamination within the dataset, and removes any sequence with evidence of PCR recombination or early cycle PCR errors, resulting in highly accurate sequence datasets. The optimized SMRT-UMI sequencing method presented here represents a highly adaptable and established starting point for accurate sequencing of diverse pathogens. These methods are illustrated through characterization of human immunodeficiency virus (HIV) quasispecies. Author SummaryThere is a great need to understand the genetic diversity of pathogens in an accurate and timely manner, but many errors can be introduced during the sample handling and sequencing steps which may prevent accurate analyses. In some cases, the errors introduced during these steps can be indistinguishable from real genetic variation and prevent analyses from identifying true sequence variation present in the pathogen population. There are established methods which can help to prevent these types of errors, but can involve many different steps and variables, all of which must be optimized and tested together to ensure the desired effect. Here we show results from testing different methods on a set of HIV+ blood plasma samples and arrive at a streamlined laboratory protocol and bioinformatic pipeline which prevents or corrects for different types of errors that can arise in sequence datasets. These methods should be an accessible starting point for anyone wanting accurate sequencing without extensive optimizations.