hd WGCNA - Preservation statistics (Zsummary & MedianRank)

Question

[AKypR]

Hello,

I am currently using the hdWCGNA package on single cell data. So I get modules on a scRNAseq of cancer cells that I project on healthy cells in order to see what are the modules that are "absent" from the healthy cells and so "exclusive" to cancer cells. I read that when doing this kind of manipulation, it is necessary to check the "Module preservation and reproducibility" so I followed the following tutorial: https://smorabit.github.io/hdWGCNA/articles/module_preservation.html

The output results are Z-statistics and MedianRank like the ones in the tutorial. I read the paper Is My Network Module Preserved and Reproducible? in which Zsummary and MedianRank are mentioned but not the different extensions Zsummary.qual, Zsummary.pres, MedianRank.qual and MedianRank.pres. While searching on internet, I found a little on the following tutorial: https://horvath.genetics.ucla.edu/html/CoexpressionNetwork/ModulePreservation/Tutorials/cholesterolPathway.pdf I thought I understood that the .qual comes from quality and .pres from preservation.

But my questions are the following:

1. What are the differences between Zsummary.qual and Zsummary.pres?
2. How to interpret these two values?
3. Is it the same for the analysis of MedianRank.qual and MedianRank.pres?

Thank you in advance for the future answers

Answer 1

[andres.firrincieli]

The differences between .qual and .pres stats are explained in Section 4 of the supplementary material Text S1

[AKypR]

Hello andres.firrincieli,

Thank you for your answer.

So if I understand correctly, unlike the Zsummary.pres which measures the preservation of modules in the other dataset, the Zsummary.qual is the quality of the different modules themselves without taking into consideration the projection of the modules on another dataset? If I'm right, then these two statistics have nothing to do with each other, right? If I'm right, then these two statistics have nothing to do with each other and do not analyze the same thing at all, right?

[andres.firrincieli]

So if I understand correctly, unlike the Zsummary.pres which measures the preservation of modules in the other dataset, the Zsummary.qual is the quality of the different modules themselves without taking into consideration the projection of the modules on another dataset?

Pretty much that. Long story short, .qual (quality) stats tells about the reproducibility of your modules in the reference dataset, while the .pres (preservation) stats tells about the reproducibility of the reference modules in the test dataset.

then these two statistics have nothing to do with each other, right?

Not really. In a certain way .qual and .pres stats are correlated; typically, very small modules tend to be less preserved but also have very low quality. Therefore, always check the .qual stats before drawing some conclusion about the preservation of a reference module in the test dataset

[AKypR]

Thank you for this clarified answer!

[AKypR]

Hello andres.firrincieli

Sorry to reopen the discussion but a small doubt persists.

From your answer, I understand that:

If Z-summary.qual is not satisfactory then in this case we consider that the module itself is not really relevant without even taking into account the projection and therefore this module should be thrown in the trash. Isn't it?

And in the case where the Z.summary.qual is satisfactory (>10) but the Z.summary.pres is bad then should we consider that the module is relevant in its reference dataset but that the projection of the latter on the other dataset is not relevant? In this case where the Z.summary.qual is ok but not the Z.summary.pres, are you agree that we can still continue to work with the module itself without taking into account the projection?

[andres.firrincieli]

If Z-summary.qual is not satisfactory then in this case we consider that the module itself is not really relevant without even taking into account the projection and therefore this module should be thrown in the trash. Isn't it?

That is correct

And in the case where the Z.summary.qual is satisfactory (>10) but the Z.summary.pres is bad then should we consider that the module is relevant in its reference dataset but that the projection of the latter on the other dataset is not relevant?

That is correct

In this case where the Z.summary.qual is ok but not the Z.summary.pres, are you agree that we can still continue to work with the module itself without taking into account the projection?

I agree.

[AKypR]

Thank you a lot!