Question: How to find differential expression across time points in RNA Seq
0
gravatar for RT
4.6 years ago by
RT10
Singapore
RT10 wrote:

Hello,

I have a time series RNA Seq data (time point 0 -> 48hours with step = 8hours). There are no replicates for each time point and also there is no control to check fold change against. But, somehow we want to be able to find genes that might have a differential expression over time point and use them for further analysis in GO enrichment and motif finding. Is there any tool that can be used or any paper where a similar analysis been done. 

Appreciate your suggestions. Thanks

PS: I could only think of finding the difference between max - min (normalized gene counts from) as fold - change for every gene. But I am not sure if this is correct.

-RT

ADD COMMENTlink modified 4.6 years ago by andrew.j.skelton735.8k • written 4.6 years ago by RT10
1
gravatar for andrew.j.skelton73
4.6 years ago by
London
andrew.j.skelton735.8k wrote:

Providing you're not interested in "novel" discovery - Then Alignment -> Read counting -> DESeq2 is probably the best way to go.

Alignment can be done using a variety of tools, specifically Tophat or RNAStar would be good places to start.

Read counting try HTSeq_Count or maybe feature counts (R wrapper package)

Differential expression can be done with DESeq2 in R. Check out this particular example https://support.bioconductor.org/p/63201/ Which shows a typical model for a time series experiment. 

ADD COMMENTlink written 4.6 years ago by andrew.j.skelton735.8k
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: 1327 users visited in the last hour