SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Significant up-regulated and downregulated genes with DESeq alisrpp Bioinformatics 8 10-29-2013 08:43 AM
how to identify differentially expressed genes from cuffdiff output? IBseq Bioinformatics 3 08-07-2013 09:06 AM
How does Cufflinks determine Differentially Regulated Genes (relCDS)? jmwhitha Bioinformatics 1 06-22-2013 02:22 PM
How to identify diff from RNA-seq data set with DNA pollution xfliwz Bioinformatics 0 01-06-2013 11:11 PM
ChIP-Seq: Tissue-specific prediction of directly regulated genes. Newsbot! Literature Watch 0 07-05-2011 03:00 AM

Reply
 
Thread Tools
Old 10-23-2013, 02:48 AM   #1
sridhar28
Member
 
Location: chennai

Join Date: May 2013
Posts: 15
Default How to identify the Up and down Regulated genes from Cuff diff output

Dear All,

How to identify the Up and down regulated genes from the diff genes, I got diffgenes using Cummerbund codes.

In my diff_gene.txt file could see most of the log2fold change values are +inf and -inf. From the list how to separate the up and down regulated genes??

Cheers

Sridhar

Last edited by sridhar28; 10-23-2013 at 08:00 AM. Reason: link not attached
sridhar28 is offline   Reply With Quote
Old 12-26-2013, 08:31 AM   #2
bassounds
Junior Member
 
Location: New York, NY

Join Date: Nov 2013
Posts: 1
Default

Hey Sridhar,

I am also getting similar results (+inf and -inf) on some of the log2fold change values. I believe this is happening because one of the conditions has a value of 0.

I am reading through cuffdiff documentation to sort out whether this is due to normalization or processing issue, or whether I actually have zero reads for some transcripts. I doubt the latter is the case, based on my data of identical cell lines with different treatments, and high coverage of the ran-seq reads.
bassounds is offline   Reply With Quote
Old 01-17-2014, 02:09 AM   #3
velt
Member
 
Location: Paris

Join Date: Jun 2013
Posts: 10
Default

Hello,

In my view, ​​-inf and inf values, are effectively due to a division by 0 when calculating the fold-change.

For example in a test_vs_control comparison, inf value for a gene means that there is many reads in test condition and 0 in the control. And conversely, -inf value for a gene means that there is 0 read for this gene in test condition and many in control condition.

Do not ignore these genes, because if the p-value is significant, they can be considered deregulated.

For me, there are two solutions, a "clean solution" and "less clean solution".

- clean solution => replace the 0 values by a non-zero value (1 for example) and recalculate the fold-change for these genes.

- less clean solution, but simpler => if the exact value of the fold-change is not important to you, and you only use it to filter the differentially expressed genes, just replace the inf & -inf values by a large value (-1000 and 1000 for example). Thus, you can select those genes that have a high fold-change (if the p-value is significant of course).

Otherwise, to select genes differentially expressed, you can filter diff_gene.txt file according to two criteria, the log2-fold change and p-value. Generally a gene is considered significant if its p-value is less than 0.05, but it is up to you to set a threshold. In addition, you can consider that a gene is truly differentially expressed if its log2 fold-change is > 1 or < -1. But again, it's up to you to set that threshold.

Thus, down-regulated genes have a fold-change <= -1 and a p-value <= 0.05 and up-regulated genes have a fold-change >= 1 and a p-value <= 0.05.

I hope these ideas will help you.
(PS: sorry if my english is bad, this language is not my native language)
velt is offline   Reply With Quote
Reply

Tags
cuffdiff 2.0.2, cummerbund, 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 04:59 PM.


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