Dataset Information

A signature of estimate-stromal-immune score-based genes associated with the prognosis of lung adenocarcinoma.

ABSTRACT:

Background

Immune and stromal component evaluation is necessary to establish accurate prognostic markers for the prediction of clinical outcomes in lung adenocarcinoma (LUAD). We aimed to develop a gene signature based on the Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE)-stromal-immune score in LUAD.

Methods

The transcriptomic profiles of patients with LUAD were obtained from The Cancer Genome Atlas (TCGA), and the immune and stromal scores were derived using the ESTIMATE algorithm. The prognostic signature genes were selected from the differentially expressed genes (DEGs) using the robust partial likelihood-based cox proportional hazards regression method. The negative log-likelihood and the Akaike Information Criterion (AIC) were used to identify the optimal gene signature. The validation was carried out in 2 independent datasets from the Gene Expression Omnibus (GSE68571 and GSE72094).

Results

Patients with high ESTIMATE, stromal, and immune scores had better overall survivals (P=0.0035, P=0.066, and P=0.0077). The expression of thirty-seven genes was related to ESTIMATE-stromal-immune score. A risk stratification model was developed based on a gene signature containing CD74, JCHAIN, and PTGDS. The ESTIMATE-stromal-immune risk score was revealed to be a prognostic factor (P=0.009) after multivariate analysis. Four groups were classified based on this risk stratification model, yielding increasing survival outcomes (log-rank test, P=0.0051). This risk stratification model and other clinicopathological factors were combined to generate a nomogram. The calibration curves showed perfect agreement between the nomogram-predicted outcomes and those actually observed. Similar observations were made in 2 independent cohorts.

Conclusions

The gene signature based on the ESTIMATE-stromal-immune score could predict the prognosis of patients with LUAD.

SUBMITTER: Ma Q

PROVIDER: S-EPMC8044489 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Background: Immune and stromal cells in the tumor microenvironment (TME) significantly contribute to the prognosis of lung adenocarcinoma; however, the TME-related immune prognostic signature is unknown. The aim of this study was to develop a novel immune prognostic model of the TME in lung adenocarcinoma. Methods: First, the immune and stromal scores among lung adenocarcinoma patients were determined using the ESTIMATE algorithm in accordance with The Cancer Genome Atlas (TCGA) database. Differentially expressed immune-related genes (IRGs) between high and low immune/stromal score groups were analyzed, and a univariate Cox regression analysis was performed to identify IRGs significantly correlated with overall survival (OS) among patients with lung adenocarcinoma. Furthermore, a least absolute shrinkage and selection operator (LASSO) regression analysis was performed to generate TME-related immune prognostic signatures. Gene set enrichment analysis was performed to analyze the mechanisms underlying these immune prognostic signatures. Finally, the functions of hub IRGs were further analyzed to delineate the potential prognostic mechanisms in comprehensive TCGA datasets. Results: In total, 702 intersecting differentially expressed IRGs (589 upregulated and 113 downregulated) were screened. Univariate Cox regression analysis revealed that 58 significant differentially expressed IRGs were correlated with patient prognosis in the training cohort, of which three IRGs (CLEC17A, INHA, and XIRP1) were identified through LASSO regression analysis. A robust prognostic model was generated on the basis of this three-IRG signature. Furthermore, functional enrichment analysis of the high-risk-score group was performed primarily on the basis of metabolic pathways, whereas analysis of the low-risk-score group was performed primarily on the basis of immunoregulation and immune cell activation. Finally, hub IRGs CLEC17A, INHA, and XIRP1 were considered novel prognostic biomarkers for lung adenocarcinoma. These hub genes had different mutation frequencies and forms in lung adenocarcinoma and participated in different signaling pathways. More importantly, these hub genes were significantly correlated with the infiltration of CD4+ T cells, CD8+ T cells, macrophages, B cells, and neutrophils. Conclusions: The robust novel TME-related immune prognostic signature effectively predicted the prognosis of patients with lung adenocarcinoma. Further studies are required to further elucidate the regulatory mechanisms of these hub IRGs in the TME and to develop new treatment strategies.

