Question: miRNA seq time course analysis tool with 2 phenotypes
gravatar for pier-luc.dudemaine
4.5 years ago by
pier-luc.dudemaine0 wrote:

Hi, I'd like to know what would be the best tool to analyze differentially expressed genes in an infection time course assay using 2 phenotypes.

I'm explaining, we did grow macrophages from 2 differents groups, healthy and sick animals (8 animals each phenotypes). Afterward, we did an infection time course (uninfected control, infection 1h, infection 4h, infection 8h, infection 12h and infection 24h).

Now I would like to visualize if the response from the cells of healthy animals differs from the one of sick animals. What would be the best way to do it? In all the case I want to keep the comparison of the infections time points with the uninfected controls, cause I already know that my controls are differents from animal to animal.

I used Deseq2 to compare each time points to my control, but I'd like to use a model where the phenotype is taken into consideration.

Any suggestions?

rna-seq analysis • 1.0k views
ADD COMMENTlink modified 4.5 years ago by Devon Ryan98k • written 4.5 years ago by pier-luc.dudemaine0
gravatar for Devon Ryan
4.5 years ago by
Devon Ryan98k
Freiburg, Germany
Devon Ryan98k wrote:

DESeq2, edgeR or limma/voom can all do that, take your pick.

ADD COMMENTlink written 4.5 years ago by Devon Ryan98k

@Devon Ryan I think DESeq2 is more suitable than other packages because the data is count and DESeq2 is built for such type of data

ADD REPLYlink written 4.5 years ago by Mo920

All three packages are meant for count data (voom is the limma method for preparing count data).

ADD REPLYlink written 4.5 years ago by Devon Ryan98k

By using the multifactor design? (I was not so sure it was suitable for that) Or you would use Interactions which seems to me to be closer to what I want, but still not sure it is the best.

The best way I can see would be to do my analysis as a Time-series experiment but I can't find how to do it when having two conditions (my two groups infected vs healthy) and comparing them as the time-series goes.

What's your thoughts?

A big thank for helping

ADD REPLYlink written 4.5 years ago by pier-luc.dudemaine0

Something of the form ~time * condition could work. You probably want time as something continuous, rather than as a factor, and you might want to modify it for some hypothesis or another. There are some examples around if you search "DESeq2 spline" or "DESeq2 time series". In short, you first need to define what you would like "comparing them as a time series" to mean, then the rest sorts itself out.

As an aside, it's best to cluster first, since that'll give you a clue about how the times should be handled.

ADD REPLYlink modified 4.5 years ago • written 4.5 years ago by Devon Ryan98k
Please log in to add an answer.


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