Identification of Core Genes and Pathways in Melanoma Metastasis via Bioinformatics Analysis

Metastasis is the leading cause of melanoma-related mortality. Current therapies are rarely curative for metastatic melanoma, revealing the urgent need to identify more effective preventive and therapeutic targets. This study aimed to screen the core genes and molecular mechanisms related to melanoma metastasis. A gene expression profile, GSE8401, including 31 primary melanoma and 52 metastatic melanoma clinical samples, was downloaded from the Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) between melanoma metastases and primary melanoma were screened using GEO2R tool. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genome (KEGG) analyses of DEGs were performed using the Database for Annotation Visualization and Integrated Discovery (DAVID). The Search Tool for the Retrieval of Interacting Genes (STRING) and Cytoscape with Molecular Complex Detection (MCODE) plug-in tools were utilized to detect the protein–protein interaction (PPI) network among DEGs. The top 10 genes with the highest degrees of the PPI network were defined as hub genes. In the results, 425 DEGs, including 60 upregulated genes and 365 downregulated genes, were identified. The upregulated genes were enriched in ECM–receptor interactions and the regulation of actin cytoskeleton, while 365 downregulated genes were enriched in amoebiasis, melanogenesis, and ECM–receptor interactions. The defined hub genes included CDK1, COL17A1, EGFR, DSG1, KRT14, FLG, CDH1, DSP, IVL, and KRT5. In addition, the mRNA and protein levels of the hub genes during melanoma metastasis were verified in the TCGA database and paired post- and premetastatic melanoma cells, respectively. Finally, KRT5-specific siRNAs were utilized to reduce the KRT5 expression in melanoma A375 cells. An MTT assay and a colony formation assay showed that KRT5 knockdown significantly promoted the proliferation of A375 cells. A Transwell assay further suggested that KRT5 knockdown significantly increased the cell migration and cell invasion of A375 cells. This bioinformatics study provided a deeper understanding of the molecular mechanisms of melanoma metastasis. The in vitro experiments showed that KRT5 played the inhibitory effects on melanoma metastasis. Therefore, KRT5 may serve important roles in melanoma metastasis.


Introduction
Melanoma, as the most common and grim malignant skin cancer, has an incidence that has unfortunately been steadily increasing in the last 40 years, and its incidence is ascending faster than that of any other solid tumor. Due to the extraordinary predisposition of melanoma to spread and its rapid progression toward metastasis, patients who develop metastasis almost always have an incurable disease, with only 6-9 months of median survival time, a 15% 3-year survival rate, and only a 4.6% 5-year survival rate [1].
Metastasis is a serious event in the clinic, leading to the majority of deaths of melanoma patients [2]. Despite the landscape of genetic alterations and elaborate molecular mechanisms discovered in melanoma, little information about the underlying biology that drives its metastasis has been fully elucidated [3]. Therefore, a more in-depth understanding of the metastatic melanoma process is urgently needed, with the aim to develop specific therapies for improving current therapy and reducing the mortality of melanoma patients.
Historically, numerous genes have been identified and are expected to be targets for preventing melanoma metastasis. For instance, PGC1A encodes PGC1α, which is a metabolic transcriptional coactivator and suppresses melanoma metastasis by protecting against oxidative stress [4]. The overexpression of KISS1 inhibits the metastasis of C8161 melanoma cells [5]. BRIC5-encoded survivin promotes melanoma metastasis through the Akt-dependent upregulation of α5 integrin [6]. Our previous work also uncovered the oncogenic functions and regulatory mechanisms of NOL7 and WDR74 in melanoma [7,8]. However, these findings were hardly sufficient to develop a complete overview of melanoma metastasis. Furthermore, the peculiarity of the high heterogeneity and clinical phenotype of melanoma cells imply a complicated regulatory mechanism of cancer metastasis [9].
Microarray-based gene expression analysis is a powerful and promising tool for functional genomics research, especially for the accurate and comprehensive analysis of complex networks involved in biological processes [10,11]. In our present study, we chose the GSE8401 profile from the Gene Expression Omnibus (GEO) database and utilized the GEO2R tool to identify the differentially expressed genes (DEGs). Subsequently, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genome (KEGG) pathway analyses were performed, and a protein-protein interaction (PPI) network was constructed. In general, the hub genes, including KRT5, IVL, and DSP, and key pathways associated with melanoma metastasis defined in the present study may provide new insights into clinical melanoma treatment. Additionally, the effects of KRT5 on cell proliferation and cell metastatic behaviors of melanoma were determined by in vitro experiments.

