Dataset Information

Construction and Verification of a Hypoxia-Stemness-Based Gene Signature for Risk Stratification in Esophageal Cancer

ABSTRACT:

Background

Numerous studies have shown that esophageal cancer (ESCA) contains areas of intertumoral hypoxia. It is widely accepted that the association of hypoxia with cancer stemness in the tumor microenvironment of ESCA is of profound clinical significance. However, reliable prognostic signatures based on hypoxia and cancer stemness are still lacking in ESCA.

Material/Methods

The t-SNE algorithm was used to estimate the hypoxia status based on the transcriptome profiles of the discovery cohort in the TCGA database. Median values of the stemness index were used to group and identify stemness-associated differentially expressed genes (DEGs). The LASSO method and Cox regression model were combined to screen for prognostic genes and to establish a genetic signature based on hypoxia-stemness. The robustness of the prognostic model was then tested in an external independent validation cohort of the GEO database.

Results

A total of 8 genes – FBLN2, IL17RB, CYP2W1, AMTN, FABP1, FOXA2, GAS1, and CTSF – were identified to construct a gene signature for ESCA risk stratification. Overall survival was significantly lower in the high-risk group than in the low-risk group in both the internal discovery set and the external validation set. The risk score was found to be an independent prognostic factor for ESCA patients. In addition, a higher risk score was significantly associated with the sensitivity of ESCA patients to gefitinib, bexarotene, dasatinib, and imatinib.

Conclusions

The hypoxia-stemness-based genetic signature established for the first time in our study could be a promising tool for ESCA cancer risk stratification.

SUBMITTER: Tang K

PROVIDER: S-EPMC8565098 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:ObjectiveDue to the molecular heterogeneity of gastric cancer, only minor patients respond to immunotherapeutic schemes. This study is aimed at developing an immune-based gene signature for risk stratification and immunotherapeutic efficacy assessment in gastric cancer.MethodsAn immune-based gene signature was developed in gastric cancer by LASSO method in the training set. The predictive performance was validated in the external datasets. KEGG pathways related to risk scores were assessed by GSEA. Based on multivariate Cox regression analysis, a nomogram was established. Sensitivity to chemotherapy drugs was evaluated between high- and low-risk samples. The relationships of risk scores with infiltration levels of immune cells, stromal scores, immune scores, immune cell subgroups, and overall response to anti-PD-L1 therapy were determined.ResultsOur results showed that high risk scores were indicative of undesirable survival outcomes both in the training set (p < 0.0001) and the validation set (p = 0.002). Moreover, this signature could independently predict patients' survival (HR: 2.656 (1.919-3.676) and p < 0.001). Subgroup analysis confirmed the sensitivity of this signature in predicting prognosis (all p < 0.05). Cancer-related pathways were primarily enriched in high-risk samples, such as MAPK and TGF-β pathways (p < 0.05). By incorporating stage and the risk score, we established a nomogram for predicting one-, three-, and five-year survival probability. Patients with high-risk scores were more sensitive to chemotherapy drugs (p < 0.05). There was heterogeneity in immune cells between high- and low-risk samples (p < 0.05). Samples with progressive disease exhibited the highest risk score, and those with complete response had the lowest risk score (p < 0.05).ConclusionThis immune-based gene signature might be representative of a promising prognostic classifier for predicting risk stratification and immunotherapeutic efficacy in gastric cancer, assisting personalized therapy and follow-up plan.

