Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • 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!

  • #2
    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.

    Comment


    • #3
      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!

      Comment


      • #4
        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').

        Comment


        • #5
          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.

          Comment

          Latest Articles

          Collapse

          • seqadmin
            Strategies for Sequencing Challenging Samples
            by seqadmin


            Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
            03-22-2024, 06:39 AM
          • seqadmin
            Techniques and Challenges in Conservation Genomics
            by seqadmin



            The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

            Avian Conservation
            Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
            03-08-2024, 10:41 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, Yesterday, 06:37 PM
          0 responses
          8 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, Yesterday, 06:07 PM
          0 responses
          8 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 03-22-2024, 10:03 AM
          0 responses
          49 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 03-21-2024, 07:32 AM
          0 responses
          66 views
          0 likes
          Last Post seqadmin  
          Working...
          X