AlgoLab / galig

A graph aligner
GNU General Public License v3.0
26 stars 12 forks source link

Novel intron retention events in genome-wide mode #17

Closed raerose01 closed 3 years ago

raerose01 commented 3 years ago

Hello, I had a more general question about running ASGAL in genome-wide mode.

I have a dataset where I know there are 3 novel retained intron events. I was wondering if, after the pre-filtering step performed with by quasi-mapping with Salmon, the reads mapping to these novel retained introns would still be included in down-stream alternative splicing analysis.

I'm asking because I have run ASGAL in genome-wide mode, but fail to detect any events. When I look at the output SAM file for a gene that should have a retained intron, there is no coverage (whereas when I look with a more typical spliced alignment to the reference genome tool, e.g. STAR, there is coverage). I have attached a screenshot of what I mean.

Screen Shot 2021-06-17 at 3 19 55 PM

Could you let me know if I am mis-understanding something, or should be running the tool differently?

Thanks!

Rachel

ldenti commented 3 years ago

Hi Rachel, unfortunately the current version of ASGAL cannot detect that kind of intron retention. Indeed, ASGAL aligns reads to the splicing graph (in other words to the known exons of a gene): in your case, there are a lot of reads mapping to an intron (as you can see from the coverage plot).

In the paper we said:

We note that the current implementation of ASGAL is not able to detect the insertion of novel exons inside an intron and intron retention events caused by the union of two exons

Your case is exactly the second one we mentioned. Maybe it isn't clear enough: in case, sorry for that.

Best, Luca

raerose01 commented 3 years ago

Hi Luca,

OK, thanks for clarifying that.

Best,

Rachel