Gene level methylation values
Entering edit mode
ag1805x ▴ 50
Last seen 4 months ago
University of Allahabad

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 • 198 views
Entering edit mode
alex.gos90 ▴ 10
Last seen 3 months ago

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


Login before adding your answer.

Traffic: 267 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6