Identification of key candidate genes involved in melanoma metastasis

Metastasis is the most lethal stage of cancer progression. The present study aimed to investigate the underlying molecular mechanisms of melanoma metastasis using bioinformatics. Using the microarray dataset GSE8401 from the Gene Expression Omnibus database, which included 52 biopsy specimens from patients with melanoma metastasis and 31 biopsy specimens from patients with primary melanoma, differentially expressed genes (DEGs) were identified, subsequent to data preprocessing with the affy package, followed by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. A protein-protein interaction (PPI) network was constructed. Mutated genes were analyzed with 80 mutated cases with melanoma from The Cancer Genome Atlas. The overall survival of key candidate DEGs, which were within a filtering of degree >30 criteria in the PPI network and involved three or more KEGG signaling pathways, and genes with a high mutation frequency were delineated. The expression analysis of key candidate DEGs, mutant genes and their associated genes were performed on UALCAN. Of the 1,187 DEGs obtained, 505 were upregulated and 682 were downregulated. ‘Extracellular exosome’ processes, the ‘amoebiasis’ pathway, the ‘ECM-receptor interaction’ pathway and the ‘focal adhesion’ signaling pathway were significantly enriched and identified as important processes or signaling pathways. The overall survival analysis of phosphoinositide-3-kinase regulator subunit 3 (PIK3R3), centromere protein M (CENPM), aurora kinase A (AURKA), laminin subunit α 1 (LAMA1), proliferating cell nuclear antigen (PCNA), adenylate cyclase 1 (ADCY1), BUB1 mitotic checkpoint serine/threonine kinase (BUB1), NDC80 kinetochore complex component (NDC80) and protein kinase C α (PRKCA) in DEGs was statistically significant. Mutation gene analysis identified that BRCA1-associated protein 1 (BAP1) had a higher mutation frequency and survival analysis, and its associated genes in the BAP1-associated PPI network, including ASXL transcriptional regulator 1 (ASXL1), proteasome 26S subunit, non-ATPase 3 (PSMD3), proteasome 26S subunit, non ATPase 11 (PSMD11) and ubiquitin C (UBC), were statistically significantly associated with the overall survival of patients with melanoma. The expression levels of PRKCA, BUB1, BAP1 and ASXL1 were significantly different between primary melanoma and metastatic melanoma. Based on the present study, ‘extracellular exosome’ processes, ‘amoebiasis’ pathways, ‘ECM-receptor interaction’ pathways and ‘focal adhesion’ signaling pathways may be important in the formation of metastases from melanoma. The involved genes, including PIK3R3, CENPM, AURKA, LAMA1, PCNA, ADCY1, BUB1, NDC80 and PRKCA, and mutation associated genes, including BAP1, ASXL1, PSMD3, PSMD11 and UBC, may serve important roles in metastases of melanoma.