Looking For Reasons Of Why The Results Of A Tophat Fusion Post Is Empty
3
0
Entering edit mode
8.6 years ago
KS ▴ 380

Hello Everyone,

I ran tophat-2.7 on the example data with --fusion-search , it ran successfully and produced the following output files: accepted hits.bam, deletions.bed, insertions.bed, junctions.bed, fusions.out

I am trying to run tophat-fusion-post, after I run it, the output is empty. Can anyone suggest the reason for it?

The command I used to run tophat (I am retuning this command under topdir directory):

tophat -o tophatMCF7 -p 8 --fusion-search /path/to/hsapiens/bowtieindex/hg19 SRR0642861.fastq SRR0642862.fastq

The structure I am using to run tophat is the following:

top_dir:

blast
blasthuman (nt and humangenomic databases)
indexes (hg19)
mcl
refGene.txt
ensGene.txt
SRR0642861.fastq 
SRR0642862.fastq
tophatMCF7 (output directory): accepted_hits.bam, deletions.bed, insertions.bed, junctions.bed, fusions.out, logs, prep_reads.info, unmapped.bam
tophat fusion • 4.8k views
ADD COMMENT
0
Entering edit mode

have you checked the log files? I did not put human blast databases inside folder blasthuman. I extracted all archives under blast directory.

ADD REPLY
1
Entering edit mode
8.5 years ago

I think I found the problem (after many errors). On line 696 of the script tophat-fusion-post, I saw:

        if string.find(file, "tophat_") != 0:

And I rename the folders tophat_ (in the above case tophat_MCF7), it seems to work! I hope that can be useful to others!

ADD COMMENT
0
Entering edit mode

Thanks, finally a ray of hope :-)

ADD REPLY
0
Entering edit mode
8.6 years ago

You're not the only one! See these SeqAnswers threads, maybe they can provide some help:

Repost: tophat-fusion outputs empty result

tophat-fusion-post result empty

I have been unable to solve the problem in my own case, and instead moved on to other fusion finding programs like FusionCatcher, ChimeraScan, FusionMap, and FusionHunter (which I have all managed to run, although not without some trial-and-error).

ADD COMMENT
0
Entering edit mode
7.3 years ago
Paul ★ 1.4k

I have got it - mcl (known gene fusions) can be obtained from previous tophatfusion version (in annotation directory), http://tophat-fusion.sourceforge.net/downloads/tophatfusion-0.1.0.tar.gz

ADD COMMENT

Login before adding your answer.

Traffic: 2991 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6