Question: Ideal settings for blasting a protein/cdna against a genome?
0
gravatar for cmdcolin
9 months ago by
cmdcolin770
United States
cmdcolin770 wrote:

Are there any ideal settings for BLAST to align a protein or cdna against a genome? Or is it perhaps even misleading to rely on BLAST this way and it's better to use something like exonerate (protein2genome/est2genome)?

blast exonerate • 528 views
ADD COMMENTlink written 9 months ago by cmdcolin770
1

Ideal settings are context-dependent: are the proteins/cDNA from the same species as the reference genome, or from a somewhat distant species? What is your purpose, finding close homologs or even distant homologs?

edit: anyway, I would use Spaln or Exonerate for proteins, and GMAP or Spaln or Exonerate for cDNA.

ADD REPLYlink modified 9 months ago • written 9 months ago by h.mon15k

It would be same species or fairly closely related (e.g. getting a danio gene, blasting against fish of interest). I don't necessarily have exact use cases but there were biologists wanting to find the location of their gene of interest in a poorly annotated genome, and they had cdna sequences from either experiment or a closely related species, and they wanted to blast against it and then look at it in the genome browser. I just noticed when trying to help that blasting a cDNA against the genome would get lots of similar paralog hits but might not have been obtaining the top matches as the right gene (or sometimes it would not even be in list) using default BLASTN of a cdna to the genome

ADD REPLYlink written 9 months ago by cmdcolin770

If that is the case, maybe the best approach is to iterate over blast algorithms: start with megablast; if nothing found, try dc-megablast; if nothing found, try blastn - of course, this is only feasible if done by hand with a few genes.

But the problem may be with the assembly, it may be too fragmented and the gene is not present or is fragmented between several contigs.

ADD REPLYlink written 9 months ago by h.mon15k
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: 733 users visited in the last hour