[BioC] (control vs treatment) in (wildtype vs mutant) microarray analysis
Thomas.H.Hampton at Dartmouth.edu
Mon May 25 17:39:22 CEST 2009
You make a great point about the problem of comparing
lists of differentially expressed genes (DEGs) under various conditions.
The lists overlap to a certain extent, but it is hard to get a
feeling for what these
overlaps mean. Certainly, the fact that a gene fails to meet a
certain level of
significance in some test does not imply that it is not
differentially expressed under
your experimental conditions.
My suggestion is that you visualize various gene sets (say, genes
affected by toxin) across all your treatments at once, A heatmap
(particularly if your gene sets are around 50-100 genes) is generally
very informative. Typically, you will find that mutant status alters
effect in interesting ways depending on the gene: some expression
will be amplified, some curtailed, etc. A consideration of the genes
affected in this
way, that is, the roles they play biologically, should yield insight.
In fact, you could begin with a heatmap of genes most affected by your
conditions in general -- say the 200 genes with most significant by
On May 25, 2009, at 11:05 AM, Cheng-Yuan Kao wrote:
> We have done affy microarrays for wildtype-control treatment,
> wildtype-toxin treatment, mutant-control treatment and mutant-toxin
> The goal is to find diffferentially expressed genes regulated by toxin
> in wildtype and then find out which of these regulation are mutant
> The first goal is typical. So we did R/ bioconductor - SAM and limma.
> Both could give us a bunch of DEGs. However, I am lost about getting
> the second aim done.
> With limma and 2x2 factorial analysis, we could find the DEGs from all
> kind of pairs, such as wildtype -Toxin/control (this answers the first
> goal) or mutant/wildtype in control treatment (this tells us how the
> mutant gene is affecting the basal expression without toxin).
> But I don't know how to find the wildtype DEGs which have regulation
> depending on the mutant gene. Say one gene is up-regulated 100 folds
> by toxin (i.e. toxin treatment/control treatment) in wildtype. Then if
> this gene is up-regulated 3 folds (toxin treatment/control treatment)
> in mutant, it will be apparently mutant gene dependent (from 100 folds
> in wildtype to 3 folds in mutant). However, this gene will be shown as
> DEG in the mutant analysis as well since it is more than the 2-fold
> cutoff. Then if I only compared the DEGs in different pair of analysis
> from limma, I will miss these kind of genes. There must be some
> dedicated way to analyze this but I could not find it. Any suggestion
> will be appreciated. Thanks a lot.
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> Search the archives: http://news.gmane.org/
More information about the Bioconductor