How to calculate "fold changes" in gene expression?
3
7
Entering edit mode
7.6 years ago
juan.rosas1 ▴ 70

Hellow I have a question on how to calculate fold changes when analyzing gene expression changes between multiple tumor and control samples per gene?

R gene • 119k views
ADD COMMENT
8
Entering edit mode

Let's say there are 50 read counts in control and 100 read counts in treatment for gene A. This means gene A is expressing twice in treatment as compared to control (100 divided by 50 =2) or fold change is 2. This works well for over expressed genes as the number directly corresponds to how many times a gene is overexpressed. But when it is other way round (i.e, treatment 50, control 100), the value of fold change will be 0.5 (all underexpressed genes will have values between 0 to 1, while overexpressed genes will have values from 1 to infinity). To make this leveled, we use log2 for expressing the fold change. I.e, log2 of 2 is 1 and log2 of 0.5 is -1.

Hope this helps

ADD REPLY
0
Entering edit mode

@arnstrm what will happen if you have the same number of replicate for both control and treated ? how do you calculate the fold change?

ADD REPLY
3
Entering edit mode
7.6 years ago

Or the bioconductor limma package if you are dealing with arrays and/or RNA-Seq to analyze your data

Limma will give you the log2 expression changes based upon statistical values

ADD COMMENT
0
Entering edit mode

But how can I use it?

ADD REPLY
0
Entering edit mode

Have you read the user's guide?

ADD REPLY
1
Entering edit mode
7.6 years ago
Ido Tamir 5.2k

You could reinvent the wheel of course, but If you ask such a question, use what pros have put a lot of thought in:

http://www.bioconductor.org/packages/release/bioc/html/DESeq2.html

http://www.bioconductor.org/packages/release/bioc/html/edgeR.html

 

ADD COMMENT
4
Entering edit mode

It should also be noted that if someone has to ask how to do this sort of thing then he/she probably shouldn't be reinventing this particular wheel...

ADD REPLY
0
Entering edit mode
6.3 years ago
Vasei ▴ 30

Fold change can also be computed in unsupervised fashion, where we don't know the class labels(like case-control or type1-type2) of the samples. In that setting we can use mean expression of a gene as the base value and compute the fold change for that gene in each sample.

ADD COMMENT
0
Entering edit mode

Can you direct to any references or studies used this method. thanks

ADD REPLY

Login before adding your answer.

Traffic: 1630 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6