Integration of protein interaction and gene co-expression information for identification of melanoma candidate genes

Cutaneous melanoma is an aggressive form of skin cancer that causes death worldwide. Although much has been learned about the molecular basis of melanoma genesis and progression, there is also increasing appreciation for the continuing discovery of melanoma genes to improve the genetic understanding of this malignancy. In the present study, melanoma candidate genes were identified by analysis of the common network from cancer type-specific RNA-Seq co-expression data and protein-protein interaction profiles. Then, an integrated network containing the known melanoma-related genes represented as seed genes and the putative genes represented as linker genes was generated using the subnetwork extraction algorithm. According to the network topology property of the putative genes, we selected seven key genes (CREB1, XPO1, SP3, TNFRSF1B, CD40LG, UBR1, and ZNF484) as candidate genes of melanoma. Subsequent analysis showed that six of these genes are melanoma-associated genes and one (ZNF484) is a cancer-associated gene on the basis of the existing literature. A signature comprising these seven key genes was developed and an overall survival analysis of 461 cutaneous melanoma cases was carried out. This seven-gene signature can accurately determine the risk profile for cutaneous melanoma tumors (log-rank P=3.27E-05) and be validated on an independent clinical cohort (log-rank P=0.028). The presented seven genes might serve as candidates for studying the molecular mechanisms and help improve the prognostic risk assessment, which have clinical implications for melanoma patients.