An In Silico Analysis Reveals an EMT-Associated Gene Signature for Predicting Recurrence of Early-Stage Lung Adenocarcinoma.

The potential micrometastasis tends to cause recurrence of lung adenocarcinoma (LUAD) after surgical resection and consequently leads to an increase in the mortality risk. Compelling evidence has suggested the underlying mechanisms of tumor metastasis could involve the activation of an epithelial-mesenchymal transition (EMT) program. Hence, the objective of this study was to develop an EMT-associated gene signature for predicting the recurrence of early-stage LUAD. The mRNA expression data of patients with early-stage LUAD were downloaded from Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) available databases. Gene Set Variation Analysis (GSVA) was first performed to provide an assessment of EMT phenotype, whereas Weighted Gene Co-expression Network Analysis (WGCNA) was constructed to determine EMT-associated key modules and genes. Based on the genes, a novel EMT-associated signature for predicting the recurrence of early-stage LUAD was identified using a least absolute shrinkage and selection operator (LASSO) algorithm and a stepwise Cox proportional hazards regression model. Kaplan-Meier survival analysis, receiver operating characteristic (ROC) curves and Cox regression analyses were used to estimate the performance of the identified gene signature. GSVA revealed diverse EMT states in the early-stage LUAD. Further correlation analyses showed that the EMT states presented high correlations with several hallmarks of cancers, tumor purity, tumor microenvironment cells, and immune checkpoint genes. More importantly, Kaplan-Meier survival analyses indicated that patients with high EMT scores had worse recurrence-free survival (RFS) and overall survival (OS) than those with low EMT scores. A novel 5-gene signature (AGL, ECM1, ENPP1, SNX7, and TSPAN12) was established based on the EMT-associated genes from WGCNA and this signature successfully predicted that the high-risk patients had a higher recurrence rate compared with the low-risk patients. In further analyses, the signature represented robust prognostic values in 2 independent validation cohorts (GEO and TCGA datasets) and a combined GEO cohort as evaluated by Kaplan-Meier survival (P-value < .0001) and ROC analysis (AUC = 0.781). Moreover, the signature was corroborated to be independent of clinical factors by univariate and multivariate Cox regression analyses. Interestingly, the combination of the signature-based recurrence risk and tumor-node-metastasis (TNM) stage showed a superior predictive ability on the recurrence of patients with early-stage LUAD. Our study suggests that patients with early-stage LUAD exhibit diverse EMT states that play a vital role in tumor recurrence. The novel and promising EMT-associated 5-gene signature identified and validated in this study may be applied to predict the recurrence of early-stage LUAD, facilitating risk stratification, recurrence monitoring, and individualized management for the patients after surgical resection.

Han Y ，Wong FC ，Wang D ，Kahlert C ... -  《-》

Clinical Significance and Immunometabolism Landscapes of a Novel Recurrence-Associated Lipid Metabolism Signature In Early-Stage Lung Adenocarcinoma: A Comprehensive Analysis.

The early-stage lung adenocarcinoma (LUAD) incidence has increased with heightened public awareness and lung cancer screening implementation. Lipid metabolism abnormalities are associated with lung cancer initiation and progression. However, the comprehensive features and clinical significance of the immunometabolism landscape and lipid metabolism-related genes (LMRGs) in cancer recurrence for early-stage LUAD remain obscure. LMRGs were extracted from Gene Set Enrichment Analysis (GSEA) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Samples from The Cancer Genome Atlas (TCGA) were used as training cohort, and samples from four Gene Expression Omnibus (GEO) datasets were used as validation cohorts. The LUAD recurrence-associated LMRG molecular pattern and signature was constructed through unsupervised consensus clustering, time-dependent receiver operating characteristic (ROC), and least absolute shrinkage and selection operator (LASSO) analyses. Kaplan-Meier, ROC, and multivariate Cox regression analyses and prognostic meta-analysis were used to test the suitability and stability of the signature. We used Gene Ontology (GO), KEGG pathway, immune cell infiltration, chemotherapy response analyses, gene set variation analysis (GSVA), and GSEA to explore molecular mechanisms and immune landscapes related to the signature and the potential of the signature to predict immunotherapy or chemotherapy response. First, two LMRG molecular patterns were established, which showed diverse prognoses and immune infiltration statuses. Then, a 12-gene signature was identified, and a risk model was built. The signature remained an independent prognostic parameter in multivariate Cox regression and prognostic meta-analysis. In addition, this signature stratified patients into high- and low-risk groups with significantly different recurrence rates and was well validated in different clinical subgroups and several independent validation cohorts. The results of GO and KEGG analyses and GSEA showed that there were differences in multiple lipid metabolism, immune response, and drug metabolism pathways between the high- and low-risk groups. Further analyses revealed that the signature-based risk model was related to distinct immune cell proportions, immune checkpoint parameters, and immunotherapy and chemotherapy response, consistent with the GO, KEGG, and GSEA results. This is the first lipid metabolism-based signature for predicting recurrence, and it could provide vital guidance to achieve optimized antitumor for immunotherapy or chemotherapy for early-stage LUAD.

Zhu M ，Zeng Q ，Fan T ，Lei Y ，Wang F ，Zheng S ，Wang X ，Zeng H ，Tan F ，Sun N ，Xue Q ，He J ... -  《Frontiers in Immunology》

