Question: How to analyze gene expression datasets with rich controls
gravatar for moxu
3.6 years ago by
moxu460 wrote:

We treat tumor cells with compound A, compound B, compound A + B mixed, and compound A-B linked respectively with different concentrations of compound(s) for each group, and we have controls (total 5 groups). We measure the gene expression for each group.

One of our goals is to find out if A-B linked & A + B mixed causes differential gene expression. If yes, what the genes are. A-B linked meaning A & B are linked together chemically and become a new compound.

A straightforward comparison is just to compare A-B linked with A + B mixed (a nested question: do a linear regression of gene expression level ~ compound concentration and compare the beta's? DESeq or edgeR might not be appropriate I am afraid). But we have experimental data for the other three groups (control, A, B). Can we make use of such datasets as well?

rna-seq R gene • 1.0k views
ADD COMMENTlink modified 3.6 years ago by zjhzwang180 • written 3.6 years ago by moxu460
gravatar for zjhzwang
3.6 years ago by
zjhzwang180 wrote:

For question 1:
I think you can use a lot of tools to do different expression analysis, such as DESeq2 ( if you data came from RNA-seq ) or limma ( if your data came from Microarray ).
For question 2:
If you want to use all experimental data, I think you can do Multivariate analysis of variance.

ADD COMMENTlink modified 3.6 years ago • written 3.6 years ago by zjhzwang180

Q1: My understanding is that DESeq or limma only does group ~ group comparison, while we have more than just two groups -- we have control (no compound), 0.1mM, 0.2mM, 0.5mM ... of compound treatments

Q2: What would be the model given the annotation in question? Like gene_expression = A + B + A*B + A-linked-B? Then the A only, B only treatments seem to be redundant.

ADD REPLYlink written 3.6 years ago by moxu460
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: 753 users visited in the last hour