Project description:BackgroundThe specific differentiation potential, unlimited proliferation, and self-renewal capacity of cancer stem cells (CSCs) are closely related to the occurrence, recurrence, and drug resistance of hepatocellular carcinoma (HCC), as well as hypoxia. Therefore, an in-depth analysis of the relationship between HCC stemness, oxygenation status, and the effectiveness of immunotherapy is necessary to improve the poor prognosis of HCC patients.MethodsThe weighted gene co-expression network analysis (WGCNA) was utilized to find hypoxia-related genes, and the stemness index (mRNAsi) was evaluated using the one-class logistic regression (OCLR) technique. Based on stemness-hypoxia-related genes (SHRGs), population subgroup categorization using NMF cluster analysis was carried out. The relationship between SHRGs and survival outcomes was determined using univariate Cox regression. The LASSO-Cox regression strategy was performed to investigate the quality and establish the classifier associated with prognosis. The main effect of risk scores on the tumor microenvironment (TME) and its response to immune checkpoint drugs was also examined. Finally, qRT-PCR was performed to explore the expression and prognostic value of the signature in clinical samples.ResultsAfter identifying tumor stemness- and hypoxia-related genes through a series of bioinformatics analyses, we constructed a prognostic stratification model based on these SHRGs, which can be effectively applied to the prognostic classification of HCC patients and the prediction of immune checkpoint inhibitors (ICIs) efficacy. Independent validation of the model in the ICGC cohort yielded good results. In addition, we also constructed hypoxic cell models in Herp3B and Huh7 cells to verify the expression of genes in the prognostic model and found that C7, CLEC1B, and CXCL6 were not only related to the tumor stemness but also related to hypoxia. Finally, we found that the constructed signature had a good prognostic value in the clinical sample.ConclusionsWe constructed and validated a stemness-hypoxia-related prognostic signature that can be used to predict the efficacy of ICIs therapy. We also verified that C7, CLEC1B, and CXCL6 are indeed associated with stemness and hypoxia through a hypoxic cell model, which may provide new ideas for individualized immunotherapy.

Project description:BackgroundIncreasing evidence showed that the clinical significance of the interaction between hypoxia and immune status in tumor microenvironment. However, reliable biomarkers based on the hypoxia and immune status in triple-negative breast cancer (TNBC) have not been well established. This study aimed to explore a gene signature based on the hypoxia and immune status for predicting prognosis, risk stratification, and individual treatment in TNBC.MethodsHypoxia-related genes (HRGs) and Immune-related genes (IRGs) were identified using the weighted gene co-expression network analysis (WGCNA) method and the single-sample gene set enrichment analysis (ssGSEA Z-score) with the transcriptomic profiles from Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) cohort. Then, prognostic hypoxia and immune based genes were identified in TNBC patients from the METABRIC (N = 221), The Cancer Genome Atlas (TCGA) (N = 142), and GSE58812 (N = 107) using univariate cox regression model. A robust hypoxia-immune based gene signature for prognosis was constructed using the least absolute shrinkage and selection operator (LASSO) method. Based on the cross-cohort prognostic hypoxia-immune related gene signature, a comprehensive index of hypoxia and immune was developed and two risk groups with distinct hypoxia-immune status were identified. The prognosis value, hypoxia and immune status, and therapeutic response in different risk groups were analyzed. Furthermore, a nomogram was constructed to predict the prognosis for individual patients, and an independent cohort from the gene expression omnibus (GEO) database was used for external validation.ResultsSix cross-cohort prognostic hypoxia-immune related genes were identified to establish the comprehensive index of hypoxia and immune. Then, patients were clustered into high- and low-risk groups based on the hypoxia-immune status. Patients in the high-risk group showed poorer prognoses to their low-risk counterparts, and the nomogram we constructed yielded favorable performance to predict survival and risk stratification. Besides, the high-risk group had a higher expression of hypoxia-related genes and correlated with hypoxia status in tumor microenvironment. The high-risk group had lower fractions of activated immune cells, and exhibited lower expression of immune checkpoint markers. Furthermore, the ratio of complete response (CR) was greatly declined, and the ratio of breast cancer related events were significantly elevated in the high-risk group.ConclusionThe hypoxia-immune based gene signature we constructed for predicting prognosis was developed and validated, which may contribute to the optimization of risk stratification for prognosis and personalized treatment in TNBC patients.

