I am currently analyzing the TCGA HNSC methylation data. I have split the probes to each gene, but I find out that for one gene, there will be several probes near that gene. I am wondering how could I get the methylation level to represent for that gene if there are mutli_probes near that gene.
I have read TCGA tutorial or other papers, they said that for genes with multi methylation probes, the probe with most anti-correlated is selected to represent the methylation level of a gene. But I am confused that, for methylation, we mostly care about the methylation in the promoter area. Shouldn't I only focus on the probes near the promoter area? So I am thinking about maybe I should get an average value for the probes near the promoter area to represent the methylation level of that gene. Is that Right?
•
written
by
gongyutang •
1 answer
It does not make biological sense to obtain the average across the entire gene, unless, perhaps, it's a single exon gene like a non-coding RNA or microRNA. You mention also obtaining the average around the promoter region. This may have greater validity; however, I still think that methylation needs to be treated at the probe-set level. For example, at a bi-directional promoter, one may have a probe on both the 5' and 3'; however, high methylation at just one of these can have consequences to downstream gene transcription. Were you to obtain the average, you may miss things like that if one probe had high methylation and the other not. Even at non-bi-directional promoters, this would be an issue.
Kevin
Log in to answer this question.
-
Human Methylation Infinum Bead - probe to gene mapping
written by shrutidesai0208 •
-
TCGA input files for COHCAP
written by dilrajkaur766 •
-
Correlating DNA methylation beta values with Gene expression levels (RNA-seq)
written by sswang25 •
-
DNA methylation probe selection for GDAC Firehose
written by octaviamdl •
-
Defining methylation level of CpG island
written by CY
-
For GEO data, how to calculate the average methylation level of the gene promoter region using R.
written by a511512345
-
Methylation at probe level and gene expression correlation
written by Will •
-
Correlation between gene expression and methylation
written by tujuchuanli
-
Interpreting level-3 TCGA methylseq data
written by mforde84
-
Methylation and RNA seq data
written by ramyar.molania •
你好，我也遇到你同样的问题，从xena下载的TCGA 450k甲基化数据会有一个基因对应多个探针id的情况，请问你后来是怎么解决的呢？
Would you please use english language? This is an international (=english speaking) community ;-)
lol, meet the same problem. I download the methylation data(450k) frome Xena and map the probe id to gene_symbol using Xena illuminaMethyl27_hg38_GDC file. Single gene map to several probe id(such as MGMT) and I don't know which probe should be selected for further analysis.
That is expected due to the fact that methylation at different parts of each gene can have different effects. You should treat each independently. It does not make sense to summarise methylation across each gene.