SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
comparing results by cuffdiff, edgeR, DESeq PFS Bioinformatics 5 03-12-2014 04:01 AM
edgeR ECHo Bioinformatics 13 04-18-2013 02:01 PM
EdgeR: Heatmaps polsum Bioinformatics 5 02-04-2012 12:47 PM
edgeR p-value and logFC alessandra85 Bioinformatics 3 08-29-2011 06:04 AM
edgeR Puva Bioinformatics 2 05-19-2011 10:04 AM

Reply
 
Thread Tools
Old 10-22-2010, 02:41 PM   #1
Balat
Member
 
Location: Australia

Join Date: May 2010
Posts: 36
Default cuffdiff vs edgeR

Hi,
I have used the new version of cuffdiff with biological replicates. After identifying differentially expressed genes, I did gene ontology based enrichment trests. I did not find any significant enriched categories. Then I used Tophat with all my libraries to get one bam file and used this in cufflinks to get one gtf file. I did this as we don't have a good annotations for the organism that I am working on. Then I used Tophat separately with each library and got bam files for each library. I used the bam file and the gtf file with bedtools and got read counts for each of library separately. I used these read count in edgeR to test differential gene expression. Gene ontology tests with differentially expressed gene showed more than 100 significantly cotegories and the up and downregulated gene categories are all relevant to the treatment I used in my experiment.

To directly compare the results from cuffdiff, I obtained read counts based on the coordinates from gene exp file from cuffdiff. I used these read counts in edgeR and used the DE genes in gene ontology enrichemnt tests. I obtained several significant gene categories and all of the enriched categories are relevant for the treatment of my experiment.

I am just wondering has anyone seen such discrepancy between cuffdiff and other read count based gene expression tests?
Balat is offline   Reply With Quote
Old 10-22-2010, 03:28 PM   #2
Cole Trapnell
Senior Member
 
Location: Boston, MA

Join Date: Nov 2008
Posts: 212
Default

Can you be a bit more detailed? How many genes were reported significantly different by cuffdiff? How did you run cuffdiff?
Cole Trapnell is offline   Reply With Quote
Old 10-22-2010, 04:42 PM   #3
Balat
Member
 
Location: Australia

Join Date: May 2010
Posts: 36
Default

Total 20700 genes were were tested with cuffdiff and 5000 genes were found to be significant between two treatments. I ran cuffdiff using gtf file obtained from cuffcompare. I had 3 replicates for two of my treatments. I used bam files from 3 replicates from each treatment together as comma separated files in cuffdiff. I extracted read counts from the gene exp file using gene coordinates with bedtools. Total 15500 genes were used for extracting read counts. As cuffdiff estimates several transcripts from same gene coordinates, the number of genes with read counts were less when I used gene coordinates for extracting gene counts. I used these read counts with edgeR.

I used top most 2000 significant genes from cuffdiff in testing for gene ontology enrichment. Similarly I used top 2000 significant genes from edgeR for testing gene category enrichment.
Balat is offline   Reply With Quote
Reply

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 07:07 AM.


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