Intercellular Molecular Crosstalk Networks within Invasive and Immunosuppressive Tumor Microenvironment Subtypes Associated with Clinical Outcomes in Four Cancer Types

Heterogeneity is a critical basis for understanding how the tumor microenvironment (TME) contributes to tumor progression. However, an understanding of the specific characteristics and functions of TME subtypes (subTMEs) in the progression of cancer is required for further investigations into single-cell resolutions. Here, we analyzed single-cell RNA sequencing data of 250 clinical samples with more than 200,000 cells analyzed in each cancer datum. Based on the construction of an intercellular infiltration model and unsupervised clustering analysis, four, three, three, and four subTMEs were revealed in breast, colorectal, esophageal, and pancreatic cancer, respectively. Among the subTMEs, the immune-suppressive subTME (subTME-IS) and matrix remodeling with malignant cells subTME (subTME-MRM) were highly enriched in tumors, whereas the immune cell infiltration subTME (subTME-ICI) and precancerous state of epithelial cells subTME (subTME-PSE) were less in tumors, compared with paracancerous tissues. We detected and compared genes encoding cytokines, chemokines, cytotoxic mediators, PD1, and PD-L1. The results showed that these genes were specifically overexpressed in different cell types, and, compared with normal tissues, they were upregulated in tumor-derived cells. In addition, compared with other subTMEs, the expression levels of PDCD1 and TGFB1 were higher in subTME-IS. The Cox proportional risk regression model was further constructed to identify possible prognostic markers in each subTME across four cancer types. Cell-cell interaction analysis revealed the distinguishing features in molecular pairs among different subTMEs. Notably, ligand–receptor gene pairs, including COL1A1-SDC1, COL6A2-SDC1, COL6A3-SDC1, and COL4A1-ITGA2 between stromal and tumor cells, associated with tumor invasion phenotypes, poor patient prognoses, and tumor advanced progression, were revealed in subTME-MRM. C5AR1-RPS19, LGALS9-HAVCR2, and SPP1-PTGER4 between macrophages and CD8+ T cells, associated with CD8+ T-cell dysfunction, immunosuppressive status, and tumor advanced progression, were revealed in subTME-IS. The spatial co-location information of cellular and molecular interactions was further verified by spatial transcriptome data from colorectal cancer clinical samples. Overall, our study revealed the heterogeneity within the TME, highlighting the potential pro-invasion and pro-immunosuppressive functions and cellular infiltration characteristics of specific subTMEs, and also identified the key cellular and molecular interactions that might be associated with the survival, invasion, immune escape, and classification of cancer patients across four cancer types.