Question: CpG Density with respect to TSS
gravatar for ResearchR
2.2 years ago by
ResearchR70 wrote:

Dear all,

I would like to replicate an analysis from the paper "Decoding the regulatory landscape of medulloblastoma using DNA methylation sequencing". In more detail, I try to reproduce the plot 1a. Does anyone of you know, how this plot has been created. I am struggling to calculate the density of CpGs with respect to the TSS location. It seems to that this is actually a density plot, with two dimensions.

Help is very well appreciated!

Thanks in advance,


ADD COMMENTlink modified 2.2 years ago • written 2.2 years ago by ResearchR70

Reads to me link like a simple percentage of CpGs per total nucleotides in a certain bin. What I would probably do is to take the center position of all your TSS in bed format. From this go 20kb in each direction, and then subdivide this 40kb windows into bins of, say, 100bp. This all can be done with bedtools makewindows. Once you have this, simply count the number of CpGs per bin, and then divide this number by 100 to get the density. This is just a technical solution. Which regions you then take depends on the question that you want to answer.

ADD REPLYlink modified 2.2 years ago • written 2.2 years ago by ATpoint38k

Thanks for the suggestion! I will give it a try!

ADD REPLYlink written 2.2 years ago by ResearchR70
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1542 users visited in the last hour