Dataset Information

An investigation to identify tumor microenvironment-related genes of prognostic value in lung squamous cell carcinoma based on The Cancer Genome Atlas

ABSTRACT:

Background

Lung squamous cell carcinoma (LUSC) is a prevalent and lethal malignancy with a poor clinical prognosis. Major constituents of the tumor microenvironment (TME) include infiltrating immune cells and stromal cells, which play a pivotal role in the progression and growth of the disease. To improve the understanding of the prognostic influence of immune and stromal cell-related genes for patients with the disease, we performed a comprehensive bioinformatics analysis to identify TME-relevant biomarkers, and investigated the potential role of these candidate signatures in LUSC.

Methods

Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE) assessed the samples of LUSC obtained from The Cancer Genome Atlas (TCGA). The samples were grouped according to their immune/stromal scores (high or low). Multivariate cox regression and receiver operating characteristic curves (ROC) were implemented to construct the risk assessment model for prognosis prediction. The co-upregulated differentially expressed genes (DEGs) in the immune and stromal groups were used for further analyses. Overall survival (OS) curves were used to determine the prognostic value of the DEGs, and the TME-related DEGs were verified with Gene Expression Omnibus (GEO) datasets. The functional assessments were performed include Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and protein–protein interaction (PPI) analyses.

Results

The immune/stromal scores calculated by ESTIMATE showed significant associations with OS (log-rank P<0.05). In addition, the prognostic risk score model based on immune and stromal scores also showed significant correlations with OS. A total of 94 TME-related genes were obviously related to poor OS. Among them, BHMT2, FES, HSPB7, NOVA2, LPAP2, and SEMA3B (BFHNLS) were confirmed using GSE4573 and GSE17710 datasets. The functional assessments exhibited those TME-related genes mostly participate in immune response, cytokine-cytokine receptor interaction, and metabolic pathways, which elucidated the probable correlation of TME with tumorigenesis in LUSC.

Conclusions

In this study, 6 potential biomarkers named BFHNLS were identified as TME-related genes with prognostic value based on immune and stromal scores of LUSC patients of TCGA, and were verified using GEO datasets, which might serve as therapeutic targets.

SUBMITTER: Liu H

PROVIDER: S-EPMC8798203 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:The host's immune system plays a pivotal role in many tumor types, including squamous cell carcinomas (SCCs). We aim to identify immunological prognosticators for lymph node metastases (LNM) and disease-specific survival (DSS) in penile SCC. For this retrospective observational cohort study, penile SCC patients (n?=?213) treated in the Netherlands Cancer Institute, were selected if sufficient formalin-fixed, paraffin-embedded tumor material was available. Analysis included previously described high-risk human papilloma virus (hrHPV) status, immunohistochemical scores for classical and non-classical human leukocyte antigen (HLA) class I, programmed death ligand-1 (PD-L1) expression, and novel data on tumor-infiltrating macrophages and cytotoxic an regulatory T-cells. Clinicopathological characteristics and extended follow-up were also included. Regression analyses investigated relationships of the immune parameters with LNM and DSS. In the total cohort, diffuse PD-L1 tumor-cell expression, CD163+ macrophage infiltration, non-classical HLA class I upregulation, and low stromal CD8+ T-cell infiltration were all associated with LNM. In the multivariable model, only tumor PD-L1 expression remained a significant predictor for LNM (odds ratio (OR) 2.8, p?=?0.05). hrHPV negativity and diffuse PD-L1 tumor-cell expression were significantly associated with poor DSS and remained so upon correction for clinical parameters [hazard ratio (HR) 9.7, p?<?0.01 and HR 2.8, p?=?0.03]. The only immune factor with different expression in HPV+ and HPV- tumors was PD-L1, with higher PD-L1 expression in the latter (p?=?0.03). In the HPV- cohort (n?=?158), LNM were associated with diffuse PD-L1 tumor-cell expression, high intratumoral CD163+ macrophage infiltration, and low number of stromal CD8+ T-cells. The first two parameters were also linked to DSS. In the multivariable regression model, diffuse PD-L1 expression remained significantly unfavorable for DSS (HR 5.0, p?<?0.01). These results emphasize the complexity of the tumor microenvironment in penile cancer and point toward several possible immunotherapy targets. Here described immune factors can aid risk-stratification and should be evaluated in clinical immunotherapy studies to ultimately lead to patient tailored treatment.