Question: Do I have to use beta-values/m-values for RRBS/WGBS data ?
0
gravatar for Tori
3.7 years ago by
Tori90
Finland
Tori90 wrote:

I am working on RRBS data. Someone suggested me to use M-values for the calculation of differential methylation. In my understanding beta-values/m-values are used to measure methylation level in array based method. In sequencing data we get single base pair resolution methylation level after methylation call. Hence, no need for calculating beta-values/m-values.

I would like to make sure that I am right. Could you please confirm that beta-value/m-value calculation is not needed for sequencing data?

wgbs rrbs methylation • 2.1k views
ADD COMMENTlink modified 3.7 years ago by Devon Ryan97k • written 3.7 years ago by Tori90
1
gravatar for Devon Ryan
3.7 years ago by
Devon Ryan97k
Freiburg, Germany
Devon Ryan97k wrote:

Every tool will have its own input requirements. Some tools will require the raw counts, other will require methylation percentages in the range 0-100 (or 0-1, it'll depend on the tool). You shouldn't be writing a tool to do this yourself, use something published and give it whatever it wants as input.

ADD COMMENTlink written 3.7 years ago by Devon Ryan97k

I have used Methylkit, RnBeads and RadMeth for analysing RRBS data, but after diving into greater detail of the result files of these methods DM regions with high p-value/qvalue does not make sense. Those statistically significant (high p-value/qvalue) DM regions do not seem significant in biological sense. Therefore, I would like to use Bismark output methylation call files in some other statistical test.

ADD REPLYlink written 3.7 years ago by Tori90
1

I've had some nice experiences with metilene, give that a try. It may be that you simply have no biologically relevant DMRs.

ADD REPLYlink written 3.7 years ago by Devon Ryan97k

Hi Devon, I tried metilene and I get thousands DMR with different length. But how can I get the DMR relation to gene. I tried bedtool, but get nothing.Maybe there is something wrong. what is you method to work this out?

ADD REPLYlink modified 3.2 years ago • written 3.2 years ago by t2g4free0

bedtools nearest should work fine.

ADD REPLYlink written 3.2 years ago by Devon Ryan97k

And I see some Question you answered before, you said not to use Fisher test. when there are two sample and no replicate, which way or software is better to get DMR? thx

ADD REPLYlink written 3.2 years ago by t2g4free0

What I was likely trying to convey is that such a design is inherently problematic. Your only option is something like a Fisher's test, though.

ADD REPLYlink written 3.2 years ago by Devon Ryan97k

How to get some result in WGBS? Hope get your answer, thx. And bison is good aligner, and some auxilary scripts is wonderful.

ADD REPLYlink written 3.2 years ago by t2g4free0
Please log in to add an answer.

Help
Access

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