- Research article
- Open Access
Genomic stratification based on microenvironment immune types and PD-L1 for tailoring therapeutic strategies in bladder cancer
BMC Cancer volume 21, Article number: 646 (2021)
The tumour microenvironment (TME) not only plays a role during tumour progression and metastasis but also profoundly influences treatment efficacy. Environment-mediated drug resistance is a result of crosstalk between tumour cells and stroma. The presence of a “stromal exhaustion” response is suggested by the T cell exhaustion signature and PD-L1 expression. The prognostic role of PD-L1 in bladder cancer has been investigated in previous studies, but the results remain inconclusive. For a more comprehensive study, we discuss potential strategies to improve effectiveness in patients with various TME statuses and PD-L1 expression levels.
We estimated the prognostic role of PD-L1 using immunohistochemistry and identified four immune subtypes according to the type of stromal immune modulation and PD-L1 expression status.
Patients in the PD-L1-low-exhausted group had the worst prognosis and showed the worst antigen-presenting cell (APC) immunosuppression status. The PD-L1-low-exhausted group showed the highest amount of infiltration by macrophage M2 cells, naïve B cells and resting mast cells. The TMB and the effectiveness of anti-PD-1 treatment were significantly increased in the PD-L1-high expression groups compared with the PD-L1-low expression groups. In the PD-L1-high groups, patients who underwent chemotherapy exhibited better overall survival rates than patients who did not undergo chemotherapy, whereas there was no significant difference in the PD-L1-low groups. We performed gene set enrichment analysis (GSEA) to explore the critical pathways that were active in the PD-L1-low-exhausted group, including the myogenesis, epithelial-mesenchymal transition and adipogenesis pathways. Copy number variations (CNVs) were related to the expression levels of differentially expressed genes upregulated in the PD-L1-low-exhausted group, including LCNL1, FBP1 and RASL11B. In addition, RASL11B played a role in predicting overall survival according to The Cancer Genome Atlas data, and this finding was validated in the PD-L1-low-exhausted group in the Gene Expression Omnibus database (GEO).
The immune environment of tumours plays an important role in the therapeutic response rate, and defining the immune groups plays a critical role in predicting disease outcome and strategy effectiveness.
Bladder cancer is the fifth most common cancer and is diagnosed in greater than 430,000 patients every year worldwide . Bladder cancer shows a sex bias with approximately 75% of patients being male and ranks fourth in incidence and eighth in mortality among cancers in males [2, 3]. Approximately 20% of bladder cancer patients are diagnosed with muscle-invasive bladder cancer (MIBC). Treatment options include radical cystectomy and chemotherapy, but the prognosis of patients with metastasis or nonorgan confinement can be poor .
Recently, immune checkpoint blockade has changed the therapy landscape in bladder cancer. Currently, anti-PD-1/PD-L1 strategies have been recommended as the first-line treatment for cisplatin-intolerant PD-1-positive patients with local late/metastatic bladder cancer (2020 NCCN v2). After several decades with no significant advanced therapeutic measures, some clinical trials have demonstrated that durable effects can be achieved in 15–25% cisplatin-resistant metastatic urothelial cancer (UC) patients treated with PD-1/PD-L1 antibodies [5,6,7,8,9,10]. Given that only a few patients benefit from PD-1/PD-L1 antibodies, there remains a critical need to explore mechanisms of resistance .
The tumour microenvironment (TME) comprises immune cells, endothelial cells and cytokines, growth factors, etc. The TME not only plays a role during tumour progression and metastasis but also has profound influences on the efficacy of treatment. Environment-mediated drug resistance is a result of crosstalk between tumour cells and stroma . “Hot tumours” filled with immune cells are often considered to be more sensitive to immunotherapy than “cold tumours”, which have fewer immune cells. The role of tumour-infiltrating lymphocytes (TILs) and cytokines in immunotherapy has been demonstrated . Nevertheless, the effect of the TME on the response to immunotherapy remains controversial in urothelial carcinoma . Efforts are currently focused on investigating the relationships between molecular subtype and immunotherapy response. Recently, the TCGA molecular subtype (k = 5) classification was applied to the IMvigor 210 atezolizumab trial data. The results indicated that the neuroendocrine-like (NE-like) and luminal subtypes were both associated with high response rates . However, the opposite pattern of NE-like tumour response to immunotherapy was observed with the Decipher classifier . Therefore, further characterization of the immune profile to identify subtypes in bladder cancer that may be relevant for therapy choice is essential.
In this study, to identify common immune subgroups, evaluate how the immune components of the TME regulate therapeutic effectiveness and explore potential strategies to improve the effectiveness of patients with various immune statuses in bladder cancer, we divided patients into four subgroups based on a combination of PD-L1 expression and the type of immune modulation, and then we performed a series of analyses. First, survival analysis was performed to assess the associations between immune subgroups and patient prognosis. We then compiled published tumour immune expression signatures and obtained scores for all samples. The leukocyte fraction, predicted neoantigens, proliferation rate, B cell receptor (BCR) expression, cancer testis antigen (CTA) expression, immunomodulator (IM) expression and DNA damage were characterized in the samples. In addition, we next estimated the mean fractions of immune cells in samples to assess the relationship between immunogenicity and immune infiltration. Finally, we sought to identify the underlying molecular mechanism governing the various immune responses across subgroups. GSEA at the RNA level and CNV analysis at the DNA level were performed to identify the key pathways and key genes in the PD-L1-low-exhausted group and provide potential strategies for patients with various immune statuses.
Within the TCGA database, we assessed 363 bladder cancer samples with mRNA expression data, CNV data and clinical data. We limited our research to data for samples with the “primary tumour” sample type (n = 357) and excluded samples with the “solid normal tissue” sample type (n = 6). For the tumour mutation burden (TMB), we obtained data from a previously published signature . The normalized level three RNA sequencing data, CNV data and clinical data were downloaded from https://genome-cancer.ucsc.edu/.
A total of 106 biopsy samples from patients with primary bladder cancer were collected at the First Hospital of China Medical University between 2009 and 2012. All the patients accepted radical surgery. The primary antibody was PD-L1 (mouse anti-human PD-L1 monoclonal, dilution 1:100, Abcam). The percentage of positive cells and staining intensity were estimated. Consistent with previous reports, staining in greater than 5% of the tumour cells was considered positive.
Clustering & definition of groups
In previous studies, “activated” stroma was characterized by an activated fibroblast state, a set of genes associated with macrophages, and the expression of other genes that point to its role in tumour promotion. The “normal” stromal factor may describe a “good” version of the stroma, whereas the “activated” stromal factor may describe an activated inflammatory stromal response to be responsible for disease progression . Based on the expression of the “activated” stroma gene signature , we divided patients into “activated” (“normal” stroma) and “exhausted” (“activated” stroma) subgroups using the consensus clustering algorithm (k = 2). Consensus clustering of the samples using “activated” stroma genes was refined using a random forest (RF) algorithm developed under the ConsensusClusterPlus R package. To investigate the clinicopathological significance of PD-L1 expression and the type of immune modulation, the cohort was divided into four subgroups: (1) “PD-L1-high-activated”, (2) “PD-L1-high-exhausted”, (3) “PD-L1-low-activated” and (4) “PD-L1-low-exhausted”.
Gene set enrichment analysis (GSEA)
GSEA was performed using a Java GSEA desktop application (www.broad.mit.edu/gsea/). The hallmark gene set “h.all.v7.0.symbols.gmt” was used in this study. A nominal P-value < 0.05 and a false discovery rate (FDR) < 0.25 were used to determine significance.
GEO validation dataset
Validation data (GES48277) were downloaded from the GEO database (www.ncbi.nlm.nih.gov/geo), and we limited our research to bladder tumour biopsy sample data. We divided patients into four subgroups according to PD-L1 expression and the type of immune modulation. To validate the role of RASL11B, we evaluated its expression in the subgroups.
We identified differentially expressed genes (DEGs) based on the criteria of a false discovery rate < 0.05 and a log fold-change > 1.0, and we performed copy number variation (CNV) analysis with the criterion of a P-value < 0.05 and estimated the infiltration of immune cells using R 3.6.1. We adopted an immunophenoscore (IPS) from The Cancer Immunome Atlas, which was based on the expression of the representative genes or gene sets comprising four categories: MHC molecules, immunomodulators, effector cells (CD8+ T cells and CD4+ T cells), and suppressor cells (Tregs and MDSCs). The IPS was indicated to be a superior predictor of response to anti-PD-1 antibodies in two independent cohorts . We adopted chemosensitivity  and immune expression signatures from a previously published signature . The χ2 test was employed to compare the clinical and pathological characteristics of the groups. Data distributions were tested using the Kolmogorov-Smirnov test. Parametric or nonparametric tests were employed according to the results. Immune expression signatures, DNA damage signatures and gene expression of IMs were compared among subgroups using ANOVA for continuous variables. Kaplan–Meier analysis for the overall survival of patients in TCGA dataset was performed using the GEPIA web server . Two-tailed P-values less than 0.05 were considered statistically significant. Statistical analyses were performed using IBM SPSS Statistics software version 24 and GraphPad Prism software 7.00.
Classification of patients into subgroups according to PD-L1 expression and the type of immune modulation
Given that immune responses can both prevent and promote cancer progression , the subtype of immune modulation should be considered. Accordingly, an “exhausted” response was suggested to predict worse survival in pancreatic ductal adenocarcinoma , whereas there was no significant difference between an “exhausted” (162/363) or an “active” response (201/363) in our study (P = 0.09) (Fig. 1a, b). As the presence of an “exhaustion” response was suggested by the T cell exhaustion signature , we analysed the PD-L1 mRNA expression level in two subgroups. The exhausted group exhibited increased PD-L1 mRNA expression than the active group (1.54 ± 1.20 vs 1.22 ± 1.16, P = 0.012) (Fig. 1c). The prognostic role of PD-L1 in bladder cancer has been investigated in previous studies, but the results remain inconclusive. For a more comprehensive study, we estimated the prognostic role of PD-L1 using immunohistochemistry. Patient biopsy samples (n = 106) were assessed for PD-L1 expression in tumour tissues (Fig. 1d). The results suggested that PD-L1 overexpression could predict worse survival outcomes in bladder cancer (P = 0.045) (Fig. 1e). Next, we divided patients into four subgroups based on a combination of PD-L1 expression and the type of immune modulation, and the Kaplan-Meier curves for overall survival are shown in Fig. 1f. The PD-L1-low-exhausted group had a worse prognosis compared with the other groups. We next investigated the clinical characteristics across subgroups and PD-L1 expression levels. The subgroups showed similar clinical characteristics. In addition, previously published bladder cancer molecular subtypes were analysed (Table 1). The basal and luminal subtypes were divided by TCGA into luminal-papillary, luminal-infiltrated, luminal, basal-squamous and neuronal subtypes . We detected the highest proportion of the luminal papillary subgroup within the PD-L1-low-activated group (chi-square, 80.9%, P < 0.0001), whereas the basal squamous subgroup harboured the highest proportion of the PD-L1-high-activated and PD-L1-high-exhausted groups (chi-square, 49.3 and 67.6%, P < 0.0001). It has been suggested that the luminal-infiltrated subtype exhibited increased PD-L1 and PD-1 expression . Interestingly, we found that approximately 53.1% of luminal-infiltrated subtype samples belong to the PD-L1-low-exhausted group and 37.5% of luminal-infiltrated subtype samples belong to the PD-L1-high-exhausted group, suggesting that a portion of the luminal-infiltrated subtype with decreased PD-L1 expression may exhibit poor anti-PD-1 treatment effectiveness.
The immune states, fraction of infiltrating immune cells and chemosensitivity of the four subgroups
Five immune expression signatures were selected to analyse the characteristics of the immune subgroups. The four immune subgroups were characterized by different immune signatures and DNA damage (Fig. 2a, b) and showed distinct immune characteristics (Fig. 2d).
The measures of DNA damage included copy number variation (CNV) burden (including number of segments and fraction of genome alterations), aneuploidy, homologous recombination deficiency (HRD), and intratumour heterogeneity (ITH). Antigen-specific BCR repertoires are critical for the recognition of tumour cells and reflect a robust antitumour response. Gene expression of IMs is critical for immunotherapy of tumours with many ongoing IM agonists and antagonist clinical trials. The leukocyte fraction (LF) was defined as the total leukocyte infiltration .
The PD-L1-high-actived group had the highest proliferation rate (Fig. 2d), the lowest CTA scores and high IFN-γ signature (vs. PD-L1-high-exhausted, P = 0.235; vs. PD-L1-low-actived, P < 0.001; vs. PD-L1-low-exhausted, P < 0.001), TGF-β signature (vs. PD-L1-high-exhausted, P < 0.001; vs. PD-L1-low-actived, P < 0.001; vs. PD-L1-low-exhausted, P < 0.001), lymphocyte infiltration (vs. PD-L1-high-exhausted, P < 0.001; vs. PD-L1-low-actived, P < 0.001; vs. PD-L1-low-exhausted, P = 0.583), and macrophage regulation (vs. PD-L1-high-exhausted, P < 0.001; vs. PD-L1-low-actived, P < 0.001; vs. PD-L1-low-exhausted, P = 0.538) and wound healing signatures (vs. PD-L1-high-exhausted, P = 0.750; vs. PD-L1-low-actived, P < 0.001; vs. PD-L1-low-exhausted, P < 0.001). This group displayed the highest number of SNV neoantigens and the highest nonsilent mutation rate.
The PD-L1-high-exhausted group exhibited the highest IFN-γ signature (vs. PD-L1-low-actived, P < 0.001; vs. PD-L1-low-exhausted, P < 0.001), TGF-β signature (vs. PD-L1-low-actived, P < 0.001; vs. PD-L1-low-exhausted, P < 0.001), lymphocyte infiltration (vs. PD-L1-low-actived, P < 0.001; vs. PD-L1-low-exhausted, P < 0.001), macrophage regulation (vs. PD-L1-low-actived, P < 0.001; vs. PD-L1-low-exhausted, P < 0.001) and wound healing signatures (vs. PD-L1-low-actived, P < 0.001; vs. PD-L1-low-exhausted, P < 0.001). This group also exhibited the highest leukocyte fraction and a high number of SNV neoantigens.
The PD-L1-low-actived group was defined by the lowest leukocyte fraction and IFN-γ signature (vs. PD-L1-low-exhausted, P = 0.040), TGF-β signature (vs. PD-L1-low-exhausted, P < 0.001), lymphocyte infiltration (vs. PD-L1-low-exhausted, P < 0.001), macrophage regulation (vs. PD-L1-low-exhausted, P < 0.001) and wound healing signatures (vs. PD-L1-low-exhausted, P < 0.230). In addition, this group exhibited the highest B cell receptor (BCR) expression, lowest aneuploidy score and lowest number of homologous recombination defects.
The PD-L1-low-exhausted group displayed the highest macrophage M2 infiltration, CTA scores, a low-to-moderate proliferation rate, IFN-γ signature, TGF-β signature, lymphocyte infiltration, macrophage regulation and wound healing signatures.
Gene expression of IMs (Fig. 2c) varied across immune subgroups. The top-ranked genes with differences between subgroups included CXCL9, CXCL10, CCL5, TNFRSF18, ITGB2, ICAM1 and all the HLA genes (consistent with their known interferon inducibility), which were most highly expressed in the PD-L1-high-exhausted group (P < 0.001).
We next estimated the mean fractions of immune cells in the four subgroups (Fig. 2e). The PD-L1-high-actived and PD-L1-high-exhausted groups exhibited the highest infiltration of M1 macrophages and CD8 T cells (P < 0.001). The PD-L1-high-actived and PD-L1-low-actived groups showed the highest infiltration of activated DCs (P < 0.001), whereas the PD-L1-low-actived group exhibited the highest infiltration of Tregs (P < 0.001). The PD-L1-low-exhausted group showed the highest infiltration of naïve B cells and resting mast cells, whereas the lowest number of activated mast cells and follicular helper T cells was also observed in this group (P < 0.001).
The TMB was significantly higher in the PD-L1-high expression groups compared with the PD-L1-low expression groups [15.06 (PD-L1-high-actived) and 22.32 (PD-L1-high-exhausted) vs. 2.87 (PD-L1-low-actived) and 5.35 (PD-L1-low-exhausted), P < 0.001], which confirmed that the effectiveness of anti-PD-1 treatment in the PD-L1-high expression groups [7.64 (PD-L1-high-actived) and 7.69 (PD-L1-high-exhausted) vs. 6.55 (PD-L1-low-actived) and 6.64 (PD-L1-low-exhausted), P < 0.001] could be related to TMB (Fig. 2f). In addition, anti-CTLA4 treatment showed no significant differences among the four subgroups.
Recently, stromal immunotypes were demonstrated to be a predictive model to identify patients who would benefit from chemotherapy in bladder cancer . We next estimated the chemosensitivity in the four subgroups. The log-transformed IC50 values are shown in Fig. 2g. Most drugs showed the highest sensitivities in the PD-L1-high-exhausted group, including nilotinib, belinostat, paclitaxel, vinblastine, doxorubicin, trametinib, bortezomib, idelalisib, cabozantinib, and bleomycin. Dabrafenib and lenalidomide showed the highest sensitivities in the PD-L1-high-activated group. Only axitinib showed sensitivity in the PD-L1-low-activated and PD-L1-low-exhausted groups. In addition, in the PD-L1-high groups, patients who underwent chemotherapy exhibited better overall survival rates compared with patients who did not undergo chemotherapy (P = 0.0005), whereas there was no significant difference in the PD-L1-low groups (P = 0.6066) (Fig. 2h). Finally, to elucidate the crosstalk between tumour cells and immune cells, we estimated the expression of several factors. FGFR3, which supported cancer-associated fibroblast (CAF) survival and activation, exhibited significantly higher expression in the PD-L1-low-actived and PD-L1-low-exhausted groups compared with the other groups. CXCL12, which led to PI3K and ERK1/2 activation in tumour cells, and CCL2, which was released by stromal cells and tumour cells to recruit myeloid cells, exhibited significantly increased expression in the exhausted groups than in the active groups (Fig. 2i) .
GSEA to predict key pathways
We performed GSEA to explore the critical pathways in the PD-L1-low-exhausted group (Fig. 3a). The upregulated mRNAs were enriched in the myogenesis (normalized enrichment score [NES] = 3.196, nominal P-value [NOM P-val] < 0.0001, enrichment score [ES] = 0.675), epithelial-mesenchymal transition (NES = 2.763, NOM P-val < 0.0001, ES = 0.594) and adipogenesis (NES = 1.731, NOM P-val < 0.0001, ES = 0.375) pathways. The downregulated mRNAs were enriched in the interferon alpha response (NES = − 3.302, NOM P-val < 0.0001, ES = − 0.784), interferon gamma response (NES = − 3.179, NOM P-val < 0.0001, ES = − 0.678) and E2F target (NES = − 2.977, NOM P-val < 0.0001, ES = − 0.630) pathways (Fig. 3b).
Exploring therapeutic targets through CNV analysis
We evaluated the CNV patterns that occurred in the PD-L1-low-exhausted group but not in the other groups. The total number of different CNV events in the PD-L1-low-exhausted group was 381 (P < 0.05), and the CNV events mainly included losses of chromosomes 3, 10 and 14 and gains of chromosomes 4 and 9 (P < 0.01) (Fig. 4a). Then, we assessed the relationship between CNVs and mRNA expression levels. The results revealed that 141 CNV-driven changes in mRNA expression of the 381 CNVs were associated with the PD-L1-low-exhausted group. To identify the expression characteristics of the PD-L1-low-exhausted group, we further investigated DEGs (Fig. 4b). A total of 6.25% (3/48) of the DEGs associated with the 141 CNVs were related to the expression level of PD-L1, including LCNL1, FBP1 and RASL11B (Fig. 4c, d). In addition, Kaplan-Meier survival curves showed that RASL11B was significant for predicting patient overall survival and disease-free survival (P = 0.026 and 0.042, respectively) (Fig. 4e).
Through analysis of the GEO validation data, we found that RASL11B (PD-L1-low-exhausted group: 10.35 ± 1.12; others: 9.99 ± 1.47; P = 0.017) showed higher expression in the PD-L1-low-exhausted group (n = 119) compared with the other groups (n = 182), which was comparable to the results in TCGA dataset.
It is well established that immune status can both support and obstruct tumour progression and therapeutic efficacy . Recently, Richard A. Moffitt et al. published a genetic “active stromo” signature and suggested that patients with an activated stroma had worse survival than patients with a normal stroma. This allowed us to explore stroma-specific subgroups with prognostic and biologic relevance . Tumour-associated macrophages (TAMs) are key regulators of the therapeutic response that produce several suppressive cytokines, thus contributing to T cell suppression. As a key factor of another path to T cell suppression, PD-L1 had different levels between the activated and exhausted groups. To explore the interaction of the two elements, we divided patients into four subgroups based on a combination of PD-L1 expression and the type of immune modulation. Interestingly, the PD-L1-low-exhausted group showed a worse prognosis than the other groups.
To analyse the characteristics of our immune subgroups, five immune expression signatures were selected, and the four subgroups were characterized by distinct immune signatures. The PD-L1-low-exhausted group had the worst prognosis and displayed composite signatures related to a low to moderate proliferation rate, IFN-γ signature, TGF-β signature lymphocyte infiltration, macrophage regulation and wound healing and had the highest M2 macrophage content and CTA scores; these findings are consistent with an immunosuppressed TME, for which a poor outcome would be expected. In contrast, the two subgroups displaying an active immune response (a high IFN-γ signature), the PD-L1-high-actived and PD-L1-high-exhausted group, had the most favourable prognosis, which is consistent with studies suggesting that a type I immune response is needed for antitumour activity . The PD-L1-high-actived group demonstrated the most SNV neoantigens and highest nonsilent mutation rate, in agreement with recent research suggesting that patients with high neoantigen number and high mutation rate exhibit the longest survival in bladder cancer. DNA damaging therapeutics such as radiotherapy and chemotherapy may help liberate neoantigens to induce T cell responses and improve the effects of immunotherapy in the group with the second highest number of SNV neoantigens, the PD-L1-high-exhausted group [26, 29, 30]. The PD-L1-low-actived group was defined by decreased leukocyte chemotaxis, similar to an immunologically quiet subtype , leading to fewer immune cells and better outcomes .
To further elucidate the mechanism, we estimated the fractions of immune cells in the four subgroups. M1 macrophages and CD8+ T cells highly infiltrated the tumours in the PD-L1-high-actived and PD-L1-high-exhausted groups. It is well established that CD8+ T cells as well as macrophages are required for synergistic curative activity in response to immune checkpoint antibodies (antibodies against CTLA-4, PD-1, and PD-L1) . High PD-L1 levels act as a target to activate high numbers of infiltrating M1 macrophages and CD8+ T cells. In the PD-L1-high subgroups, patients who underwent chemotherapy exhibited better overall survival rates than patients who did not undergo chemotherapy, which may have inhibited high chemosensitivity in these two subgroups.
It should be noted that most of the chemotherapeutic drugs tested in this study had the highest sensitivities in the PD-L1-high-exhausted group, but chemotherapy suppresses the antitumour response of CD8+ T cells [32, 33] and induces IL-10-producing M2 macrophages, which may lead to transformation towards a PD-L1-low-exhausted status .
The highest infiltration of Tregs and the lowest infiltration of CD8+ T cells were found in the PD-L1-low-actived group. It was reported that STAT3 inhibition in combination with radiation improved tumour growth delay, decreased Tregs and enhanced effector T cells and M1 macrophages, which may improve curative effects . In addition, the PD-L1-low-actived group exhibited the lowest infiltration of M1 macrophages, which may cause low chemosensitivity. TAMs protect CSCs from chemotherapy and secrete cathepsins that may blunt the chemosensitivity of cancer cells [36, 37]. An agonist CD40 antibody was found to modify macrophages to a tumouricidal phenotype in pancreatic cancer . Moreover, an agonistic CD40 monoclonal antibody (CP-870,893) in combination with gemcitabine showed partial clinical effects in advanced pancreatic cancer patients . Therefore, a CD40 antibody may improve the chemosensitivity of patients in the PD-L1-low-actived group. The results revealed that axitinib showed sensitivity in the PD-L1-low-exhausted group. Combination therapies of immune checkpoint inhibitors with targeted agents have demonstrated both safety and efficacy and have received FDA approval for the first-line treatment of advanced RCC . For their use in the treatment of bladder cancer, more clinical trials are needed. The PD-L1-low-exhausted group showed high infiltration of M2 macrophages and low infiltration of CD8+ T cells. Therefore, the key to transformation to the tumouricidal type is CD8+ T cell proliferation and activation. However, the poor antigen-presenting cell (APC) status of the PD-L1-low-exhausted group leads to inefficient presentation of antigen to CD8+ T cells and includes high infiltration of naïve B cells and low T follicular helper cell and DC activation. T follicular helper cells help B cells differentiate into plasma and memory B cells .
In conclusion, to reverse the poor APC status, improving T follicular helper cell and DC activation may be helpful. Treatment with FLT3 L/poly I:C, which induces the activation of CD103+ cDC1s at tumour sites and induce proliferation of naïve CD8 T cells and generated CTLs, has been shown to enhance antitumour responses [42,43,44].
We performed GSEA to explore the critical pathways in the PD-L1-low-exhausted group. Interestingly, the upregulated mRNAs were enriched in the epithelial mesenchymal transition (EMT) pathway. Recently, Li Wang et al. demonstrated that in metastatic urothelial cancer patients treated with a PD-1 antibody, nivolumab, higher expression of EMT-related genes was associated with lower response rates and shorter PFS and OS. Therefore, cotargeting PD-1 and stromal elements may reverse immune resistance.
The upregulated mRNAs in the other group were enriched in the E2F target pathway. E2F target pathway activation leads to an increase in the levels of proteins required for DNA synthesis, eventually allowing the cell to duplicate itself. For the treatment of other groups, it is well established that the effects of CDK4/6 inhibitors reduce the activity of the E2F target DNA methyltransferase 1 and suppress the proliferation of tumour cells and regulatory T cells .
The genome-wide CNV study revealed 381 gene-level areas of gains or deletions. The results revealed that 141 CNV-driven changes in mRNA expression of the 381 CNVs were associated with the PD-L1-low-exhausted group. DEG testing gave a set of 48 gene-level transcripts consistently associated with the PD-L1-low-exhausted group. From these results, the CNV analysis provided translational benefits in DEGs including LCNL1, FBP1 and RASL11B. RASL11B was significant for predicting patient overall survival and disease-free survival. RASL11B is a small GTPase belonging to a Ras subfamily of putative tumour suppressor genes. Upregulated RASL11B inhibits cell proliferation, invasion, and migration, induces G0/G1 cell cycle arrest and promotes cell apoptosis in clear cell renal cell carcinoma (CCRCC) . RASL11B has promise as a therapeutic target.
In summary, four stable immune subgroups were found in bladder cancer (Fig. 5). These subgroups were found to be associated with prognosis and genetic and immune modulatory alterations that shape immune environments types. Patients in the PD-L1-high-actived and PD-L1-high-exhausted groups were found to have a significantly higher overall survival and better response to chemotherapy and immunotherapy, suggesting that cancer therapeutics can be tailored based on the immune subgroups of the tumour microenvironment. Regarding the lower number of SNV neoantigens in the PD-L1-high-exhausted subgroup, radiotherapy and chemotherapy in combination with anti-PD-1/PD-L1 therapy may be useful in these cases. In the immunologically quiet stage, PD-L1-low-actived patients showed the highest infiltration of Tregs and the lowest infiltration of M1 macrophages and CD8+ T cells and thus may benefit from STAT3 inhibition in combination with radiation or an agonistic CD40 antibody. In addition, pembrolizumab or avelumab with axitinib has demonstrated both safety and efficacy in the first-line treatment of advanced RCC and is recommended for PD-L1-low-actived-type patients, who need more clinical trials. PD-L1-low-exhausted patients showed the worst APC immunosuppression status and had the worst prognosis, and therapeutic approaches targeting RASL11B to improve T cell follicular helper and DC activation (FLT3 L/poly I:C) may be helpful. With an increasing understanding that the immune environment of tumours plays an important role in the response to therapy, defining the immune subtype plays a critical role in predicting disease outcome and strategy effectiveness.
The role of the microenvironment in tumours along with various cell types along with crosstalk was assessed. To elucidate the crosstalk between tumour cells and stroma, we estimated the expression of some factors. Soluble factors secreted by tumour cells and stroma promote angiogenesis and support tumour cell survival. The interactions between tumour cells and stroma mediated through CXCL12 lead to activation of the PI3K and ERK1/2 pathways, which are associated with tumour viability and the active EMT pathway, promoting cell survival and metastasis. CCL2 released by CAFs and tumour cells recruits myeloid-derived suppressor cells (MDSCs) and macrophages. TAMs, which promote tumour cell viability, secrete FGF, which binds with FGFR3 and supports CAF survival and activation. Furthermore, TAMs and MDSCs suppress T-cell function in the TME .
Several limitations should be addressed in this study. First, we adopted the CD274 mRNA expression level as a surrogate for PD-L1 expression. Second, the retrospective nature of TCGA causes bias from variations in follow-up information. Furthermore, our findings need to be validated at the protein level.
In this study, we identified common immune subgroups based on a combination of PD-L1 expression and the type of immune modulation and performed a series of analyses. This suggests the need for prospective studies in immune subgroup patients, especially to better address drug design and understand the outcome of gene-specific therapies.
The PD-L1-low-exhausted group exhibited a worse prognosis than the other groups and showed tumouricidal features (high infiltration of M2 macrophages and low infiltration of CD8+ T cells) and poor APC immune status and activated CAF status. Targeting the immune status of the PD-L1-low-exhausted group would have a tremendous impact on individual treatment.
Copy number variations (CNVs) and differentially expressed genes upregulated in the PD-L1-low-exhausted group, such as RASL11B, played a role in predicting overall survival in TCGA data and was validated in the PD-L1-low-exhausted group in the GEO database. This important result also highlights RASL11B as a therapeutic target, but further investigations are needed.
We believe that reporting data from a large patient cohort might be extremely valuable for many aspects, such as providing patients with a better understanding of tumour immune mechanisms and therapy.
Availability of data and materials
The datasets generated and/or analysed during the current study are available in the Genomic Data Commons Data Portal (https://portal.gdc.cancer.gov/). The public access to the database is open. The authors did not have special access privileges.
Kamat AM, Hahn NM, Efstathiou JA, Lerner SP, Malmstrom PU, Choi W, et al. Bladder cancer. Lancet. 2016;388(10061):2796–810.
Antoni S, Ferlay J, Soerjomataram I, Znaor A, Jemal A, Bray F. Bladder Cancer incidence and mortality: a global overview and recent trends. Eur Urol. 2017;71(1):96–108.
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020. CA Cancer J Clin. 2020;70(1):7–30.
Sanli O, Dobruch J, Knowles MA, Burger M, Alemozaffar M, Nielsen ME, et al. Bladder cancer. Nat Rev Dis Primers. 2017;3:17022.
Bellmunt J, de Wit R, Vaughn DJ, Fradet Y, Lee JL, Fong L, et al. Pembrolizumab as second-line therapy for advanced Urothelial carcinoma. N Engl J Med. 2017;376(11):1015–26.
Sharma P, Retz M, Siefker-Radtke A, Baron A, Necchi A, Bedke J, et al. Nivolumab in metastatic urothelial carcinoma after platinum therapy (CheckMate 275): a multicentre, single-arm, phase 2 trial. Lancet Oncol. 2017;18(3):312–22.
Apolo AB, Infante JR, Balmanoukian A, Patel MR, Wang D, Kelly K, et al. Avelumab, an anti-programmed death-ligand 1 antibody, in patients with refractory metastatic Urothelial carcinoma: results from a multicenter, Phase Ib Study. J Clin Oncol. 2017;35(19):2117–24.
Massard C, Gordon MS, Sharma S, Rafii S, Wainberg ZA, Luke J, et al. 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. 2016;34(26):3119–25.
Balar AV, Galsky MD, Rosenberg JE, Powles T, Petrylak DP, Bellmunt J, et al. Atezolizumab as first-line treatment in cisplatin-ineligible patients with locally advanced and metastatic urothelial carcinoma: a single-arm, multicentre, phase 2 trial. Lancet. 2017;389(10064):67–76.
Rosenberg JE, Hoffman-Censits J, Powles T, van der Heijden MS, Balar AV, Necchi A, et al. Atezolizumab in patients with locally advanced and metastatic urothelial carcinoma who have progressed following treatment with platinum-based chemotherapy: a single-arm, multicentre, phase 2 trial. Lancet. 2016;387(10031):1909–20.
Wang L, Saci A, Szabo PM, Chasalow SD, Castillo-Martin M, Domingo-Domenech J, et al. EMT- and stroma-related gene expression and resistance to PD-1 blockade in urothelial cancer. Nat Commun. 2018;9(1):3503.
Wu T, Dai Y. Tumor microenvironment and therapeutic response. Cancer Lett. 2017;387:61–8.
Zeng D, Li M, Zhou R, Zhang J, Sun H, Shi M, et al. Tumor microenvironment characterization in gastric Cancer identifies prognostic and Immunotherapeutically relevant gene signatures. Cancer Immunol Res. 2019;7(5):737–50.
Powles T, Kockx M, Rodriguez-Vida A, Duran I, Crabb SJ, Van Der Heijden MS, et al. Clinical efficacy and biomarker analysis of neoadjuvant atezolizumab in operable urothelial carcinoma in the ABACUS trial. Nat Med. 2019;25(11):1706–14.
Kim J, Kwiatkowski D, McConkey DJ, Meeks JJ, Freeman SS, Bellmunt J, et al. The Cancer genome atlas expression subtypes stratify response to checkpoint inhibition in advanced Urothelial Cancer and identify a subset of patients with high survival probability. Eur Urol. 2019;75(6):961–4.
Fong MHY, Feng M, McConkey DJ, Choi W. Update on bladder cancer molecular subtypes. Transl Androl Urol. 2020;9(6):2881–9.
Rooney MS, Shukla SA, Wu CJ, Getz G, Hacohen N. Molecular and genetic properties of tumors associated with local immune cytolytic activity. Cell. 2015;160(1–2):48–61.
Amankulor NM, Kim Y, Arora S, Kargl J, Szulzewsky F, Hanke M, et al. Mutant IDH1 regulates the tumor-associated immune system in gliomas. Genes Dev. 2017;31(8):774–86.
Moffitt RA, Marayati R, Flate EL, Volmar KE, Loeza SG, Hoadley KA, et al. Virtual microdissection identifies distinct tumor- and stroma-specific subtypes of pancreatic ductal adenocarcinoma. Nat Genet. 2015;47(10):1168–78.
Charoentong P, Finotello F, Angelova M, Mayer C, Efremova M, Rieder D, et al. Pan-cancer Immunogenomic analyses reveal genotype-Immunophenotype relationships and predictors of response to checkpoint blockade. Cell Rep. 2017;18(1):248–62.
Wang Y, Wang Z, Xu J, Li J, Li S, Zhang M, et al. Systematic identification of non-coding pharmacogenomic landscape in cancer. Nat Commun. 2018;9(1):3192.
Thorsson V, Gibbs DL, Brown SD, Wolf D, Bortone DS, Ou Yang TH, et al. The immune landscape of Cancer. Immunity. 2018;48(4):812–30 e814.
Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98–W102.
Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144(5):646–74.
Sia D, Jiao Y, Martinez-Quetglas I, Kuchuk O, Villacorta-Martin C, Castro de Moura M, et al. Identification of an immune-specific class of hepatocellular carcinoma, Based on Molecular Features. Gastroenterology. 2017;153(3):812–26.
Robertson AG, Kim J, Al-Ahmadie H, Bellmunt J, Guo G, Cherniack AD, et al. Comprehensive molecular characterization of muscle-invasive bladder Cancer. Cell. 2017;171(3):540–56 e525.
Fu H, Zhu Y, Wang Y, Liu Z, Zhang J, Xie H, et al. Identification and validation of stromal Immunotype predict survival and benefit from adjuvant chemotherapy in patients with muscle-invasive bladder Cancer. Clin Cancer Res. 2018;24(13):3069–78.
Galon J, Angell HK, Bedognetti D, Marincola FM. The continuum of cancer immunosurveillance: prognostic, predictive, and mechanistic signatures. Immunity. 2013;39(1):11–26.
Brown JS, Sundar R, Lopez J. Combining DNA damaging therapeutics with immunotherapy: more haste, less speed. Br J Cancer. 2018;118(3):312–24.
Chen YP, Zhang Y, Lv JW, Li YQ, Wang YQ, He QM, et al. Genomic analysis of tumor microenvironment immune types across 14 solid Cancer types: immunotherapeutic implications. Theranostics. 2017;7(14):3585–94.
Saha D, Martuza RL, Rabkin SD. Macrophage polarization contributes to Glioblastoma eradication by combination Immunovirotherapy and immune checkpoint blockade. Cancer Cell. 2017;32(2):253–67 e255.
DeNardo DG, Brennan DJ, Rexhepaj E, Ruffell B, Shiao SL, Madden SF, et al. Leukocyte complexity predicts breast cancer survival and functionally regulates response to chemotherapy. Cancer Discov. 2011;1(1):54–67.
Nakasone ES, Askautrud HA, Kees T, Park JH, Plaks V, Ewald AJ, et al. Imaging tumor-stroma interactions during chemotherapy reveals contributions of the microenvironment to resistance. Cancer Cell. 2012;21(4):488–503.
Dijkgraaf EM, Heusinkveld M, Tummers B, Vogelpoel LT, Goedemans R, Jha V, et al. Chemotherapy alters monocyte differentiation to favor generation of cancer-supporting M2 macrophages in the tumor microenvironment. Cancer Res. 2013;73(8):2480–92.
Oweida AJ, Darragh L, Phan A, Binder D, Bhatia S, Mueller A, et al. STAT3 modulation of regulatory T cells in response to radiation therapy in head and neck Cancer. J Natl Cancer Inst. 2019;111(12):1339–49.
Jinushi M, Chiba S, Yoshiyama H, Masutomi K, Kinoshita I, Dosaka-Akita H, et al. Tumor-associated macrophages regulate tumorigenicity and anticancer drug responses of cancer stem/initiating cells. Proc Natl Acad Sci U S A. 2011;108(30):12425–30.
Shree T, Olson OC, Elie BT, Kester JC, Garfall AL, Simpson K, et al. Macrophages and cathepsin proteases blunt chemotherapeutic response in breast cancer. Genes Dev. 2011;25(23):2465–79.
Beatty GL, Chiorean EG, Fishman MP, Saboury B, Teitelbaum UR, Sun W, et al. CD40 agonists alter tumor stroma and show efficacy against pancreatic carcinoma in mice and humans. Science. 2011;331(6024):1612–6.
Beatty GL, Torigian DA, Chiorean EG, Saboury B, Brothers A, Alavi A, et al. A phase I study of an agonist CD40 monoclonal antibody (CP-870,893) in combination with gemcitabine in patients with advanced pancreatic ductal adenocarcinoma. Clin Cancer Res. 2013;19(22):6286–95.
Zarrabi K, Paroya A, Wu S. Emerging therapeutic agents for genitourinary cancers. J Hematol Oncol. 2019;12(1):89.
Shannack MM, Linterman MA. T Follicular Helper Cells. In: Bradshaw RA, Stahl PD, editors. Encyclopedia of Cell Biology. Waltham: Academic Press; 2016. p. 541–9.
Fu C, Jiang A. Dendritic cells and CD8 T cell immunity in tumor microenvironment. Front Immunol. 2018;9:3059.
Salmon H, Idoyaga J, Rahman A, Leboeuf M, Remark R, Jordan S, et al. Expansion and activation of CD103(+) dendritic cell progenitors at the tumor site enhances tumor responses to therapeutic PD-L1 and BRAF inhibition. Immunity. 2016;44(4):924–38.
Broz ML, Binnewies M, Boldajipour B, Nelson AE, Pollack JL, Erle DJ, et al. Dissecting the tumor myeloid compartment reveals rare activating antigen-presenting cells critical for T cell immunity. Cancer Cell. 2014;26(5):638–52.
Goel S, DeCristo MJ, Watt AC, BrinJones H, Sceneay J, Li BB, et al. CDK4/6 inhibition triggers anti-tumour immunity. Nature. 2017;548(7668):471–5.
He H, Dai J, Zhuo R, Zhao J, Wang H, Sun F, et al. Study on the mechanism behind lncRNA MEG3 affecting clear cell renal cell carcinoma by regulating miR-7/RASL11B signaling. J Cell Physiol. 2018;233(12):9503–15.
This work is supported by Shenyang Bureau of Science and Technology Project No. 17–230–9-38. The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
The biopsy samples were collected from the patients after obtaining written informed consent. The study was approved by the First Hospital of China Medical University Institutional Ethics Committee. All data was anonymized, and all identifiable information and biological samples were stored according to the local guidelines.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Lyu, X., Wang, P., Qiao, Q. et al. Genomic stratification based on microenvironment immune types and PD-L1 for tailoring therapeutic strategies in bladder cancer. BMC Cancer 21, 646 (2021). https://doi.org/10.1186/s12885-021-08350-1
- Immune subtypes
- Bladder cancer
- Therapeutic strategies