Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • RNA-Seq: TopHat with Roche/454 reads

    Hellow

    I have got Roche/454 RNA-Seq reads from human samples.
    It seems that TopHat/Cufflinks is one of the most used pipeline for analyzing RNA-Seq data, however, for short reads only (non-454 reads).
    The program Bowtie, which is called by TopHat, cannot map 454 reads.

    Still I would like to take adavantage of TopHat's features if possible.
    Has anybody tried to use 454 reads with TopHat, implemented a workaround, or are there any alternatives to TopHat/Bowtie that yield good results for 454 reads?

    Any help, comments are appreciated!

  • #2
    Here is an item from the release notes for TopHat v 1.2.0, just released on 1/18/2011

    Version 1.2.0 also addresses the following issues:

    * A problem with alignment of reads of mixed length has been fixed. Users with a wide range of read lengths, such as trimmed 454 data sets, are strongly encouraged to upgrade.
    So it would appear the 454 reads can be used with TopHat.

    Comment


    • #3
      I have used the latest version. But Bowtie is still not able to map to long 454 reads in the first mapping round. Surely, for the IUM (initially unmapped reads) which comprise the majority of the 454 reads, TopHat chops them into smaller parts that can be mapped by bowtied. But I think that it would be better if one could use the full length 454 reads.

      Comment


      • #4
        Is there any reason you don't want to use the 454 software? Newbler is very good.

        Comment


        • #5
          Can newbler detect new splice variants ?

          Comment


          • #6
            A cDNA de-novo assembly gives an output of isotigs (essentially mRNAs) grouped into isogroups (representing genes) by traversing the contig graph (each contig theoretically represents an exon) which gives you a list of the contigs/exons per isogroup/gene and the combinations of them that form the isotigs/mRNAs. Of course 454 errors make this all much less straightforward than it sounds.

            Not sure if there is an easy way to get to new splice variants from here, but I do know that gsMapper allows you to use a fully annotated genome as a reference including gene names, exon positions and known SNPs which may allow for contigs to be mapped to exons.

            Comment


            • #7
              I have met the same problem. I have also had a look at G-Mo.R-Se, a Tophat alternative. Again, the aligner it uses -- SOAP, only works with short reads...

              Still looking for a splicing junction discovery tool which is independent of the aligners, i.e. works with SAM/BAM files directly.

              Comment


              • #8
                It seems that "MapSplice" fits my needs very well. It splits the raw reads into fragments first, so it could deal with both short and long reads.

                Comment


                • #9
                  The long reads of the 454 are very informative for exon junction sites. I would try to find a method that makes use of this information rather than discarding it by splitting up the reads.

                  Comment


                  • #10
                    Originally posted by Jeremy View Post
                    The long reads of the 454 are very informative for exon junction sites. I would try to find a method that makes use of this information rather than discarding it by splitting up the reads.
                    MapSplice splits the raw reads into fragments first, but the fragments are not processed independently. After each fragment is aligned, the alignments of all fragments of a read will be assembled to form alignments as long as possible

                    Comment


                    • #11
                      Help with Roche454 data and Bowtie!

                      I am working with Roche 454data and am trying to align sequences using Bowtie.
                      $ bowtie -k 2 -a BRCA1_coding -q 1FR.fastq
                      But I am getting No alignments.
                      # reads processed: 1911
                      # reads with at least one reported alignment: 0 (0.00%)
                      # reads that failed to align: 1911 (100.00%)
                      No alignments
                      Is there any way to obtain atleast a few alignments?
                      Or does Bowtie not support 454 data alignment??

                      Comment


                      • #12
                        Our mapper (Genomatix) has a local spliced alignment mode that allows to align complete transcripts to the genome without prior split of the reads (it works fine for 454 reads). Attached is a screenshot of assembled (velvet) RNA-Seq reads mapped to a reference genome.

                        By the way there is another thread discussing a similar topic:
                        Discussion of next-gen sequencing related bioinformatics: resources, algorithms, open source efforts, etc
                        Attached Files

                        Comment


                        • #13
                          Roche454

                          Dear all
                          i an beginner in NGS pls help me.
                          i hv roche 454 data in 10 set in all different condition and i hv to do trancriptome analysis . i m clc genomic software so pls help me how to proceed the task

                          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, 03-27-2024, 06:37 PM
                          0 responses
                          12 views
                          0 likes
                          Last Post seqadmin  
                          Started by seqadmin, 03-27-2024, 06:07 PM
                          0 responses
                          11 views
                          0 likes
                          Last Post seqadmin  
                          Started by seqadmin, 03-22-2024, 10:03 AM
                          0 responses
                          52 views
                          0 likes
                          Last Post seqadmin  
                          Started by seqadmin, 03-21-2024, 07:32 AM
                          0 responses
                          68 views
                          0 likes
                          Last Post seqadmin  
                          Working...
                          X