Identification of DEGs
The GSE8401 profile, including 31 primary melanoma samples and 52 melanoma metastasis samples, was submitted to the GEO2R online analysis tool. Using an adjusted p value < 0.05 and |log2 (Fold change)| > 1.5 as the cutoff criteria, a total of 425 DEGs, including 60 upregulated genes and 365 downregulated genes, were picked out. The heat map of the full range of genes and volcano plot of DEGs are shown in Figure 1. A list of DEGs is shown in Table S1. The top 15 significantly upregulated and downregulated genes are listed in Table 1. PSPH was the most significantly upregulated gene, and S100A7 was the most significantly downregulated gene.

GO Function and KEGG Pathway Enrichment Analyses of DEGs
The top five enriched terms identified in each GO category using Database for Annotation Visualization and Integrated Discovery (DAVID) software are shown in Table S2. The upregulated genes were significantly found to participate in the formation of cellular components, including condensed chromosome kinetochore, chromosome and centromeric region; condensed chromosome and centromeric region; and condensed chromosome, chromosomal region. The biological processes included sister chromatid segregation, mitotic cell cycle process, nuclear chromosome segregation, mitotic nuclear division, and cell cycle process. The molecular functions included small molecule binding, adenyl ribonucleotide binding, adenyl nucleotide binding, ATP binding, and nucleotide binding. In addition, the downregulated genes were mainly involved in cellular components, such as the extracellular region, extracellular region part, extracellular exosome, extracellular vesicle, and extracellular organelle; biological processes including epidermis development, skin development, keratinocyte differentiation, epidermal cell differentiation, and epithelial cell differentiation; and molecular functions including structural molecule activity, endopeptidase inhibitor activity, structural constituent of cytoskeleton, endopeptidase regulator activity, and peptidase inhibitor activity.  A KEGG pathway analysis was then conducted. As shown in Table 2, the upregulated genes were enriched in ECM-receptor interactions, progesterone-mediated oocyte maturation, regulation of actin cytoskeleton, and metabolic pathways, while the downregulated genes were enriched in amoebiasis, melanogenesis, and ECM-receptor interactions.

Hub Genes and Module Screening from PPI Network
Based on the information in the Search Tool for the Retrieval of Interacting Genes (STRING) protein query from the public databases of the DEGs, the top 10 genes with the highest degrees of connectivity were selected and defined as hub genes ( Figure 2). Additionally, 250 nodes and 751 edges in the PPI network among those DEGs were analyzed using the Molecular Complex Detection (MCODE) plug-in in Cytoscape, and the top 2 significant modules were filtered ( Figure 3). Based on the GO function and KEGG pathway analyses, these two modules were principally associated with cell cycle, ECM-receptor interactions, focal adhesion, and the PI3K-Akt signaling pathway.

Validation of Hub Gene Expression between Primary Melanoma and Metastatic Site in TCGA Database
To confirm the reliability and accuracy of the results through the above bioinformatics analysis, we verified the mRNA level of the hub genes in the TCGA database. As shown in Figure 4, The results showed that gene expression levels of COL17A1, DSG1, KRT14, FLG, CDH1, DSP, IVL, and KRT5 had significant decreases in metastatic melanoma (n = 71) compared with primary melanoma (n = 105). However, the gene expression level of CDK1 and EGFR had no significant differences between metastatic melanoma and primary melanoma.

Validation of Protein Expression of Hub Genes in Paired Premetastatic and Postmetastatic Melanoma Cells
Based on the above data ( Figure 4), the gene expressions of CDH1, KRT5, COL17A1, KRT14, IVL, DSP, DSG1, and FLG show the consistent downregulation in metastatic melanoma compared with primary melanoma. To further confirm the protein expression levels of these eight hub genes in the development of melanoma metastasis, paired premetastatic melanoma cells (A375) and postmetastatic melanoma cells (A375M) were established through an experimental animal model of melanoma metastasis (see the Materials and Methods section). Next, western blotting was carried out to measure the protein level of those eight hub genes. As shown in Figure 5, the expression of these eight hub genes was significantly downregulated in metastatic melanoma cells compared with the primary melanoma cells (p < 0.01).

