methylKit icon indicating copy to clipboard operation
methylKit copied to clipboard

MethylationStats shows one peak

Open gevro opened this issue 11 months ago • 4 comments

Hi, My methylationstats profile doesn't look like what is in the vignette. Is this a computational issue or does this suggest a specific sample issue?

Thanks Screenshot 2024-03-13 at 8 29 01 PM

gevro avatar Mar 14 '24 03:03 gevro

Hi @gevro,

if there were no errors when importing your data I assume it is a sample-specific issue. Could you provide further details about the data source and the sample?

You should probably also consult the person who generated the data, to see if this is expected.

Best Alex

alexg9010 avatar Mar 14 '24 10:03 alexg9010

It is a bisulfite RRBS WGS sample. Is this simply an indication of low conversion efficiency?

gevro avatar Mar 14 '24 13:03 gevro

Hi @gevro

The methylation percentage pattern for RRBS is different than the one for WGBS. RRBS targets sites with high CpG density, e.g. CpG islands, which are usually highly methylated. WGBS covers the whole genome including areas with low CpG density, which are usually lowly methylated.

See https://www.illumina.com/science/sequencing-method-explorer/kits-and-arrays/rrbs-seq-scrrbs.html

Best Alex

alexg9010 avatar Mar 15 '24 10:03 alexg9010

Sorry, I meant to write it is WGBS. Is this pattern abnormal for WGBS?

gevro avatar Mar 15 '24 14:03 gevro

Hi @gevro,

This depends on the organism that you are studying. For mammals, the normal pattern would be a bimodal distribution for CpG and very low methylation for non-CpG sites, like in the vignette. On the other hand, for plants, you would expect a different pattern where also non-CpG sites are methylated (see here).

Best Alex

alexg9010 avatar Mar 18 '24 13:03 alexg9010

Thanks!

gevro avatar Mar 18 '24 13:03 gevro