Question: WGCNA on methylation data
gravatar for enora.fremy
4 months ago by
enora.fremy0 wrote:

Hello everyone,

I am working on methylation data and I want to use WGCNA on them. The data are the beta values of 400K CpG probes and ~30 samples which were transformed by log2(beta+1) (The data then have a nice scale-free topology profile when the soft threshold is calculated, I think it is a good sign? ). So, if someone already use WGCNA for this kind of data, I would have some questions: Firstly, are my data properly normalized? And secondly, how can I choose the maximal number of blocks for the blockWiseModule function ? I have access to a cluster, can I deduct the maxBlockSize when I know the available RAM?

Thanks in advance,

Have a nice day !


ADD COMMENTlink modified 4 months ago by Peter Langfelder2.2k • written 4 months ago by enora.fremy0
Answer: WGCNA on methylation data
gravatar for Peter Langfelder
4 months ago by
United States
Peter Langfelder2.2k wrote:

I have done a few WGCNA analyses on Illumina 450k methylation data. I used beta values and filtered the data down to about 300k most variant probes. I haven't thought deeply about which transformation is best for methylation data but using beta values worked fairly well. I would not use a log2(beta +1) transformation since it will exaggerate differences in CpGs with low methylation; this transformation is useful for count data.

maxBlockSize indeed depends mainly on available RAM. I would say you should be able to use maxBlockSize around 30k if you have 64GB of RAM, perhaps even 32GB. Increase the maxBlockSize by a factor of 1.4 for every doubling of RAM; with 256GB you should be able to use block size of ~60k, 1TB would let you use ~120k etc. You can play with the block sizes to some degree; in my experience the actual maximum block size usable on a system depends also on how the system is configured, how much swap space there is etc.


ADD COMMENTlink written 4 months ago by Peter Langfelder2.2k

Thanks a lot for your quick reply and your great job on WGCNA !

After several tests, I used a maxBlockSize of 60K with 128GB of RAM with the beta values. It worked out fine :)


ADD REPLYlink written 4 months ago by enora.fremy0
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 342 users visited in the last hour