Question: Fusion-gene finding - Is longer read-length always better?
2
gravatar for ljq
4.9 years ago by
ljq20
Singapore
ljq20 wrote:

Hi all,

I am very new to fusion-gene discovery and I am tasked along with some fellow colleagues to apply bioinformatic methods on clinical data for the detection of novel fusions among subtypes/groups of lymphomas patients.

We have a few sets of 2x100bp RNA-seq (~2x75Million reads) dataset and we applied FusionCatcher onto them. After some FastQC on the input reads, we found out that the first ~15 bases of each reads may be contaminated by 'not-so-random' priming or could be even be from sequencing the adapter.

We subsequently applied FusionCatcher to the hard-trimmed-first-20bps (lets call the candidate fusion from this set - B) and hard-trimmed-first-20bp-back-30bp datasets (result - C) and got different sets of results. Finally, we call the results on applying FusionCatcher with the original 2x100bp reads as A.

The absolute counts of candidate fusions in the results decrease with decreasing read-lengths. C is a power-subset of B. However, there are some detected candidate which are exclusive to B as compared with A.

After seeing this, i am wondering if there is a 'perfect' read-length to work with. (The fusions exclusive to B are under validation and i will update everyone with the results in the time to come)

What are your experiences on dealing with RNA-seq and fusion-gene finding?

--jq

rna-seq fusion • 1.8k views
ADD COMMENTlink modified 4.9 years ago by Devon Ryan88k • written 4.9 years ago by ljq20
1

When you say "hard-trimmed-first-20bps", are you referring to the reads after lopping off the first 20 bases or to the first 20 bases that show bias? I assume the former, but what you wrote technically means the latter.

BTW, why not just run things through a standard read trimmer to remove adapters and use the results? The non-random priming just shows that reads aren't completely randomly selected, but the biased part is still part of the legitimate sequence.

ADD REPLYlink written 4.9 years ago by Devon Ryan88k

I was applying the trim to all the reads; regardless of them being biased or not.

The program which i was using, FusionCatcher, remove adaptors automatically (i am not sure of the removal performance). Devon, I will look into Trim-Galore and Trimmomatic. You are right about the non-random priming of sequences. Thanks!

ADD REPLYlink written 4.9 years ago by ljq20

FusionCatcher is using this to find automatically the adapter and remove them:

http://code.google.com/p/fusioncatcher/source/browse/remove_adapter.py

ADD REPLYlink written 4.9 years ago by bstrs0

cross posted: http://seqanswers.com/forums/showthread.php?t=43244

ADD REPLYlink written 4.9 years ago by Pierre Lindenbaum117k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2329 users visited in the last hour