Gene level methylation values
1
0
Entering edit mode
Arindam ▴ 80
@ag1805x-15215
Last seen 4 weeks ago
University of Eastern Finland

How to extract gene level methylation values?

I used Bismark to align BS-Seq data followed by methylKit for differential methylation identification. By default, the analysis considers CpG sites, but by using suitable annotation BED file, we can perform region level (promoters, exons, TSS) analysis. I was wondering how to obtain gene level methylation values?

methylation methylKit bsseq • 1.6k views
ADD COMMENT
0
Entering edit mode
alex.gos90 ▴ 10
@alexgos90-13597
Last seen 3 months ago
Germany

Hi @ag1805x,

dependent on what you are interested in there may be multiple ways to get gene level values. You may perform a gene-level (aka regional) analysis by using the regionCounts() function to extract a methylRaw/methylRawList object for each region of interest as described here in the vignette and proceed as usual with the downstream analysis. Or if you want to keep the single CpGs statistics you could first run the analysis and then extract CpG sites overlapping your genes of interest using the selectByOverlap() function as described in the vignette. Finally, you could convert your analyzed methylDiff object into a GRanges object and proceed using GenomicRanges functions.

Best, Alex

ADD COMMENT

Login before adding your answer.

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