Skip to content

Reference-based, heuristic algorithm for identification of potential contaminants in mothur-formatted OTU tables

License

Notifications You must be signed in to change notification settings

Arkadiy-Garber/Taxonsluice

Repository files navigation

taxonsluice

Heuristic algorithm for removing likely contaminants from microbiome datasets. The currently accepted input to this program is a mothur-formated OTU table. However, we are currently working on making this tool compatibale with QIIME output and amplicon sequence variants (ASVs).

Please see the Wiki for more information on the algorithm, its installation, and a mock community example.

Dependencies:

-Python3

-BLAST

Usage:

This script takes as input an OTU table that is created by mothur (using the make.shared command: https://www.mothur.org/wiki/MiSeq_SOP#OTUs), as well as a tab-delimited file that maps each sample to a corresponding blank (see sample_blank_map.txt for specific format). The latter file is necessary. However, if you only have one blank for all your samples, you can make this file in the same format, but the second column will only have one value. See sample_blank_map-1blank.txt as an example of how to make the blank map if you only have one blank.

This script will run with only these two inputs, and output two files:

-an OTU table with only the OTUs whose abundance in the sample-specific and non-specific-blanks is 10% or less of the abundance in the environmental samples.

-an OTU table containing sequences that were flagged due to presence within blanks

However, if you provide the script with the location of the SILVA database, and a fasta file containing the OTU sequences, an additional output will include an annotated summary of the OTUs that were flagged. This output will include the closest match to the flagged OTUs in SILVA, the study in which those matches originated, and the source of isolation.

The SILVA database can be downloaded from https://www.arb-silva.de/no_cache/download/archive/current/Exports/. The file you should get is "SILVA_132_SSUParc_tax_silva.fasta.gz". As of April 2018, release 132 is the latest release.

sample command (with the representative 16S sequences in FASTA format and SILVA database provided for classification of potentially-rare OTUs)

python3 taxonsluice.py -blank_map mySamplesToBlanks.txt -otu_table myOTUs.txt -seq_file myOTUs.fasta -silva_DB SILVA_128_SSURef_tax_silva.fasta -rare 5 -t 4 -silva_aln 10 -out_folder /path/to/output/directory/

sample command (simple version without 16S or SILVA database provided)

python3 taxonsluice.py -blank_map mySamplesToBlanks.txt -otu_table myOTUs.txt -rare 5 -out_folder /path/to/output/directory/

Citing taxonsluice:

Ramírez, G. A., Garber, A. I., Lecoeuvre, A., Angelo, T. D., Wheat, C. G., Orcutt, B. N. (2019). Ecology of Subseafloor Crustal Biofilms, 10(August), 1–17. https://doi.org/10.3389/fmicb.2019.01983

Please also cite various dependencies used by taxonsluice.

About

Reference-based, heuristic algorithm for identification of potential contaminants in mothur-formatted OTU tables

Topics

Resources

License

Stars

Watchers

Forks

Packages

No packages published

Languages