Toward standards in clinical microbiome studies: comparison of DNA extraction methods and bioinformatic pipelines

When studying the microbiome using next generation sequencing, DNA extraction method, sequencing procedures and bioinformatic processing are crucial to obtain reliable data. Method choice has been demonstrated to strongly affect the final biological interpretation. We assessed the performance of three DNA extraction methods and two bioinformatic pipelines for bacterial microbiota profiling through 16S rRNA gene sequencing, using positive and negative controls for DNA extraction and sequencing and eight different types of high- or low-biomass samples. Performance was evaluated based on quality control passing, DNA yield, richness, diversity and compositional profiles. All DNA extraction methods retrieved the theoretical relative bacterial abundance with maximum three-fold change, and library preparation and sequencing induced little variation. Bioinformatic pipelines showed different results for estimating richness, but diversity and compositional profiles were comparable. DNA extraction methods were successful for feces and oral swabs and variation induced by DNA extraction methods was lower than inter-subject variation. For low-biomass samples, a mixture of genera present in negative controls and sample-specific genera, possibly representing biological signal, were observed. We conclude that the tested bioinformatic pipelines perform equally with pipeline-specific advantages and disadvantages. Two out of three extraction methods performed equally well, while one method was less accurate. Lastly, we demonstrate the importance of including negative controls when analyzing low bacterial biomass samples.