Hi Di,

Thank you very much for you email.

My major challenges with identifying differentially expressed genes is the
microarray data are from different platforms (Illumina and Affymetrix), and
those are the only data available for my project. In addition, my question
does not necessarily to find differentially expressed genes of each cell
type, but the *expressed genes* of each cell type are more interested. I
hope to find a way that avoids direct comparison between cell type and cell
type. I tried to rank the gene expression for each cell type and set a
cutoff for expressed and unexpressed genes, but the cutoff is arbitrary and
affects the downstream analysis. In this case, would you have any
suggestions? Any advice on obtaining differentially express ed genes for
microarray data from different platforms is also appreciated.

By the way, would you mind sending me the title of the paper that you
mentioned.

Thank you very much,
Wendy




On 4 November 2011 17:21, Wu, Di <dwu@fas.harvard.edu> wrote:

> Hi Wendy,
>
> I am not sure whether using a right gene set test is your current problem.
> It seems you want to find the signature genes for each of the cell types.
>  Therefore, for this question, it seems a differential expression problem
> to me.
>
> I understand, when you have data from several cell types, you probably
> don't have one particular cell type as a control group to all other cell
> types. I had the similar problem in the mammary gland cell type data (Lim
> 2010, Nature Medicine). What I have done is to compare the cell type A to
> each of the other three cell types, then get the overlapped up (or down)
> regulated genes in the three comparisons. These genes are the signature
> genes (expressed genes or lower-expressed genes) for the cell type A. The
> same thing can be done for the other cell types.
>
> Regarding gene set tests,  testing which pathways, GO terms or other gene
> lists are enriched in your gene list,  there are different ways. Some
> required the raw data (our "roast" and "romer" functions in limma among
> others ). The geneSetTest function in limma only used the ranks of genes.
>
> I will be happy to discuss with you more about gene set tests if that is
> actually what you face to or if you need to use them later.
>
> Hope this help,
> Di
>
>
> ----
> Di Wu
> Postdoctoral fellow
> Harvard University, Statistics Department
> Harvard Medical School
> Science Center, 1 Oxford Street, Cambridge, MA 02138-2901 USA
>
> ________________________________________
> From: bioconductor-bounces@r-project.org [
> bioconductor-bounces@r-project.org] On Behalf Of Wendy Qiao [
> wendy2.qiao@gmail.com]
> Sent: Friday, November 04, 2011 5:01 PM
> To: bioconductor@r-project.org
> Subject: [BioC] gene enrichment analysis without a control sample
>
> Hi all,
>
> I have a microarray dataset compiled from several sources, so I am facing
> some challenges with identifying the expressed genes of each cell type. I
> am thinking to use the enriched gene sets of each cell type as the
> expressed genes of that cell type. However, the gene set enrichment
> analysis (http://www.broadinstitute.org/gsea/index.jsp) needs both control
> and sample data. I am wondering if there is gene set enrichment tool for
> the analysis of one cell type only.
>
> Thank you, Wendy
>
>         [[alternative HTML version deleted]]
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor@r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives:
> http://news.gmane.org/gmane.science.biology.informatics.conductor
>

	[[alternative HTML version deleted]]

