How to find fusion genes using bioinformatics, Need help?

[ketan_bnf]

Hi all,

I am new to this forum, actually i have 454 sequences and want to find fusion genes (chromosomal translocation) may be causing cancer, right now i am using BLAT, collecting sequences having min. 100bp alignment and above 93% identity.

1) Is there any tool or pipeline that directly predict fusion genes?

2) Is there any perl script to extract result from BLAT output using my parameters?

3) Am i going in right way?

Regards.

[malachig]

Is your data from genomic DNA or RNA? There are probably others, but here are some possibilities off the top:

DeFuse, BreakDancer, Trans-ABySS, NextGENe, some relevant papers here and here, and here.

We have used DeFuse successfully at our center so that would be my first recommendation if you have RNA-seq data...

[ketan_bnf]

Thanks for your reply.


My data is a transcriptome data. I didn't know about BFAST and SAM tools, i'll try it and confirm whether it's apporproate for my data.


Regards,

[townway]

would you share some your experience of defuse?

when I modified config.txt like
# Directory where you have decompressed your dataset
dataset_dir = /gs1/users/defuse-0.3.3/
dataset_prefix = /gs1/users/defuse-0.3.3/Homo_sapiens.GRCh37.59

it always showed the error "does not exist the file /gs1/users/tangwei/defuse-0.3.3/Homo_sapiens.GRCh37.59"


Thank you

Wei

Is your data from genomic DNA or RNA? There are probably others, but here are some possibilities off the top:

DeFuse, BreakDancer, Trans-ABySS, NextGENe, some relevant papers here and here, and here.

We have used DeFuse successfully at our center so that would be my first recommendation if you have RNA-seq data...

[andrewm]

Do you have a space after Homo_sapiens.GRCh37.59
If you would like to email me the config file directly, i can take a look.

Andrew

amcphers
at
bccrc
dot
ca

[townway]

Thank you Andrew.
The main problem for me as you said is the text editor messing everything up.
now it is running.

Do you have a space after Homo_sapiens.GRCh37.59
If you would like to email me the config file directly, i can take a look.

Andrew

amcphers
at
bccrc
dot
ca

[Brajbio]

Did some one has worked on FusionSeq or FusionHunter. How are these tools comparatively with DeFuse in terms of performance and usability???

[Brajbio]

can tophat be used for detection of gene fusion?

[jochensupper]

Brajbio,

most gene fusion tools use paired-end information to detect gene fusion candidates. So in principle you could map with tophat but depending on your downstream gene fusion detection tool, the "split" alignments might not be taken into account. Some tools, however, re-align reads to the detected breakpoint region to determine the exact breakpoint.

can tophat be used for detection of gene fusion?

We (Genomatix) have a gene fusion detection tool that will use "split" reads to detect gene fusion candidates.

[kissthefuture]

would you share some your experience of defuse?

when I modified config.txt like
# Directory where you have decompressed your dataset
dataset_dir = /gs1/users/defuse-0.3.3/
dataset_prefix = /gs1/users/defuse-0.3.3/Homo_sapiens.GRCh37.59

it always showed the error "does not exist the file /gs1/users/tangwei/defuse-0.3.3/Homo_sapiens.GRCh37.59"


Thank you

Wei

hi, I find that your defuse version is so old.Please use 0.3.6 or 0.6.2

[malachig]

I would also suggest you update to 0.4.2. The way these things are specified in the config file have changed slightly. The new example may be more clear. After you have tried that, if you are still having difficulty, email your config file to the author or post it here and we can take a look....

[GenoMax]

Mapsplice may also be an option you can try: http://www.netlab.uky.edu/p/bioinfo/MapSplice

[srividya]

deFuse error.

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

[andrewm]

Your config.txt appears to be fine. However, the perl function abs_path is returning null. What operating system are you using? And which version of perl?

If you want, you can reach me quicker by email at the address listed on this page:

http://defuse.sourceforge.net

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

[mjn138]

FusionMap

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.