Question: Differential Expression and Co-expression analysis from Cell type-specific datasets
0
gravatar for volvicpellegrino
7 months ago by
volvicpellegrino0 wrote:

I use translating ribosome affinity purification (TRAP - the same principle as RiboTag) to extract translating mRNA from a specific cell type, where ribosomes are GFP-tagged. I have an RNA-seq dataset from TRAP-processed samples, that are enriched in genetic markers specific to this cell type. I want to study the effects of different conditions on this cell type through differential gene expression (DESeq2) and co-expression analysis (WGCNA).

Using DESeq2 on kallisto counts: Out of ~23,000 genes, ~4000 are significantly enriched in my TRAP RNA compared to RNA from the whole tissue input (Total RNA). ~5000 genes are significantly depleted in TRAP RNA compared to Total RNA (these include specific markers of non-target cell types).

My question: Should I filter genes before performing TRAP sample vs TRAP sample comparisons under different conditions, to remove the 5000 depleted genes?

A paper discussing analysis of TRAP datasets talks about filtering microarray probes with intensities lower than that of specific markers of non-target cell types. Would this make sense in an RNAseq scenario?

trap deseq2 wgcna ribotag • 331 views
ADD COMMENTlink modified 7 months ago by benformatics1.1k • written 7 months ago by volvicpellegrino0
0
gravatar for benformatics
7 months ago by
benformatics1.1k
ETH Zurich
benformatics1.1k wrote:

The thing is that just because those genes aren't being actively translated (i.e. present in your TRAP data) doesn't mean they aren't present/transcribed at biologically significant levels in the cells.

TRAP sounds like an enrichment method of known targets so it make sense to use the true negatives to set a threshold for expression (as mentioned for the microarray analysis). In RNA-Seq you don't have that luxury.

-

You should however filter out "unexpressed" genes (e.g. genes below a certain expression threshold in the total RNA). This is standard protocol for RNA-Seq DGE (differential gene expression).

ADD COMMENTlink written 7 months ago by benformatics1.1k

Thanks for your reply, however I don't understand why the microarray threshold can't be used in RNA-Seq?

TRAP involves immunoprecipitating GFP tagged ribosomes that are conditionally expressed in a target cell type. There is some non-specific RNA from other cell types, such as astrocytes. I could use the abundance value for a cell marker (e.g. Gfap) as a threshold to filter genes with read counts below this?

ADD REPLYlink written 7 months ago by volvicpellegrino0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 973 users visited in the last hour