Hi,

I am dealing with Whole Genome Bisulfite Data (WGBS), form a time series experiment (4 stages), and I am wondering I f I could do Differentially Methylated Regions (DMRs) in a time series continous procedure. I mean, I used to find DMRs between One condition and Another, but not in a time course (DMRs that change sequentially between stage 1 and 4).

My previous thinking was to find those DMRs between stage 1 and 2, 2 and 3 and 3 and 4, and then overlap them using bedtools to find those one that are changing sequentially

DO you know a more straightforward strategy?

Thanks in advance