Project description:BackgroundLung adenocarcinoma (LUAD) is the most common form of lung cancer, accounting for 30% of all cases and 40% of all non-small-cell lung cancer cases. Immune-related genes play a significant role in predicting the overall survival and monitoring the status of the cancer immune microenvironment. The present study was aimed at finding an immune-related gene signature for predicting LUAD patient outcomes.MethodsFirst, we chose the TCGA-LUAD project in the TCGA database as the training cohort for model training. For model validating, we found the datasets of GSE72094 and GSE68465 in the GEO database and took them as the candidate cohorts. We obtained 1793 immune-related genes from the ImmPort database and put them into a univariate Cox proportional hazard model to initially look for the genes with potential prognostic ability using the data of the training cohort. These identified genes then entered into a random survival forests-variable hunting algorithm for the best combination of genes for prognosis. In addition, the LASSO Cox regression model tested whether the gene combination can be further shrinkage, thereby constructing a gene signature. The Kaplan-Meier, Cox model, and ROC curve were deployed to examine the gene signature's prognosis in both cohorts. We conducted GSEA analysis to study further the mechanisms and pathways that involved the gene signature. Finally, we performed integrating analyses about the 22 TICs, fully interpreted the relationship between our signature and each TIC, and highlighted some TICs playing vital roles in the signature's prognostic ability.ResultsA nine-gene signature was produced from the data of the training cohort. The Kaplan-Meier estimator, Cox proportional hazard model, and ROC curve confirmed the independence and predictive ability of the signature, using the data from the validation cohort. The GSEA analysis results illustrated the gene signature's mechanism and emphasized the importance of immune-related pathways for the gene signature. 22 TICs immune infiltration analysis revealed resting mast cells' key roles in contributing to gene signature's prognostic ability.ConclusionsThis study discovered a novel immune-related nine-gene signature (BTK, CCR6, S100A10, SEMA3C, GPI, SCG2, TNFRSF11A, CCL20, and DKK1) that predicts LUAD prognosis precisely and associates with resting mast cells strongly.

Project description:Background: Lung adenocarcinoma (LUAD) is an exceedingly diverse disease, making prognostication difficult. Inflammatory responses in the tumor or the tumor microenvironment can alter prognosis in the process of the ongoing cross-talk between the host and the tumor. Nonetheless, Inflammatory response-related genes' prognostic significance in LUAD, on the other hand, has yet to be determined. Materials and Methods: The clinical data as well as the mRNA expression patterns of LUAD patients were obtained from a public dataset for this investigation. In the TCGA group, a multigene prognostic signature was built utilizing LASSO Cox analysis. Validation was executed on LUAD patients from the GEO cohort. The overall survival (OS) of low- and high-risk cohorts was compared utilizing the Kaplan-Meier analysis. The assessment of independent predictors of OS was carried out utilizing multivariate and univariate Cox analyses. The immune-associated pathway activity and immune cell infiltration score were computed utilizing single-sample gene set enrichment analysis. GO keywords and KEGG pathways were explored utilizing gene set enrichment analysis. Results: LASSO Cox regression analysis was employed to create an inflammatory response-related gene signature model. The high-risk cohort patients exhibited a considerably shorter OS as opposed to those in the low-risk cohort. The prognostic gene signature's predictive ability was demonstrated using receiver operating characteristic curve analysis. The risk score was found to be an independent predictor of OS using multivariate Cox analysis. The functional analysis illustrated that the immune status and cancer-related pathways for the two-risk cohorts were clearly different. The tumor stage and kind of immune infiltrate were found to be substantially linked with the risk score. Furthermore, the cancer cells' susceptibility to anti-tumor medication was substantially associated with the prognostic genes expression levels. Conclusion: In LUAD, a new signature made up of 8 inflammatory response-related genes may be utilized to forecast prognosis and influence immunological state. Inhibition of these genes could also be used as a treatment option.

Project description:Background: Lung cancer has a high incidence and mortality rate worldwide. Vasculogenic mimicry (VM) is a specific modality of tumor angiogenesis that could potentially be a new target for tumor therapy. The purpose of this study was to explore the role of VM-related genes in assessing the prognosis and immune landscape of lung cancer. Methods: VM-related genes were obtained from previous studies, and the expression data and clinical data of lung adenocarcinoma (LUAD) patients were obtained from the TCGA database and GEO database. We performed enrichment analysis of 24 VM-related genes and screened hub genes by constructing a protein-protein interaction network and using Cytoscape software. Subsequently, we developed the VM score based on univariate Cox regression analysis and Lasso analysis and validated the VM score on the GSE72094 dataset. In addition, we constructed a nomogram based on the VM score in the TCGA cohort. Finally, we explored the correlation between the VM score and the tumor microenvironment, immune cell infiltration, immune checkpoints, and drug sensitivity. Results: Enrichment analysis revealed that VM-related genes were associated with the HIF signaling pathway and angiogenic pathway. We developed a VM score based on 3 genes (EPHA2, LAMC2 and LOXL2) in LUAD patients. Kaplan-Meier analysis showed that the VM score was associated with poor prognosis in LUAD patients. The receiver operating characteristic curve suggested that the VM score and nomogram are valid predictors for the overall survival of LUAD patients. The VM score was significantly correlated with immune cell infiltration, such as naïve B cells, neutrophils, and eosinophils, and there was a difference in the TME between the high VM score group and the low VM score group. LUAD patients in the high VM score group were more sensitive to antitumor drugs. Conclusion: In summary, the VM score developed in this study is a valuable indicator for evaluating the prognosis and immune landscape of LUAD patients. VM may be a potential target for antitumor therapy in lung cancer.