Project description:BackgroundInflammation is known to have an intricate relationship with tumorigenesis and tumor progression while it is also closely related to tumor immune microenvironment. Whereas the role of inflammation-related genes (IRGs) in lung squamous carcinoma (LUSC) is barely understood. Herein, we recognized IRGs associated with overall survival (OS), built an IRGs signature for risk stratification and explored the impact of IRGs on immune infiltration landscape of LUSC patients.MethodsThe RNA-sequencing and clinicopathological data of LUSC patients were downloaded from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) database, which were defined as training and validation cohorts. Cox regression and least absolute shrinkage and selection operator analyses were performed to build an IRG signature. CIBERSORT, microenvironment cell populations-counter and tumor immune dysfunction and rejection (TIDE) algorithm were used to perform immune infiltration analysis.ResultsA two-IRG signature consisting of KLF6 and SGMS2 was identified according to the training set, which could categorize patients into two different risk groups with distinct OS. Patients in the low-risk group had more anti-tumor immune cells infiltrated while patient with high-risk had lower TIDE score and higher levels of immune checkpoint molecules expressed. The IRG signature was further identified as an independent prognostic factor of OS. Subsequently, a prognostic nomogram including IRG signature, age, and cancer stage was constructed for predicting individualized OS, whose concordance index values were 0.610 (95% CI: 0.568-0.651) in the training set and 0.652 (95% CI: 0.580-0.724) in validation set. Time-dependent receiver operator characteristic curves revealed that the nomogram had higher prediction accuracy compared with the traditional tumor stage alone.ConclusionThe IRG signature was a predictor for patients with LUSC and might serve as a potential indicator of the efficacy of immunotherapy. The nomogram based on the IRG signature showed a relatively good predictive performance in survival.

Project description:Background and purpose: Radioresistance remains a major reason of radiotherapeutic failure in esophageal squamous cell carcinoma (ESCC). Our study is to screen the immune-related long non-coding RNA (ir-lncRNAs) of radiation-resistant ESCC (rr-ESCC) via Gene Expression Omnibus (GEO) database and to construct a prognostic risk model. Methods: Microarray data (GSE45670) related to radioresistance of ESCC was downloaded from GEO. Based on pathologic responses after chemoradiotherapy, patients were divided into a non-responder (17 samples) and responder group (11 samples), and the difference in expression profiles of ir-lncRNAs were compared therein. Ir-lncRNA pairs were constructed for the differentially expressed lncRNAs as prognostic variables, and the microarray dataset (GSE53625) was downloaded from GEO to verify the effect of ir-lncRNA pairs on the long-term survival of ESCC. After modelling, patients are divided into high- and low-risk groups according to prognostic risk scores, and the outcomes were compared within groups based on the COX proportional hazards model. The different expression of ir-lncRNAs were validated using ECA 109 and ECA 109R cell lines via RT-qPCR. Results: 26 ir-lncRNA genes were screened in the GSE45670 dataset with differential expression, and 180 ir-lncRNA pairs were constructed. After matching with ir-lncRNA pairs constructed by GSE53625, six ir-lncRNA pairs had a significant impact on the prognosis of ESCC from univariate analysis model, of which three ir-lncRNA pairs were significantly associated with prognosis in multivariate COX analysis. These three lncRNA pairs were used as prognostic indicators to construct a prognostic risk model, and the predicted risk scores were calculated. With a median value of 2.371, the patients were divided into two groups. The overall survival (OS) in the high-risk group was significantly worse than that in the low-risk group (p < 0.001). The 1-, 2-, and 3-year prediction performance of this risk-model was 0.666, 0.702, and 0.686, respectively. In the validation setting, three ir-lncRNAs were significantly up-regulated, while two ir-lncRNAs were obviouly down-regulated in the responder group. Conclusion: Ir-lncRNAs may be involved in the biological regulation of radioresistance in patients with ESCC; and the prognostic risk-model, established by three ir-lncRNAs pairs has important clinical value in predicting the prognosis of patients with rr-ESCC.