GuΓ©guen et al. benchmarked 38 combinations of sequence-resolution algorithms (DADA2 standalone; LotuS3 with DADA2 or UPARSE), taxonomic classifiers (Kraken2, DECIPHER, RDP, MMseqs2, Lambda, Metaxa2), and databases (SILVA, GreenGenes2, RDP, RefSeq, Metaxa2) across seven 16S hypervariable regions and three mock communities (Zymo 8 genera, ZIEL1 13, ZIEL2 18; 29 unique genera), evaluated by F1, AUPRC, and CLR-based Euclidean distance. Notably, this is the first benchmark to include the 2024 updates of RDP v2.14, GreenGenes2, and LotuS3 .
Strengths: systematic coverage of the full workflow, transparent code and data (GitHub/Zenodo; PRJNA674596), and honest reporting that identification and quantification diverge. The GC2 underperformance and the L2-pseudo-count caveat are candidly discussed. Limitations the authors acknowledge: only three mocks with one species per genus, a single primer pair, incomplete combination space, and the Escherichia-Shigella merge potentially favoring Metaxa2. BGPT adds: external evidence supports region-dependence of results (V1-V2 vs V3-V4 biases for taxa like Akkermansia and Bifidobacterium in human gut profiling ), and complex 227-strain mocks showed substantial reference-miss rates even for top denoisers , so the 29-genus mocks here likely overstate real-world accuracy. GreenGenes2's failure is plausibly an artifact of its V4-heavy 23M-sequence expansion, which the authors infer but do not test directly. What would falsify the headline: an independent rerun on PRJNA674596 or additional mocks where V4 or a GG2 combination beats V1-V3 DADA2+MMseqs2+SILVA.
Default: V1-V3 + DADA2+MMseqs2+SILVA (pooled or per-sample) or LotuS3+DADA2+RDP; DADA2+Metaxa2 only when quantification of genera stable since the 2012 SILVA 111 reference matters. Future work should control 16S copy-number variation for quantification.
Know what changed, what holds up, and what remains uncertain. Every Friday. No ads.