How to find differential expression across time points in RNA Seq
Entering edit mode
7.3 years ago
RT ▴ 10


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.


differential-expression RNA-Seq time-series • 4.3k views
Entering edit mode
7.3 years ago

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 Which shows a typical model for a time series experiment.


Login before adding your answer.

Traffic: 2197 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6