Cellular dynamics in tumour microenvironment along with lung cancer progression underscore spatial and evolutionary heterogeneity of neutrophil.

BACKGROUND: The cellular dynamics in the tumour microenvironment (TME) along with non-small cell lung cancer (NSCLC) progression remain unclear. METHODS: Multiplex immunofluorescence test detecting 10 immune-related markers on 553 primary tumour (PT) samples of NSCLC was conducted and spatial information in TME was assessed by the StarDist depth learning model. The single-cell transcriptomic atlas of PT (n = 4) and paired tumour-draining lymph nodes (TDLNs) (n = 5 for tumour-invaded, n = 3 for tumour-free) microenvironment was profiled. Various bioinformatics analyses based on Gene Expression Omnibus, TCGA and Array-Express databases were also used to validate the discoveries. RESULTS: Spatial distances of CD4+ T cells-CD38+ T cells, CD4+ T cells-neutrophils and CD38+ T cells-neutrophils prolonged and they were replaced by CD163+ macrophages in PT along with tumour progression. Neutrophils showed unique stage and location-dependent prognostic effects. A high abundance of stromal neutrophils improved disease-free survival in the early-stage, whereas high intratumoural neutrophil infiltrates predicted poor prognosis in the mid-to-late-stage. Significant molecular and functional reprogramming in PT and TDLN microenvironments was observed. Diverse interaction networks mediated by neutrophils were found between positive and negative TDLNs. Five phenotypically and functionally heterogeneous subtypes of tumour-associated neutrophil (TAN) were further identified by pseudotime analysis, including TAN-0 with antigen-presenting function, TAN-1 with strong expression of interferon (IFN)-stimulated genes, the pro-tumour TAN-2 subcluster, the classical subset (TAN-3) and the pro-inflammatory subtype (TAN-4). Loss of IFN-stimulated signature and growing angiogenesis activity were discovered along the transitional trajectory. Eventually, a robust six neutrophil differentiation relevant genes-based model was established, showing that low-risk patients had longer overall survival time and may respond better to immunotherapy. CONCLUSIONS: The cellular composition, spatial location, molecular and functional changes in PT and TDLN microenvironments along with NSCLC progression were deciphered, highlighting the immunoregulatory roles and evolutionary heterogeneity of TANs.

Multiplex immunofluorescence and single-cell transcriptomic profiling reveal the spatial cell interaction networks in the non-small cell lung cancer microenvironment.

BACKGROUND: Conventional immunohistochemistry technologies were limited by the inability to simultaneously detect multiple markers and the lack of identifying spatial relationships among cells, hindering understanding of the biological processes in cancer immunology. METHODS: Tissue slices of primary tumours from 553 IA∼IIIB non-small cell lung cancer (NSCLC) cases were stained by multiplex immunofluorescence (mIF) assay for 10 markers, including CD4, CD38, CD20, FOXP3, CD66b, CD8, CD68, PD-L1, CD133 and CD163, evaluating the amounts of 26 phenotypes of cells in tumour nest and tumour stroma. StarDist depth learning model was utilised to determine the spatial location of cells based on mIF graphs. Single-cell RNA sequencing (scRNA-seq) on four primary NSCLC cases was conducted to investigate the putative cell interaction networks. RESULTS: Spatial proximity among CD20+ B cells, CD4+ T cells and CD38+ T cells (r2  = 0.41) was observed, whereas the distribution of regulatory T cells was associated with decreased infiltration levels of CD20+ B cells and CD38+ T cells (r2  = -0.45). Univariate Cox analyses identified closer proximity between CD8+ T cells predicted longer disease-free survival (DFS). In contrast, closer proximity between CD133+ cancer stem cells (CSCs), longer distances between CD4+ T cells and CD20+ B cells, CD4+ T cells and neutrophils, and CD20+ B cells and neutrophils were correlated with dismal DFS. Data from scRNA-seq further showed that spatially adjacent N1-like neutrophils could boost the proliferation and activation of T and B lymphocytes, whereas spatially neighbouring M2-like macrophages showed negative effects. An immune-related risk score (IRRS) system aggregating robust quantitative and spatial prognosticators showed that high-IRRS patients had significantly worse DFS than low-IRRS ones (HR 2.72, 95% CI 1.87-3.94, p < .001). CONCLUSIONS: We developed a framework to analyse the cell interaction networks in tumour microenvironment, revealing the spatial architecture and intricate interplays between immune and tumour cells.