Published July 19, 2020
| Version Version 3
Dataset
Open
DADA2 formatted 16S rRNA gene sequences for both bacteria & archaea
Contributors
Contact person:
Data curator:
Other:
Researcher:
- 1. Trend Laboratory, Curtin University of Technology
- 2. Commonwealth Scientific and Industrial Research Organisation (CSIRO)
- 3. WA Human Microbiome Collaboration Centre (WAHMCC), Trend Laboratory, Curtin University
Description
These two combined bacterial and archaeal 16S rRNA gene sequence databases were collated from various sources and formatted for the purpose of using the "assignTaxonomy" command within the DADA2 pipeline.
- RefSeq+RDP: This database contains 14676 bacterial & 660 archaea full 16S rRNA gene sequences. It was compiled in 14/05/2018 from predominantly the NCBI RefSeq 16S rrna database (https://www.ncbi.nlm.nih.gov/refseq/targetedloci/16S_process/) and was supplemented with extra sequences from the RDP database (https://rdp.cme.msu.edu/misc/resources.jsp).
- Genome Taxonomy Database (GTDB): The new version of our dada2 formatted GTDB reference sequences now contains 21965 bacteria and 1126 archaea full 16S rRNA gene sequences. If you wonder why there are fewer species with 16S rRNA, that is because some metagenomics assembled genomes (MAGs) lack the 16S gene and thus cannot be extracted. The database was downloaded from https://data.ace.uq.edu.au/public/gtdb/data/releases/release95/ on 19/07/2020. Please read the release notes and file descriptions.
The formatting to DADA2 format of the databases was done using a simple awk bash scripts. The script takes as input a fasta file as provided by the core databases creators and then it outputs a fasta file with all 7 taxonomy ranks separated by ";" as required for DADA2 compatibility. Additionally, we have concatenated the unique sequence ID be it NCBI/RDP or GTDB ID to the species entry. We see this as an important QC step to highlight the issues/confidence associated with short read taxonomy assignment at the more finer rank levels.
Notes
The RefSeq+RDP database was updated due to a quotation mark bug that was wrongly placed in front of some of the species names. A file with all the affected species names has been uploaded to review. This shouldn't affect any assignments but might have caused some issues reading into R.
The GTDB was updated due to a new release with taxonomy changes has been made available. The core GTDB team advises that everyone using the GTDB to convert to the release 95. I have also formatted all the 16S rRNA sequences in the GTDBr95 that have passed QC. If anyone finds a need for them I can share outside of here because I do not want to confuse anyone.
awk script can be provided on request.
Files
Version2AffectedSeqs.txt
Additional details
References
- Parks, D. H., et al. (2018). "A standardized bacterial taxonomy based on genome phylogeny substantially revises the tree of life." Nature Biotechnology.
- Cole, J. R., Q. Wang, J. A. Fish, B. Chai, D. M. McGarrell, Y. Sun, C. T. Brown, A. Porras-Alfaro, C. R. Kuske, and J. M. Tiedje. 2014. Ribosomal Database Project: data and tools for high throughput rRNA analysis Nucl. Acids Res. 42(Database issue):D633-D642; doi: 10.1093/nar/gkt1244 [PMID: 24288368]
- NCBI 16S RefSeq Nucleotide sequence records: https://www.ncbi.nlm.nih.gov/nuccore?term=33175%5BBioProject%5D+OR+33317%5BBioProject%5D