Objective: To identify genes with aberrant promoter methylation for developing novel diagnostic markers and therapeutic targets against primary colorectal cancer (CRC).
Methods: Two paired CRC and adjacent normal tissues were collected from two CRC patients. A Resi: MBD2b protein-sepharose-4B column was used to enrich the methylated DNA fragments. Difference in the average methylation level of each DNA methylation region between the tumor and control samples was determined by log2 fold change (FC) in each patient to screen the differentially methylated DNA regions. Genes with log2FC value ≥4 or ≤-4 were identified to be hypermethylated and hypomethylated, respectively. Then, the underlying functions of methylated genes were speculated by Gene Ontology database and pathway enrichment analyses. Furthermore, a protein-protein interaction network was built using Search Tool for the Retrieval of Interacting Genes/Proteins database, and the transcription factor binding sites were screened via the Encyclopedia of DNA Elements (ENCODE) database.
Results: Totally, 2,284 and 1,142 genes were predicted to have aberrant promoter hypermethylation or hypomethylation, respectively. MAP3K5, MAP3K8, MAPK14, and MAPK9 with promoter hypermethylation functioned via MAPK signaling pathway, focal adhesion, or Wnt signaling pathway, whereas MAP2K1, MAPK3, MAPK11, and MAPK7 with promoter hypomethylation functioned via TGF-beta signaling pathway, neurotrophin signaling pathway, and chemokine signaling pathway. CREBBP, PIK3R1, MAPK14, APP, ESR1, MAPK3, and HRAS were the seven hubs in the constructed protein-protein interaction network. RPL22, RPL36, RPLP2, RPS7, and RPS9 were commonly regulated by transcription factors, and YY1 and IRF4 were hypermethylated.
Conclusion: MAPK14, MAPK3, HRAS, YY1, and IRF4 may be considered as potential biomarkers for early diagnosis and therapy of CRC.
Keywords: aberrant DNA methylation; microarray analysis; pathway enrichment analysis; primary colorectal cancer; transcription factor.