A prognostic nomogram for lung adenocarcinoma based on immune-infiltrating Treg-related genes: from bench to bedside

Accumulating evidence suggests that lymphocyte infiltration in the tumor microenvironment is positively correlated with tumorigenesis and development, while the role of Tregs (regulatory T cells) has been controversial. Therefore, we attempted to discover the possible value of Tregs for lung adenocarcinoma (LUAD).

The gene-sequencing data of LUAD were applied from three Gene Expression Omnibus (GEO) datasets-GSE10072, GSE32863 and GSE43458; the corresponding fractions of tumor-infiltrating immune cells were extracted from the CIBERSORTx portal. Weighted gene coexpression network analysis (WGCNA) and protein-protein interaction (PPI) network analysis were conducted to identify the significant module and candidate genes related to Tregs. The role of candidate genes in LUAD was further verified using data from The Cancer Genome Atlas (TCGA) database. Finally, we constructed a nomogram model to predict the prognosis of LUAD by plotting Kaplan-Meier (K-M), receiver operating characteristic (ROC) and calibration curves, which elucidated the performance of the nomogram.

In total, 10,047 genes in 333 samples (196 tumor and 137 normal samples) from the GEO database were included. By WGCNA and PPI analysis, we identified a significant black module and 36 candidate genes related to Treg. Next, the candidate genes were verified using TCGA data by Cox regression analysis to screen 13 hub genes that stratified LUAD patients into low- or high-risk groups. Low-risk patients showed a significantly longer overall survival (OS) than high-risk patients (3-year OS: 70.2% vs. 35.2%; 5-year OS: 36.6% vs. 0; P=1.651E-09), and the areas under the ROC curves (AUCs) showed good (3-year AUC: 0.733; 5-year AUC: 0.777). Next, we constructed a survival nomogram combining the hub genes and clinical parameters; the low-risk patients still showed a favorable prognosis compared with that of the high-risk patients (P=7.073E-13), and the AUCs were better (3-year AUC: 0.763; 5-year AUC: 0.873).

We revealed the role of immune-infiltrating Treg-related genes in LUAD and constructed a prognostic nomogram, which may help clinicians make optimal therapeutic decisions and help patients obtain better outcomes.