Use of StringTie to assemby of unannotated intron retention transcripts
CBenetti opened this issue · 1 comments
CBenetti commented
The aim is to be able to obtain the full transcript sequence in case of unannotated intron retention events. Following the StringTie workflow, this would mean considering retained introns (which have been previously identified with IRFinder) as potential new exons, as to reconstruct the full transcript and quantify them trough the classical pipeline, while still being able to use existing annotations.
A) The question is: is there an existing setting which I could implement to obtain such result, or is it something I should costumize?
B) If it isn't an existing feature of StringTie, could the following be a solution:
- When generating BAM files, adding in STAR splice junctions corresponding to the exon-intron junction as to generate XS tags for those reads spanning the exon - intron junction for both files
- Add to gtf annotation file coordinates of retained introns and gene of appartenece
- Running StringTie with custom annotations and custom BAM files
Thank you in advance
CuteGold0407 commented
hi,have you solved this issue?