This is a test version of Biostars. For the public version, visit https://www.biostars.org.
My PCA plot looks good but i don't have differentially expressed genes

Dear All,

I have the following PCA plot from microarray data:

enter image description here

It looks quite good because the green samples and red samples, which i am interested in, are well separated along PC1. However, when i ran differential expression analysis, there are no significant differentially expressed genes and the histogram looks like this:

enter image description here

Is this result normal, or to be expected?

Thanks!

pca expression differential
• 1,994 views • link
written by druggable • druggable

You should use the raw p-values, not the adjusted p-values, to inspect the histogram

• 0 views • link
written by Papyrus Papyrus

2 answers

The PCA plot shows what seems to be a massive batch effect (left to right) and also one outlier (purple point at the bottom). The green vs red group comparison is confounded with the batch effect, so getting results will be difficult.

Can't make any comments about your DE results without knowing how you analysed the data. There are ways to account for batch effects and outliers.

• 0 views • link
written by Gordon Smyth Gordon Smyth

Your sample size is small, which may be one reason you don't detect a lot of differential expression. That being said, I would still expect to see some signal, so it may come down to how you're processing your data? It's almost impossible to tell from the information given.

Gordon's point about batch effects is also a good one -- what does each color represent in the PCA?

• 0 views • link
written by Jautis Jautis

The colors represent phenotypes: healthy (pink), obese + severe diabetes (green), obese + well controled diabetes (blue), and obese only (purple). Unfortunately, there is no information on the batch effect.

I started with raw cel files and did the following:

data = ReadAffy(celfile.path=celpath)

data.rma.norm = affy::rma(data)

rma = exprs(data.rma.norm)

Then I used prcomp for the PCA.

• 0 views • link
written by druggable • druggable

Log in to answer this question.

More posts like this