Annotation for EPICv2
1
0
Entering edit mode
nurkenber • 0
@bba1982e
Last seen 5 hours ago
Kazakhstan

I am trying to get a more complete annotation for EPICv2. The annotation from Illumina "Infinium MethylationEPIC v2.0" (link below), has 614365 out of 937055 probes not having information on associated gene and a lot of "unknowns" when it comes to gene features such as TSS1500 or 5UTR, as seen on picture below. When I check the position of CpG on UCSC I can find related gene. For example, cg06979118 is overlapped with SHANK2 gene, but has no gene annotation in "Infinium MethylationEPIC v2.0".

enter image description here

I followed ChAMP tutorial for analysis, that uses ChAMPdata (v2.31.1) package, which loads up the annotation based on "Infinium MethylationEPIC v2.0". The newer and developmental ChAMPdata versions do not have annotation for EPICv2 at all.

I have few questions

1) Has anyone used ChAMP and was able to load more complete annotation for EPICv2, meaning to find ChAMPdata package with complete annotation for EPICv2, since functions in ChAMP call data("probe.features.epicv2") even in champ.GSEA,

2) I have a half solution to my problem, I can annotate DMP after I find them but even then I am still not sure where to get a complete annotation. I can overlap probes used for EPICv1 with probes used for EPICv2. But then I am not sure how to perform GSEA analysis in this case since champ.GSEA calls data("probe.features.epicv2"), even though I provide it with DMPs and DMRs. Maybe someone could suggest the pipeline for that?

Infinium MethylationEPIC v2.0: https://support.illumina.com/array/array_kits/infinium-methylationepic-beadchip-kit/downloads.html

Tutorial for ChAMP: https://github.com/YuanTian1991/ChAMP-DemoRun/blob/main/EPICv2/illumina_demo_data_iScan/main.md

EpicV2 ChAMPdata • 59 views
ADD COMMENT
0
Entering edit mode

Don't know if any of these are of use or if these are what ChAMP uses under the hood:

IlluminaHumanMethylationEPICv2anno.20a1.hg38

IlluminaHumanMethylationEPICv2manifest

EPICv2manifest

ADD REPLY
0
Entering edit mode
@james-w-macdonald-5106
Last seen 18 minutes ago
United States

The CpG you point out is in the intronic region of SHANK2, so will not be in any of the locations you have in your barplot.

enter image description here

The annotations we have for this array come directly from Illumina, and we (and here by 'we' I mean Zuguang Gu who submitted the annotations) simply process their data to put into packages that are useful for minfi and ChAMP to use. If you want more sophisticated location information you could always use ChIPseeker or ChIPPeakAnno or the functions in GenomicRanges that those packages are based on.

Login before adding your answer.

Traffic: 847 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