Characterization of the Immune Microenvironmental Landscape of Lung Squamous Cell Carcinoma with Immune Cell Infiltration

Background. Increasing evidence supports that immune cell infiltration (ICI) patterns play a key role in the tumor progression of lung squamous cell carcinoma (LUSC). However, to date, the immune infiltration picture of LUSC has not been elucidated. Method. TCGA was used to download multiomics data from LUSC samples. At the same time, we included two datasets on lung squamous cell carcinoma, GSE17710 and GSE157010. To reveal the landscape of tumor immune microenvironment (TIME), the ESTIMATE algorithm, ssGSEA approach, and CIBERSORT analysis are used. To quantify the ICI pattern in a single tumor, consistent clustering is used to determine the LUSC subtype based on the ICI pattern, and principal component analysis (PCA) is used to obtain the ICI score. The prognostic value of the Kaplan-Meier curves is confirmed. GSEA (Gene Set Enrichment Analysis) was used to perform functional annotation. To investigate the immunotherapeutic effects of the ICI score, the immunophenotyping score (IPS) is used. Finally, analyze the mutation data with the “maftools” R package. Results. We identified four different immune infiltration patterns with different prognosis and biological characteristics in 792 LUSC samples. The identification of ICI patterns in individual tumors developed under ICI-related characteristic genes based on the ICI score helps to analyze the biological process, clinical results, immune cell infiltration, immunotherapy effects, and genetic variation. Immune failure is indicated by a high ICI score subtype marked by immunosuppression. Patients with low ICI scores have an abundance of efficient immune cells, which corresponds to the immunological activation phenotype and may have therapeutic benefits. The immunophenotypic score was used as a surrogate indicator of immunotherapy results, and samples with low ICI scores obtained significantly higher immunophenotypic scores. Finally, the relationship between the ICI score and tumor mutation burden (TMB) was proven. Conclusion. This study fully clarified the indispensable role of the ICI model in the complexity and diversity of TIME. The quantitative identification of ICI patterns in a single tumor will help draw the picture of TIME and further optimize precision immunotherapy.