BackgroundAutophagy has an important association with tumorigenesis, progression and prognosis. However, the mechanism of autophagy-regulated genes on the risk prognosis of bladder cancer (BC) patients has not been fully elucidated yet. In this study, we created a prognostic model of BC risk based on autophagyrelated genes, which further illustrates the value of genes associated with autophagy in the treatment of BC.
MethodsWe rst downloaded human autophagy-associated genes and BC datasets from Human Autophagy Database and The Cancer Genome Atlas (TCGA) database, and nally obtained differential prognosisassociated genes for autophagy by univariate regression analysis and differential analysis of cancer versus normal tissues. Subsequently, we downloaded two datasets from Gene Expression Omnibus (GEO), GSE31684 and GSE15307, to expand the total number of samples. Based on these genes, we distinguished the molecular subtypes (C1, C2) and gene classes (A, B) of BC by consistent clustering analysis. Using the genes merged from TCGA and the two GEO datasets, we conducted least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression analysis to obtain risk genes and construct autophagy-related risk prediction models. The accuracy of this risk prediction model was assessed by Receiver operating characteristic (ROC) and calibration curves, and then Nomograms were constructed to predict the survival of bladder cancer patients at 1, 3, and 5 years, respectively. According to the median value of the risk score, we divided BC samples into high and low risk groups. Kaplan-Meier (K-M) survival analysis was performed to compare survival differences between subgroups. Then, we used single sample gene set enrichment analysis (ssGSEA) for immune cell in ltration abundance, immune checkpoint genes, immunotherapy response, gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis, and tumor mutation burden (TMB) analysis for different subgroups. We also applied quantitative real-time polymerase chain reaction (PCR) and immunohistochemistry (IHC) techniques to verify the expression of these six genes in the model. Finally, we chose the IMvigor210 dataset for external validation.
ResultsSix risk genes associated with autophagy (SPOCD1, FKBP10, NAT8B, LDLR, STM3 and ANXA2) were nally screened by lasso regression algorithm and multivariate COX regression analysis. ROC and calibration curves showed that the model established was accurate and reliable. Univariate and multivariate regression analyses were used to verify that the risk model was an independent predictor. K-M survival analysis indicated that patients in the high-risk group had signi cantly worse overall survival than those in the low-risk group. Analysis by algorithms such as correlation analysis, gene set variation Page 3/29 analysis (GSVA) and ssGSEA showed that differences in immune microenvironment, enrichment of multiple biologically active pathways, TMB, immune checkpoint genes, and human leukocyte antigens (HLA) wer...