Muscle-invasive urothelial bladder cancer is a common malignancy with poor outcomes for which immune checkpoint blockade is now showing promise. Despite clinical activity of PD-1/PD-L1–targeted therapy in this disease, most patients do not benefit and resistance mechanisms remain unknown. The non–T-cell-inflamed tumor microenvironment correlates with poor prognosis and resistance to immunotherapies. In this study, we determined tumor-oncogenic pathways correlating with T-cell exclusion. We first establish in this report that T-cell–inflamed bladder tumors can be identified by immune gene expression profiling with concordance with CD8+ T-cell infiltration. Upregulation of genes encoding immune checkpoint proteins PD-L1, IDO, FOXP3, TIM3, and LAG3 was associated with T-cell–inflamed tumors, suggesting potential for sensitivity to checkpoint blockade. β-Catenin, PPAR-γ, and FGFR3 pathways were activated in non–T-cell-inflamed tumors. No difference was seen in overall somatic mutational density between groups. The three pathways identified represent targetable potential pathways of tumor-intrinsic immunotherapy resistance. Cancer Immunol Res; 4(7); 563–8. ©2016 AACR.

Urothelial bladder cancer is the most common malignancy of the urinary system and estimated to cause 165,084 annual deaths worldwide (1). Despite chemotherapy, the outcomes for metastatic urothelial cancer remain poor with an overall 5-year survival rate of 15% (2). Given the need for novel approaches in this disease, the clinical exploration of immunotherapeutic checkpoint blockade is occurring rapidly (3–5). Inhibition of PD-1/PD-L1 interactions in bladder cancer has shown an objective response rate (ORR) of 52% in patients with tumor-infiltrating immune cells showing PD-L1 expression (6). Despite these promising data, the ORR rate in PD-L1–negative patients was just 11%. Thus, the majority of patients with urothelial cancer do not benefit from immune checkpoint blockade, and mechanisms of resistance remain incompletely understood.

The presence of an existing antitumor T-cell response in all cancers is important for the activity of immunotherapies, including checkpoint blockade (7, 8). Tumors with a T-cell–inflamed microenvironment are characterized by infiltration of CD8+ T cells, chemokines, and an interferon signature, which are associated with improved survival and response to immunotherapies (9–12). In urothelial bladder cancer, increased T-cell infiltration has been correlated with longer patient survival (13). Efforts have now been focused on understanding mechanisms driving T-cell exclusion. Recently, β-catenin activation has been identified in melanoma as the first tumor-oncogenic pathway leading to a non–T-cell-inflamed phenotype (14). Stabilization of β-catenin leads to ATF3-dependent transcriptional repression of CCL4, which ablates recruitment of Batf3-lineage dendritic cells into the tumor microenvironment. This, in turn, leads to failed T-cell priming, defective T-cell recruitment, and loss of therapeutic response to anti–PD-L1 and anti–CTLA-4 therapies in a murine model. In other cancer types, it is not yet known whether the same or additional tumor-intrinsic signaling pathways affect the development of a T-cell–inflamed phenotype. We therefore pursued the identification of pathways and mutations associated with non–T-cell-inflamed urothelial bladder tumors to discover new candidate therapeutic targets for improving immunotherapy efficacy.

Identification of T-cell–inflamed and non–T-cell-inflamed subtypes in TCGA

