Biostar Beta. Not for public use.
How To Annotate Methylation 450K?
0
Entering edit mode
2.3 years ago
liu4gre • 200
United States

I have the file from encode, with the calculated beta values and probe id. I am wondering how to map these probes to genes symbol? And how to determine their locations, such as promoter, 3'UTR or gene body?

Thanks

methylation • 9.3k views
ADD COMMENTlink
0
Entering edit mode

I am wonder do you have any success over your question !! i am at same situation .

ADD REPLYlink
1
Entering edit mode
13 months ago
Freiburg, Germany

You could start with the IlluminaHumanMethylation450k Bioconductor package, which contains mappings for at least part of what you're looking for (maybe all of it, I haven't checked).

ADD COMMENTlink
0
Entering edit mode

I think this is probably sufficient. I know IMA provides region-level calculations for the categories that you listed, but that isn't a Bioconductor package and might be using something else.

If not, I know for certain they are included in the .bpm file that you can download here:

https://support.illumina.com/array/array_kits/infinium_humanmethylation450_beadchip_kit/downloads.ilmn

The website was being slow for me, but it should work. If not, you can just use the copy of the .bpm file included with the COHCAP demo dataset (the standalone version of COHCAP also adds these labels at the CpG site level, but this is probably not the most efficient way to accomplish this single task):

http://sourceforge.net/projects/cohcap/

Once you delete some of early and final lines (or write code to ignore those lines), this is just a tab-delimited text file that you can parse.

ADD REPLYlink
0
Entering edit mode

Hi Charles, Can you give me some introduce to .bmp files (methylation 450K, such as GPL13534_HumanMethylation450_15017482_v.1.2.bpm.gz)? It is binary files, I didn't find too much introduction on the interenet.

ADD REPLYlink
0
Entering edit mode

If you don't want to parse the Illumina annotation files, you can use these commands for the Bioconductor code:

library("IlluminaHumanMethylation450kanno.ilmn12.hg19")
data("IlluminaHumanMethylation450kanno.ilmn12.hg19")
annotation.table = getAnnotation(IlluminaHumanMethylation450kanno.ilmn12.hg19)
ADD REPLYlink
0
Entering edit mode
17 months ago
Researcher • 30

Were anyone able to get this 450K annotation file? Please do share with me too.

ADD COMMENTlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.1