Development and validation of polyamines metabolism-associated gene signatures to predict prognosis and immunotherapy response in lung adenocarcinoma.

Polyamines metabolism is closely related to tumor development and progression, as well as tumor microenvironment (TME). In this study, we focused on exploring whether polyamines metabolism-associated genes would provide prognosis and immunotherapy response prediction in lung adenocarcinoma (LUAD). The expression profile data of polyamines metabolism-associated genes were acquired from the Cancer Genome Atlas (TCGA) database. Utilizing the least absolute shrinkage and selection operator (LASSO) algorithm, we created a risk score model according to polyamines metabolism-associated gene signatures. Meanwhile, an independent cohort (GSE72094) was employed to validate this model. Through the univariate and multivariate Cox regression analyses, the independent prognostic factors were identified. Subsequently, quantitative real-time polymerase chain reaction (qRT-PCR) was performed to detect their expression in LUAD cells. By consensus clustering analysis, polyamines metabolism-associated subgroups were determined in LUAD patients, with differential gene expression, prognosis, and immune characteristics analyses explored. A total of 59 polyamines metabolism genes were collected for this study, of which 14 genes were identified for the construction of risk score model using LASSO method. High- and low- risk groups of LUAD patients in TCGA cohort were distinguished via this model, and high-risk group presented dismal clinical outcomes. The same prognostic prediction of this model had been also validated in GSE72094 cohort. Meanwhile, three independent prognostic factors (PSMC6, SMOX, SMS) were determined for constructing the nomogram, and they were all upregulated in LUAD cells. In addition, two distinct subgroups (C1 and C2) were identified in LUAD patients. Comparing the two subgroups, 291 differentially expressed genes (DEGs) were acquired, mainly enriching in organelle fission, nuclear division, and cell cycle. Comparing to C1 subgroup, the patients in C2 subgroup had favorable clinical outcomes, increased immune cells infiltration, and effective immunotherapy response. This study identified polyamines metabolism-associated gene signatures for predicting the patients' survival, and they were also linked to immune cells infiltration and immunotherapy response in LUAD patients.

Wang N ，Chai M ，Zhu L ，Liu J ，Yu C ，Huang X ... -  《Frontiers in Immunology》

Development of a copper metabolism-related gene signature in lung adenocarcinoma.

The dysregulation of copper metabolism is closely related to the occurrence and progression of cancer. This study aims to investigate the prognostic value of copper metabolism-related genes (CMRGs) in lung adenocarcinoma (LUAD) and its characterization in the tumor microenvironment (TME). The differentially expressed CMRGs were identified in The Cancer Genome Atlas (TCGA) of LUAD. The least absolute shrinkage and selection operator regression (LASSO) and multivariate Cox regression analysis were used to establish the copper metabolism-related gene signature (CMRGs), which was also validated in Gene Expression Omnibus (GEO) database (GSE72094). The expression of key genes was verified by quantitative real-time PCR (qRT-PCR). Then, the CMRGS was used to develop a nomogram to predict the 1-year, 3-year, and 5-year overall survival (OS). In addition, differences in tumor mutation burden (TMB), biological characteristics and immune cell infiltration between high-risk and low-risk groups were systematically analyzed. Immunophenoscore (IPS) and an anti-PD-L1 immunotherapy cohort (IMvigor210) were used to verify whether CMRGS can predict the response to immunotherapy in LUAD. 34 differentially expressed CMRGs were identified in the TCGA dataset, 11 of which were associated with OS. The CMRGS composed of 3 key genes (LOXL2, SLC31A2 and SOD3) had showed good clinical value and stratification ability in the prognostic assessment of LUAD patients. The results of qRT-PCR confirmed the expression of key CMRGs in LUAD and normal tissues. Then, all LUAD patients were divided into low-risk and high-risk groups based on median risk score. Those in the low-risk group had a significantly longer OS than those in the high-risk group (P<0.0001). The area under curve (AUC) values of the nomogram at 1, 3, and 5 years were 0.734, 0.735, and 0.720, respectively. Calibration curves comparing predicted and actual OS were close to ideal model, indicating a good consistency between prediction and actual observation. Functional enrichment analysis showed that the low-risk group was enriched in a large number of immune pathways. The results of immune infiltration analysis also confirmed that there were a variety of immune cell infiltration in the low-risk group. In addition, multiple immune checkpoints were highly expressed in the low-risk group and may benefit better from immunotherapy. CMRGS is a promising biomarker to assess the prognosis of LUAD patients and may be serve as a guidance on immunotherapy.

Chang W ，Li H ，Zhong L ，Zhu T ，Chang Z ，Ou W ，Wang S ... -  《Frontiers in Immunology》

Identification of immune activation-related gene signature for predicting prognosis and immunotherapy efficacy in lung adenocarcinoma.

