Flexible taxonomic assignment of ambiguous sequencing reads
Open Access
- 7 January 2011
- journal article
- Published by Springer Science and Business Media LLC in BMC Bioinformatics
- Vol. 12 (1), 8
- https://doi.org/10.1186/1471-2105-12-8
Abstract
To characterize the diversity of bacterial populations in metagenomic studies, sequencing reads need to be accurately assigned to taxonomic units in a given reference taxonomy. Reads that cannot be reliably assigned to a unique leaf in the taxonomy (ambiguous reads) are typically assigned to the lowest common ancestor of the set of species that match it. This introduces a potentially severe error in the estimation of bacteria present in the sample due to false positives, since all species in the subtree rooted at the ancestor are implicitly assigned to the read even though many of them may not match it. We present a method that maps each read to a node in the taxonomy that minimizes a penalty score while balancing the relevance of precision and recall in the assignment through a parameter q. This mapping can be obtained in time linear in the number of matching sequences, because LCA queries to the reference taxonomy take constant time. When applied to six different metagenomic datasets, our algorithm produces different taxonomic distributions depending on whether coverage or precision is maximized. Including information on the quality of the reads reduces the number of unassigned reads but increases the number of ambiguous reads, stressing the relevance of our method. Finally, two measures of performance are described and results with a set of artificially generated datasets are discussed. The assignment strategy of sequencing reads introduced in this paper is a versatile and a quick method to study bacterial communities. The bacterial composition of the analyzed samples can vary significantly depending on how ambiguous reads are assigned depending on the value of the q parameter. Validation of our results in an artificial dataset confirm that a combination of values of q produces the most accurate results.Keywords
This publication has 57 references indexed in Scilit:
- QIIME allows analysis of high-throughput community sequencing dataNature Methods, 2010
- Forensic identification using skin bacterial communitiesProceedings of the National Academy of Sciences of the United States of America, 2010
- A core gut microbiome in obese and lean twinsNature, 2008
- Next-generation DNA sequencingNature Biotechnology, 2008
- Accurate taxonomy assignments from 16S rRNA sequences produced by highly parallel pyrosequencersNucleic Acids Research, 2008
- A comparison of random sequence reads versus 16S rDNA sequences for estimating the biodiversity of a metagenomic libraryNucleic Acids Research, 2008
- Short pyrosequencing reads suffice for accurate microbial community analysisNucleic Acids Research, 2007
- Endophytic Bacterial Diversity in Rice (Oryza sativa L.) Roots Estimated by 16S rDNA Sequence AnalysisMicrobial Ecology, 2007
- Microbial diversity in the deep sea and the underexplored “rare biosphere”Proceedings of the National Academy of Sciences of the United States of America, 2006
- Database resources of the National Center for Biotechnology InformationNucleic Acids Research, 2000