This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Identify the isoform of a transcript using RNA Seq data

Hi,

my favorite transcript is supposed to consist of 12 exons (it was only described in germ cells, not in somatic tissues). Thus, I want to check what isoforms are actually expressed in my (somatic) tissue of interest.

I have paired-end RNA Seq data from a larger number of samples. I would like to use this dataset to shed light onto the isoform.

Thus far, I have aligned my reads using STAR, and counted exons using featureCounts. Seemingly, only a few exons are expressed, but quite consistent over all samples. There is no described isoform of the transcript that looks like what I have.

How would you suggest to go from here to obtain a consensus sequence? Do I need to assemble the whole transcriptome, or are there any shortcuts I can take, as I am only interested in a single transcript?

Many thanks!

rna-seq
• 3,888 views • link
updated by EagleEye EagleEye • written by aggregatibacter aggregatibacter

3 answers

I'd subset the alignments to contain only that regions and then run stringTie.

• 1 views • link
written by Devon Ryan Devon Ryan

That sounds like a great idea, many thanks, Devon!

• 1 views • link
written by aggregatibacter aggregatibacter

Did you try running Cufflinks2 RABT assembly; I would suggest you to run in both mode; default and RABT and see if you can catch your desired isoform

• 1 views • link
written by Manvendra Singh Manvendra Singh

These posts might give you some idea,

A: How to determine alternative splicing read counts

A: workflow/tool for alternative splicing

• 1 views • link
written by EagleEye EagleEye

Log in to answer this question.

More posts like this