[This article was first published on Chitka, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here)
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.
Nowadays, regions with relatively lower CpG density are gaining importance in DNA methylation studies. This is based on the fact that a majority of CpG rich regions (CpG islands) are non-dynamic and less variant in terms of methylation status probed across a variety of tissues and cell populations (Irizarry 2009, Ziller 2013). It is now proven that methylation is more dynamic along the CpG shores (< 2kb flanking CpG Islands) and CpG shelves (<2kb flanking outwards from a CpG shore). While it is easy to retreive the genomic co-ordinates of CpG Islands from UCSC browser, public resources for retrieving the genomic co-ordinates for shores and shelfs were missing or scarce. Here, I am displaying the code (in R using Bioconductor package GenomicRanges) I use for generating objects for CpG islands, CpG island shores and CpG island shelves. The resulting objects are GRanges objects which could be used in multiple downstream applications in Bioconductor related packages.
To leave a comment for the author, please follow the link and comment on their blog: Chitka.
R-bloggers.com offers daily e-mail updates about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job.
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.