Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • #16
    Originally posted by Brajbio View Post
    can tophat be used for detection of gene fusion?
    Several people asked about TopHat. Since last summer (2011), we have released TopHat-Fusion, which can detect fusion transcripts from either single reads or paired end reads. It is very fast and highly effective at filtering out the numerous false positives that plague these types of tools.

    Comment


    • #17
      I recommend MapSplice for detecting gene fusions in PE data. I use the full-running, cluster, MapPer, and filter-fusion-by-repeat options. I like how it reports putative gene fusions in an easily annotated bed file.

      Comment


      • #18
        Originally posted by salzberg View Post
        Several people asked about TopHat. Since last summer (2011), we have released TopHat-Fusion, which can detect fusion transcripts from either single reads or paired end reads. It is very fast and highly effective at filtering out the numerous false positives that plague these types of tools.
        Hi,
        I am sure Tophat-fusion is awesome, but it would be great to get some feedback on the problems observed by several users and reported here :

        Discussion of next-gen sequencing related bioinformatics: resources, algorithms, open source efforts, etc


        Thanks,

        Nicolas
        PS: I've sent an email to [email protected] but with no success....

        Comment


        • #19
          Also I forgot to mention that TopHat-Fusion can handle SOLiD reads (as can TopHat).

          Comment


          • #20
            Problems with running DeFuse

            Problem solved
            Last edited by JMath; 03-21-2012, 01:36 AM.

            Comment


            • #21
              This bug has now been fixed with defuse version 0.5.0

              Originally posted by srividya View Post
              Hello,

              I am trying to use defuse. I tried to create the 2bit reference genome file. I got some errors. Can anyone check if there is some mistake with the config file that i changed?

              /Apps/serial/defuse/defuse-0.4.3/scripts/create_reference_dataset.pl -c config.txt

              This is the error that I got when I ran the create_reference_dataset.pl,

              Use of uninitialized value in concatenation (.) or string at /Apps/serial/defuse/defuse-0.4.3/scripts/cmdrunner.pm line 39.

              Thanks,
              Srividya

              Comment


              • #22
                I'm curious about a comparison of these fusion gene discovery algorthms !
                For different types of genome sizes and genuses ? and for different type of data ?

                Anyone is aware about a paper talking about that ?

                Comment


                • #23
                  Originally posted by salzberg View Post
                  Several people asked about TopHat. Since last summer (2011), we have released TopHat-Fusion, which can detect fusion transcripts from either single reads or paired end reads. It is very fast and highly effective at filtering out the numerous false positives that plague these types of tools.
                  Unless you have a master list of all actual fusion events in the cell lines MCF7 and SKBR3, the claim that TopHat-Fusion has "fewer false positives" entails that the known fusions are the ONLY FUSIONS, which is a claim that is both is unsubstantiated and misleading: how do you differentiate between tophat-fusion being over-fitted for known data from defuse being non-specific?

                  Furthermore, why didn't you guys do a PCR for the 42 novel fusions and 42 randomly selected deFuse fusions to at least demonstrate the sensitivity and specificity in discovering new fusions? At the very least, rerun TopHat fusion with the same data as Table 3 in the deFuse paper and make use of their PCR results.

                  How many of the 42 Tophat-Fusion novel fusions agreed with the 1670 deFuse novel fusions?

                  edit: I am in no way affiliated with either group; I just cannot stand lopsided claims in the literature.
                  Last edited by Ichinichi; 04-18-2012, 11:46 AM.

                  Comment


                  • #24
                    Ichinichi, if you think there are really 1670 unknown fusion events in one of these cell lines, by all means spend some time trying to validate them.

                    We don't make "lopsided claims" in our papers but I'm not sure what you're referring to by your statement. We stand by the results in the TopHat-Fusion paper.

                    Comment


                    • #25
                      tophat-fusion deFuse comparison

                      As the deFuse author Ill add my comments into this discussion.

                      For the tophat-fusion paper the authors set the anchor length for deFuse to 13bp. I never intended this parameter to be user modified, and as such the parameter isnt propagated to everywhere in the pipeline, and will have a strange effect on calculating probabilities for the fusions. The result will be more false predictions with probability that exceeds the threshold (default 0.5).

                      I reran the MCF7 and SKBR3 libraries setting anchor length to 4bp, and also setting discord_read_trim for MCF7 to 42 as suggested by a warning message produced by deFuse.

                      MCF7 and SKBR3 produce 66 and 126 fusions with > 0.5 probability respectively.

                      Generally, i think its a strength of deFuse that we report more, including events such as novel 5’ exons, or unexpected intergenic splicing. If people are not interested in these they can filter them. People can also select their own sensitivity / specificity by thresholding on the probability score. I think this is preferable to having your results filtered for you.

                      Finally, I think tophat-fusion is an important contribution, especially since as it says in the paper it doesnt depend on annotations. It would have been nice to then have seen examples of fusions predicted by tophat-fusion and missed by deFuse because of its reliance on annotated genes. Instead the authors focused on the specificity of deFuse, using results from an incorrect run of the software. A quick email would have resolved this issue. As a rule I believe comparisons should be done with full cooperation of those being compared to, as it would certainly add some reliability to the results.

                      Comment


                      • #26
                        Regarding the comparison to DeFuse in our paper: we ran that at the insistence of the reviewers. As many people know, reviewers often ask for these kinds of "bake-off" evaluations even when not appropriate or necessary. We did not include these in our original manuscript. However I must disagree with andrewm that comparisons "should be done with full cooperation of those being compared to" - this almost never happens. I cannot count the number of times programs of my group have been used in publications as a comparison to some new piece of software, and I don't think I've ever been informed in advance, much less asked to "cooperate."

                        Comment


                        • #27
                          problem with tophat-fusion-post

                          Hi Salzberg,

                          I've been having a lot of trouble trying to get any kind of nontrivial output from tophat-fusion-post, as per my many posts on Seqanswers, even on the MCF7 example given on the webpage. I am trying to get this example to work so that I can try and find some reliable fusion events from purely single-end data and it's been a really frustrating experience.

                          At this stage, I would very much appreciate a detailed command history to make that example work or maybe some tips that you as a seasoned tophat-fusion users could provide on why the output of tophat-fusion-post 2.0.4 (or 2.0.3 or earlier) may be empty despite no problems on alignment with tophat-fusion and there being obvious candidates passing filtration from fusions.out. Assuming I eventually overcome these problems (hopefully!), can you also advise me on how to proceed with purely SE data? E,g,m num-fusion-pairs = 0 must be set, etc.?

                          thanks alot!

                          tm


                          Originally posted by salzberg View Post
                          Several people asked about TopHat. Since last summer (2011), we have released TopHat-Fusion, which can detect fusion transcripts from either single reads or paired end reads. It is very fast and highly effective at filtering out the numerous false positives that plague these types of tools.

                          Comment


                          • #28
                            FusionCatcher

                            Have you tried FusionCatcher?

                            Finder of Somatic Fusion Genes in RNA-seq data. Contribute to ndaniel/fusioncatcher development by creating an account on GitHub.


                            It has found many novel fusion genes in RNA-seq data which were validated later using RT-PCR.

                            FusionCatcher has found plenty of fusion genes in BT-474, SKBR3, KPL4, and MCF7. Some of those fusion genes found by FusionCatcher to be even alternatively spliced and later they were also validated with RT-PCR!

                            FusionCatcher has been used for finding novel and known fusion genes in the following articles:
                            - S. Kangaspeska, S. Hultsch, H. Edgren, D. Nicorici, A. Murumägi, O.P. Kallioniemi, Reanalysis of RNA-sequencing data reveals several additional fusion genes with multiple isoforms, PLOS One, Oct. 2012. http://dx.plos.org/10.1371/journal.pone.0048745
                            - H. Edgren, A. Murumagi, S. Kangaspeska, D. Nicorici, V. Hongisto, K. Kleivi, I.H. Rye, S. Nyberg, M. Wolf, A.L. Borresen-Dale, O.P. Kallioniemi, Identification of fusion genes in breast cancer by paired-end RNA-sequencing, Genome Biology, Vol. 12, Jan. 2011. http://genomebiology.com/2011/12/1/R6

                            Originally posted by salzberg View Post
                            Regarding the comparison to DeFuse in our paper: we ran that at the insistence of the reviewers. As many people know, reviewers often ask for these kinds of "bake-off" evaluations even when not appropriate or necessary. We did not include these in our original manuscript. However I must disagree with andrewm that comparisons "should be done with full cooperation of those being compared to" - this almost never happens. I cannot count the number of times programs of my group have been used in publications as a comparison to some new piece of software, and I don't think I've ever been informed in advance, much less asked to "cooperate."
                            Last edited by ndaniel; 02-03-2013, 03:02 AM. Reason: adding info

                            Comment


                            • #29
                              Originally posted by mjn138 View Post
                              You may also try out FusionMap, from Omicsoft. http://omicsoft.com/fusionmap/

                              It works with single end reads as well. The paper is at http://bioinformatics.oxfordjournals...tr310.abstract

                              It runs on both windows and linux.
                              How can I use FusionMap with hg38? It does not seem to be possible.
                              Any help is highly appreciated

                              Comment


                              • #30
                                Just change the reference library to Human.B38 or Human.hg38. You can see a list of precompiled libraries and gene models here: http://www.arrayserver.com/wiki/inde..._from_OmicSoft

                                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
                                12 views
                                0 likes
                                Last Post seqadmin  
                                Started by seqadmin, Yesterday, 06:07 PM
                                0 responses
                                10 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