This is a test version of Biostars. For the public version, visit https://www.biostars.org.
From WGBS-seq to DMR calling and annotation

Hi, friends, I am new in Whole Genome Bisulfite Sequencing (WGBS-seq) analysis. About Differentially Methylated Regions (DMRs) calling, I have read some papers directly using Fisher Test to detect DMRs and there are also some R packages built to do the detection. Do you have any advice on DMR calling? Also for the annotation (genes, transposable elements), any advice on carrying out the annotation? Thanks so much in advance!

sequencing r

For DMR calling, you could use (among many packages) methylKit. It has a solid documentation, is actively being maintained, authors respond to issues on Github, and has the option to aggregate CpG counts over user-specified regions in order to call DMRs. Once you have the DMRs, you could use (again among many) annotatR.

You're very welcome.

1 answer

Great! Thank you so much for your advice!

Log in to answer this question.