title: Pathoscope: species identification and strain attribution with unassembled sequencing data.
abstract: Emerging next-generation sequencing technologies have revolutionized the collection of genomic data for applications in bioforensics, biosurveillance, and for use in clinical settings. However, to make the most of these new data, new methodology needs to be developed that can accommodate large volumes of genetic data in a computationally efficient manner. We present a statistical framework to analyze raw next-generation sequence reads from purified or mixed environmental or targeted infected tissue samples for rapid species identification and strain attribution against a robust database of known biological agents. Our method, Pathoscope, capitalizes on a Bayesian statistical framework that accommodates information on sequence quality, mapping quality, and provides posterior probabilities of matches to a known database of target genomes. Importantly, our approach also incorporates the possibility that multiple species can be present in the sample and considers cases when the sample species/strain is not in the reference database. Furthermore, our approach can accurately discriminate between very closely related strains of the same species with very little coverage of the genome and without the need for multiple alignment steps, extensive homology searches, or genome assembly--which are time-consuming and labor-intensive steps. We demonstrate the utility of our approach on genomic data from purified and in silico "environmental" samples from known bacterial agents impacting human health for accuracy assessment and comparison with other approaches.
authors.affiliation: Department of Statistics, Brigham Young University, Provo, Utah 84602, USA;
authors.forename: Owen E
authors.initials: OE
authors.lastname: Francis
authors.forename: Matthew
authors.initials: M
authors.lastname: Bendall
authors.forename: Solaiappan
authors.initials: S
authors.lastname: Manimaran
authors.forename: Changjin
authors.initials: C
authors.lastname: Hong
authors.forename: Nathan L
authors.initials: NL
authors.lastname: Clement
authors.forename: Eduardo
authors.initials: E
authors.lastname: Castro-Nallar
authors.forename: Quinn
authors.initials: Q
authors.lastname: Snell
authors.forename: G Bruce
authors.initials: GB
authors.lastname: Schaalje
authors.forename: Mark J
authors.initials: MJ
authors.lastname: Clement
authors.forename: Keith A
authors.initials: KA
authors.lastname: Crandall
authors.forename: W Evan
authors.initials: WE
authors.lastname: Johnson
mesh_terms: D000465:Algorithms; D001408:Bacillus anthracis / Q000235:genetics; D001419:Bacteria / Q000145:classification* / Q000235:genetics*; D001499:Bayes Theorem; D023482:Bioterrorism; D042726:Burkholderia mallei / Q000235:genetics; D016957:Burkholderia pseudomallei / Q000235:genetics; D003014:Clostridium botulinum / Q000235:genetics; D019295:Computational Biology / Q000379:methods*; D030541:Databases, Genetic; D004926:Escherichia coli / Q000235:genetics; D004927:Escherichia coli Infections / Q000382:microbiology; D005060:Europe; D005604:Francisella tularensis / Q000235:genetics; D016680:Genome, Bacterial; D023281:Genomics; D059014:High-Throughput Nucleotide Sequencing; D006801:Humans; D017422:Sequence Analysis, DNA; D012984:Software; D013045:Species Specificity; D015010:Yersinia pestis / Q000235:genetics