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?
•
written
by
wolfgang.rumpf •
1 answer
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.
Log in to answer this question.
-
Novel Splice Variant Discovery from RNAseq Data
written by saipra003 •
-
How to filter reads only with single SNV
written by sunhyunchang •
-
Detection of aberrant mRNA splicing
written by ÁngelMG •
-
HISAT2 not reporting all splice sites?
written by kapoozy •
-
Clarification of good pipeline for transcriptome assembly of RNA-seq data
written by nash.claire
-
Strategies to shortlist True LoF mutations
written by kautilya
-
Which software do you use for RNA-seq data quality control?
written by Christian
-
Finding Novel Transcripts From A De Novo Transcriptome Assembly
written by ilovepython •
-
Where Can I Find Data Sets Of Cancer Publicly Available?
written by Pascal
-
Counting Short-Read Gene Alignments With Bowtie/Tophat.
written by Chris Warth
