This is a limited proof of concept to search for research data, not a production system.

Search the MIT Libraries

Title: DADA2 formatted 16S rRNA gene sequences for both bacteria & archaea

Type Dataset Ali Alishum (2019): DADA2 formatted 16S rRNA gene sequences for both bacteria & archaea. Zenodo. Dataset. https://zenodo.org/record/3188334

Authors: Ali Alishum (Trend Laboratory, Curtin University of Technology) ; Ali Alishum (Trend Laboratory, Curtin University of Technology) ; Seersholm Frederik (Trend Laboratory, Curtin University of Technology) ; Greenfield Paul (Commonwealth Scientific and Industrial Research Organisation (CSIRO)) ; Christophersen Claus (WA Human Microbiome Collaboration Centre (WAHMCC), Trend Laboratory, Curtin University) ;

Links

Summary

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): our dada2 formatted GTDB reference sequence set contains 20486 bacteria and 1073 archaea full 16S rRNA gene sequences. The database was downloaded from (http://gtdb.ecogenomic.org/downloads) on 20/11/2018.

The formatting to DADA2 format of the databases was done using a locally written python 2.7 script. The script takes as input a taxonomy .txt file and a fasta file as provided by the core databases creators and then these two files are matched according to a unique sequence identifier available in both files. 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.

More information

  • DOI: 10.5281/zenodo.3188334

Subjects

  • DADA2 format, 16S rRNA, Bacterial, Archaeal

Dates

  • Publication date: 2019
  • Issued: January 16, 2019

Notes

Other: 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. Python script can be provided on request. Other: {"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"]}

Rights

Much of the data past this point we don't have good examples of yet. Please share in #rdi slack if you have good examples for anything that appears below. Thanks!

Format

electronic resource

Relateditems

DescriptionItem typeRelationshipUri
IsVersionOfhttps://doi.org/10.5281/zenodo.2541238
IsPartOfhttps://zenodo.org/communities/rmg
IsPartOfhttps://zenodo.org/communities/zenodo