Question: Selecting Msa For Phylogenetic Tree
0
gravatar for Pappu
6.4 years ago by
Pappu1.9k
Pappu1.9k wrote:

I searched the uniprot database by hmmsearch using a custom alignment built from selected sequence. Now I am wondering if it will be a good idea to realign the hits using Muscle or just use the hmmsearch results for building phylogenetic tree. Also should I remove the columns with >33% gap in the MSA before building the phylogenetic tree?

phylogenetics • 2.6k views
ADD COMMENTlink modified 6.4 years ago • written 6.4 years ago by Pappu1.9k
2
gravatar for Pawel Szczesny
6.4 years ago by
Pawel Szczesny3.2k
Poland
Pawel Szczesny3.2k wrote:

What you have from your search is a multiple pairwise alignment (sequences aligned to a reference, one-vs-all relationship). Phylogenetic trees basically by definition require multiple sequence alignment (they represent all-vs-all relationship), so yes, you need to re-align the hits using MUSCLE, MAFFT or whatever else is there for calculation of MSA.

ADD COMMENTlink written 6.4 years ago by Pawel Szczesny3.2k

Thanks for your prompt response. Should I realign the MSA obtained from hmmsearch or get the sequences present in the MSA and align them using Muscle/MAFFT?

ADD REPLYlink written 6.4 years ago by Pappu1.9k
1

It all depends on the question you are asking (and whether your sequence was a fragment, single-domain protein, multi-domain protein, etc.). In the most simple case, if you just used hmmsearch to find remote homologs of your query protein, take out full sequences of the hits and realign them with Muscle/MAFFT. And _again_: what you got from hmmsearch is _not_ MSA.

ADD REPLYlink written 6.4 years ago by Pawel Szczesny3.2k

Thanks. I meant that I can make alignment from hmmalign and it is also recommended in FastTree website.

ADD REPLYlink written 6.4 years ago by Pappu1.9k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1766 users visited in the last hour