Question: Time Course Rna-Seq Analysis
9
gravatar for kanwarjag
6.7 years ago by
kanwarjag1.0k
United States
kanwarjag1.0k wrote:

I have a data with three replicates for 4 time points and underlying aim is to find signature genes for a particular phenotype which will occur after 3rd time point and goes down at 4th time point. I was wondering how can I analyze this, Should DEseq package (count based ) method be of some use I have not seen any time course study with RNA-seq using Deseq. Any other suggestion or pointers will be helpful.

Thanks

replicates rna-seq • 7.6k views
ADD COMMENTlink modified 5.9 years ago by Charles Warden7.4k • written 6.7 years ago by kanwarjag1.0k
1

are you testing the hypothesis that there is a non-zero change in expression (slope) or that the slopes differ over time between two groups?

ADD REPLYlink written 6.7 years ago by Dave Bridges1.3k

My aim is to find out the genes which are specifically associated with particular phenotype over time. Yes in other word genes which are nonzero over time course

ADD REPLYlink modified 6.7 years ago • written 6.7 years ago by kanwarjag1.0k
7
gravatar for Dave Bridges
6.7 years ago by
Dave Bridges1.3k
Ann Arbor, MI
Dave Bridges1.3k wrote:

You should be able to use fitNbinomGLM and nbinomGLMTest from the DESeq package:

library(DESeq)
cds <- makeExampleCountDataSet()[ 1:100, ]
cds <- estimateSizeFactors( cds )
cds <- estimateDispersions( cds, method="pooled" ) #need to use either pooled or blind
time.course <- c(0,2,5,15,24)
fit1 <- fitNbinomGLMs( cds, count ~ time.course ) #fits to the formula that the time.course is a covariate
fit0 <- fitNbinomGLMs( cds, count ~ 1 ) #the null hypothesis
nbinomGLMTest( fit1, fit0 ) #chi squared test for each gene

I am not 100% sure, but i think you will also need to adjust the results of nbinomGLMTest for multiple observations.

If you need a time course plus a condition then just modify the formula to add that variable

ADD COMMENTlink modified 6.7 years ago • written 6.7 years ago by Dave Bridges1.3k
2
gravatar for Charles Warden
5.9 years ago by
Charles Warden7.4k
Duarte, CA
Charles Warden7.4k wrote:

If you use log2(RPKM + rounding factor), you can use the timecourse package:

http://www.bioconductor.org/packages/2.12/bioc/html/timecourse.html

This would generally be my preferred strategy (with rounding factor = 0.1). However, it looks like you only care about changes between the 3rd and 4th timepoint. If this is the case, you only need to define two groups (pre and most) and any differential expression package could probably be used. If the sensitivity is too high (because the "pre" group is so large), you could use a combination of strategies: 1) compare 4rd versus 3rd timepoint (using DESeq, for example) and 2) filter for genes not differentially expressed between 1st to 3rd timepoints (using timecourse)

ADD COMMENTlink written 5.9 years ago by Charles Warden7.4k
0
gravatar for ashutoshdhingra
5.9 years ago by
ashutoshdhingra0 wrote:

I am in a similar situtation like you. For differentiation expression during time course control vs mutats. Can you tell me which Bioconductor package and script did you use for DE? I have neural differentiation time line of control vs mutants (2 replicates for 5 time points).

Your reply will be highly appreciared.

Best regards, Ashutosh

ADD COMMENTlink written 5.9 years ago by ashutoshdhingra0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1774 users visited in the last hour