Interpretations of microbial community studies are biased by the selected 16S rRNA gene amplicon sequencing pipeline.

Abstract

One of the major methods to identify microbial community composition, to unravel microbial population dynamics, and to explore microbial diversity in environmental samples is DNA- or RNA-based 16S rRNA (gene) amplicon sequencing. Subsequent bioinformatics analyses are required to extract valuable information from the high-throughput sequencing approach. However, manifold bioinformatics tools complicate their choice and might cause differences in data interpretation, making the selection of the pipeline a crucial step. Here, we compared the performance of most widely used 16S rRNA gene amplicon sequencing analysis tools (i.e. Mothur, QIIME1, QIIME2, and MEGAN) using mock datasets and environmental samples from contrasting terrestrial and freshwater sites. Our results showed that QIIME2 outcompeted all other investigated tools in sequence recovery (>10 times less false positives), taxonomic assignments (>22% better F-score) and diversity estimates (>5% better assessment), while there was still room for improvement e.g. imperfect sequence recovery (recall up to 87%) or detection of additional false sequences (precision up to 72%). Furthermore, we found that microbial diversity estimates and highest abundant taxa varied among analysis pipelines (i.e. only one in five genera was shared among all analysis tools) when analyzing environmental samples, which might skew biological conclusions. Our findings were subsequently implemented in a high-performance computing conformant workflow following the FAIR (Findable, Accessible, Interoperable, and Re-usable) principle, allowing reproducible 16S rRNA gene amplicon sequence analysis starting from raw sequence files. Our presented workflow can be utilized for future studies, thereby facilitating the analysis of high-throughput DNA- or RNA-based 16S rRNA (gene) sequencing data substantially. ImportanceMicroorganisms play an essential role in biogeochemical cycling events across the globe. Phylogenetic marker gene analysis is a widely used method to explore microbial community dynamics in space and time, to predict the ecological relevance of microbial populations, or to identify microbial key players in biogeochemical cycles. Several computational analysis methods were developed to aid 16S rRNA gene analysis but choosing the best method is not trivial. In this study, we compared popular analysis methods (i.e. Mothur, QIIME1 and 2, and MEGAN) using samples with known microbial composition (i.e. mock community samples) and environmental samples from contrasting habitats (i.e. groundwater, soil, sediment, and river water). Our findings provide guidance for choosing the currently optimal 16S rRNA gene sequencing analysis method and we implemented our recommended pipeline into a reproducible workflow, which follows highest bioinformatics standards and is open source and free to use.