This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Which normalization method to use FPKM/TPM?

Hello, If I have three conditions with 3 replicates each. And If I want to check the expression of a particular Gane A across these conditions, what normalization technique would be better FPKM/TPM? And can we use these normalized counts FPKM/TPM counts for performing differential expression study, with DESeq2 / limma?

rna-seq chip-seq
• 7,782 views • link
updated by Kristoffer Vitting-Seerup Kristoffer Vitting-Seerup • written by glady glady

I agree that TPM is prefered - but generally I would say that is due to the normalization to effective transcript length (instead of annotated) and other sequncing biases.

Furthermore FPKM/RPKM are only unstrustworthy in cases where you have a global shift in the length distribution of transcripts - which in my experience is quite rare.

• 0 views • link
written by Kristoffer Vitting-Seerup Kristoffer Vitting-Seerup

could you clarify what do you mean by that?

• 0 views • link
written by Santosh Anand Santosh Anand

In fact, it is FPKM / RPKM that are rarely suitable for differential expression analysis because there is no cross-sample normalisation performed when deriving these units..

Please read this: A comprehensive evaluation of normalization methods for Illumina high-throughput RNA sequencing data analysis

The Total Count and RPKM [FPKM] normalization methods, both of which are still widely in use, are ineffective and should be definitively abandoned in the context of differential analysis.

Also, by Harold Pimental: What the FPKM? A review of RNA-Seq expression units

The first thing one should remember is that without between sample normalization (a topic for a later post), NONE of these units are comparable across experiments. This is a result of RNA-Seq being a relative measurement, not an absolute one.

• 0 views • link
written by Kevin Blighe Kevin Blighe

3 answers

If I were you, I would normalise the raw counts in DESeq2 and then compare the conditions in a pairwise fashion, and also via ANOVA / Likelihood Ratio Test. So, with 3 conditions, you would have 4 results tables. So:

  • A vs B
  • B vs C
  • A vs C
  • A vs B vs C

Doing that will already give you a p-value and log [base 2] fold change for your gene, and every other gene.

You can also compare expression levels of the gene visually by using the normalised counts or the regularised log-transformed or variance stabilised transformed counts

Kevin

• 1 views • link
written by Kevin Blighe Kevin Blighe

I didn't understand the fourth condition. How different is the fourth condition (A vs B vs C) from the above three? What different results should we except from this? How would I interpret its results? Because we might already get the list of DE genes & isoforms from the first three condition.

• 1 views • link
written by glady glady

A vs B will find genes that differ (statistically significantly) between condition A and condition B. These statistically significantly differentially expressed genes say nothing about condition C. So, it is reasonable to assume that a proportion of these genes will have equivalent expression levels in C as they do in either A or B.

A vs B vs C, however, will, generally-speaking, find those genes that differ between all conditions, i.e., it is essentially an ANOVA / Likelihood Ratio Test

• 1 views • link
written by Kevin Blighe Kevin Blighe

Okay, thank you for the explanation.

• 1 views • link
written by glady glady

TPM is probably a good way to go about getting normalized counts. However, using normalised counts as a starting point for DESeq2 is a big NO. You could optionally use the fpkm function in DESeq to get FPKM values.

• 1 views • link
written by vj vj

For miRNA samples, is there a need to consider the gene length while calculating the normalized counts like rpkm? Since, the all the miRNA have almost the same length.

• 1 views • link
written by glady glady

Perhaps not. In that case you can use RPM but if you are going to use the term "RPKM" then by definition it requires length normalisation. In this case, obviously RPKM ~ n * RPM, where "n" is just a constant, making pretty much no difference.

Vijay

• 1 views • link
written by vj vj

Please keep in mind that neither FPKM nor RPKM are suitable for differential expression analysis...

• 1 views • link
written by Kevin Blighe Kevin Blighe

For studying the differentially expressed transcripts/isoforms, can we use the raw counts. Because most of the algorithms give a strange output when isoform raw counts are used as inputs for differential expression study. Eg, DESeq2, limma.

I was planning to use tximport on the isoform raw counts and then go for a DESeq2 analysis. Is this the right way?

• 1 views • link
written by glady glady

If you want to gauge differential expression over isoforms, then EdgeR's diffSpliceDGE() would be better than DESeq2. DEXSeq is another option.

If you want differential expression over genes, then, yes, use tximport to input the raw counts to DESeq2, where they will be normalised. This is outlined in the DESeq2 Vignette

• 1 views • link
written by Kevin Blighe Kevin Blighe

Thank you for the comments.

• 1 views • link
written by glady glady

Should the sum of the RPM column be precisely 1 M, for each miRNA sample? Because now when I calculate the sum of RPM, it is around 0.5 million

• 1 views • link
written by glady glady

Yes, by definition it should be 1M ! (If you sequence a million reads, how will it be distributed among your miRNA)

• 1 views • link
written by vj vj

I divided the read counts by total no. of reads in the sample, instead of "total number of mapped reads". I guess, this was the mistake

• 1 views • link
written by glady glady

I guess by RPM, you mean RPKM.

Whatever you will do, the normalization will not yield 1M. The reason is that there is a double normalization (division) by library size and gene-length. If it has been only lib-size, they would add up to 1M. But then you also divide by gene-length and that will change (lower) the sum from 1M. And this is exactly the problem with RPKM as the sum will vary in different samples according to both the size and expression level of genes.

A rough estimate is this: You are getting 0.5 as RPM. You need to multiply everything by 2 to make 1. That means, the library mapping rate is 0.5 only, and that's a very low mapping rate for RNAseq data, and could be a problem in the data.

• 1 views • link
written by Santosh Anand Santosh Anand

You should use the raw counts for DE analysis - see this well written workflow.

• 0 views • link
written by Kristoffer Vitting-Seerup Kristoffer Vitting-Seerup

Log in to answer this question.

More posts like this