Construction of an endoplasmic reticulum stress-related signature in lung adenocarcinoma by comprehensive bioinformatics analysis

Lung Adenocarcinoma (LUAD) is a major component of lung cancer. Endoplasmic reticulum stress (ERS) has emerged as a new target for some tumor treatments.

The expression and clinical data of LUAD samples were downloaded from The Cancer Genome Atlas (TCGA) and The Gene Expression Omnibus (GEO) database, followed by acquiring ERS-related genes (ERSGs) from the GeneCards database. Differentially expressed endoplasmic reticulum stress-related genes (DE-ERSGs) were screened and used to construct a risk model by Cox regression analysis. Kaplan-Meier (K-M) curves and receiver operating characteristic (ROC) curves were plotted to determine the risk validity of the model. Moreover, enrichment analysis of differentially expressed genes (DEGs) between the high- and low- risk groups was conducted to investigate the functions related to the risk model. Furthermore, the differences in ERS status, vascular-related genes, tumor mutation burden (TMB), immunotherapy response, chemotherapy drug sensitivity and other indicators between the high- and low- risk groups were studied. Finally, quantitative real-time polymerase chain reaction (qRT-PCR) was used to validate the mRNA expression levels of prognostic model genes.

A total of 81 DE-ERSGs were identified in the TCGA-LUAD dataset, and a risk model, including HSPD1, PCSK9, GRIA1, MAOB, COL1A1, and CAV1, was constructed by Cox regression analysis. K-M and ROC analyses showed that the high-risk group had a low survival, and the Area Under Curve (AUC) of ROC curves of 1-, 3- and 5-years overall survival was all greater than 0.6. In addition, functional enrichment analysis suggested that the risk model was related to collagen and extracellular matrix. Furthermore, differential analysis showed vascular-related genes FLT1, TMB, neoantigen, PD-L1 protein (CD274), Tumor Immune Dysfunction and Exclusion (TIDE), and T cell exclusion score were significantly different between the high- and low-risk groups. Finally, qRT-PCR results showed that the mRNA expression levels of 6 prognostic genes were consistent with the analysis.

A novel ERS-related risk model, including HSPD1, PCSK9, GRIA1, MAOB, COL1A1, and CAV1, was developed and validated, which provided a theoretical basis and reference value for ERS-related fields in the study and treatment of LUAD.