Three Prognostic Biomarkers Correlate with Immune Checkpoint Blockade Response in Bladder Urothelial Carcinoma

Guo, Ya; Zhang, Yin Bin; Li, Yi; Su, Wang Hui; He, Shan; Pan, Shu Pei; Xu, Kun; Kou, Wei Hua

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

International Journal of Genomics

On this page

Abstract Background Methods Results Discussion Conclusions Abbreviations Data Availability Ethical Approval Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

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

Three Prognostic Biomarkers Correlate with Immune Checkpoint Blockade Response in Bladder Urothelial Carcinoma

Ya Guo,¹Yin Bin Zhang,²Yi Li,¹Wang Hui Su,¹Shan He,³Shu Pei Pan,¹Kun Xu,¹and Wei Hua Kou¹

Academic Editor: Luca Falzone

Received15 Mar 2022

Accepted09 May 2022

Published26 May 2022

Abstract

Aim. We aim to develop a signature that could accurately predict prognosis and evaluate the response to immune checkpoint blockade (ICB) in bladder urothelial carcinoma (BLCA). Methods. Based on comprehensive analysis of public database, we identified prognosis-related hub genes and investigated their predictive values for the ICB response in BLCA. Results. Among 69 common DEGs, three genes (AURKA, BIRC5, and CKS1B) were associated with poor prognosis, and which were related to histological subtypes, TP53 mutation status, and the C2 (IFN-gamma dominant) subtype. Three genes and their related risk model can effectively predict the response of immunotherapy. Their related drugs were identified through analysis of drug bank database. Conclusions. Three genes could predict prognosis and evaluate the response to ICB in BLCA.

1. Background

Bladder cancer is the ninth most common malignant tumor worldwide, and bladder urothelial carcinoma (BLCA) is its most frequent pathological type [1]. Despite the development of diagnostic and treatment techniques, the 5-year survival rate of patients with BLCA varies from 5% to 70%, and their prognosis remains unfavorable [2]. Detection of biomarkers can indicate a particular disease state and can be used in screening for differential prognosis, evaluation of treatment response, and monitoring of disease progression [3]. As demonstrated in different studies, the availability of several cancer databases and a comprehensive bioinformatics analysis has allowed for the accurate identification of key biomarkers for early diagnosis of malignancies or prediction of prognosis and cancer recurrence [4–7]. Therefore, it is reliable to develop useful biomarkers for predicting BLCA prognosis based on systematic bioinformatics analysis.

BLCA is the 10th most common malignant tumor worldwide, which is divided into nonmuscle-invasive BLCA (NMIBC) and muscle-invasive (MIBC). MIBC has an unfavorable prognosis with a 5-year overall survival (OS) of approximately 50% [8, 9]. In recent years, cancer immunotherapies by immune checkpoint blockade (ICB) have attracted considerable attention for its influence on the treatment of locally advanced and metastatic BLCA [10]. Recent studies reported that immune checkpoint blockade (ICB) has promising improved survival rates in individuals with advanced BLCA with a high tumor mutation burden (TMB) and neoantigens. Patients with tumors that overexpression of programmed cell death 1 receptor (PD-1) or its ligand (PD-L1) appear to benefit most from this therapy [9, 11]. However, only a proportion of patients respond to treatment with ICB, the results of immunotherapy are still not satisfactory [12, 13]. Therefore, immunotherapy needs to continue to improve in terms of improvement in their effectiveness in the treatment of BLCA.

Previous studies have reported that highly infiltrating lymphocytes are related to the prognosis of BLCA [14]. The evidence indicates that CD8+ T cells are involved in tumor adaptive immunity, and their infiltration is associated with prognosis [15, 16]. Programmed death ligand 1 (PD-L1), also known as B7-H1 or CD274, is involved in inhibiting T cell-mediated antitumor immunity through interaction with PD-1 [17, 18]. Previous studies have reported that high expression of PD-L1 is associated with worse cancer outcomes in various malignancies [19]. CD8+ T cell infiltration and tumor mutation burden (TMB) have been reported to be correlated with response to atezolizumab (anti-PD-L1) in metastatic urothelial cancer (mUC) [20–22]. However, none of these factors is sufficient to achieve accurate outcome prediction, and identification of ICB response biomarkers is a critical challenge in the field [23]. In general, there is an urgent need to identify key biomarkers that can be used to effectively evaluate the response to immunotherapy in BLCA patients.

In this study, we first identified DEGs from the GEPIA (Gene Expression Profiling Interactive Analysis) and Oncomine databases and then identified the overlapping DEGs between them using a Venn diagram. We further performed gene enrichment and protein–protein interaction (PPI) analyses to select hub genes. Moreover, prognosis-related hub genes were identified using comprehensive bioinformatics analysis and confirmed in three online databases. We then explored the expression of the key prognosis-related genes with different clinical factors using the UALCAN and TISIDB database. Finally, we evaluated whether the risk model based on three key genes could be used to predict immunotherapy response in BLCA.

2. Methods

2.1. Identification of DEGs

