Identification of Genes and Construction of Prognostic Model of Lung Adenocarcinoma Based on Propionate Metabolism-Related Genes

TL;DR

This study identifies genes related to propionate metabolism that help predict the survival of lung adenocarcinoma patients.

Contribution

A novel prognostic model using 16 propionate metabolism-related genes for predicting LUAD patient outcomes is developed.

Findings

The model showed AUC values of 0.722, 0.696, and 0.700 for 1-, 3-, and 5-year overall survival.

GSVA and GSEA revealed significant enrichment in G2/M checkpoint, glycolysis, HIF-1, IL-17, and p53 pathways.

qRT-PCR confirmed differential expression of ADIPOQ, CYP27A1, GCDH, and SERPINE1 between LUAD and normal cells.

Abstract

Dysregulation of propionate metabolism can enhance the invasive properties of lung adenocarcinoma (LUAD) cells and increase their metastatic potential. Therefore, we constructed a predictive model based on propionate metabolism-related genes (PMRGs) to evaluate the prognosis of patients with LUAD.

mRNA expression profiles and clinical data of LUAD patients were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). The predictive model was constructed using least absolute shrinkage and selection operator (LASSO). The associations between the risk score and tumor prognosis, immune infiltration, drug sensitivity, signaling pathways, and clinical features were evaluated using the CIBERSORT algorithm, the Genomics of Drug Sensitivity in Cancer (GDSC) database, gene set variation analysis (GSVA), gene set enrichment analysis (GSEA), and nomogram analyses. For hub genes, motif enrichment, genome-wide association study (GWAS) analysis and single-cell analysis were performed. Expression differences between LUAD cells and normal lung epithelial cells were validated using quantitative real-time polymerase chain reaction (qRT-PCR).

We established a LUAD prognostic model containing 16 hub genes. The area under the receiver operating characteristic curve (AUC) values for the model were 0.722 (1-year overall survival (OS)), 0.696 (3-year OS), and 0.700 (5-year OS). GSVA revealed significant enrichment in the G2/M checkpoint, E2F targets, and glycolysis pathways. GSEA indicated enrichment of the hypoxia-inducible factor-1 (HIF-1), interleukin (IL)-17, and p53 signaling pathways. Transcription factor analysis identified the motif cisbp-M4287 as the most significantly enriched, with a normalized enrichment score (NES) of 5.9. qRT-PCR results showed that the expression of ADIPOQ, CYP27A1, and GCDH was downregulated, while the expression of SERPINE1 was upregulated. No statistically significant differences were observed in the expression of CYP17A1 and EHHADH.

A novel prognostic model based on PMRGs was established to predict overall survival in LUAD patients.

Figures