Question: arrayExpress differential gene express analysis
gravatar for jadepinket
9 months ago by
jadepinket0 wrote:

I have a dataset that looks like below:

Gene ID Gene Name p-value1 log2foldchange1 p-value2 log2foldchange2 p-value3
ACEGIKM00000000001 AABR07013255.1 NA 0 NA 0 NA
ACEGIKM00000000007 Gad1 NA 0 NA 0 NA
ACEGIKM00000000008 Alx4 NA 0 NA 0 NA
ACEGIKM00000000009 Tmco5b NA 0 NA 0 NA
ACEGIKM00000000010 Cbln1 NA 0 NA 0 NA
ACEGIKM00000000012 Tcf15 NA 0 NA 0 NA
ACEGIKM00000000017 Steap1 0.293657 0.2 0.176462 0.2 0.08213
ACEGIKM00000000021 AABR07061902.1 0.058899 -0.3 0.919169 0 0.95051
ACEGIKM00000000024 Hebp1 0.904233 0 0.589132 0.1 0.637529
ACEGIKM00000000033 Tmcc2 NA 0 NA -0.1 NA
ACEGIKM00000000034 Nuak2 0.580938 -0.1 0.882088 0



I want to Identify differentially expressed genes (for example, using p-value or fold change or both) for each treatment that includes direction of change, then I want to Identify the most important pathways impacted by the treatment. I also want to do visualizations that show the changes in gene expression as a function of treatment.

I also have a counts file, but since I already have this file with fold change and pvalues, I was hoping I could get my answers from this.

Any suggestions are appreciated.


ADD COMMENTlink modified 9 months ago by Aaron Lun18k • written 9 months ago by jadepinket0
gravatar for Aaron Lun
9 months ago by
Aaron Lun18k
Cambridge, United Kingdom
Aaron Lun18k wrote:

If the p-values are already BH-adjusted, you can just define DE genes as those with p-values below a desired threshold, e.g., 5%. If not, you'll have to adjust them yourself using p.adjust with method="BH". After you identify DE genes, you can identify enriched pathways/processes using goana or kegga from the limma package. Check whether your organism is supported, though; the IDs aren't familiar to me.

What you mean by "visualizations" is vague. If the log-fold changes describe the effect of treatment, you can just plot them. They look a bit weird, though, usually there's more than 1 decimal point.

ADD COMMENTlink modified 9 months ago • written 9 months ago by Aaron Lun18k
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.2.0
Traffic: 190 users visited in the last hour