Open Access

Luise Rauer, Amedeo De Tomassi, Christian L. Müller, Claudia Hülpüsch, Claudia Traidl-Hoffmann, Matthias Reiger, Avidan U. Neumann

• Introduction Microbiome amplicon sequencing data are distorted by multiple protocol-dependent biases from bacterial DNA extraction, contamination, sequence errors, and chimeras, hindering clinical microbiome applications. In particular, extraction bias is a major confounder in sequencing-based microbiome analyses, with no correction method available to date. Here, we suggest using mock community controls to computationally correct extraction bias based on bacterial morphological properties. Methods We compared dilution series of 3 cell mock communities with an even or staggered composition. DNA of these mock, and additional skin microbiome samples, was extracted with 8 different extraction protocols (2 buffers, 2 extraction kits, 2 lysis conditions). Extracted DNA was sequenced (V1–V3 16S rRNA gene) together with corresponding DNA mocks. Results Microbiome composition was significantly different between extraction kits and lysis conditions, but not between buffers. IndependentIntroduction Microbiome amplicon sequencing data are distorted by multiple protocol-dependent biases from bacterial DNA extraction, contamination, sequence errors, and chimeras, hindering clinical microbiome applications. In particular, extraction bias is a major confounder in sequencing-based microbiome analyses, with no correction method available to date. Here, we suggest using mock community controls to computationally correct extraction bias based on bacterial morphological properties. Methods We compared dilution series of 3 cell mock communities with an even or staggered composition. DNA of these mock, and additional skin microbiome samples, was extracted with 8 different extraction protocols (2 buffers, 2 extraction kits, 2 lysis conditions). Extracted DNA was sequenced (V1–V3 16S rRNA gene) together with corresponding DNA mocks. Results Microbiome composition was significantly different between extraction kits and lysis conditions, but not between buffers. Independent of the extraction protocol, chimera formation increased with higher input cell numbers. Contaminants originated mostly from buffers, and considerable cross-contamination was observed in low-input samples. Comparing the microbiome composition of the cell mocks to corresponding DNA mocks revealed taxon-specific protocol-dependent extraction bias. Strikingly, this extraction bias per species was predictable by bacterial cell morphology. Morphology-based computational correction of extraction bias significantly improved resulting microbial compositions when applied to different mock samples, even with different taxa. Equivalent correction of the skin samples showed a substantial impact on microbiome compositions. Conclusions Our results indicate that higher DNA density increases chimera formation during PCR amplification. Furthermore, we show that computational correction of extraction bias based on bacterial cell morphology would be feasible using appropriate positive controls, thus constituting an important step toward overcoming protocol biases in microbiome analysis.…