- Open Access
Construction and validation of a gene signature related to bladder urothelial carcinoma based on immune gene analysis
BMC Cancer volume 22, Article number: 926 (2022)
This study developed a gene signature associated with a malignant and common tumor of the urinary system, the Bladder Urothelial Carcinoma (BLCA).
The Cancer Genome Atlas (TCGA) database was searched to obtain 414 BLCA samples and the expression spectra of 19 normal samples. Single-sample Gene Set Enrichment Analysis (ssGSEA) was conducted to determine the enrichment levels in the BLCA samples of the 29 immune genes. Unsupervised hierarchical clustering, gene set enrichment analysis (GSEA), single-factor Cox analysis, least absolute shrinkage and selection operator (LASSO) regression models, and GEO queues were used to determine the BLCA immune gene subtype, analyze the biological pathway differences between immune gene subtypes, determine the characteristic genes of BLCA associated with prognosis, identify the BLCA-related genes, and verify the gene signature, respectively.
We identified two immune gene subtypes (immunity_L and immunity_H). The latter was significantly related to receptors, JAK STAT signaling pathways, leukocyte interleukin 6 generation, and cell membrane signal receptor complexes. Four characteristic genes (RBP1, OAS1, LRP1, and AGER) were identified and constituted the gene signature. Significant survival advantages, higher mutation frequency, and superior immunotherapy were observed in the low-risk group patients. The gene signature had good predictive ability. The results of the validation group were consistent with TCGA queue results.
We constructed a 4-gene signature that helps monitor BLCA occurrence and prognosis, providing an important basis for developing personalized BLCA immunotherapy.
Bladder cancer is amongst the ten most common forms of cancers worldwide. In 2020, approximately 213,000 deaths from approximately 573,000 new cases of bladder cancer were estimated to have occurred . The incidence of bladder cancer in men is approximately three to four times that in women, but most women are already in an advanced stage of the disease when they are diagnosed with bladder cancer, thus leading to a worse prognosis. There are several subtypes of bladder cancer, including bladder urothelial carcinoma (BLCA), squamous cell carcinoma, adenocarcinoma, small cell carcinoma, and sarcoma, among which the first is the main one . There has been a rising trend in the incidence and mortality rates of BLCA over the years. With the aging of the population and the increase in environmental pollution, the public health burden of BLCA will also increase [3, 4]. The usual procedures used for the treatment of BLCA are transurethral resection, radical cystectomy, radiotherapy, and chemotherapy. However, BLCA is aggressive and highly recurrent and has a poor prognosis. Therefore, the diagnosis of BLCA, its treatment, and the five-year survival rates have remained basically unchanged [2, 5,6,7]. Traditional therapies can no longer meet the current status of BLCA treatment. In recent years, with the development of tumor immunology and the introduction of immunosuppressive agents in cancer treatment, the treatment of bladder cancer is expected to progress . A large number of studies has analyzed the immune characteristics of BLCA patients from the perspective of immune infiltration [9,10,11]. However, there are still few studies on exploring new potential prognostic and immunotherapeutic biomarkers of BLCA with respect to tumor immune interaction.
The number of studies on the role played by the tumor immune microenvironment in tumor development and tumorigenesis has been increasing with the rapid advancement of the next-generation sequencing technology [12, 13]. The tumor microenvironment has a huge impact on the occurrence and development of BLCA and also helps to predict the response of patients with BLCA to immunotherapy [14, 15]. Immunotherapies such as immune checkpoint blockade, cytokine therapy, and cell therapy for the treatment of cancer have evolved in recent years [16, 17]. Immunotherapy is a potential neoadjuvant therapy for the treatment of BLCA. However, only a few patients benefit from this treatment [18,19,20]. The increase in immune cell infiltration in BLCA is an effective indicator of the response to immunotherapy [21, 22]. Significant correlation between immunotherapy and Human Leukocyte Antigen (HLA) and PD-L1 expression was observed [23,24,25]. The accumulation of somatic mutations can promote the occurrence and development of tumors and is conducive to the expression of neoantigens . Previous cumulative studies have also shown that the tumor mutational burden (TMB) can predict the prognosis of patients with tumors [27, 28], and recent studies have begun to explore the correlation between the TMB and immunotherapy response [29, 30]. Therefore, it is necessary to understand the key factors affecting immunotherapy in BLCA.
In this study, the abundance level of immune cells in BLCA were quantified applying single-sample Gene Set Enrichment Analysis (ssGSEA). We further determined the immune subtype of BLCA according to the ssGSEA score. The differential genes between the BLCA immune subtypes were determined, and the intersection with the immune gene set was used to obtain BLCA-related immune genes related to the prognosis. Next, a new gene signature was constructed based on these immune genes related to the prognosis of BLCA to calculate the survival risk of the patient. In addition, according to the clinicopathological parameters and the signature, a nomogram to predict the 1-, 3-, and 5-year overall survival rates was constructed. Our research provides a new potential gene signature and immunotherapy classification method for the treatment of BLCA to promote the development of a personalized and precise treatment approach.
Materials and methods
Downloading and processing the research data
We downloaded the research data of this study (TCGA-BLCA) from The Cancer Genome Atlas (TCGA) database, including the transcriptome data and clinical information associated with 414 BLCA samples and 19 normal samples. We used "bladder cancer” (title) and "series" (entry type) and "Homo sapiens" (organism) as keywords in the GEO database ( http://www.ncbi.nlm.nih.gov/geo/) to obtain two datasets (GSE32894 and GSE31684) and series matrix files of related microarray platforms. Then, the two datasets were merged, and the expression data of that dataset were batch standardized using the "sva" package in R, so as to obtain a normalized gene expression matrix file for external validation analysis. The TCGA database provided the mutation data for the BLCA samples. In addition, the 2,483 immune genes were obtained from the Immunology Database and Analysis Portal (ImmPort (https://www.immport.org/). The imvigor210 immunotherapy dataset was taken from a study published by Sanjeev Mariathasan et al. .
Cluster analysis of BLCA samples
To obtain the ssGSEA score subsequent to quantifying the level of enrichment of the immune cells in every BLCA sample, single-sample gene enrichment analysis (ssGSEA) was applied in the TCGA-BLCA cohort. Next, to determine the immune subtype of BLCA, an unsupervised cluster analysis of BLCA afflicted patients was conducted based on the ssGSEA score.
Analysis of the correlation between BLCA immune genome subtypes and molecular characteristics
ESTIMATE is an algorithm used to estimate the immune scores, matrix scores, and ESTIMATE scores in tumor microenvironments. Using the ESTIMATE algorithm, the tumor purity, matrix content, and level of immune cell infiltration of each BLCA sample were calculated. To establish the relationship between the different BLCA immune infiltration and immune subtypes, the stromal, immune, and ESTIMATE scores of the various immune subtypes were compared. The differences in HLA genes in the BLCA immune subtypes were assessed, and their levels of expression in each subgroup determined. In addition, for comparing the ratio of immune cell subgroups between the BLCA immune subtypes, and for evaluating tumor-infiltrating immune cells in the BLCA samples, the CIBERSORT deconvolution algorithm was used . The CIBERSORT platform (https://cibersortx.stanford.edu/) provided the gene expression characteristic matrix of the 22 tumor-infiltrating immune cells.
Functional analysis between different immune subgroups and risk groups
Gene enrichment analysis (GSEA) was used to clarify differences in molecular mechanisms between different immune gene subtypes and between different risk groups. Gene Ontology (GO) analysis was based on c5.go.v7.symbols.gmt and consisted of three parts: biological process (BP), molecular function (MF), and cellular component (CC). Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis was performed using c2.cp.kegg.v7.symbols.gmt [33,34,35]. Statistically significant values were considered to have an FDR < 0.25 and a P < 0.05.
Identification of differentially expressed genes (DEGs) between different immune subgroups
The ‘limma’ package in R was used for screening the DEGs amongst the various immune subtypes, with a False Discovery Rate (FDR) < 0.05 criterion.
Identification of the significance of immune genes in BLCA prognosis
We used the intersection between DEGs and immune gene sets to identify immune genes related to BLCA. Using a P < 0.001 as the screening criterion, Univariate Cox Regression Analysis was applied.
Construction and verification of the BLCA gene signature
We used the least absolute shrinkage and selection operator (LASSO) to analyze and construct the BLCA-related gene signature. The LASSO algorithm reduces the number of variables by constructing a penalty function to effectively avoid overfitting. The penalty parameters of the model are based on the minimum partial likelihood deviation λ. This value was determined by tenfold cross-validation. The characteristic genes and corresponding coefficients contained in the risk signature were determined by the λ-value. The formula of the BLCA gene signature was as follows: BLCA risk score = gene1 × coef1 + gene2 × coef2 + gene3 × coef3 + …… + gene(n) × coef(n), with the regression coefficient being represented as ‘coef’. The BLCA samples from the TCGA cohort were segregated into high-risk and low-risk categories based on the median values of the risk scores. To evaluate the BLCA gene signature value and the prognostic ability, multivariate and univariate Cox regression analyses, survival curve, and Receiver Operating Characteristic (ROC) curve were applied. In addition, we used a test cohort and GEO queue as verification queues, and verified the stability and reliability of the gene signature using the same method.
Analyzing the correlation between the built gene signature, the immune microenvironment, clinical characteristics, and the immune therapy response
A heatmap was used to analyze the correlation between the gene signature and the immune cell genes. We also explored the correlation between BLCA-related genes and clinical characteristics, including age, sex, grading, and stage. The TCGA database provided the clinical feature information. The immunophenoscore (IPS) is a good predictor of CTLA-4 and PD-1 reactivity, and immune checkpoint inhibitor immunoformidal scoring files were obtained from the Cancer Immunogenic Group database [36,37,38]. Based on the IPS files, we analyzed and predicted the differences in the responses in the group to immunotherapy using CTLA-4 and PD-1 blockers. Moreover, we also predicted the survival difference between high-risk and low-risk groups based on imvigor210.
Development and verification of the nomogram
Based on the risk score and the independent clinicopathological characteristics, a nomogram was developed for evaluating the survival probability for the 1-, 3-, and 5-year periods of each BLCA patient. A calibration chart and a time-dependent ROC curve were used to verify the nomogram performance.
R v4.0.3 (R Foundation for Statistical Computing, Vienna, Austria) was used to conduct all statistical analyses. The “SurvivalROC” package was used to calculate the area below the ROC curve (AUC). An AUC ≥ 0.70 and an AUC ≥ 0.6 were defined as significant and acceptable predictors, respectively. The Kaplan–Meier method and order-to-peer test were used for survival analysis. Statistical significance was considered in case of P-value < 0.05.
Identification of the BLCA immune gene subgroup
In this study, we used ssGSEA to quantify the concentration levels of 29 immune-related genes in BLCA. We identified 414 BLCA tumor samples as belonging to two immune gene subtypes (immunity_L and immunity_H) based on unsupervised clustering and on the ssGSEA scores (Fig. 1A). Figure 1B shows that the ssGSEA score can be used to distinguish immunity_L from immunity_H BLCA patients. The heat map shows the expression levels of 29 immune-related differential genes between the two different subgroups of immune genes (Fig. 1C).
Correlation analysis of the BLCA immunogene subgroups, the tumor microenvironment, and the HLA gene
We applied the ESTIMATE algorithm to determine the ESTIMATE, stromal, and immune scores of the BLCA tumor microenvironment. Figure 2A shows the immunity scores of the immunity_H group, which had higher stromal and ESTIMATE scores than the immunity_L group. In contrast, a lower level of tumor purity was observed for the immunity_H than the immunity_L group. In addition, there were significant differences in the immune, stromal, and ESTIMATE scores of the two BCLA immune gene subtypes (P < 0.001; Fig. 2B). Furthermore, we systematically estimated the proportions of 22 tumor immune cells in BLCA using the CIBERSORT algorithm. Compared to in the immunity_H group, higher levels of activation of the monocytes and dendritic cells were observed in the immunity L group (P < 0.05; Fig. 2C). Conversely, higher levels of expression of the macrophages M0 and M2 were observed (P < 0.05) in the immunity H group than in the immunity L group.
We also analyzed differences in the expression of HLA genes between the two BLCA immune gene subgroups (Fig. 2D). The HLA gene expression level between the immunity_L and immunity_H groups was found to significantly differ, whereas the level of expression of HLA genes was found to be higher in the immunity_H group than in the immunity_L group (P < 0.001).
Gene function enrichment analysis and identification of differentially expressed immune genes between the BLCA immune gene subtypes
To understand the biological pathway differences between the two different immune gene subtypes of BLCA, we performed GSEA using c2.cp.kegg.v7.symbols.gmt and c5.go.v7.symbols.gmt. The KEGG enrichment analysis results showed that cytokine receptor interactions, B-cell receptor signaling pathways, apoptosis, Toll-like receptor signaling pathways, and JAK STAT signaling pathways were significantly active in the immunity_H group (Fig. 3A). The results of GO enrichment analysis showed that in BP, pathways such as leukocyte interleukin 6 production, cell polarization, and positive regulation of cell adhesion were significantly active in the immunity_H group. In the case of MF, pathways such as cytokine binding, cytokine receptor binding, and extracellular matrix binding were significantly active in the immunity_H group. Regarding the CC, the extracellular matrix containing collagen, the membrane signal receptor complex, and the protein complex involved in cell adhesion were significantly related to the immunity_H group (Table 1).
To screen for the DEGs associated with immunity in the two BLCA immunogen subtypes, we first identified gene expression differences between the immunity_L and immunity_H groups. 3,340 DEGs were identified in total. The expression levels of these genes are presented in a heat map (Fig. 3B), while a volcanic map shows the differences in the expression of these genes between the two immune gene subtypes (Fig. 3C). Next, we used the previously obtained differences to express the intersection of genes and immune genes. Finally, 478 BLCA immune-related DEGs were obtained (Fig. 3D).
Construction of the immunity gene signature related to BLCA
We performed Univariate Cox analysis of the immune genes based on their differential expression to identify the immune genes associated with BLCA prognosis, using P < 0.01 as the threshold. A total of 21 immune genes associated with the prognosis of patients with BLCA were obtained (Fig. 4A). The LASSO regression model was used to screen the BLCA signature genes for the development of a new BLCA gene signature (Fig. 4B, C). Finally, four BLCA signature genes were used to develop the gene signature. The formula used for developing the BLCA gene signature is as follows: risk score = 0.1601 × expression of RBP1 + (-0.2471) × expression of OAS1 + (0.2570) × expression of LRP1 + (-0.4403) × expression of AGER.
The TCGA cohort was randomly segregated into test and training groups. Based on the BLCA gene signature, the patients with BLCA in the training group were further segregated into low-risk and high-risk groups (Fig. 5A). An increase in the risk score was associated with a decrease in the survival time of the BLCA patients (Fig. 5B). The overall survival of the patients in the low-risk group was observed to be longer than that of patients in the high-risk group, based on the Kaplan–Meier survival curve results (P < 0.001, Fig. 5C). To evaluate the predictive power of the BLCA gene signature, the total survival projections for the 1-, 3-, and 5-year periods based on the gene signature were 0.678, 0.697, and 0.704, respectively (Fig. 5D). This indicated that the BLCA gene signature has good predictive power. In the patients with BLCA, risk score, TNM stage, stage, and age were associated with poor survival in the Univariate analysis results. The potential use of the risk score of the gene signature as independent prognostic factor is illustrated by the multivariate analysis results (HR: 1.76; 95% CI: 1.43–2.16; P < 0.001; Tables 2 and 3).
Functional verification of the BLCA-related gene signature
To verify the reliability of the application of the BLCA gene signature to different populations, we calculated the risk scores of patients with BLCA in the test group and GEO queue using the same formula. The test group and GEO queue were used to verify that the results of the queue are consistent with the previous results. Depending on the risk scores of the gene signature, the patients with BLCA in the validation queue were segregated into low-risk and high-risk categories (Fig. 5E, I). An increase in the risk score was associated with lower survival outcomes in patients (Fig. 5F, J). Moreover, compared with the survival results of the high-risk group, the results of the low-risk group were superior, as seen from the results of the Kaplan–Meier survival curve (test group: P = 0.001; GEO queue: P = 0.032; Fig. 5G, K). In the test group, the total survival projections for one, three, and five years based on the AUC results were 0.658, 0.686, and 0.714, respectively (Fig. 5H). In the GEO cohort, the total survival projections for one, three, and five years based on the AUC results were 0.574, 0.639, and 0.638, respectively (Fig. 5L). The potential to accurately predict the prognosis of the patients with BLCA and the applicability of the gene signature in different populations was confirmed by these results. The relationship between the risk score and immunophenotype was analyzed and visualized using an alluvial map. Lower risk scores were observed in the immunity L group(Fig. 5M).
Correlation between gene signature and immunophenotyping, immune cells, and somatic mutation
Compared to that in the immunity_L group, the expression of RBP1 and LRP1 in the immunity_H group was higher than that in the immunity_L group. Conversely, the expression of AGER in the immunity H group was lower than that in the immunity L group (Fig. 6A). Correlation between the four characteristic genes and the immune cells significant for the development of the BLCA. The results showed that AGER was significantly positively correlated with regulatory T cells, plasma cells, and B cells. LRP1 was significantly positively correlated with resting mast cells and macrophages M2, while it was significantly negatively correlated with follicular helper T cells, dendritic cells, and CD8T cells. OAS1 was significantly positively correlated with CD8T cells, memory CD4T cells, and dendritic cells, whereas RBP1 was significantly negatively correlated with plasma and B cells (Fig. 6B).
Studies have shown a correlation between mutations in somatic cells in tumor tissues and immune therapy response. We analyzed the distribution patterns of the tumor mutation load (TMB) in the low- and high-risk groups. The frequency of mutations in the low-risk group was higher than that in the high-risk group (Fig. 6C). Compared with the patients with a low mutation frequency, the survival results of the patients with a high mutation frequency were found to be markedly better, as revealed by the survival analysis results (Fig. 6D). Moreover, in both the low and high mutation groups, the patients in the low-risk group indicated a clear advantage of survival (Fig. 6E).
Predictive value of the gene signature for immunotherapy effects
An analysis of the expression of the Tumor Immune Dysfunction and Exclusion (TIDE) in the high-risk and the low-risk groups revealed higher incidence of TIDE in the low-risk group than in the high-risk group (P < 0.05; Fig. 7A). In the CTLA4 treatment group, patients in the low-risk group exhibited better outcomes (P = 7.7e − 07; Fig. 7B). More importantly, we used the imvigor210 dataset to evaluate the ability of our model to predict the effect of immunotherapy based on four characteristic genes. Compared to the patients in the high-risk group, the patients in the low-risk group were associated with better survival outcomes in the Kaplan Meier analysis (Fig. 7C).
Nomogram prediction model construction and evaluation
We built a line chart based on clinical characteristics, such as the BLCA-related genes, age, sex, and stage, to predict the total survival rates of patients with BLCA over one, three, and five years (Fig. 8A). The total survival projections for one, three, and five years based on time-dependent ROCs were 0.768, 0.758, and 0.778, respectively (Fig. 8B). We also used calibration charts to evaluate the predictive effect of this line chart on the total lifetime of patients with BLCA (Fig. 8C). These results suggest that the genes that combine clinical characteristics have potential clinical benefits.
Functional enrichment analysis between high- and low-risk groups
The GSEA results revealed differences in biological pathways between the high- and low-risk groups. The high-risk group was mainly enriched in cell adhesion molecules, cytokine receptor interaction, ecm receptor interaction, focal adhesion, and regulation of actin cytoskeleton (Fig. 9A). Conversely, the low-risk group was mainly enriched in metabolism of cytochrome P450 to heterologous substances, oxidative phosphorylation, mutual transformation of pentose and glucuronic acid, porphyrin and chlorophyll metabolism, and ribosome (Fig. 9B).
BLCA is a common malignant genitourinary tract tumor . However, progress in bladder cancer treatment has been limited. Although surgery is still the preferred treatment for BLCA, patients have a high rate of postoperative recurrence . Hence, new strategies for treatment need to be developed in which the identification of BLCA biomarkers is imperative. In the present study, the relative abundance of the 29 immune genes in the BLCA samples was calculated using ssGSEA. Two different BLCA immune gene subtypes (immunity_L and immunity_H) were determined using unsupervised clustering methods based on ssGSEA scores. We screened differentially expressed immune genes between the two immune gene subtypes and identified the immune genes associated with prognosis using univariate cox analysis. Finally, we built a gene signature related with BLCA based on the immune genes associated with prognosis. The BLCA-related gene signature could divide patients into low- and high-risk groups, and the patients in the low-risk group had significant survival advantages and showed significant immunotherapy effects. In summary, we aimed to provide a new gene signature to promote the personalized treatment of BLCA.
We observed that the expression levels of monocytes and dendritic cells activated in the immunity_L group were higher than those in the immunity_H group, whereas the expression levels of macrophages M0 and M2 in the immunity_L group were lower than those in the immunity_H group. Besides being mononuclear phagocytes, the monocytes were found to be significant regulators of cancer development and progression. Contrasting roles were played by the different subsets in the prevention of cancer cell metastasis and diffusion, along with the promotion of tumor growth [41, 42]. Recent research by Wang et al. revealed that intermediate monocytes can be induced by IFN-γ to inhibit cancer metastasis by promoting natural killer cell activation by FOXO1 and interleukin 27 . Significant roles are played by the macrophages in inflammation and innate immunity . M2 macrophages are a subtype of macrophages that can inhibit inflammatory response and promote angiogenesis and tissue remodeling and repair . Previous studies have shown that bladder cancer cells can stimulate BMP4 to induce macrophages and polarize the M2 phenotype, thereby promoting bladder cancer development . The levels of RBP1 and LRP1 were significantly higher in the immunity H group than in the immunity_L group.
We analyzed the biological pathway differences between the two immune gene subgroups of BLCA using GSEA. Cytokine receptor interactions, receptor signaling pathways, JAK STAT signaling pathways, leukocyte interleukin 6 generation, and mass membrane signal receptor complexes were more active in the Immunity_H group. Having a crucial role in the immune system, Toll-like Receptors (TLRs) are members of the pattern recognition receptor family. TLRs can mediate the identification of the pathogen-related molecular patterns derived from pathogens and damage-related molecular patterns from damaged cells, further driving host innate and adaptive immunity activation [47,48,49]. Chronic inflammation can promote tumor development. TLR is an effective activator of the inflammatory response . TLR expression in cancer cells can trigger an NF-B signaling cascade, the production of inflammatory cytokines, and anti-apoptosis proteins, thereby promoting tumor growth . In several cancers, TLR4 has been indicated to play a significant role. TLR4 activation increases the expression of VEGF and TGF-β1 in prostate cancer cells, which promotes tumor development . Overexpression of TLR4 has been associated with poor prognosis in breast and colon cancers [53, 54]. In addition, TLR4 has been found to have antitumor activity in skin cancer . Previous studies have found that TLR4 is less expressed in non-myo-inundated bladder cancer . TLR4 may have a significant impact on BLCA development. The JAK/STAT signaling pathway is an important component of functional reactions. Three major proteins, a cell surface receptor, JAK, and STAT, are involved in JAK/STAT signaling . An abnormal activation of the JAK/STAT signaling pathway may be a key factor promoting tumor development . STAT3 promotes cell cycle processes by promoting the activation of cell cycle protein-dependent kinases. STAT5 protects cells from apoptosis by activating the transcription of Bcl-x to produce the anti-apoptosis protein Bcl-XL. The inhibition of STAT activity reduces cell proliferation and increases apoptosis . Therefore, antagonist JAK/STAT signaling may help to inhibit the development of BLCA.
In the present study, we constructed a gene signature containing four BLCA signature genes (RBP1, OAS1, LRP1, and AGER). Retinol-binding protein type 1 (RBP1) is a cytosolic carrier responsible for regulating retinol homeostasis in human tissues [60, 61]. Previous studies have observed abnormal expression of RBP1 in a variety of human cancers [62,63,64,65,66]. RBP1 is responsible for regulating intracellular retinoic acid homeostasis, which is related to morphogenesis, cell proliferation, and differentiation. Therefore, the deletion of RBP1 may lead to tumorigenesis in bladder cancer. The abnormal CpG hypermethylation of the RBP1 gene promoter may be related to bladder cancer development . DNA damage is an important step in cancer development . Poly (ADP-Ribose) is a quick synthesis at the DNA damage site that facilitates repair. The high expression of OAS1 promotes the ability of cancer cells to survive DNA damage by reducing PAR synthesis, thus preventing cell death . Low-density lipoprotein receptor-related protein (LRP) is a transmembrane receptor that plays a role in cell swallowing, cell signaling, and other cell protein transport . LRP1 can act as a receptor for medium-term growth factors, signaling through ERK1/2 to induce the expression of internalized metalloproteinase, and promote the survival of micro-metastasis, thereby promoting cancer progression [70,71,72,73]. AGER is a member of the superfamily of cell surface receptors for immunoglobulins. A study by Wang et al. has shown that AGER overexpression, which is associated with age, reduces the proliferation, invasion, and migration of lung cancer cell lines and increases apoptosis . However, the mechanisms underlying the actions of these characteristic genes in BLCA require further studies.
The survival of patients in the TCGA-BLCA queue was measured to determine the prognostic value of the gene signature. The gene signature could divide the BLCA patients into high- and low-risk groups based on the risk scores, and patients in the low-risk group had significant survival advantages. The immunity_L group was more associated with the low-risk group. In addition, we monitored the progress of BLCA through the risk curve generated based on the risk score. The ROC curve showed that the gene signature has good predictive ability. All results were validated in the test group and GEO queue. This suggests that our gene signature can be applied to other individuals in queues. The risk score, being an independent predictor of the total survival of the BLCA patients, was confirmed further by multivariate analysis. Moreover, our genetic signature risk scores were significantly associated with immune response predictors, such as immune scores, and TIDE. Good response to immunotherapy was displayed by the low-risk group patients. This gene signature facilitates the personalized treatment for bladder cancer. As a result, our gene signature is expected to be a valuable tool for evaluating the prognosis of patients with BLCA.
The correlation between the genes that build immune cells and the gene signature was also analyzed. The four characteristic genes were significantly associated with immune cells, such as regulatory T cells, plasma cells, B cells, resting mast cells, macrophages M2, follicle-assisted T cells, and dendritic cells. Previous studies have shown that tumor-related macrophages can promote and may even be necessary for angiogenesis in bladder cancer [75, 76]. Macrophages M2 are phenotype phenocytes and may promote BLCA development by inducing immunosuppression, angiogenesis, and cell invasion . Dendritic cells are considered the most effective anti-cancer immunosuppressors . Previous studies have found the presence of immature and minimally activated dendritic cells in the urine and bladder of BLCA patients . We believe that the reduction or inhibition of dendritic cells in the tumor microenvironment may be a key factor promoting BLCA development. Regulatory T cells (Tregs) play an important role in cancer immune escape . Poor prognosis in a variety of tumor patients was revealed by the high levels of infusion of Treg in the associated tumor tissues . We believe that the removal of Treg cells from the BLCA tumor microenvironment has the potential to stimulate and enhance antitumor immune response in patients with BLCA.
The present study had some limitations. First, we only analyzed data from public databases. Second, we need to examine larger queues of BLCA patients receiving immunotherapy in the future to verify our results.
In the present study, we developed and validated a gene signature containing four genes. The gene signature could be used as an independent prognostic factor for BLCA patients. The gene signature may help monitor the occurrence and prognosis of BLCA. Moreover, the clinical response of patients with BLCA to immunotherapy can be predicted using the gene signature. Our findings thus provide new insights into personalized BLCA immunotherapy.
Availability of data and materials
The BLCA data used for analysis in this study were obtained from the TCGA database (https://tcga-data.nci.nih.gov/tcga/) and the GEO database (http://www.ncbi.nlm.nih.gov/geo/).
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, Bray F. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin. 2021;71(3):209–49.
Berdik C. Unlocking bladder cancer. Nature. 2017;551(7679):S34–5.
Yang Y, Cheng Z, Jia X, Shi N, Xia Z, Zhang W, Shi X. Mortality trends of bladder cancer in China from 1991 to 2015: an age-period-cohort analysis. Cancer Manage Res. 2019;11:3043–51.
Ebrahimi H, Amini E, Pishgar F, Moghaddam SS, Nabavizadeh B, Rostamabadi Y, Aminorroaya A, Fitzmaurice C, Farzadfar F, Nowroozi MR, et al. Global, Regional and National Burden of Bladder Cancer, 1990 to 2016: Results from the GBD Study. 2016. J Urol. 2019;201(5):893–901.
Kirkali Z, Chan T, Manoharan M, Algaba F, Busch C, Cheng L, Kiemeney L, Kriegmair M, Montironi R, Murphy WM. Bladder cancer: epidemiology, staging and grading, and diagnosis. Urology. 2005;66(6-suppl-S1):4–34.
Humphrey PA, Moch H, Cubilla AL, Ulbright TM, Reuter VE. The 2016 WHO Classification of Tumours of the Urinary System and Male Genital Organs-Part B: Prostate and Bladder Tumours. Eur Urol. 2016;70(1):106–19.
Goebell PJ, Legal W, Weiss C, Fietkau R, Wullich B, Krause S. [Multimodal therapy for bladder sparing with high grade bladder tumors]. Der Urologe Ausg A. 2008;47(7):838, 840–32, 844–35.
Hogendoorn PC, Athanasou, Bielack, Alava D, Tos DAP. Bladder cancer: ESMO Practice Guidelines for diagnosis, treatment and follow-up. Ann Oncol. 2014;40(5):584–91.
Cao R, Yuan L, Ma B, Wang G, Tian Y. Immune-related long non-coding RNA signature identified prognosis and immunotherapeutic efficiency in bladder cancer (BLCA). Cancer Cell Int. 2020;20:276.
Wang Z, Tu L, Chen M, Tong S. Identification of a tumor microenvironment-related seven-gene signature for predicting prognosis in bladder cancer. BMC Cancer. 2021;21(1):692.
Zhang LH, Li LQ, Zhan YH, Zhu ZW, Zhang XP. Identification of an IRGP Signature to Predict Prognosis and Immunotherapeutic Efficiency in Bladder Cancer. Front Mol Biosci. 2021;8:607090.
Smith CC, Beckermann KE, Bortone DS, Cubas A, Bixby LM, Lee SJ, Panda A, Ganesan S, Bhanot G, Wallen EM. Endogenous retroviral signatures predict immunotherapy response in clear cell renal cell carcinoma. J Clin Invest. 2018;128(11):4804–20.
Lu X, Horner JW, Paul E, Shang X, Troncoso P, Deng P, Jiang S, Chang Q, Spring DJ, Sharma P, et al. Effective combinatorial immunotherapy for castration-resistant prostate cancer. Nature. 2017;543(7647):728–32.
Li F, Guo H, Wang Y, Liu B, Zhou H. Profiles of tumor-infiltrating immune cells and prognostic genes associated with the microenvironment of bladder cancer. Int Immunopharmacol. 2020;85:106641.
Carosella ED, Ploussard G, Lemaoult J, Desgrandchamps F. A Systematic Review of Immunotherapy in Urologic Cancer: Evolving Roles for Targeting of CTLA-4, PD-1/PD-L1, and HLA-G. Eur Urol. 2015;68(2):267–79.
Del Paggio JC. Immunotherapy: Cancer immunotherapy and the value of cure. Nat Rev Clin Oncol. 2018;15(5):268–70.
Christofi T, Baritaki S, Falzone L, Libra M, Zaravinos A. Current Perspectives in Cancer Immunotherapy. Cancers. 2019;11(10):1472.
Morales A, Eidinger D, Bruce AW. Intracavitary Bacillus Calmette-Guerin in the treatment of superficial bladder tumors. J Urol. 1976;167(2):891–4.
Necchi A, Anichini A, Raggi D, Briganti A, Massa S, Lucianò R, Colecchia M, Giannatempo P, Mortarini R, Bianchi M. Pembrolizumab as Neoadjuvant Therapy Before Radical Cystectomy in Patients With Muscle-Invasive Urothelial Bladder Carcinoma (PURE-01): An Open-Label, Single-Arm, Phase II Study. J Clin Oncol. 2018;36(34):JCO1801148.
Braun DA, Burke KP, Allen EV. Genomic Approaches to Understanding Response and Resistance to Immunotherapy. Clin Cancer Res. 2016;22(23):5642–50.
Cheng W, Fu D, Xu F, Zhang Z. Unwrapping the genomic characteristics of urothelial bladder cancer and successes with immune checkpoint blockade therapy. Oncogenesis. 2018;7(1):2.
Sharma P, Allison JP. The future of immune checkpoint therapy. Science. 2015;348(6230):56.
Kim HS, Seo HK. Immune checkpoint inhibitors for urothelial carcinoma. Invest Clin Urol. 2018;59(5):285–96.
Patel SP, Kurzrock R. PD-L1 Expression as a Predictive Biomarker in Cancer Immunotherapy. Mol Cancer Ther. 2015;14(4):847–56.
Zhang B, Yu W, Feng X, Zhao Z, Fan Y, Meng Y, Hu S, Cui Y, He Q, Zhang H. Prognostic significance of PD-L1 expression on tumor cells and tumor-infiltrating mononuclear cells in upper tract urothelial carcinoma. Med Oncol. 2017;34(5):94.
Schreiber RD, Gubin MM, Mardis ER, Artyomov. Tumor neoantigens: building a framework for personalized cancer immunotherapy. J Clin Invest. 2015;125(9):3413–21.
Chen H, Chong W, Wu Q, Yao Y, Mao M, Wang X. Association of LRP1B Mutation With Tumor Mutation Burden and Outcomes in Melanoma and Non-small Cell Lung Cancer Patients Treated With Immune Check-Point Blockades. Front Immunol. 2019;10:1113.
Chen Y, Liu Q, Chen Z, Wang Y, Zhang H: PD-L1 expression and tumor mutational burden status for prediction of response to chemotherapy and targeted therapy in non-small cell lung cancer. J Exp Clin Cancer Res. 2019, 38(1).
Goodman AM, Kato S, Bazhenova L, Patel SP, Frampton GM, Miller V, Stephens PJ, Daniels GA, Kurzrock R. Tumor Mutational Burden as an Independent Predictor of Response to Immunotherapy in Diverse Cancers. Mol Cancer Ther. 2017;16(11):2598–608.
Rizvi NA, Hellmann MD, Snyder A, Kvistborg P, Makarov V, Havel JJ, Lee W, Yuan J, Wong P, Ho TS, et al. Cancer immunology. Mutational landscape determines sensitivity to PD-1 blockade in non-small cell lung cancer. Science. 2015;348(6230):124–8.
Mariathasan S, Turley SJ, Nickles D, Castiglioni A, Yuen K, Wang Y, Kadel EE III, Koeppen H, Astarita JL, Cubas R, et al. TGFβ attenuates tumour response to PD-L1 blockade by contributing to exclusion of T cells. Nature. 2018;554(7693):544–8.
Chen B, Khodadoust MS, Liu CL, Newman AM, Alizadeh AA. Profiling Tumor Infiltrating Immune Cells with CIBERSORT. Methods Mol Biol (Clifton, NJ). 2018;1711:243–59.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Prot Sci. 2019;28(11):1947–51.
Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(D1):D545-d551.
Charoentong P, Finotello F, Angelova M, Mayer C, Efremova M, Rieder D, Hackl H, Trajanoski Z. Pan-cancer Immunogenomic Analyses Reveal Genotype-Immunophenotype Relationships and Predictors of Response to Checkpoint Blockade. Cell Rep. 2017;18(1):248–62.
Eliezer M, Allen V, Diana, Miao, Bastian, Schilling, Sachet A, Shukla, Christian, Blank: Genomic correlates of response to CTLA-4 blockade in metastatic melanoma. Science (New York, NY) 2015.
Hugo W, Zaretsky J, Sun L, Song C, Moreno B, Hu-Lieskovan S, Berent-Maoz B, Pang J, Chmielowski B, Cherry G. Genomic and Transcriptomic Features of Response to Anti-PD-1 Therapy in Metastatic Melanoma. Cell. 2016;165(1):35–44.
Jemal A, Bray F, Center MM, Ferlay J, Forman D. Global cancerstatistics. CA Cancer J Clin. 2011;6(2):169–90.
Lightfoot AJ, Breyer BN, Rosevear HM, Erickson BA, Konety BR, O'Donnell MA. Multi-institutional analysis of sequential intravesical gemcitabine and mitomycin C chemotherapy for non-muscle invasive bladder cancer. Urol Oncol. 2014;32(1):35.e15–39.
Qian BZ, Li J, Hui Z, Kitamura T, Zhang J, Campion LR, Kaiser EA, Snyder LA, Pollard JW. CCL2 recruits inflammatory monocytes to facilitate breast-tumour metastasis. Nature. 2012;475(7355):222–5.
Richard N, Hanna, Caglar, Cekic, Duygu, Sag, Robert, Tacke, Graham D, Thomas: Patrolling monocytes control tumor metastasis to the lung. Science (New York, NY) 2015.
Wang R, Bao W, Pal M, Liu Y, Yazdanbakhsh K, Zhong H. Intermediate monocytes induced by IFN-γ inhibit cancer metastasis by promoting NK cell activation through FOXO1 and interleukin-27. J Immunotherapy Cancer. 2022;10(1).
Gordon S, Martinez FO. Alternative activation of macrophages: mechanism and functions. Immunity. 2010;32(5):593–604.
Sica A, Mantovani A. Macrophage plasticity and polarization: in vivo veritas. Macrophage plasticity and polarization: in vivo veritas. J Clin Invest. 2012;122(3):787–95.
Martínez VG, Rubio C, Martínez-Fernández M, Segovia C, López-Calderón F, Garín MI, Teijeira A, Munera-Maravilla E, Varas A, Sacedón R, et al. BMP4 Induces M2 Macrophage Polarization and Favors Tumor Progression in Bladder Cancer. Clin Cancer Res. 2017;23(23):7388–99.
Takeda K, Akira S. Role of toll-like receptor in innate immunity. Tanpakushitsu kakusan koso Protein, nucleic acid, enzyme. 2002;47(16 Suppl):2097–102.
Takeda K, Akira S. Toll-like receptors. Curr Protoc Immunol. 2015;109:14.12.11-14.12.10.
Takeda K. Toll-like receptors in innate immunity. Int Immunol. 2004;17(1):1.
Chung-Yi C, Chiu-Li K, Chi-Ming L. The Cancer Prevention, Anti-Inflammatory and Anti-Oxidation of Bioactive Phytochemicals Targeting the TLR4 Signaling Pathway. Int J Mol Ences. 2018;19(9):2729.
Solmaz O. Moghadam, Mohammad, Reza, Nowroozi: Toll-like receptors: The role in bladder cancer development, progression and immunotherapy. Scand J Immunol. 2019;90(6): e12818.
Zhao S, Zhang Y, Zhang Q, Wang F, Zhang D. Toll-like receptors and prostate cancer. Front Immunol. 2014;5:352.
Yang H, Zhou H, Feng P, Zhou X, Wen H, Xie X, Shen H, Zhu X. Reduced expression of Toll-like receptor 4 inhibits human breast cancer cells proliferation and inflammatory cytokines secretion. J Exp Clin Cancer Res. 2010;29(1):92.
Moradi-Marjaneh R, Hassanian SM, Fiuji H, Soleimanpour S, Ferns GA, Avan A, Khazaei M: Toll like receptor signaling pathway as a potential therapeutic target in colorectal cancer. J Cell Physiol. 2018.
Yusuf N, Nasti TH, Long JA, Naseemuddin M, Lucas AP, Xu H, Elmets CA. Protective role of Toll-like receptor 4 during the initiation stage of cutaneous chemical carcinogenesis. Can Res. 2008;68(2):615–22.
Ayari C, Bergeron A, Larue H, Ménard C, Fradet Y. Toll-like receptors in normal and malignant human bladders. J Urol. 2011;185(5):1915–21.
Bose S, Banerjee S, Mondal A, Chakraborty U, Pumarol J, Croley CR, Bishayee A. Targeting the JAK/STAT Signaling Pathway Using Phytocompounds for Cancer Prevention and Therapy. Cells. 2020;9(6).
Arumuggam N, Bhowmick NA, Rupasinghe H. A Review: Phytochemicals Targeting JAK/STAT Signaling and IDO Expression in Cancer. Phytother Res. 2015;29(6):805–17.
Thomas SJ, Snowden JA, Zeidler MP, Danson SJ. The role of JAK/STAT signalling in the pathogenesis, prognosis and treatment of solid tumours. Br J Cancer. 2015;113(3):365–71.
Napoli J. L: Cellular retinoid binding-proteins, CRBP, CRABP, FABP5: Effects on retinoid metabolism, function and related diseases. Pharmacol Ther. 2017;173:19–33.
Silvaroli JA, Arne JM, Chelstowska S, Kiser PD, Banerjee S, Golczak M. Ligand Binding Induces Conformational Changes in Human Cellular Retinol-binding Protein 1 (CRBP1) Revealed by Atomic Resolution Crystal Structures. J Biol Chem. 2016;291(16):8528–40.
Kuppumbatti YS, Bleiweiss IJ, Mandeli JP, Samuel W, Rafael MYL. Cellular Retinol-Binding Protein Expression and Breast Cancer. J Natl Cancer Inst. 2000;6:475.
Esteller M, Guo M, Moreno V, Peinado MA, Herman JG. Hypermethylation-associated Inactivation of the Cellular Retinol-Binding-Protein 1 Gene in Human Cancer. Can Res. 2002;62(20):5902.
Jerónimo C, Henrique R, Oliveira J, Lobo F, Lopes C. Aberrant cellular retinol binding protein (CRBP1) gene expression and promoter methylation in prostate cancer. J Clin Pathol. 2004;57(8):872–6.
Palan PR, Romney SL. Cellular binding proteins for vitamin A in human carcinomas and in normal tissues. Can Res. 1980;40(11):4221.
Toki K, Enokida H, Kawakami K, Chiyomaru T, Nakagawa M. CpG hypermethylation of cellular retinol-binding protein 1 contributes to cell proliferation and migration in bladder cancer. Int J Oncol. 2010;37(6):1379.
Ashis B. DNA Damage, Mutagenesis and Cancer. Int J Mol Sci. 2018;19(4):970.
Kondratova AA, Cheon H, Dong B, Holvey-Bates EG, Hasipek M, Taran I, Gaughan C, Jha BK, Silverman RH, Stark GR. Suppressing PAR ylation by 2′,5′,ligoadenylate synthetase 1 inhibits DNA damage- -nduced cell death. EMBO J. 2020;39(11):e101573.
Gonias SL, Karimi-Mostowfi N, Murray SS, Mantuano E, Gilder AS. Expression of LDL receptor-related proteins (LRPs) in common solid malignancies correlates with patient survival. PLoS ONE. 2017;12(10):e0186649.
Chen S, Bu G, Takei Y, Sakamoto K, Ikematsu S, Muramatsu T, Kadomatsu K. Midkine and LDL-receptor-related protein 1 contribute to the anchorage-independent cell growth of cancer cells. J Cell Sci. 2007;120(22):4009–15.
Bart VG, Stéphane D, Hervé E, Roebroek A. The Matricellular Receptor LRP1 Forms an Interface for Signaling and Endocytosis in Modulation of the Extracellular Tumor Environment. Front Pharmacol. 2015;6:271.
Song H, Li Y, Lee J, Schwartz AL, Bu G. Low-Density Lipoprotein Receptor-Related Protein 1 Promotes Cancer Cell Migration and Invasion by Inducing the Expression of Matrix Metalloproteinases 2 and 9. Can Res. 2009;69(3):879–86.
Montel V, Gaultier A, Lester RD, Campana WM, Gonias SL. The low-density lipoprotein receptor-related protein regulates cancer cell survival and metastasis development. Can Res. 2007;67(20):9817.
Wang Q, Zhu W, Xiao G, Ding M, Chang J, Liao H. Effect of AGER on the biological behavior of non-small cell lung cancer H1299 cells. Mol Med Rep. 2020;22(2):810–8.
Van D, Bos L, Van D. Epithelial Plasticity, Cancer Stem Cells, and the Tumor-Supportive Stroma in Bladder Carcinoma. Mol Cancer Res. 2012;10(8):995–1009.
Hanada T, Nakagawa M, Emoto A, Nomura T, Nomura Y. Prognostic value of tumor-associated macrophage count in human bladder cancer. Int J Urol. 2010;7(7):263–9.
Luís L. Daniela, Oliveira, Ana, Tavares, Teresina, Amaro, Ricardo, Cruz: The predominance of M2-polarized macrophages in the stroma of low-hypoxic bladder tumors is associated with BCG immunotherapy failure - ScienceDirect. Urol Oncol. 2014;32(4):449–57.
Nishiyama T, Tachibana M, Horiguchi Y, Nakamura K, Ikeda Y, Takesako K, Murai M. Immunotherapy of bladder cancer using autologous dendritic cells pulsed with human lymphocyte antigen-A24-specific MAGE-3 peptide. Clin Cancer Res. 2001;7(1):23.
Beatty JD, Islam S, North ME, Knight SC, Ogden CW. Urine dendritic cells: a noninvasive probe for immune activity in bladder cancer? BJU Int. 2015;94(9):1377–83.
Horn T, Grab J, Schusdziarra J, Schmid S, Maurer T, Nawroth R, Wolf P, Pritsch M, Gschwend JE, Kübler H. Antitumor T cell responses in bladder cancer are directed against a limited set of antigens and are modulated by regulatory T cells and routine treatment approaches. Int J Cancer J Int Du Cancer. 2013;133(9):2145–56.
Fridman WH, Pagès F, Sautès-Fridman C, Galon J. The immune contexture in human tumours: impact on clinical outcome. Nat Rev Cancer. 2012;12(4):298–306.
This work was supported by the Natural Science Foundation of Liaoning Province (2019-KF-01–12).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Xing, P., Jiang, Z. & Liu, Y. Construction and validation of a gene signature related to bladder urothelial carcinoma based on immune gene analysis. BMC Cancer 22, 926 (2022). https://doi.org/10.1186/s12885-022-09794-9
- Gene signature