Rna Seq Cell Type Specificity
2
2
Entering edit mode
11.6 years ago
k.nirmalraman ★ 1.1k

Hello,

I now have RNA-Seq data from different cell types and ideally the idea was to perform DE. However, I was wondering if it would be possible to identify cell type specific marker genes for these cell types.

What do you think is the ideal way to approach this question.

I am using DESeq for differential analysis.

Thanks in advance!

rna rna-seq deseq bioconductor • 5.7k views
ADD COMMENT
0
Entering edit mode

I don't really understand what is the question here. The answer is in your question. Performing DE analysis will give you cell type specific markers for the cell type you compare.

ADD REPLY
0
Entering edit mode

I have atleast 5 different cell types and wanted to know if there is any other kind of approach/better approach can one take other than DESEq

ADD REPLY
1
Entering edit mode

The definition of differential expression is pretty clear. However, the definition of cell-type-specific is not. If you can come up with a clear definition of that concept for your system, it might be easier to answer your question.

ADD REPLY
1
Entering edit mode

Cell Type Specific: Genes that are predominantly expressed in a particular cell type and not in other (of course when performing DE Analysis, one would get this contrast with respect to two cell types.) But in this case, when one has many cell types, my approach would be to contrast all combinations.

My question is to find out if there are any direct method that is available to do this.

In other words, to identify a set of features that explicitly identify a cell type ie a sample from a given count data.... like clustering! Hope this helps to understand!

ADD REPLY
0
Entering edit mode

See my answer below.

ADD REPLY
0
Entering edit mode

Hi,

I am looking to solve similar problem. I have been referring to the CellCODE enter link description here package but haven't been able to utilize it for my dataset.

Please let me know if you were able to figure it out and what approach did you use?

Thanks.

ADD REPLY
2
Entering edit mode
11.6 years ago
David ▴ 740

There are more than one way to go for it but the first step is to extract the 5 DE gene lists from your 5 cell types.

As a quick test I would do a Venn diagram. If you are lucky you will obtain genes that you declared "DE" and that are specific for the cell types.

After that you can try a PCA, for example, and looks at the gene ordering in the component capturing the most variance citation

This is one way to go for it.

ADD COMMENT
1
Entering edit mode
11.6 years ago

Assuming that what you are looking for is a set of genes that is differentially expressed in only one condition compared to all others, you could:

  1. Apply voom from the limma package to transform your counts
  2. Perform DE testing in the standard way using limma
  3. Use the decideTests function in limma, which results in each gene being classified as up-or-down regulated while "conditioning" on other tests. There are some parameters to control how this is done, so you'll want to read the docs.

OR

  1. Perform DE testing in edgeR
  2. Use decideTestsDGE analogously to #3 above.
ADD COMMENT

Login before adding your answer.

Traffic: 1543 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6