Comprehensive analysis of 33 human cancers reveals clinical implications and immunotherapeutic value of the solute carrier family 35 member A2

Abstract

Background: Solute carrier family 35 member A2 (SLC35A2), which belongs to the SLC35 solute carrier family of human nucleoside sugar transporters, has shown regulatory roles in various tumors and neoplasms. However, the function of SLC35A2 across human cancers remains to be systematically assessed. Insights into the prediction ability of SLC35A2 in clinical practice and immunotherapy response remains limited.

Materials and methods: We obtained the gene expression and protein levels of SLC35A2 in a variety of tumors from Molecular Taxonomy of Breast Cancer International Consortium, The Cancer Genome Atlas, Gene Expression Omnibus, Chinese Glioma Genome Atlas, and Human Protein Atlas databases. The SLC35A2 level was validated by immunohistochemistry. The predictive value for prognosis was evaluated by Kaplan-Meier survival and Cox regression analyses. Correlations between SLC35A2 expression and DNA methylation, genetic alterations, tumor mutation burden (TMB), microsatellite instability (MSI), and tumor microenvironment were performed using Spearman's correlation analysis. The possible downstream pathways of SLC35A2 in different human cancers were explored using gene set variation analysis. The potential role of SLC35A2 in the tumor immune microenvironment was evaluated via EPIC, CIBERSORT, MCP-counter, CIBERSORT-ABS, quanTIseq, TIMER, and xCell algorithms. The difference in the immunotherapeutic response of SLC35A2 under different expression conditions was evaluated by the tumor immune dysfunction and exclusion (TIDE) score as well as four independent immunotherapy cohorts, which includes patients with bladder urothelial carcinoma (BLCA, N = 299), non-small cell lung cancer (NSCLC, N = 72 and N = 36) and skin cutaneous melanoma (SKCM, N = 25). Potential drugs were identified using the CellMiner database and molecular docking.

Results: SLC35A2 exhibited abnormally high or low expression in 23 cancers and was significantly associated with the prognosis. In various cancers, SLC35A2 expression and mammalian target of rapamycin complex 1 signaling were positively correlated. Multiple algorithmic immune infiltration analyses suggested an inverse relation between SLC35A2 expression and infiltrating immune cells, which includes CD4+T cells, CD8+T cells, B cells, and natural killer cells (NK) in various tumors. Furthermore, SLC35A2 expression was significantly correlated with pan-cancer immune checkpoints, TMB, MSI, and TIDE genes. SLC35A2 showed significant predictive value for the immunotherapy response of patients with diverse cancers. Two drugs, vismodegib and abiraterone, were identified, and the free binding energy of cytochrome P17 with abiraterone was higher than that of SLC35A2 with abiraterone.

Conclusion: Our study revealed that SLC35A2 is upregulated in 20 types of cancer, including lung adenocarcinoma (LUAD), breast invasive carcinoma (BRCA), colon adenocarcinoma (COAD), and lung squamous cell carcinoma (LUSC). The upregulated SLC35A2 in five cancer types indicates a poor prognosis. Furthermore, there was a positive correlation between the overexpression of SLC35A2 and reduced lymphocyte infiltration in 13 cancer types, including BRCA and COAD. Based on data from several clinical trials, patients with LUAD, LUSC, SKCM, and BLCA who exhibited high SLC35A2 expression may experience improved immunotherapy response. Therefore, SLC35A2 could be considered a potential predictive biomarker for the prognosis and immunotherapy efficacy of various tumors. Our study provides a theoretical basis for further investigating its prognostic and therapeutic potentials.

Keywords: SLC35A2; glycosylation; immunotherapy; pan-cancer analysis; tumor microenvironment.