Transcriptome Analyses Identify a Metabolic Gene Signature Indicative of Antitumor Immunosuppression of EGFR Wild Type Lung Cancers With Low PD-L1 Expression

Abstract

PurposeWith the development and application of targeted therapies like tyrosine kinase inhibitors (TKIs) and immune checkpoint inhibitors (ICIs), non-small cell lung cancer (NSCLC) patients have achieved remarkable survival benefits in recent years. However, epidermal growth factor receptor (EGFR) wild-type and low expression of programmed death-ligand 1 (PD-L1) NSCLCs remain unmanageable. Few treatments for these patients exist, and more side effects with combination therapies have been observed. We intended to generate a metabolic gene signature that could successfully identify high-risk patients and reveal its underlying molecular immunology characteristics.MethodsBy identifying the bottom 50% PD-L1 expression level as PD-L1 low expression and removing EGFR mutant samples, a total of 640 lung adenocarcinoma (LUAD) and lung squamous carcinoma (LUSC) tumor samples and 93 adjacent non-tumor samples were finally extracted from The Cancer Genome Atlas (TCGA). We identified differentially expressed metabolic genes (DEMGs) by R package limma and the prognostic genes by Univariate Cox proportional hazards regression analyses. The intersect genes between DEMGs and prognostic genes were put into the least absolute shrinkage and selection operator (LASSO) penalty Cox regression analysis. The metabolic gene signature contained 18 metabolic genes generated and successfully stratified LUAD and LUSC patients into the high-risk and low-risk groups, which was also validated by the Gene Expression Omnibus (GEO) database. Its accuracy was proved by the time-dependent Receiver Operating Characteristic (ROC) curve, Principal Components Analysis (PCA), and nomogram. Furthermore, the Single-sample Gene Set Enrichment Analysis (ssGSEA) and diverse acknowledged methods include XCELL, TIMER, QUANTISEQ, MCPcounter, EPIC, CIBERSORT-ABS, and CIBERSORT revealed its underlying antitumor immunosuppressive status. Besides, its relationship with somatic copy number alterations (SCNAs) and tumor mutational burden (TMB) was also discussed.ResultsIt is noteworthy that metabolism reprogramming is associated with the survival of the double-negative LUAD and LUSC patients. The SCNAs and TMB of critical metabolic genes can inhibit the antitumor immune process, which might be a promising therapeutic target.

Prognostic characteristics of a six-gene signature based on ssGSEA in sarcoma

“…Considering TIME’s importance for tumor cell growth, we correlated the risk score with TIME. Firstly, the correlation between risk score and immune cell infiltration status was examined [ 23 ]. There are a number of immune cells that are involved in the analysis, including aDC, B cells, cytotoxic cells, DC, eosinophils, iDC, mast cells, NK cells, pDC, T cells, Tcm, Tgd, and Th2 cells.…”

Transcriptome analysis reveals the prognostic and immune infiltration characteristics of glycolysis and hypoxia in head and neck squamous cell carcinoma

“…Therefore, we conducted a correlation analysis between the risk score and TIME. Firstly, multiple analysis methods including TIMER, CIBERSORT, CIBERSORT-ABS, QUANTISEQ, MCPCOUNTER, XCELL and EPIC were used to analyze the correlations between risk scores and immune cell infiltration status [ 24 ]. The immune cells involved in the analysis included CD4 T cells, CD8 T cells, B cells, neutrophils, and so on.…”