Question: CpG site not annotated with annotatr
gravatar for soli5
13 months ago by
soli50 wrote:


I am using annotatr to annotate differentially methylated cytosines and I have a list of CpGs not included in the exons, promoters, introns or intergenic results. Where could they be classified if it's not in one of those regions. I checked one of the non annotated CpGs and I found it 4Kb downstream a gene. Is there an explanation for CpGs not being present in those regions?



cpg annotatr • 195 views
ADD COMMENTlink modified 13 months ago by rcavalca130 • written 13 months ago by soli50
Answer: CpG site not annotated with annotatr
gravatar for rcavalca
13 months ago by
United States
rcavalca130 wrote:


Thanks for using annotatr. If it were 4kb upstream of a gene, I would answer that it falls in the 1to5kb regions, which is distinct from the promoters. However, you said 4kb downstream which suggests that it should probably be intergenic. To help me figure this out, could you give me the location of the CpG and what genome it's with respect to, and I'll have a look on the Genome Browser to see what could be going on.



ADD COMMENTlink written 13 months ago by rcavalca130

Hi Raymond,

Thanks for your reply and help. I actually meant upstream. In any case, in the list I got there are different situations. For example chr12:9818891 falls in a region that to me could be promoter, 1to5Kb or intergenic,  but it was not annotated as promoter or intergenic, would it be in a 1to5Kb region?. I did not perform that annotation.
I give you a small list of examples (from this particular comparison I got 1195 CpGs not in those 4 categories), it would really help understand where these CpGs would go.

Another question, I see overlaps between exon and promoter lists for example, I understand that this is due to isoform configuration (one cytosine could be in the exon of an isoform and in the promoter of another isoform). But other than exon, promoter and intron annotations, are the other categories supposed to be mutually exclusive?






I hope my questions are clear. Thank you


ADD REPLYlink written 13 months ago by soli50

Hi Sol,

Since it is upstream, my guess would be that it would show up in the 1to5kb annotation. The promoter annotation only extends 1kb upstream from the TSS. I verified in the UCSC Genome Browser for hg19.

As for the list of other CpGs, I looked at a handful at random and they seem to primarily be in that 1to5kb region upstream of a TSS. I would recommend adding that annotation and seeing if that does the trick.



ADD REPLYlink modified 13 months ago • written 13 months ago by rcavalca130
Please log in to add an answer.


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