KRT5 Knockdown Promotes Cell Proliferation, Migration, and Invasion of Melanoma
To further prove the reliability and accuracy of this bioinformatics analysis, one of the hub gene, KRT5, was selected for further biological experiments analysis. The reduction efficiency was determined by siRNA knockdown of KRT5, followed by RT-qPCR and western blotting ( Figure 6A,B). The MTT experiments demonstrated that the proliferative rate of A375 cells was significantly increased upon KRT5 knockdown ( Figure 6C). The colony formation experiments revealed that the number of colonies in the KRT5-siRNA transfection group was significantly higher compared with that in the control group ( Figure 6D). The Transwell assay indicated that the migration and invasion abilities were significantly enhanced following KRT5 knockdown in A375 cells ( Figure 6E,F). These results suggested the inhibitory effect of KRT5 in melanoma metastasis.

Discussion
Melanoma is the most malignant type of skin cancer, and its metastasis remains essentially incurable because the mutated genes and underlying molecular mechanisms are poorly uncovered [12]. In the present study, the gene expression profile of GSE8401 was analyzed, and 425 DEGs, including 60 upregulated genes and 365 downregulated genes, were filtered between 31 primary melanoma samples and 52 metastatic melanoma samples, which were obtained from the clinical melanoma patients.
The GO function analysis revealed that upregulated DEGs mainly participated in small molecule binding, nucleotide binding, and cell cycle processes, while the downregulated DEGs were involved in extracellular region, extracellular region part, extracellular organelle, exosome, and vesicles. Based on these results, cell mitosis and malignant proliferation were activated, whereas the interaction with the extracellular environment was suppressed during the metastatic process.
The KEGG pathway analysis revealed that upregulated DEGs were mainly enriched in ECM-receptor interactions, progesterone-mediated oocyte maturation, regulation of the actin cytoskeleton, and metabolic pathways and that downregulated DEGs were enriched in amoebiasis, melanogenesis, and ECM-receptor interactions. Both upregulated genes and downregulated genes were enriched in ECM-receptor interactions, indicating that this pathway is pronouncedly valuable to melanoma metastasis. Cross-talk between the ECM and melanoma metastasis is well known to be commonly elaborate in previous reports [13,14]. Cancer metastasis involves multiple complex processes that are critically influenced by ECM components [15]. Various ECM-related proteins are significantly dysregulated during the progression of cancer, causing both biochemical and biomechanical changes that together promote cancer metastasis [16]. Herein, the expressions of ITGA4, SPP1, and ITGB3 were increased, while the COL6A2, LAMC2, LAMB3, SDC1, ITGB4, LAMC3, COMP, and LAMA3 expressions were decreased. It prompted the genes ITGA4, SPP1, and ITGB3 to promote the interactions between melanoma cells and the ECM and thus facilitated melanoma metastasis. On the contrary, the genes COL6A2, LAMC2, LAMB3, SDC1, ITGB4, LAMC3, COMP, and LAMA3 might play opposite effects since the regulation of actin cytoskeleton, metabolic pathways, and melanogenesis are the key pathways involved in melanoma metastasis [17][18][19]. The regulation of action cytoskeleton is mainly related to cell migration. The dynamic actin cytoskeleton spatially and temporally mediates protrusion, adhesion, contraction, and retraction from the active cell. The variability of action cytoskeleton of cancer cells confers the aggressive phenotype, such as the EMT. Several genes, such as FAK and c-fos, were attested to regulate the actin cytoskeleton to promote melanoma metastasis [20,21]. Dysregulation of metabolic pathways in cancer cells results in the metabolic reprogramming in cancer, leading to enhanced substance uptake to supply the energy production and biosynthesis. It was declared that metabolic remodeling is pivotal for melanoma cells to adapt to tumor microenvironment and to maintain the growth and dissemination of melanoma cells [22]. Obviously, melanogenesis is essential for melanoma development. Deregulated melanogenesis potentially contributes to more aggressive behaviors of melanoma cells [23].
Direct evidence of the participation of progesterone-mediated oocyte maturation and amoebiasis in melanoma metastasis has so far not been illustrated. Rare studies in the literature have shown the interaction between cancer metastasis and progesteronemediated oocyte maturation/amoebiasis. Mood et al. showed that progesterone-mediated Xenopus oocytes maturation pathway was related to the G2/M transition of oocytes [24]. In the current study, CDK1, MAD2L1, and BUB1, which are associated with the cell cycle, were identified to enrich progesterone-mediated oocyte maturation. This implies that the progesterone-mediated oocyte maturation pathway might participate in melanoma metastasis by influencing the cell cycle. In addition, Nelis et al. presented that cancer cells and amebic trophozoites follow the same metastatic route for the liver and other organs [25]. Cancer-related genes that are adhesion/migration-related genes, including LAMB3, LAMA3, LAMC2, and LAMC3, and immune-related genes, including SERPINB3, SERPINB4, and IL1R2, as well as the proliferation-related genes, including ARG1 and SERPINB2, which were enriched in the amoebiasis pathway in this analysis. This implies that the amoebiasis pathway might be involved in melanoma metastasis through these genes. These pathways screened in this study provided promising targets for new drug intervention to fight melanoma.
Most importantly, 10 hub genes were defined. We verified the mRNA expression of those hub genes in various samples of primary melanoma and metastatic melanoma, which were obtained from the TCGA dataset. The overall direction in the TCGA dataset yielded concordance with the results of bioinformatics analysis upon GSE8401, except for EGFR and CDK1. It appears that CDH1, KRT5, COL17A1, KRT14, IVL, DSP, DSG1, and FLG may be potentially useful biomarkers in melanoma metastasis. To further confirm this hypothesis, we tested the protein levels of these eight genes in paired metastatic melanoma cells and primary melanoma cells. Consistently, the protein expression levels of CDH1, KRT5, COL17A1, KRT14, IVL, DSP, DSG1, and FLG were significantly downregulated in metastatic melanoma compared with primary melanoma.
Among these hub genes, we identified the attenuated expressions of CDH1, COL17A1, DSG1, KRT14, and FLG in melanoma metastases compared with primary melanoma via bioinformatics analysis upon the GSE8401 and TCGA datasets. This implies that these genes may suppress melanoma metastasis. There have been several articles that revealed the expression pattern and functional role of those genes in melanoma development. CDH1 encodes E-cadherin, a key biomarker of epithelial-mesenchymal transition (EMT) in cancer cells and a mediator of cell-cell adhesion in epithelial tissues [26]. Loss of Ecadherin promotes the aggressive behaviors of melanoma cells via constitutively active snail expression during the metastasis process [20,27,28]. COL17A1 encodes the collagen alpha-1 (XVII) chain, which is responsible for the adhesion of cells and matrix. Kai Tao et al. found that COL17A1 was decreased in melanoma compared with normal tissues through bioinformatics analysis, while an inverse result was observed in the qRT-PCR validation [29]. DSG1-encoded desmoglein 1 is the component of intercellular desmosome junctions and is involved in cell-cell adhesion [30]. Ji et al. declared that DSG1 might have a potential value for the prognosis and treatment of melanoma, and Herlyn et al. detected that inhibiting DSG1 contributed to melanoma metastasis [31,32]. KRT14 encodes keratin 14, which is an epithelial proliferative marker [33]. It was reported that KRT14 potentially facilitates melanoma tumorigenesis [34,35]. FLG encodes Filaggrin, which is a functional protein in the epidermis and plays a critical role in skin homeostasis. FLG was observed to participate in melanoma development [31,36].
It is noticeable that the functional characteristics of KRT5, IVL, and DSP in melanoma have not been fully explored, especially in melanoma metastasis. KRT5 encodes keratin 5, which belongs to the keratin family, which are intermediate filament proteins. Keratin 5 is responsible for the structural integrity of epithelial cells and contributes to cell polarization, cytoskeleton regulation, and protein translation [37]. Previous research has confirmed that keratin 5 is a stem cell marker in breast cancer and is associated with cancer recurrence and chemotherapy resistance in ovarian cancer [37,38]. IVL encoding the involucrin is a marker of keratinocyte terminal differentiation and maintains the morphological characteristics of the epidermis [39]. Previous findings have suggested that involucrin is a biomarker of EMT in squamous cell carcinoma [40]. Another group also demonstrated that involucrin might be involved in breast cancer, cervical cancer, and oral cancer [41,42]. DSP encodes desmoplakin, which is a major high molecular weight protein of desmosomes. It was shown that DSP functions as a tumor suppressor in cancer migration [43,44], whereas researchers also indicated that the deficiency of desmoplakin induced loss of the epithelial phenotype and acquisition of aggressive phenotype and thus facilitated melanoma metastasis [45]. The explicit function and regulatory mechanism of DSP in melanoma metastasis are not yet understood.
In the present study, the downregulation of the KRT5, IVL, and DSP expressions implied inhibitory effects of these genes in melanoma metastasis. To further test this hypothesis, KRT5 was selected for further biological analysis as there are few studies on KRT5 in the literature and the mechanism related to melanoma metastasis is unclear. We suppressed the KRT expression in melanoma A375 cells through KRT5-specfic siRNAs and then performed the in vitro experiments to clarify the effects of KRT5 on cell proliferation, migration, and invasion of melanoma. The experimental results showed that KRT5 knockdown can significantly enhance the cell metastasis capacities in melanoma, including cell proliferation, migration, and invasion.
In this study, we highlighted that KRT5, IVL, and DSP may serve as novel targets to suppress melanoma metastasis. The possible effect of KRT5 on melanoma metastasis were preliminarily investigated through in vitro biological experiments. Further studies about the detailed biological function and regulatory mechanisms of these genes in melanoma metastasis are needed in the future.

