This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Distance for gene expression samples

Hi,

Imagine I have 200 samples, and for each, I have expression values of 500 genes (so it's a matrix with 500 rows and 200 columns). Expression values are normalized, it’s TPM. I need to compute a distance between the samples. What is the better way of doing it? Should I take log(TPM+1), then scale with z-transformation, and do Euclidean distance? Or there is a better way? And should I scale rows (calculate z-scores within genes but across samples) or columns?

Thanks

rna-seq distance
• 4,920 views • link
updated by Kevin Blighe Kevin Blighe • written by chipolino • chipolino

Hi, maybe try a tutorial first to understand the basics. With google I find many RNA-seq tutorials, like this one for example. Good luck.

• 0 views • link
written by Benn Benn

Well, I do know the basics. And I know that people usually calculate the distance with log CPM. However, my question is about TPM. What did I miss in the basics in your opinion?

• 0 views • link
written by chipolino • chipolino

I mean the basics to understand the difference between row or column z-scores

• 0 views • link
written by Benn Benn

Row z-score is to make sure that genes with higher expression don't have much influence on what samples are similar to each other. I just wanted to make sure that it makes sense and I didn't miss anything

• 0 views • link
written by chipolino • chipolino

Yes, this is the default way that pheatmap and heatmap.2 do it, i.e., scale by row. It is also the way that I usually scale my data to Z-scale. Assuming you have samples as columns and genes as rows, this can be done via:

t(scale(t(x)))

I show that each produces the same output in this proof, here: A: cannot replicate the pheatmap scale function

• 0 views • link
written by Kevin Blighe Kevin Blighe

1 answer

Hey, there is actually no standard way. For example, either pheatmap() or heatmap.2() (gplots) (cannot remember which) will scale your data to Z-scores, first, and then perform clustering on these [Z-scores] and produce a heatmap using the same; whereas, the other function will perform clustering on the un-scaled data and then present the Z-scaled data in the actual heatmap.

You can have complete control over your own scaling by switching off the scaling feature in both functions (or whatever other function you're using).

Your logged TPM+1 data should already be on a 'presentable' distribution. I see no issue further scaling this to Z-scores and performing both clustering and generating the heatmap on that Z-scaled data. Euclidean distance would be fine, or 1 minus Pearson correlation. Ward's linkage (ward.D2) as the linkage method will then give an evenly distributed tree layout.

Kevin

• 0 views • link
written by Kevin Blighe Kevin Blighe

Thank you for your awesome answer!

• 0 views • link
written by chipolino • chipolino

Log in to answer this question.

More posts like this