Identification and Validation of an Individualized Prognostic Signature of Bladder Cancer Based on Seven Immune Related Genes

Background: There has been no report of prognostic signature based on immune-related genes (IRGs). This study aimed to develop an IRG-based prognostic signature that could stratify patients with bladder cancer (BLCA).

Methods: RNA-seq data along with clinical information on BLCA were retrieved from the Cancer Genome Atlas (TCGA) and gene expression omnibus (GEO). Based on TCGA dataset, differentially expressed IRGs were identified via Wilcoxon test. Among these genes, prognostic IRGs were identified using univariate Cox regression analysis. Subsequently, we split TCGA dataset into the training (n = 284) and test datasets (n = 119). Based on the training dataset, we built a least absolute shrinkage and selection operator (LASSO) penalized Cox proportional hazards regression model with multiple prognostic IRGs. It was validated in the training dataset, test dataset, and external dataset GSE13507 (n = 165). Additionally, we accessed the six types of tumor-infiltrating immune cells from Tumor Immune Estimation Resource (TIMER) website and analyzed the difference between risk groups. Further, we constructed and validated a nomogram to tailor treatment for patients with BLCA.

Results: A set of 47 prognostic IRGs was identified. LASSO regression and identified seven BLCA-specific prognostic IRGs, i.e., RBP7, PDGFRA, AHNAK, OAS1, RAC3, EDNRA, and SH3BP2. We developed an IRG-based prognostic signature that stratify BLCA patients into two subgroups with statistically different survival outcomes [hazard ratio (HR) = 10, 95% confidence interval (CI) = 5.6–19, P < 0.001]. The ROC curve analysis showed acceptable discrimination with AUCs of 0.711, 0.754, and 0.772 at 1-, 3-, and 5-year follow-up respectively. The predictive performance was validated in the train set, test set, and external dataset GSE13507. Besides, the increased infiltration of CD4⁺ T cells, CD8+ T cells, macrophage, neutrophil, and dendritic cells in the high-risk group (as defined by the signature) indicated chronic inflammation may reduce the survival chances of BLCA patients. The nomogram demonstrated to be clinically-relevant and effective with accurate prediction and positive net benefit.

Conclusion: The present immune-related signature can effectively classify BLCA patients into high-risk and low-risk groups in terms of survival rate, which may help select high-risk BLCA patients for more intensive treatment.

Introduction

Bladder cancer (BLCA) is the most common malignancy of the urinary system with high morbidity and mortality rates (Bray et al., 2018). Approximately 25% of BLCA patients are diagnosed with muscle-invasive or metastatic disease during the early stages of prognosis (Akbani et al., 2017). Meanwhile patients with non-muscle-invasive BLCA continue to suffer from the high progression rates (Cambier et al., 2016; Chen et al., 2019). Overall, the 5-year survival rate at all stages of bladder cancer remains no more than 20% (Aggen and Drake, 2017). Once the tumor progresses to locally advanced or metastatic stage, standard treatment for BLCA with combination chemotherapy are insufficient with low response and survival rates (Von et al., 2000; Maase et al., 2006). With the emergence of immune checkpoint therapy including programmed cell death protein (PD-1) and programmed death-ligand 1 (PD-L1), the treatment of advanced BLCA patients with durable response has become possible (Koshkin and Grivas, 2018). However, most BLCA patients do not adequately respond to PD-1 or PD-L1-targeted therapy; and hence, it is imperative to develop prognostic biomarkers to closely monitor progression and shed light on treatment stratification.

The most representative of BLCA type is urothelial cancer (UC). Up to 30% of urothelial cancer (UC) specimens have demonstrated differential expression in PD-L1, which is associated with increased all-cause mortality (Inman et al., 2007; Boorjian et al., 2008; Faraj et al., 2015). However, standalone PD-L1 expression acts as an unviable biomarker since significant heterogeneity of association were observed between PD-L1 staining and clinical results of BLCA patients (Rosenberg et al., 2016; Massard et al., 2017; Powles et al., 2017; Sharma et al., 2017; Bellmunt et al., 2017). Since PD-L1 expression is subject to immunohistochemistry (IHC) score; the expression of immune-related gene (IRG) may serve as a better biomarker as it can be quantified from multiple cell types within a sample (Sweis et al., 2016). In this context, the IRG-based prognostic signatures have been proposed for patients diagnosed with nonsquamous non–small cell lung cancer (Li et al., 2017) and papillary thyroid cancer (Lin et al., 2019), which show significant prognostic values. However, the clinical relevance and prognostic significance of IRGs-based signature in BLCA remains unknown.

The purpose of the present study was to investigate the clinical implications of IRGs on prognostic stratification and their potential as biomarkers for targeted BLCA therapy. In the manuscript, we performed integrated analysis using IRG expression profiles and clinical information of patients with BLCA retrospectively. Individualized prognostic signature based on IRGs was developed and validated in independent datasets, while underlying mechanisms were explored using bioinformatics analysis.

Materials and Methods

Data Source

From The Cancer Genome Atlas (TCGA), normalized RNA sequencing (RNA-Seq) data sets with estimation of Fragments Per Kilobase of transcript per Million mapped reads (FPKM) from 414 tumor samples and 19 non-tumor samples were retrieved. Clinical data were also derived from TCGA into integrated analysis. Gene expression data and clinical information in GSE13507 based on the Illumina human-6 v2.0 expression BeadChip platform were downloaded from gene expression omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/). A total of 568 patients were involved in the development and validation of the prognostic signature, i.e., 403 patients in TCGA dataset and 165 patients in GSE13507 cohort.

