Hi,
I am trying to write a script in python to analyse my qPCR data, however I get confused with terminology how to report and plot results. I'd appreciate some advice on how to do this properly
What is the best way to report the results? In my lab they usually report mRNA expression of the Gene of Interest (GOI) relative to the reference gene, but the numbers you get seem quite irrelevant from my point of view; if you plot them, however, you clearly see the differences. I am also aware of the fold change, which seems more relevant if you just report the value, but when you plot them I think there's the problem that it's skewed to upregulated genes (which can go from one to infinity, whereas the downregulated ones can just be from 0 to 1).
To solve the latter, I have seen recommendations to plot log2 fold change, but that seems weird, because isn't that just the same value as ddCt?
If you do recommend reporting the fold change, what kind of error bars should I use when making the bar plot? Are the distributions still normal?
Thank you very much!
Hi,
What exactly is the purpose of the plot you are trying to make? Show the raw Ct values? To show the normalized (compared to ref gene) values? The fold change values between two conditions?
All options are possible, it depends on what your purpose is of the plot.
So my final purpose is to study differences in expression of a gene of interest in treated vs control (or disease vs healthy). For that, yes, I want to have the values normalised to the reference gene. So I wonder what is the best way to show that, I guess it's partly a matter of opinions, but I'd like to make an informed decision :)
In that case I would choose for log2 transform fold change values, especially when both up- and down regulated genes are put in one plot. Log2 transformation makes up- and down regulation (ratios) symmetrical.
Aha, thanks. However, isn't the log 2 of the fold change equal to ddCt? Also, in the case of plotting that, what kind of error bars do you use?
Could you please put an example of such analysis and display of data?
Thanks a lot!
I think ddCt is similar to log2 FC.
You can check this package in R for further explanation:
https://www.bioconductor.org/packages/devel/bioc/vignettes/ddCt/inst/doc/rtPCR.pdf