Immune and stromal scoring system associated with tumor microenvironment and prognosis: a gene-based multi-cancer analysis.
Ontology highlight
ABSTRACT: Background
Tumor microenvironment (TME) is associated with tumor progression and prognosis. Previous studies provided tools to estimate immune and stromal cell infiltration in TME. However, there is still a lack of single index to reflect both immune and stromal status associated with prognosis and immunotherapy responses.Methods
A novel immune and stromal scoring system named ISTMEscore was developed. A total of 15 datasets were used to train and validate this system, containing 2965 samples from lung adenocarcinoma, skin cutaneous melanoma and head and neck squamous cell carcinoma.Results
The patients with high immune and low stromal scores (HL) were associated with low ratio of T cell co-inhibitory/stimulatory molecules and low levels of angiogenesis markers, while the patients with low immune and high stromal scores (LH) had the opposite characteristics. The HL patients had immune-centered networks, while the patients with low immune and low stromal scores (LL) had desert-like networks. Moreover, copy number alteration burden was decreased in the HL patients. For the clinical characteristics, our TME classification was an independent prognostic factor. In the 5 cohorts with immunotherapy, the LH patients were linked to the lowest response rate.Conclusions
ISTMEscore system could reflect the TME status and predict the prognosis. Compared to previous TME scores, our ISTMEscore was superior in the prediction of prognosis and immunotherapy response.
SUBMITTER: Zeng Z
PROVIDER: S-EPMC8336334 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Objective: Infiltrating immune and stromal cells are essential for osteosarcoma progression. This study set out to analyze immune-stromal score-based gene signature and molecular subtypes in osteosarcoma. Methods: The immune and stromal scores of osteosarcoma specimens from the TARGET cohort were determined by the ESTIMATE algorithm. Then, immune-stromal score-based differentially expressed genes (DEGs) were screened, followed by univariate Cox regression analysis. A LASSO regression analysis was applied for establishing a prognostic model. The predictive efficacy was verified in the GSE21257 dataset. Associations between the risk scores and chemotherapy drug sensitivity, immune/stromal scores, PD-1/PD-L1 expression, immune cell infiltrations were assessed in the TARGET cohort. NMF clustering analysis was employed for characterizing distinct molecular subtypes based on immune-stromal score-based DEGs. Results: High immune/stromal scores exhibited the prolonged survival duration of osteosarcoma patients. Based on 85 prognosis-related stromal-immune score-based DEGs, a nine-gene signature was established. High-risk scores indicated undesirable prognosis of osteosarcoma patients. The AUCs of overall survival were 0.881 and 0.849 in the TARGET cohort and GSE21257 dataset, confirming the well predictive performance of this signature. High-risk patients were more sensitive to doxorubicin and low-risk patients exhibited higher immune/stromal scores, PD-L1 expression, and immune cell infiltrations. Three molecular subtypes were characterized, with distinct clinical outcomes and tumor immune microenvironment. Conclusion: This study developed a robust prognostic gene signature as a risk stratification tool and characterized three distinct molecular subtypes for osteosarcoma patients based on immune-stromal score-based DEGs, which may assist decision-making concerning individualized therapy and follow-up project.
| S-EPMC8495166 | biostudies-literature
Project description:Hypoxia, a typical hallmark of numerous tumors, indicates poor infiltration of antitumor lymphocytes, as well as facilitates the development, progression, and drug resistance of malignant cells. Here, the present research was performed to identify novel hypoxia-related molecular markers and their correlation to the tumor immune microenvironment (TIME) in colon cancer. The expression of hypoxia-related gene signature was extracted from The Cancer Genome Atlas (TCGA) COAD cohort. Based on this signature, a risk score model was constructed using the Lasso regression model. Its discrimination ability and stability were validated in another independent cohort (GSE17536) from Gene Expression Omnibus (GEO) database. Moreover, molecular biology experiments (quantitative real-time PCR and multiple immunohistochemistry) were performed to validate the results of bioinformatics analyses. Three hub genes, including PPFIA4, SERPINE1, and STC2, were chosen to build the risk score model. All of these genes were increasingly expressed in the hypoxia subgroup (HS). Compared with the normoxia subgroup (NS), HS had worse pathological features (T, N, M, and stage) and overall survival (OS), more expression of immune checkpoint molecules, poorer infiltration of some pro-inflammation immune cells (CD4+ T cells and CD8+ T cells), and enriched infiltration of M0/M2 macrophages. After the risk model was proven to be valuable and stable, a nomogram was built based on this model and some clinicopathological factors. Moreover, it had been identified that three hub genes were all increasingly expressed in hypoxic conditions by quantitative real-time PCR (qPCR). The results of multiple immunohistochemistry (mIHC) also showed that higher expression of hub genes was associated with poorer infiltration of pro-inflammation immune cells (CD8+ T cells and M1 macrophages) and richer infiltration of anti-inflammation immune cells (Treg cells and M2 macrophages). In conclusion, the present study uncovered the relations among hypoxia, TIME, and clinicopathological features of colon cancer. It might provide new insight and a potential therapeutic target for immunotherapy.
| S-EPMC9208084 | biostudies-literature
Project description:Tumor-associated neutrophils (TANs) can promote tumor progression. This study aimed to investigate the molecular signature that predict the prognosis and immune response of breast cancer (BRCA) based on TAN-related gene (TANRG) expression data. The RNA-seq data of BRCA were gathered from The Cancer Genome Atlas (TCGA) and gene expression omnibus (GEO) datasets. Univariate Cox regression analysis and the least absolute shrinkage and selection operator for selecting prognostic genes. A neo-TAN-related risk signature was constructed by multivariate Cox regression analysis. Time-dependent receiver operating characteristic (ROC) curve analyses and Kaplan-Meier analyses were performed to validate the signature in GEO cohorts and the triple-negative breast cancer (TNBC) subtype. We constructed an independent prognostic factor model with 11 TANRGs. The areas under the ROC curve (AUCs) of the TCGA training cohorts for 3-, 5-, and 7-year overall survival were 0.72, 0.73, and 0.73, respectively. The AUCs of the GEO test cohorts for 3-, 5-, and 7-year overall survival were 0.83, 0.89, and 0.94 (GSE25066) and 0.67, 0.69, and 0.73 (GSE58812), respectively. The proportion of immune subtypes differed among the different risk groups. The IC50 values differed significantly between risk groups and can be used as a guide for systemic therapy. The prognostic model developed by TANRGs has excellent predictive performance in BRCA patients. In addition, this feature is closely related to the prediction of survival, immune activity and treatment response in BRCA patients.
| S-EPMC10912776 | biostudies-literature
Project description:BackgroundThe glioma-associated stromal cell (GASC) is a recently identified type of cell in the glioma microenvironment and may be a prognostic marker for glioma. However, the potential mechanisms of GASCs in the glioma microenvironment remain largely unknown. In this work, we aimed to explore the mechanisms of GASCs in gliomas, particularly in high-grade gliomas (HGG).MethodsWe used glioma datasets from The Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA). We utilized the Single-sample Gene Set Enrichment Analysis (ssGSEA) algorithm to discriminate between patients with high or low GASC composition. The xCELL and CIBERSORT algorithms were used to analyze the composition of stromal cells and immune cells. Risk score and a nomogram model were constructed for prognostic prediction of glioma.ResultsWe observed for the first time that the levels of M2 macrophages and immune checkpoints (PD-1, PD-L1, PD-L2, TIM3, Galectin-9, CTLA-4, CD80, CD86, CD155, and CIITA) were significantly higher in the high GASC group and showed positive correlation with the GASC score in all glioma population and the HGG population. Copy number variations of DR3 and CIITA were higher in the high-GASC group. THY1, one of the GASC markers, exhibited lower methylation in the high GASC group. The constructed risk score was an independent predictor of glioma prognostics. Finally, a credible nomogram based on the risk score was established.ConclusionsGASCs stimulate glioma malignancy through the M2 macrophage, and are associated with the level of immune checkpoints in the glioma microenvironment. The methylation of THY1 could be used as prognostic indicator and treatment target for glioma. However, further studies are required to verify these findings.
| S-EPMC8117153 | biostudies-literature
Project description:BackgroundThe tumor microenvironment (TME) involves infiltration of multiple immune cell subsets, which could influence the prognosis and clinical characteristics. The increasing evidence on the role of tumor-infiltrating lymphocytes (TILs) in primary and metastatic melanomas supports that the immune system is involved in the progression and outcomes of melanoma. However, the immune infiltration landscape in melanoma has not been systematically elucidated.MethodsIn this study, we used CIBERSORT and ESTIMATE algorithms to analyze immune infiltration pattern of 993 melanoma samples. Then we screened differential expression genes (DEGs) related to immune subtypes and survival. The immune cell infiltration (ICI) score was constructed by using principal-component analysis (PCA) based on immune signature genes from DGEs. Gene set enrichment analysis (GSEA) was applied to explore high and low ICI score related pathways. Finally, the predictive ability of ICI score was evaluated in survival prognosis and immunotherapy benefit.ResultWe identified three ICI clusters and three gene clusters associated with different immune subtypes and survival outcomes. Then the ICI score was constructed, and we found that high ICI score exhibited activated immune characteristics and better prognosis. High ICI score was significantly enriched in immune pathways and highly expressed immune signature genes. More importantly, we confirmed that melanoma patients with high ICI score had longer overall survival and rate of response to immunotherapy.ConclusionWe presented a comprehensive immune infiltration landscape in melanoma. Our results will facilitate understanding of the melanoma tumor microenvironment and provide a new immune therapy strategy.
| S-EPMC8134682 | biostudies-literature
Project description:BackgroundMelanoma is a cancer of the skin with potential to spread to other organs and is responsible for most deaths due to skin cancer. It is imperative to identify immune biomarkers for early melanoma diagnosis and treatment.Results63 immune-related genes of the total 1039 unique IRGs retrieved were associated with overall survival of melanoma. A multi-IRGs classifier constructed using eight IRGs showed a powerful predictive ability. The classifier had better predictive power compared with the current clinical data. GSEA analysis showed multiple signaling differences between high and low risk score group. Furthermore, biomarker was associated with multiple immune cells and immune infiltration in tumor microenvironment.ConclusionsThe immune-related genes prognosis biomarker is an effective potential prognostic classifier in the immunotherapies and surveillance of melanoma.MethodsMelanoma samples of genes were retrieved from TCGA and GEO databases while the immune-related genes (IRGs) were retrieved from the ImmPort database. WGCNA, Cox regression analysis and LASSO analysis were used to classify melanoma prognosis. ESTIMATE and CIBERSORT algorithms were used to explore the relationship between risk score and tumor immune microenvironment. GSEA analysis was performed to explore the biological signaling pathway.
| S-EPMC7202520 | biostudies-literature
Project description:Malignant pleural mesothelioma (MPM) is a rare but highly aggressive thoracic malignancy. ESTIMATE algorithm-derived immune scores are commonly used to quantify the immune and stromal components in tumors. Thus, this algorithm may help determine the tumor microenvironment (TME)-related gene expression profile associated with tumor immunity. This study aimed at mining public databases to determine a potential correlation between differentially expressed genes (DEGs) and survival in patients with MPM. We categorized patients from the Gene Expression Omnibus database according to their immune/stromal scores into high- and low-score groups. Functional enrichment analysis and the construction of protein-protein interaction networks showed that the DEGs identified were primarily involved in the TME. Furthermore, we validated these genes in an independent cohort of patients with MPM from The Cancer Genome Atlas database. DEG analysis showed that 29 DEGs were cancer driver genes. Subsequently, 14 TME-related genes, which have been previously neglected, were shown to exhibit significant prognostic potential in MPM. In conclusion, immune/stromal scores are novel predictors of a poor prognosis in patients with MPM. We identified DEGs that are involved in immunity against MPM and may contribute to patient survival. Owing to their potential as prognostic factors for MPM, these 14 TME-related genes need to be studied in detail in the future.
| S-EPMC7526499 | biostudies-literature
Project description:The phenotype of pyroptosis has been extensively studied in a variety of tumors, but the relationship between pyroptosis and esophageal squamous cell carcinoma (ESCC) remains unclear. Here, 22 pyroptosis genes were downloaded from the website of Gene Set Enrichment Analysis (GSEA), 79 esophageal squamous cell carcinoma samples and GSE53625 containing 179 pairs of esophageal squamous cell carcinoma samples were collected from the Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO), respectively. Then, pyroptosis subtypes of esophageal squamous cell carcinoma were obtained by cluster analysis according to the expression difference of pyroptosis genes, and a pyroptosis scoring model was constructed by the pyroptosis-related genes screened from different pyroptosis subtypes. Time-dependent receiver operator characteristic (timeROC) curves and the area under the curve (AUC) values were used to evaluate the prognostic predictive accuracy of the pyroptosis scoring model. Kaplan-Meier method with log-rank test were conducted to analyze the impact of the pyroptosis scoring model on overall survival (OS) of patients with esophageal squamous cell carcinoma. Nomogram models and calibration curves were used to further confirm the effect of the pyroptosis scoring model on prognosis. Meanwhile, CIBERSORTx and ESTIMATE algorithm were applied to calculate the influence of the pyroptosis scoring model on esophageal squamous cell carcinoma immune microenvironment. Our findings revealed that the pyroptosis scoring model established by the pyroptosis-related genes was associated with the prognosis and immune microenvironment of esophageal squamous cell carcinoma, which can be used as a biomarker to predict the prognosis and act as a potential target for the treatment of esophageal squamous cell carcinoma.
| S-EPMC9845255 | biostudies-literature
Project description:Background: Pancreatic cancer (PC) is one of the most lethal types of cancer with extremely poor diagnosis and prognosis, and the tumor microenvironment plays a pivotal role during PC progression. Poor prognosis is closely associated with the unsatisfactory results of currently available treatments, which are largely due to the unique pancreatic tumor microenvironment (TME). Methods: In this study, a total of 177 patients with PC from The Cancer Genome Atlas (TCGA) cohort and 65 patients with PC from the GSE62452 cohort in Gene Expression Omnibus (GEO) were included. Based on the proportions of 22 types of infiltrated immune cell subpopulations calculated by cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT), the TME was classified by K-means clustering and differentially expressed genes (DEGs) were determined. A combination of the elbow method and the gap statistic was used to explore the likely number of distinct clusters in the data. The ConsensusClusterPlus package was utilized to identify radiomics clusters, and the samples were divided into two subtypes. Result: Survival analysis showed that the patients with TMEscore-high phenotype had better prognosis. In addition, the TMEscore-high had better inhibitory effect on the immune checkpoint. A total of 10 miRNAs, 311 DEGs, and 68 methylation sites related to survival were obtained, which could be biomarkers to evaluate the prognosis of patients with pancreatic cancer. Conclusions: Therefore, a comprehensive description of TME characteristics of pancreatic cancer can help explain the response of pancreatic cancer to immunotherapy and provide a new strategy for cancer treatment.
| S-EPMC8217872 | biostudies-literature
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.
| S-EPMC10777728 | biostudies-literature