Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Immunol.

Sec. Cancer Immunity and Immunotherapy

This article is part of the Research TopicUrothelial Neoplasms: An Integrated Approach to Prevention, Diagnostics, and Personalized TherapyView all 18 articles

Revealing Intra-Group Immunotherapy Response Heterogeneity in Metastatic Urothelial Carcinoma through Interpretable Feature Extraction and Spectral Clustering

Provisionally accepted
Yoshiyuki  NagumoYoshiyuki NagumoXiucai  YeXiucai YeTianyi  ShiTianyi ShiBryan  J MathisBryan J MathisTetsuya  SakuraiTetsuya SakuraiHiroyuki  NishiyamaHiroyuki Nishiyama*
  • University of Tsukuba, Tsukuba, Japan

The final, formatted version of the article will be published soon.

Introduction: Immune checkpoint inhibitors (ICIs) have improved outcomes in metastatic urothelial carcinoma (mUC) but clinical responses remain highly heterogenous. Traditional binary classification of response overlooks clinically relevant variability within each group but a more detailed understanding of intra-group heterogeneity may support subclass-specific therapeutic strategies. Methods: We developed a novel analysis framework that integrates interpretable feature extraction and spectral clustering to identify patient subclasses associated with heterogeneous responses to ICIs. This method was applied to tumor transcriptomic data from the IMvigor210 cohort (n = 298), comprising mUC patients treated with atezolizumab. Interpretable features based on SHapley Additive exPlanations (SHAP) were computed from a response classification model to quantify patient-level gene contributions, which were then used for spectral clustering. An independent cohort (GSE176307, n = 88) was used for external validation. Results: This approach identified four patient clusters with distinct immune phenotypes and response patterns. Cluster 3 (92.3% responders) showed an inflamed phenotype with high PD-L1 expression, T cell activation, and TP53 mutations. Cluster 1 (100% non-responders) displayed an immune-desert phenotype with FGFR3 mutations and elevated TGF-β signaling. Cluster 2 was more heterogeneous, containing two subgroups (Sub 1 and Sub 2) with differing immune activity and immunosuppressive gene expression, corresponding to response rates of 23.2% and 77.3%, respectively. Similar patterns were observed in the validation cohort. Conclusions: Our framework, which combines SHAP-based interpretable feature extraction with spectral clustering, revealed subclass-level heterogeneity in ICI response, highlighting biologically distinct immune subclasses. This approach may facilitate the development of subclass-specific therapeutic strategies.

Keywords: urothelial carcinoma, Immune checkpoint inhibitor, response, biomarker, Gene Expression, Shap, spectral clustering

Received: 15 May 2025; Accepted: 11 Dec 2025.

Copyright: © 2025 Nagumo, Ye, Shi, Mathis, Sakurai and Nishiyama. 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) or licensor 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: Hiroyuki Nishiyama

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.