Question: How to use TPM from RNA seq data analysis for differential gene expression analysis? which statistical methods are reuired to be performed.
gravatar for amnahsiddiqa
24 months ago by
amnahsiddiqa0 wrote:

i am a newbie in RNA seq data analysis. I want to check a gene as DEG in a dataset of RNA-chip seq experiment. The data provided is in the form of a single column for each treatment type and lists the expression level of each gene normalized to transcripts per million (TPM).i need help on using what type of analysis using R to find out the DEGs.

The snapshot of data is like this.

Gene EntrezID Normal_TPM Diabetes_TPM
Apoa1 11806 14668.15 2875.06
Mup3 17842 9992.58 1697.63
Serpina3k 20714 8031.3 2849.67
Alb 11657 6801.26 6912.08
Cyp2e1 13106 6580.8 7816.79

Thanks in advance

rna-seq • 3.8k views
ADD COMMENTlink modified 24 months ago by Govardhan Anande100 • written 24 months ago by amnahsiddiqa0
gravatar for Govardhan Anande
24 months ago by
Govardhan Anande100 wrote:

To analyse differential expression analysis of genes in R, you can use DESeq, DESeq2 or edgeR.

Formula for TPM is here, so if you can get total reads aligned for each sample then you can find out aligned reads freq, which you can use as input for above programs and can perform differential expression analysis.

According to your snapshot, it looks like your data is already analysed for differential expression.

Correct me if I am wrong.


ADD COMMENTlink modified 24 months ago • written 24 months ago by Govardhan Anande100

edgeR works with raw counts, so maybe EBseq could be a better advice for TPMs?

ADD REPLYlink written 24 months ago by cpcantalapiedra140

Hi Govardhan, I think you're mixing up CPM (counts per million) with TPM (transcripts per million). See here how it's computed.

ADD REPLYlink written 24 months ago by michael.ante2.4k

yeah, so you can get TPM formula here then.

ADD REPLYlink written 24 months ago by Govardhan Anande100
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.3.0
Traffic: 1172 users visited in the last hour