Differential expression in edgeR with cuffdiff count output
2
0
Entering edit mode
7.0 years ago
mjg ▴ 30

Hello,

I have rna-seq data from paired samples and want to do isoform-level differential expression. 

My understanding is that is not possible to explicitly specify paired design in cuffdiff, so I considered the option of using the isoforms.read_group_tracking file and extract the 'raw_frags' column, which I *think* corresponds to counts. Then use this as input in edgeR where I can specify paired-sample design.

Does this sound sensible?

If not, my question is: which package can I use to do isoform-level DE which takes into consideration paired-sample designs.

 

Thanks,

Maria

 

RNA-Seq edgeR cuffdiff • 2.4k views
ADD COMMENT
1
Entering edit mode
6.9 years ago
mark.ziemann ★ 1.4k

DEX-Seq Estimates exon usage coefficients from the fitted terms of the GLM and may be just what you're looking for.

http://www.bioconductor.org/packages/release/bioc/manuals/DEXSeq/man/DEXSeq.pdf

ADD COMMENT
0
Entering edit mode

Thank you for your suggestion! I did not know about DEX-Seq, I am running it now and will compare it with my previous results. 

Maria
 

ADD REPLY
1
Entering edit mode
6.9 years ago
Lior Pachter ▴ 610

You can use sleuth.

ADD COMMENT
0
Entering edit mode

I have played sleuth. But if I am interested with DE genes but not DE transcripts in downstream analysis, how could I do it with sleuth?

e.g. I got the DE transcripts table with Q<0.05 with associated gene name. But as you know, one gene could have multiple transcripts.

like

target_id          Q                        ext_gene
ENST00000257570    0                        OASL
ENST00000339275    1.08823325997074e-138    OASL
ENST00000543677    4.50865532164411e-20     OASL

Could I use it for finding DE genes like edgeR/DESeq? how?

Thank you!

ADD REPLY

Login before adding your answer.

Traffic: 1833 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