- Research article
- Open Access
Screening of gene markers related to the prognosis of metastatic skin cutaneous melanoma based on Logit regression and survival analysis
BMC Medical Genomics volume 14, Article number: 96 (2021)
Bioinformatics was used to analyze the skin cutaneous melanoma (SKCM) gene expression profile to provide a theoretical basis for further studying the mechanism underlying metastatic SKCM and the clinical prognosis.
We downloaded the gene expression profiles of 358 metastatic and 102 primary (nonmetastatic) CM samples from The Cancer Genome Atlas (TCGA) database as a training dataset and the GSE65904 dataset from the National Center for Biotechnology Information database as a validation dataset. Differentially expressed genes (DEGs) were screened using the limma package of R3.4.1, and prognosis-related feature DEGs were screened using Logit regression (LR) and survival analyses. We also used the STRING online database, Cytoscape software, and Database for Annotation, Visualization and Integrated Discovery software for protein–protein interaction network, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses based on the screened DEGs.
Of the 876 DEGs selected, 11 (ZNF750, NLRP6, TGM3, KRTDAP, CAMSAP3, KRT6C, CALML5, SPRR2E, CD3G, RTP5, and FAM83C) were screened using LR analysis. The survival prognosis of nonmetastatic group was better compared to the metastatic group between the TCGA training and validation datasets. The 11 DEGs were involved in 9 KEGG signaling pathways, and of these 11 DEGs, CALML5 was a feature DEG involved in the melanogenesis pathway, 12 targets of which were collected.
The feature DEGs screened, such as CALML5, are related to the prognosis of metastatic CM according to LR. Our results provide new ideas for exploring the molecular mechanism underlying CM metastasis and finding new diagnostic prognostic markers.
Skin cutaneous melanoma (SKCM) is a common skin malignancy with poor prognosis due to aggressiveness and metastasis . In recent years, the incidence of SKCM has significantly increased, and the survival rate of patients remains poor; the 5-year survival rate of metastatic SKCM patients is less than 5% . As metastasis is an important cause of poor prognosis of SKCM, it is necessary to identify its underlying molecular mechanisms and determin its molecular biomarkers.
Previous studies have focused on molecular markers related to SKCM metastasis and prognosis. For example, S100 has been clearly identified as a diagnostic marker of SKCM, and in terms of SKCM prognosis, the most studied diagnostic marker is S100B. Other studies have reported that high S100B expression indicates tumor recurrence and metastasis [3, 4]. Da Forno et al.  showed that high expression of Wnt-5a indicates an increase in SKCM aggressiveness, distant metastases, and poor prognosis. Melanoma inhibitory protein is also considered a diagnostic marker of SKCM metastasis and poor prognosis . Ci et al.  showed that CDCA8 overexpression promotes the malignant progression of SKCM and leads to poor prognosis. Yang et al.  found that STK26, KCNT2, and CASP12 expression is correlated with the prognosis of SKCM using weighted gene co-expression network analysis. However, the mechanism underlying metastasis of a nevus into SKCM is still unclear, and further research is urgently needed.
Logit regression (LR) is popular among medical practitioners because of its interpretability and ease of application without the need for a computer . Zejnullahu et al.  examined the prevalence and risk factors of postpartum depression using LR. Sufriyana et al.  used LR to explore the predictive performances for pregnancy care to inform clinicians’ decision making. Sokou et al. developed and validated a prediction model for clinical variables. These reports suggested that LR model for predicting the risk assessment is regularly used in medicine to guide management decisions.
In this study, we downloaded gene expression detection data of 460 (358 metastatic and 102 primary [nonmetastatic]) SKCM patients from The Cancer Genome Atlas (TCGA) database. The limma package of R3.4.1 was used to screen for differentially expressed genes (DEGs), the survival package of R3.4.1 was used for single- and multifactor Cox regression analysis, and feature DEGs were selected using LR. Next, we constructed a protein–protein interaction (PPI) network, followed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses.
Gene expression profile data and clinical information
We downloaded 472 SKCM patients’ gene expression profile data using Illumina HiSeq 2000 from the TCGA database (https://gdc-portal.nci.nih.gov/). According to the clinical information, 460 tumor samples had information about metastasis, of which 358 metastatic and 102 nonmetastatic SKCM tumor samples were used as a training dataset. Table 1 showed all the clinical information present in the dataset; there was significant association in pathologic N status, tumor recurrence, and death between metastatic SKCM and nonmetastatic SKCM. In addition, we downloaded another data set of SKCM patient gene expression profiles (GSE65904) from the National Center for Biotechnology Information database (NCBI) (https://www.ncbi.nlm.nih.gov/), which contained 214 SKCM tumor samples. The detection platform used was the Illumina HumanHT-12 V4.0 expression beadchip. According to the clinical information, 150 tumor samples had information about metastasis, of which 83 metastatic and 67 nonmetastatic SKCM tumor samples were used as a validation dataset.
Data preprocessing and DEG screening
Tumor samples in the training dataset were divided into metastatic and nonmetastatic groups according to the clinical information. DEGs in both groups were analyzed using the limma package (https://bioconductor.org/packages/release/bioc/html/limma.html) version 3.34.7 of R3.4.1. FDR < 0.05 and |log2FC|> 1 were used as the cut-off threshold.
Next, we performed two-way hierarchical clustering based on the centered Pearson correlation algorithm based on the DEG expression in the training dataset using pheatmap version 1.0.8 (https://cran.r-project.org/web/packages/pheatmap/index.html) of R3.4.1, and then performed the following analysis:
To determine whether the clinical information about the tumor samples in different clusters was significantly different, we ran the chisq.test function (http://www.bioconductor.org/help/search/index.html?q=chisq.test/) of R3.4.1 on the clinical information.
According to the survival prognosis information about tumor samples in different clusters, we used the Kaplan–Meier (KM) curve in the survival package version 2.41–1 (http://bioconductor.org/packages/survivalr/) of R3.4.1 to calculate the correlation between clusters and survival prognosis.
Screening of independent prognostic DEGs
On the basis of the clinical prognosis information about the SKCM tumor samples included in this analysis, we univariate Cox regression analysis in the survival package version 2.41–1 of R3.4.1  to screen DEGs with significant differences in prognosis. We further used multivariate Cox regression analysis to screen independent prognosis-related DEGs, and log-rank P value < 0.05 was selected as the threshold for significant correlation.
Screening of feature DEGs
On the basis of previously screened DEGs independently associated with survival prognosis, we used the glm function of R3.4.1 to perform LR [14, 15] to screen out feature DEGs and classify metastatic and nonmetastatic SKCM tumor samples. All genes with P < 0.05 were considered feature DEGs. Next, on the basis of the expression of the feature DEGs, we calculated the discriminant accuracy and classified all tumor samples in the training and validation datasets into metastatic and nonmetastatic groups. The KM survival curve in the survival package version 2.41-1 of R3.4.1 was used to perform correlation analysis of the actual survival prognosis of the metastatic and nonmetastatic SKCM tumor samples obtained from the LR classification.
Cluster analysis of feature DEGs
For the training and validation datasets, we used pheatmap version 1.0.8  of R3.4.1 to perform two-way hierarchical clustering based on the centered Pearson correlation algorithm based on feature DEG expression. On the basis of the survival prognosis information about samples in different clusters, we used the KM survival curve in the survival package version 2.41-1 (http://bioconductor.org/packages/survivalr/) of R3.4.1 to determine the correlation between clusters and survival prognosis in the training and validation datasets. The expression levels of feature DEGs in different clusters were then demonstrated.
PPI network construction function analysis
We constructed a PPI network to evaluate the interactions between proteins encoded by the DEGs. The PPI network was based on the STRING database version 10.0 (http://string-db.org/)  and was visualized using Cytoscape software version 3.6.1 (http://www.cytoscape.org/) . Next, we used the Database for Annotation, Visualization and Integrated Discovery (DAVID) version 6.8 (https://david.ncifcrf.gov/) to conduct GO function annotation, including Biology Process, and KEGG signaling pathway enrichment analysis [19, 20]. Finally, significantly relevant GO functions were screened using Fisher’s exact test with P < 0.05.
On the basis of the clinical information, SKCM tumor samples in the TCGA training dataset were divided into metastatic (n = 358) and nonmetastatic (n = 102) groups. We screened 876 DEGs, including 353 significantly downregulated and 523 significantly upregulated DEGs. Figure 1a showed the test volcano map. The two-dimensional hierarchical clustering heat map (Fig. 1b) showed that samples with similar gene expression patterns were hierarchically clustered into the same group, indicating that the screened DEGs can well distinguish tumor samples with different prognoses. The two clusters of the clustering chart contained 247 and 213 tumor samples, respectively. Comparing the clinical information of the tumor samples in the two clusters showed that the samples were significantly different in age, pathologic T status, pathologic stage, radiotherapy, tumor metastasis, and survival time (Table 2). In addition, comparing the prognostic levels of the tumor samples showed that the survival ratio of cluster 1 was significantly better than cluster 2 (P = 9.141e−09) (Fig. 1c).
Screening of independent prognostic DEGs
On the basis of the clinical prognosis information about the 460 SKCM tumor samples, we analyzed the 876 screened DEGs by univariate single-factor Cox regression analysis using the survival package version 2.41-1 of R3.4.1 and screened 435 DEGs significantly associated with prognosis. These 435 DEGs were further analyzed using multivariate Cox regression analysis, and 61 independent prognostic DEGs were screened.
Screening of feature DEGs by LR
The 61 independent prognostic DEGs were further screened using LR to obtain 11 feature DEGs: ZNF750, NLRP6, transglutaminase 3 (TGM3), KRTDAP, CAMSAP3, KRT6C, CALML5, SPRR2E, CD3G, RTP5, and FAM83C (Table 3). LR was also used to classify the tumor samples in the training and validation datasets on the basis of the expression levels of these 11 feature DEGs (Table 4). Metastatic SKCM accounted for 98% of the tumor samples in the training dataset and 87% of the tumor samples in the validation dataset. In addition, KM curve analysis was performed on the basis of LR classification prediction results (Fig. 2a, b). In the TCGA training and GSE65904 validation datasets, the predicted survival prognosis of the nonmetastatic group was better than the metastatic group (P = 1.578e−08 and 2.786e−08, respectively).
Cluster analysis of feature DEGs
In the TCGA training and GSE65904 validation datasets, we performed two-way hierarchical clustering based on the centered Pearson correlation algorithm based on the expression of the 11 feature DEGs (Fig. 3a). We also performed KM curve analysis on the tumor samples in clusters 1 and 2 on the basis of the expression of the 11 DEGs. In TCGA the training and GSE65904 validation datasets, the survival prognosis in cluster 1 was better compared to cluster 2 (P = 7.268e−06 and 1.199e−03, respectively) (Fig. 3b). The expression levels of the 11 feature DEGs in the two clusters of the TCGA training and GSE65904 validation datasets are shown in Fig. 4. Tumor samples with similar gene expression patterns were hierarchically clustered into the same group, indicating that the selected DEGs can well distinguish tumor samples with different prognoses. We identified eight genes (CALML5, CAMSAP3, FAM83C, KRTDAP, SPRR2E, TGM3, and ZNF750) with high expression associated with poor prognosis and three genes (CD3G, NLRP6, and RTP5) with low expression associated with poor prognosis (Fig. 4a). In the GSE65904 validation dataset, we obtained similar results: three genes (CD3G, NLRP6, and RTP5) with low expression were associated with poor prognosis (Fig. 4b).
We constructed a PPI network based on the 876 screened DEGs using the STRING online database. After screening the 11 feature DEGs using LR, we only focused on their PPI network, which contained 107 pairs of interaction connections related to the 11 feature DEGs. Cytoscape software was used to construct the PPI network (Fig. 5) containing 78 nodes (65 downregulated and 13 upregulated) and 107 edges. We collected 12 targets of CALML5, including TGM3, PDE1C, CASP14, and CASP14, and 26 targets of SPRR2E, including SPRR2A, IVL, and LCE2A.
GO and KEGG signaling pathway enrichment analysis
We obtained 21 significantly related GO biological processes and 9 KEGG signaling pathways using the online tool DAVID (Table 5). For GO functions, the DEGs were primarily enriched in “epidermis development,” “keratinization,” “keratinocyte differentiation,” and “establishment of skin barrier,” which involved many genes, such as CALML5, SPRR2A, and HRNR. For KEGG signaling pathways, the DEGs were primarily enriched in “calcium signal pathway,” “GnRH signal pathway,” and “melanin production pathway.” Of the 9 KEGG signaling pathways, in the “hsa04916: Melanogenesis” pathway, CALML3 and CALML5 deserved attention, of which CALML5 was one of important feature DEGs.
SKCM is the most malignant skin tumor type derived from melanocytes, with high mortality, high metastasis, and difficulty in treating. As SKCM grows, cancer cells penetrate the skin and mucous membranes, eventually reach blood vessels or lymphatic channels, and quickly spread throughout the body and major organs . Metastasis is an important cause of poor prognosis in SKCM patients . The 5-year survival rate of metastatic SKCM patients is only 4.6%. Therefore, there is an urgent need to identify potential prognosis-related genes for SKCM and provide more powerful gene markers for the diagnosis of metastatic prognosis in SKCM patients.
In this study, we performed a comprehensive bioinformatics analysis based on the gene expression profile datasets (TCGA and GSE65904 databases) of SKCM patients. By classifying and comparing the gene expression data in the training dataset, we found 876 DEGs, of which 353 were downregulated and 523 were upregulated. With univariate and multivariate Cox regression analysis in the survival package of R3.4.1, we screened 61 independent prognostic DEGs, and then using LR and cluster analysis, we successfully screened 11 feature DEGs: ZNF750, NLRP6, TGM3, KRTDAP, CAMSAP3, KRT6C, CALML5, SPRR2E, CD3G, RTP5, and FAM83C. From among these 11 feature DEGs, we screened feature DEGs that significantly increased or decreased between cluster 1 and cluster 2 in the training dataset with P < 0.001. To understand the biological functions of the screened DEGs, we used the STRING database to search for interactions between the products of screened DEGs and constructed a PPI network. Then, we used the DAVID online tool perform GO function and KEGG signaling pathway enrichment analysis of genes on the PPI network. Of the screened 9 KEGG signaling pathways, one is worth pointing out: “hsa04916: Melanogenesis.” There are two genes involved in this pathway, CALML3 and CALML5, of which CALML5 is one of the important feature DEGs.
CALML5 encodes the CALML5 protein with 146 amino acids. Unlike the generally expressed calmodulin, CALML5 expression is limited to the epidermis and other layered epithelial tissues, and it is highly expressed in the differentiated epidermis . In this study, we found that CALML5 is involved in epidermal development, the calcium signaling pathway, and vascular smooth muscle contraction. Previous studies have reported that CALML5 is involved in terminal differentiation of keratinocytes and encodes a calcium-binding protein expressed in the epidermis . In addition, CALML5, as a calmodulin-like protein, participates in not only epidermal differentiation but also intracellular signal transduction. Therefore, CALML5 plays an important role in cell proliferation, differentiation, apoptosis, and migration . Kurozumi et al. showed that CALML5 is a key gene for lymphatic vascular infiltration in early breast cancer and has potential prognostic value . In addition, CALML5 is involved in the “hsa04916: Melanogenesis” pathway. Mac et al.  showed that calmodulin-like protein can the SKCM cell growth. In addition, Ke et al.  found that CALML5 may be considered a novel biomarker for lung adenocarcinoma diagnosis, which is useful for predicting the risk of lung adenocarcinoma. Misawa et al.  found that CALML5 has a high predictive ability as an emerging biomarker for a validation set, capable of discriminating between the plasma of patients and healthy individuals. However, a few studies have reported an association between CALML5 and SKCM, although the mechanism is unclear. Our results will provide a novel view for the prognosis of SKCM.
In this study, we collected 12 targets of CALML5, including TGM3, PDE1C, and CASP14. Smirnov et al.  showed that TGM3 is absent in melanocytes as well as SKCM samples and that the expression pattern of TGM3 renders it a potential specific marker for basal cell carcinoma diagnosis. Hu et al.  reported that TGM3 controls multiple oncogenic pathways in hepatocellular carcinoma (HCC), contributing to increased cell proliferation and epithelial–mesenchymal transition (EMT). TGM3 also potentially enhances HCC metastasis. TGM3 may serve as a novel therapeutic target in HCC. Wu et al.  found that PDE1C is associated with SKCM development. Shimizu et al. found that PDE1C messenger RNA is expressed and may play an important role in human malignant SKCM melanoma-associated antigen cells. Chen et al.  reported that CASP14 might be a potential biomarker for gastric cancer diagnosis and an independent prognostic factor of gastric cancer. These results indicate that TGM3 may be an important target gene for SKCM treatment. Thus, CALML5 may play an important role in SKCM metastasis and promote SKCM occurrence and development, so CALML5 may become a potential target for the treatment of SKCM metastasis in the future.
Another calmodulin-like protein, CALML3, is a tumor suppressor gene, which significantly inhibits liver cancer growth and lung metastasis. CALML3 is a gene with “early warning” value for liver cancer and lung metastasis, and it is expected to become a new marker for early diagnosis of lung metastasis of liver cancer and a new target for inhibiting liver cancer growth and lung metastasis . These findings show that CALML3 may be closely related to SKCM metastasis and prognosis, and that it has certain clinical significance for the prediction of SKCM and its prognosis.
Although we identified 11 feature DEGs related to SKCM metastasis, the detailed mechanisms have not yet been explored. For example, a further accurate classification with a large sample size and clinical information needs to identify SKCM metastatic and nonmetastatic patients. In addition, whether the 11 feature DEGs are involved in several molecular pathways, such as hsa04916: Melanogenesis, needs to be investigated. Corresponding experimental research is also needed to verify feature gene functions.
We identified 876 DEGs (353 downregulated and 523 upregulated) in the TCGA training dataset. In addition, 11 important prognostic-related feature DEGs, such as ZNF750, NLRP6, TGM3, CALML5, CD3G, and RTP5, may play an important role in SKCM metastasis. These DEGs are involved in 9 KEGG signaling pathways, such as the “hsa04916: Melanogenesis” pathway. Of the 11 feature DEGs, CALML5 is involved in the “hsa04916: Melanogenesis” pathway, 12 targets of which were collected, such as TGM3, PDE1C, and CASP14. This study provides new ideas for exploring the molecular mechanism underlying SKCM metastasis and finding new diagnostic prognostic markers.
Availability of data and materials
The datasets supporting the conclusions of this article are available in the [GSE65904] and [TCGA] repository, [https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE65904] and [https://xenabrowser.net/datapages/?cohort=GDC%20TCGA%20Melanoma%20(SKCM)&removeHub=https%3A%2F%2Fxena.treehouse.gi.ucsc.edu%3A443].
Differentially expressed genes
Kyoto Encyclopedia of Genes and Genomes
Skin cutaneous melanoma
The Cancer Genome Atlas
Shah DJ, Dronca RS. Latest advances in chemotherapeutic, targeted, and immune approaches in the treatment of metastatic melanoma. Mayo Clin Proc. 2014;89(4):504–19.
Leonardi GC, Falzone L, Salemi R, Zanghì A, Spandidos DA, McCubrey JA, et al. Cutaneous melanoma: from pathogenesis to therapy (Review). Int J Oncol. 2018;52(4):1071–80.
Xiong TF, Pan FQ, Li D. Expression and clinical significance of S100 family genes in patients with melanoma. Melanoma Res. 2019;29(1):23–9.
Wagner NB, Forschner A, Leiter U, Garbe C, Eigentler TK. S100B and LDH as early prognostic markers for response and overall survival in melanoma patients treated with anti-PD-1 or combined anti-PD-1 plus anti-CTLA-4 antibodies. Br J Cancer. 2018;119(3):339–46.
Da Forno PD, Pringle JH, Hutchinson P, Osborn J, Huang Q, Potter L, et al. WNT5A expression increases during melanoma progression and correlates with outcome. Clin Cancer Res. 2008;14(18):5825–32.
Bosserhoff AK. Melanoma inhibitory activity (MIA): an important molecule in melanoma development and progression. Pigment Cell Res. 2005;18(6):411–6.
Ci C, Tang B, Lyu D, Liu W, Qiang D, Ji X, et al. Overexpression of CDCA8 promotes the malignant progression of cutaneous melanoma and leads to poor prognosis. Int J Mol Med. 2019;43(1):404–12.
Yang L, Xu Y, Yan Y, Luo P, Chen S, Zheng B, et al. Common nevus and skin cutaneous melanoma: prognostic genes identified by gene co-expression network analysis. Genes (Basel). 2019;10(10):747.
Bisaso KR, Karungi SA, Kiragga A, Mukonzo JK, Castelnuovo B. A comparative study of logistic regression based machine learning techniques for prediction of early virological suppression in antiretroviral initiating HIV patients. BMC Med Inform Decis Mak. 2018;18(1):77.
Zejnullahu VA, Ukella-Lleshi D, Zejnullahu VA, Miftari E, Govori V. Prevalence of postpartum depression at the clinic for obstetrics and gynecology in Kosovo teaching hospital: Demographic, obstetric and psychosocial risk factors. Eur J Obstet Gynecol Reprod Biol. 2021;256:215–20.
Sufriyana H, Husnayain A, Chen YL, Kuo CY, Singh O, Yeh TY, et al. Comparison of multivariable logistic regression and other machine learning algorithms for prognostic prediction studies in pregnancy care: systematic review and meta-analysis. JMIR Med Inform. 2020;8(11):e16503.
Sokou R, Piovani D, Konstantinidi A, Tsantes AG, Parastatidou S, Lampridou M, et al. A risk score for predicting the incidence of hemorrhage in critically ill neonates: development and validation study. Thromb Haemost. 2021;121:131–9.
Wang P, Wang Y, Hang B, Zou X, Mao JH. A novel gene expression-based prognostic scoring system to predict survival in gastric cancer. Oncotarget. 2016;7(34):55343–51.
Lee JW, Lee JB, Park M, Song SH. An extensive comparison of recent classification tools applied to microarray data. Comput Stat Data Anal. 2005;48(4):869–85.
Zhu J, Hastie T. Classification of gene microarrays by penalized logistic regression. Biostatistics (Oxford, England). 2004;5(3):427–43.
Wang L, Cao C, Ma Q, Zeng Q, Wang H, Cheng Z, et al. RNA-seq analyses of multiple meristems of soybean: novel and alternative transcripts, evolutionary and functional implications. BMC Plant Biol. 2014;14:169.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.
da Huang W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37(1):1–13.
Zhang X, Ding C, Tian H, Dong X, Meng X, Zhu W, et al. ZNF23 suppresses cutaneous melanoma cell malignancy via mitochondria-dependent pathway. Cellular Physiol Biochem. 2017;43(1):147–57.
Cohen I, Birnbaum RY, Leibson K, Taube R, Sivan S, Birk OS. ZNF750 is expressed in differentiated keratinocytes and regulates epidermal late differentiation genes. PLoS ONE. 2012;7(8):e42628.
Babini E, Bertini I, Capozzi F, Chirivino E, Luchinat C. A structural and dynamic characterization of the EF-hand protein CLSP. Structure. 2006;14(6):1029–38.
Donovan M, Ambach A, Thomas-Collignon A, Prado C, Bernard D, Jammayrac O, et al. Calmodulin-like skin protein level increases in the differentiated epidermal layers in atopic dermatitis. Exp Dermatol. 2013;22(12):836–7.
Villalobo A, González-Muñoz M, Berchtold MW. Proteins with calmodulin-like domains: structures and functional roles. Cellular and molecular life sciences : CMLS. 2019;76(12):2299–328.
Kurozumi S, Joseph C, Sonbul S, Alsaeed S, Kariri Y, Aljohani A, et al. A key genomic subtype associated with lymphovascular invasion in invasive breast cancer. Br J Cancer. 2019;120(12):1129–36.
Mac Neil S, Walker SW, Senior HJ, Bleehen SS, Tomlinson S. Effects of extracellular calmodulin and calmodulin antagonists on B16 melanoma cell growth. J Invest Dermatol. 1984;83(1):15–9.
Ke H, Wu Y, Wang R, Wu X. Creation of a prognostic risk prediction model for lung adenocarcinoma based on gene expression, methylation, and clinical characteristics. Med Sci Monit. 2020;26:e925833.
Misawa K, Imai A, Matsui H, Kanai A, Misawa Y, Mochizuki D, et al. Identification of novel methylation markers in HPV-associated oropharyngeal cancer: genome-wide discovery, tissue verification and validation testing in ctDNA. Oncogene. 2020;39(24):4741–55.
Smirnov A, Anemona L, Montanaro M, Mauriello A, Annicchiarico-Petruzzelli M, Campione E, et al. Transglutaminase 3 is expressed in basal cell carcinoma of the skin. Eur J Dermatol. 2019;29(5):477–83.
Hu JW, Yang ZF, Li J, Hu B, Luo CB, Zhu K, et al. TGM3 promotes epithelial-mesenchymal transition and hepatocellular carcinogenesis and predicts poor prognosis for patients after curative resection. Dig Liver Dis. 2020;52(6):668–76.
Wu W, Liu H, Song F, Chen LS, Kraft P, Wei Q, et al. Associations between smoking behavior-related alleles and the risk of melanoma. Oncotarget. 2016;7(30):47366–75.
Chen Q, Huang X, Dong X, Wu J, Teng F, Xu H. Long non-coding RNA ERICH3-AS1 is an unfavorable prognostic factor for gastric cancer. PeerJ. 2020;8:e8050.
Yang B, Li M, Tang W, Liu W, Zhang S, Chen L, et al. Dynamic network biomarker indicates pulmonary metastasis at the tipping point of hepatocellular carcinoma. Nat Commun. 2018;9(1):678.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Jia, G., Song, Z., Xu, Z. et al. Screening of gene markers related to the prognosis of metastatic skin cutaneous melanoma based on Logit regression and survival analysis. BMC Med Genomics 14, 96 (2021). https://doi.org/10.1186/s12920-021-00923-0
- Cutaneous melanoma
- Differentially expressed genes