This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Kraken2 database curation

Hi all, I'm working on mouse gut microbiome samples and want to use kraken to get their taxonomic profile. I'm using kraken2 with the databases nt and bacteria (plus some others). The problem is that there are bacterial sequences integrated in some genomes in nt, they can be easy to track like weird mammals (I assume no bat entered the lab) but perhaps some parasites have bacterial DNA or fungi and these might be relevant. My question is, is there a neat way to remove those pseudo-bacterial sequences from the database or do some post analysis to remove these unspecific mappings?

Thanks

metagenomics kraken
• 6,210 views • link
written by Asaf Asaf

1 answer

So, five months later I'm happy to introduce domain_classifier which is a pretty simple naive-Bayes classifier to tell if a sequence is prokaryote or eukaryote. I wrote a civet pipeline, which is a pipeline management system internal to the Jackson Laboratory but also available on github to build the kraken2 database

What this package does is first predict PFAM domains on predicted ORFs and then use these domains to classify into a taxonomic domain. To filter the kraken DB I simply remove DNA sequences that strongly disagree with the reported taxonomy. This also removes mitochondrial and chloroplast genomes.

• 0 views • link
written by Asaf Asaf

Log in to answer this question.

More posts like this