Microbial community resemblance methods differ in their ability to detect biologically relevant patterns

Abstract

A comparison of ordination analysis methods used to reveal gradients or clusters in sequence data from microbial communities reveals which methods are best suited for which community structure at different sequencing depth. High-throughput sequencing methods enable characterization of microbial communities in a wide range of environments on an unprecedented scale. However, insight into microbial community composition is limited by our ability to detect patterns in this flood of sequences. Here we compare the performance of 51 analysis techniques using real and simulated bacterial 16S rRNA pyrosequencing datasets containing either clustered samples or samples arrayed across environmental gradients. We found that many diversity patterns were evident with severely undersampled communities and that methods varied widely in their ability to detect gradients and clusters. Chi-squared distances and Pearson correlation distances performed especially well for detecting gradients, whereas Gower and Canberra distances performed especially well for detecting clusters. These results also provide a basis for understanding tradeoffs between number of samples and depth of coverage, tradeoffs that are important to consider when designing studies to characterize microbial communities.