SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Cufflinks, differentially expressed genes statsteam Bioinformatics 5 11-15-2013 11:28 AM
Getting differentially expressed genes based on RPKM values casshyr RNA Sequencing 3 04-30-2012 09:18 AM
DESeq and EdgeR: too many differentially expressed genes!?!? cutcopy11 Bioinformatics 5 12-08-2011 12:14 AM
Comparing mouse and human differentially expressed genes stephenhart General 3 11-16-2011 01:14 AM
Detecting differentially expressed genes using aligner outputs questioner Bioinformatics 6 11-03-2011 07:15 AM

Reply
 
Thread Tools
Old 04-08-2012, 10:06 AM   #1
Amative
Member
 
Location: USA

Join Date: Dec 2011
Posts: 45
Default Finding diferentially expressed genes

Hi All,

I have the reads count for the following RNA-Seq samples:
(x1.1, x1.2, x2.1, x2.2, x3.1, x3.2, y1.1, y1.2, y2.1, y2.2, z1.1, z1.2, z2.1, z2.2)

What i need is to find differentially expressed genes in a group of the samples above, example: differentially expressed genes in all x samples but not in the y nor z.

Not: x1.1 & x1.2 are biological replicates.

Please let me know of any way/tool to do this job.

Thank you!
Amative is offline   Reply With Quote
Old 04-08-2012, 03:25 PM   #2
chadn737
Senior Member
 
Location: US

Join Date: Jan 2009
Posts: 392
Default

DESeq
EdgeR
DEGseq
BaySeq


Theres a few for you to check out.
chadn737 is offline   Reply With Quote
Old 04-09-2012, 11:41 AM   #3
Amative
Member
 
Location: USA

Join Date: Dec 2011
Posts: 45
Default

Thanks,
I tried the a couple of them before, but only to find differntiallly expressed genes between only two samples (x1.1, x1.2, y1.1, y1.2)
Now the case is a bit more complicated, i want it as i described above to find those DEG is a group of samples that is not expressed in the remaining. If any of these packages provide this can you please let me know which one and if you have links to any examples that looks like mine I would really appreciate that.
Amative is offline   Reply With Quote
Old 04-09-2012, 11:52 AM   #4
chadn737
Senior Member
 
Location: US

Join Date: Jan 2009
Posts: 392
Default

As far as I know, nobody has yet released (I know of groups who are working on such approaches) a package to deal with data in that fashion. Although if anyone here knows of one, I would also be interested in knowing.

Otherwise, as far as I know, all approaches require a pairwise comparison between different treatments.

So the approach I would use would be to do the pairwise comparisons and then filter accordingly.

So find your differentially expressed genes between x-y, x-z, y-z then find genes that are differentially expressed in both x-y and x-z and filter out genes that are differentially expressed in y-z. You should be left with a genes that are only differentially expressed between x and both y and z conditions, but which are unchanged in between your y and z conditions.

We actually used a similar approach to this when comparing transcriptomes between 3 different varieties of the same species, but looking for differences in only one of the conditions that could explain a phenotype.
chadn737 is offline   Reply With Quote
Old 04-09-2012, 06:06 PM   #5
Amative
Member
 
Location: USA

Join Date: Dec 2011
Posts: 45
Default

Thanks for the details

When you say "find your differentially expressed genes between x-y", do you mean finding differentially expressed genes in all x conditions (x1.1, x1.2, x2.1, x2.2, x3.1, x3.2) vs all y conditions (y1.1, y1.2, y2.1, y2.2) ???

Or do I have first to find out those differentially expressed between x conditions:
(x1.1,x1.2) vs (x2.1 vs x2.2)
(x1.1,x1.2) vs (x3.1 vs x3.2)
(x2.1 vs x2.2) vs (x3.1 vs x3.2)
and then the same for y conditions
(y1.1, y1.2) vs (y2.1, y2.2).

Up to my knowledge i can only use the packages above with only two conditions at a time
(x1.1,x1.2) vs (x2.1 vs x2.2)

I am not sure if there is a way to do two samples with more than one condition at the same time, as follow,
(x1.1, x1.2, x2.1, x2.2, x3.1, x3.2) vs (y1.1, y1.2, y2.1, y2.2) ???

Or that is considered same as my first question!!

Thanks
Amative is offline   Reply With Quote
Old 04-09-2012, 06:19 PM   #6
chadn737
Senior Member
 
Location: US

Join Date: Jan 2009
Posts: 392
Default

Quote:
Originally Posted by Amative View Post
Thanks for the details

When you say "find your differentially expressed genes between x-y", do you mean finding differentially expressed genes in all x conditions (x1.1, x1.2, x2.1, x2.2, x3.1, x3.2) vs all y conditions (y1.1, y1.2, y2.1, y2.2) ???

Or do I have first to find out those differentially expressed between x conditions:
(x1.1,x1.2) vs (x2.1 vs x2.2)
(x1.1,x1.2) vs (x3.1 vs x3.2)
(x2.1 vs x2.2) vs (x3.1 vs x3.2)
and then the same for y conditions
(y1.1, y1.2) vs (y2.1, y2.2).

Up to my knowledge i can only use the packages above with only two conditions at a time
(x1.1,x1.2) vs (x2.1 vs x2.2)

I am not sure if there is a way to do two samples with more than one condition at the same time, as follow,
(x1.1, x1.2, x2.1, x2.2, x3.1, x3.2) vs (y1.1, y1.2, y2.1, y2.2) ???

Or that is considered same as my first question!!

Thanks
According to your opening post, all your x samples are biological reps, all your y samples are biological reps, and dido for all your z samples.

You NEED biological reps to do differential expression. So what I mean by doing pairwise comparisons of the conditions is all biological reps of x versus all biological reps of y and so on. There is no point in separating these replicates as they are what allows you to estimate the biological variance between samples.
chadn737 is offline   Reply With Quote
Old 04-09-2012, 08:01 PM   #7
Amative
Member
 
Location: USA

Join Date: Dec 2011
Posts: 45
Default

I am sorry, I was not clear enough before,

All x samples are from the same tissue but different stages
(x1.1,x1.2) are biological replicates of tissue x stage 1
(x2.1,x2.2) are biological replicates of tissue x stage 2
.
.

(y1.1,y1.2) are biological replicates of tissue y stage 1
(y2.1,y2.2) are biological replicates of tissue y stage 2
.
.
.

I hope it is more clear now!
Amative is offline   Reply With Quote
Old 04-09-2012, 08:13 PM   #8
chadn737
Senior Member
 
Location: US

Join Date: Jan 2009
Posts: 392
Default

Ah, this vastly complicates things. Can you clarify then exactly what your aim is. Since X1 is different from X2 are you interested only in a specific stage or a specific condition regardless of stage?
chadn737 is offline   Reply With Quote
Old 04-09-2012, 09:44 PM   #9
Amative
Member
 
Location: USA

Join Date: Dec 2011
Posts: 45
Default

Yea, I know it is complicated !!

Actually, I need both,
Find those differentially expressed between any two conditions
(x1.1,x1.2) vs (y1.1,y1.2)
which can be easily done with the packages listed above.

The one I am still confused how to do it, is finding those genes differentially expressed genes in all x conditions but not in any of y conditions.

1) I'm thinking of finding the DEG between each possible pair of conditions from sample x
then
2) I'm thinking of finding the DEG between each possible pair of conditions from sample y

then
3) Union all results from 1
4) Union all results from 2

find out those genes in 3 and not 4


what do you think?!! make sense?
Amative 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 06:56 PM.


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