This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Interpreting Differential Methylation Values

I'm having some trouble interpreting some results here on Differential Methylation that I'm hoping someone can clear up.

I've run some biseq results through an R package called Methylsig. The input file looks something like this:

chrBase chr base strand coverage freqC freqT
chr21.43008527 chr21 43008527 F 32 100.00 0.00
chr21.43008531 chr21 43008531 F 32 96.88 3.12
chr21.43008543 chr21 43008543 F 32 90.62 9.38
chr21.43008674 chr21 43008674 R 27 100.00 0.00
chr21.43008710 chr21 43008710 R 67 94.03 5.97
chr21.43008720 chr21 43008720 R 67 92.54 7.46

The result is a data table like this:

methylSigDiff object with 3,260 rows
--------------------------
chr start end strand pvalue qvalue meth.diff logLikRatio
1 chr21 43008527 43008527 + 0.3219880 0.6664426 -2.4167475 1.27484057
2 chr21 43008531 43008531 + 0.8718018 1.0000000 0.4158129 0.02957797
3 chr21 43008543 43008543 + 0.5583045 0.8746145 -2.8071020 0.40674783
4 chr21 43014041 43014041 + 0.3670413 0.7056865 2.3585424 0.98274395
5 chr21 43014044 43014044 + 0.9166589 1.0000000 -0.6133413 0.01210852

My problem here is understanding what exactly "meth.diff" is. I get that -2.4167475 showing less methylation vs the control and 2.3585424 is showing increased methylation vs the control. But I'm not sure exactly how to interpret that.. Ideally I would like to be able to use the Diff.Methylation value to calculate FreqC and FreqT backwards.

If anyone could simplify this I would greatly appreciate it.

sequencing methylation bi-seq

My guess is that meth.diff is the percent change. So -2.4167475 would be a -2.4% methylation change. Just check one position in all of the samples and that should let you see if such a value makes sense.

0 answers

No answers yet.

Log in to answer this question.