AUTHOR=Hu Rixin , Tao Tao , Yu Lu , Ding Qiuxia , Zhu Guanghui , Peng Guoyu , Zheng Shiwen , Yang Leyun , Wu Song TITLE=Multi-Omics Characterization of Tumor Microenvironment Heterogeneity and Immunotherapy Resistance Through Cell States–Based Subtyping in Bladder Cancer JOURNAL=Frontiers in Cell and Developmental Biology VOLUME=Volume 9 - 2021 YEAR=2022 URL=https://www.frontiersin.org/journals/cell-and-developmental-biology/articles/10.3389/fcell.2021.809588 DOI=10.3389/fcell.2021.809588 ISSN=2296-634X ABSTRACT=Due to the strong heterogeneity of bladder cancer (BC), there is often a substantial variation in the prognosis and efficiency of immunotherapy among BC patients. For the precision treatment and assessment of prognosis, the subtyping of BC plays a critical role. Despite various subtyping methods proposed previously, most of them are based on a limited number of molecules and none of them is developed on the basis of cell states. In this study, by integrating single-cell RNA-seq and bulk RNA-seq data, we constructed a single-cell atlas and identified the absolute proportion of 22 different cell states in BC, including immune and nonimmune cell states derived from tumor tissues. To explore the heterogeneity of BC, BC was identified into four different subtypes in multiple cohorts using an improved consensus clustering algorithm based on cell states. Among the four subtypes, C1 had median prognosis and best overall response rate (ORR), which characterized an immunosuppressive tumor microenvironment (TME). C2 was enriched in EMT/invasion, angiogenesis, immunosuppression, and immune exhaustion. Surely C2 performed worst in prognosis and ORR. C3 with worse ORR than C2 was enriched in angiogenesis and almost nonimmune exhaustion. Displaying an immune effective environment, C4 performed best in prognosis and ORR. We found that patients with just immunosuppressive environments are suitable for immunotherapy, but patients with immunosuppressive environments accompanied by immune exhaustion or angiogenesis may resist to immunotherapy. Furthermore, we conducted exploration into the heterogeneity of transcriptome, mutational profiles, and somatic copy-number alterations (SCNAs) in four subtypes, which could explain the significant differences in prognosis and ORR related to cell states. We also found PD-1 in immune cells and tumor cells could both influence ORR in BC. The level of TGFβ in a cell state can be opposite to the overall level of tissues and the level in a specific cell state could predict ORR more accurately. Thus, our work will further the understanding of the heterogeneity and immunotherapy resistance in BC which is expected to assist clinical practice and serve as a supplement to the current subtyping method from a novel perspective of cell states.