Hi there,
I am using CHAMP to perform differential methylation analysis, including DMP and DMR. Regarding DMR, I used DMRcate algorithm and I am a little bit confused about the result.
In the derived table of DMRcate embedded in CHAMP, I got a column of overlapping.promoters. For example:
seqnames start end width strand no.cpgs minfdr Stouffer maxbetafc meanbetafc
chr20:61049813-61051915 chr20 61049813 61051915 2103 * 27 0.000000e+00 0.000000e+00 0.4770680 0.35455081
chr20:57424521-57431303 chr20 57424521 57431303 6783 * 77 0.000000e+00 2.725637e-268 -0.2084268 0.07728631
chr20:24448859-24452131 chr20 24448859 24452131 3273 * 21 0.000000e+00 1.434465e-236 0.4263522 0.28291317
chr20:21491781-21498921 chr20 21491781 21498921 7141 * 26 3.342271e-208 4.594347e-231 0.4385002 0.25698803
chr20:61806628-61810795 chr20 61806628 61810795 4168 * 23 0.000000e+00 8.953525e-215 0.4182034 0.24114120
chr20:37351945-37358077 chr20 37351945 37358077 6133 * 23 1.003282e-251 2.186334e-208 0.3632555 0.20988842
overlapping.promoters
chr20:61049813-61051915 GATA5-001, RP13-379O24.3-001
chr20:57424521-57431303 GNAS-050, GNAS-037, GNAS-058, GNAS-001, GNAS-009, GNAS-049, GNAS-AS1-001, GNAS-036, GNAS-059, GNAS-012, GNAS-AS1-004
chr20:24448859-24452131 SYNDIG1-001
chr20:21491781-21498921 NKX2-2-001, NKX2-2-AS1-001
chr20:61806628-61810795 MIR124-3-201, RP5-963E22.4-001
chr20:37351945-37358077 SLC32A1-001
And what confusing me is why there is a number coming after each gene, like GATA5-001, what is that supposed to mean? And how can I select the significant DMR, just refer to mindfr or stouffer? What about the threshold? I didn't see the threshold in its R help document.
Hope someone could give me a hint! Many thanks advanced.
Did you find an answer to this? I had the same question.
Thanks.