Decoding tumor immune microenvironment heterogeneity by single-cell and spatial multi-omics: From immunotherapy resistance to translational biomarkers

Immune checkpoint blockade has transformed cancer therapy, yet primary and acquired resistance remain major clinical challenges. Increasing evidence indicates that immunotherapy resistance cannot be fully explained by tumor-intrinsic alterations or conventional biomarkers such as PD-L1 expression, tumor mutational burden, or microsatellite instability. Instead, therapeutic response is shaped by the tumor immune microenvironment (TIME) as a heterogeneous, spatially organized, and dynamically evolving ecosystem. Single-cell omics has revealed diverse immune and stromal cell states, including progenitor and terminally exhausted T cells, suppressive myeloid programs, B-cell/TLS-associated immune-reactive states, and CAF-mediated exclusion phenotypes. Spatial transcriptomics, spatial proteomics, and imaging-based approaches further demonstrate that these cell states assemble into distinct immune niches, including immune-inflamed, T-cell-excluded, myeloid-suppressive, metabolic/hypoxic, and TLS-associated niches. These spatial ecosystems determine whether antitumor immune cells can access malignant cells, receive antigen-presenting support, or become restrained by stromal, vascular, metabolic, and myeloid barriers. In this review, we summarize how single-cell and spatial multi-omics redefine TIME heterogeneity in immunotherapy resistance, highlight ligand-receptor communication networks linking cell states to spatial immune dysfunction, and discuss emerging translational biomarkers for patient stratification. We further propose that future immunotherapy biomarkers should evolve from static single-marker assays toward longitudinal, spatially resolved, and interpretable multi-omics models that guide precision combination immunotherapy.