I have a question regarding the RNA-Seq analysis.

We have a time course data regarding a mouse model wt and mutant treated with a drug (10uM) and the taken down at different time points: 8 time points and did RNA-seq on these samples currently we have tpm values and raw gene counts from rna seq. Also we have only one sample for each time point. we would like to see for changes in genes across different time point but not sure statistical method should we use to determine this, should we consider for example a single time point wt vs mt, or see trends across wt and mt individually I know that there are some bioc packages for NGS analysis (edgeR, DESeq,etc). Is any package or method good for one-sample-per-time-point data?

  • $\begingroup$ For "one-sample-per-time-point data" there isn't any good statistical test available. There are tricks (such as described in edgeR manual page 21), but there is no sound statistical ways to analyze these poor designs. I wonder why you choose such a poor design? And not ask a statistician or bioinformatician first before you do the (expensive) NGS? $\endgroup$
    – benn
    Commented Aug 2, 2017 at 7:20
  • $\begingroup$ I would recommend to do at least MDA plots to see if there is visible condition / time signal. However, I agree that 1 replicate is very bad idea. Is there chance to get more replicates? $\endgroup$ Commented Aug 2, 2017 at 9:02
  • $\begingroup$ @IanSudbery another poor soul like me I guess, definitely a duplicate need to up my search skills $\endgroup$
    – sbradbio
    Commented Aug 2, 2017 at 11:48
  • $\begingroup$ There is a discussion here about how to approach this analysis using limma: support.bioconductor.org/p/33783 $\endgroup$
    – Jay Moore
    Commented Feb 15, 2019 at 23:07