Question: How can I filter my RNA-seq transcript data to display ONLY novel isoforms and variants?
0
gravatar for wolfgang.rumpf
7 months ago by
wolfgang.rumpf10 wrote:

I have some RNA-seq data from (human) cancer cells and am trying to find non-canonical transcripts (e.g. splice variants, fusion products, etc.). I have already created a SAM/BAM file of the transcripts and retrieved accession information (using MagicBLAST), but now what I'd like to do is filter out the "known" (canonical) transcripts and leave only the novel transcript variants. It seems like there should be an easy way to do this, but I'm at a loss as to where to find the a single collection of canonical transcripts to filter against....help?

ADD COMMENTlink modified 6 months ago by geek_y10.0k • written 7 months ago by wolfgang.rumpf10
0
gravatar for geek_y
6 months ago by
geek_y10.0k
Barcelona
geek_y10.0k wrote:

Take all the Exon-exon junctions coordinates (i.e intron coordinates i.e retained introns) and check if they are in latest Gencode and FANTOM CAT GTF files. If not, pull out the transcripts that contain those junctions and create a new GTF file.

This will give you all novel transcripts that can have novel exon skipping/inclusion events and also alternate 3' and 5' splice sites.

ADD COMMENTlink modified 6 months ago • written 6 months ago by geek_y10.0k
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: 1021 users visited in the last hour