Abstract
Immune checkpoint inhibitor (ICI) treatment has been used to treat advanced urothelial cancer. Molecular markers might improve risk stratification and prediction of ICI benefit for urothelial cancer patients. We analyzed 406 cases of bladder urothelial cancer from The Cancer Genome Atlas (TCGA) data set and identified 161 messenger RNAs (mRNAs) as differentially expressed immunity genes (DEIGs). Using the LASSO Cox regression model, an eight-mRNA-based risk signature was built. We validated the prognostic and predictive accuracy of this immune-related risk signature in 348 metastatic urothelial cancer (mUC) samples treated with anti-PD-L1 (atezolizumab) from IMvigor210. We built an immune-related risk signature based on the eight mRNAs: ANXA1, IL22, IL9R, KLRK1, LRP1, NRG3, SEMA6D, and STAP2. The eight-mRNA-based risk signature successfully categorizes patients into high-risk and low-risk groups. Overall survival was significantly different between these groups, regardless if the initial TCGA training set, the internal TCGA testing set, all TCGA set, or the ICI treatment set. The hazard ratio (HR) of the high-risk group to the low-risk group was 3.65 (p < 0.0001), 2.56 (p < 0.0001), 3.36 (p < 0.0001), and 2.42 (p = 0.0009). The risk signature was an independent prognostic factor for prediction survival. Moreover, the risk signature was related to immunity characteristics. In different tumor mutational burden (TMB) subgroups, it successfully categorizes patients into high-risk and low-risk groups, with significant differences of clinical outcome. Our eight-mRNA-based risk signature is a stable biomarker for urothelial cancer and might be able to predict which patients benefit from ICI treatment. It might play a role in precision individualized immunotherapy.
Introduction
Bladder cancer is a common tumor of the urinary system; 90% of the pathological types are urothelial cancer (UC). For bladder cancer patients with local progression or distant metastasis, cisplatin combined with gemcitabine is the first choice. However, the effect is not satisfactory. The median survival time of patients is only 15 months, and the 5 years survival rate is difficult to reach 15% (Griffiths et al., 2011). In recent years, with the rapid development of tumor immune checkpoint inhibitors (ICIs), especially the rapid development of programmed cell death molecule 1 (PD-1)/programmed cell death molecule ligand 1 (PD-L1) inhibitors, the treatment of bladder cancer patients has brought new options. In patients with metastatic UC (mUC) who cannot receive cisplatin chemotherapy and PD-L1 positive, ICIs can already be used as a first-line treatment (Nadal and Bellmunt, 2019).
Although it is proved the efficacy of PD-1/PD-L1 inhibitors is better than that of traditional platinum-based chemotherapy (Bellmunt et al., 2017; Sidaway, 2017), studies have confirmed that only about 20% of solid tumor patients can benefit from the treatment (Braun et al., 2016). Therefore, it is becoming more and more important to identify and verify biomarkers that can accurately predict ICI treatment efficacy. There have been some clinical studies exploring the corresponding biomarkers, such as PD-L1 expression (Powles et al., 2014), CD8+ T cell (Ghatalia and Plimack, 2019), tumor mutational burden (TMB) (Yarchoan et al., 2017), and microsatellite instability (MSI) (Dudley et al., 2016). However, these biomarkers have shortcomings in clinical application. In addition, multi-factor joint prediction may be able to provide better prediction results (Mazzaschi et al., 2020).
In this study, we analyzed the messenger RNA (mRNA) transcriptome data of 406 bladder cancer patients from The Cancer Genome Atlas (TCGA) data combined with immune-related genes to establish an immune-related risk signature, and we verified the 348 mUC patients receiving anti-PD-L1 therapy from IMvigor210 study. We emphasized the strong predictive power of the risk score in selecting patients with good response to atezolizumab and verified its role in ICI treatment.
Materials and Methods
Clinical Cohorts and Data Sets
The gene expression sequence matrix and clinical characteristics of 406 bladder cancer patients can be downloaded from TCGA data set1. The immune gene sets come from the ImmPort2 and InnateDB3 data sets (Breuer et al., 2013; Bhattacharya et al., 2014). Clinical information and gene transcription information of 348 patients with mUC who received ICI treatment are downloaded from4 (Mariathasan et al., 2018). The infiltration of 22 immune cells was downloaded from the TIMER database5 and Dongqiang Zeng’s research (Li T. et al., 2020; Zeng et al., 2020).
Bioinformatic Analysis
R package DESeq2 was used for gene expression differential analysis, and R package clusterProfiler for Gene Ontology/Kyoto Encyclopedia of Genes and Genomes/Gene Set Enrichment Analysis (GO/KEGG/GSEA) function enrichment analysis and visualization (Subramanian et al., 2005; Love et al., 2014). The ggstatsplot package was used to evaluate the relationship between risk score, TMB, and immunophenotype.
Establishment and Evaluation of Risk Prediction Model
We randomly divide the samples in TCGA cohort into training/validation (3:1) groups to identify and evaluate predictors. The “glmnet” R package was used for LASSO analysis, and 14 immune-related genes were identified (Ternes et al., 2016). Then we conducted multiple Cox regression analysis to establish an eight-mRNA-based risk prediction model. Use the formula to generate the risk score for each patient: risk score = EXP 1 ∗ β 1 EXP2 ∗ β 2 … EXPÑ ∗ β Ñ, where “EXP” represents the expression level of key genes and β is the corresponding regression coefficient. The “timeROC” package was used to establish the receiver operating characteristic (ROC) curve and verify the area value under it (AUC). Draw a Kaplan–Meier curve to show the association of risk scores and potential prognostic genes with patient survival.
Statistical Analysis
The Kaplan–Meier method was used to analyze the correlation between relate risk factors with patient survival. Statistical tests were performed using R software, version 3.6. 2 (R Foundation for Statistical Computing; Vienna, Austria). Values of p < 0.05 were considered statistically significant.
Results
Describe Immune-Related Gene Features and Construct Immune-Related Prediction Models
Based on the Cibersort algorithm, the infiltration levels of 22 immune cells in 406 bladder cancer patients in TCGA data set were evaluated. The fuzzy clustering algorithm divided the samples into two categories: 210 samples were clustered into cluster 1 and 196 samples were clustered into cluster 2 (Figure 1). Cluster 1 had high-level immune characteristics (high-immunity group), and cluster 2 had low-level immune characteristics (low-immunity group). Missing clinical data are shown as blank on the top of the heatmap.
FIGURE 1
We were trying to study immune-related genes stratified by immune phenotype and their prognostic potential and to establish a good immunological prediction model, which can make more accurate individual risk stratification and prognosis prediction for UC patients. After preliminary screening through single-factor regression analysis and difference analysis based on immunity clustering, a total of 1,976 genes were identified as differentially expressed genes (DEGs). All DEG expression levels are shown in Figure 2A, and the log2 and p-values of all DEGs are shown in Figure 2B. Subsequently, 161 genes were identified as differentially expressed immunity genes (DEIGs) based on ImmPort and InnateDB databases (Figure 2C). In order to further explore the prognostic significance of DEIGs, 161 important genes were used for multiple LASSO regression and multiple Cox regression analysis (Figure 2D), and eight key DEIGs were identified. Distribution of each DEIG in TCGA set and ICI treatment set is shown in Supplementary Figure 1. The Kaplan–Meier curve of each DEIG in TCGA set is shown in Supplementary Figure 2 [ANXA1: p < 0.0001, hazard ratio (HR) = 2.02; IL22: p < 0.0001, HR = 0.50; IL9R: p < 0.0001, HR = 0.71; KLRK1: p < 0.0001, HR = 0.54; LRP1: p = 0.0051, HR = 1.96; NRG3: p = 0.00065, HR = 1.40; SEMA6D: p = 0.01, HR = 1.65; STAP2: p = 0.00029, HR = 0.68]. Finally, according to the relative coefficient in the multiple regression analysis, the risk score was calculated according to the following formula: (0.144493498 ∗ ANXA1) + (−0.387235675 ∗ IL22) + (−0.103863619 ∗ IL9R) + (−0.291034924 ∗ KLRK1) + (0.137267967 ∗ LRP1) + (0.090089369 ∗ NRG3) + (0.091881101 ∗ SEMA6D) + (−0.138073124 ∗ STAP2) (Table 1).
FIGURE 2
TABLE 1
| coef | HR | HR.95L | p-value | |
| ANXA1 | 0.1444935 | 1.16 | 1.04–1.29 | 0.0085 |
| IL22 | –0.3872357 | 0.68 | 0.47–0.98 | 0.0365 |
| IL9R | –0.1038636 | 0.9 | 0.82–0.99 | 0.0325 |
| KLRK1 | –0.2910349 | 0.75 | 0.66–0.85 | <0.0001 |
| LRP1 | 0.13726797 | 1.15 | 0.96–1.37 | 0.128 |
| NRG3 | 0.09008937 | 1.09 | 0.99–1.21 | 0.0771 |
| SEMA6D | 0.0918811 | 1.1 | 0.97–1.24 | 0.1539 |
| STAP2 | –0.1380731 | 0.87 | 0.74–1.02 | 0.0878 |
Characteristics of differentially expressed immunity genes (DEIGs) in the risk signature.
Survival Analysis, Prognostic Value, and Immune Infiltration Verification of the Risk Signature
The prognostic value of eight DEIG signatures was further evaluated in three verification sets (TCGA test set, all TCGA set, and independent ICI treatment set). We calculated the risk score of each patient using the same formula, and we divided them into high-risk and low-risk groups by 1.54 as a cutoff. Consistent with the results of TCGA training set, the prognosis of high-risk patients in the three validation sets was worse than that of patients in the low-risk group (Figures 3A–D, left; TCGA training set: p < 0.0001, HR = 3.65; TCGA test set: p < 0.0001, HR = 2.56; all TCGA set: p < 0.0001, HR = 3.36; ICI treatment set: p = 0.0009, HR = 2.42). The results of the time-dependent ROC curve analysis verified the predictive value of the established risk model (Figures 3A–D, right), suggesting that the prognosis prediction for 3–5 years was more robust. The univariate and multivariate Cox analyses of TCGA set showed that the risk signature can be used as an independent prognostic factor (Supplementary Table 1). We examined the correlation between risk signature and the bladder cancer immune microenvironment. Both TCGA and ICI data sets showed a relatively consistent trend of infiltration. In terms of immune cell infiltration, such as T.cell.CD8.positive, T.cell.CD4.activated, and Macrophage.M0, the infiltration trend was the same in the two data sets. The differences were statistically significant (Figures 3E,F).
FIGURE 3
Verification and Comparison of the Correlation Between the Risk Signature and Immune Checkpoint Inhibitor Treatment Efficacy
There was a significant correlation between the risk signature and TMB, regardless if in TCGA or ICI treatment set (Figures 4A,B). At the same time, the immune subtypes classified according to CD8 cell infiltration (desert, excluded, and inflamed) were also obviously related to the risk signature (Figure 4C). TCGA type II subgroup had the lowest risk score (Figure 4D). This is consistent with the previous results (Mariathasan et al., 2018).
FIGURE 4
We divided TMB into high-risk and low-risk groups and then subdivided them into subgroups based on the risk scores level (TCGA set: p < 0.0001, HR = 0.49; ICI treatment set: p = 0.0006, HR = 0.54). The results suggested that even in the TMB subgroup, the risk signature still remained its prognostic ability, regardless if in TCGA or ICI treatment set (Figure 5, TCGA high-TMB group: p < 0.0001, HR = 4.26; TCGA low-TMB group: p < 0.0001, HR = 2.83; ICI treatment high-TMB group: p = 0.0034, HR = 10.64; ICI treatment low-TMB group: p = 0.031, HR = 2.56). Multivariate risk regression also confirmed this result (Table 2).
FIGURE 5
TABLE 2
| HR | 95%CI | p-value | |
| Risk score | 4.83 | 2.14−10.93 | 0.0002 |
| TMB | 0.96 | 0.94−0.99 | 0.0012 |
Multivariate Cox regression of risk scores and tumor mutational burden (TMB) in immune checkpoint inhibitor (ICI) treatment set.
Gene Ontology/Kyoto Encyclopedia of Genes and Genomes/Gene Set Enrichment Analysis
In order to further explore the molecular mechanisms related to risk scores, we divided TCGA cohort patients into high-risk and low-risk groups. The results of GO and KEGG suggested that the risk signature was related to the extracellular matrix and energy metabolism changes in the tumor microenvironment (Figures 6A,B). GSEA results suggested that the high-risk group is positively correlated with steroid metabolism and YP450 metabolism, while pathways such as cytokine interaction and immune response are positively correlated with low-scoring risks (Figure 6C).
FIGURE 6
Discussion
As a new treatment method, ICI has initially proven efficacy and safety in the treatment of UC. Unfortunately, not all patients with cancer respond to ICI treatment. In previous studies, the use of risk signature derived from the gene transcriptome to monitor the immune status of tumors and guide individualized treatment has proven to be meaningful (Li et al., 2017; Cristescu et al., 2018; Wang et al., 2019). Therefore, the development of meaningful genetic markers to monitor the immune status of patients not only can monitor the prognosis of patients but also can screen out potential ICI responding patients, avoiding the waste of medical resources and overtreatment. In this study, we validated immune-related gene risk model based on eight DEIGs, which proved to be a reliable indicator of favorable ICI efficacy and can identify bladder cancer patients with poor prognosis. Our results showed that the prognosis is worse if the risk score is higher. At the same time, the time-dependent ROC curve results suggested that the 3–5 years’ prognosis prediction for UC patients was more robust. In terms of tumor immune cell infiltration, whether in TCGA group or the ICI treatment group, the differences in T.cell.CD8.positive, T.cell.CD4.activated, and Macrophage.M0 were statistically significant. Patients in the high-risk group had significantly lower representation of T.cell.CD8.positive and T.cell.CD4.activated and significantly higher abundance of Macrophage.M0. This is also consistent with the results of other studies (Li W. et al., 2020; Lin et al., 2020).
PD-L1 is currently the most mature and in-depth biomarker, but the results obtained in different ICI studies are not consistent. It may be due to the different monitoring methods and positive standards between different platforms, and the evaluation is subjective. At present, due to technical requirements, it is difficult to apply TMB to routine clinical practice (Fenizia et al., 2018). In addition, TMB has not yet proven its predictive or prognostic value for overall survival (Addeo et al., 2019). According to the immune microenvironment, most solid tumors can be divided into three different immunological phenotypes: immune inflamed, immune excluded, or immune desert (Chen and Mellman, 2013; Hegde et al., 2016). Studies have shown that immune inflamed subtypes have the best response to ICI treatments (such as anti-PD-1 and anti-CTLA-4) (Mariathasan et al., 2018; Galon and Bruni, 2019). The immune-related risk signature we established was significantly correlated with TMB and immunophenotype. A lower risk score means a higher TMB, a better response, and a better prognosis. At last, the K-M curve of TMB subtype showed that the risk signature was able to be used as a stable predictor. In order to further clarify the mechanism of immune risk score, we subsequently used TCGA data set to conduct GO, KEGG, and GSEA. The results of GO and KEGG show that the risk score is related to the energy metabolism and synthesis of the tumor microenvironment, and the formation and activation of extracellular matrix. In the analysis of GSEA results, we can see that the synthesis and metabolism of steroids are positively correlated with high-scoring risks, while pathways such as cytokine interaction and immune response are positively correlated with low-scoring risks. This is consistent with the results of other studies. Zeng et al. (2020) found that the defect of M1 macrophage function is related to poor prognosis of UC immunotherapy, and it is also positively related to steroid synthesis and metabolism.
Among the eight DEIGs, there are few studies in bladder cancer, but some of their interactions with immunity have been explored and verified in other researches. ANXA1 can enhance the function of regulatory T cells (Tregs) and reduce the survival rate of patients with breast cancer (Bai et al., 2020). IL-22 producing T cells in colorectal cancer enhance T cell function by recruiting neutrophils, thereby enhancing immune response (Tosti et al., 2020). Th9 cells promote the expansion of CD8+ T cells in an IL-9R-dependent manner in colorectal cancer (Wang et al., 2020). CIK cells can target lung cancer cells expressing NKG2D/KLRK1 ligand, and the killing effect can be partially blocked by NKG2D/KLRK1 ligand inhibitors (Yin et al., 2017). The correlation between LRP1 mRNA expression and patient survival was observed in bladder urothelial carcinoma. At the same time, the LRP1 protein can regulate the immune function by regulating the movement and adhesion of T cells (Gonias et al., 2017; Panezai et al., 2017). STAP2 maintains the cytotoxicity of functional memory CD8+ T cells by controlling cytokine signaling inhibitor 3 (Muraoka et al., 2017). SEMA6D act as a modulator in the late stage of the primary immune response (O’Connor et al., 2008).
Although the risk signature based on eight DEIGs embodies a powerful predictive function in selecting patients with good response to atezolizumab, its accuracy and effectiveness should be further verified in a prospective cohort study receiving immunotherapy. In addition, the molecular mechanism of the protein encoded by DEIGs in UC still needs to be explored in vitro and in vivo.
The risk signature is a stable biomarker that can be used to predict immunotherapy efficacy and immunophenotype determination, and it can be used as a supplement to TMB.
Statements
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: The gene expression sequence matrix and clinical characteristics of 406 bladder cancer patients was downloaded from the TCGA data set (https://portal.gdc.cancer.gov). The immune gene sets were downloaded from the Immport (https://s3.immport.org/release/genelists/GeneList.txt?download=true) and InnateDB (https://www.innatedb.com/download/innatedb_curated_genes.xls) data sets. Clinical information and gene transcription information of 348 patients with mUC who received ICI treatment was downloaded from http://research- pub.gene.com/IMvigor210CoreBiologies. The infiltration of 22 immune cells was downloaded from the TIMER database (http:// timer.cistrome.org/infiltration_estimation_for_tcga.csv.gz) and Dongqiang Zeng’s research (https://www.thno.org/v10/p7002/thnov10p7002s2.xlsx).
Author contributions
JLu and WC designed the study. PL, SH, YY, YT, LT, CG, and JW obtained and assembled data. PL, SH, JLi, JW, YH, ZL, and JC analyzed and interpreted the data. PL, JX, MZ, XC, and WC wrote the report. All authors approved the final version. JC, JLu, and WC are the guarantors.
Funding
This study was supported by grants from the National Natural Science Foundation of China (Award Nos. 81725016, 81872094, 81772718, 81602219, 81972376, and 81902576) and the Guangdong Provincial Science and Technology Foundation of China (Award Nos. 2017B020227004 and 2017A030313538).
Acknowledgments
We thank TCGA, ImmPort, InnateDB, TIMER, study IMvigor210, and Dongqiang Zeng’s research for their efforts and for providing the data.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.646982/full#supplementary-material
Supplementary Figure 1Distribution of each DEIG in TCGA set and ICI treatment set.
Supplementary Figure 2Kaplan–Meier curve of each DEIG in TCGA set.
Supplementary Table 1Univariate and multivariate Cox regression analysis of clinical pathologic features in TCGA set.
References
1
AddeoA.BannaG. L.WeissG. J. (2019). Tumor mutation burden-from hopes to doubts.JAMA Oncol.5934–935. 10.1001/jamaoncol.2019.0626
2
BaiF.ZhangP.FuY.ChenH.ZhangM.HuangQ.et al (2020). Targeting ANXA1 abrogates Treg-mediated immune suppression in triple-negative breast cancer.J. Immunother. Cancer8:e000169. 10.1136/jitc-2019-000169
3
BellmuntJ.de WitR.VaughnD. J.FradetY.LeeJ. L.FongL.et al (2017). Pembrolizumab as second-line therapy for advanced urothelial carcinoma.N. Engl. J. Med.3761015–1026. 10.1056/NEJMoa1613683
4
BhattacharyaS.AndorfS.GomesL.DunnP.SchaeferH.PontiusJ.et al (2014). ImmPort: disseminating data to the public for the future of immunology.Immunol. Res.58234–239. 10.1007/s12026-014-8516-1
5
BraunD. A.BurkeK. P.Van AllenE. M. (2016). Genomic approaches to understanding response and resistance to immunotherapy.Clin. Cancer Res.225642–5650. 10.1158/1078-0432.CCR-16-0066
6
BreuerK.ForoushaniA. K.LairdM. R.ChenC.SribnaiaA.LoR.et al (2013). InnateDB: systems biology of innate immunity and beyond–recent updates and continuing curation.Nucleic Acids Res.41D1228–D1233. 10.1093/nar/gks1147
7
ChenD. S.MellmanI. (2013). Oncology meets immunology: the cancer-immunity cycle.Immunity391–10. 10.1016/j.immuni.2013.07.012
8
CristescuR.MoggR.AyersM.AlbrightA.MurphyE.YearleyJ.et al (2018). Pan-tumor genomic biomarkers for PD-1 checkpoint blockade-based immunotherapy.Science362eaar3593. 10.1126/science.aar3593
9
DudleyJ. C.LinM. T.LeD. T.EshlemanJ. R. (2016). Microsatellite instability as a biomarker for PD-1 blockade.Clin. Cancer Res.22813–820. 10.1158/1078-0432.CCR-15-1678
10
FeniziaF.PasqualeR.RomaC.BergantinoF.IannacconeA.NormannoN. (2018). Measuring tumor mutation burden in non-small cell lung cancer: tissue versus liquid biopsy.Transl. Lung Cancer Res.7668–677. 10.21037/tlcr.2018.09.23
11
GalonJ.BruniD. (2019). Approaches to treat immune hot, altered and cold tumours with combination immunotherapies.Nat. Rev. Drug Discov.18197–218. 10.1038/s41573-018-0007-y
12
GhataliaP.PlimackE. (2019). Biomarkers for neoadjuvant checkpoint blockade response in urothelial cancer.Nat. Med.251650–1651. 10.1038/s41591-019-0645-6
13
GoniasS. L.Karimi-MostowfiN.MurrayS. S.MantuanoE.GilderA. S. (2017). Expression of LDL receptor-related proteins (LRPs) in common solid malignancies correlates with patient survival.PLoS One12:e0186649. 10.1371/journal.pone.0186649
14
GriffithsG.HallR.SylvesterR.RaghavanD.ParmarM. K. (2011). International phase III trial assessing neoadjuvant cisplatin, methotrexate, and vinblastine chemotherapy for muscle-invasive bladder cancer: long-term results of the BA06 30894 trial.J. Clin. Oncol.292171–2177. 10.1200/JCO.2010.32.3139
15
HegdeP. S.KaranikasV.EversS. (2016). The where, the when, and the how of immune monitoring for cancer immunotherapies in the era of checkpoint inhibition.Clin. Cancer. Res.221865–1874. 10.1158/1078-0432.CCR-15-1507
16
LiB.CuiY.DiehnM.LiR. (2017). Development and validation of an individualized immune prognostic signature in early-stage nonsquamous non-small cell lung cancer.JAMA Oncol.31529–1537. 10.1001/jamaoncol.2017.1609
17
LiT.FuJ.ZengZ.CohenD.LiJ.ChenQ.et al (2020). TIMER2.0 for analysis of tumor-infiltrating immune cells.Nucleic Acids Res.48W509–W514. 10.1093/nar/gkaa407
18
LiW.ZengJ.LuoB.MaoY.LiangY.ZhaoW.et al (2020). [High expression of activated CD4(+) memory T cells and CD8(+) T cells and low expression of M0 macrophage are associated with better clinical prognosis in bladder cancer patients].Xi Bao Yu Fen Zi Mian Yi Xue Za Zhi3697–103.
19
LinJ.YangJ.XuX.WangY.YuM.ZhuY. (2020). A robust 11-genes prognostic model can predict overall survival in bladder cancer patients based on five cohorts.Cancer Cell Int.20:402. 10.1186/s12935-020-01491-6
20
LoveM. I.HuberW.AndersS. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2.Genome Biol.15:550. 10.1186/s13059-014-0550-8
21
MariathasanS.TurleyS. J.NicklesD.CastiglioniA.YuenK.WangY.et al (2018). TGFbeta attenuates tumour response to PD-L1 blockade by contributing to exclusion of T cells.Nature554544–548. 10.1038/nature25501
22
MazzaschiG.MinariR.ZeccaA.CavazzoniA.FerriV.MoriC.et al (2020). Soluble PD-L1 and circulating CD8+PD-1+ and NK cells enclose a prognostic and predictive immune effector score in immunotherapy treated NSCLC patients.Lung Cancer1481–11. 10.1016/j.lungcan.2020.07.028
23
MuraokaD.SeoN.HayashiT.Hyuga-AmaikeC.OkamoriK.TawaraI.et al (2017). Signal-transducing adaptor protein-2 promotes generation of functional long-term memory CD8+ T cells by preventing terminal effector differentiation.Oncotarget830766–30780. 10.18632/oncotarget.15403
24
NadalR.BellmuntJ. (2019). Management of metastatic bladder cancer.Cancer Treat. Rev.7610–21. 10.1016/j.ctrv.2019.04.002
25
O’ConnorB. P.EunS. Y.YeZ.ZozulyaA. L.LichJ. D.MooreC. B.et al (2008). Semaphorin 6D regulates the late phase of CD4+ T cell primary immune responses.Proc. Natl. Acad. Sci. U.S.A.10513015–13020. 10.1073/pnas.0803386105
26
PanezaiJ.BergdahlE.SundqvistK. G. (2017). T-cell regulation through a basic suppressive mechanism targeting low-density lipoprotein receptor-related protein 1.Immunology152308–327. 10.1111/imm.12770
27
PowlesT.EderJ. P.FineG. D.BraitehF. S.LoriotY.CruzC.et al (2014). MPDL3280A (anti-PD-L1) treatment leads to clinical activity in metastatic bladder cancer.Nature515558–562. 10.1038/nature13904
28
SidawayP. (2017). Bladder cancer: pembrolizumab is superior to chemotherapy.Nat. Rev. Urol.14:261. 10.1038/nrurol.2017.38
29
SubramanianA.TamayoP.MoothaV. K.MukherjeeS.EbertB. L.GilletteM. A.et al (2005). Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles.Proc. Natl. Acad. Sci. U.S.A.10215545–15550. 10.1073/pnas.0506580102
30
TernesN.RotoloF.MichielsS. (2016). Empirical extensions of the lasso penalty to reduce the false discovery rate in high-dimensional Cox regression models.Stat. Med.352561–2573. 10.1002/sim.6927
31
TostiN.CremonesiE.GovernaV.BassoC.KancherlaV.Coto-LlerenaM.et al (2020). Infiltration by interleukin-22 producing T cells promotes neutrophil recruitment and predicts favorable clinical outcome in human colorectal cancer.Cancer Immunol. Res.81452–1462. 10.1158/2326-6066.CIR-19-0934
32
WangC.LuY.ChenL.GaoT.YangQ.ZhuC.et al (2020). Th9 cells are subjected to PD-1/PD-L1-mediated inhibition and are capable of promoting CD8 T cell expansion through IL-9R in colorectal cancer.Int. Immunopharmacol.78:106019. 10.1016/j.intimp.2019.106019
33
WangZ.SongQ.YangZ.ChenJ.ShangJ.JuW. (2019). Construction of immune-related risk signature for renal papillary cell carcinoma.Cancer Med.8289–304. 10.1002/cam4.1905
34
YarchoanM.HopkinsA.JaffeeE. M. (2017). Tumor mutational burden and response rate to PD-1 inhibition.N. Engl. J. Med.3772500–2501. 10.1056/NEJMc1713444
35
YinX.LuX.XiuwenZ.MinZ.XiaoR.MaoZ.et al (2017). Role of NKG2D in cytokine-induced killer cells against lung cancer.Oncol. Lett.133139–3143. 10.3892/ol.2017.5800
36
ZengD.YeZ.WuJ.ZhouR.FanX.WangG.et al (2020). Macrophage correlates with immunophenotype and predicts anti-PD-L1 response of urothelial cancer.Theranostics107002–7014. 10.7150/thno.46176
Summary
Keywords
immune-related risk signature, immunity gene, immune checkpoint inhibitor, urothelial cancer, tumor microenvironment
Citation
Li P, Hao S, Ye Y, Wei J, Tang Y, Tan L, Liao Z, Zhang M, Li J, Gui C, Xiao J, Huang Y, Chen X, Cao J, Luo J and Chen W (2021) Identification of an Immune-Related Risk Signature Correlates With Immunophenotype and Predicts Anti-PD-L1 Efficacy of Urothelial Cancer. Front. Cell Dev. Biol. 9:646982. doi: 10.3389/fcell.2021.646982
Received
28 December 2020
Accepted
03 February 2021
Published
18 March 2021
Volume
9 - 2021
Edited by
Wei Zhao, Chengdu Medical College, China
Reviewed by
An Zhao, University of Chinese Academy of Sciences, China; Qing Zhang, Nanjing Drum Tower Hospital, China
Updates
Copyright
© 2021 Li, Hao, Ye, Wei, Tang, Tan, Liao, Zhang, Li, Gui, Xiao, Huang, Chen, Cao, Luo and Chen.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jiazheng Cao, 530694776@qq.comJunhang Luo, luojunh@mail.sysu.edu.cnWei Chen, chenw3@mail.sysu.edu.cn
†These authors have contributed equally to this work
This article was submitted to Molecular and Cellular Oncology, a section of the journal Frontiers in Cell and Developmental Biology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.