This is a test version of Biostars. For the public version, visit https://www.biostars.org.
By What Criterion Can We Screen The Blast Results?

I used to select the first several sequences as the ideal sequences I want to get. But there are cases that a lot of sequences with high scores and low evalues will retrieve after blasting. When it happens I was suggested to pick up some sequences randomly and exclude the sequences from same organism. Is it ok? Or Does any other better solutions exist?

blast
• 2,964 views • link
updated by Jan Kosinski Jan Kosinski • written by Zhizhong Zhizhong

1 answer

It depends on biological question you want to answer with your BLAST search, and probably whether you deal with DNA or protein sequences.

I can tell you what I do to create multiple sequence alignments of a query protein family based on BLAST result.

  1. retrieve all BLAST hits up to very bad E-value
  2. cluster hits using CLANS IMPORTANT: use only regions of hits that are aligned to query, otherwise clustering can get wrong due to other domains in hits
  3. Select those hits which appear to belong to my family based on the clustering.
  4. Retrieve full sequences of hits
  5. Align them using mafft or other alignment building program
  6. For further analysis, like phylogenetic analysis, use alignments that are filtered versions of the above. Depending on your problem you need different filtering method.

Blindly excluding sequences from the same organism is dangerous, they may be paralogs, but again it depends what you want to answer.

• 0 views • link
written by Jan Kosinski Jan Kosinski

Thanks for your useful answer. when blasting some short sequences such as 16s RNA of bacteriums, the results can be that all the sequence retrieved with indentical E-value and score. If I want to build phylogenetic trees, what should I do with the sequences?

• 0 views • link
written by Zhizhong Zhizhong

I would analyze the whole set manually, for example: 1) align ALL sequences, 2) group by sequence similarity, most of sequence alignment programs do that by default 3) for every group of sequences referring to the same GeneID take one, with a length most close to average length of seq in the alignment, and without regions clearly dissimilar to other seqs in the alignment 4) for every group of nearly identical sequences from the same species, take one only from representative strain (e.g. from all E.coli strains take only those from K12 strain).

• 0 views • link
written by Jan Kosinski Jan Kosinski

Log in to answer this question.

More posts like this