A framework for assessing 16S rRNA marker-gene survey data analysis methods using mixtures.
Open Access
- 13 March 2020
- journal article
- research article
- Published by Springer Science and Business Media LLC in Microbiome
- Vol. 8 (1), 1-18
- https://doi.org/10.1186/s40168-020-00812-1
Abstract
There are a variety of bioinformatic pipelines and downstream analysis methods for analyzing 16S rRNA marker-gene surveys. However, appropriate assessment datasets and metrics are needed as there is limited guidance to decide between available analysis methods. Mixtures of environmental samples are useful for assessing analysis methods as one can evaluate methods based on calculated expected values using unmixed sample measurements and the mixture design. Previous studies have used mixtures of environmental samples to assess other sequencing methods such as RNAseq. But no studies have used mixtures of environmental to assess 16S rRNA sequencing. We developed a framework for assessing 16S rRNA sequencing analysis methods which utilizes a novel two-sample titration mixture dataset and metrics to evaluate qualitative and quantitative characteristics of count tables. Our qualitative assessment evaluates feature presence/absence exploiting features only present in unmixed samples or titrations by testing if random sampling can account for their observed relative abundance. Our quantitative assessment evaluates feature relative and differential abundance by comparing observed and expected values. We demonstrated the framework by evaluating count tables generated with three commonly used bioinformatic pipelines: (i) DADA2 a sequence inference method, (ii) Mothur a de novo clustering method, and (iii) QIIME an open-reference clustering method. The qualitative assessment results indicated that the majority of Mothur and QIIME features only present in unmixed samples or titrations were accounted for by random sampling alone, but this was not the case for DADA2 features. Combined with count table sparsity (proportion of zero-valued cells in a count table), these results indicate DADA2 has a higher false-negative rate whereas Mothur and QIIME have higher false-positive rates. The quantitative assessment results indicated the observed relative abundance and differential abundance values were consistent with expected values for all three pipelines. We developed a novel framework for assessing 16S rRNA marker-gene survey methods and demonstrated the framework by evaluating count tables generated with three bioinformatic pipelines. This framework is a valuable community resource for assessing 16S rRNA marker-gene survey bioinformatic methods and will help scientists identify appropriate analysis methods for their marker-gene surveys.Keywords
This publication has 41 references indexed in Scilit:
- The SILVA ribosomal RNA gene database project: improved data processing and web-based toolsNucleic Acids Research, 2012
- Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studiesNucleic Acids Research, 2012
- Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variationNucleic Acids Research, 2012
- UCHIME improves sensitivity and speed of chimera detectionBioinformatics, 2011
- Search and clustering orders of magnitude faster than BLASTBioinformatics, 2010
- Ironing out the wrinkles in the rare biosphere through improved OTU clusteringEnvironmental Microbiology, 2010
- QIIME allows analysis of high-throughput community sequencing dataNature Methods, 2010
- Introducing mothur: Open-Source, Platform-Independent, Community-Supported Software for Describing and Comparing Microbial CommunitiesApplied and Environmental Microbiology, 2009
- edgeR: a Bioconductor package for differential expression analysis of digital gene expression dataBioinformatics, 2009
- Naïve Bayesian Classifier for Rapid Assignment of rRNA Sequences into the New Bacterial TaxonomyApplied and Environmental Microbiology, 2007