Abstract:Background Assessments of the soil microbiome provide valuable insight to ecosystem function due to the integral role microorganisms play in biogeochemical cycling of carbon and nutrients. For example, treatment effects on nitrogen cycling functional groups are often presented alongside one another to demonstrate how agricultural management practices affect various nitrogen cycling processes. However, the functional groups commonly evaluated in nitrogen cycling microbiome studies range from phylogenetically narrow (e.g., N-fixation, nitrification) to broad [e.g., denitrification, dissimilatory nitrate reduction to ammonium (DNRA)]. The bioinformatics methods used in such studies were developed for 16S rRNA gene sequence data, and how these tools perform across functional genes of different phylogenetic diversity has not been established. For example, an OTU clustering method that can accurately characterize sequences harboring comparatively little diversity may not accurately resolve the diversity within a gene comprised of a large number of clades. This study uses two nitrogen cycling genes, nifH , a gene which segregates into only three distinct clades, and nrfA , a gene which is comprised of at least eighteen clades, to investigate differences which may arise when using heuristic OTU clustering (abundance-based greedy clustering, AGC) vs. true hierarchical OTU clustering (Matthews Correlation Coefficient optimizing algorithm, Opti-MCC). Detection of treatment differences for each gene were evaluated to demonstrate how conclusions drawn from a given dataset may differ depending on clustering method used. Results The heuristic and hierarchical methods performed comparably for the more conserved gene, nifH . The hierarchical method outperformed the heuristic method for the more diverse gene, nrfA ; this included both the ability to detect treatment differences using PERMANOVA, as well as higher resolution in taxonomic classification. The difference in performance between the two methods may be traced to the AGC method’s preferential assignment of sequences to the most abundant OTUs: when analysis was limited to only the largest 100 OTUs, results from the AGC-assembled OTU table more closely resembled those of the Opti-MCC OTU table. Additionally, both AGC and Opti-MCC OTU tables detected comparable treatment differences using the rank-based ANOSIM test. This demonstrates that treatment differences were preserved using both clustering methods but were structured differently within the OTU tables produced using each method. Conclusion For questions which can be answered using tests agnostic to clustering method (e.g., ANOSIM), or for genes of relatively low phylogenetic diversity (e.g., nifH ), most upstream processing methods should lead to similar conclusions from downstream analyses. For studies involving more diverse genes, however, care should be exercised to choose methods that ensure accurate clustering for all genes. This will mitigate the risk of introducing Type II errors by allowing for detection of comparable treatment differences for all genes assessed, rather than disproportionately detecting treatment differences in only low-diversity genes.

SILVA, RDP, Greengenes, NCBI and OTT — how do these taxonomies compare?

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

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

A comparison between Greengenes, SILVA, RDP, and NCBI reference databases in four published microbiota datasets

Comparing bioinformatic pipelines for microbial 16S rRNA amplicon sequencing

To compare the performance of prokaryotic taxonomy classifiers using curated 16S full-length rRNA sequences

Biootu: an Improved Method for Simultaneous Taxonomic Assignments and Operational Taxonomic Units Clustering of 16s Rrna Gene Sequences.

A Comparison of Methods for Clustering 16s Rrna Sequences into Otus

16S-ITGDB: An Integrated Database for Improving Species Classification of Prokaryotic 16S Ribosomal RNA Sequences

The SILVA and "All-species Living Tree Project (LTP)" taxonomic frameworks

Uniting the classification of cultured and uncultured bacteria and archaea using 16S rRNA gene sequences

Comparison of Mothur and QIIME for the Analysis of Rumen Microbiota Composition Based on 16S rRNA Amplicon Sequences

Embracing Ambiguity in the Taxonomic Classification of Microbiome Sequencing Data

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

From reads to operational taxonomic units: an ensemble processing pipeline for MiSeq amplicon sequencing data

Improving Species Level‐taxonomic Assignment from 16S rRNA Sequencing Technologies

Optimising high-throughput sequencing data analysis, from gene database selection to the analysis of compositional data: a case study on tropical soil nematodes

Nitrogen Cycling Microbial Diversity and Operational Taxonomic Unit Clustering: When to Prioritize Accuracy Over Speed

Evaluating the accuracy of amplicon-based microbiome computational pipelines on simulated human gut microbial communities

TreeOTU: Operational Taxonomic Unit Classification Based on Phylogenetic Trees

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