Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Alternative to blast?

    Is it just me or does blast seem increasingly to be out of date and a major bottleneck for RNA-seq applications?

    The most popular aligners currently (e.g. bowtie, blast) trade off speed for low memory use, but now memory is cheap. There are very fast, memory-intensive aligners for some problems (e.g. star, mummer) but I don't yet know of one that can replace blast for basic problems such as annotating transcripts against a protein database. This basic operation takes us sometimes two weeks using blastx on a 24-cpu machine, which isn't really sustainable for RNA-seq processing.

    So my question is, does anyone know of a better aligner for this problem, and does anyone else agree that someone *should* create an aligner that is more adapted to current hardware costs?

  • #2
    Have you used blat? http://genome.ucsc.edu/FAQ/FAQblat.html If you are looking for homologous matches this may be an option (but not against a huge db like genpept but if you are going against a proteome it would be fine).

    You should also specify what DB you are using to blastx against for the 2 week (24 cores?) run. Are you using some kind of parallel method for that search or is it a serial job?

    Comment


    • #3
      Blat is definitely better since it keeps the index in memory but it doesn't have the blastx mode (does it??) and it could still be much faster e.g. by using a suffix tree.

      What we usually do is blast say 100k or 200k transcripts against the Uniprot taxonomic subsets and some smaller databases. The bacterial Uniprot is the biggest and the one that takes the longest. Usually we allocate one CPU to each database target, so we could improve on that by also threading the larger targets using blast's own threading.

      But it doesn't change the fact that blast is built on an indexing strategy which economizes memory more than necessary, with consequent reduction in speed. I would not be surprised if 100x speedup is easy to achieve with very practical memory use.

      Comment


      • #4
        Originally posted by Will Nelson View Post
        Blat is definitely better since it keeps the index in memory but it doesn't have the blastx mode (does it??) and it could still be much faster e.g. by using a suffix tree..........
        Have a look below into the blat options ...

        Code:
        options:
           -t=type     Database type.  Type is one of:
                         dna - DNA sequence
                         prot - protein sequence
                         dnax - DNA sequence translated in six frames to protein
                       The default is dna
           -q=type     Query type.  Type is one of:
                         dna - DNA sequence
                         rna - RNA sequence
                         prot - protein sequence
                         dnax - DNA sequence translated in six frames to protein
                         rnax - DNA sequence translated in three frames to protein

        Comment


        • #5
          True, true...I haven't messed with blast for a while.

          But look: this is a 10+ year-old program which has not been updated in forever. It doesn't thread. If you want to use it in blastx mode, then the proteins have to be the *query*, meaning they are streamed and not indexed, which is extremely inefficient for search a large protein DB. Moreover blat uses a seed index rather than the more efficient suffix tree....again trading off time for memory.

          Blat is more scalable than blast, or it would be if the two problems above were addressed, but it certainly is nowhere near the best one can do, either for standalone usage, or much less as the engine of a large-scale cloud annotation service.

          Next-gen sequencing needs a next-gen alignment solution. One of the groups with serious experience at this needs to step up and build something better.

          Comment


          • #6
            LAST is good for homologous sequences.

            Comment


            • #7
              For a blast alternative, how about usearch?

              However, it sounds to me like OP isn't setting up his blasts properly. 200k queries against some subset of uniprot (or even the whole thing) with 24 cores shouldn't take even one day given sufficient RAM..
              Last edited by rhinoceros; 11-22-2013, 12:55 AM.
              savetherhino.org

              Comment


              • #8
                For BLASTX use pauda

                Comment


                • #9
                  Originally posted by Will Nelson View Post
                  True, true...I haven't messed with blast for a while.

                  But look: this is a 10+ year-old program which has not been updated in forever. It doesn't thread.
                  Wrong on both parts. Blast does thread. And improvements to it are on-going. Just because a program was created 10+ years ago does not make it obsolete.

                  Being current and multi-threaded doesn't neccessarily make Blast the best solution however I agree with 'rhinoceros' -- 200K queries vs uniprot using 24 core should not take too long. I routinely annotate large rnaSeq results via Blast. This gives at least a 'first-pass' level of annotation. What I have given up on is Blast2Go; that program is way too slow for a large number of reads.

                  Comment


                  • #10
                    Pauda might be very interesting for you as suggested by jimmybee

                    Otherwise Gblast was just mentioned here
                    Discussion of next-gen sequencing related bioinformatics: resources, algorithms, open source efforts, etc

                    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
                    10 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, Yesterday, 06:07 PM
                    0 responses
                    9 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 03-22-2024, 10:03 AM
                    0 responses
                    51 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 03-21-2024, 07:32 AM
                    0 responses
                    67 views
                    0 likes
                    Last Post seqadmin  
                    Working...
                    X