Hi,
I work on alternative splicing with RNA-Seq data. Recently, while mapping reads using Tophat version 1.31 in order to characterize exon skipping events, I found a peculiar problem. I have tried to summarize the issue with this figure:
The problem is this: Most of the times, when the last segment is less than 10 bases (during segment match), tophat allows for a default of up to 2 mismatches for each segment and so this small segment also gets 2 mismatches. Consider the scenario where the last segment is 4 bases = GTTG, as shown in the figure. The read clearly matches the correct exon (with no skipping), but in both cases, tophat maps with 1 and 2 mismatches leading to an exon-skipping event. In this particular case, I had around 1000 such reads exon skipped out of about 10000+ reads that mapped normally across these exons (no skipping).
1) Is this a bug? Have you experienced this issue? If so, how did you get around this issue?
I have two things in mind to get around this issue:
1) I am planning to take all these "special" events or junctions that are not in the annotation and then look for the start and end junctions (as to which genes they map to) and then look if the "current" end matches any of the other exons in the current gene. The good thing is that this can help with any software. However, its laborious as you have to check fasta sequences for every such read.
2) Another way is to run tophat initially with 0 segment mismatches and then obtain the unmapped reads and convert them to FASTQ again and then map again with 1 mismatch and then do the same for 2 mismatches. I am not sure if this will help though.
I'd appreciate any thoughts on this. Thank you.
I work on alternative splicing with RNA-Seq data. Recently, while mapping reads using Tophat version 1.31 in order to characterize exon skipping events, I found a peculiar problem. I have tried to summarize the issue with this figure:
The problem is this: Most of the times, when the last segment is less than 10 bases (during segment match), tophat allows for a default of up to 2 mismatches for each segment and so this small segment also gets 2 mismatches. Consider the scenario where the last segment is 4 bases = GTTG, as shown in the figure. The read clearly matches the correct exon (with no skipping), but in both cases, tophat maps with 1 and 2 mismatches leading to an exon-skipping event. In this particular case, I had around 1000 such reads exon skipped out of about 10000+ reads that mapped normally across these exons (no skipping).
1) Is this a bug? Have you experienced this issue? If so, how did you get around this issue?
I have two things in mind to get around this issue:
1) I am planning to take all these "special" events or junctions that are not in the annotation and then look for the start and end junctions (as to which genes they map to) and then look if the "current" end matches any of the other exons in the current gene. The good thing is that this can help with any software. However, its laborious as you have to check fasta sequences for every such read.
2) Another way is to run tophat initially with 0 segment mismatches and then obtain the unmapped reads and convert them to FASTQ again and then map again with 1 mismatch and then do the same for 2 mismatches. I am not sure if this will help though.
I'd appreciate any thoughts on this. Thank you.
Comment