Dealing with controls
During DNA extraction, library preparation and sequencing, technical blanks, negative controls and positive controls may be added to detect contamination, confirm successful amplification of target genes, detect false results, and overall, verify that the full workflow (from extraction to bioinformatics) behaves as expected[1]. Field blanks (e.g., empty tubes left opened at the site or extraction of sample storage buffer) may also be included to enable removal of site or sampling material contaminants a posteriori[2].
This section documents how we analyse control samples and translate those results into practical actions in the data set. Using the R package decontam, the workflow enables identification and removal of contaminant DNA sequences based on patterns observed in blanks and negative controls[3]. In addition, positive controls can be used to evaluate whether the observed community composition matches the expected taxonomic profile[2].
The two downloadable scripts below (003_positive_control.R and 003_negative_control.R) implement a standard approach for dealing with negative and positive controls in a phyloseq object.
For a step-by-step explanation of how each control type is handled, see the sections below: