The Overexpression of SLC25A13 Predicts Poor Prognosis and Is Correlated with Immune Cell Infiltration in Patients with Skin Cutaneous Melanoma

lv, Yue; Yuan, Chun-hui; Han, Lu-yao; Huang, Gao-ru; Ju, Ling-ce; Chen, Ling-hui; Han, Hai-ying; Zhang, Chong; Zeng, Ling-hui

doi:https://doi.org/10.1155/2022/4091978

Disease Markers

On this page

Abstract Introduction Methods Results Discussion Conclusions Data Availability Conflicts of Interest Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 4091978 | https://doi.org/10.1155/2022/4091978

The Overexpression of SLC25A13 Predicts Poor Prognosis and Is Correlated with Immune Cell Infiltration in Patients with Skin Cutaneous Melanoma

Yue lv,¹Chun-hui Yuan,²Lu-yao Han,¹Gao-ru Huang,¹Ling-ce Ju,¹Ling-hui Chen,³Hai-ying Han,¹Chong Zhang,^1,2and Ling-hui Zeng^1,2

Academic Editor: Liu Jinhui

Received14 Feb 2022

Revised07 Apr 2022

Accepted09 Apr 2022

Published14 May 2022

Abstract

Purpose. Skin cutaneous melanoma (SKCM) is one of the most malignant and aggressive cancers with poor prognosis due to its rapid progression towards metastasis. Thus, finding clinically relevant biomarkers for early diagnosis, prognosis, and therapy prediction is essential. This study focused on the identification of SLC25A13 as a novel biomarker for SKCM and is aimed at investigating the biological functions of solute carrier family 25 member 13 (SLC25A13) in the development of SKCM. Methods. GEPIA was used to analyze the diagnostic and prognostic values of SLC25A13 in SKCM using the TCGA dataset. PrognoScan was used to validate the prognostic value of SLC25A13 and its coexpressed genes in SKCM. TISIDB was established to reveal the relationship between the expression of SLC25A13 and immune infiltration in SKCM. The protein expression of SLC25A13 in SKCM was evaluated by the Human Protein Atlas. The signaling pathways and biological functions of SLC25A13 in SKCM were analyzed by LinkOmics. Metascape was applied to analyze the functional enrichment analysis of SLC25A13. Protein-protein interaction analysis of SLC25A13 was performed by GeneMANIA. Results. The mRNA and protein levels of SLC25A13 in the SKCM were much higher than those in the normal tissue. Furthermore, the overexpression of SLC25A13 predicts worse outcomes of SKCM patients. Moreover, the SLC25A13 expression was negatively correlated with the immune infiltration level of SKCM. The overexpression of SLC25A13 coexpressed genes, such as ACLY and AFG3L2, and SCL25A13 interacting genes also predicted the unfavorable prognosis of SKCM patients. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of SLC25A13 coexpressed genes showed that these genes are enriched in ATPase activity, cell cycle, mTOR, and VEGFA-VEGFR2 signaling pathways, which were relevant to tumor development and angiogenesis. Gene set enrichment analysis (GSEA) demonstrated that the SLC25A13 expression was related to infiltrating immune cells in SKCM. Conclusion. Our findings revealed that SLC25A13 might be a potential prognostic and therapeutic biomarker for SKCM.

1. Introduction

The incidence of skin cutaneous melanoma (SKCM) is increasing in both males and females, and SKCM caused approximate 72% of deaths in skin carcinoma [1]. SKCM lacks effective treatment except early surgical resection, and it usually results in a poor prognosis and shockingly high mortality [2]. Meanwhile, SKCM is highly heterogeneous, which makes personalized treatment difficult [3]. Immune checkpoint inhibitors, such as cytotoxic T-lymphocyte antigen-4 (CTLA-4) and programmed cell death protein 1 (PD-1) inhibitors, are the main treatment measures for advanced-stage SKCM, but the efficacy of immunotherapy is limited due to different response rates [4]. In addition, primary and secondary resistance and the absence of predictive markers of response are biological challenges and clinical dilemmas in immunotherapy for SKCM patients [5]. Thus, it is urgent to identify a novel predictive biomarker for the immunotherapy of SKCM patients.

Solute carrier family 25 member 13 (SLC25A13), a calcium-binding/stimulated aspartate–glutamate carrier, plays a key role in metabolic pathways including aerobic glycolysis, gluconeogenesis, urea cycle, and synthesis of proteins and nucleotides [6]. The genetic deficiency of SLC25A13 causes adult-onset type II citrullinemia [7]. Furthermore, pathogenic germline variants correlated with increased tumor mutational burden have been found in 7 genes, including SLC25A13 and TP53, and germline variants can predict the response to immune checkpoint inhibitors [8]. In addition, SLC25A13 is correlated with tumor aggressiveness and poorer prognosis of colorectal cancer [9]. Thus, it is of importance to reveal the role of SLC25A13 in cancer development and evaluate the prognostic value of immunotherapy of SKCM patients.

