This is my first time dealing with methylation data. I would like to know if my approach is appropriate or biased?
I have a list of differential variable cpgs, and I have to perform Gene ontology analysis. I read multiple published papers and majority of them have used "gometh" function "missMethyl" package. Instead of using "gometh" function, I obtained the Entrez Gene IDs that the significant CpGs are mapped to using "getMappedEntrezIDs" function. Further, performed GO Enrichment Analysis on this gene set using "enrichgo" function of clusterProfiler.
Is this approach correct or is it baised?
1 answer
I would prefer the gometh approach because it has been specificially designed to adress the biases related to methylation data. Indeed as it is very well explained here (https://genomebiology.biomedcentral.com/articles/10.1186/s13059-021-02388-x), the assumption of independently measured genes is not respected because some CpGs are assigned to several genes. Additionally, some genes are more likely to occur in your DMP list because they inherently contain more CpGs measured by the array. For all these reasons, you should prefer gometh over other methods.
Log in to answer this question.
-
missMethyl probe number bias plot
written by ana •
-
getMappedEntrezIDs from missmethyl package R
written by Rut_MDLV •
-
Illumina methylation EPIC V2 array
written by juliviglino •
-
Extract geneID following GO enrichment analysis
written by hellokwmin •
-
Making plots for enrichment and KEGG results
written by ymervenuryavuz •
-
Error using gometh from missMethyl
written by olivia.grant •
-
Adding geneset to missMethyl
written by solveigbrunstad •
-
How to calculate methylation FC for genes based on CpG FC data?
written by maria2019
-
genlist for clusterprofiler
written by Omics data mining
-
R bioconductor in analysying methylation data
written by preethibagopi5 •
