Identifying a Risk Signature of Methylation-Driven Genes as a Predictor of Survival Outcome for Colon Cancer Patients

Aberrant expression of gene is driven by its promoter methylation and is the key molecular basis of carcinogenic processes. This study aimed at identifying a risk signature of methylation-driven (MD) genes and evaluating its prognostic value for colon cancer (CC) patients. The expression profiles of methylation and mRNA in CC samples were obtained from the TCGA database, and the MethylMix algorithm was used to identify MD genes. The relationships between their expression levels and overall survival (OS) of CC patients were analyzed, and a prognostic signature of MD genes was established. The risk score of gene signature was calculated, and the median was used to divide all patients into high (H) and low (L) risk groups. The prognostic value of gene signature was tested by the TCGA cohort and an independent validation cohort (GSE17538 dataset). In total, 69 MD genes were identified, and 7 were associated with OS of CC patients. Ultimately, 4 (TWIST1, LDOC1, EPHX3, and STC2) were screened out to establish a risk signature. The H-risk patients (>0.923) had a worse OS than L-risk patients (≤0.923) in both the TCGA (5-year cumulative survival: 52.9% vs 72.0%, P=0.005) and GSE17538 cohort (49.4% vs 69.3%, P=0.004). The AUC values of MD genes signature for the prediction of 3- and 5-year OS were 0.648 and 0.643 in the TCGA dataset and 0.634 and 0.624 in the GSE17538 dataset, respectively. The risk signature of four MD genes was identified as an independent predictor of OS for CC patients (HR for TCGA dataset: 2.071, 95% CI=1.196-3.586, P=0.009; HR for GSE17538 dataset: 2.021, 95% CI=1.290-3.166, P=0.002). The risk signature of four MD genes might be a useful prognostic tool and help doctors improve the clinical management of CC patients.