Lung adenocarcinoma (LUAD) is a major subtype of non-small cell lung cancer (NSCLC) with a highly heterogeneous tumor microenvironment. Immune checkpoint inhibitors (ICIs) are more effective in tumors with a pre-activated immune status. However, the potential of the immune activation-associated gene (IAG) signature for prognosis prediction and immunotherapy response assessment in LUAD has not been established. Therefore, it is critical to explore such gene signatures. RNA sequencing profiles and corresponding clinical parameters of LUAD were extracted from the TCGA and GEO databases. Unsupervised consistency clustering analysis based on immune activation-related genes was performed on the enrolled samples. Subsequently, prognostic models based on genes associated with prognosis were built using the last absolute shrinkage and selection operator (LASSO) method and univariate Cox regression. The expression levels of four immune activation related gene index (IARGI) related genes were validated in 12 pairs of LUAD tumor and normal tissue samples using qPCR. Using the ESTIMATE, TIMER, and ssGSEA algorithms, immune cell infiltration analysis was carried out for different groups, and the tumor immune dysfunction and rejection (TIDE) score was used to evaluate the effectiveness of immunotherapy. Based on the expression patterns of IAGs, the TCGA LUAD cohort was classified into two clusters, with those in the IAG-high pattern demonstrating significantly better survival outcomes and immune cell infiltration compared to those in the IAG-low pattern. Then, we developed an IARGI model that effectively stratified patients into different risk groups, revealing differences in prognosis, mutation profiles, and immune cell infiltration within the tumor microenvironment between the high and low-risk groups. Notably, significant disparities in TIDE score between the two groups suggest that the low-risk group may exhibit better responses to ICIs therapy. The IARGI risk model was validated across multiple datasets and demonstrated exceptional performance in predicting overall survival in LUAD, and an IARGI-integrated nomogram was established as a quantitative tool for clinical practice. The IARGI can serve as valuable biomarkers for evaluating the tumor microenvironment and predicting the prognosis of LUAD patients. Furthermore, these genes probably provide valuable guidance for establishing effective immunotherapy regimens for LUAD patients.

Zeng W ，Wang J ，Yang J ，Chen Z ，Cui Y ，Li Q ，Luo G ，Ding H ，Ju S ，Li B ，Chen J ，Xie Y ，Tong X ，Liu M ，Zhao J ... -  《Frontiers in Immunology》

Prognosis and immunotherapy significances of a cancer-associated fibroblasts-related gene signature in lung adenocarcinoma.

Luo Y ，Zhang S ，Xie H ，Su Q ，He S ，Lei Z ... -  《-》

Identification of novel gene signature for lung adenocarcinoma by machine learning to predict immunotherapy and prognosis.

Lung adenocarcinoma (LUAD) as a frequent type of lung cancer has a 5-year overall survival rate of lower than 20% among patients with advanced lung cancer. This study aims to construct a risk model to guide immunotherapy in LUAD patients effectively. LUAD Bulk RNA-seq data for the construction of a model, single-cell RNA sequencing (scRNA-seq) data (GSE203360) for cell cluster analysis, and microarray data (GSE31210) for validation were collected from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database. We used the Seurat R package to filter and process scRNA-seq data. Sample clustering was performed in the ConsensusClusterPlus R package. Differentially expressed genes (DEGs) between two groups were mined by the Limma R package. MCP-counter, CIBERSORT, ssGSEA, and ESTIMATE were employed to evaluate immune characteristics. Stepwise multivariate analysis, Univariate Cox analysis, and Lasso regression analysis were conducted to identify key prognostic genes and were used to construct the risk model. Key prognostic gene expressions were explored by RT-qPCR and Western blot assay. A total of 27 immune cell marker genes associated with prognosis were identified for subtyping LUAD samples into clusters C3, C2, and C1. C1 had the longest overall survival and highest immune infiltration among them, followed by C2 and C3. Oncogenic pathways such as VEGF, EFGR, and MAPK were more activated in C3 compared to the other two clusters. Based on the DEGs among clusters, we confirmed seven key prognostic genes including CPA3, S100P, PTTG1, LOXL2, MELTF, PKP2, and TMPRSS11E. Two risk groups defined by the seven-gene risk model presented distinct responses to immunotherapy and chemotherapy, immune infiltration, and prognosis. The mRNA and protein level of CPA3 was decreased, while the remaining six gene levels were increased in clinical tumor tissues. Immune cell markers are effective in clustering LUAD samples into different subtypes, and they play important roles in regulating the immune microenvironment and cancer development. In addition, the seven-gene risk model may serve as a guide for assisting in personalized treatment in LUAD patients.

Shu J ，Jiang J ，Zhao G 《Frontiers in Immunology》