Project description:The neutrophils exhibit both anti-tumor and pro-tumor effects in cancers. The correlation between neutrophils and tumor development in lung adenocarcinoma (LUAD) is still uncertain, possibly due to a lack of specific neutrophil infiltration evaluation methods. In this study, we identified 30 hub genes that were significantly associated with neutrophil infiltration in LUAD through data mining, survival analysis, and multiple tumor-infiltrating immune cells (TICs) analysis, including TIMER, CIBERSORT, QUANTISEQ, XCELL, and MCPCOUNTER. Consensus clustering analysis showed that these 30 hub genes were correlated with clinical features in LUAD. We further developed a neutrophil scoring system based on these hub genes. The neutrophil score was significantly correlated with prognosis and tumor immune microenvironment (TIME) in LUAD. It was also positively associated with PD-L1 expression and negatively associated with tumor mutational burden (TMB). When combined with the neutrophil score, the predictive capacity of PD-L1 and TMB for prognosis was significantly improved. Thus, the 30 hub genes might play an essential role in the interaction of neutrophils and LUAD, and the neutrophil scoring system might effectually assess the infiltration of neutrophils. Furthermore, we verified the expression of these 30 genes in the LUAD tumor tissues collected from our department. We further found that overexpressed TNFAIP6 and TLR6 and downregulated P2RY13, SCARF1, DPEP2, PRAM1, CYP27A1, CFP, GPX3, and NCF1 in LUAD tissue might be potentially associated with neutrophils pro-tumor effects. The following in vitro experiments demonstrated that TNFAIP6 and TLR6 were significantly overexpressed, and P2RY13 and CYP27A1 were significantly downregulated in LUAD cell lines, compared to BEAS-2B cells. Knocking down TNFAIP6 in A549 and PC9 resulted in the upregulation of FAS, CCL3, and ICAM-1, and the downregulation of CCL2, CXCR4, and VEGF-A in neutrophils when co-culturing with the conditioned medium (CM) from LUAD cells. Knocking down TNFAIP6 in LUAD also led to an elevated early apoptosis rate of neutrophils. Therefore, overexpressed TNFAIP6 in LUAD cancer cells might lead to neutrophils "N2" polarization, which exhibited pro-tumor effects. Further research based on the genes identified in this pilot study might shed light on neutrophils' effects on LUAD in the future.

Project description:BackgroundLung cancer has the highest morbidity and mortality of cancers worldwide. Lung adenocarcinoma (LUAD) is the most common pathological subtype of lung cancer and surgery is its most common treatment. The dysregulated expression of DNA repair genes is found in a variety of cancers and has been shown to affect the origin and progression of these diseases. However, the function of DNA repair genes in surgically-treated LUAD is unclear.MethodsWe sought to determine the association between the signature of DNA repair genes for patients with surgical LUAD and their overall prognosis. We obtained gene expression data and corresponding clinical information of LUAD from The Cancer Genome Atlas (TCGA) database. The differently expressed DNA repair genes of surgically-treated LUAD and normal tissues were identified using the Wilcoxon rank-sum test. We used uni- and multivariate Cox regression analyses to shrink the aberrantly expressed genes, which were then used to construct the prognostic signature and the risk score formula associated with the independent prognosis of surgically-treated LUAD. We used Kaplan-Meier and Cox hazard ratio analyses to confirm the diagnostic and prognostic roles. Two validation sets (GSE31210 and GSE37745) were downloaded from the Gene Expression Omnibus (GEO) and were used to externally verify the prognostic value of the signature. OSluca online database verifies the hazard ratio for the DNA repair genes by which the signature was constructed. We investigated the correlation between the signature of the DNA repair genes and the clinical parameters. The potential molecular mechanisms and pathways of the prognostic signature were explored using Gene Set Enrichment Analysis (GSEA).ResultsWe determined the prognostic signature based on six DNA repair genes (PLK1, FOXM1, PTTG1, CCNO, HIST3H2A, and BLM) and calculated the risk score based on this formula. Patients with surgically-treated LUAD were divided into high-risk and low-risk groups according to the median risk score. The high-risk group showed poorer overall survival than the low-risk group; the signature was used as an independent prognostic indicator and had a greater prognostic value in surgically-treated LUAD. The prognostic value was replicated in GSE31210 and GSE37745. OSluca online database analysis shows that six DNA repair genes were associated with poor prognosis in most lung cancer datasets. The prognostic signature risk score correlated with the pathological stage and smoking status in surgically-treated LUAD. The GSEA of the risk signature in high-risk patients showed pathways associated with the cell cycle, oocyte meiosis, mismatch repair, homologous recombination, and nucleotide excision repair.ConclusionsA six-DNA repair gene signature was determined using TCGA data mining and GEO data verification. The gene signature may serve as a novel prognostic biomarker and therapeutic target for surgically-treated LUAD.