0
$\begingroup$

I have case control studies,4 different condition where I have 1 control and against which I'm have 3 comparisons to make.

I did WGCNA for individuals samples .Here I have the kept the number of genes equal in all the cases .

Now my question is whether the number of modules which are obtained from WGCNA would be similar in number if I run module preservation between two condition .

Lets say i get 15 modules in Control condition and 16 treated condition if I run module preservation the number of modules which would come would be 15+16= 31 or its independent .

The power I have kept same for both WGCNA and module preservation analysis. But Still I get different module numbers.

Any suggestion would be highly appreciated

Improve this question
$\endgroup$

1 Answer 1

Reset to default
2
$\begingroup$

Your question is not clear, you seem to mix the relevant terminology...

First, you cannot run WGCNA for individual samples, only for individual conditions (with or without WT, depending on what question you're trying to answer and assuming you have enough samples in each condition). If that's what you have done, fine. The number of modules in each independent analysis will generally be different, that's because each condition is different. The module identification built into WGCNA as the default approach does not allow the user to specify the number of modules beforehand, but if you insist on getting a pre-specified number of modules, you could run WGCNA using the "manual" way (for example, Tutorial I, section 2b) and substitute another clustering method (e.g., PAM) that allows you to specify the number of clusters. You may find the results inferior though.

Second, module preservation does not result in new modules. Module preservation evaluates whether modules found in data set A are preserved in a different data set B. To get new, common modules, you could run a consensus network analysis on pairs of conditions or all 3 conditions together. The number of modules you would find in a consensus analysis cannot be predicted from the number of modules you find in individual data sets.

Improve this answer
$\endgroup$
4
  • $\begingroup$ "assuming you have enough samples in each condition)." sorry yes i mixed up and yes i have control have like 55 samples and treatment like 50 samples .Thank you for the clarification now it seems bit more clear I will read the tutorial a bit more $\endgroup$
    – kcm
    Commented Oct 18, 2019 at 19:36
  • $\begingroup$ Conceptual clarification I have done as I said Control sample WGCNA as well as treatment ones with numerous biological replicate .I get various modules and yes "umber of modules in each independent analysis will generally be different, that's because each condition is different." yes this is what i get .So if do module preservation analysis where I want to see if my "data set A are preserved in a different data set B" then the number of resulted module would be same as what I did WGCNA for Control sample where I would see w.r.t to my reference. Am I explaining it correctly ? $\endgroup$
    – kcm
    Commented Oct 19, 2019 at 9:31
  • 1
    $\begingroup$ Sorry, I still don't understand your question. Module preservation outputs, for each module in a reference network (A), preservation statistics of the module in a test network (B). Module preservation does not change the number of modules. Some modules from A may be preserved in B and some not; sometimes all modules are preserved, sometimes none of the modules are preserved. The number of modules from A that are preserved in B may be smaller or larger than the number of modules found in a separate network analysis of B. $\endgroup$
    – Peter Langfelder
    Commented Oct 20, 2019 at 18:17
  • $\begingroup$ Thank you for the clarification $\endgroup$
    – kcm
    Commented Oct 22, 2019 at 19:44

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.