Question: How do I find out regions of methylation enrichment using MEDIPS bioconductor package?
gravatar for lakhujanivijay
4.3 years ago by
lakhujanivijay5.4k wrote:

I have a single MEDIP-seq sample and I am only interested in finding out regions of methylation. How do I go about it? Since, I am not interested in differentially methylated regions (DMR's) as I don't have multiple samples, how could I achieve the same using MEDIPS?


I know there are a lot of tools available but this one is highly cited and used. There are other tools which in turn use MEDIPS at the backend (example - MeQA).

What has been already done?

Since I am working with a genome which is not already present in the MEDIPS custom genome package (BSgenome), I was able to create one for myself based on below guidelines and forge it into BSgenome.

Here is the Forging_custom_genome_tutorial

Now, I am looking at this tutorial.

The problem:

The methodology has been beautifully explained for the demo dataset called "MeDIPSdata" provided within MEDIPS package (which consists of a couple of BAM files),however, it is not clear to me as to how to proceed with a single BAM file (since I have a single sample)?

  1. Which steps are important and cannot be skipped (QC, normalization)?
  2. How to find regions of methylation for a single sample?

PS: Respond only if you can help me with regards to MEDIPS package. No other package information would be useful for me.


methylation medips medip-seq • 1.2k views
ADD COMMENTlink modified 4.3 years ago • written 4.3 years ago by lakhujanivijay5.4k
Please log in to add an answer.


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