This study compared the gene transcription levels between SKCM and Paracancerous tissue using RNA sequencing (data not shown), and the top 100 genes overexpressed in SKCM were considered as potential oncogenes involved in cancer process and verified by the cancer genome atlas (TCGA) data sets of SKCM. In this study, we investigated the diagnostic and prognostic value of SLC25A13 for SKCM patients and clarified the correlation between the expression of SLC25A13 and immune cell infiltration in SKCM patients; meanwhile, we also predicted the potential gene function of SLC25A13 by enrichment analysis of its neighbor genes in SKCM (Figure 1).

2. Methods

2.1. Analysis of the Expression of SLC25A13 in SKCM and Normal Skin Tissues

Gene Expression Profiling Interactive Analysis (GEPIA, http://gepia.cancer-pku.cn/) is a web-based tool, and The Cancer Genome Atlas (TCGA) and Genotype Tissue Expression (GTEx) database provide fast and customizable functionalities [10]. In this study, the expression of SLC25A13 in SKCM and normal skin tissues was compared and represented by box plots via the GEPIA. In addition, the overall survival and disease-free survival of SKCM patients based on the expression of SLC25A13 were analyzed using GEPIA.

Human Protein Atlas (https://www.proteinatlas.org) provides the expression profiles in human tissues of genes both on the mRNA and protein level [11]. The protein levels of SLC25A13 were compared in SKCM and normal skin tissues by Human Protein Atlas.

2.2. Survival Analyses in SKCM Patients

PROGgene (http://www.compbio.iupui.edu/proggene) serves as a tool to perform survival analysis on single genes, multiple genes as a signature, ratio of expression of two genes, and curated/published gene signatures [12]. The overall survival of SKCM patients based on the expression of SLC25A13 was analyzed using PROGgene.

PrognoScan (http://dna00.bio.kyutech.ac.jp/PrognoScan/) shows the prognostic values of genes in multiple microarray datasets with meta-analysis [13]. PrognoScan was used to explore the prognostic value of interacting genes of SLC25A13 in SKCM patients.

2.3. Clinical Prediction Value of the Established Prognostic Model

The TCGA database (https://portal.gdc.cancer.gov/) and the UCSC Xena database (https://xena.ucsc.edu/) were used to obtain SKCM data. All subsequent analyses were performed via version 4.1.3. To test whether the expression of SLC25A13 was the effective prognostic indicator, the univariate and multivariate Cox regressions were performed [14, 15]. Then, diagnostic ROC and PR analysis of the SLC25A13 expression in normal samples and SKCM samples was performed by package “ROCR” [16]. Finally, to assess the risk of death for SKCM, a nomogram was constructed by combining the expression of SLC25A13 and clinical data. The patient “TCGA-EE-A2A2-06A” was illustrated.

2.4. Correlation Analysis of Immune Mutation

The mutation analysis of SKCM was used the package “maftools” [17, 18]. In the analysis process, a sample “TCGA-FW-A3R5-06A-11D-A23B-08 31097” was deleted because its mutation count is an outlier.

2.5. Correlation Analysis between the Expression of SLC25A13 and Immune Features in SKCM Patients

TISIDB (http://cis.hku.hk/TISIDB/index.php) is a website portal for interaction between tumor and immune systems in 28 types of tumor-infiltrating lymphocytes across human cancers [19]. The correlation between the abundance of tumor-infiltrating lymphocytes and the expression of SLC25A13 was investigated by TISIDB, and the relationship between abundance of immunomodulators and the expression of SLC25A13 was determined using TISIDB.

2.6. Analysis of Coexpressed Genes of SLC25A13 in SKCM

The multiomics data from all 32 TCGA cancer types and 10 Clinical Proteomics Tumor Analysis Consortium (CPTAC) cancer cohorts are available at LinkedOmics (http://www.linkedomics.org) [20]. In this study, LinkFinder module was used to investigate the coexpression genes of SLC25A13, and the LinkInterpreter module was employed to perform enrichment analysis of co-expression genes associated with SLC25A13 in SKCM.

The STRING database (http://string-db.org/) aims to predict protein-protein associations including direct (physical) and indirect (functional) relations from knowledge transfer between organisms and interactions aggregated from other databases [21]. Gene Ontology (GO) and KEGG analysis of SLC25A13 coexpressed genes were evaluated by the STRING database to predict the biological functional of SLC25A13.

Cytoscape is an open source software project for integrating biomolecular interaction networks with high-throughput expression data and other molecular states into a unified conceptual framework [22]. Cytoscape and MCODE were established to estimate the most indispensable module in the protein-protein interaction networks.

2.7. Functional Enrichment Analysis and Gene Set Enrichment Analysis (GSEA) of SLC25A13

GeneMANIA (http://www.genemania.org), a flexible web interface, is utilized for generating hypotheses about gene function, analyzing gene lists, and prioritizing genes for functional assays [23]. In the present study, the protein-protein interaction network of SLC25A13 and its interacting genes was conducted by GeneMANIA database.

Metascape (https://metascape.org/) is a web site that combines functional enrichment, interactome analysis, gene annotation, and membership search to leverage more than 40 independent knowledge bases within one integrated portal [24]. In this study, Metascape was used to demonstrate the functional enrichment analysis of SLC25A13 and its neighboring genes.

GSEA is a powerful analytical method to perform genome-wide expression profiling of two phenotypes and evaluate cumulative changes in the expression of groups of multiple genes defined based on previous biological knowledge [25]. In this study, GSEA was applied to evaluate the biological functions and signaling pathways regulated by SLC25A13.

3. Results

3.1. SLC25A13 Is Overexpressed in SKCM Compared with Normal Skin Tissues

Firstly, the levels of SLC25A13 in SKCM and normal skin tissues were analyzed by immunohistochemistry data collected from Human Protein Atlas [11]. As shown in Figure 2(a), high levels of SLC25A13 were detected in SKCM samples, whereas the expression of SLC25A13 could not or barely be detected in normal skin tissues. Statistical analysis manifested that the expression of SLC25A13 in the SKCM was much higher than that in normal skin tissues (Figure 2(b)). In line with the protein expression, the mRNA levels of SLC25A13 were also upregulated in SKCM compared with normal skin tissues (Figure 2(c)) [10].

(a)

(b)

(c)

3.2. High Levels of SLC25A13 Predict Poor Disease Outcome in SKCM

To find whether the SLC25A13 expression is associated with overall survival (OS) and disease-free survival (DFS), we compared the OS and DFS of SKCM patients with high or low SLC25A13 expression. As shown in Figures 3(a) and 3(b), patients with low levels of SLC25A13 had longer OS () and DFS () [10]. Furthermore, this was also validated using PROGgene based on GSE19234 and GSE22153 datasets. As shown in Figures 3(c) and 3(d), the overexpression of SLC25A13 was associated with worse OS in SKCM patients [13, 26, 27]. To explore the prognostic value of SLC25A13 gene in SKCM, univariate COX analysis was performed, and SLC25A13 was found to be an independent prognostic factor (Figure 3(e), ). SLC25A13 was also found to be an independent prognostic factor by multivariate COX analysis (Figure 3(f), ). To explore the diagnostic value of SLC25A13 gene in SKCM, diagnostic ROC analysis of SLC25A13 gene in 556 normal subjects and 468 SKCM patients was performed. ROC analysis showed that the area under the curve (AUC) was 0.849 (Figure 3(g)). Meanwhile, the area under the curve (AUC) was 0.894 according to PRC analysis. These results indicated that SLC25A13 might be a good diagnostic molecule for SKCM. To further determine the survival of SKCM patients, we drew a nomogram combining the risk value and clinical characteristics of the model [28]. As shown in Figure 4, we found that the 1-, 3-, and 5-year mortality rates of patient “TCGA-EE-A2A2-06” were 0.102, 0.478 and 0.64, respectively.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

Figure 3

The overexpression of SLC25A13 predicts poor disease outcome in SKCM patients. (a, b) High levels of SLC25A13 predicted poor disease outcome in SKCM, and the data were obtained by GEPIA. (c, d) Overexpressed SLC25A13 predicted poor disease outcome in SKCM, and the data were obtained by PROGgene. (e) Univariate COX regression showed that SLC25A13 was an independent prognostic factor for SKCM (). (f) Multivariate COX regression showed that SLC25A13 was an independent prognostic factor for SKCM (). (g) ROC analysis showed that the area under the curve (AUC) was 0.849. Meanwhile, the area under the curve (AUC) was 0.894 according to PRC analysis. These results suggest that this gene may be a good diagnostic molecule for SKCM.

3.3. The SLC25A13 Expression Is Associated with Immune Infiltration of Tumor Immune Cells and Immunomodulators in SKCM

The TISIDB database was used to evaluate the relation between the expression of SLC25A13 and immune infiltration of tumors [19]. As shown in Figure 5(a), the level of SLC25A13 was negatively correlated with the level of 28 immune cell types in numerous cancer types. Specifically, significant negative correlations between the expression of SLC25A13 and the immune infiltration levels of multiple immune cell types in SKCM and the proportions of Th17 (, ), Th1 (, ), Tem CD8 (, ), and MDSC (, ) were lower in the high SLC25A13 expression group than in the low SLC25A13 expression group (Figure 5(b)). Besides, the SLC25A13 expression was also negatively correlated with multiple immune stimulators in SKCM (Figure 5(c)). As shown in Figure 5(d), the top-ranked genes include the following: TNFRSF25 (, ), ICOSLG (, ), TNFRSF18 (, ), and C10orf54 (, ). For the MHC molecules (Figure 5(e)), SLC25A13 was mostly negatively correlated with HLA-DMA (, ), HLA-DMB (, ), HLA-DRA (, ), and HLA-DOA (, ) (Figure 5(f)). These results indicated that SLC25A13 was indeed a determinant of immune infiltration and immunomodulators in SKCM.

(a)

(b)

(c)

(d)

(e)

(f)

3.4. Analysis of Coexpressed Genes Correlated with SLC25A13 in SKCM

To investigate the SLC25A13 coexpressed genes, we identified the potential candidate genes using LinkedOmics [20] and measured their Pearson coefficient of correlation with SLC25A13. The top 50 significant genes that were negatively (Figure 6(a)) and positively (Figure 6(b)) correlated with SLC25A13 were presented in the heat map. As shown in Figure 6(c), 8959 genes (with red label) showed statistically significant positive correlations with SLC25A13, while 10769 genes (with green label) showed statistically significant negative correlation with SLC25A13.

(a)

(b)

(c)

(d)

3.5. The Prognostic Value of SLC25A13 Associated Genes for SKCM Patients

A total of 295 coexpressed genes with strong positive correlations with SLC25A13 were used to construct a protein-protein interaction network, and the vital module in the network was highlighted in yellow (Figure 6(d)). Furthermore, the prognostic values of SLC25A13 coexpressed genes in SKCM were investigated using PrognoScan [13]. As shown in Figures 7(a)–7(i), high level of SLC25A13 coexpressed genes predicted poor clinical outcome of SKCM patients (ACLY: [1.07-6.97], ; AFG3L2: [1.16-8.73], ; ASL: [1.10-7.46], ; CHCHD3: [2.60-47.70], ; CTH: [1.41-7.19], ; GFM1: [2.11-29.49], ; MDH2: [1.52-9.46], ; NDUFS1: [1.06-11.22], ; SLC25A15: [1.71-14.69], ). Thus, patients with higher level of SLC25A13 and its core coexpressed genes in SKCM have poorer prognosis relative to those SKCM patients with normal SLC25A13.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

To better understand the biological functions of SLC25A13, the protein-protein interaction network of SLC25A13 and its interacting genes was constructed using GeneMANIA (Figure 8(a)) [23]. Meanwhile, the prognostic values of the interacting genes of SLC25A13 were analyzed using PrognoScan [13]. As shown in Figures 8(b)–8(h), the overexpression of SLC25A13 interacting genes predicted poor clinical outcome of SKCM patients (SLC25A13: [1.64-10.05], ; GOT1: [1.16-7.03], ; GOT2: [1.14-12.13], ; RNF31 [1.07-9.79], ; SIRT7: [1.05-7.98], ; TUBG1: [1.77-17.11], ; ELAVL1: [1.40-10.99], ). Altogether, the overexpression of SCL25A13 interacting genes tended to cause poor prognosis in SKCM patients.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

3.6. GO and KEGG Analysis of SLC25A13 Coexpressed Genes

We input 649 coexpressed genes of SLC25A13 into the STRING database to perform functional enrichment analysis [29]. GO analysis of biological processes revealed that these genes were enriched in mitotic spindle midzone assembly, spindle midzone assembly, DNA replication checkpoint, and pigment biosynthetic process (Figure 9(a)). These genes were located in melanosome membrane, melanosome, and mitochondrial nucleoids (Figure 9(c)) and associated with the molecular functions of retromer complex binding, DNA replication origin binding, ATPase activity, etc. (Figure 9(b)). The results from KEGG analysis indicated that pathways including cell cycle and mTOR signaling pathway were particularly enriched, which were closely related to cancer progression (Figure 9(d)).

(a)

(b)

(c)

(d)

3.7. Functional Enrichment Analysis of SLC25A13 in Patients with SKCM

Metascape analysis showed that SLC25A13 and its neighboring genes were enriched in glucose metabolism, amino acid metabolism, and cell cycle (Supplementary Figure 1A-1B) [24]. Moreover, they were also enriched in pathways such as VEGFA-VEGFR2 signaling pathway and urea cycle and associated pathways, which were relevant to tumor development and angiogenesis. Meanwhile, SLC25A13 and its neighboring genes were enriched in transcription factor target genes, including MEF2C, NCOA2, SNAI1, and SOX10 target genes, which were participated in the cancer progression (Supplementary Figure 1 C-D).

GSEA of SLC25A13 was utilized to further investigate the biological function of SLC25A13 in SKCM by LinkedOmics [20]. GO analysis of molecular function showed that SLC25A13 positively participated in ATPase activity, while negatively regulated antigen binding and immunoglobulin binding (Supplementary Figure 2). GO analysis of the biological process showed that SLC25A13 positively participated in pigment metabolic and negatively regulated T cell activation, lymphocyte mediated immunity, humoral immune response, response to tumor necrosis factor, and leukocyte proliferation (Supplementary Figure 3). GO analysis of cell components displayed that SLC25A13 positively regulated pigment granules but negatively regulated immunological synapse and MHC protein complexes in SKCM (Supplementary Figure 4). KEGG pathway analysis showed that SLC25A13 positively regulated Fanconi anemia pathway, cell cycle, DNA replication, basal transcription factors and negatively regulated natural killer cell-mediated cytotoxicity and antigen processing and presentation (Supplementary Figure 5).

4. Discussion

SKCM is a highly aggressive skin tumor worldwide due its metastatic nature, and the incidence rate of SKCM keeps increasing in recent years [30]. Once SKCM spreads through the dermis, it easily spreads to vital organs and lymph nodes and causes poor prognosis [31]. Therefore, it is urgent to identify biomarkers to predict the prognosis of SKCM. This study indicated that the mRNA level and protein expression of SLC25A13 was significantly higher in SKCM than that in normal skin tissues. Moreover, Cox regression analysis showed that the SLC25A13 overexpression was an independent risk factor for poor outcome of SKCM patients. Furthermore, the overexpression of SCL25A13 associated genes including ACLY, AFG3L2, ASL, CHCHD3, CTH, GFM1, MDH2, NDUFS1, SLC25A15, GOT1, GOT2, RNF31, SIRT7, TUBG1, and ELAVL1 were also related to poor prognosis of SKCM patients. These findings revealed that SLC25A13 is a promising and prognostic biomarker in SKCM.

Multiple members of the mitochondrial carrier family (SLC25) are involved in the invasion and metastasis of many cancers. For example, Wong et al. reported that SLC25A22 promotes the proliferation and migration of colorectal cancer cells, and high levels of SLC25A22 predict poor prognosis in colorectal cancer [32]. SLC25A13 was also found to be pivotal in colorectal cancer aggressiveness [9]. Although the function of SCL25A13 in SKCM development may need further investigation, the function enrichment analysis revealed that SLC25A13 and its neighboring genes were enriched in cell cycle, VEGFA-VEGFR2 signaling, and mTOR pathway, which were related to the development of SKCM. SKCM develops from mutation of pigment-producing cells in the skin, and there is uncontrolled growth of pigment-producing cells during SKCM development [33]. In this study, the function enrichment analysis of SLC25A13 and its neighboring genes revealed that SLC25A13 was involved in pigment granules, pigment biosynthetic, and metabolic processes. Our results suggested that SLC25A13 might be critically involved in the progression of SKCM. It is well known that cancer cells have a higher metabolic rate than normal cells as a result of environmental selection and the accumulation of mutations [34]. Then, we downloaded the mutation data of SKCM through mafTools package and showed the landscape map of the top 10 genes with the highest mutation frequency, among which TTN and MUC16 genes are more prone to mutation (Supplementary Figure 6A). However, there is no relationship between SLC25A13 and tumor mutation load (TMB) through correlation analysis (Supplementary Figure 6B, , ). Thus, SLC25A13 was not associated with TMB of SKCM.

SLC25A13 plays a pivotal role under glucose-deprived conditions and is correlated with tumor aggressiveness [9]. The presence and amount of glucose have a great impact on the immune cell functions of the immune system, resulting in the tumor pathological stage [35]. Thus, we were interested in investigating the relationship between the expression of SLC25A13 and immune cell infiltration of SKCM. The data of the current study revealed that the SLC25A13 expression was negatively correlated with immune cell infiltration, including effector memory CD8+ and CD4+ T cells, Th1 cells, Th2 cells, and Treg cells. Furthermore, the function enrichment analysis revealed that SLC25A13 negatively regulates T cell activation, lymphocyte mediated immunity, humoral immune response, natural killer cell-mediated cytotoxicity, etc. These data indicated that the SLC25A13 expression was negatively correlated with immune cell infiltration in SKCM. Furthermore, SKCM treatment has been revolutionized with the approval of immune checkpoint inhibitors, which have a significant impact on the prognosis of patients with SKCM [36]. Nevertheless, only a small percentage of SKCM patients can benefit from immunotherapy [37]. Tumor immune cell infiltration can affect the sensitivity of immunotherapy, and higher immune cell infiltration scores showed a significant immune therapeutic advantage and clinical benefit [38]. Thus, the expression of SLC25A13 might be a potential biomarker useful for screening suitable SKCM patients for immunotherapy. Comprehensive psychological intervention can remarkably improve the negative emotions, social support degree, sleep quality, and immune functions of cancer patients receiving conventional nursing intervention during the perioperative period [39]. As SLC25A13 could be a biomarker for tumor immune cell infiltration in SKCM, the feasibility of exploring SLC25A13 as a biomarker of psychological intervention in SKCM treatment might need further investigation.

5. Conclusions

This work demonstrated that the SLC25A13 expression was increased in SKCM compared with normal skin tissues, and the high expression of SLC25A13 and its neighbor genes were associated with a poor prognosis in SKCM patients. In addition, the expression of SLC25A13 was negatively correlated with the tumor-infiltrating cells in the tumor immune microenvironment. Functional enrichment analysis of SLC25A13 showed that it was associated with tumor development and progression. These findings suggest that SLC25A13 may be a potential prognostic biomarker for the diagnosis and prognosis of SKCM patients.

Data Availability

All data generated or analyzed used to support the findings of this study are included within the article and supplementary materials.

Conflicts of Interest

The authors report no conflicts of interest in this work.

Acknowledgments

The study was funded by the National Natural Science Foundation of China (no. 32100521 to CHY), Natural Science Foundation of Zhejiang Province (no. LTY21H160001 to LHZ, LQ22C060001 to CHY, and LQ20H160023 to LHC), Huadong Medicine Joint Funds of the Zhejiang Provincial Natural Science Foundation of China (no. LHDMY22H160001 to CZ), Public-Service Technology Research Plan of Zhejiang Province (no. LGF21H310002 to CZ), and Xin-Miao Talent Plan of Zhejiang Province (2022-10 to YL).

Supplementary Materials

Supplementary Figure 1: functional enrichment analysis of SLC25A13 and its neighboring genes in patients with SKCM. (a) Heatmap of GO enriched terms stained by values. (b) GO enrichment network painted by cluster ID, where nodes that share the same cluster ID are typically mutural close. (c, d) Summary of enrichment analysis in transcription factor targets. Supplementary Figure 2: GO analysis of molecular function of SLC25A13 in SKCM using LinkedOmics. (a) GSEA of SLC25A13 was utilized to further investigate the molecular function of SLC25A13 in SKCM. The correlations between SLC25A13 expression and the presence of antigen binding (b), ATPase activity (c), and immunoglobulin binding (d) were analyzed by GSEA using LinkedOmics. Supplementary Figure 3: GO analysis of biological process of SLC25A13 in SKCM using LinkedOmics. (a) GSEA of SLC25A13 was utilized to make further investigation on the biological process of SLC25A13 in SKCM. The pertinence between SLC25A13 expression and the presence of pigment metabolic (b), T cell activation (c), lymphocyte mediated immunity (d), humoral immune response (e), response to tumor necrosis factor (f), and leukocyte proliferation (g) were analyzed by GSEA using LinkedOmics. Supplementary Figure 4: GO analysis of cell component of SLC25A13 in SKCM using LinkedOmics. (a) GSEA of SLC25A13 was utilized to further investigate the cell component of SLC25A13 in SKCM. The correlations between the SLC25A13 expression and the presence of pigment granule (b), immunological synapse (c), and MHC protein complex (d) were analyzed by GSEA using LinkedOmics. Supplementary Figure 5: KEGG pathways analysis of SLC25A13 in SKCM using LinkedOmics. (a) GSEA of SLC25A13 was utilized to further investigate KEGG pathways of SLC25A13 in SKCM. The correlations between SLC25A13 expression and the presence of fanconi anemia pathway (b), cell cycle (c), basal transcription factors (d), natural killer cell mediated cytotoxicity (e), DNA replication (f), and antigen processing (g) were analyzed by GSEA using LinkedOmics. Supplementary Figure 6: mutation analysis of SKCM. (a) The landscape map of the top 10 genes with the highest mutation frequency, among which TTN and MUC16 genes are more prone to mutation. (b) SLC25A13 was not associated with tumor mutation load (TMB). (, ). (Supplementary Materials)

References

A. Jemal, M. Saraiya, P. Patel et al., “Recent trends in cutaneous melanoma incidence and death rates in the United States, 1992-2006,” Journal of the American Academy of Dermatology, vol. 65, no. 5, pp. S17.e1–S17.e11, 2011.
View at: Publisher Site | Google Scholar
A. Ju, J. Tang, S. Chen, Y. Fu, and Y. Luo, “Pyroptosis-related gene signatures can robustly diagnose skin cutaneous melanoma and predict the prognosis,” Frontiers in oncology, vol. 11, p. 709077, 2021.
View at: Publisher Site | Google Scholar
S. Zhou, Y. Han, J. Li et al., “New prognostic biomarkers and drug targets for skin cutaneous melanoma via comprehensive Bioinformatic analysis and validation,” Frontiers in oncology, vol. 11, p. 745384, 2021.
View at: Publisher Site | Google Scholar
S. Zhang, K. Chen, H. Liu et al., “PMEL as a prognostic biomarker and negatively associated with immune infiltration in skin cutaneous melanoma (SKCM),” Journal of Immunotherapy, vol. 44, no. 6, pp. 214–223, 2021.
View at: Publisher Site | Google Scholar
G. C. Leonardi, S. Candido, L. Falzone, D. A. Spandidos, and M. Libra, “Cutaneous melanoma and the immunotherapy revolution (review),” International journal of oncology, vol. 57, no. 3, pp. 609–618, 2020.
View at: Publisher Site | Google Scholar
Z. H. Zhang, W. X. Lin, M. Deng, X. J. Zhao, and Y. Z. Song, “Molecular analysis of _SLC25A13_ gene in human peripheral blood lymphocytes: marked transcript diversity, and the feasibility of cDNA cloning as a diagnostic tool for citrin deficiency,” Gene, vol. 511, no. 2, pp. 227–234, 2012.
View at: Publisher Site | Google Scholar
M. Komatsu, T. Kimura, M. Yazaki et al., “Steatogenesis in adult-onset type II citrullinemia is associated with down- regulation of PPARα,” Biochimica et Biophysica Acta, vol. 1852, no. 3, pp. 473–481, 2015.
View at: Publisher Site | Google Scholar
A. Chatrath, A. Ratan, and A. Dutta, “Germline variants predictive of tumor mutational burden and immune checkpoint inhibitor efficacy,” iScience, vol. 24, no. 3, p. 102248, 2021.
View at: Publisher Site | Google Scholar
M. Miyo, M. Konno, N. Nishida et al., “Metabolic adaptation to nutritional stress in human colorectal cancer,” Scientific Reports, vol. 6, no. 1, p. 38415, 2016.
View at: Publisher Site | Google Scholar
Z. Tang, C. Li, B. Kang, G. Gao, C. Li, and Z. Zhang, “GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses,” Nucleic acids research, vol. 45, no. W1, pp. W98–W102, 2017.
View at: Publisher Site | Google Scholar
M. Uhlen, L. Fagerberg, B. M. Hallstrom et al., “Tissue-based map of the human proteome,” Science, vol. 347, no. 6220, p. 1260419, 2015.
View at: Publisher Site | Google Scholar
C. P. Goswami and H. Nakshatri, “PROGgeneV2: enhancements on the existing database,” BMC Cancer, vol. 14, no. 1, p. 970, 2014.
View at: Publisher Site | Google Scholar
H. Mizuno, K. Kitada, K. Nakai, and A. Sarai, “PrognoScan: a new database for meta-analysis of the prognostic value of genes,” BMC Medical Genomics, vol. 2, no. 1, p. 18, 2009.
View at: Publisher Site | Google Scholar
Y. Cao, J. Xie, L. Chen et al., “Construction and validation of a novel pyroptosis-related gene signature to predict the prognosis of uveal melanoma,” Frontiers in cell and developmental biology, vol. 9, p. 761350, 2021.
View at: Publisher Site | Google Scholar
L. Xue, J. Liu, J. Xie, and J. Luo, “Prognostic value of SLC16A3 (MCT4) in lung adenocarcinoma and its clinical significance,” International journal of general medicine, vol. Volume 14, pp. 8413–8425, 2021.
View at: Publisher Site | Google Scholar
T. Sing, O. Sander, N. Beerenwinkel, and T. Lengauer, “ROCR: visualizing classifier performance in R,” Bioinformatics, vol. 21, no. 20, pp. 3940-3941, 2005.
View at: Publisher Site | Google Scholar
A. Mayakonda, D. C. Lin, Y. Assenov, C. Plass, and H. P. Koeffler, “Maftools: efficient and comprehensive analysis of somatic variants in cancer,” Genome Research, vol. 28, no. 11, pp. 1747–1756, 2018.
View at: Publisher Site | Google Scholar
K. Ellrott, M. H. Bailey, G. Saksena et al., “Scalable open science approach for mutation calling of tumor exomes using multiple genomic pipelines,” Cell systems, vol. 6, no. 3, pp. 271–281.e7, 2018.
View at: Publisher Site | Google Scholar
B. Ru, C. N. Wong, Y. Tong et al., “TISIDB: an integrated repository portal for tumor-immune system interactions,” Bioinformatics, vol. 35, no. 20, pp. 4200–4202, 2019.
View at: Publisher Site | Google Scholar
S. V. Vasaikar, P. Straub, J. Wang, and B. Zhang, “LinkedOmics: analyzing multi-omics data within and across 32 cancer types,” Nucleic acids research, vol. 46, no. D1, pp. D956–D963, 2018.
View at: Publisher Site | Google Scholar
D. Szklarczyk, J. H. Morris, H. Cook et al., “The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible,” Nucleic acids research, vol. 45, no. D1, pp. D362–D368, 2017.
View at: Publisher Site | Google Scholar
P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, no. 11, pp. 2498–2504, 2003.
View at: Publisher Site | Google Scholar
D. Warde-Farley, S. L. Donaldson, O. Comes et al., “The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function,” Nucleic acids research, vol. 38, supplement_2, pp. W214–W220, 2010.
View at: Publisher Site | Google Scholar
Y. Zhou, B. Zhou, L. Pache et al., “Metascape provides a biologist-oriented resource for the analysis of systems-level datasets,” Nature Communications, vol. 10, no. 1, p. 1523, 2019.
View at: Publisher Site | Google Scholar
X. J. Qin, X. Lin, G. Xue et al., “CXCL10 is a potential biomarker and associated with immune infiltration in human papillary thyroid cancer,” Bioscience Reports, vol. 41, no. 1, 2021.
View at: Publisher Site | Google Scholar
D. Bogunovic, D. W. O'Neill, I. Belitskaya-Levy et al., “Immune profile and mitotic index of metastatic melanoma lesions enhance clinical staging in predicting patient survival,” Proceedings of the National Academy of Sciences of the United States of America, vol. 106, no. 48, pp. 20429–20434, 2009.
View at: Publisher Site | Google Scholar
G. Jonsson, C. Busch, S. Knappskog et al., “Gene expression profiling-based identification of molecular subtypes in stage IV melanomas with different clinical outcome,” Clinical cancer research: an official journal of the American Association for Cancer Research, vol. 16, no. 13, pp. 3356–3367, 2010.
View at: Publisher Site | Google Scholar
J. Xie, L. Chen, Q. Tang et al., “A necroptosis-related prognostic model of uveal melanoma was constructed by single-cell sequencing analysis and weighted co-expression network analysis based on public databases,” Frontiers in Immunology, vol. 13, p. 847624, 2022.
View at: Publisher Site | Google Scholar
D. Szklarczyk, A. L. Gable, K. C. Nastou et al., “The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets,” Nucleic acids research, vol. 49, no. D1, pp. D605–D612, 2021.
View at: Publisher Site | Google Scholar
X. Zhou, M. Peng, Y. He et al., “CXC chemokines as therapeutic targets and prognostic biomarkers in skin cutaneous melanoma microenvironment,” Frontiers in oncology, vol. 11, p. 619003, 2021.
View at: Publisher Site | Google Scholar
B. Huang, W. Han, Z. F. Sheng, and G. L. Shen, “Identification of immune-related biomarkers associated with tumorigenesis and prognosis in cutaneous melanoma patients,” Cancer cell international, vol. 20, no. 1, p. 195, 2020.
View at: Publisher Site | Google Scholar
C. C. Wong, Y. Qian, X. Li et al., “SLC25A22 promotes proliferation and survival of colorectal cancer cells with KRAS mutations and xenograft tumor progression in mice via intracellular synthesis of aspartate,” Gastroenterology, vol. 151, no. 5, pp. 945–960.e6, 2016.
View at: Publisher Site | Google Scholar
P. M. Szczypinski and K. Sprawka, “Orthorectification of skin nevi images by means of 3D model of the human body,” Sensors, vol. 21, no. 24, p. 8367, 2021.
View at: Publisher Site | Google Scholar
J. H. Park, W. Y. Pyun, and H. W. Park, “Cancer metabolism: phenotype, signaling and therapeutic targets,” Cell, vol. 9, no. 10, p. 2308, 2020.
View at: Publisher Site | Google Scholar
A. Indini, F. Grossi, M. Mandala, D. Taverna, and V. Audrito, “Metabolic Interplay between the Immune System and Melanoma Cells: Therapeutic Implications,” Biomedicines, vol. 9, no. 6, p. 607, 2021.
View at: Publisher Site | Google Scholar
G. C. Leonardi, L. Falzone, R. Salemi et al., “Cutaneous melanoma: from pathogenesis to therapy (review),” International journal of oncology, vol. 52, no. 4, pp. 1071–1080, 2018.
View at: Publisher Site | Google Scholar
J. Meng, X. Huang, Y. Qiu, M. Yu, J. Lu, and J. Yao, “Characterization of m6A-related genes landscape in skin cutaneous melanoma to aid immunotherapy and assess prognosis,” International journal of general medicine, vol. Volume 14, pp. 5345–5361, 2021.
View at: Publisher Site | Google Scholar
Y. Liao, D. He, and F. Wen, “Analyzing the characteristics of immune cell infiltration in lung adenocarcinoma via bioinformatics to predict the effect of immunotherapy,” Immunogenetics, vol. 73, no. 5, pp. 369–380, 2021.
View at: Publisher Site | Google Scholar
X. Shi, L. Ma, J. Hao, and W. Yan, “Regulatory effects of comprehensive psychological intervention on adverse emotions and immune status of cervical cancer patients during the perioperative period,” American journal of translational research, vol. 13, no. 6, pp. 6362–6371, 2021.
View at: Google Scholar

Copyright

Copyright © 2022 Yue lv et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

532

Downloads

635

Citations