How Can I Find Transcripts Which Omit Specific Exons Of A Gene?
1
1
Entering edit mode
10.1 years ago
Dan D 7.3k

Most of my work has been on NGS QC, but lately I've been getting approached with queries for downstream analysis. An investigator was interested in evaluating differential expression of a small set of genes. I ran the investigator's RNASeq datasets through the Tuxedo suite and also through a BWA -> DESeq pipeline to get familiar with established techniques.

Now the investigator wants to focus on expression of a particular murine gene, but only on transcripts which are missing two specific exons. I have the consensus sequences of these exons.

What I'm not sure about is the general strategy to obtain this particular expression information. Having some sort of idea about which tools to investigate would be especially helpful.

rna-seq transcript expression exon • 2.5k views
ADD COMMENT
1
Entering edit mode

Just a clarifying question: are you trying to identify genes that include/exclude exons between two conditions? Or are you saying that you know the gene you want to look at which expressed several isoforms, and you are trying to deconvolve the expression of each isoform from the "total" expression of the gene?

ADD REPLY
0
Entering edit mode

Good question, Steve. In my case, it's the latter of the two possibilities. Thanks!

ADD REPLY
3
Entering edit mode
10.1 years ago
JC 13k

You can try to infer the transcript expressions for each allele with Cufflinks or check differential exon expression with DEXseq in R:Bioconductor.

ADD COMMENT
0
Entering edit mode

Thanks! Looks like DEXseq may be the ticket. I'll have a look right now.

ADD REPLY
1
Entering edit mode

FWIW, if you are really looking to do the "latter" part of what I said as a comment to your original question, then DEXSeq won't be the ticket. DEXSeq is appropriate when you want to identify exons that are differentially expressed (the first part of what I said in the comment)

ADD REPLY
0
Entering edit mode

Cool, thanks for the heads up. I'm deep in the documentation right now. I'm also reading through the Tuxedo suite docs, but I haven't yet determined a strategy.

ADD REPLY
1
Entering edit mode

if you want to deconvolve isoform expressions, then Cufflinks it's the ticket

ADD REPLY
0
Entering edit mode

Excellent. I got sidetracked on other projects but I will revisit the cufflinks output to see what I need to do. I'll report back on this thread when I've made progress. Many thanks!

ADD REPLY
0
Entering edit mode

I have the same problem, DEXseq will tell which exons are different but not which junctions are different. Other than manually looking at the bam files in IGV, is there a computational way to see how many reads span 1-3 rather than 1-2-3?

ADD REPLY
1
Entering edit mode

maybe you can extract the reads in the region and separating them by the CIGAR information with some script

ADD REPLY

Login before adding your answer.

Traffic: 829 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6