A novel epithelial-mesenchymal transition (EMT)-related gene signature of predictive value for the survival outcomes in lung adenocarcinoma.

Lung adenocarcinoma (LUAD) is a remarkably heterogeneous and aggressive disease with dismal prognosis of patients. The identification of promising prognostic biomarkers might enable effective diagnosis and treatment of LUAD. Aberrant activation of epithelial-mesenchymal transition (EMT) is required for LUAD initiation, progression and metastasis. With the purpose of identifying a robust EMT-related gene signature (E-signature) to monitor the survival outcomes of LUAD patients. In The Cancer Genome Atlas (TCGA) database, least absolute shrinkage and selection operator (LASSO) analysis and cox regression analysis were conducted to acquire prognostic and EMT-related genes. A 4 EMT-related and prognostic gene signature, comprising dickkopf-like protein 1 (DKK1), lysyl oxidase-like 2 (LOXL2), matrix Gla protein (MGP) and slit guidance ligand 3 (SLIT3), was identified. By the usage of datum derived from TCGA database and Western blotting analysis, compared with adjacent tissue samples, DKK1 and LOXL2 protein expression in LUAD tissue samples were significantly higher, whereas the trend of MGP and SLIT3 expression were opposite. Concurrent with upregulation of epithelial markers and downregulation of mesenchymal markers, knockdown of DKK1 and LOXL2 impeded the migration and invasion of LUAD cells. Simultaneously, MGP and SLIT3 silencing promoted metastasis and induce EMT of LUAD cells. In the TCGA-LUAD set, receiver operating characteristic (ROC) analysis indicated that our risk model based on the identified E-signature was superior to those reported in literatures. Additionally, the E-signature carried robust prognostic significance. The validity of prediction in the E-signature was validated by the three independent datasets obtained from Gene Expression Omnibus (GEO) database. The probabilistic nomogram including the E-signature, pathological T stage and N stage was constructed and the nomogram demonstrated satisfactory discrimination and calibration. In LUAD patients, the E-signature risk score was associated with T stage, N stage, M stage and TNM stage. GSEA (gene set enrichment analysis) analysis indicated that the E-signature might be linked to the pathways including GLYCOLYSIS, MYC TARGETS, DNA REPAIR and so on. In conclusion, our study explored an innovative EMT based prognostic signature that might serve as a potential target for personalized and precision medicine.

Cui Y ，Wang X ，Zhang L ，Liu W ，Ning J ，Gu R ，Cui Y ，Cai L ，Xing Y ... -  《Frontiers in Oncology》

Development and validation of a robust immune-related prognostic signature in early-stage lung adenocarcinoma.

Wu P ，Zheng Y ，Wang Y ，Wang Y ，Liang N ... -  《Journal of Translational Medicine》

Identification of a three-gene expression signature and construction of a prognostic nomogram predicting overall survival in lung adenocarcinoma based on TCGA and GEO databases.

Lung adenocarcinoma (LUAD) is the major cause of cancer mortality. Traditional prognostic factors have limited importance after including other parameters. Thus, developing a more credible prognostic model combined with genes and clinical parameters is necessary. The messenger RNA (mRNA) expression and clinical information from The Cancer Genome Atlas (TCGA)-LUAD datasets and microarray data from three Gene Expression Omnibus (GEO) databases were obtained. We identified differentially-expressed genes (DEGs) between lung tumor and normal tissues through integrated analysis of the three GEO datasets. Univariate and multivariate Cox regression analyses were conducted to select survival-associated DEGs and to establish a prognostic gene signature which was associated with overall survival (OS). The expression of gene proteins was assessed in 180 LUAD tissue microarrays (TMAs) by immunohistochemistry (IHC). We verified its predictive performance with a Kaplan-Meier (KM) curve, receiver operating characteristic (ROC) curve, and Harrell's concordance index (C-index) and validated it in external GEO databases. Multivariate Cox regression analysis was performed to identify the significant prognostic indicators in LUAD. Furthermore, we established a prognostic nomogram based on TCGA-LUAD dataset. A three-gene signature was constructed to predict the OS of LUAD patients. The KM analysis, ROC curve, and C-index present a good predictive ability of the gene signature in TCGA dataset [P<0.0001; C-index 0.6375; 95% confidence interval (CI): 0.5632-0.7118; area under the ROC curve (AUC) 0.674] and the external GEO datasets (P=0.05, 0.004, and 0.04, respectively). Univariate and multivariate Cox regression analyses also verified that LUAD patients with low-risk scores had a decreased risk of death compared to those with a high-risk score in TCGA database [hazard ratio (HR) =0.3898; 95% CI: 0.1938-0.7842; P<0.05]. Finally, we constructed a nomogram integrating the gene signature and clinicopathological parameters (P<0.0001; C-index 0.762; 95% CI: 0.714-0.845; AUC 0.8136). Compared with conventional staging, a nomogram can effectively improve prognosis prediction. The nomogram is closely associated to the OS of LUAD patients. This consequence may be beneficial to individualized treatment and clinical decision-making.

Zhou Y ，Gao S ，Yang R ，Du C ，Wang Y ，Wu Y ... -  《-》