SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
reciprocal blast renesh Bioinformatics 7 02-10-2017 08:33 AM
Reciprocal blast help Shorash Bioinformatics 5 07-21-2014 10:39 PM
detect orthologs using reciprocal BLAST nicedad Bioinformatics 2 06-07-2012 07:37 AM
Reciprocal Blast With Multiple Contigs Per Gene jwhittall Bioinformatics 0 07-21-2010 01:20 PM
Need Help Regarding Reciprocal best hit in BLAST Bharat Bioinformatics 1 03-11-2010 03:04 AM

Reply
 
Thread Tools
Old 10-28-2014, 12:48 PM   #1
James88
Junior Member
 
Location: Leicester

Join Date: Oct 2014
Posts: 8
Default Reciprocal Blast Guidance...

Hello,

I'm working on creating a core set of orthologous genes in ~40 c.diff strains in order to create a binary tree at the end.

I was told to use reciprocal blast searches. I understand this and have both a script to perform them and handle the output.

My issue is... When doing this reciprocal blast in a pairwise fashion how do I decide which strains to compare with each other? For example, should I compare strain1 with all other strains and then strain2 with all others etc. or should I compare each one with a sensible reference genome (for example, that which the assembly of these genomes was originally carried out)?

Obviously the first method would involve an incredible number of runs. The latter would take fewer, but which is more sensible and will give meaningful results?

I'm a masters student by the way so whilst I think I know everything, I possibly don't...

Thanks!
James88 is offline   Reply With Quote
Old 10-29-2014, 12:11 AM   #2
LeightonP
Member
 
Location: Scotland

Join Date: Feb 2011
Posts: 29
Default

Generally, you need to compare all strains with all other strains. You only have reciprocal BLAST hits between two strains A and B if the BLAST has been performed reciprocally: A vs B, and B vs A.

That's 40 x 40 - 40 = 1560 comparisons, for 40 strains.

I can imagine circumstances where, if the only differences between a pair of isolates were SNPs, you could make assumptions and avoid the comparison. I'd typically just run the comparison, anyway.
LeightonP is offline   Reply With Quote
Old 10-29-2014, 08:43 AM   #3
James88
Junior Member
 
Location: Leicester

Join Date: Oct 2014
Posts: 8
Default

Thank you for the reply/advice. I agree... That certainly felt like the most sensible thing to do. I had better get writing a little script!
James88 is offline   Reply With Quote
Old 10-29-2014, 12:19 PM   #4
LeightonP
Member
 
Location: Scotland

Join Date: Feb 2011
Posts: 29
Default

It's a common task, so you might find that someone has already written a tool/script that would suit your purposes, in a language that suits your workflow. There's also a tool in the galaxy toolshed https://toolshed.g2.bx.psu.edu/ ('blast_rbh').
LeightonP is offline   Reply With Quote
Old 11-05-2014, 02:45 AM   #5
James88
Junior Member
 
Location: Leicester

Join Date: Oct 2014
Posts: 8
Default

I had already... acquired one from the Kostas lab but have since been told that I apparently don't need it as we're using OrthAgogue now which allegedly does the reciprocal blast from an input of an all-vs-all blast output of all my protein sequences.

The joys of bioinformatics!

Again, thank you for the help, it's appreciated.
James88 is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 02:34 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO