SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
Q values in cuffdiff with 2 biological replicates nr23 Bioinformatics 1 02-20-2014 06:47 PM
Different RPKM values in same dataset using Cufflinks or Cuffdiff (v1.3.0) hlwright Bioinformatics 3 11-09-2012 03:51 AM
Differential gene expression: Can Cufflinks/Cuffcompare handle biological replicates? marcora Bioinformatics 38 12-14-2010 03:57 PM
cufflinks problem: one gene have multiple rpkm values calliopsis Bioinformatics 1 05-25-2010 08:18 AM
Differential gene expression: Can Cufflinks/Cuffcompare handle biological replicates? marcora Bioinformatics 0 05-19-2010 01:11 AM

Reply
 
Thread Tools
Old 02-11-2014, 11:56 AM   #1
thanhhoang
Member
 
Location: Ohio, USA

Join Date: Jul 2013
Posts: 16
Question Variation in RPKM values of different biological replicates produced by Cufflinks

Hi all,
I want to rank gene expression based on the transcript abundance to see the relative expressions of different genes in my RNA-seq sample. I ran Cufflinks for 3 biological replicates for that condition to get the FPKM values. However, the FPKM values are very different within the 3 biological replicates, which I am not sure why and how to deal with this.
In Cufflinks manual, It said that Cufflinks can only use classic-fpkm normalization on one library at a time. Should I use Cuffdiff2 to normalize FPKMs and fragment counts across all libraries (This policy is identical to the one used in DESeq) and take the average of these normalized RPKMs?. Or do you think it is ok to get the FPKMs from biological replicates produced by Cufflinks and feed them into DESeq to normalize them?

Thanks for your opinions.
Thanh
thanhhoang is offline   Reply With Quote
Old 02-12-2014, 05:30 AM   #2
mbblack
Senior Member
 
Location: Research Triangle Park, NC

Join Date: Aug 2009
Posts: 245
Default

The FPKM values would not be expected to be the same or even necessarily very similar amongst biological replicates. That is the point of using biological replicates - large differences in relative transcript abundance are often due (at least in part) to innate biological variation in expression.

For some conditions/systems/instances the actual individual variation in expression may be the single largest source of variation in expression amongst the samples. For example, in human exposure to arsenic, gene expression estimates show the intra-individual variance is far greater than the inter-individual variance. PCA plots of gene expression estimates for individual humans will not show any overlap even out to 3 or 4 standard deviations - each individual's expression is largely unique (at least for some tissues).

That's the problem with trying to rank or compare individual expression estimates - you have single point estimates of expression per transcript, and those are whatever they are.
__________________
Michael Black, Ph.D.
ScitoVation LLC. RTP, N.C.

Last edited by mbblack; 02-12-2014 at 05:34 AM.
mbblack is offline   Reply With Quote
Reply

Tags
cufflinks, deseq, fpkm, rna-seq

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 06:45 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO