Identification of aberrantly methylated differentially expressed genes and pro-tumorigenic role of KIF2C in melanoma

Background: Skin Cutaneous Melanoma (SKCM) is known as an aggressive malignant cancer, which could be directly derived from melanocytic nevi. However, the molecular mechanisms underlying the malignant transformation of melanocytes and melanoma tumor progression still remain unclear. Increasing research showed significant roles of epigenetic modifications, especially DNA methylation, in melanoma. This study focused on the identification and analysis of methylation-regulated differentially expressed genes (MeDEGs) between melanocytic nevus and malignant melanoma in genome-wide profiles. Methods: The gene expression profiling datasets (GSE3189 and GSE114445) and gene methylation profiling datasets (GSE86355 and GSE120878) were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) and differentially methylated genes (DMGs) were identified via GEO2R. MeDEGs were obtained by integrating the DEGs and DMGs. Then, a functional enrichment analysis of MeDEGs was performed. STRING and Cytoscape were used to describe the protein-protein interaction (PPI) network. Furthermore, survival analysis was implemented to select the prognostic hub genes. Next, we conducted gene set enrichment analysis (GSEA) of hub genes. To validate, SKCM cell culture and lentivirus infection was performed to reveal the expression and behavior pattern of KIF2C. Patients and specimens were collected and then immunohistochemistry (IHC) staining was conducted. Results: We identified 237 hypomethylated, upregulated genes and 182 hypermethylated, downregulated genes. Hypomethylation-upregulated genes were enriched in biological processes of the oxidation-reduction process, cell proliferation, cell division, phosphorylation, extracellular matrix disassembly and protein sumoylation. Pathway enrichment showed selenocompound metabolism, small cell lung cancer and lysosome. Hypermethylation-downregulated genes were enriched in biological processes of positive regulation of transcription from RNA polymerase II promoter, cell adhesion, cell proliferation, positive regulation of transcription, DNA-templated and angiogenesis. The most significantly enriched pathways involved the transcriptional misregulation in cancer, circadian rhythm, tight junction, protein digestion and absorption and Hippo signaling pathway. After PPI establishment and survival analysis, seven prognostic hub genes were CKS2, DTL, KIF2C, KPNA2, MYBL2, TPX2, and FBL. Moreover, the most involved hallmarks obtained by GSEA were E2F targets, G2M checkpoint and mitotic spindle. Importantly, among the 7 hub genes, we found that down-regulated level of KIF2C expression significantly inhibited the proliferative ability of SKCM cells and suppressed the metastasis capacity of SKCM cells. Conclusions: Our study identified potential aberrantly methylated-differentially expressed genes participating in the process of malignant transformation from nevus to melanoma tissues based on comprehensive genomic profiles. Transcription profiles of CKS2, DTL, KIF2C, KPNA2, MYBL2, TPX2, and FBL provided clues of aberrantly methylation-based biomarkers, which might improve the development of precision medicine. KIF2C plays a pro-tumorigenic role and potentially inhibited the proliferative ability in SKCM.