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 high-throughput DNA- or RNA-based 16S rRNA (gene) amplicon sequencing in combination with bioinformatics analyses. However, focusing on environmental samples from contrasting habitats, it was not systematically evaluated (i) which analysis methods provide results that reflect reality most accurately, (ii) how the interpretations of microbial community studies are biased by different analysis methods and (iii) if the most optimal analysis workflow can be implemented in an easy-to-use pipeline. Here, we compared the performance of 16S rRNA (gene) amplicon sequencing analysis tools (i.e., Mothur, QIIME1, QIIME2, and MEGAN) using three mock datasets with known microbial community composition that differed in sequencing quality, species number and abundance distribution (i.e., even or uneven), and phylogenetic diversity (i.e., closely related or well-separated amplicon sequences). Our results showed that QIIME2 outcompeted all other investigated tools in sequence recovery (>10 times fewer false positives), taxonomic assignments (>22% better F-score) and diversity estimates (>5% better assessment), suggesting that this approach is able to reflect the in situ microbial community most accurately. Further analysis of 24 environmental datasets obtained from four contrasting terrestrial and freshwater sites revealed dramatic differences in the resulting microbial community composition for all pipelines at genus level. For instance, at the investigated river water sites Sphaerotilus was only reported when using QIIME1 (8% abundance) and Agitococcus with QIIME1 or QIIME2 (2 or 3% abundance, respectively), but both genera remained undetected when analyzed with Mothur or MEGAN. Since these abundant taxa probably have implications for important biogeochemical cycles (e.g., nitrate and sulfate reduction) at these sites, their detection and semi-quantitative enumeration is crucial for valid interpretations. A high-performance computing conformant workflow was constructed to allow FAIR (Findable, Accessible, Interoperable, and Re-usable) 16S rRNA (gene) amplicon sequence analysis starting from raw sequence files, using the most optimal methods identified in our study. Our presented workflow should be considered for future studies, thereby facilitating the analysis of high-throughput 16S rRNA (gene) sequencing data substantially, while maximizing reliability and confidence in microbial community data analysis.

Comparative analysis of size-fractional eukaryotic microbes in subtropical riverine systems inferred from 18S rRNA gene V4 and V9 regions

Comparison of Microbial Diversity Determined with the Same Variable Tag Sequence Extracted from Two Different PCR Amplicons

Unveiling the overlooked small-sized microbiome in river ecosystems

Comparison of Deep-Sea Picoeukaryotic Composition Estimated from the V4 and V9 Regions of 18S Rrna Gene with a Focus on the Hadal Zone of the Mariana Trench

A comparison of two gene regions for assessing community composition of eukaryotic marine microalgae from coastal ecosystems

A multi-amplicon 16S rRNA sequencing and analysis method for improved taxonomic profiling of bacterial communities

Determining Microeukaryotic Plankton Community Around Xiamen Island, Southeast China, Using Illumina Miseq And Pcr-Dgge Techniques

Use of high throughput sequencing and light microscopy show contrasting results in a study of phytoplankton occurrence in a freshwater environment.

Microbial Eukaryotic Diversity and Distribution in a River Plume and Cyclonic Eddy-Influenced Ecosystem in the South China Sea.

Genetic diversity patterns of microbial communities in a subtropical riverine ecosystem (Jiulong River, southeast China)

Patterns in the Composition of Microbial Communities from a Subtropical River: Effects of Environmental, Spatial and Temporal Factors

Comparison of the Levels of Bacterial Diversity in Freshwater, Intertidal Wetland, and Marine Sediments by Using Millions of Illumina Tags

18S rRNA gene sequencing reveals significant influence of anthropogenic effects on microeukaryote diversity and composition along a river-to-estuary gradient ecosystem

Interpretations of Environmental Microbial Community Studies Are Biased by the Selected 16S rRNA (Gene) Amplicon Sequencing Pipeline

[Strategy of Selecting 16S Rrna Hypervariable Regions for Metagenome-Phylogenetic Marker Genes Based Analysis].

Taxonomic Precision of Different Hypervariable Regions of 16S Rrna Gene and Annotation Methods for Functional Bacterial Groups in Biological Wastewater Treatment.

Size-fractionated microbiome observed during an eight-month long sampling in Jiaozhou Bay and the Yellow Sea

Sequencing Insights into Microbial Communities in the Water and Sediments of Fenghe River, China.

Microeukaryotic Community Shifting Along a Lentic-Lotic Continuum

Evaluation of Different Primers of the 18S Rrna Gene to Profile Amoeba Communities in Environmental Samples

High-resolution Single-Molecule Long-Fragment Rrna Gene Amplicon Sequencing of Bacterial and Eukaryotic Microbial Communities