Molecular Medicine Reports Molecular Medicine Reports 1791-2997 1791-3004 D.A. Spandidos 10.3892/mmr.2017.8219 mmr-17-02-2907 Articles A multigene support vector machine predictor for metastasis of cutaneous melanoma WeiDong Department of Plastic and Esthetic Surgeries, Sichuan Academy of Medical Sciences and Sichuan Provincial People's Hospital, Chengdu, Sichuan 610072, P.R. China Correspondence to: Dr Dong Wei, Department of Plastic and Esthetic Surgeries, Sichuan Academy of Medical Sciences and Sichuan Provincial People's Hospital, 32 West Second Section First Ring Road, Chengdu, Sichuan 610072, P.R. China, E-mail: dongdhhdhdh@aliyun.com 022018 07122017 17 2 2907 2914 09032017 21112017 Copyright: © Wei et al. 2018 This is an open access article distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non-commercial and no modifications or adaptations are made.

Gene expression profiles of cutaneous melanoma were analyzed to identify critical genes associated with metastasis. Two gene expression datasets were downloaded from Gene Expression Omnibus (GEO) and another dataset was obtained from The Cancer Genome Atlas (TCGA). Differentially expression genes (DEGs) between metastatic and non-metastatic melanoma were identified by meta-analysis. A protein-protein interaction (PPI) network was constructed for the DEGs using information from BioGRID, HPRD and DIP. Betweenness centrality (BC) was calculated for each node in the network and the top feature genes ranked by BC were selected to construct the support vector machine (SVM) classifier using the training set. The SVM classifier was then validated in another independent dataset. Pathway enrichment analysis was performed for the feature genes using Fisher's exact test. A total of 798 DEGs were identified and a PPI network including 337 nodes and 466 edges was then constructed. Top 110 feature genes ranked by BC were included in the SVM classifier. The prediction accuracies for the three datasets were 96.8, 100 and 94.4%, respectively. A total of 11 KEGG pathways and 13 GO biological pathways were significantly over-represented in the 110 feature genes, including endometrial cancer, regulation of actin cytoskeleton, focal adhesion, ubiquitin mediated proteolysis, regulation of apoptosis and regulation of cell proliferation. A SVM classifier of high prediction accuracy was acquired. Several critical genes implicated in melanoms metastasis were also revealed. These results may advance understanding of the molecular mechanisms underlying metastasis, and also provide potential therapeutic targets.

 cutaneous melanoma metastasis differentially expressed genes protein-protein interaction interaction SVM classifier Introduction

Cutaneous melanoma is the most dangerous type of skin cancer. In the United States, the five-year survival rates is 98% among those with localized disease and 17% among those in whom spread has occurred. Apparently, metastatis is the most important factor affecting prognosis of melanoma patients. Study on melanoma metastasis mechanism would beneficial to the development of novel therapeutic strategies to extend the overall survival of patients.

Certain signaling pathways have been involved in metastasis of cutaneous melanoma, such as integrin alpha V beta 3 signaling (1) and signal transducer and activator of transcription 3 (STAT3) signaling (2). Several indicators and biomarkers associated with metastasis have also been identified in melanoma. LINE1 hypermethylation in peripheral blood of cutaneous melanoma patients is associated with metastasis (3). Increased level of circulating U2 small nuclear RNA fragments indicates metastasis in melanoma patients (4). Chen et al found that Tip60 may regulate melanoma metastasis and could be a potential therapeutic target (5). Thang et al reported that deltex-3-like (DTX3L) stimulates metastasis of melanoma (6). The study by De Semir et al (7) indicated that pleckstrin homology domain-interacting protein (PHIP) is a marker and mediator of melanoma metastasis. Serpin family E member 1 (SERPINE1) expression discriminates site-specific metastasis in human melanoma (8). Galectin-3 expression favors metastasis in murine melanoma (9).

