4.3.6. GO and KEGG Enrichment Analysis of DMR-Related Genes

GO enrichment analysis of DMR-related genes was performed using the GOseq R package [98], which also corrects for gene length bias. A GO term was considered to be significantly enriched in DMR-related genes at a corrected *p*-value threshold of 0.05. KEGG analysis (http://www.genome.jp/ kegg/), an approach for understanding high-level functions and relationships in biological systems, was applied to uncover the pathway enrichment of DMR-associated genes [99]. KOBAS software [100] was used to test for statistical enrichment of DMR-related genes, which were then subdivided into all, hypermethylation, and hypomethylation categories and assigned to KEGG pathways.
