Multi-Omics Reveals Activated Fibroblasts with Dual Functional Roles in Repair and Negative Barrier in the CD8 + T Cell Cytotoxic Niche for Lung Cancer Neoadjuvant Therapy

ABSTRACT A major challenge in elucidating immune activation and tolerance is that single-omics technologies are inherently limited. Since single-cell transcriptomics lacks spatial information and spatial transcriptomics lacks resolution or depth, neither can adequately model the multidimensional tumor immune microenvironment (TiME) features. We present NICHE (Niche Integrated Cellular Heterogeneity Elucidator), an integrative multi-omics framework. By aligning physical cell-cell interactions with single-cell ligand-receptor (L-R) expression, NICHE leverages single-cell transcriptomics, spatial proteomics, and AI modeling to systematically decode the composition, spatial interactions, and communications within immune functional niches. To validate NICHE’s capabilities, we applied it to human tonsil tissue. The framework successfully identified key structural and functional units and deeply elucidated the cell-cell interactions and molecular crosstalk within them. Furthermore, these results were validated through orthogonal methods, including single-cell spatial transcriptomics and multiplexed immunofluorescence, which confirmed the accuracy and reproducibility of our analyses. In non-small cell lung cancer (NSCLC) patients undergoing neoadjuvant immunotherapy, CD8 + T cell cytotoxic niches mediate tumor clearance with variable outcomes, to define the cellular and molecular drivers behind this functional spectrum, we utilized NICHE to analyze the TiME from patients with complete versus partial pathological response, both pre- and post-treatment. We found that neoadjuvant immunotherapy may induce the activation of CXCL12 + , NECTIN2 + , POSTN + , and COL6A1 + fibroblasts through the activation of cytotoxic T cells. Paradoxically, these fibroblasts support tissue repair but also establish an immunosuppressive niche that protects residual tumor cells, ultimately attenuating treatment efficacy. Our work establishes a multidimensional framework for dissecting dynamic immune activity, with direct implications for understanding tumor immunology and identifying novel therapeutic targets.