Annotation using Pfam

Question

[Md Moinuddin]

Hi Community,

I have 17 metagenomically assembled genomes (MAGs). I want to annotate the protein of these genome using Pfam.

I want a step by step guideline for the annotations in Linux bash.

My idea is: I need to download the Pfam database but there are lot of files and I am confused which one should I download. Next I need to use HMM for the annotations.

Thank you so much!

Answer 1

[Mensur Dlakic]

This is not going to be step by step, but hopefully still useful. First download and unpack Pfam HMM database:

wget https://ftp.ebi.ac.uk/pub/databases/Pfam/current_release/Pfam-A.hmm.gz
gunzip Pfam-A.hmm.gz
hmmpress Pfam-A.hmm

Assuming your protein sequences are in proteome.faa, do the search:

hmmscan -o annotations.txt --tblout annotations.tbl -E 0.01 --cpu 4 Pfam-A.hmm proteome.faa

Programs from the HMMer package are described here.

[Md Moinuddin]

Thank you so much