SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Can I use DESeq and edgeR for mixed ANOVA abebe Illumina/Solexa 18 07-13-2017 05:45 AM
Can edgeR/DESeq have more than one covariate? arrchi Bioinformatics 8 10-28-2013 03:37 PM
DESeq and edgeR up/down regulation murphycj Bioinformatics 7 09-21-2011 08:09 AM
BaySeq vs GLM in EdgeR and DESeq Hilary April Smith Bioinformatics 0 08-03-2011 11:14 AM
edgeR vs DESeq vs bayseq Azazel Bioinformatics 1 10-07-2010 08:11 AM

Reply
 
Thread Tools
Old 08-10-2011, 07:35 AM   #1
PFS
Member
 
Location: USA

Join Date: Mar 2010
Posts: 55
Question comparing results by cuffdiff, edgeR, DESeq

Has anyone tried to compare the results from the various tools that offers differential expression analyses for RNASEQ data?

I understand they have different underlying models and assumptions, but I would expect some overlap. At a first glance, when I tried to compare the results on my data, I got completely different DE genesets and I am puzzled.
PFS is offline   Reply With Quote
Old 08-12-2011, 06:17 AM   #2
chenyao
Member
 
Location: Beijing

Join Date: Jul 2011
Posts: 74
Default

Quote:
Originally Posted by PFS View Post
Has anyone tried to compare the results from the various tools that offers differential expression analyses for RNASEQ data?

I understand they have different underlying models and assumptions, but I would expect some overlap. At a first glance, when I tried to compare the results on my data, I got completely different DE genesets and I am puzzled.
This is quite interesting question?

I also want to know the answer. With my experence on microarray, different statistcal method may generate different DEG list, but they should have a lot overlap.
chenyao is offline   Reply With Quote
Old 08-12-2011, 11:52 AM   #3
mbblack
Senior Member
 
Location: Research Triangle Park, NC

Join Date: Aug 2009
Posts: 245
Default

Quote:
Originally Posted by PFS View Post
Has anyone tried to compare the results from the various tools that offers differential expression analyses for RNASEQ data?

I understand they have different underlying models and assumptions, but I would expect some overlap. At a first glance, when I tried to compare the results on my data, I got completely different DE genesets and I am puzzled.
I analyzed the same RNAseq data set using cuffdiff, DESeq and EdgeR. The data was a whole transcriptome SOLiD 4 fragment run for 3 control and 3 treatment animals (mice), so the final mapped reads ranged from approx. 71M reads, to a high of approx. 171M reads (because of the layout of the barcoding on the slides they were not all equivalent). Raw counts were extracted from BioScope mapped reads using BAMtools and a UCSC RefGene bed file as reference.

Going solely by FDR < 0.05 as a cutoff:

EdgeR - only 217 significantly differentially expressed genes
DESeq - 337 significantly differentially expressed genes
There were 198 genes in common in those two lists.

Cuffdiff gave 202 significantly differentially expressed genes by q-value, but I don't know offhand how many of those are genes common to the other two.

The issue that has me really tearing my hair out is that I have microarray (affy) data for these same animals. A LIMMA analysis of those 6 arrays gives some 3000 significant genes by FDR < 0.05. With the coverage I have with the RNAseq data, I should at least be comparable with the array data, not left with a difference in significant genes of 1000s.

We were expecting to approach or exceed array results with about one tenth of the mapped reads I have for these initial runs.
mbblack is offline   Reply With Quote
Old 10-08-2011, 02:05 AM   #4
songyj
Member
 
Location: china

Join Date: Sep 2011
Posts: 15
Default how wonderful

Quote:
Originally Posted by mbblack View Post
I analyzed the same RNAseq data set using cuffdiff, DESeq and EdgeR. The data was a whole transcriptome SOLiD 4 fragment run for 3 control and 3 treatment animals (mice), so the final mapped reads ranged from approx. 71M reads, to a high of approx. 171M reads (because of the layout of the barcoding on the slides they were not all equivalent). Raw counts were extracted from BioScope mapped reads using BAMtools and a UCSC RefGene bed file as reference.

Going solely by FDR < 0.05 as a cutoff:

EdgeR - only 217 significantly differentially expressed genes
DESeq - 337 significantly differentially expressed genes
There were 198 genes in common in those two lists.

Cuffdiff gave 202 significantly differentially expressed genes by q-value, but I don't know offhand how many of those are genes common to the other two.

The issue that has me really tearing my hair out is that I have microarray (affy) data for these same animals. A LIMMA analysis of those 6 arrays gives some 3000 significant genes by FDR < 0.05. With the coverage I have with the RNAseq data, I should at least be comparable with the array data, not left with a difference in significant genes of 1000s.

We were expecting to approach or exceed array results with about one tenth of the mapped reads I have for these initial runs.
Wow..
I only used Tophat/cufflinks ,2 ,only 2 significantly differentially expressed genes in isoform_exp.diff
seems like maybe I could try DESeq, hope it will give me more differentially expressed genes
songyj is offline   Reply With Quote
Old 03-08-2014, 06:09 AM   #5
super0925
Senior Member
 
Location: UK

Join Date: Feb 2014
Posts: 206
Default

Hi All I am a rookie in RNA-seq.
I found a problem that gene ID output from Cuffdiff and egdeR/DESeq are different so that I cannot find the common DE genes.
I use the data set with the Drosophila melanogaster genome.
The edgeR output gene ID (e.g. FBgn0000370,FBgn0000500,…)
But the Cuffdiff output gene ID is something like (XLOC_000028,XLOC_000038,…) and gene symbol (e.g.,KH1,RpLP1,…).
Could you please recommend any software or command to translate them automatically? Could you have command line for me to use to uniform or join them?
super0925 is offline   Reply With Quote
Old 03-12-2014, 04:01 AM   #6
Thias
Member
 
Location: Münster, Germany

Join Date: Mar 2013
Posts: 44
Default

BioMart is an option. I frequently use it as R package ("biomaRt"), but there is also an easy to use online tool.
Thias is offline   Reply With Quote
Reply

Tags
cuffdiff, degseq, deseq, edger, rnaseq

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 02:12 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO