Question: Cufflinks and fpkm
gravatar for qudrat
14 months ago by
qudrat70 wrote:

Hello guys, Can somebody suggest me how to choose cut off fpkm value for a list transcripts generated from Cufflinks.

rna-seq gene • 579 views
ADD COMMENTlink modified 14 months ago by Kevin Blighe42k • written 14 months ago by qudrat70

Here we go again:

You should know that the old 'Tuxedo' pipeline of Tophat(2) and Cufflinks is no longer the "advisable" tool for RNA-seq analysis. The software is deprecated/ in low maintenance and should be replaced by HISAT2, StringTie and ballgown. See this paper: Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. (If you can't get access to that publication, let me know and I'll -cough- help you.) There are also other alternatives, including alignment with STAR and bbmap, or pseudo-alignment using salmon.

ADD REPLYlink written 14 months ago by WouterDeCoster39k
gravatar for Kevin Blighe
14 months ago by
Kevin Blighe42k
Republic of Ireland
Kevin Blighe42k wrote:

There's no correct answer. Also, you should not be using Cufflinks - it has been updated to StringTie'. Moreover, FPKM is not an ideal expression unit to use in terms of comparisons across samples. If you are analysing TCGA data, you neither have to use the FPKM counts, in most cases, as the raw count HTseq files are available for download.

With that, my swift answer is to eliminate genes whose mean FPKM value are below 10. However, due to the fact that there is no cross-sample normalisation performed when deriving FPKM counts, a value of 10 means different things in different samples.

ADD COMMENTlink modified 7 months ago • written 14 months ago by Kevin Blighe42k

An update (6th October 2018):

You should abandon RPKM / FPKM. They are not ideal where cross-sample differential expression analysis is your aim; indeed, they render samples incomparable via differential expression analysis:

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.

ADD REPLYlink modified 3 months ago • written 9 months ago by Kevin Blighe42k
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: 1001 users visited in the last hour