The Cancer Genome Atlas (TCGA) Bladder Urothelial Carcinoma (BLCA) RNA-seq gene expression data were downloaded from Broad Genome Data Analysis Center (http://gdac.broadinstitute.org, release date July 15, 2014) as previously described (14). The upper quartile-normalized and log2-transformed RNA-seq by expectation maximization (RSEM) values for 267 samples was summarized at gene level (15). Genes expressed in <80% of the samples were filtered. Exome somatic mutation data were downloaded from the TCGA portal (November 16, 2014) for 238 tumor-normal sample pairs with variants from “BI Automated Mutation Calling.”

Unsupervised hierarchical clustering using R function hclust v3.1.0 with K equal to 10 was conducted on 16,197 genes. A 725-gene cluster contained 12 from our previously described T-cell signature (CD8A, CCL2, CCL3, CCL4, CXCL9, CXCL10, ICOS, GZMK, HLA-DMA, HLA-DMB, HLADOA, and HLA-DOB). Consensus clustering of the primary tumor samples was performed using BioConductor package ConsensusClusterPlus v1.18.0 (16) based on expression of the 725 genes, with 2,000 resampling cycles, using hierarchical clustering and Euclidean distance. Three distinct tumor subtypes were identified, with low, moderate, and high expression level of the T-cell genes (Fig. 1A). Principal component analysis (PCA) mapping showed separation of subtypes.

Figure 1.

A, T-cell–inflamed and non–T-cell-inflamed bladder tumors can be distinguished by immune gene expression profiling. Eighty-eight (33%) of the tumors show minimal expression of T-cell–related immune genes (non–T-cell-inflamed), whereas 95 (36%) show overexpression (T-cell–inflamed). The remainders show a mixed expression pattern of T-cell–related genes. B, representative examples of CD8 IHC showing an absence of intratumoral T cells in a non–T-cell-inflamed tumor (left) and marked T-cell infiltration in a T-cell–inflamed tumor (right). C, immune subtypes identified by immune gene expression profiling show significant association with the presence of T cells by IHC, defined as having >1 CD8+ T cells per 0.1 mm2 field (P = 0.01, Fisher exact test, n = 19).

Figure 1.

A, T-cell–inflamed and non–T-cell-inflamed bladder tumors can be distinguished by immune gene expression profiling. Eighty-eight (33%) of the tumors show minimal expression of T-cell–related immune genes (non–T-cell-inflamed), whereas 95 (36%) show overexpression (T-cell–inflamed). The remainders show a mixed expression pattern of T-cell–related genes. B, representative examples of CD8 IHC showing an absence of intratumoral T cells in a non–T-cell-inflamed tumor (left) and marked T-cell infiltration in a T-cell–inflamed tumor (right). C, immune subtypes identified by immune gene expression profiling show significant association with the presence of T cells by IHC, defined as having >1 CD8+ T cells per 0.1 mm2 field (P = 0.01, Fisher exact test, n = 19).

Close modal

Characterization of activated pathways and differentially mutated genes

Genes differentially expressed between tumor groups were detected from the global set of 16,197 genes by using ANOVA and filtered by FDR q value < 0.01 and fold change ≥ 2.0 (or ≤ −2.0). Canonical pathways significantly enriched in the genes of interest were identified by Ingenuity Pathways Analysis (Ingenuity Systems; http://www.ingenuity.com) based on experimental evidence from the Ingenuity Knowledge Base (March 23, 2014). A total of 78,762 somatic variants were converted to variant cell format (VCF) and annotated with exonic function, population frequency, dbSNP, COSMIC database, and deleterious function prediction. Synonymous variants, those with no annotation available, and those with a population frequency >1%, were filtered. The resulting 55,996 variants were kept for further analysis. The variants were then summarized at the gene and patient level for comparison of mutation profiles between groups. Single-nucleotide variants (SNV) located in selected genes were analyzed using the Sorting Tolerant From Intolerant (SIFT) prediction algorithm (http://sift.jcvi.org/; ref. 17).

Immunohistochemistry

The immunohistochemical staining (IHC) on human bladder tumors was performed by the Human Tissue Resource Center (HTRC) at the University of Chicago using tissue from cystectomy specimens from muscle-invasive bladder cancer, corresponding to specimens previously submitted for inclusion in the TCGA. Staining was performed using a CD8-specific monoclonal antibody (CD8 clone C8/144B; NeoMarkers) and a β-catenin monoclonal antibody (clone CAT-5H1; Life Technologies) in combination with a secondary goat anti-mouse immunoglobulin G conjugated to an alkaline phosphatase (Biocare Medical). Slides were scanned using a CRi Panoramic Scan Whole Slide Scanner and viewed with Pannoramic Viewer 1.15.4 (3DHISTECH). CD8+ cells from three 0.1 mm2 fields were counted manually and averaged. Tumors were considered positive for intratumoral T cells if >1 CD8+ cell was present. For β-catenin, tumors with any nuclear staining were considered positive. A urologic pathologist (G.P. Paner) scored samples for intensity and distribution in a blinded fashion.

Classifying subtypes of urothelial bladder cancer by immune gene profiling

We first aimed to determine whether T-cell–inflamed and non–T-cell-inflamed bladder tumors could be identified using immune gene expression profiling. Genes from a previously derived T-cell signature were used to classify 267 samples of urothelial bladder cancer analyzed from the TCGA (7, 10). First, unsupervised hierarchical clustering was conducted on 16,197 genes. A cluster of 725 genes containing 12 genes from the T-cell signature was used to perform consensus clustering of tumor samples, which identified groups with high, intermediate, and low expression of the immune genes (Fig. 1A). Those groups were compared to identify pathways linked with the non–T-cell-inflamed tumor microenvironment. The non–T-cell-inflamed group comprised 88 samples (33%), whereas the T-cell–inflamed group represented 95 samples (36%). The remainder of samples showed intermediate expression of the T-cell signature genes. Thus, like metastatic melanoma, bladder cancers could be segregated based on the degree of expression of genes indicative of a T-cell–inflamed tumor microenvironment (14).

T-cell–inflamed cancers overexpress inhibitory molecules and were infiltrated by T cells

In order to validate the correlation between the T-cell–inflamed tumor microenvironment gene signature and actual presence of intratumoral T cells, a subset of 19 TCGA samples derived from our institution was analyzed by IHC for CD8. We found that 37% of samples showed a notable absence of intratumoral CD8+ T cells, whereas 63% showed significant CD8+ T-cell infiltration (Fig. 1B). T-cell–infiltrated tumors contained a median of 14 CD8+ T cells per high power field (range, 2–47). Comparison between IHC and gene expression profiling revealed that the presence of intratumoral CD8+ T cells was strongly associated with the T-cell–inflamed gene expression signature (P = 0.01; Fig. 1C).

In other cancers, expression of immune-inhibitory molecules has been linked with a T-cell–inflamed phenotype (9). Therefore, we evaluated the expression of inhibitory molecules in bladder cancers and examined them for correlation with the T-cell signature. Expression of PD-L1 mRNA was positively correlated with that of CD8A (P < 0.0001; Fig. 2A), and the T-cell–inflamed tumors had higher expression of both genes as compared with the non–T-cell-inflamed tumors. A similar positive correlation was observed with other immune-inhibitory markers IDO, FOXP3, TIM3, and LAG3 (Fig. 2B). As in melanoma, T-cell–inflamed bladder cancers are characterized by the highest expression of genes that reflect negative feedback pathways that limit ongoing T-cell activation.

Figure 2.

A, expression of PD-L1 is positively correlated with expression of CD8A (P < 0.0001, Pearson correlation). T-cell–inflamed tumors (red) show higher expression of both genes compared with non–T-cell-inflamed tumors (blue). B, additional immune-inhibitory molecules, IDO1, FOXP3, TIM3, and LAG3 show a similarly significant correlation (all P < 0.0001, Pearson correlation).

Figure 2.

A, expression of PD-L1 is positively correlated with expression of CD8A (P < 0.0001, Pearson correlation). T-cell–inflamed tumors (red) show higher expression of both genes compared with non–T-cell-inflamed tumors (blue). B, additional immune-inhibitory molecules, IDO1, FOXP3, TIM3, and LAG3 show a similarly significant correlation (all P < 0.0001, Pearson correlation).

Close modal

Non–T-cell-inflamed tumors have activated β-catenin and PPAR-γ

Identification of molecular pathways active in tumor cells that might mediate T-cell exclusion from the tumor microenvironment has become a high priority as we move toward developing new therapeutic interventions to expand the fraction of patients responding to immunotherapies. Based on the recent characterization of the Wnt/β-catenin pathway in melanoma as being causal for preventing T-cell activation and trafficking into the tumor microenvironment (14), we pursued a similar analysis in bladder cancer. To identify pathways linked to the non–T-cell-inflamed phenotype, we compared gene expression patterns between T-cell–inflamed and non–T-cell-inflamed bladder cancers. We found 3,112 genes differentially expressed between the two groups (Q < 0.01; Supplementary Table S1), with 730 being preferentially expressed in the non–T-cell-inflamed tumors. We focused on the latter set of genes and evaluated for pathway enrichment using Ingenuity Pathways Analysis. The top activated upstream regulators identified were β-catenin and PPAR-γ (both P = 0.003). Thirty genes related to β-catenin signaling and 22 genes related to PPARG-γ were overexpressed in the non–T-cell-inflamed tumors (Fig. 3A and B).

Figure 3.

A and B, cellular localization of molecules upregulated in the non–T-cell-inflamed group and their relationship with (A) β-catenin and (B) PPARG. Arrows point from the regulator toward a molecule being regulated. Intensity of red color is proportional to level of upregulation. Dotted lines indicate indirect relationships. Lines highlighted in orange show predicted activation. C, all bladder tumors showed cytomembranous staining for β-catenin (left); however, a fraction of non–T-cell-inflamed tumors showed nuclear staining for β-catenin (right). D, nuclear β-catenin staining was only found in tumors with an absence of CD8+ T cells (P = 0.036, Fisher exact test, n = 19).

Figure 3.

A and B, cellular localization of molecules upregulated in the non–T-cell-inflamed group and their relationship with (A) β-catenin and (B) PPARG. Arrows point from the regulator toward a molecule being regulated. Intensity of red color is proportional to level of upregulation. Dotted lines indicate indirect relationships. Lines highlighted in orange show predicted activation. C, all bladder tumors showed cytomembranous staining for β-catenin (left); however, a fraction of non–T-cell-inflamed tumors showed nuclear staining for β-catenin (right). D, nuclear β-catenin staining was only found in tumors with an absence of CD8+ T cells (P = 0.036, Fisher exact test, n = 19).

Close modal

We further explored the two pathways by evaluating gene expression of PPARG and of WNT7B, the primary ligand associated with β-catenin signaling in melanoma (14). Both genes showed an inverse correlation with expression of CD8A (P < 0.0001), further suggesting a link between activation of those pathways and the absence of intratumoral T cells. To confirm this pattern at the protein level, we used IHC to examine expression of stabilized β-catenin versus extent of CD8+ T-cell infiltration. Indeed, an inverse relationship was observed between nuclear staining of β-catenin and the presence of intratumoral CD8+ T cells (P = 0.036; Fig. 3C and D). Nuclear staining was only found in samples without intratumoral CD8+ T cells, although it only accounted for 43% of the non–T-cell-inflamed tumors. These results are consistent with what was observed in melanoma, in which 48% of the non–T-cell-inflamed tumors showed β-catenin pathway activation (14).

To assess a potential underlying mechanism of T-cell exclusion in this tumor type, we evaluated expression of the chemokine CCL4 and the transcription factor BATF3, which are associated with recruitment of the critical population of dendritic cells necessary for T-cell activation in the melanoma context (14). CD8 transcripts were correlated with CCL4 and BATF3 expression within the tumor microenvironment (Supplementary Fig. S1A). We confirmed the expected inverse correlation with WNT7B expression, but also found that PPARG showed the same relative expression pattern (Supplementary Fig. S1B). Thus, the recruitment of BATF3-lineage dendritic cells necessary for cross presentation of tumor antigens might also be negatively influenced by PPARG. Lastly, of the 29 β-catenin target genes overexpressed in our dataset, only IHH and GATA2 overlapped with the 22 target genes of PPARG. Despite this lack of overlap, target gene expression for both pathways was strongly correlated (P < 0.0001), indicating a potential synergistic mechanism of immune exclusion.

Mutated FGFR3 and coactivators of PPAR-γ in non–T-cell-inflamed cancers

We next interrogated somatic mutations between T-cell–inflamed and non–T-cell-inflamed subtypes to identify mutations in pathways that might be exclusive to the non–T-cell-inflamed subset. Analysis of exome sequencing data revealed no difference in the overall number of nonsynonymous mutations per patient for each group (P = 0.80). Thus, a lower mutational density alone could not explain the absence of T cells in a major subset of tumors. We therefore reviewed differential mutations at the individual gene level. We found 891 genes mutated in two or more samples exclusive to noninflamed tumors (Supplementary Table S2). FGFR3 was the most common mutation exclusive to noninflamed tumors (P < 0.0001). Fourteen separate FGFR3 mutations were detected in 11 of the non–T-cell-inflamed samples, with no mutations detected in the T-cell–inflamed samples. Nine unique mutations were identified (Table 1), and two were recurrent (G370C in three samples and Y373C in four samples). The two most common FGFR3 mutations, Y373C and G370C, have been reported previously in bladder cancer and lead to activation of the receptor through increased ligand-independent dimerization and phosphorylation (18). In addition, recurrent in-frame activating FGFR3–TACC3 fusions were found exclusively in non–T-cell-inflamed tumors (three samples). We also analyzed FGFR3 gene expression and found an inverse correlation with expression of CD8A (P < 0.0001). Non–T-cell-inflamed tumors showed higher expression of FGFR3 as compared with the T-cell–inflamed tumors.

Table 1.

FGFR3 mutations found in non–T-cell-inflamed bladder cancers

TCGA sample IDMutationAmino acid substitution
TCGA-DK-A3IS-01 chr4:1803377:G:A E216K 
TCGA-DK-A3IS-01 chr4:1803395:G:A D222N 
TCGA-DK-A3IS-01 chr4:1803435:G:A G235D 
TCGA-FJ-A3ZF-01 chr4:1803738:G:A V306I 
TCGA-GU-A42R-01 chr4:1805533:C:T H349Y 
TCGA-UY-A78N-01 chr4:1805561:C:T P358L 
TCGA-CF-A5U8-01   
TCGA-E7-A3Y1-01 chr4:1806089:G:T G370C 
TCGA-E7-A7DU-01   
TCGA-BT-A42C-01   
TCGA-CF-A47W-01   
TCGA-DK-A6B0-01 chr4:1806099:A:G Y373C 
TCGA-E7-A5KF-01   
TCGA-DK-A3IS-01 chr4:1808937:C:-: A790fs 
Fusions 
TCGA-CF-A3MF-01   
TCGA-CF-A3MG-01 FGFR3-TACC3  
TCGA-CF-A3MH-01   
TCGA sample IDMutationAmino acid substitution
TCGA-DK-A3IS-01 chr4:1803377:G:A E216K 
TCGA-DK-A3IS-01 chr4:1803395:G:A D222N 
TCGA-DK-A3IS-01 chr4:1803435:G:A G235D 
TCGA-FJ-A3ZF-01 chr4:1803738:G:A V306I 
TCGA-GU-A42R-01 chr4:1805533:C:T H349Y 
TCGA-UY-A78N-01 chr4:1805561:C:T P358L 
TCGA-CF-A5U8-01   
TCGA-E7-A3Y1-01 chr4:1806089:G:T G370C 
TCGA-E7-A7DU-01   
TCGA-BT-A42C-01   
TCGA-CF-A47W-01   
TCGA-DK-A6B0-01 chr4:1806099:A:G Y373C 
TCGA-E7-A5KF-01   
TCGA-DK-A3IS-01 chr4:1808937:C:-: A790fs 
Fusions 
TCGA-CF-A3MF-01   
TCGA-CF-A3MG-01 FGFR3-TACC3  
TCGA-CF-A3MH-01   

Finally, we also discovered that PPARGC1A and PPARGC1B were exclusively mutated in non–T-cell-inflamed tumors. Both are coactivators of PPAR-γ, which was one of the top activated regulators in our gene expression pathway analysis. PPARGC1A mutations included A728D, C726Y, K663*, R657Q, R485W, and R408T. PPARGC1B mutations included P607L, V703I, E787K, and S994L. These mutations were evaluated by a prediction algorithm, SIFT, and 70% were predicted to affect protein function (17). Together, these data identify the β-catenin pathway, the PPAR-γ pathway, and active FGFR3 as three potential cancer-intrinsic pathways that could contribute to T-cell exclusion in bladder cancer.

Using our immune gene signature, we found that urothelial bladder cancer segregated into T-cell–inflamed and non–T-cell-inflamed subtypes, and the latter phenotype correlates with an absence of CD8+ T cells by IHC. Unsupervised clustering analyses of muscle-invasive bladder cancer based on global gene expression segregate tumors into two to four intrinsic subtypes (19–21), but previous studies did not evaluate segregation based on T-cell infiltration in bladder cancer, as we have done. Our strategy was based on the presence of immune cells in the tumor microenvironment, rather than clustering by tumor-derived genomic pathways. Although the latter approach might better identify chemotherapeutic resistance pathways, the former is more likely to yield insight into immunotherapy resistance. PPAR-γ pathway activation and FGFR3 mutations are associated with the “luminal” intrinsic bladder subtype, and we found both were also linked with the non–T-cell-inflamed phenotype. Further bioinformatic and experimental analyses are needed to determine the relationship of the non–T-cell-inflamed subtype and previously described bladder cancer–intrinsic subtypes.

Our use of immune gene expression profiling identified three molecular pathways linked with the non–T-cell-inflamed subtype. Our methodology selected for genes overexpressed in tumors with an absence of immune-infiltrating cells, which by inference were highly likely to arise from tumor-intrinsic activated molecular pathways. Modulation of these targets could prove useful to overcome primary resistance to immune checkpoint blockade such as PD-1/PD-L1–targeted therapies. One limitation to this analysis is its correlative nature. However, similar methods have led to the discovery of a causal link between activated β-catenin and intratumoral T-cell restriction via ATF3-mediated CCL4 repression in melanoma, which was, indeed, confirmed through careful experimental studies in genetically engineered mouse models (14). Our study independently validated this pathway as active in non–T-cell-inflamed bladder cancers. We also identified two other molecular pathways not previously found in melanoma. The effect of modulation of these pathways on the induction of a T-cell–inflamed phenotype will be an attractive area of research to pursue in future studies. The availability of FGFR3 inhibitors for clinical testing provides a rationale and feasibility for combination with monoclonal antibodies to PD-1 or PD-L1.

An understanding of the non–T-cell-inflamed phenotype in urothelial bladder, as well as other cancers, is critical to improving response rates for future immunotherapeutic development (22, 23). Three immune checkpoint inhibitors have been approved by the FDA since 2011 for melanoma, renal, and non–small cell lung cancers (23–25). It is critical to continue advancing our biologic understanding of tumors with primary resistance to expand the proportion of patients responding to these immunotherapies. A better understanding of resistance will catalyze the development of new clinical trial strategies with combination therapies that overcome resistance in urothelial bladder cancer and, ultimately, other tumor types.

W.M. Stadler reports receiving commercial research support from Johnson & Johnson and is a consultant/advisory board member for Astra-Zeneca. No potential conflicts of interest were disclosed by the other authors.

Conception and design: R.F. Sweis, S. Spranger, G. Steinberg, T.F. Gajewski

Development of methodology: R.F. Sweis, S. Spranger, R. Bao, T.F. Gajewski

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): R.F. Sweis, G.P. Paner, G. Steinberg

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): R.F. Sweis, S. Spranger, R. Bao, T.F. Gajewski

