Functional enrichment analysis of LYSET and identification of related hub gene signatures as novel biomarkers to predict prognosis and immune infiltration status of clear cell renal cell carcinoma

Abstract

Abstract Purpose The latest research shows that the lysosomal enzyme trafficking factor (LYSET) encoded by TMEM251 is a key regulator of the amino acid metabolism reprogramming (AAMR) and related pathways significantly correlate with the progression of some tumors. The purpose of this study was to explore the potential pathways of the TMEM251 in clear cell renal cell carcinoma (ccRCC) and establish related predictive models based on the hub genes in these pathways for prognosis and tumor immune microenvironment (TIME). Methods We obtained mRNA expression data and clinical information of ccRCC samples from The Cancer Genome Atlas (TCGA), E-MATE-1980, and immunotherapy cohorts. Single-cell sequencing data (GSE152938) were downloaded from the Gene Expression Omnibus (GEO) database. We explored biological pathways of the LYSET by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of TMEM251-coexpression genes. The correlation of LYSET-related pathways with the prognosis was conducted by Gene Set Variation Analysis (GSVA) and unsupervised cluster analysis. The least absolute shrinkage and selection operator (LASSO) and Cox regression were used to identify hub prognostic genes and construct the risk score. Immune infiltration analysis was conducted by CIBERSORTx and Tumor Immune Estimation Resource (TIMER) databases. The predictive value of the risk score and hub prognostic genes on immunotherapy responsiveness was analyzed through the tumor mutation burden (TMB) score, immune checkpoint expression, and survival analysis. Immunohistochemistry (IHC) was finally used to verify the expressions of hub prognostic genes. Results The TMEM251 was found to be significantly correlated with some AAMR pathways. AAGAB, ENTR1, SCYL2, and WDR72 in LYSET-related pathways were finally identified to construct a risk score model. Immune infiltration analysis showed that LYSET-related gene signatures significantly influenced the infiltration of some vital immune cells such as CD4 + cells, NK cells, M2 macrophages, and so on. In addition, the constructed risk score was found to be positively correlated with TMB and some common immune checkpoint expressions. Different predictive values of these signatures for Nivolumab therapy responsiveness were also uncovered in immunotherapy cohorts. Finally, based on single-cell sequencing analysis, the TMEM251 and the hub gene signatures were found to be expressed in tumor cells and some immune cells. Interestingly, IHC verification showed a potential dual role of four hub genes in ccRCC progression. Conclusion The novel predictive biomarkers we built may benefit clinical decision-making for ccRCC. Our study may provide some evidence that LYSET-related gene signatures could be novel potential targets for treating ccRCC and improving immunotherapy efficacy. Our nomogram might be beneficial to clinical choices, but the results need more experimental verifications in the future.