Integrated multi-omic analysis and experiment reveals the role of endoplasmic reticulum stress in lung adenocarcinoma

Background Lung cancer is a highly prevalent malignancy worldwide and is associated with high mortality rates. While the involvement of endoplasmic reticulum (ER) stress in the development of lung adenocarcinoma (LUAD) has been established, the underlying mechanism remains unclear. Methods In this study, we utilized data from The Cancer Genome Atlas (TCGA) to identify differentially expressed endoplasmic reticulum stress-related genes (ERSRGs) between LUAD and normal tissues. We performed various bioinformatics analyses to investigate the biological functions of these ERSRGs. Using LASSO analysis and multivariate stepwise regression, we constructed a novel prognostic model based on the ERSRGs. We further validated the performance of the model using two independent datasets from the Gene Expression Omnibus (GEO). Additionally, we conducted functional enrichment analysis, immune checkpoint analysis, and immune infiltration analysis and drug sensitivity analysis of LUAD patients to explore the potential biological function of the model. Furthermore, we conducted a battery of experiments to verify the expression of ERSRGs in a real-world cohort. Results We identified 106 ERSRGs associated with LUAD, which allowed us to classify LUAD patients into two subtypes based on gene expression differences. Using six prognostic genes (NUPR1, RHBDD2, VCP, BAK1, EIF2AK3, MBTPS2), we constructed a prognostic model that exhibited excellent predictive performance in the training dataset and was successfully validated in two independent external datasets. The risk score derived from this model emerged as an independent prognostic factor for LUAD. Confirmation of the linkage between this risk model and immune infiltration was affirmed through the utilization of Gene Set Enrichment Analysis (GSEA), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. The q-PCR results verified significant differences in the expression of prognostic genes between cancer and paracancer tissues. Notably, the protein expression of NUPR1, as determined by immunohistochemistry (IHC), exhibited an opposite pattern compared to the mRNA expression patterns. Conclusion This study establishes a novel prognostic model for LUAD based on six ER stress-related genes, facilitating the prediction of LUAD prognosis. Additionally, NUPR1 was identified as a potential regulator of stress in LUAD.