This is a beta test.
Question: How to map CpG ids to nearest genes?
0
Entering edit mode

I am working in dna methylation data and I am having a .csv file containing cpg ids. I want to map the cpg sites to nearest gene. How to do that in R?

Thanks in advance

1
Entering edit mode

Load the relevant txdb database and use the nearest() command, presumably with the transcripts.

ADD COMMENTlink 2.2 years ago Devon Ryan 90k
Entering edit mode
0

Hi, Devon. Could you give a detail? I also have a question. When two genes are near the CpG, how to decide the nearest one. Should the TSS, promoter and gene body be considered? @Devon Ryan

ADD REPLYlink 2.0 years ago
Kangli Wang
• 0
Entering edit mode
0

Honestly it's a bit arbitrary what one considers. Biologically speaking, it'd make sense to consider genes downstream (when considering the gene's strand) or containing a SNP, but I'd have to check and see if that's possible with nearest().

ADD REPLYlink 2.0 years ago
Devon Ryan
90k
0
Entering edit mode

If you are working with Illumina 450k or similar, you can use the manifest table to annotate your probes. See here: https://bioconductor.org/packages/release/data/annotation/html/IlluminaHumanMethylation450kmanifest.html

ADD COMMENTlink 2.2 years ago nandobonf • 0
0
Entering edit mode

I try to install IlluminaHumanMethylation450k.db package for annotation but I came across this error.

Error: package or namespace load failed for 'IlluminaHumanMethylation450k.db':
 .onLoad failed in loadNamespace() for 'IlluminaHumanMethylation450k.db', details:
  call: NULL
  error: 'fun' is defunct.
Use 'FDb.InfiniumMethylation.hg19' instead.
Use 'FDb.InfiniumMethylation.hg18' instead.
Use 'mapToGenome() function in minfi or methylumi' instead.
See help("Defunct")
Error: loading failed
Execution halted
ERROR: loading failed
* removing 'C:/R/R-3.4.1/library/IlluminaHumanMethylation450k.db'
ADD COMMENTlink 2.2 years ago preethibagopi5 • 0 • updated 2.2 years ago Devon Ryan 90k
Entering edit mode
0

Please note the "ADD COMMENT" button, which can be used to reply to specific posts.

ADD REPLYlink 2.2 years ago
Devon Ryan
90k

Login before adding your answer.

Powered by the version 1.6