Gene expression profiling is a powerful tool to unveil genes implicated in metastasis of melanoma (10,11). In present study, those gene expression data were combined and more differentially expressed genes (DEGs) were identified via meta-analysis. Feature genes were revealed via support vector machine (SVM) classification. Meanwhile, a SVM classifier was acquired and validated. These findings could advance the understanding about melanoma metastasis and also provide potential therapeutic targets.

 Materials and methods <sec> <title>Gene expression data and pre-treatment

Two gene expression datasets (GSE46517 and GSE7553) were retrieved from Gene Expression Omnibus (GEO) with key words ‘cutaneous melanoma’, ‘homo sapiens’ and ‘metastasis’ by the end of May 13th, 2016. Dataset GSE46517 contained 73 metastatic melanoma samples and 31 primary melanoma samples. Dataset GSE7553 contained 40 metastatic melanoma samples and 14 primary melanoma samples. Moreover, gene expression profiles of 481 skin cutaneous melanoma specimens were downloaded from The Cancer Genome Atlas (TCGA, https://tcga-data.nci.nih.gov/tcga) with keyword ‘cutaneous melanoma’, out of which 183 were metastatic melanoma samples, 71 primary melanoma samples and others not specified.

Gene expression dataset GSE7553 was acquired with Affymetrix platform. Background correction and normalization were performed with package oligo (12) of R. Missing values were filled with median method. Backgroud correction was carried out with MSA method. Normalization was achieved with quantiles method. As for dataset GSE46517 from Agilent platform, probes were mapped into genes. For probes corresponding to a same gene, the averaged value was considered as the expression value of the gene. Log2 transformation and normalization with the median were then applied. Level 3 gene expression data from TCGA were normalized with package limma (13) of R.

 Screening of differentially expressed genes (DEGs)

Meta-analysis of the three gene expression dataset was performed with package MetaDE (14) of R to identify differentially expressed genes (DEGs). This method firstly tested heterogeneity of gene expression value in various platforms with three statistic parameters: tau2, Q value and Qpval. Then it tested differential expression of genes between different groups with statistic parameters P-value and false positive rate (FDR). To ensure the homogeneity of feature genes, tau2=0, Qpval >0.05 and FDR <0.05 were set as the cut-offs.

 Construction of protein-protein interaction (PPI) network

The PPI information was downloaded from Database of Protein, Chemical, and Genetic Interactions (BioGRID, http://thebiogrid.org/) (15), Human Protein Reference Database (HPRD, http://www.hprd.org/) (16) and Database of Interacting Proteins (DIP, http://dip.doe-mbi.ucla.edu/dip/Main.cgi) (17). The protein products of the DEGs were mapped into the whole network and the PPI network for the DEGs was then acquired. The network was then visualized with Cytoscape (18).

 Calculation of betweenness centrality (BC)

Feature genes were screened out from the DEGs using betweenness centrality (BC), which reflected hubness of the node in the PPI network. The BC was calculated as follow:

CB(v)=∑t≠v≠u∈Vσst(v)σst

Where σst is the shortest path from s to t; σst(v) is number of the shortest path from s to t passing through node v; BC score is between 0 and 1, and greater BC score indicates higher degree of hubness.

 Training and validation of SVM classifier

Gene expresson data from TCGA were chosen as the training set. Genes were ranked based upon BC value and top 10 genes were selected out to train SVM classifier. An increment of 10 genes were added into the classifier until metastatic melanomas could be totally separated from primary melanomas. These DEGs were regarded as feature genes. Gene expression datasets GSE46517 and GSE7553 were used as the validation set. Sensitivity (Se), specificity (Sp), positive predictive value (PPV), negative predictive value (NPV) and area under curve (AUC) were calculated to evaluate the SVM classifier.

 Pathway enrichment analysis

Gene ontology (GO) biological pathways and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways related to feature genes were identified using runHyperGO and runHyperKEGG from package EMA of R. The significance was calculated using Fisher's exact test as follow. P-value <0.05 was set as the threshold.

p=1–∑i=0x–1(Mi)(N–MK–i)(NK)

Where N is the total number of gene; M is number of gene in the pathway; K is number of feature gene.

 Results <sec> <title>Differentially expressed genes

According to the criteria, a total of 798 DEGs were revealed from the three gene expression datasets. Heat map of the expression levels of the 798 DEGs is shown in Fig. 1.

 PPI network

Based upon information from BioGRID, HPRD and DIP, a PPI network containing 337 nodes and 466 edges was obtained (Fig. 2). Distribution of degree is shown in Fig. 3. Most genes (232 genes) showed small degree (Log transformed degree <1) while only 1 gene had Log transformed degree >4. Therefore, like most biological networks, this PPI network exhibited scale-free property. These genes with high degree were hub genes and might play important roles in the development of diseases.

 SVM classifier

Top 10 DEGs ranked by BC are listed in Table I: Integrin subunit alpha 4 (ITGA4), activating transcription factor 2 (ATF2), SRY-box 2 (SOX2), keratin 31 (KRT31), cullin 2 (CUL2), replication protein A1 (RPA1), calmodulin 1 (CALM1), FUS RNA binding protein (FUS), angiomotin like 2 (AMOTL2) and ras-related C3 botulinum toxin substrate 1 (RAC1).

As shown in Fig. 4A, the accuracy increased from 85.4 to 96.8% while more genes were included in the classifier. Details are shown in Table II. The accuracy rate reached 96.8% while top 110 DEGs were included in the classifier, and it didn't increase though more genes were added in. Therefore, the top 110 genes were used to construct the SVM classifier. The scatter plot of prediction result for the training dataset is shown in Fig. 4B and the receiver operating characteristic (ROC) curve is shown in Fig. 4C.

The SVM classifier was validated using dataset GSE46517 and GSE7553. The accuracy rates in dataset GSE46517 (Fig. 5A) and GSE7553 (Fig. 5B) were 100 and 94.4%, respectively. Three non-metastatic melanoma samples from dataset GSE7753 were predicted wrong. The accuracy, sensitivity, specificity, PPV, NPV and AUC are listed in Table III. The receiver operating characteristic (ROC) curves are shown in Fig. 5C and D.

 Biological pathways

A total of 11 KEGG pathways (Fig. 6, top) and 13 GO biological pathways (Fig. 6, bottom) were identified in the 110 feature genes. Endometrial cancer, regulation of actin cytoskeleton, focal adhesion, ubiquitin mediated proteolysis, regulation of apoptosis, and regulation of cell proliferation were significantly over-represented in the feature genes.

 Discussion

Three gene expression datasets were downloaded from GEO or TCGA. A total of 798 DEGs were identified in metastatic melanoma compared with non-metastatic melanoma via meta-analysis. A PPI network containing 337 nodes and 466 edges was constructed. BC was calculated for each gene in the network. Top 110 feature genes ranked by BC were included in the SVM classifier. The prediction accuracies for the three datasets were 96.8, 100, and 94.4%, respectively. A total of 11 KEGG pathways and 13 GO biological pathways were significantly over-represented in the 110 feature genes, such as endometrial cancer, regulation of actin cytoskeleton, focal adhesion, ubiquitin mediated proteolysis, regulation of apoptosis, and regulation of cell proliferation. Among the 110 feature genes, several genes have been involved in metastasis of melanoma.

ATF2 is a transcription factor and a member of the leucine zipper family of DNA binding proteins. It confers radiation resistance to human melanoma cells (19). Lau et al further pointed out that transcriptional repression of IFNβ1 by ATF2 confers melanoma resistance to therapy (20). It also promotes melanoma metastasis by suppressing protein fucosylation (21). However, ATF2-derived peptides can inhibit melanoma growth and metastasis (22). Therefore, ATF2 is a critical gene in melanoma metastasis and could be a therapeutic target.

SOX2 is a member of the SRY-related HMG-box (SOX) family of transcription factors involved in the regulation of embryonic development and in the determination of cell fate. It regulates self-renewal and tumorigenicity of human melanoma-initiating cells (23). It also contributes to melanoma cell invasion (24). Downregulation of SOX2 in cancer stem cells suppresses growth and metastasis of lung cancer (25). We speculated that it could be a potential therapeutic target for metastatic melanoma.

RAC1 is a GTPase which belongs to the RAS superfamily of small GTP-binding proteins, which regulate a diverse array of cellular events, including the control of cell growth, cytoskeletal reorganization, and the activation of protein kinases. RAC1 is a critical genes associated with melanoma metastasis. A previous study has reported that a RAC1 mutant (Pro29 to serine, P29S) promotes melanocyte proliferation and migration (26). A recent study indicates that RAC1 P29S regulates PD-L1 expression in melanoma (27). Besides, matrix metallopeptidase 14 modulates melanoma cell dissemination and metastasis through activation of MMP2 and RAC1 (28). Immunity-related GTPase family M member 1 (IRGM1) enhances B16 melanoma cell metastasis through PI3K-RAC1 mediated epithelial mesenchymal transition (29).

Certain genes take parts in metastasis of various cancers and their roles in melanoma metastasis remain to be uncovered. AMOTL2 is related to angiomotin and is a member of the motin protein family. It disrupts apical-basal cell polarity and promotes tumour invasion (30). Raf-1 proto-oncogene (RAF1) activation in gastrointestinal carcinoid cells decreases tumor cell adhesion and thus promotes metastasis (31). Vav guanine nucleotide exchange factor 1 (VAV1) is a member of the VAV gene family. Inhibition of VAV1 can attenuate metastasis of pancreatic cancer (32). Besides, RPA1 is regarded as a useful prognostic indicator in colon cancer (33) and esophageal carcinoma (34). We thought that it might exert prognostic effect in melanoma. Further studies on these genes could expand the knowledge about melanoma metastasis.

In the present study, the biological pathways that were significantly enriched in the metastasis related feature genes in melanoma were screened. To find out the pathways involved in the feature genes in metastatic melanoma, Valsesia et al also conducted a network-guided investigation (35). They connected the genes of somatic copy number alterations in metastatic melanoma with significant biological pathways. The results showed the frequently altered pathways in melanoma, including angiogenesis and melanogenesis, while we found that regulation of apoptosis and regulation of cell proliferation participate in melanoma.

Overall, a multigene SVM predictor of high accuracy was obtained for metastatic melanoma. Several critical genes linked with melanoma metastasis were also revealed. Further researches on the feature genes could improve the understanding about the molecular mechanisms underlying metastasis and provide potential therapeutic targets.

 References BaiJZhangJWuJShenLZengJDingJWuYGongZLiAXuSJWA regulates melanoma metastasis by integrin alphaVbeta3 signalingOncogene2912271237201010.1038/onc.2009.40819946336 CaoHHChuJHKwanHYSuTYuHChengCYFuXQGuoHLiTTseAKInhibition of the STAT3 signaling pathway contributes to apigenin-mediated anti-metastatic effect in melanomaSci Rep621731201610.1038/srep21731269118384766576 De AraújoÉSKashiwabaraAYAchatzMIMoredoLFDe SáBCDupratJPRosenbergCCarraroDMKrepischiACLINE-1 hypermethylation in peripheral blood of cutaneous melanoma patients is associated with metastasisMelanoma Res25173177201510.1097/CMR.0000000000000141256477374352068 KuhlmannJDWimbergerPWilschKFluckMSuterLBrunnerGIncreased level of circulating U2 small nuclear RNA fragments indicates metastasis in melanoma patientsClin Chem Lab Med53605611201510.1515/cclm-2014-106425741740 ChenGChengYTangYMartinkaMLiGRole of Tip60 in human melanoma cell migration, metastasis, and patient survivalJ Invest Dermatol13226322641201210.1038/jid.2012.19322673729 ThangNDYajimaIKumasakaMYIidaMSuzukiTKatoMDeltex-3-like (DTX3L) stimulates metastasis of melanoma through FAK/PI3K/AKT but not MEK/ERK pathwayOncotarget61429014299201510.18632/oncotarget.3742260334504546467 De SemirDNosratiMBezrookoveVDarAAFedermanSBienvenuGVennaSRangelJClimentJMeyer TamgüneyTMPleckstrin homology domain-interacting protein (PHIP) as a marker and mediator of melanoma metastasisProc Natl Acad Sci USA10970677072201210.1073/pnas.1119949109225117203344978 KleinRMBernsteinDHigginsSPHigginsCEHigginsPJSERPINE1 expression discriminates site-specific metastasis in human melanomaExp Dermatol21551554201210.1111/j.1600-0625.2012.01523.x227162553383094 ComodoANLacerda BachiALSoaresMFFrancoMde Paulo Castro TeixeiraVGalectin-3 expression favors metastasis in murine melanomaAdv Biosci Biotechnol045562201310.4236/abb.2013.410A3007 KabbarahONogueiraCFengBNazarianRMBosenbergMWuMScottKLKwongLNXiaoYCordon-CardoCIntegrative genome comparison of primary and metastatic melanomasPLoS One5e10770201010.1371/journal.pone.0010770205207182875381 RikerAIEnkemannSAFodstadOLiuSRenSMorrisCXiYHowellPMetgeBSamantRSThe gene expression profiles of primary and metastatic melanoma yields a transition point of tumor progression and metastasisBMC Med Genomics113200810.1186/1755-8794-1-13184424022408576 CarvalhoBBengtssonHSpeedTPIrizarryRAExploration, normalization, and genotype calls of high-density oligonucleotide SNP array dataBiostatistics8485499200710.1093/biostatistics/kxl04217189563 SmythGKLimma: Linear models for microarray dataStat Biol Health14202005 LiJTsengGCAn adaptively weighted statistic for detecting differential gene expression when combining multiple transcriptomic studiesAnn Appl Stat59941019201110.1214/10-AOAS393 Chatr-AryamontriABreitkreutzBJOughtredRBoucherLHeinickeSChenDStarkCBreitkreutzAKolasNO'DonnellLThe BioGRID interaction database: 2015 updateNucleic Acids Res43Database issueD470DS478201510.1093/nar/gku120425428363 Keshava PrasadTSGoelRKandasamyKKeerthikumarSKumarSMathivananSTelikicherlaDRajuRShafreenBVenugopalAHuman protein reference database-2009 updateNucleic Acids Res37Database issueD767D772200910.1093/nar/gkn89218988627 XenariosIRiceDWSalwinskiLBaronMKMarcotteEMEisenbergDDIP: The database of interacting proteinsNucleic Acids Res28289291200010.1093/nar/28.1.28910592249102387 SmootMEOnoKRuscheinskiJWangPLIdekerTCytoscape 2.8: New features for data integration and network visualizationBioinformatics27431432201110.1093/bioinformatics/btq67521149340 RonaiZYangYMFuchsSYAdlerVSardanaMHerlynMATF2 confers radiation resistance to human melanoma cellsOncogene16523531199810.1038/sj.onc.12015669484842 LauESedyJSanderCShawMAFengYScortegagnaMClapsGRobinsonSChengPSrivasRTranscriptional repression of IFNβ1 by ATF2 confers melanoma resistance to therapyOncogene3457395748201510.1038/onc.2015.22257286764558399 LauEFengYClapsGFukudaMNPerlinaADonnDJilaveanuLKlugerHFreezeHHRonaiZAThe transcription factor ATF2 promotes melanoma metastasis by suppressing protein fucosylationSci Signal8ra124201510.1126/scisignal.aac6479266455814818095 BhoumikAGangiLRonaiZInhibition of melanoma growth and metastasis by ATF2-derived peptidesCancer Res6482228230200410.1158/0008-5472.CAN-04-071415548688 SantiniRPietrobonoSPandolfiSMontagnaniVD'AmicoMPenachioniJYVinciMCBorgognoniLSteccaBSOX2 regulates self-renewal and tumorigenicity of human melanoma-initiating cellsOncogene3346974708201410.1038/onc.2014.71246819554180644 GirouardSDLagaACMihmMCScolyerRAThompsonJFZhanQWidlundHRLeeCWMurphyGFSOX2 contributes to melanoma cell invasionLab Invest92362370201210.1038/labinvest.2011.18822184093 XiangRLiaoDChengTZhouHShiQChuangTSMarkowitzDReisfeldRALuoYDownregulation of transcription factor SOX2 in cancer stem cells suppresses growth and metastasis of lung cancerBr J Cancer10414101417201110.1038/bjc.2011.210214680473101944 KrauthammerMKongYHaBHEvansPBacchiocchiAMcCuskerJPChengEDavisMJGohGChoiMExome sequencing identifies recurrent somatic RAC1 mutations in melanomaNat Genet4410061014201210.1038/ng.2359228422283432702 VuHLRosenbaumSPurwinTJDaviesMAAplinAERAC1 P29S regulates PD-L1 expression in melanomaPigment Cell Melanoma Res28590598201510.1111/pcmr.12392261767074675336 ShaverdashviliKWongPMaJZhangKOsmanIBedogniBMT1-MMP modulates melanoma cell dissemination and metastasis through activation of MMP2 and RAC1Pigment Cell Melanoma Res27287296201410.1111/pcmr.1220124387669 TianLLiLXingWLiRPeiCDongXFuYGuCGuoXJiaYIRGM1 enhances B16 melanoma cell metastasis through PI3K-Rac1 mediated epithelial mesenchymal transitionSci Rep512357201510.1038/srep12357262029104512008 MojallalMZhengYHultinSAudebertSvan HarnTJohnssonPLenanderCFritzNMiethCCorcoranMAmotL2 disrupts apical-basal cell polarity and promotes tumour invasionNat Commun54557201410.1038/ncomms555725080976 GreenblattDYKunnimalaiyaanMChenHRaf-1 activation in gastrointestinal carcinoid cells decreases tumor cell adhesionAm J Surg193331335200710.1016/j.amjsurg.2006.09.016173205291838566 RazidloGLMagnineCSlettenACHurleyRMAlmadaLLFernandez-ZapicoMEJiBMcNivenMATargeting pancreatic cancer metastasis by inhibition of Vav1, a driver of tumor cell invasionCancer Res7529072915201510.1158/0008-5472.CAN-14-3103259773354506209 GivalosNGakiopoulouHSkliriMBousboukeaKKonstantinidouAEKorkolopoulouPLeloudaMKouraklisGPatsourisEKaratzasGReplication protein A is an independent prognostic indicator with potential therapeutic implications in colon cancerMod Pathol20159166200710.1038/modpathol.380071917361204 DahaiYSanyuanSHongLDiZChongZA relationship between replication protein A and occurrence and prognosis of esophageal carcinomaCell Biochem Biophys67175180201310.1007/s12013-013-9530-y23420526 ValsesiaARimoldiDMartinetDIbbersonMBenaglioPQuadroniMWaridelPGaillardMPidouxMRapinBNetwork-guided analysis of genes with altered somatic copy number and gene expression reveals pathways commonly perturbed in metastatic melanomaPLoS One6e18369201110.1371/journal.pone.0018369214946573072964

Heat map of the expression levels of the 798 DEGs. ‘M’ means metastatic cutaneous melanoma; ‘N’ means non-metastatic cutaneous melanoma. X-axis is samples and Y-axis is genes expression level. Green color means higher gene expression level; red color means lower gene expression level. DEGs, differentially expressed genes.

The protein-protein interaction network of the differentially expressed genes (DEGs). Up-regulated genes in metastatic melanoma are in red, down-regulated genes in metastatic melanoma are in green.

Distribution of degree. X-axis is Log transformed degree and Y-axis is number of node.

Accuracy rate with various numbers of differentially expressed genes (A) for the training dataset. Scatter plot showing prediction result of the training dataset (B). Metastatic cutaneous melanoma samples are in red and non-metastatic cutaneous melanoma samples are in black. The receiver operating characteristic (ROC) curve of for the training dataset (C).

Scatter plots and ROC curves for dataset GSE46517 (A and C) and GSE7553 (B and D). Metastatic cutaneous melanoma samples are in red and non-metastatic cutaneous melanoma samples are in black. ROC, receiver operating characteristic.

KEGG pathways and GO biological pathways significantly over-represented in the 110 feature genes. X-axis indicates number of gene and color of the bar indicates-lg (P-value). KEGG, Kyoto Encyclopedia of Genes and Genomes; GO, gene ontology.

Top 10 genes ranked by betweenness centrality.

Gene BC Exp Degree P-value FDR Q Qp tau2
ITGA4 2.06×10−1 1 25 1.00×10−20 9.22×10−19 1.39E+00 5.00×10−1 0
ATF2 1.47×10−1 1 16 2.46×10−6 1.78×10−4 1.54E+00 4.64×10−1 0
SOX2 1.43×10−1 0 15 3.73×10−4 6.85×10−3 5.21×10−2 9.74×10−1 0
KRT31 1.36×10−1 0 10 1.00×10−20 9.22×10−19 4.96×10−1 7.80×10−1 0
CUL2 1.22×10−1 1 16 2.69×10−4 5.49×10−3 1.19×10−1 9.42×10−1 0
RPA1 1.10×10−1 0 16 4.07×10−3 3.69×10−2 6.28×10−1 7.31×10−1 0
CALM1 9.78×10−2 0 13 7.63×10−4 1.12×10−2 8.76×10−1 6.45×10−1 0
FUS 9.26×10−2 0 12 6.57×10−6 3.27×10−4 1.95E+00 3.77×10−1 0
AMOTL2 9.25×10−2 0   4 9.86×10−6 4.24×10−4 8.40×10−1 6.57×10−1 0
RAC1 8.70×10−2 0   7 2.38×10−3 2.54×10−2 1.13E+00 5.69×10−1 0

BC, betweenness centrality; Exp, expression; FDR, false positive rate. ‘1’ in column Exp indicates up-regulation in metastatic melanoma compared with non-metastatic melanoma while ‘0’ indicates down-regulation in metastatic melanoma.

Prediction results of the SVM classifier in the training dataset with varying numbers of differentially expressed genes.

SVM classifier Correct prediction Wrong prediction Accuracy rate (%) Error rate (%)
Top10 217 37 85.4 14.6
Top20 225 29 88.6 11.4
Top30 230 24 90.6   9.4
Top40 231 23 90.9   9.1
Top50 232 22 91.3   8.7
Top60 234 20 92.1   7.9
Top70 236 18 92.9   7.1
Top80 238 16 93.7   6.3
Top90 240 14 94.5   5.5
Top100 242 12 95.3   4.7
Top110 246   8 96.9   3.1
Top120 246   8 96.9   3.1
Top130 246   8 96.9   3.1

Prediction results of the SVM classifier in the 3 independent datasets.

Dataset No. of sample Accuracy Sensitivity Specificity PPV NPV AUROC
TCGA 254 0.968 1.000 0.987 0.958 1.000 0.991
GSE46517 128 1.000 1.000 1.000 1.000 1.000 1.000
GSE7753 54 0.944 1.000 0.886 0.942 1.000 0.947

PPV, positive predictive value; NPV, negative predictive value; AUROC, area under receiver operating characteristic curve.