Normalize depth before calculating TPM
1
0
Entering edit mode
colore • 0
@colore-14431
Last seen 5.0 years ago

Hi

We have RNA-seq data (three biological replicates in a single condition) and we want to calculate TPM with the sole purpose of ranking genes by expression level (in quantiles) .

I understand that, while TPM uses the total reads number to normalize sequencing depth and library size, for DE  is taken a more sophisticated approach (TMM in edgeR or "median ratio method" in DEseq2).

I have two questions:

1) Should I normalize my counts with a "DE method" instead of relying on the total reads system?

 

2) Which is the best way to collapse the three replicates? 

I was thinking of taking the median TPM for each gene, but there might be some smarter way.

 

Thank you

 

rna-seq TPM library size factor • 1.6k views
ADD COMMENT
0
Entering edit mode
chris86 ▴ 400
@chris86-8408
Last seen 3.0 years ago
UCL, United Kingdom

TPM is for comparing expression strength of different genes, this is what you want. The DESEQ normalisation method is for doing differential expression between groups, not comparing expression strength. I'd just take the mean for the replicates, but the median is reasonable.

ADD COMMENT

Login before adding your answer.

Traffic: 455 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