Hello everyone
I have a question: I did transcriptome analysis on a cassava genotype harvested at 24 hr and 72 hr, compared with non-infested healthy plants. I had 3 biological replicates for each time-point, analysed with tophat/cufflinks/cuffdiff. The 3 biological replicates RNA-seq data from the cufflinks was used to test correlation on R and on excel between the replicates, however, I get a negative correlation, because there is a high difference in the FPKM values. Can you advise as to how I can calculate correlation between the replicates?
OR any other method?
Molemi.
•
written
by
rauwaneme •
2 answers
You can try to log your values, this sometimes helps when there are some highly expressed genes that vary a lot and therefore drive the negative correlation... Can you maybe share a correlation plot? This could give us a bit more information and probably we could help better then.
Please try a different summarization units other than FPKM, e.g. TPM or CPM in addition and try a MDS (multidimensional scaling) plot. We have had many posts about the flawed concept of FPKM and this might be another effect. In my experience, samples in FPKM values cluster according to technical batches more than biological similarity. If there still is low correlation between biological replicates, you can start investigating them further.
Log in to answer this question.
-
Different results for differentially expressed genes in RNA-seq analysis with cuffdiff and R::genef…
written by strangephone1988 •
-
2d- heatmap for RNA-Seq study
written by Sujata •
-
cuffdiff "value_1" & "value_2", are they mean FPKM values from replicates for each condition?
written by vivekananthrp •
-
From cuffmerge to cuffdiff: What input for cuffdiff?
written by BioinfGuru
-
How to handle technical and biological replicates using cufflinks
written by ejmcmaho •
-
How to get an 'average' gene expression levels of bilological replicates?
written by hana
-
Cufflinks & Cuffdiff Novel Isoform Detection
written by andrew.j.skelton73
-
Large Variation In Rpkm Values Of Different Biological Replicates From Cufflinks
written by hoangtv •
-
Fpkm Values Differ Between Different Cuffdiff Runs Even When Same Data Is Used For One Of The Condi…
written by Nick
-
Tophat, Cufflinks/Cuffdiff Pipeline With Normalization
written by daniel.soronellas
Did you use rank correlation?
If I understand correctly, you've calculated Pearson's correlation coefficient between replicates and you obtained negative values. Although one generally expects replicates to have high positive correlation, this is not always the case. Low correlation between replicates can be used as a quality control criterion e.g. to remove bad/failed replicates.
Have you looked at the data i.e. plotted replicates against each other ? This would make it easier to spot if the problem is caused by outliers. You may want to use raw counts for this.