Microarray Data
The gene expression dataset GSE8401 (GDS3966) (https://www.ncbi.nlm.nih.gov/ geo/query/acc.cgi?acc=GSE8401, 3 January 2022), based on the Agilent GPL96 platform ((HG-U133A) Affymetrix Human Genome U133A Array), was obtained from GEO, a free and publicly available database [46]. The GSE8401 dataset includes 83 clinical samples, containing 31 primary melanomas and 52 melanoma metastases. The heat map of the full expression range of genes was acquired using the GEO online data analysis tool (https://www.ncbi.nlm.nih.gov/sites/GDSbrowser, 3 January 2022). The volcano plot of the full expression range of genes was established using Microsoft Excel.

Identification of the DEGs
GEO2R (https://www.ncbi.nlm.nih.gov/geo/geo2r/, 3 January 2022), an interactive web tool, allows users to compare different groups of clinical cancer samples in a GEO series to identify genes that were differentially expressed across different experimental conditions [47]. We applied this tool to detect the DEGs between primary melanomas and melanoma metastases in GSE8410. The results are presented in a table of genes ordered by significance. The adjusted p-value was calculated to reduce the false positive rate, applying the Benjamini and Hochberg false discovery rate method by default. Genes within the cutoff criteria of an adjusted p-value < 0.05 and |log2 (Fold change)| ≥ 1.5 were designated DEGs.

GO Function and KEGG Pathway Analyses of DEGs
Gene ontology (GO) analysis is a commonly used and productive method for annotating genes and gene products and for identifying the biological characteristics of high-throughput genome or transcriptome data [48]. The Kyoto Encyclopedia of Genes and Genomes (KEGG) is an open and collective database integrating genomes, biological pathways, diseases, drugs, and chemical substances. DEGs were subjected to the Database for Annotation Visualization and Integrated Discovery (DAVID, https://david.ncifcrf.gov/, 3 January 2022), an online bioinformatics tool, to interpret the GO functions and enriched KEGG pathways and to visualize the biological processes (BP), molecular functions (MF), cellular components (CC) and pathways of those DEGs [49]. A p-value < 0.05 and FDR < 0.05 were set as the cutoff criteria.

Construction of PPI Network and Module Analysis
The Search Tool for the Retrieval of Interacting Genes (STRING) is an online tool that assesses protein-protein interaction (PPI) network information [50]. STRING (version 10.5) was used to evaluate the potential PPI relationships among those DEGs. Only experimentally validated interactions with a combined score ≥0.4 were selected as significance. The PPI network was constructed and visualized using Cytoscape software 3.6.0 [51]. The molecular complex detection (MCODE) plug-in in Cytoscape was used to screen the modules of the PPI network. The inferred modules used the default settings with the degree cutoff = 2, node score cutoff = 0.2, K-core = 2, and max depth = 100. Additionally, the hub genes were mapped into STRING with a confidence score ≥0.4 and a maximum number of interactors ≤ 5. The KEGG pathway analysis of the genes in each module was performed using DAVID.

Definitions of Hub Genes
Based on the information in the STRING protein query and degree analysis of the PPI among DEGs using Cytoscope software, the top 10 genes with the highest degrees were defined as hub genes.

Validation of Hub Gene Expression in TCGA Database
The expression level of identified hub genes between primary melanoma and metastatic site was validated using TCGA data, which contains 105 primary melanomas and 71 metastatic melanomas (https://tcga-data.nci.nih.gov/tcga/, 3 January 2022). The comparison between the two datasets was performed with the t-test. A p-value < 0.05 was considered significant.

Cell, Cell Culture, and Postmetastatic Cell Line Establishment
The human melanoma A375 cell line was purchased from the Cell Resource Center of Shanghai Institute for Biological Sciences (Chinese Academy of Sciences, Shanghai, China) and used as the premetastatic parental cell line. The post-metastatic melanoma cell line, A375M, was derived from the pulmonary metastatic nodules of pre-metastasized parental A375 cells via trypsinization as described previously [7,52]. Briefly, A375 cells (3 × 10 5 cells in 0.1 mL saline solution per mouse) were intravenously injected into the tails of 6-8-week-old BALB/C nude mice (Slac Animal Inc., Shanghai, China). After 7 weeks, all mice were sacrificed, and the pulmonary metastatic nodules were stripped to obtain monoplast via trypsinization. Subsequently, these monoplasts were cultured in vitro to establish the A375M cell line.
All cell lines were regularly subjected to mycoplasma testing. The A375 and A375M cell lines were characterized by Genetic Testing Biotechnology Corporation (Suzhou, China) using short tandem repeat (STR) markers. The animal studies were performed following the animal protocol and procedures approved by the Institutional Animal Care and Use Committee (IACUC) of Fuzhou University. The identification code of this animal study is protocol #2019-SG-014, which are consistent with the AAALAS guidelines. All possible efforts were made to minimize animal suffering and sacrifice.

siRNA-Mediated Knockdown of KRT5 in A375 Cells
A375M cells at 70-80% confluency were transfected in OPTI-MEM medium (Invitrogen, Carlsbad, CA, USA) with the indicated siRNA duplexes using lipofectamine 3000 (Invitrogen, Carlsbad, CA, USA). After 6 h of incubation, the transfection medium was removed. Cells were washed with phosphate-buffered solution (PBS) and then cultured in complete medium for 24 to 48 h before further experiments. The reduction efficiency of KRT5 in A375M cells was detected by RT-qPCR and western blotting. The KRT5-specific siRNA duplexes were synthesized by Sangon Biotech (Shanghai, China), and the sequences were listed as follows, 5 -3 : CAUCUCUGUUGUCACAAGCAGUGUU.

Cell Proliferation Assay
For cell proliferative rate assay, target cells were plated in 96-well plates (5000 cells per well). Then, cells were incubated with MTT solution (JT343, Dingguochangsheng Bio, Bei-jing, China) for 4 h, and absorbance at 490 nm was measured by a microplate reader (Tecan Infifinite1 200Pro, Hombrechtikon, Switzerland) daily for 5 days. For cell colony formation assay, target cells (10,000 cells/well) were plated in 6-well plates and cultured for 14 days. Cell colonies were staining with crystal violet solution (DC079, Dingguochangsheng Bio, Beijing, China) and quantified by a Canon scanner (Shanghai, China).

Transwell Assay
Transwell assay was performed as described previously with minor modifications [8,54]. Target cells were suspended in culture medium containing 1% BSA but without FBS. These cells (10 5 cells/200 µL culture medium per group) were plated into the top chamber of Transwell chambers (for migration) or chamber coated with Matrigel (Corning, NY, USA, for invasion). A total of 800 µL culture mediums containing 20% FBS was added into the lower chamber of Transwell chambers. Twenty-four hours later, the number of migrated or invaded cells was evaluated.

Statistical Analysis
The data for all experiments were managed using GraphPad Prism software 8.0 and are represented as the means ± s.d. A paired t-test was used for two-group comparisons. A p-value < 0.05 was considered statistically significant.

Conclusions
Collectively, our bioinformatics analysis presented here identified the DEGs and hub genes involved in melanoma metastasis, which might have important roles in melanoma metastasis. A total of 425 DEGs and 10 hub genes were defined, and the enrichment analysis suggested that an interaction with ECM may play a dominant role in melanoma metastasis. Overall, KRT5, IVL, and DSP might represent potential functions for the prevention and treatment of melanoma metastasis. Notably, KRT5 was confirmed to play inhibitory effect on melanoma metastasis. The findings of this study may contribute to the more profound elucidation of mechanisms of melanoma metastasis. However, further verification experiments are necessary to confirm the results of these analyses.

Conflicts of Interest:
The authors declare no conflict of interest.