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

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

Methods: 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.

Results: 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.

Conclusion: 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.