MeDIP-seq and Seqmonk

Grenouille · 06-05-2013, 03:50 AM

Hi everyone,

I am looking for some advice regarding the analysis of a MeDIP-seq dataset. I generated paired-end sequences of two biological replicates of the same genotype, while just one Input sequence is available.
I was planning to use Seqmonk to have a first glance but despite it's very user-friendly I am having a hard time understanding which tools I should use; in particular, I do not see at which step the Input dataset is taken into account to correct the IP datasets. The only option that makes the distinction between IP and Input seems to be the use of MACS to define probes but I was wondering if the MACS algorithm is good at finding what I expect to be rather wide peaks.
I hope my queston is not too confused.
Thanks in advance.

turnersd · 06-05-2013, 05:12 AM

Haven't used SeqMonk. Take a look at MEDIPS. What are others using for MeDIP-Seq analysis with biological replicates?

Grenouille · 06-05-2013, 05:31 AM

Thank you for the reply. I am already using MEDIPS but just to compare different samples, apparently it is not appropriate for the identification of sample-wise local enrichment. I was told that to look for methylation-rich regions it is better to use a regular peak finder tool.

turnersd · 06-05-2013, 05:39 AM

I've used MACS and ZINBA on MEDIPS data, but this data doesn't really look "peaky" the way other xIP experiments might, and not sure how methods like these that were mainly designed for ChIP-seq are performing on MeDIP.

Grenouille · 06-05-2013, 05:59 AM

That's exactly my point. I would say that MeDIP peaks shouldn't look too different from ChIP peaks for histone modifications (such as H3K27me3) but as I am a beginner I wanted to ask for the opinion of more experienced people!

paolo.kunder · 06-05-2013, 07:06 AM

Hi,
I am also using MEDIPS for sample-wise enrichment, from the manual:

Quote:

In principle, it is possible to calculate genome wide coverage and methylation
proles for only one MEDIPS SET or for only one group of MEDIPS SETs
using the function
MEDIPS.meth

But in ther meantime I used also MACS to calls for peaks refining the bandwitch and results are consistent. Then I called peaks with some other programs and extracted a consensus peakste with R package DiffBind

Anyway, I am really interesting in sharing knowledge with you guys..

Grenouille · 06-10-2013, 08:42 AM

Hi Paolo,

thanks for replying.
Did you include Input samples in your analysis? I know that the usefulness of such data is controversial but in my case I think I am really obliged to take Inputs into account; I am working with several different genotypes, while only one reference genome is available. That means that, due to structural variants, for some samples I could have some regions that show no coverage because the sequence is absent in that particular genotype, and here is where the Input becomes essential.
Are you working with similar data?

tanaybhatt · 04-25-2017, 03:58 AM

Can you direct me to a published paper where DMR analysis has been done without taking the input in consideration (or without input subtraction/normalization)?

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
meDIP-chip/meDIP-seq/MethylCap-seq using mouse DNA	Jones90	Epigenetics	6	09-13-2012 11:00 PM
Chip-seq or Medip-seq not for C57BL/6 but Balb/c	gatapishi	General	1	09-06-2012 10:54 PM
SeqMonk: Export features (e.g. CpG Islands) from Ensembl for import into SeqMonk?	jjw14	Bioinformatics	8	06-07-2012 11:35 PM
WGA in MeDip-Seq	mugiu17	Epigenetics	3	10-12-2011 05:13 AM
MeDIP-seq versus ChIP-seq analyses?	jjw14	Bioinformatics	8	08-24-2010 08:22 AM

06-05-2013, 03:50 AM	#1
Grenouille Junior Member Location: France Join Date: Dec 2012 Posts: 4	MeDIP-seq and Seqmonk Hi everyone, I am looking for some advice regarding the analysis of a MeDIP-seq dataset. I generated paired-end sequences of two biological replicates of the same genotype, while just one Input sequence is available. I was planning to use Seqmonk to have a first glance but despite it's very user-friendly I am having a hard time understanding which tools I should use; in particular, I do not see at which step the Input dataset is taken into account to correct the IP datasets. The only option that makes the distinction between IP and Input seems to be the use of MACS to define probes but I was wondering if the MACS algorithm is good at finding what I expect to be rather wide peaks. I hope my queston is not too confused. Thanks in advance. Last edited by Grenouille; 06-05-2013 at 03:58 AM.

06-05-2013, 05:12 AM	#2
turnersd Senior Member Location: Charlottesville, VA Join Date: May 2011 Posts: 115	Haven't used SeqMonk. Take a look at MEDIPS. What are others using for MeDIP-Seq analysis with biological replicates?

06-05-2013, 05:31 AM	#3
Grenouille Junior Member Location: France Join Date: Dec 2012 Posts: 4	Thank you for the reply. I am already using MEDIPS but just to compare different samples, apparently it is not appropriate for the identification of sample-wise local enrichment. I was told that to look for methylation-rich regions it is better to use a regular peak finder tool.

06-05-2013, 05:39 AM	#4
turnersd Senior Member Location: Charlottesville, VA Join Date: May 2011 Posts: 115	I've used MACS and ZINBA on MEDIPS data, but this data doesn't really look "peaky" the way other xIP experiments might, and not sure how methods like these that were mainly designed for ChIP-seq are performing on MeDIP.

06-05-2013, 05:59 AM	#5
Grenouille Junior Member Location: France Join Date: Dec 2012 Posts: 4	That's exactly my point. I would say that MeDIP peaks shouldn't look too different from ChIP peaks for histone modifications (such as H3K27me3) but as I am a beginner I wanted to ask for the opinion of more experienced people!

06-10-2013, 08:42 AM	#7
Grenouille Junior Member Location: France Join Date: Dec 2012 Posts: 4	Hi Paolo, thanks for replying. Did you include Input samples in your analysis? I know that the usefulness of such data is controversial but in my case I think I am really obliged to take Inputs into account; I am working with several different genotypes, while only one reference genome is available. That means that, due to structural variants, for some samples I could have some regions that show no coverage because the sequence is absent in that particular genotype, and here is where the Input becomes essential. Are you working with similar data?

04-25-2017, 03:58 AM	#8
tanaybhatt Junior Member Location: Bangalore Join Date: Apr 2017 Posts: 3	Can you direct me to a published paper where DMR analysis has been done without taking the input in consideration (or without input subtraction/normalization)?