Abstract:16S rRNA is the universal gene of microbes, and it is often used as a target gene to obtain profiles of microbial communities via next-generation sequencing (NGS) technology. Traditionally, sequences are clustered into operational taxonomic units (OTUs) at a 97% threshold based on the taxonomic standard using 16S rRNA, and methods for the reduction of sequencing errors are bypassed, which may lead to false classification units. Several denoising algorithms have been published to solve this problem, such as DADA2 and Deblur, which can correct sequencing errors at single-nucleotide resolution by generating amplicon sequence variants (ASVs). As high-resolution ASVs are becoming more popular than OTUs and only one analysis method is usually selected in a particular study, there is a need for a thorough comparison of OTU clustering and denoising pipelines. In this study, three of the most widely used 16S rRNA methods (two denoising algorithms, DADA2 and Deblur, along with de novo OTU clustering) were thoroughly compared using 16S rRNA amplification sequencing data generated from 358 clinical stool samples from the Colorectal Cancer (CRC) Screening Cohort. Our findings indicated that all approaches led to similar taxonomic profiles (with P > 0.05 in PERMNAOVA and P <0.001 in the Mantel test), although the number of ASVs/OTUs and the alpha-diversity indices varied considerably. Despite considerable differences in disease-related markers identified, disease-related analysis showed that all methods could result in similar conclusions. Fusobacterium, Streptococcus, Peptostreptococcus, Parvimonas, Gemella , and Haemophilus were identified by all three methods as enriched in the CRC group, while Roseburia, Faecalibacterium, Butyricicoccus , and Blautia were identified by all three methods as enriched in the healthy group. In addition, disease-diagnostic models generated using machine learning algorithms based on the data from these different methods all achieved good diagnostic efficiency (AUC: 0.87–0.89), with the model based on DADA2 producing the highest AUC (0.8944 and 0.8907 in the training set and test set, respectively). However, there was no significant difference in performance between the models ( P >0.05). In conclusion, this study demonstrates that DADA2, Deblur, and de novo OTU clustering display similar power levels in taxa assignment and can produce similar conclusions in the case of the CRC cohort.

Clustering 16S Rrna for OTU Prediction: a Method of Unsupervised Bayesian Clustering.

Two-Stage Clustering (Tsc): A Pipeline For Selecting Operational Taxonomic Units For The High-Throughput Sequencing Of Pcr Amplicons

Subsampled Open-Reference Clustering Creates Consistent, Comprehensive OTU Definitions and Scales to Billions of Sequences

Comparison of Methods for Picking the Operational Taxonomic Units From Amplicon Sequences

Interpreting 16S metagenomic data without clustering to achieve sub-OTU resolution

M-pick, a modularity-based method for OTU picking of 16S rRNA sequences

ASV vs OTUs clustering: Effects on alpha, beta, and gamma diversities in microbiome metabarcoding studies

Clustering on Human Microbiome Sequencing Data: A Distance-Based Unsupervised Learning Model

Defining Loci in Restriction-Based Reduced Representation Genomic Data from Nonmodel Species: Sources of Bias and Diagnostics for Optimal Clustering

Updating the 97% identity threshold for 16S ribosomal RNA OTUs

TreeOTU: Operational Taxonomic Unit Classification Based on Phylogenetic Trees

Using Mothur to Determine Bacterial Community Composition and Structure in 16S Ribosomal RNA Datasets

TACO: Taxonomic Prediction of Unknown OTUs Through OTU Co-Abundance Networks.

Otu Analysis Using Metagenomic Shotgun Sequencing Data

Nonparametric clustering of RNA-sequencing data

Bayesian Nonparametric Clustering with Feature Selection for Spatially Resolved Transcriptomics Data

Clustering pipeline for determining consensus sequences in targeted next-generation sequencing

Hierarchical Clustering Supported by Reciprocal Nearest Neighbors

A New Unsupervised Binning Approachfor Metagenomic Sequences Based onN-grams and Automatic Feature Weighting.

Omada: robust clustering of transcriptomes through multiple testing

An independent evaluation in a CRC patient cohort of microbiome 16S rRNA sequence analysis methods: OTU clustering, DADA2, and Deblur