Annotating Illumna epic probes
1
2
Entering edit mode
11 days ago

I think I have a very easy question but unfortunately I haven't found a good answer yet, it's basically that I would like to annotate, using R, illumina epic probe sites according to their location, We did an experiment using Illumina Epic probes, and I now have my matrix of average beta values differences between two conditions for each probe, and I am trying to create a graph such as (E) and (F) here, where the Y axis is the average beta value difference for the probe and the X-axis corresponds to the location of this probe (transcription factor binding site, open chromatin, DNase I hypersensitive sites, etc). For this I need to annotate the probes, how can I do that?

I've been looking at the IlluminaMethylationManifest but it is not clear at all how to use it for this.

Any suggestions appreciated.

R illumina EPIC • 378 views
ADD COMMENT
0
Entering edit mode

Just in case, check if you're using EPIC v1 and EPIC v2, because the annotations are different (see this and this packages)

ADD REPLY
0
Entering edit mode

Thanks, I am using EPIC v2, yet, it is not clear to me how to get right information from the manifest.

ADD REPLY
0
Entering edit mode

If you install and load these packages, you can get a data.frame within R with the information of probe names, coordinates, genes mapped to, etc., as such:

library(IlluminaHumanMethylationEPICv2manifest)
library(IlluminaHumanMethylationEPICv2anno.20a1.hg38)

annotationEPIC.v2 <- getAnnotation(IlluminaHumanMethylationEPICv2anno.20a1.hg38)
annotationEPIC.v2 <- as.data.frame(annotationEPIC.v2)
ADD REPLY
0
Entering edit mode
11 days ago

Using minfi, this should be useful:

library(minfi)
library(minfiDataEPIC)
library(GenomicRanges)

epic.anno <- getAnnotation(minfiDataEPIC::MsetEPIC)
epic.anno <- epic.anno[order(epic.anno$chr, epic.anno$pos), ]

cpgs.ranges <- GRanges(epic.anno$chr, IRanges(epic.anno$pos, epic.anno$pos, names = epic.anno$Name))
cpgs.ranges@elementMetadata <- epic.anno
ADD COMMENT
0
Entering edit mode

This is great! although it seems that it is for Illumina EPIC V1, and I am using V2

ADD REPLY

Login before adding your answer.

Traffic: 2379 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6