<p>EPIC data can be processed in the same manner as the previous iteration of methylation array data from Illumina (450k). This means that starting with .idat files, normalization should be performed (for example, via the minfi package). A <a href="http://biorxiv.org/content/early/2016/07/23/065490" rel="noreferrer">recent paper</a> from the creators of minfi is particularly helpful because it makes clear that normalized EPIC data from their package can be immediately compared against, for example, <a href="https://wiki.nci.nih.gov/display/TCGA/Data+level" rel="noreferrer">level 3 TCGA data</a>.</p>

<p>After that, I suggest using the manifest to attach genomic coordinates to your probes and segregating them into functional regions. By testing differential methylation in only regulatory regions, for example, you can increase the statistical power by reducing the overall number of tests to the ones you expect to yield major differences.</p>

<p>There are existing packages out there for differential methylation analysis, but without knowing your replicate structure or aims, it is difficult to point you in the right direction.</p>

CC BY-SA 3.0
2017-06-06T23:10:14.527Z
2017-06-06T23:10:14.527Z
Ben D.
<p>I recently got some <a href="https://www.illumina.com/products/by-type/microarray-kits/infinium-methylation-epic.html" rel="noreferrer">EPIC DNA methylation</a> data and I was wondering what are some good practices to follow? 
I am interested in knowing about normalization and differential analysis. Thank you.</p>

CC BY-SA 3.0
2017-06-06T21:49:12.547Z
Hi deepseas, welcome and thank you for posting a question to Bioinformatics StackExchange. I've added a link to the Infinium EPIC web page for people who might be curious about what EPIC DNA methylation data is. Your question is quite broad, and would be substantially improved by adding more of a story behind why you want to do this analysis (e.g. is this a single individual study? do you want to compare with RNASeq data?). Bioinformatics is a very large subject area, and there are many different ways to do [almost] the same thing, so any help with clarifying questions is greatly appreciated.
gringer
CC BY-SA 3.0
2017-06-06T23:27:08.023Z
Thanks for the response! This data has 16 samples from different individuals with 5 groups (5+6+1+2+2). It will be difficult to get power for groups with 1 and 2. but the main aim is to get an idea about certain known genes of interest.
deepseas
CC BY-SA 3.0
2017-06-08T13:10:51.933Z
It will be difficult to make contrasts like edgeR,limma or deseq2 makes with low number of replicates per condition. 2 might still work but 1 in a single condition will not hold good. Ideally when you work with EPIC arrays or 450k both will work with any standard 450k analysis tool. Try tools like DimMER/Rnbeads. However things hold better with higher replicates per groups. All the best!
ivivek_ngs
<r><normalization><methylation><microarray>
What are some good practices to follow during EPIC DNA methylation data analysis?
# What are some good practices to follow during EPIC DNA methylation data analysis?

(Asked by: username_0 on 2017-06-06)

I recently got some [EPIC DNA methylation](https://www.illumina.com/products/by-type/microarray-kits/infinium-methylation-epic.html) data and I was wondering what are some good practices to follow?
I am interested in knowing about normalization and differential analysis. Thank you.

**Tags:** r, normalization, methylation, microarray
> **Comments:**
> **username_1**: Hi deepseas, welcome and thank you for posting a question to Bioinformatics StackExchange. I've added a link to the Infinium EPIC web page for people who might be curious about what EPIC DNA methylation data is. Your question is quite broad, and would be substantially improved by adding more of a story behind why you want to do this analysis (e.g. is this a single individual study? do you want to compare with RNASeq data?). Bioinformatics is a very large subject area, and there are many different ways to do [almost] the same thing, so any help with clarifying questions is greatly appreciated.  
> **username_0**: Thanks for the response! This data has 16 samples from different individuals with 5 groups (5+6+1+2+2). It will be difficult to get power for groups with 1 and 2. but the main aim is to get an idea about certain known genes of interest.  
> **username_2**: It will be difficult to make contrasts like edgeR,limma or deseq2 makes with low number of replicates per condition. 2 might still work but 1 in a single condition will not hold good. Ideally when you work with EPIC arrays or 450k both will work with any standard 450k analysis tool. Try tools like DimMER/Rnbeads. However things hold better with higher replicates per groups. All the best!  

---

## 1 answers

### Answer by username_3 on 2017-06-06. Score: 6 | Accepted answer

EPIC data can be processed in the same manner as the previous iteration of methylation array data from Illumina (450k). This means that starting with .idat files, normalization should be performed (for example, via the minfi package). A [recent paper](http://biorxiv.org/content/early/2016/07/23/065490) from the creators of minfi is particularly helpful because it makes clear that normalized EPIC data from their package can be immediately compared against, for example, [level 3 TCGA data](https://wiki.nci.nih.gov/display/TCGA/Data+level).

After that, I suggest using the manifest to attach genomic coordinates to your probes and segregating them into functional regions. By testing differential methylation in only regulatory regions, for example, you can increase the statistical power by reducing the overall number of tests to the ones you expect to yield major differences.

There are existing packages out there for differential methylation analysis, but without knowing your replicate structure or aims, it is difficult to point you in the right direction.


---