Writing, review, and/or revision of the manuscript: R.F. Sweis, R. Bao, W.M. Stadler, G. Steinberg, T.F. Gajewski

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): R.F. Sweis, G.P. Paner, W.M. Stadler

Study supervision: R.F. Sweis, G. Steinberg, T.F. Gajewski

The authors are extremely grateful for a generous donation from Diana and Andrew Da Miano supporting this research. The authors also thank Terri Li and the University of Chicago HTRC for their assistance with IHC.

This study was supported by NIH grants T32 CA009566 and T32 GM007019 (to R.F. Sweis).

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
Ferlay
J
,
Soerjomataram
II
,
Dikshit
R
,
Eser
S
,
Mathers
C
,
Rebelo
M
, et al
Cancer incidence and mortality worldwide: Sources, methods and major patterns in GLOBOCAN 2012
.
Int J Cancer
2015
;
136
:
E359
86
.
2.
von der Maase
H
,
Sengelov
L
,
Roberts
JT
,
Ricci
S
,
Dogliotti
L
,
Oliver
T
, et al
Long-term survival results of a randomized trial comparing gemcitabine plus cisplatin, with methotrexate, vinblastine, doxorubicin, plus cisplatin in patients with bladder cancer
.
J Clin Oncol
2005
;
23
:
4602
8
.
3.
Bracarda
S
,
Altavilla
A
,
Hamzaj
A
,
Sisani
M
,
Marrocolo
F
,
Del Buono
S
, et al
Immunologic checkpoints blockade in renal cell, prostate, and urothelial malignancies
.
Semin Oncol
2015
;
42
:
495
505
.
4.
Carosella
ED
,
Ploussard
G
,
LeMaoult
J
,
Desgrandchamps
F
. 
A Systematic review of immunotherapy in urologic cancer: Evolving roles for targeting of CTLA-4, PD-1/PD-L1, and HLA-G
.
Eur Urol
2015
;
68
:
267
79
.
5.
Kim
JW
,
Tomita
Y
,
Trepel
J
,
Apolo
AB
. 
Emerging immunotherapies for bladder cancer
.
Curr Opin Oncol
2015
;
27
:
191
200
.
6.
Powles
T
,
Eder
JP
,
Fine
GD
,
Braiteh
FS
,
Loriot
Y
,
Cruz
C
, et al
MPDL3280A (anti-PD-L1) treatment leads to clinical activity in metastatic bladder cancer
.
Nature
2014
;
515
:
558
62
.
7.
Harlin
H
,
Meng
Y
,
Peterson
AC
,
Zha
Y
,
Tretiakova
M
,
Slingluff
C
, et al
Chemokine expression in melanoma metastases associated with CD8+ T-cell recruitment
.
Cancer Res
2009
;
69
:
3077
85
.
8.
Ji
RR
,
Chasalow
SD
,
Wang
L
,
Hamid
O
,
Schmidt
H
,
Cogswell
J
, et al
An immune-active tumor microenvironment favors clinical response to ipilimumab
.
Cancer Immunol Immunother
2012
;
61
:
1019
31
.
9.
Spranger
S
,
Spaapen
RM
,
Zha
Y
,
Williams
J
,
Meng
Y
,
Ha
TT
, et al
Up-regulation of PD-L1, IDO, and T(regs) in the melanoma tumor microenvironment is driven by CD8(+) T cells
.
Sci Transl Med
2013
;
5
:
200ra116
.
10.
Gajewski
TF
,
Louahed
J
,
Brichard
VG
. 
Gene signature in melanoma associated with clinical activity: A potential clue to unlock cancer immunotherapy
.
Cancer J
2010
;
16
:
399
403
.
11.
Gajewski
TF
,
Schreiber
H
,
Fu
YX
. 
Innate and adaptive immune cells in the tumor microenvironment
.
Nat Immunol
2013
;
14
:
1014
22
.
12.
Gajewski
TF
,
Woo
SR
,
Zha
Y
,
Spaapen
R
,
Zheng
Y
,
Corrales
L
, et al
Cancer immunotherapy strategies based on overcoming barriers within the tumor microenvironment
.
Curr Opin Immunol
2013
;
25
:
268
76
.
13.
Sharma
P
,
Shen
Y
,
Wen
S
,
Yamada
S
,
Jungbluth
AA
,
Gnjatic
S
, et al
CD8 tumor-infiltrating lymphocytes are predictive of survival in muscle-invasive urothelial carcinoma
.
Proc Natl Acad Sci U S A
2007
;
104
:
3967
72
.
14.
Spranger
S
,
Bao
R
,
Gajewski
TF
. 
Melanoma-intrinsic β-catenin signalling prevents anti-tumour immunity
.
Nature
2015
;
523
:
231
5
.
15.
Li
B
,
Dewey
CN
. 
RSEM: Accurate transcript quantification from RNA-Seq data with or without a reference genome
.
BMC Bioinformatics
2011
;
12
:
323
.
16.
Wilkerson
MD
,
Hayes
DN
. 
ConsensusClusterPlus: A class discovery tool with confidence assessments and item tracking
.
Bioinformatics
2010
;
26
:
1572
3
.
17.
Kumar
P
,
Henikoff
S
,
Ng
PC
. 
Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm
.
Nat Protoc
2009
;
4
:
1073
81
.
18.
Iyer
G
,
Milowsky
MI
. 
Fibroblast growth factor receptor-3 in urothelial tumorigenesis
.
Urol Oncol
2013
;
31
:
303
11
.
19.
Choi
W
,
Porten
S
,
Kim
S
,
Willis
D
,
Plimack
ER
,
Hoffman-Censits
J
, et al
Identification of distinct basal and luminal subtypes of muscle-invasive bladder cancer with different sensitivities to frontline chemotherapy
.
Cancer Cell
2014
;
25
:
152
65
.
20.
Cancer Genome Atlas Research Network
. 
Comprehensive molecular characterization of urothelial bladder carcinoma
.
Nature
2014
;
507
:
315
22
.
21.
Damrauer
JS
,
Hoadley
KA
,
Chism
DD
,
Fan
C
,
Tiganelli
CJ
,
Wobker
SE
, et al
Intrinsic subtypes of high-grade bladder cancer reflect the hallmarks of breast cancer biology
.
Proc Natl Acad Sci U S A
2014
;
111
:
3110
5
.
22.
Gajewski
TF
. 
The next hurdle in cancer immunotherapy: Overcoming the non-T cell-inflamed tumor microenvironment
.
Semin Oncol
2015
;
42
:
663
71
.
23.
Hamid
O
,
Robert
C
,
Daud
A
,
Hodi
FS
,
Hwu
WJ
,
Kefford
R
, et al
Safety and tumor responses with lambrolizumab (anti-PD-1) in melanoma
.
N Engl J Med
2013
;
369
:
134
44
.
24.
Hodi
FS
,
O'Day
SJ
,
McDermott
DF
,
Weber
RW
,
Sosman
JA
,
Haanen
JB
, et al
Improved survival with ipilimumab in patients with metastatic melanoma
.
N Engl J Med
2010
;
363
:
711
23
.
25.
Brahmer
J
,
Reckamp
KL
,
Baas
P
,
Crinò
L
,
Eberhardt
WE
,
Poddubskaya
E
, et al
Nivolumab versus docetaxel in advanced squamous-cell non-small-cell lung cancer
.
N Engl J Med
2015
;
373
:
123
35
.