I am trying to figure out which isoform of a protein is expressed in my cells. I have reads in exons suggesting that either multiple isoforms are expressed, or a yet undiscovered isoform is expressed. To test this I need a tool that uses paired end reads across different exons in order to 'link' those exons together, so that I know that these exons are expressed as one isoform.
To elaborate on the situation:
Isoform1: |---[exon1]-----------------[exon3]---|
Isoform2: |---[exon1]---[exon2]-----------------|
my track: |---[exon1]---[exon2]---[exon3]---|
So there is a chance that isoform 1 and 2 are both expressed, or that my track is actually a new isoform. If I have reads with one end in exon 1 and the other end in exon 2, I know that isoform 2 is expressed. If one end is in exon 2 and the other in exon 3 I know that a yet unknown isoform is expressed.
Can someone suggest a tool to do this type of analysis?
•
written
by
thomasterTW •
2 answers
Cufflinks is a bit obvious so you have probably tried it already.
I would recommend a full de novo analysis using eg Trinity / CLC and map back to the genome with GMAP. Then you can see any new isoforms quite easily.
Log in to answer this question.
-
How to identify first codon (Methionine) misssense variants?
written by lincaijin1994 •
-
How to interpret the exon numbering for minus strand?
written by sebastian.swirski •
-
Run featureCounts with Gencode GTF
written by tianleivv •
-
Will pseudoaligners detect rearrangements based on known transcripts?
written by b10hazard •
-
De novo RNA-seq
written by blawney •
-
Alternative splicing: explanation about delta psi values
written by user31888
-
GFF3 ID field: "Official" nomenclature
written by user230613
-
Converting gtf format to bed format
written by thjnant
-
Differential expression level of different exons from RNA-seq indicates what?
written by biomedicineman1 •
-
Generating Internal Exon Based On Known Mrna Isoform
written by Puriney