Data Preprocessing and Differential Analysis

In the R software, differentially expressed genes (DEGs) between BLCA and normal tissues were identified using Wilcoxon test after within-array replicate probes were replaced with their average via limma package (Ritchie et al., 2015; Yue et al., 2019). The p-value was adjusted with the false discovery rate (FDR) (Benjamini and Hochberg, 1995). FDR < 0.05 and |log2(FC)| value > 1 was considered significant.

The Kyoto Encyclopedia of Genes and Genomes (KEGG) (Kanehisa and Goto, 2000) pathway enrichment were analyzed with the DEGs using the clusterProfiler R package (Yu et al., 2012). P < 0.05 was considered statistically significant.

Development and Validation of a Prognostic Signature

By accessing the Immunology Database and Analysis Portal (IMMPORT) (Bhattacharya et al., 2014) website (https://www.immport.org), we retrieved a latest list of immune‐related genes, out of which we identified BLCA-specific immune‐related genes (IRGs) after matching the DEGs. Survival-associated IRGs were identified using univariate Cox regression analysis with a threshold value of p < 0.01.

Patients in TCGA dataset was randomly assigned in a 7:3 ratio to a training set and test set with the same proportion of each BLCA stage. With expression profiles of the identified survival-associated IRGs, we conducted least absolute shrinkage and selection operator (LASSO) regression analysis in the training set. Subsequently we calculated the individualized risk score with coefficients and constructed a prognostic signature which separates the high-risk BLCA patients from the low-risk group. Clinical relevance was validated using survival analysis between groups with thresholds of p < 0.05 using the R software survival and survminer package; whereas, the receiver operating characteristic (ROC) analysis was performed (via the survival ROC package), and the area under the curve (AUC) was calculated at multiple time-point to evaluate the discrimination (Heagerty et al., 2000).

Clinical characteristics including age, gender, stage, and tumor-node-metastasis (TNM) status were collected from TCGA database and integrated with transcriptome profile derived from TCGA dataset. Multivariate cox regression analysis was performed using clinical data and risk scores to see if the prognostic value of risk scores was independent of clinical characteristics. A value of p < 0.05 was considered significant statistically.

External Validation of the Prognostic Signature in the Test Set and GSE13507 Cohort

The prognostic signature with the same risk score formula and cutoff value was then validated in the test set and external dataset GSE13507 respectively. We chose the GSE13507 as it involved the largest sample size of BLCA patients with survival data. Likewise, the prognostic model was presented as a risk plot in each dataset that encompassed the expression level of the included genes, distribution of risk score, and survival status of individuals.

Gene Set Enrichment Analysis Analysis

We performed Gene Set Enrichment Analysis (GSEA) between high-risk and low-risk group as separated by the 7-IRG signature via clusterProfiler and enrichplot R package. Two functions (gseGO and gseKEGG) were applied to identify the enriched terms in Gene Ontology (GO) and KEGG with a false discovery rate (FDR) value < 0.05 (Subramanian et al., 2005).

Difference of Tumor-Infiltrating Immune Cells Between Groups Defined by the Signature in Bladder Cancer

Six types of tumor-infiltrating immune cells were retrieved from Tumor Immune Estimation Resource (TIMER) (https://cistrome.shinyapps.io/timer/), (Li et al., 2017) also known as a web server for comprehensive analysis of tumor-infiltrating immune cells. The abundance of immune cells was tested one by one to detect the differences between the prognostic classified risk groups using Wilcoxon test.

Construction of a Nomogram Based on the Immune-Related Gene Signature

A nomogram encompassing the risk score based on expression of prognostic IRGs and clinicopathological factors was constructed using the rms R package. Discrimination of the nomogram was validated using ROC analysis at 3- and 5-year follow-up, and predictive accuracy was tested by presenting the difference between predicted survival and actual survival using calibration plot. Further, decision curve analysis (DCA) was performed to examine the clinical utility of the nomogram by quantifying the net benefits at different threshold probabilities.

Results

Bladder Cancer Specific Immune-Related Genes

We managed to obtain 4,876 DEGs base on the transcription profile in TCGA dataset, out of which 3,453 genes were upregulated and 1,423 downregulated. The DEGs were presented in Supplementary Material S1. The KEGG analysis (Figure 1) indicated that the genes were mainly involved in PI3K−Akt and MAPK signaling pathway, which are pivotal in the regulation of immune responses (Liu et al., 2007; Weichhart and Säemann, 2008). Specific IRGs were identified with intersection of the immune-related genes and DEGs in BLCA. As shown in Figure 2, we identified 120 upregulated, and 140 downregulated BLCA-specific IRGs.

FIGURE 1

Figure 1 The top 30 most enriched Kyoto Encyclopedia of Genes and Genomes pathways.

FIGURE 2

Figure 2 Differentially expressed immune‐related genes. Heatmap (A) and a volcano plot (B).

Development and Internal Validation of a Prognostic Signature

A total of 403 patients with BLCA in TCGA dataset was recorded with follow-up time >0. Based on the data in these patients, totally 47 survival-associated IRGs were identified using univariate Cox regression. Table 1 shows the general profile of survival-associated IRGs in BLCA. Subsequently, we split TCGA dataset randomly into training set (n = 284) and test set (n = 119). With expression profiles of survival-associated IRGs, we conducted a LASSO regression in the training set and identified seven BLCA-specific prognostic IRGs (Figures 3A, B). The boxplot of expression levels of these genes between BLCA and normal samples in TCGA dataset along with results of Wilcoxon test is provided in Supplementary Material S2.

TABLE 1

Table 1 Characteristics of bladder cancer‐specific immune‐related genes (univariate cox analysis).

FIGURE 3

Figure 3 Development of the prognostic signature based on seven immune‐related genes (IRGs) in the training set. (A, B) LASSO regression identified 7 IRGs (C) Time-dependent ROC curve of the 7-IRG prognostic signature (D) Survival analysis between signature-defined risk groups (E) Heatmap of expression profiles of included IRGs (F) Distribution of groups based on the signature (G) Survival status of patients in different groups.

After extracting the coefficient values, we calculated individualized risk scores with coefficient‐weighted expression levels of seven IRGs with the following formula:

$$Riskscore=expressionlevelofRBP7\ast 0.0044+PDGFRA\ast 0.0021+AHNAK\ast 0.0054+OAS1\ast (-0.0021)+RAC3\ast 0.0043+EDNRA\ast 0.0142+SH3BP2\ast (-0.0019)$$

The median of the risk score is 0.24. Individuals with a risk score higher than 0.24 were classified as a high-risk group while the other as a low-risk group; consequently, a prognostic signature based on seven IRGs was developed. The survival analysis indicated that the survival rate was remarkably lower in the high‐risk group as opposed to low‐risk group [hazard ratio (HR) = 10, 95% confidence interval (CI) = 5.6–19, p-value <0.001, Figure 3C]; whereas, the ROC curve analysis (Figure 3D) showed acceptable discrimination with AUCs of 0.711, 0.754, and 0.772 at 1-, 3-, and 5-year follow-up respectively. Figures 3E–G represents the risk plot encompassing distribution of groups based on the signature, survival status of individuals between groups, and the expression level of included IRGs. It shows a clear separation of survival status between risk groups with red dots being death and blue ones alive. While a large amount of death occurred in the high-risk group, most of the patients in low-risk group stayed alive at follow-up.

We calculated individual risk score with the aforementioned formula and classified the patients in the test set into high-risk and low-risk groups. Similarly, we validated the clinical utility and discrimination in both datasets. Figure 4 summarizes the results of validation in the test set. A significant separation was shown in the Kaplan-Meier survival curve in the test set (HR = 6.9, 95% CI = 1.8–27, p-value = 0.01, Figure 4A). ROC curve analysis demonstrated acceptable discrimination with an AUC of 0.68 in predicting 5-year overall survival (Figure 4B). Meanwhile, the risk plot shows markedly different survival status between groups (Figures 4C–E).

FIGURE 4

Figure 4 Validation of the the seven immune‐related gene (7-IRG) prognostic signature in the test set. (A) Time-dependent receiver operating characteristic (ROC) curve of the 7-IRG prognostic signature. (B) Survival analysis between signature-defined risk groups. (C) Heatmap of expression profiles of included IRGs. (D) Distribution of groups based on the signature. (E) Survival status of patients in different groups. (F) Multivariable analyses of the risk score, age, gender, tumor stage, and T/N/M, where confidence interval (CI) and hazard ratio (HR) stand for a confidence interval, and hazard ratio, respectively.

Multivariate Cox regression analysis with TCGA dataset suggested that the 7-IRG based prognostic signature could be an independent predictor after other variables including age, gender, stage, and information coded in TNM concerning status of tumor, lymph node, and distant metastasis were adjusted (Figure 4F).

External Validation in the Test and GSE13507 Datasets

In GSE13507 cohort, we validated the IRG-based prognostic signature for overall survival, cancer specific survival and progression-free survival using the same formula and cutoff value. The results were presented in Figure 5 with Kaplan-Meier curve and time-dependent ROC curve. Consistent with above findings, significantly different survival outcomes were observed between risk groups in overall survival (OS) (HR = 73,000, 95% CI = 0.53–10^10, p-value = 0.01, Figure 5A), cancer-specific survival (CSS) (HR = 8.4*10^9, 95% CI = 2,700–2.6*10^16, p-value <0.001, Figure 5B), and progression-free survival (PFS) (HR = 1.3*10^8, 95% CI = 11-1.7*10^15, p-value = 0.02, Figure 5C). Similarly, ROC curve analysis indicated that the prognostic signature can effectively predict 65% of the overall survival, 71% of the cancer specific survival, and 71% of the progression-free survival at 1- and 3-year follow-up (Figures 5D, E). Risk plots for different survival indicators were presented in Supplementary Material S3.

FIGURE 5

Figure 5 Validation of the the seven immune‐related gene (7-IRG) prognostic signature in GSE13507 cohort. (A–C) Survival analysis between signature-defined risk groups [overall survival (OS) (A), cancer-specific survival (CSS) (B), and progression-free survival (PFS) (C)]. (D–F) Time-dependent receiver operating characteristic (ROC) curve for predicting survival of the 7-IRG prognostic signature [OS (D), CSS (E), and PFS (F)].

Gene Set Enrichment Analysis

The predictive performance of the 7-IRG signature could be associated with the biological function of these IRGs in BLCA. To explore the underlying mechanism, we performed GSEA between high-risk and low-risk groups based on the signature to identify the enriched GO term as well as KEGG pathway. Figure 6 shows the results of enrichment analysis. GO terms (Figure 6A) including collagen−activated tyrosine kinase receptor signaling pathway, dendritic cell antigen processing and presentation, extracellular matrix component, as well as leukocyte migration and chemotaxis involved in inflammatory response were significantly enriched. In the case of KEGG pathway (Figure 6B), we found enhanced activity of several immune-related pathways in the high-risk group, such as chemokine signaling pathway, cytokine−cytokine receptor interaction, ECM−receptor interaction, IL−17 signaling pathway, leukocyte transendothelial migration, and PI3K−Akt signaling pathway. Particularly, upregulated PI3K−Akt pathway was in consistency with the KEGG enrichment of DEGs between tumor and control specimens in BLCA (Figure 2).

FIGURE 6

Figure 6 Gene Set Enrichment Analysis (GSEA) analysis between risk groups as classified by the immune-related gene (IRG)-based signature. (A) GSEA in gene ontology (GO) terms. (B) GSEA in Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways.

Difference of Tumor-Infiltrating Immune Cells Between Risk Groups

We analyzed the difference in tumor-infiltrating immune cells in TCGA samples between the risk groups to explore the relationship between the present IRG-based prognostic signature and the tumor immune microenvironment. The results showed that abundance of six types of tumor-infiltrating immune cells except B cells were significantly enriched (P < 0.001) in the high risk group as compared to the low risk group by Wilcoxon test. By contrast, higher B cell level was observed in the low-risk group (P < 0.001) (Figure 7).

FIGURE 7

Figure 7 The difference of tumor-infiltrating immune cells among risk groups as defined by the seven immune‐related gene (7-IRG) prognostic signature.

Construction of a Nomogram Based on the Seven Immune-Related Gene Signature

We developed a nomogram to predict 1-, 3-, and 5-year overall survival using the 7-IRG signature and the aforementioned clinical factors of colorectal cancer (Figure 8A). The ROC analysis (Figure 8B) showed adequate discrimination with an AUC of 0.759, and 0.783 at 3- and 5-year follow-up, indicating that the nomogram could distinguish over 75% of survival outcome at these time-points. In addition, the calibration plot (Figure 8C-D) demonstrated optimal predictive accuracy with predicted survival rate approximately equivalent to actual survival. Further, results of DCA (Figure 8E-F) showed that most of the dashed curve were above the two solid lines (black and gray), suggesting positive net benefits. In other words, clinical decision made upon the nomogram would be favorable than treat-none or treat-all scheme.

FIGURE 8

Figure 8 Construction of a nomogram based on the seven immune‐related gene (7-IRG) signature. (A) A nomogram based on the signature and clinical information. (B) Time-dependent receiver operating characteristic (ROC) curve for predicting overall survival (OS) of the nomogram. (C, D) Calibration plot evaluating the predictive accuracy of the nomogram [at 3-year survival (C) at 5-year survival (D)]. (E, F) Decision curve analysis evaluating the clinical utility of the nomogram [at 3-year survival (E) at 5-year survival (F)].

Discussion

BLCA patients are at substantial risk for recurrence and metastasis. Over the past several decades, treatments for metastatic BLCA remained unsatisfied with platinum chemotherapy until immune checkpoint therapy came up (Bellmunt et al., 2014; Milowsky et al., 2016; Ghasemzadeh et al., 2016). With the growing popularity of immunotherapy, attention has been shifting to developing novel biomarkers related to tumor immune milieus for estimates of treatment response and survival outcome (Aggen and Drake, 2017).

Till date, no report on IRG-based prognostic signatures exists, although considerable efforts have been made to develop prognostic signatures based on differentially expressed genes (Gao et al., 2019; Yan et al., 2019; Zhang et al., 2019). We identified seven prognostic IRGs based on a comprehensive analysis which could serve as valuable biomarkers in the clinical setting. Further, the prognostic signature based on the seven IRGs could categorize BLCA patients into two subgroups with statistically different survival outcomes, which was validated in both TCGA and GSE13507 datasets. Multivariate Cox analysis indicated that the IRG-based prognostic signature was independent of clinical factors. Further, it presented discrimination in overall survival rate, cancer-specific survival and progression-free survival in GSE13507 cohort, indicating the 7-IRG prognostic signature was independent of datasets. Using the 7-IRG signature along with age, gender and TNM_stage, we built a nomogram to predict overall survival for patients with BLCA. The nomogram promised to be clinically-relevant and credible in predictive performance. Additionally, we explored the underlying mechanisms using GSEA and immune cell analysis between risk groups. These findings support the potential translation of the present IRG-based prognostic signatures into clinical practice.

Our IRG-based signature highlighted seven IRGs, i.e., RBP7, PDGFRA, AHNAK, OAS1, RAC3, EDNRA, and SH3BP2. Expression of PDGFRA was reported in BLCA specimens; however, such report lacked sufficient experimental study (Terada, 2009; Terada, 2013). AHNAK was identified as a unique intracellular protein with different expression level and subcellular localization between BLCA sample and control (Lee et al., 2018); whereas, its prognostic value was observed in other studies (Chu et al., 2018; Wu et al., 2019). High expression of EDNRA was found to be associated with poor outcome in patients with advanced BLCA in a bioinformatics report (Laurberg et al., 2014). By contrast, no reports concerning RBP7, OAS1, RAC3, and SH3BP2 were published in BLCA, and therefore, the role of these IRGs in BLCA requires further investigation.

To explore the mechanisms by which the IRG-based signature effectively stratifies BLCA patients, GSEA between risk groups as classified by the signature demonstrated significant activity of multiple immune-related pathways in the high-risk patients. DEGs between risk groups were significantly enriched in chemokine signaling pathway and cytokine-receptor interaction, which are involved in chemotaxis, angiogenesis, as well as inflammatory processes (Lippitz, 2013; Li et al., 2017). The cytokines play a crucial role in the immune response of cytotoxic T lymphocytes (CTL), and the regulation of cellular differentiation (Agarwal et al., 2006). Besides, an enhanced inflammatory milieu is reported to be associated with tumor progression and poor prognosis (Lin and Karin, 2007; Lippitz, 2013). To further elucidate the role of tumor milieu associated with the 7-IRG signature, we analyzed the estimations of six types of immune cells and observed increased abundance of CD4+ T cells, CD8+ T cells, macrophage, neutrophil and dendritic cells in the high-risk group. T cells infiltration was reported to promote tumor invasion and metastasis via the androgen receptor (Tao et al., 2016) and estrogen receptor signaling (Tao et al., 2018) among BLCA patients. This aspect aligned with a clinical study that demonstrated the prognostic role of a specific subset of CD4+ T cells (Th17) (Chugh et al., 2013). The infiltration of tumor-associated macrophages is recognized to facilitate tumor progression in BLCA via tumorigenesis, angiogenesis, and disruption of adaptive immune response (Miyake et al., 2016). Moreover, the recruitment of neutrophils could increase the level of human neutrophil peptides, which cause tumor angiogenesis and growth (Thompson et al., 2015). An increase of tumor-infiltrating neutrophils was linked to immunosuppression in muscle invasion of BLCA patients that lead to inadequate treatment response and prognosis (Zhou et al., 2017). Additionally, the presence of tumor-infiltrated dendritic cells is associated with tumor progression and poor prognosis. By contrast, higher B cell level in the low-risk group indicated the tumor-suppressive role of B cell infiltration, which was consistent with a previous study (Jiang et al., 2019). In summary, the misregulation of tumor immune microenvironment may be responsible for the difference in survival outcome observed between groups as defined by the prognostic signature.

In the present study, we reported a prognostic signature based on IRG expression for predicting survival rate in BLCA patients, which was observed to be clinically relevant and effective in different datasets. To our best knowledge, this has been the first reported IRG-based signature in BLCA. Nevertheless, our results consisted of several limitations. First, as non-tumor samples were less than BLCA specimen, the results were biased to an extent. Second, the molecular mechanisms of BLCA could not be fully elucidated due to lack of in vitro and in vivo studies. Further studies are therefore warranted.

Conclusions

We developed and validated a first-ever IRG-based prognostic signature that stratify BLCA patients into two subgroups with statistically different survival outcomes, for which misregulation of tumor immune microenvironment might be responsible. These findings may provide insight on development of novel immune biomarkers and target therapy.

Data Availability Statement

Expression profile can be accessed in The Cancer Genome Atlas (TCGA) (https://portal.gdc.cancer.gov/) and gene expression omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/), while immune related gene list can be retrieved from the Immunology Database and Analysis Portal (IMMPORT) website (https://www.immport.org). Abundance of tumor-infiltrating immune cells were retrieved from the Tumor Immune Estimation Resource (TIMER) website (https://cistrome.shinyapps.io/timer/). The data sets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Author Contributions

HQ, XH, CH, YL, and BY contributed to data analysis, interpretation, and drafting. HQ, XH, CH, YL, and JL contributed to study design, study supervision, and final approval of the manuscript. All authors read and approved the final manuscript.

Funding

The study was supported by The Introduced Project of Suzhou Clinical Medical Expert Team (SZYJTD201725, recipient: JL) and Key Project of Jiangsu Provincial Department of Science and Technology (BE2017007-5, recipient: YL).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We would like to thank Richard Shi from Nanjing Medical University for his valuable lessons in medical English writing.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.00012/full#supplementary-material

Supplementary Table S1 | Differentially expressed genes between BLCA samples and non-tumor samples.

Supplementary Figure S2 | Boxplot of expression levels in these genes between BLCA and normal samples in TCGA dataset with results of Wilcoxon test.

Supplementary Figure S3 | Riskplots for the 7-IRG signature in GSE13507 [(1)overall survival; (2) cancer-specific survival; (3) progression-free survival] (A) Heatmap of expression profiles of included IRGs (B) Distribution of groups based on the signature (C) Survival status of patients in different groups.

References

Agarwal, A., Verma, S., Burra, U., Murthy, N., Mohanty, N., Saxena, S. (2006). Flow cytometric analysis of Th1 and Th2 cytokines in PBMCs as a parameter of immunological dysfunction in patients of superficial transitional cell carcinoma of bladder. Cancer Immunol. Immunother. 55 (6), 734–743. doi: 10.1007/s00262-005-0045-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Aggen, D. H., Drake, C. G. (2017). Biomarkers for immunotherapy in bladder cancer: a moving target. J. Immunother. Cancer 5 (1), 94. doi: 10.1186/s40425-017-0299-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Aggen, D. H., Drake, C. G. (2017). Biomarkers for immunotherapy in bladder cancer: a moving target. J. Immunother. Cancer 5 (1), 94. doi: 10.1186/s40425-017-0299-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Akbani, R., Al-Ahmadie, H., Albert, M., Alexopoulou, I., Ally, A., Antic, T., et al. (2017). Comprehensive molecular characterization of muscle-invasive bladder cancer. Cell 171 (3), págs. 540–556. doi: 10.1016/j.cell.2017.09.007

CrossRef Full Text | Google Scholar

Bellmunt, J., Orsola, A., Leow, J., Wiegel, T., De Santis, M., Horwich, A. (2014). Bladder cancer: esmo practice guidelines for diagnosis, treatment and follow-up. Ann. Oncol. 25 (suppl_3), iii40–iii48. doi: 10.1093/annonc/mdu223

PubMed Abstract | CrossRef Full Text | Google Scholar

Bellmunt, J., De Wit, R., Vaughn, D. J., Fradet, Y., Lee, J.-L., Fong, L., et al. (2017). Pembrolizumab as second-line therapy for advanced urothelial carcinoma. New Engl. J. Med. 376 (11), 1015–1026. doi: 10.1056/NEJMoa1613683

CrossRef Full Text | Google Scholar

Benjamini, Y., Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc.: Ser. B. (Methodological) 57 (1), 289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x

CrossRef Full Text | Google Scholar

Bhattacharya, S., Andorf, S., Gomes, L., Dunn, P., Schaefer, H., Pontius, J., et al. (2014). ImmPort: disseminating data to the public for the future of immunology. Immunol. Res. 58 (2-3), 234–239. doi: 10.1007/s12026-014-8516-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Boorjian, S. A., Sheinin, Y., Crispen, P. L., Farmer, S. A., Lohse, C. M., Kuntz, S. M., et al. (2008). T-cell coregulatory molecule expression in urothelial cell carcinoma: clinicopathologic correlations and association with survival. Clin. Cancer Res. Off. J. Am. Assoc. Cancer Res. 14 (15), 4800. doi: 10.1158/1078-0432.CCR-08-0731

CrossRef Full Text | Google Scholar

Bray, F., Ferlay, J., Soerjomataram, I., Siegel, R. L., Torre, L. A., Jemal, A. (2018). Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA: A Cancer J. Clin. 68 (6), 394–424. doi: 10.3322/caac.21492

CrossRef Full Text | Google Scholar

Cambier, S., Sylvester, R. J., Collette, L., Gontero, P., Brausi, M. A., Van Andel, G., et al. (2016). EORTC nomograms and risk groups for predicting recurrence, progression, and disease-specific and overall survival in non–muscle-invasive stage Ta–T1 urothelial bladder cancer patients treated with 1–3 years of maintenance bacillus Calmette-Guérin. Eur. Urol. 69 (1), 60–69. doi: 10.1016/j.eururo.2015.06.045

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Z., He, S., Zhan, Y., He, A., Fang, D., Gong, Y., et al. (2019). TGF-β-Induced transgelin promotes bladder cancer metastasis by regulating epithelial-mesenchymal transition and invadopodia formation. EBioMedicine 47, 208–220. doi: 10.1016/j.ebiom.2019.08.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Chu, J., Li, N., Li, F. (2018). A risk score staging system based on the expression of seven genes predicts the outcome of bladder cancer. Oncol. Lett. 16 (2), 2091–2096. doi: 10.3892/ol.2018.8904

PubMed Abstract | CrossRef Full Text | Google Scholar

Chugh, S., Anand, V., Swaroop, L., Sharma, M., Seth, A., Sharma, A. (2013). Involvement of Th17 cells in patients of urothelial carcinoma of bladder. Hum. Immunol. 74 (10), 1258–1262. doi: 10.1016/j.humimm.2013.06.032

PubMed Abstract | CrossRef Full Text | Google Scholar

Faraj, S. F., Munari, E., Guner, G., Taube, J., Anders, R., Hicks, J., et al. (2015). Assessment of Tumoral PD-L1 expression and intratumoral CD8+ TCells in Urothelial Carcinoma. Urology 85 (3), 703.e1–703.e6. doi: 10.1016/j.urology.2014.10.020

CrossRef Full Text | Google Scholar

Gao, X., Zhang, S., Chen, Y., Wen, X., Chen, M., Wang, S., et al. (2019). Development of a novel six-long noncoding RNA signature predicting survival of patients with bladder urothelial carcinoma. J. Cell. Biochem 120 (12), 19796–19809. doi: 10.1002/jcb.29285

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghasemzadeh, A., Bivalacqua, T. J., Hahn, N. M., Drake, C. G. (2016). New strategies in bladder cancer: a second coming for immunotherapy. Clin. Cancer Res. 22 (4), 793–801. doi: 10.1158/1078-0432.CCR-15-1135

PubMed Abstract | CrossRef Full Text | Google Scholar

Heagerty, P. J., Lumley, T., Pepe, M. S. (2000). Time-dependent ROC curves for censored survival data and a diagnostic marker. Biometrics 56 (2), 337–344. doi: 10.1111/j.0006-341X.2000.00337.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Inman, B. A., Sebo, T. J., Frigola, X., Dong, H., Bergstralh, E. J., Frank, I., et al. (2007). PD-L1 (B7-H1) expression by urothelial carcinoma of the bladder and BCG-induced granulomata: associations with localized stage progression. Cancer: Interdiscip. Int. J. Am. Cancer Soc. 109 (8), 1499–1505. doi: 10.1002/cncr.22588

CrossRef Full Text | Google Scholar

Jiang, Q., Fu, Q., Chang, Y., Liu, Z., Zhang, J., Xu, L., et al. (2019). CD19+ tumor-infiltrating B-cells prime CD4+ T-cell immunity and predict platinum-based chemotherapy efficacy in muscle-invasive bladder cancer. Cancer Immunol. Immunother. 68 (1), 45–56. doi: 10.1007/s00262-018-2250-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanehisa, M., Goto, S. (2000). KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28 (1), 27–30. doi: 10.1093/nar/28.1.27

PubMed Abstract | CrossRef Full Text | Google Scholar

Koshkin, V. S., Grivas, P. (2018). Emerging role of immunotherapy in advanced urothelial carcinoma. Curr. Oncol. Rep. 20 (6), 48. doi: 10.1007/s11912-018-0693-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Laurberg, J. R., Jensen, J. B., Schepeler, T., Borre, M., Ørntoft, T. F., Dyrskjøt, L. (2014). High expression of GEM and EDNRA is associated with metastasis and poor outcome in patients with advanced bladder cancer. BMC Cancer 14 (1), 638. doi: 10.1186/1471-2407-14-638

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, H., Kim, K., Woo, J., Park, J., Kim, H., Lee, K. E., et al. (2018). Quantitative proteomic analysis identifies AHNAK (Neuroblast Differentiation-associated Protein AHNAK) as a novel candidate biomarker for bladder urothelial carcinoma diagnosis by liquid-based cytology. Mol. Cell. Proteomics Mcp 17 (9), 1788. doi: 10.1074/mcp.RA118.000562

CrossRef Full Text | Google Scholar

Li, B., Cui, Y., Diehn, M., Li, R. (2017). Development and validation of an individualized immune prognostic signature in early-stage nonsquamous non-small cell lung cancer. JAMA Oncol. 3 (11), 1529. doi: 10.1001/jamaoncol.2017.1609

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, T., Fan, J., Wang, B., Traugh, N., Chen, Q., Liu, J. S., et al. (2017). TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 77 (21), e108–e110. doi: 10.1158/0008-5472.CAN-17-0307

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, W.-W., Karin, M. (2007). A cytokine-mediated link between innate immunity, inflammation, and cancer. J. Clin. Invest. 117 (5), 1175–1183. doi: 10.1172/JCI31537

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, P., Guo, Y.-n., Shi, L., Li, X.-j., Yang, H., He, Y., et al. (2019). Development of a prognostic index based on an immunogenomic landscape analysis of papillary thyroid cancer. Aging (Albany NY) 11 (2), 480. doi: 10.18632/aging.101754

PubMed Abstract | CrossRef Full Text | Google Scholar

Lippitz, B. E. (2013). Cytokine patterns in patients with cancer: a systematic review. Lancet Oncol. 14 (6), e218–e228. doi: 10.1016/S1470-2045(12)70582-X

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y., Shepherd, E. G., Nelin, L. D. (2007). MAPK phosphatases—regulating the immune response. Nat. Rev. Immunol. 7 (3), 202. doi: 10.1038/nri2035

PubMed Abstract | CrossRef Full Text | Google Scholar

Maase, H. V. D., Sengelov, L., Roberts, J. T., Ricci, S., Dogliotti, L., Oliver, T., et al. (2006). Long-term survival results of a randomized trial comparing gemcitabine plus cisplatin, with methotrexate, vinblastine, doxorubicin, plus cisplatin in patients with bladder cancer. J. Urol. 175 (2), 482–482. doi: 10.1016/S0022-5347(05)00370-8

CrossRef Full Text | Google Scholar

Massard, C., Gordon, M. S., Sharma, S., Rafii, S., Wainberg, Z. A., Luke, J., et al. (2017). Safety and efficacy of durvalumab (MEDI4736), an anti-programmed cell death ligand-1 immune checkpoint inhibitor, in patients with advanced urothelial bladder cancer. J. Clin. Oncol. 34 (26), 3119–3125. doi: 10.1200/JCO.2016.67.9761

CrossRef Full Text | Google Scholar

Milowsky, M. I., Rumble, R. B., Booth, C. M., Gilligan, T., Eapen, L. J., Hauke, R. J., et al. (2016). Guideline on muscle-invasive and metastatic bladder cancer (European Association of Urology guideline): American society of clinical oncology clinical practice guideline endorsement. J. Clin. Oncol. 34 (16), 1945–1952. doi: 10.1200/JCO.2015.65.9797

PubMed Abstract | CrossRef Full Text | Google Scholar

Miyake, M., Hori, S., Morizawa, Y., Tatsumi, Y., Nakai, Y., Anai, S., et al. (2016). CXCL1-mediated interaction of cancer cells with tumor-associated macrophages and cancer-associated fibroblasts promotes tumor progression in human bladder cancer. Neoplasia 18 (10), 636–646. doi: 10.1016/j.neo.2016.08.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Powles, T., O’Donnell, P. H., Massard, C., Arkenau, H.-T., Friedlander, T. W., Hoimes, C., et al. (2017). Updated efficacy and tolerability of durvalumab in locally advanced or metastatic urothelial carcinoma. J. Clin. Oncol. 35 (6), 286–286. doi: 10.1200/JCO.2017.35.6_suppl.286

CrossRef Full Text | Google Scholar

Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., et al. (2015). limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 43 (7), e47–e47. doi: 10.1093/nar/gkv007

PubMed Abstract | CrossRef Full Text | Google Scholar

Rosenberg, J. E., Hoffman-Censits, J., Powles, T., Van Der Heijden, M. S., Balar, A. V., Necchi, A., et al. (2016). 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 387 (10031), 1909–1920. doi: 10.1016/S0140-6736(16)00561-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, P., Retz, M., Siefker-Radtke, A., Baron, A., Necchi, A., Bedke, J., et al. (2017). Nivolumab in metastatic urothelial carcinoma after platinum therapy (CheckMate 275): a multicentre, single-arm, phase 2 trial. Lancet Oncol. 18 (3), 312–322. doi: 10.1016/S1470-2045(17)30065-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Subramanian, A., Tamayo, P., Mootha, V. K., Mukherjee, S., Ebert, B. L., Gillette, M. A., et al. (2005). Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. 102 (43), 15545–15550. doi: 10.1073/pnas.0506580102

CrossRef Full Text | Google Scholar

Sweis, R. F., Spranger, S., Bao, R., Paner, G. P., Stadler, W. M., Steinberg, G., et al. (2016). Molecular drivers of the non-T-cell-inflamed tumor microenvironment in urothelial bladder cancer. Cancer Immunol. Res. 4 (7), 563–568. doi: 10.1158/2326-6066.CIR-15-0274

PubMed Abstract | CrossRef Full Text | Google Scholar

Tao, L., Qiu, J., Jiang, M., Song, W., Yeh, S., Yu, H., et al. (2016). Infiltrated T cells promote bladder cancer progression via increasing IL-1→androgen receptor (AR)→HIF-1α→VEGFa signals. Mol. Cancer Ther. 15 (8), 1943. doi: 10.1158/1535-7163.MCT-15-0306

PubMed Abstract | CrossRef Full Text | Google Scholar

Tao, L., Qiu, J., Slavin, S., Ou, Z., Liu, Z., Ge, J., et al. (2018). Recruited T cells promote the bladder cancer metastasis via up-regulation of the estrogen receptor β/IL-1/c-MET signals. Cancer Lett. 430, 215–223. doi: 10.1016/j.canlet.2018.03.045

PubMed Abstract | CrossRef Full Text | Google Scholar

Terada, T. (2009). Autopsy case of primary small cell carcinoma of the urinary bladder: KIT and PDGFRA expression and mutations. Pathol. Int. 59 (4), 247–250. doi: 10.1111/j.1440-1827.2009.02358.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Terada, T. (2013). Urinary bladder urothelial carcinoma with expression of KIT and PDGFRA and showing diverse differentiations into plasmacytoid, clear cell, acantholytic, nested, and spindle variants, and into adenocarcinoma, signet-ring cell carcinoma, small cell carcinoma, large cell carcinoma, and pleomorphic carcinoma. Int. J. Clin. Exp. Pathol. 6 (6), 1150.

PubMed Abstract | Google Scholar

Thompson, D. B., Siref, L. E., Feloney, M. P., Hauke, R. J., Agrawal, D. K. (2015). Immunological basis in the pathogenesis and treatment of bladder cancer. Expert Rev. Clin. Immunol. 11 (2), 265. doi: 10.1586/1744666X.2015.983082

PubMed Abstract | CrossRef Full Text | Google Scholar

Von, d. M. H., Hansen, S. W., Roberts, J. T., Dogliotti, L., Oliver, T., Moore, M. J., et al. (2000). Gemcitabine and cisplatin versus methotrexate, vinblastine, doxorubicin, and cisplatin in advanced or metastatic bladder cancer: results of a large, randomized, multinational, multicenter, phase III study. J. Clin. Oncol. Off. J. Am. Soc. Clin. Oncol. 18 (17), 3068. doi: 10.1200/JCO.2000.18.17.3068

CrossRef Full Text | Google Scholar

Weichhart, T., Säemann, M. (2008). The PI3K/Akt/mTOR pathway in innate immune cells: emerging therapeutic applications. Ann. Rheumatic Dis. 67 (Suppl 3), iii70–iii74. doi: 10.1136/ard.2008.098459

CrossRef Full Text | Google Scholar

Wu, Q., Fu, C., Li, M., Li, J., Li, Z., Qi, L., et al. (2019). CINP is a novel cofactor of KLF5 required for its role in the promotion of cell proliferation, survival and tumor growth. Int. J. Cancer 144 (3), 582–594. doi: 10.1002/ijc.31908

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, X., Liu, X.-P., Guo, Z.-X., Liu, T.-Z., Li, S. (2019). Identification of hub genes associated with progression and prognosis in patients with bladder cancer. Front. In Genet. 10, 408. doi: 10.3389/fgene.2019.00408

CrossRef Full Text | Google Scholar

Yu, G., Wang, L.-G., Han, Y., He, Q.-Y. (2012). clusterProfiler: an R package for comparing biological themes among gene clusters. Omics: A J. Integr. Biol. 16 (5), 284–287. doi: 10.1089/omi.2011.0118

CrossRef Full Text | Google Scholar

Yue, C., Ma, H., Zhou, Y. (2019). Identification of prognostic gene signature associated with microenvironment of lung adenocarcinoma. PeerJ 7, e8128. doi: 10.7717/peerj.8128

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, C., Li, Z., Hu, J., Qi, F., Li, X., Luo, J. (2019). Identification of five long noncoding RNAs signature and risk score for prognosis of bladder urothelial carcinoma. J. Cell. Biochem. 121 (1), 856–866. doi: 10.1002/jcb.29330

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, L., Xu, L., Chen, L., Fu, Q., Liu, Z., Chang, Y., et al. (2017). Tumor-infiltrating neutrophils predict benefit from adjuvant chemotherapy in patients with muscle invasive bladder cancer. Oncoimmunology 6 (4), 00–00. doi: 10.1080/2162402X.2017.1293211

CrossRef Full Text | Google Scholar