GEPIA (http://gepia.cancer-pku.cn/) is an online network tool based on data from The Cancer Genome Atlas (TCGA) and GTEx, which can be used to study interactions between DEGs, as well as for survival analysis, profile plotting, and detection of similar genes. Oncomine is an online resource containing microarray data (https://www.Oncomine.org) [24]. In this study, we first used the GEPIA and Oncomine databases to identify DEGs by comparison of tumor samples with normal samples. RNA-Seq data from the TCGA-based (The Cancer Genome Atlas) GEPIA database was performed to identify differentially expressed genes (DEGs) between 404 tumor and 19 normal tissues. ANOVA method was used to obtain the differently expressed genes (DEGs), mRNAs with and were selected as DEGs. In the Oncomine database, 288 samples from TCGA gene expression dataset were used to screen DEGs between tumor and normal groups. The selection criteria for DEGs are , , and . Then, we identified the overlapping DEGs between them using a Venn diagram (http://bioinformatics.psb.ugent.be/webtools/Venn/) [25].

2.2. Functional Analysis and Pathway Enrichment Analysis

Metascape (http://metascape.org/) is an online resource for gene annotation and analysis [26]. In the present study, Metascape was used to perform gene ontology (GO) and pathway analyses of 69 common hub genes. Pathway and process enrichment analyses were conducted based on several sources, including GO biological processes, The Kyoto Encyclopedia of Genes and Genomes pathways, reactome gene sets, and CORUM. Terms with a value less than 0.01, a minimum count of 3, and an enrichment factor greater than 1.5 were considered to represent significant processes or pathways.

2.3. Construction of PPI Network and Identification of Hub Genes

We evaluated PPI information of common genes using the STRING online database (https://string-db.org/cgi/input.pl) and then visualized the resulting interaction network using Cytoscape software (http://www.cytoscape.org/) [27, 28]. A confidence score greater than 0.4 was defined as significant. The Molecular Complex Detection (MCODE) plugin in Cytoscape was used to further screen key genes in the PPI network with , , and .

2.4. Development of Prognosis-Related Model

We downloaded gene expression profile and clinical data from TCGA (https://portal.gdc.cancer.gov/). To reduce statistical bias, BLCA patients were excluded if clinical information or overall survival (OS) was missing from their records. The prognostic value of 11 identified hub DEGs that were analyzed using the R survival package [29]. The DEGs with significant prognostic value were selected for further analysis. Based on these prognosis-related DEGs, least absolute shrinkage and selection operator (LASSO) Cox regression analysis was applied to establish prognostic model. Patients were then divided into high- and low-risk groups according to median risk score. Receiver operating characteristic (ROC) curve was used to assess the predictive accuracy of each gene and risk score. Univariate and multivariate Cox regression analysis was performed for independent analysis with other clinical characteristics. Nomogram was then used to assess 1-year, 3-year, and 5-year overall survival.

2.5. Validation of Prognostic Value of Three Key DEGs

PROGgenesV2 (http://genomics.jefferson.edu/proggene/filter.php) is a web resource that allows researchers to study the correlations between genes and overall survival (OS) in multiple cancers based on TCGA and GEO data [30]. PrognoScan (http://www.prognoscan.org/) was used to evaluate the associations between gene expression and patient prognosis, according to measures including OS and disease-free survival (DFS), across a large collection of publicly available cancer microarray datasets [31]. The OSblca database (http://bioinfo.henu.edu.cn/BLCA/BLCAList.jsp) provides a useful tool to assess novel prognostic biomarkers in bladder cancer, based on data from 1,075 bladder cancer patients, including OS, disease-specific survival (DSS), disease-free interval, and progression-free interval [32]. In this study, we further confirmed the prognostic value of key genes based on the abovementioned three databases. The hub genes identified in this way were defined as key prognosis-related genes.

2.6. Association between Prognosis-Related Key Genes and Clinical Characteristics

The UALCAN database (http://ualcan.path.uab.edu/) is a comprehensive web resource for analyzing cancer OMICS data (TCGA and MET500) [33, 34]. TISIDB (http://cis.hku.hk/TISIDB) is a publicly available resource that allows the user to explore the function of a gene and its role in tumor-immune features. TISIDB consists of 10 modules: function, literature, screening, immunotherapy, lymphocyte, immunomodulator, chemokine, subtype, clinical, and drug [35]. A previous study reported that the expression and prognostic values of DEGs were associated with clinical characteristics, including TNM stage, smoking history, lymph invasion, histological type, and immune subtype [5, 36]. Therefore, we explored the relationship between three key prognosis-related genes and clinical characteristics using UALCAN and TISIDB database.

2.7. Immune Cell Infiltration Analysis

TIMER is a user-friendly web portal for the systematic analysis of immune infiltrates across different types of cancer (https://cistrome. shinyapps.io/timer/) [37]. In this study, we used TIMER and TISIDB to analyze the associations between the three key genes and immune cell infiltration. was considered as statistically significant.

2.8. Evaluation of the Value of Prognostic-Related Genes in Response to Immune Checkpoint Blockade

Previous studies have suggested that TMB and PD-L1 expressions are correlated with response to atezolizumab in mUC [20]. In this study, the correlation between prognostic-related gene expression and TMB score was calculated using Spearman’s correlation [24, 38, 39]. The immunomodulator module of TISIDB was used to examine the associations between PDL1 and selected genes. Moreover, we used the screening module of TISIDB to explore whether the expression of prognosis-related genes showed significant differences between responders and nonresponders to immunotherapy. Tumor Immune Dysfunction and Exclusion (TIDE, http://tide.dfci.harvard.edu) algorithm was then performed to estimate custom biomarker predictive power of response outcome and overall survival [23, 40].

2.9. Exploration of the Model in the Tumor Immune Microenvironment and Immunotherapeutic Treatment

RNA sequence transcription data, mutation data, and relevant clinical information of BLCA patients were obtained from the TCGA (https://cancergenome.nih.gov/) database [41]. Immune function was analyzed based on three prognostic gene models using R package limma, GSVA, GSEABase, pheatmap, and reshape2 [42]. We used the TIDE algorithm to predict immunotherapy response [43].

2.10. Correlation between Prognostic-Related Genes and Their Target Drug

Gene Set Cancer Analysis (GSCALite, http://bioinfo.life.hust.edu.cn/web/GSCALite/) is a web server for the analysis a set of genes in cancers with different function modules [44]. In this study, we analyzed the drugs in the Drug Bank database targeting these three genes using the drug module of TISIDB. We further applied GSCALite to analyze the drug sensitivity of key ceRNA signatures.

3. Results

3.1. Identification of Hub Genes

A total of 750 DEGs were identified from the GEPIA database, and 1,881 DEGs were identified from Oncomine. Sixty-nine common genes were screened out using the Venn diagram (Supplementary 1, Supplementary 2, Figure 1(a)). We then performed GO and pathway enrichment analyses for the common genes using Metascape. The results showed that these common genes were involved in 20 main GO terms and pathways, of which the top five were mitotic cell cycle, cGMP-PKG signaling pathway, extracellular matrix organization, muscle contraction, and response to hydrogen peroxide (Figure 1(b)). Finally, a PPI network of DEGs was constructed using the STRING and Cytoscape software, containing 69 nodes and 83 edges (Figure 1(c)). One significant module was identified using MCODE. This module contained 11 genes, Aurora-A kinase (AURKA), BIRC5, CENPA, CKS1B, ECT2, MYBL2, NUF2, RRM2, TK1, TPX2, and UBE2C, which were defined as hub genes (Figure 1(d)).

(a)

(b)

(c)

(d)

Figure 1

Identification of 11 hub genes. (a) Identification of common genes between GEPIA and Oncomine by Venn diagram. (b) Enriched terms of common genes identified by Metascape. Network of enriched terms colored by cluster ID. (c) PPI network of DEGs constructed with STRING software: nodes represent proteins; continuous lines represent direct interactions (physical), while indirect ones (functional) are represented by interrupted lines; and line thickness indicates the strength of data support. (d) Identification of hub genes using MCODE. Upregulated genes are represented by red nodes, while downregulated genes are denoted by green nodes. Node size is positively correlated with value. The line color is determined by the combined score provided by STRING.

3.2. Construction of Prognostic Gene Model

The association between the expression of 11 hub genes and overall survival was evaluated using the R survival package. High expression of three genes (AURKA, BIRC5, and CKS1B) was related to an unfavorable prognosis (Figures 2(a)–2(k)). LASSO analysis was applied to establish a prognostic gene model based on these three prognostic DEGs (Figures 3(a) and 3(b)). The . The BLCA patients were divided into high- and low-risk score group based on risk score. Figure 3(c) displayed the risk score distribution, survival status, and the expression of these genes. BLCA patients with high-risk score had an unfavorable prognosis than those with low-risk score (Figure 3(d)). ROC curve analysis indicated that the AUCs of the 1-, 3-, and 5-year prognosis models were 0.593, 0.556, and 0.54, respectively (Figure 3(e)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(a)

(b)

(c)

(d)

(e)

Figure 3

Construction of a prognostic model based on three key gene. (a) LASSO coefficient profiles of the three key genes. (b) Plots of the ten-fold cross-validation error rates. (c) Distribution of risk score, survival status, and the expression of three prognostic genes in BLCA. (d) Overall survival curves for BLCA patients in the high-/low-risk group. (e) ROC analysis was performed to measure the predictive value. BLCA: bladder urothelial carcinoma; LASSO: least absolute shrinkage and selection operator; ROC: receiver operating characteristic curve; HR: hazard ratio; CI: confidence interval.

3.3. Construction of Predictive Nomogram

According to the clinicopathologic features and three prognostic genes, we constructed a predictive nomogram to predict the survival probability (Figures 4(a) and 4(b)). The C-index of the nomogram was 0.624 (95% CI, 0.582-1). Calibration curves also showed a favorable predictive power of the nomogram (Figures 4(c) and 4(d)). Furthermore, we validated the prognostic value of the three genes in BLCA using the PROGgenesV2, PrognoScan, and OSblca databases. Our results showed that BLCA patients with higher expression levels of the three hub genes exhibited poorer OS, DFS, and DSS, indicating that the three hub genes may be associated with unfavorable prognosis (Table 1). In summary, our data suggested that the three key genes could serve as biomarkers of poor prognosis.

(a)

(b)

(c)

(d)

3.4. Association between Three Key Genes and Clinical Parameters in BLCA

A previous study reported that the OS of BLCA patients was significantly associated with clinical characteristics, including TNM stage, smoking history, lymph invasion, and histological type [32]. The three hub genes identified here were associated with various clinical characteristics including age, histological subtype, molecular subtype, nodal metastasis status, sample type, smoking, cancer stage, and TP53 mutation status (Table 2). Most importantly, overexpression of the three hub genes was positively correlated with histological subtypes (Figures 5(a)–5(c)). The expression of the three hub genes was higher in the basal squamous and neuronal subtypes than in the luminal subtype (Figures 5(d)–5(f)). BLCA patients with TP53 mutations also showed high expression of the three hub genes (Figures 5(g)–5(i)). Moreover, we found that the three prognosis-related genes had the highest expression levels in the C2 (IFN-gamma dominant) subtype and the lowest in the C3 (inflammatory) subtype (Figures 5(j)–5(l)). Taken together, these results suggest that increased expression of the three key genes might predict poor prognosis in patients with BLCA.

Table 2

Relationships of three prognostic genes with clinical characteristics. Clinical characteristics included age, histological subtype (papillary or nonpapillary tumor), molecular subtype (luminal papillary, luminal infiltrated, luminal, basal squamous, and neuronal), nodal metastasis status (N0: no regional lymph node metastasis; N1: metastases in 1–3 axillary lymph nodes; N2: metastases in 4–9 axillary lymph nodes; N3: metastases in 10 or more axillary lymph nodes), sample type, smoking, cancer stage, and TP53 mutation status.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

Figure 5

Evaluation of the association between three prognosis-related key genes and clinical factors. Expression of three key prognosis-related genes based on different sample types, according to (a)–(c) histological subtypes, (d)–(f) molecular subtypes, and (g)–(i) TP53 mutation status. (j)–(l) Associations between expression of three prognosis-related genes and immune subtypes across BLCA via TISIDB database: (j) AURKA, (k) BIRC5, and (l) CKS1B. Act_CD8: activated CD8 T cell; Tcm_CD8: central memory CD8 T cell; Tem_CD8: effector memory CD8 T; C1: wound healing; C2: IFN-γ dominant; C3: inflammatory; C4: lymphocyte depleted; C5: immunologically quiet; C6: TGF-β dominant. Kruskal–Wallis test was used to evaluate the statistical significance of differential expression.

3.5. CD8+ T Cell Infiltration Predicts Poor Prognosis in BLCA

The TIMER and TISIDB databases were used to explore the relationship between the three prognosis-related genes and tumor-infiltrating immune cells. The three genes were positively associated with levels of infiltrating CD8+ T cells, neutrophils, and dendritic cells. Expression of BIRC5 was negatively correlated with infiltration of B cells (Table 3, Figures 6(a)–6(f)). High levels of infiltration of CD8+ T cells were associated with poor prognosis (Figure 6(g)). These results suggest that these genes may affect prognosis via regulation of CD8+ T cells.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

3.6. Three Key Genes Correlated with ICB Response

Our results indicated that the expression levels of the three hub genes were positively correlated with infiltration of PD-L1 (CD274) expression and TMB (Figures 7(a)–7(f)). We also found that high infiltration of , high expression of CD274, and high TMB were associated with prolonged overall survival after anti-PD-L1 therapy in BLCA (Figures 7(g)–7(i)). The association score of ICB survival outcome showed that three genes were correlated with ICB benefit in different cancer, especially in BLCA (Figures 8(a)–8(d)). Finally, our results showed that the three hub genes exhibited a significant difference in expression between responders and nonresponders to atezolizumab in urothelial cancer via TISIDB analysis (Table 4). Taken together, these results suggest that the impact of the three hub genes on response to immunotherapy in BLCA may be associated with TMB and PD-L1 expression. Our identified three key genes might serve as an indicator to evaluate response to ICB immunotherapy.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

Figure 8

Estimation of the TME and cancer immunotherapy response and evaluation target drugs of these genes. Association between three genes and immunosuppressive indices (columns), including T cell dysfunction score and ICB survival outcome. (b)–(d) Survival differences between groups with high and low expression of biomarkers after anti-PDL1. (b) AURKA, (c) BIRC5, and (d) CKS1B. (e) Immune function difference between high-risk groups and low-risk group. (f) The difference in immunotherapy response between high- and low-risk groups based on the TIDE score. (g)–(i) Drugs targeting three key genes were collected from the DrugBank database: (g) AURKA, (h) BIRC5, and (i) CKS1B. Red rectangle represents the current gene, blue rectangle represents a drug, and green rectangle indicates other targets. (j, k) Correlation of three prognostic genes expression and drug sensitivity from GSCA database. (j) CTRP drug sensitivity and expression correlation. (k) GDSC drug sensitivity and expression correlation. GSCA: Gene Set Cancer Analysis; CTRP: Cancer Therapeutics Response Portal; GDSC: Genomics of Drug Sensitivity in Cancer; ICB: immune checkpoint blockade; OS: overall survival; TME: tumor immune microenvironment.

3.7. Evaluation of the Tumor Immune Microenvironment and Immunotherapy Response Based on Prognostic-Related Model

To explore the underlying molecular mechanisms of three gene-related model, we performed immune function enrichment analysis. The low-risk and high-risk groups displayed significant differences in the expression of immune indicators, including type II IFN response and MHC class I (Figure 8(e)). We further evaluated whether the risk model based on three key genes could be used to predict immunotherapy response in BLCA. Our result showed that the TIDE score of the high-risk group is less than that of the low-risk group, indicating that compared with the low-risk group, the high-risk group can benefit more from immunotherapy (Figure 8(f)).

3.8. Drug–Gene Interaction Network

Drugs targeting three key genes were collected from the Drug Bank database. AURKA and 18 other targets were correlated with 15 drugs. BIRC5 and two targets were related to 3 drugs. CKS1B and another 3 targets interacted with 2 drugs (Figures 8(g)–8(i)). Based on GSCALite analysis, we found that most drugs were effective in association with increased expression of AURKA both in CTRP and GDSC database, while BIRC5 was positively regulated by 2 drugs and negatively regulated by 2 drugs in GDSC database. Specifically, these molecules could be exploited as potential therapeutic drug targets for BLCA (Figures 8(j) and 8(k)).

4. Discussion

BLCA is one of the most common urinary cancers in the world, with high recurrence and mortality rates that limit the efficacy of treatment [45]. It is therefore essential to understand the molecular mechanism of BLCA. With the development of bioinformatics technology, increasing numbers of studies are using bioinformatics analysis to develop biomarkers and explore the molecular mechanism of BLCA [46, 47]. However, there are still no reliable biomarkers associated with prognosis of BLCA patients. In recent years, immunotherapy has attracted considerable attention owing to its influence on the treatment of locally advanced and metastatic BLCA, but the objective response rate of BLCA to immune checkpoint inhibitors (ICIs) was low [12, 48]. Therefore, it is critical to identify satisfactory signatures to effectively predict prognosis and evaluate the benefit of immunotherapy in BLCA patients.

AURKA is a member of the serine/threonine kinase family and has a role in regulation of the cell cycle [49]. Accumulating evidence indicates that AURKA is overexpressed in various cancers, including breast cancer, head and neck cancer, esophagus cancer, hematological malignancies, colorectal cancer, stomach cancer, pancreatic cancer, and ovarian and prostate cancers [50, 51]. Pathological overexpression of AURKA is correlated with shorter survival of cancer patients. According to previous reports, high expression of Aurora-A in tumor cells is closely related to poor prognosis [46, 51]. BIRC5 is a member of the inhibitor of apoptosis gene family, which has dual roles in promoting cell proliferation and preventing apoptosis [50]. Overexpression of BIRC5 has been reported in several malignancies, and higher BIRC5 expression was also found to be associated with decreased survival [52]. CKS1B is an oncogene that has been reported to show increased expression in various tumors [53]. In accordance with previous studies, our results demonstrated that high expression of three hub genes was associated with poor prognosis (Figures 2(a), 2(b), and 2(d)). Further analysis revealed that these three key genes were overexpressed in nonpapillary tumors, the basal squamous subtype, TP53 mutation patients, and C2 (IFN-gamma dominant) subtype (Figures 5(a)–5(l)). Three hub genes were highly correlated with TMB and PD-L1 expressions (Figures 7(a)–7(f)). Previous studies have reported that TP53 mutation is associated with poor prognosis [54]. TMB has prognostic roles in various cancer types, including BLCA [55]. High expression of PD-L1 is associated with worse cancer outcomes [19]. A previous study reported that the C2 subtype had the highest levels of , as well as having less favorable outcomes [36]. These results indicated that the expression of these three key genes may represent a marker of poor prognosis in BLCA, as well as in nonpapillary tumors, the basal squamous subtype, TP53 mutation patients, tumor with high TMB, and the C2 subtype.

Tumor immunity is an extremely complex biological process, and factors affecting the efficacy of immune checkpoint inhibitors (IMCIs) are associated with PD-L1 expression level, tumor-infiltrating lymphocytes (TILs), and tumor mutational burden (TMB) [13]. TMB has been investigated in various malignancies and found to be correlated with response to atezolizumab in mUC [20–22]. Increased infiltration has been reported to be correlated with better immunotherapeutic effect [56]. A previous study showed that biomarkers related to infiltration could facilitate the monitoring of immunotherapy response and the exploration of the immune infiltration mechanism in clear cell renal cell carcinoma. A previous study showed that biomarkers related to infiltration could facilitate the monitoring of immunotherapy response and the exploration of the immune infiltration mechanism in clear cell renal cell carcinoma. AURKA is overexpressed in cancer cells but not in normal tissues, making it a potential target for immunotherapy. AURKA-specific CD8(+) T cells can selectively lyse leukemia cells [21]. (CTLs) generated in vitro can recognize AURKA epitope (YLILEYAPL). In addition, these CTLs can kill leukemia cells endogenously expressing AURKA, indicating that the homologous epitopes are naturally processed and presented at a level sufficient for immunotherapeutic applications [57]. BIRC5 is a member of the apoptosis inhibitor gene family, which regulates several cancers by activating cell apoptosis process [58]. BIRC5 is correlated with T cell survival and proliferation, which can increase the accumulation and persistence of CD8(+) T cells following an encounter with Ag [59]. CKS1B is associated with cell cycle. The association between CKS1B and ICI is rare.

To further explore the molecular mechanism of three genes related model. We found that immune function was different between high- and low-risk groups (Figure 8(e)). Our result indicated that three hub genes were highly correlated with , TMB, and PD-L1 expression (Figures 6(a)–6(f) and 7(a)–7(f)). We further found that , TMB, and PD-L1 expression were positively correlated with OS after anti-PDL1 therapy (Figures 7(g)–7(i)). The association score of ICB survival outcome showed that three genes were correlated with ICB benefit (Figures 8(a)–8(d)). In addition, we found that three hub genes exhibited a significant difference in expression between responders and nonresponders to atezolizumab in urothelial cancer via TISIDB analysis (Table 4). Studies showed that TIDE score can effectively predict the response of immunotherapy [23, 42]. In this study, TIDE score analysis suggested that patients with high-risk score get a better response to immunotherapy. Taken together, we can infer that the expression of three key genes and their related model may provide reliable biomarkers to evaluate the response to immunotherapy, and the mechanism of three genes in ICB response remains to be further studied. Finally, we developed potential drug targets that interact with three genes (Figures 8(g)–8(k)). Collectively, these data suggest that three prognostic genes were correlated with ICB response in BLCA, which may be associated with , TMB, and PD-L1 expression. Drug–gene interaction network analysis indicated that these genes and their related drugs could be used in the development of new targets for BLCA immunotherapy.

5. Conclusions

In summary, three key genes in BLCA were found to be correlated with poor prognosis and immunotherapy response in BLCA. Three prognostic genes and their related drugs may help to develop new targets for improving BLCA immunotherapy. Combination of three genes’ inhibitor and anti-PDL1 may provide new insights for improving effectiveness of immunotherapy. Nevertheless, the present study has certain limitations. For example, our experimental design mainly focuses on the computational nature. In fact, no validation analysis has been performed based on BLCA cell lines or clinical samples, which greatly limit the impact of our results. Moreover, the biological mechanism of three gene-related model has not been fully elucidated. We will perform external experiments based on BLCA cell lines or clinical samples to support our results and investigate the underlying molecular mechanism of three genes in prediction immunotherapy response in our following work.

Abbreviations

BLCA:	Bladder urothelial carcinoma
DFS:	Disease-free survival
DSS:	Disease-specific survival
GEPIA:	Gene expression profiling interactive analysis
GO:	Gene ontology
IMCIs:	Immune checkpoint inhibitors
mUC:	Metastatic urothelial cancer
OS:	Overall survival
PPI:	Protein–protein interaction
TIIC:	Tumor-infiltrating immune cell
TMB:	Tumor mutation burden
TIDE:	Tumor immune dysfunction and exclusion
TME:	Tumor immune microenvironment.

Data Availability

The data that support the findings of this study are included in the article/Supplementary Material, and further inquiries can be directed to the corresponding author.

Ethical Approval

This research was supported by the Independent Ethics Committee (IEC) of the Second Affiliated Hospital of Xi’an Jiao Tong University.

Conflicts of Interest

The authors declare that there are no competing interests associated with this manuscript.

Authors’ Contributions

G. wrote the main text of the article and designed the experiments. Y.B.Z. reviewed and edited the article. Y. L, W.H.S., and S.P.P. prepared Figures 1–4. S.H. analyzed part of the data. K.X. and W.H.K. prepared Figures 5–8 and table. All authors reviewed the manuscript. Ya Guo and Yin Bin Zhang are co-first authors.

Acknowledgments

Special thanks are due to Zhong Guo Liang for bioinformatics-related technical support. This work was supported by the Shaanxi Provincial Natural Science Foundation (grant number 2017JQ8057).

Supplementary Materials

Supplementary 1: different expression genes from GEPIA database. Supplementary 2: different expression genes from Oncomine database. Supplementary 3: original source data: the relevant code or R software. (Supplementary Materials)

References

M. Kang, K.-H. Lee, H. S. Lee et al., “Concurrent autophagy inhibition overcomes the resistance of epidermal growth factor receptor tyrosine kinase inhibitors in human bladder cancer cells,” International Journal of Molecular Sciences, vol. 18, no. 2, p. 321, 2017.
View at: Publisher Site | Google Scholar
D. Liu, B. Zhou, and R. Liu, “An RNA-sequencing-based transcriptome for a significantly prognostic novel driver signature identification in bladder urothelial carcinoma,” Peer J, vol. 8, article e9422, 2020.
View at: Publisher Site | Google Scholar
A. Hentschel, E. E. V. D. Toom, A. N. Vis et al., “A systematic review on mutation markers for bladder cancer diagnosis in urine,” BJU International, vol. 127, no. 1, pp. 12–27, 2021.
View at: Publisher Site | Google Scholar
K. Jin, S. Qiu, D. Jin et al., “Development of prognostic signature based on immune-related genes in muscle-invasive bladder cancer: bioinformatics analysis of TCGA database,” Aging, vol. 13, no. 2, pp. 1859–1871, 2021.
View at: Publisher Site | Google Scholar
J. L. Deng, Y. H. Xu, and G. Wang, “Identification of potential crucial genes and key pathways in breast cancer using bioinformatic analysis,” Frontiers in Genetics, vol. 10, p. 695, 2019.
View at: Publisher Site | Google Scholar
D.-P. Huang, Y.-H. Zeng, W.-Q. Yuan et al., “Bioinformatics analyses of potential miRNA-mRNA regulatory axis in HBV-related hepatocellular carcinoma,” International Journal of Medical Sciences, vol. 18, no. 2, pp. 335–346, 2021.
View at: Publisher Site | Google Scholar
L. Falzone, S. Candido, R. Salemi et al., “Computational identification of microRNAs associated to both epithelial to mesenchymal transition and NGAL/MMP-9 pathways in bladder cancer,” Oncotarget, vol. 7, no. 45, pp. 72758–72766, 2016.
View at: Publisher Site | Google Scholar
Y. Fu, S. Sun, J. Bi, C. Kong, and L. Yin, “A novel immune-related gene pair prognostic signature for predicting overall survival in bladder cancer,” BMC Cancer, vol. 21, no. 1, p. 810, 2021.
View at: Publisher Site | Google Scholar
J. W. Van Hattum, B.-M. De Ruiter, J. R. Oddens, M. C. C. M. Hulshof, T. M. De Reijke, and A. D. Bins, “Bladder-sparing chemoradiotherapy combined with immune checkpoint inhibition for locally advanced urothelial bladder cancer-a review,” Cancers, vol. 14, no. 1, p. 38, 2022.
View at: Publisher Site | Google Scholar
D. A. Braun, K. P. Burke, and E. M. V. Allen, “Genomic approaches to understanding response and resistance to immunotherapy,” Clinical Cancer Research, vol. 22, no. 23, pp. 5642–5650, 2016.
View at: Publisher Site | Google Scholar
L. Duan, X. Liu, Z. Luo et al., “G-protein subunit gamma 4 as a potential biomarker for predicting the response of chemotherapy and immunotherapy in bladder cancer,” Genes, vol. 13, no. 4, p. 693, 2022.
View at: Publisher Site | Google Scholar
X. Chen, Y. Jin, L. Gong et al., “Bioinformatics analysis finds immune gene markers related to the prognosis of bladder cancer,” Frontiers in Genetics, vol. 11, p. 607, 2020.
View at: Publisher Site | Google Scholar
S. Pan, Y. Zhan, X. Chen, B. Wu, and B. Liu, “Bladder cancer exhibiting high immune infiltration shows the lowest response rate to immune checkpoint inhibitors,” Frontiers in Oncology, vol. 9, p. 1101, 2019.
View at: Publisher Site | Google Scholar
B. Lin, T. Zhang, X. Ye, and H. Yang, “High expression of EMP1 predicts a poor prognosis and correlates with immune infiltrates in bladder urothelial carcinoma,” Oncology Letters, vol. 20, no. 3, pp. 2840–2854, 2020.
View at: Publisher Site | Google Scholar
S. M. A. Mahmoud, E. C. Paish, D. G. Powe et al., “Tumor-infiltrating CD8+ lymphocytes predict clinical outcome in breast cancer,” Journal of Clinical Oncology : Official Journal of the American Society of Clinical Oncology, vol. 29, no. 15, pp. 1949–1955, 2011.
View at: Publisher Site | Google Scholar
P. Sharma, Y. Shen, S. Wen et al., “CD8 tumor-infiltrating lymphocytes are predictive of survival in muscle-invasive urothelial carcinoma,” Proceedings of the National Academy of Sciences of the United States of America, vol. 104, no. 10, pp. 3967–3972, 2007.
View at: Publisher Site | Google Scholar
L. B. Jilaveanu, B. Shuch, C. R. Zito et al., “PD-L1 expression in clear cell renal cell carcinoma: an analysis of nephrectomy and sites of metastases,” Journal of Cancer, vol. 5, no. 3, pp. 166–172, 2014.
View at: Publisher Site | Google Scholar
S. L. Topalian, C. G. Drake, and D. M. Pardoll, “Immune checkpoint blockade: a common denominator approach to cancer therapy,” Cancer Cell, vol. 27, no. 4, pp. 450–461, 2015.
View at: Publisher Site | Google Scholar
E. Xylinas, B. D. Robinson, L. A. Kluth et al., “Association of T-cell co-regulatory protein expression with clinical outcomes following radical cystectomy for urothelial carcinoma of the bladder,” European journal of surgical oncology : the journal of the European Society of Surgical Oncology and the British Association of Surgical Oncology, vol. 40, no. 1, pp. 121–127, 2014.
View at: Publisher Site | Google Scholar
A. Balar, M. Galsky, J. Rosenberg et al., “Atezolizumab as first-line treatment in cisplatin-ineligible patients with locally advanced and metastatic urothelial carcinoma: a single-arm, multicentre, phase 2 trial,” Lancet (London, England), vol. 389, no. 10064, pp. 67–76, 2017.
View at: Publisher Site | Google Scholar
J. Lin, M. Yu, X. Xu et al., “Identification of biomarkers related to CD8 T cell infiltration with gene co-expression network in clear cell renal cell carcinoma,” Aging, vol. 12, no. 4, pp. 3694–3712, 2020.
View at: Publisher Site | Google Scholar
J. Rosenberg, J. Hoffman-Censits, T. Powles et al., “Atezolizumab in patients with locally advanced and metastatic urothelial carcinoma who have progressed following treatment with platinum-based chemotherapy: a single-arm, multicentre, phase 2 trial,” Lancet (London, England), vol. 387, no. 10031, pp. 1909–1920, 2016.
View at: Publisher Site | Google Scholar
P. Jiang, S. Gu, D. Pan et al., “Signatures of T cell dysfunction and exclusion predict cancer immunotherapy response,” Nature Medicine, vol. 24, no. 10, pp. 1550–1558, 2018.
View at: Publisher Site | Google Scholar
M. K. Hassan, D. Kumar, M. Naik, and M. Dixit, “The expression profile and prognostic significance of eukaryotic translation elongation factors in different cancers,” PLoS One, vol. 13, no. 1, article e0191377, 2018.
View at: Publisher Site | Google Scholar
Y. Niu, B. Hu, X. Li et al., “Comparative digital gene expression analysis of tissue-cultured plantlets of highly resistant and susceptible banana cultivarsin response to Fusarium oxysporum,” International Journal of Molecular Sciences, vol. 19, no. 2, p. 350, 2018.
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
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. Szklarczyk, A. Franceschini, S. Wyder et al., “STRING v10: protein-protein interaction networks, integrated over the tree of life,” Nucleic Acids Research, vol. 43, no. D1, pp. D447–D452, 2015.
View at: Publisher Site | Google Scholar
J. Chen, X. Zhao, Y. Yuan, and J. J. Jing, “The expression patterns and the diagnostic/prognostic roles of PTPN family members in digestive tract cancers,” Cancer Cell International, vol. 20, no. 1, p. 238, 2020.
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
Y.-Z. Chen, M. Duo Zuo, H.-L. Ren, A.-J. Fan, and G. Ying, “Bioinformatics analysis of expression and alterations of BARD1 in breast cancer,” Technology in Cancer Research & Treatment, vol. 18, article 153303381989226, 2019.
View at: Publisher Site | Google Scholar
G. Zhang, Q. Wang, M. Yang et al., “OSblca: a web server for investigating prognostic biomarkers of bladder cancer patients,” Frontiers in Oncology, vol. 9, p. 466, 2019.
View at: Publisher Site | Google Scholar
F. Chen, D. S. Chandrashekar, S. Varambally, and C. J. Creighton, “Pan-cancer molecular subtypes revealed by mass-spectrometry-based proteomic characterization of more than 500 human cancers,” Nature Communications, vol. 10, no. 1, p. 5679, 2019.
View at: Publisher Site | Google Scholar
D. S. Chandrashekar, B. Bashel, S. A. Balasubramanya et al., “UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses,” Neoplasia, vol. 19, no. 8, pp. 649–658, 2017.
View at: Google Scholar
B. Ru, C. N. Wong, Y. Tong et al., “TISIDB: an integrated repository portal for tumor-immune system interactions,” Bioinformatics (Oxford, England), vol. 35, no. 20, pp. 4200–4202, 2019.
View at: Publisher Site | Google Scholar
V. Thorsson, D. Gibbs, S. Brown et al., “The immune landscape of cancer,” Immunity, vol. 51, no. 2, pp. 411-412, 2019.
View at: Publisher Site | Google Scholar
T. Li, J. Fan, B. Wang et al., “TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells,” Cancer Research, vol. 77, no. 21, pp. e108–e110, 2017.
View at: Publisher Site | Google Scholar
M. Goldman, B. Craft, M. Hastie et al., “Visualizing and interpreting cancer genomics data via the Xena platform,” Nature Biotechnology, vol. 38, no. 6, pp. 675–678, 2020.
View at: Publisher Site | Google Scholar
J. Chen, A. Wang, J. Ji et al., “An innovative prognostic model based on four genes in Asian patient with gastric cancer,” Cancer Research and Treatment, vol. 53, no. 1, pp. 148–161, 2021.
View at: Publisher Site | Google Scholar
J. Fu, K. Li, W. Zhang et al., “Large-scale public data reuse to model immunotherapy response and resistance,” Genome Medicine, vol. 12, no. 1, p. 21, 2020.
View at: Publisher Site | Google Scholar
C. Yue, Y. Ren, H. Ge et al., “Comprehensive analysis of potential prognostic genes for the construction of a competing endogenous RNA regulatory network in hepatocellular carcinoma,” Oncotargets and Therapy, vol. Volume 12, pp. 561–576, 2019.
View at: Publisher Site | Google Scholar
F. Xu, X. Huang, Y. Li, Y. Chen, and L. Lin, “m⁶A-related lncRNAs are potential biomarkers for predicting prognoses and immune responses in patients with LUAD,” Molecular therapy. Nucleic acids, vol. 24, pp. 780–791, 2021.
View at: Publisher Site | Google Scholar
F. Xu, X. Zhan, X. Zheng et al., “A signature of immune-related gene pairs predicts oncologic outcomes and response to immunotherapy in lung adenocarcinoma,” Genomics, vol. 112, no. 6, pp. 4675–4683, 2020.
View at: Publisher Site | Google Scholar
C.-J. Liu, F.-F. Hu, M. Xia, L. Han, Q. Zhang, and A. Y. Guo, “GSCALite: a web server for gene set cancer analysis,” Bioinformatics, vol. 34, no. 21, pp. 3771-3772, 2018.
View at: Publisher Site | Google Scholar
Z. O. Ashtiani, G. Pourmand, S. A. Salami, M. Ayati, and J. Tavakkoly-Bazzaz, “Dysregulated expression of long intergenic non-coding RNAs (LincRNAs) in urothelial bladder carcinoma,” International Journal of Molecular and Cellular Medicine, vol. 6, no. 4, pp. 212–221, 2017.
View at: Publisher Site | Google Scholar
S. Dutertre, M. Cazales, M. Quaranta et al., “Phosphorylation of CDC25B by Aurora-A at the centrosome contributes to the G2-M transition,” Journal of Cell Science, vol. 117, no. 12, pp. 2523–2531, 2004.
View at: Publisher Site | Google Scholar
T. Ochi, H. Fujiwara, and M. Yasukawa, “Aurora-A kinase: a novel target both for cellular immunotherapy and molecular target therapy against human leukemia,” Expert Opinion on Therapeutic Targets, vol. 13, no. 12, pp. 1399–1410, 2009.
View at: Publisher Site | Google Scholar
H. Zhang, J. Song, J. Dong et al., “Tumor microenvironment analysis identified subtypes associated with the prognosis and the tumor response to immunotherapy in bladder cancer,” Frontiers in Genetics, vol. 12, article 551605, 2021.
View at: Publisher Site | Google Scholar
T. Ikezoe, J. Yang, C. Nishioka et al., “A novel treatment strategy targeting Aurora kinases in acute myelogenous leukemia,” Molecular Cancer Therapeutics, vol. 6, no. 6, pp. 1851–1857, 2007.
View at: Publisher Site | Google Scholar
F. Li, G. Ambrosini, E. Chu et al., “Control of apoptosis and mitotic spindle checkpoint by survivin,” Nature, vol. 396, no. 6711, pp. 580–584, 1998.
View at: Publisher Site | Google Scholar
E. Tanaka, Y. Hashimoto, T. Ito et al., “The clinical significance of Aurora-A/STK15/BTAK expression in human esophageal squamous cell carcinoma,” Clinical cancer research : an official journal of the American Association for Cancer Research, vol. 11, no. 5, pp. 1827–1834, 2005.
View at: Publisher Site | Google Scholar
J. Liang, W. Zhao, P. Tong et al., “Comprehensive molecular characterization of inhibitors of apoptosis proteins (IAPs) for therapeutic targeting in cancer,” BMC Medical Genomics, vol. 13, no. 1, p. 7, 2020.
View at: Publisher Site | Google Scholar
A. Lapitz Dambolenea, A. Arbelaiz Cossio, C. J. O’Rourke et al., “Patients with cholangiocarcinoma present specific RNA profiles in serum and urine extracellular vesicles mirroring the tumor expression: novel liquid biopsy biomarkers for disease diagnosis,” Cell, vol. 9, no. 3, p. 721, 2020.
View at: Publisher Site | Google Scholar
Y. Luo, W. Huang, H. Zhang, and G. Liu, “Prognostic significance of CD117 expression and <i>TP53</i> missense mutations in triple-negative breast cancer,” Oncology Letters, vol. 15, no. 5, pp. 6161–6170, 2018.
View at: Publisher Site | Google Scholar
C. Zhang, L. Shen, F. Qi, J. Wang, and J. Luo, “Multi-omics analysis of tumor mutation burden combined with immune infiltrates in bladder urothelial carcinoma,” Journal of Cellular Physiology, vol. 235, no. 4, pp. 3849–3863, 2020.
View at: Publisher Site | Google Scholar
W. H. Fridman, F. Pagès, C. Sautès-Fridman, and J. Galon, “The immune contexture in human tumours: impact on clinical outcome,” Nature Reviews. Cancer, vol. 12, no. 4, pp. 298–306, 2012.
View at: Publisher Site | Google Scholar
K. Kuzushima, M. Yasukawa, T. Ochi et al., “Aurora-A kinase: a novel target of cellular immunotherapy for leukemia,” Blood, vol. 113, no. 1, pp. 66–74, 2009.
View at: Publisher Site | Google Scholar
J. Shen, R. Wang, Y. Chen et al., “An immune-related signature predicted survival in patients with kidney papillary cell carcinoma,” Frontiers in Oncology, vol. 11, article 670047, 2021.
View at: Publisher Site | Google Scholar
B. Zhao, A. Song, R. Haque et al., “Cooperation between molecular targets of costimulation in promoting T cell persistence and tumor regression,” The Journal of Immunology, vol. 182, no. 11, pp. 6744–6752, 2009.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Ya Guo 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

450

Downloads

725

Citations