An Oxidative Stress-Associated Seven-Gene Prognostic Signature in Lung Adenocarcinoma: Integrative Transcriptomic Analysis Across Public Cohorts

Abstract

Lung adenocarcinoma is molecularly heterogeneous, and oxidative-stress programs can support either tumor restraint or tumor adaptation depending on cellular context. This study integrated public lung adenocarcinoma transcriptomic cohorts to identify oxidative-stress-associated expression features and evaluate their prognostic relevance. Expression profiles from The Cancer Genome Atlas, Genotype-Tissue Expression project, and GEO series GSE31210, GSE40791, and GSE30219 were analyzed. Differential expression, weighted gene co-expression network analysis, functional enrichment, univariable Cox regression, and least absolute shrinkage and selection operator Cox modeling were combined to derive a risk signature. Immune-cell enrichment, gene set enrichment analysis, gene set variation analysis, and pan-cancer analyses were used for biological characterization. A total of 1,305 genes differed between tumor and control samples, including 498 upregulated and 807 downregulated genes. Intersection of differentially expressed genes, the oxidative-stress-associated co-expression module, and the oxidative-stress gene set yielded 44 genes enriched in responses to reactive oxygen species and hydrogen peroxide, antioxidant and peroxidase activities, focal adhesion, Rap1 signaling, and PI3K-Akt signaling. A seven-gene signature comprising FBLN5, HBB, FYN, HGF, TFAP2A, PLIN5, and F2RL1 stratified the 523-sample training cohort and the 207-sample internal validation cohort into groups with different overall survival. Time-dependent areas under the receiver operating characteristic curve at 1, 3, and 5 years were 0.677, 0.622, and 0.649 in training and 0.613, 0.691, and 0.706 in internal validation. In the 85-case GSE30219 external cohort, corresponding values were 0.588, 0.661, and 0.631; survival separation followed the expected direction but did not reach statistical significance (log-rank P = 0.100). Seventeen immune-cell signatures differed between risk groups, while high-risk tumors were enriched for cell-cycle, DNA-replication, mismatch-repair, glycolytic, E2F, G2M-checkpoint, MYC-target, and mTORC1-related programs. The signature therefore captures reproducible oxidative-stress-associated transcriptional variation with moderate prognostic discrimination. Its clinical utility requires prospective evaluation, complete clinical adjustment, and experimental validation.