This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Analysing MIRA-seq methylation data

Hi,

I have a MIRA-seq data that I want to analyze. In the past I have used Bisulfite data for methylation analysis. There are a number of options available for Bisulfite data (e.g. aligners e.t.c. Aligning with Bismark etc and downstream analysis with methylkit/methpipe e.t.c.) I want to know what are the options or pipeline available for MIRA-seq data. Also what are the best practices for DNA methylation analysis with MIRA-seq.

methylation-analysis mira-seq
• 3,799 views • link
updated by Ram Ram • written by Saad Khan Saad Khan

1 answer

MIRA is a MeDIP-based strategy. In general, you'll want to use a normal aligner (such as BWA or Bowtie) and a ChIP-Seq/ChIP-chip peak-calling algorithm (MATS, MACS, etc.). I've used MATS on MIRA tiling array data, but I can't absolutely promise that is the best strategy.

• 0 views • link
updated by Ram Ram • written by Charles Warden Charles Warden

Can MEDIPS be used for the same purpose or do we have to use MACS/MATS?

Thanks

• 0 views • link
written by smk • smk

I haven't tried it, but I Googled the MEDIPS Bioconductor package and I think it looks appropriate. There are lots of valid peak calling algorithms out there - I just think MATS/MACS are the most commonly used.

• 1 views • link
updated by Ram Ram • written by Charles Warden Charles Warden

How do you do differential methylation analysis once you do peak calling with MATS/MACS

• 1 views • link
written by smk • smk

The peak calling compares the signal between the two groups (which perhaps has been normalized against an input reference). So, peaks called due to a significant difference in signal between the groups are most likely to be differentially methylated.

• 0 views • link
written by Charles Warden Charles Warden

Log in to answer this question.

More posts like this