To define dominant molecular and cellular features associated with PD-1/PD-L1 blockade resistance in metastatic urothelial cancer.

We pursued an unbiased approach using bulk RNA sequencing data from two clinical trials to discover (IMvigor 210) and validate (CheckMate 275) pretreatment molecular features associated with resistance to PD-1/PD-L1 blockade in metastatic urothelial cancer. We then generated single-cell RNA sequencing (scRNA-seq) data from muscle-invasive bladder cancer specimens to dissect the cellular composition underlying the identified gene signatures.

We identified an adaptive immune response gene signature associated with response and a protumorigenic inflammation gene signature associated with resistance to PD-1/PD-L1 blockade. The adaptive immune response:protumorigenic inflammation signature expression ratio, coined the 2IR score, best correlated with clinical outcomes, and was externally validated. Mapping these bulk gene signatures onto scRNA-seq data uncovered their underlying cellular diversity, with prominent expression of the protumorigenic inflammation signature by myeloid phagocytic cells. However, heterogeneity in expression of adaptive immune and protumorigenic inflammation genes was observed among single myeloid phagocytic cells, quantified as the myeloid single cell immune:protumorigenic inflammation ratio (M_sc2IR) score. Single myeloid phagocytic cells with low M_sc2IR scores demonstrated upregulation of proinflammatory cytokines/chemokines and downregulation of antigen presentation genes, were unrelated to M1 versus M2 polarization, and were enriched in pretreatment blood samples from patients with PD-L1 blockade–resistant metastatic urothelial cancer.

The balance of adaptive immunity and protumorigenic inflammation in individual tumor microenvironments is associated with PD-1/PD-L1 resistance in urothelial cancer with the latter linked to a proinflammatory cellular state of myeloid phagocytic cells detectable in tumor and blood.

See related commentary by Drake, p. 4139