Tumor Expression Profile Analysis Developed and Validated a Prognostic Model Based on Immune-Related Genes in Bladder Cancer

A risk signature of necroptosis-related lncRNA to predict prognosis and probe molecular characteristics for male with bladder cancer

Bladder cancer (BC) is a frequently diagnosed cancer with high mortality. Male patients have a higher risk of developing BC than female patients. As a type of caspase-independent cell death, necroptosis plays a significant role in the occurrence and progression of BC. The aberrant function of long non-coding RNAs (lncRNAs) plays an indispensable role in GI. However, the relationship between lncRNA and necroptosis in male patients with BC remains unclear. The clinical information and RNA-sequencing profiles of all BC patients were retrieved from The Cancer Genome Atlas Program. A total of 300 male participants were selected for the study. We conducted to identify the necroptosis-related lncRNAs (NRLs) by Pearson correlation analysis. Subsequently, least absolute shrinkage and selection operator Cox regression were conducted to establish a risk signature with overall survival-related NRLs in the training set and to validate it in the testing set. Finally, we verified the effectiveness of the 15-NRLs signature in prognostic prediction and therapy via survival analysis, receiver operating characteristic curve analysis, and Cox regression. Furthermore, we analyzed the correlation between the signature risk score and pathway enrichment analysis, immune cell infiltration, anticancer drug sensitivity, and somatic gene mutations. We developed 15-NRLs (AC009974.1, AC140118.2, LINC00323, LINC02872, PCAT19, AC017104.1, AC134312.5, AC147067.2, AL139351.1, AL355922.1, LINC00844, AC069503.1, AP003721.1, DUBR, LINC02863) signature, and divided patients into a high-risk group and low-risk group through the median risk score. Kaplan-Meier and receiver operating characteristic curves showed that the prognosis prediction had satisfactory accuracy. Cox regression analysis indicated that the 15-NRLs signature was a risk factor independent of various clinical parameters. Additionally, immune cell infiltration, half-maximal inhibitory concentration, and somatic gene mutations differed significantly among different risk subsets, implying that the signature could assess the clinical efficacy of chemotherapy and immunotherapy. This 15-NRLs risk signature may be helpful in assessing the prognosis and molecular features of male patients with BC and improve treatment modalities, thus can be further applied clinically.

Clinical Significance and Tumor Microenvironment Characterization of a Novel Immune-Related Gene Signature in Bladder Cancer

Cancer immunotherapy plays a crucial role in bladder cancer (BC) progression. Increasing evidence has elucidated the clinicopathologic significance of the tumor microenvironment (TME) in predicting outcomes and therapeutic efficacy. This study sought to establish a comprehensive analysis of the immune-gene signature combined with TME to assist in BC prognosis. We selected sixteen immune-related genes (IRGs) after a weighted gene co-expression network and survival analysis. Enrichment analysis revealed that these IRGs were actively involved in Mitophagy and Renin secretion pathways. After multivariable COX analysis, the IRGPI comprising NCAM1, CNTN1, PTGIS, ADRB3, and ANLN was established to predict the overall survival of BC, which was validated in both TCGA and GSE13507 cohorts. In addition, a TME gene signature was developed for molecular and prognosis subtyping with unsupervised clustering, followed by a panoramic landscape characterization of BC. In summary, the IRGPI model developed in our study provided a valuable tool with an improved prognosis for BC.

Construction and validation of a bladder cancer risk model based on autophagy-related genes

Autophagy 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 autophagy-related genes, which further illustrates the value of genes associated with autophagy in the treatment of BC. We first downloaded human autophagy-associated genes and BC datasets from Human Autophagy Database and The Cancer Genome Atlas (TCGA) database, and finally obtained differential prognosis-associated 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 the 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 infiltration 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. Six risk genes associated with autophagy (SPOCD1, FKBP10, NAT8B, LDLR, STMN3, and ANXA2) were finally 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 significantly worse overall survival than those in the low-risk group. Analysis by algorithms such as correlation analysis, gene set variation analysis (GSVA), and ssGSEA showed that differences in immune microenvironment, enrichment of multiple biologically active pathways, TMB, immune checkpoint genes, and human leukocyte antigens (HLAs) were observed in the different risk groups. Then, we constructed nomograms that predicted the 1-, 3-, and 5-year survival rates of different BC patients. In addition, we screened nine sensitive chemotherapeutic drugs using the correlation between the obtained expression status of risk genes and drug sensitivity results. Finally, the external dataset IMvigor210 verified that the model is reliable and efficient. We established an autophagy-related risk prognostic model that is accurate and reliable, which lays the foundation for future personalized treatment of bladder cancer.

Identification of Immune-Related Genes Associated With Bladder Cancer Based on Immunological Characteristics and Their Correlation With the Prognosis

BACKGROUND

To identify the immune-related genes of bladder cancer (BLCA) based on immunological characteristics and explore their correlation with the prognosis.

METHODS

We downloaded the gene and clinical data of BLCA from the Cancer Genome Atlas (TCGA) as the training group, and obtained immune-related genes from the Immport database. We downloaded GSE31684 and GSE39281 from the Gene Expression Omnibus (GEO) as the external validation group. R (version 4.0.5) and Perl were used to analyze all data.

RESULT

Univariate Cox regression analysis and Lasso regression analysis revealed that 9 prognosis-related immunity genes (PIMGs) of differentially expressed immune genes (DEIGs) were significantly associated with the survival of BLCA patients (p < 0.01), of which 5 genes, including NPR2, PDGFRA, VIM, RBP1, RBP1 and TNC, increased the risk of the prognosis, while the rest, including CD3D, GNLY, LCK, and ZAP70, decreased the risk of the prognosis. Then, we used these genes to establish a prognostic model. We drew receiver operator characteristic (ROC) curves in the training group, and estimated the area under the curve (AUC) of 1-, 3- and 5-year survival for this model, which were 0.688, 0.719, and 0.706, respectively. The accuracy of the prognostic model was verified by the calibration chart. Combining clinical factors, we established a nomogram. The ROC curve in the external validation group showed that the nomogram had a good predictive ability for the survival rate, with a high accuracy, and the AUC values of 1-, 3-, and 5-year survival were 0.744, 0.770, and 0.782, respectively. The calibration chart indicated that the nomogram performed similarly with the ideal model.

CONCLUSION

We had identified nine genes, including PDGFRA, VIM, RBP1, RBP1, TNC, CD3D, GNLY, LCK, and ZAP70, which played important roles in the occurrence and development of BLCA. The prognostic model based on these genes had good accuracy in predicting the OS of patients and might be promising candidates of therapeutic targets. This study may provide a new insight for the diagnosis, treatment and prognosis of BLCA from the perspective of immunology. However, further experimental studies are necessary to reveal the underlying mechanisms by which these genes mediate the progression of BLCA.