This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Any idea about extracting bacteria related proteins(or protein clusters) from uniref90 database?

Hi guys,

I want to do a blastp to align some microbial DNA sequences(fastq) against uniref90 protein catalog, in order to answer some questions about protein abundance in my data. The thing is that uniref90 fasta file downloaded from UniProt website has amino acid sequences from all kind of organism including prokaryote and eukaryote. But I only want to focus on bacteria. The uniref90 has already provided a taxid for the lowest common taxon, this could be Genus, Family but they do not list all taxonomic ranks to tell that it is a bacteria.

Some is there any about how to do this extraction?

Thanks in advance

uniref uniprot
• 3,274 views • link
updated by mario_de_manner • mario_de_manner • written by boaty boaty

Thanks, I will try this

• 0 views • link
written by boaty boaty

Were you able to find a running solution, @boaty? I am interested in the exact same question

• 0 views • link
written by mario_de_manner • mario_de_manner

1 answer

I don't think this is necessarily a great idea, but here goes. First you download all bacterial sequences from taxonomic divisions (sprot and trembl .dat files, convert them to fasta with esl-reformat from HMMer). Then you can use either cd-hit or MMseqs2 to trim them to 90% redundancy. It should not take too long for 90% but you will need a computer with a healthy amount of memory and disk space.

Yet another way - and also not great - would be to download all bacterial proteins as above, extract their IDs. Next download UniRef90, format it using makeblastdb and extract bacterial proteins with blastdbcmd and those IDs from the previous step. I think the first approach will be faster.

• 0 views • link
written by Mensur Dlakic Mensur Dlakic

Thanks I will try the second one, I am not sure that the in-house 90% identity file will be the same as uniref90

• 0 views • link
written by boaty boaty

Log in to answer this question.

More posts like this