A Prognostic Risk Signature of Two Autophagy-Related Genes for Predicting Triple-Negative Breast Cancer Outcomes.

Background: Triple-negative breast cancer (TNBC) is recognized as the most aggressive molecular subtype of breast cancer. Recent studies have highlighted the complex role of autophagy in the pathogenesis of TNBC. Methods: In this study, we evaluated 18,330 genes, including 1111 autophagy-related genes, (ARGs), across 579 TNBC samples from online databases. Differentially expressed ARGs in TNBC were identified using high-throughput RNA-seq data from the Cancer Genome Atlas (TCGA). Prognostic factors were examined through Cox regression and multivariate Cox analyses, with predictive efficacy assessed using receiver operating characteristic (ROC) curves. A nomogram integrating the risk signature with clinicopathological factors, such as TNM stage, was developed. Immunohistochemical analysis of clinical samples was also conducted. Results: EIF4EBP1 and NPAS3 were significantly correlated with prognostic outcomes in patients with TNBC. Multivariate Cox regression analysis demonstrated that the expression levels of these two genes were accurate predictors of disease progression in TNBC samples from TCGA and the GSE31519 dataset. The efficacy of this predictive model was validated using ROC curve analysis and calibration plots, confirming its ability to accurately estimate the 1-, 2-, and 3-year survival rates for individuals with TNBC. Additionally, EIF4EBP1 and NPAS3 expression influenced drug sensitivity in TNBC cell lines, with notably lower NPAS3 expression in TNBC tissues, particularly in Stage III cases. This study is the first to report NPAS3 expression in patients with TNBC. Conclusion: The autophagy-related genes EIF4EBP1 and NPAS3 may serve as independent prognostic factors for individuals with TNBC.

Identification and Verification of a Prognostic Risk Signature in Oral Squamous Cell Carcinoma.

INTRODUCTION: Oral squamous cell carcinoma (OSCC) is a prevalent malignant condition. This study aimed to investigate the role of mTORC1 signaling and develop a prognostic model for OSCC. MATERIALS AND METHODS: The single-sample gene set enrichment analysis (ssGSEA) algorithm was utilized to calculate the Z-Score of Hallmarks in OSCC, followed by univariate Cox regression analysis to identify processes associated with prognosis. Weighted gene co-expression network analysis (WGCNA) was performed using transcriptomic data from the cancer genome atlas (TCGA) cohort to identify genes correlated with mTORC1 signaling. A six-gene prognostic model was constructed using multifactorial Cox regression analysis and validated using an external dataset. RESULTS: The study uncovered a strong linkage between mTORC1, glycolysis, hypoxia, and the prognosis of OSCC. mTORC1 signaling emerged as the most significant risk factor, negatively impacting patient survival. Additionally, a six-gene prognostic risk score model was developed which provided a quantitative measure of patients' survival probabilities. Interestingly, within the context of these findings, TP53 gene mutations were predominantly observed in the high-risk group, potentially underlining the genetic complexity of this patient subgroup. Additionally, differential immune cell infiltration and an integrated nomogram were also reported. CONCLUSION: This study highlights the importance of mTORC1 signaling in OSCC prognosis and presents a robust prognostic model for predicting patient outcomes.

Construction and validation of a SASP-related prognostic signature in patients with acute myeloid leukaemia.

Acute myeloid leukaemia (AML) is a common and highly aggressive haematological malignancy in adults. Senescence-associated secretory phenotype (SASP) plays important roles in tumorigenesis and progression of tumour. However, the prognostic value of SASP in patients with AML has not been clarified. The present study aims to explore the prognostic value of SASP and develop a prognostic risk signature for AML. The RNA-sequencing data was collected from the TCGA, GTEx and TARGET databases. Subsequently, differentially expressed gene analysis, univariate Cox regression and LASSO regression were applied to identified prognostic SASP-related genes and construct a prognostic risk-scoring model. The risk score of each patient were calculated and patients were divided into high- or low-risk groups by the median risk score. This novel prognostic signature included 11 genes: G6PD, CDK4, RPS6KA1, UBC, H2BC12, KIR2DL4, HSF1, IFIT3, PIM1, RUNX3 and TRIM21. The patients with AML in the high-risk group had shorter OS, demonstrating that the risk score acted as a prognostic predictor, which was validated in the TAGET-AML dataset. Univariate and multivariate analysis revealed the risk score was an independent prognostic factor in patients with AML. Furthermore, the present study revealed that the risk score was associated with immune landscape, immune checkpoint gene expression and chemotherapeutic efficacy. In the present study, we constructed and validated a unique SASP-related prognostic model to assess therapeutic effect and prognosis in patients with AML, which might contribute to understanding the role of SASP in AML and guiding the treatment for AML.

Immune infiltration and prognosis in gastric cancer: role of NAD+ metabolism-related markers.

Background: This study endeavored to develop a nicotinamide adenine dinucleotide (NAD+) metabolism-related biomarkers in gastric cancer (GC), which could provide a theoretical foundation for prognosis and therapy of GC patients. Methods: In this study, differentially expressed genes (DEGs1) between GC and paraneoplastic tissues were overlapped with NAD+ metabolism-related genes (NMRGs) to identify differentially expressed NMRGs (DE-NMRGs). Then, GC patients were divided into high and low score groups by gene set variation analysis (GSVA) algorithm for differential expression analysis to obtain DEGs2, which was overlapped with DEGs1 for identification of intersection genes. These genes were further analyzed using univariate Cox and least absolute shrinkage and selection operator (LASSO) regression analyses to obtain prognostic genes for constructing a risk model. Enrichment and immune infiltration analyses further investigated investigate the different risk groups, and qRT-PCR validated the prognostic genes. Results: Initially, we identified DE-NMRGs involved in NAD biosynthesis, with seven (DNAJB13, CST2, THPO, CIDEA, ONECUT1, UPK1B and SNCG) showing prognostic significance in GC. Subsequent, a prognostic model was constructed in which the risk score, derived from the expression profiles of these genes, along with gender, emerged as robust independent predictors of patient outcomes in GC. Enrichment analysis linked high-risk patients to synaptic membrane pathways and low-risk to the CMG complex pathway. Tumor immune infiltration analysis revealed correlations between risk scores and immune cell abundance, suggesting a relationship between NAD+ metabolism and immune response in GC. The prognostic significance of our identified genes was validated by qRT-PCR, which confirmed their upregulated expression in GC tissue samples. Conclusion: In this study, seven NAD+ metabolism-related markers were established, which is of great significance for the development of prognostic molecular biomarkers and clinical prognosis prediction for gastric cancer patients.

A novel amino acid metabolism-related gene signature to predict the overall survival of esophageal squamous cell carcinoma patients.

Background: Esophageal squamous cell carcinoma (ESCC) has a poor early detection rate, prognosis, and survival rate. Effective prognostic markers are urgently needed to assist in the prediction of ESCC treatment outcomes. There is accumulating evidence of a strong relationship between cancer cell growth and amino acid metabolism. This study aims to determine the relationship between amino acid metabolism and ESCC prognosis. Methods: This study comprehensively evaluates the association between amino acid metabolism-related gene (AAMRG) expression profiles and the prognosis of ESCC patients based on data from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) databases. Real-time quantitative polymerase chain reaction (RT-qPCR) was used to verify the expression of prognosis-related genes. Results: A univariate Cox regression analysis of TCGA data identified 18 prognosis-related AAMRGs. The gene expression profiles of 90 ESCC tumor and normal tissues were obtained from the GSE20347 and GSE67269 datasets. Two differently expressed genes (DEGs) were considered as ESCC prognosis-related genes; and they were branched-chain amino acid transaminase 1 (BCAT1) and methylmalonic aciduria and homocystinuria type C protein (MMACHC). These two AAMRGs were used to develop a novel AAMRG-related gene signature to predict 1- and 2-year prognostic risk in ESCC patients. Both BCAT1 and MMACHC expression were verified by RT-qPCR. A prognostic nomogram that incorporated clinical factors and BCAT1 and MMACHC gene expression was constructed, and the calibration plots showed that it had good prognostic performance. Conclusions: The AAMRG signature established in our study is efficient and could be used in clinical settings to predict the early prognosis of ESCC patients.

Construction of a 5-Gene super-enhancer-related signature for osteosarcoma prognosis and the regulatory role of TNFRSF11B in osteosarcoma.

Osteosarcoma, one of the most common primary malignancies in children and adolescents, has the primary characteristics of a poor prognosis and high rate of metastasis. This study used super-enhancer-related genes derived from two different cell lines to construct five novel super-enhancer-related gene prognostic models for patients with osteosarcoma. The training and testing datasets were used to confirm the prognostic models of the five super-enhancer-related genes, which resulted in an impartial predictive element for osteosarcoma. The immunotherapy and prediction of the response to anticancer drugs have shown that the risk signature of the five super-enhancer-related genes positively correlate with chemosensitivity. Furthermore, functional analysis of the risk signature genes revealed a significant relationship between gene groups and the malignant characteristics of tumours. TNF Receptor Superfamily Member 11b (TNFRSF11B) was selected for functional verification. Silencing of TNFRSF11B suppressed the proliferation, migration, and invasion of osteosarcoma cells in vitro and suppressed osteosarcoma growth in vivo. Moreover, transcriptome sequencing was performed on MG-63 cells to study the regulatory mechanism of TNFRSF11B in osteosarcoma cells, and it was discovered that TNFRSF11B is involved in the development of osteosarcoma via the phosphoinositide 3-kinase signalling pathway. Following the identification of TNFRSF11B as a key gene, we selected an inhibitor that specifically targeted this gene and performed molecular docking simulations. In addition, risedronic acid inhibited osteosarcoma growth at both cellular and molecular levels. In conclusion, the super-enhancer-related gene signature is a viable therapeutic tool for osteosarcoma prognosis and treatment.

Development of a propionate metabolism-related gene-based molecular subtypes and scoring system for predicting prognosis in bladder cancer.

PURPOSE: Bladder cancer (BLCA) is a prevalent malignancy. Dysregulated propionate metabolism, a key cancer factor, suggests a potential target for treating metastatic cancer. However, a complete understanding of the link between propionate metabolism-related genes (PMRGs) and bladder cancer is lacking. METHODS: From the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases, we gathered BLCA patient data, which was classified into distinct subgroups using non-negative matrix factorization (NMF). Survival and pathway analyses were conducted between these clusters. The PMRGs model, created through univariate Cox and least absolute shrinkage and selection operator (LASSO) analyses, was assessed for prognostic significance using Kaplan-Meier and receiver operating characteristic (ROC) curves. A comprehensive evaluation included clinical, tumor microenvironment (TME), drug sensitivity, and immunotherapy analyses. Finally, the expression of HSD17B1 essential genes was confirmed via quantitative real-time polymerase chain reaction (qRT-PCR), with further validation through Transwell, wound healing, colony-formation, and EDU assays. RESULTS: We discovered two distinct subcategories (CA and CB) within BLCA using NMF analysis, with CA demonstrating significantly better overall survival compared to CB. Additionally, six PMRGs emerged as critical factors associated with propionate metabolism and prognosis. Kaplan-Meier analysis revealed that high-risk PMRGs were correlated with a poorer prognosis in BLCA patients. Moreover, significant differences were observed between the two groups in terms of infiltrated immune cells, immune checkpoint expression, TME scores, and drug sensitivity. Notably, we found that suppressing HSD17B1 gene expression inhibited the invasion of bladder cancer cells. CONCLUSION: Our study proposes molecular subtypes and a PMRG-based score as promising prognostic indicators in BLCA. Additionally, cellular experiments underscore the pivotal role of HSD17B1 in bladder cancer metastasis and invasion, suggesting its potential as a novel therapeutic target.

Characterization of tumor endothelial cells (TEC) in gastric cancer and development of a TEC-based risk signature using single-cell RNA-seq and bulk RNA-seq data.

BACKGROUND: Tumor endothelial cells (TECs) are essential participants in tumorigenesis. This study is focused on elucidating the TEC traits in gastric cancer (GC) and constructing a prognostic risk model to predict the clinical outcome of GC patients. METHODS: Single-cell RNA sequencing (scRNA-seq) data were obtained from the GEO database. Using specific markers, the Seurat R package aided in processing scRNA-seq data and identifying TEC clusters. Based on TEC cluster-associated genes identified by Pearson correlation analysis, TEC-related prognostic genes were screened by lasso-Cox regression analysis, thereby constructing a risk signature. A nomogram was created by combining the risk signature with clinicopathological features. RESULTS: Based on the scRNA-seq data, 5 TEC clusters were discovered in GC, with 3 of them showing prognostic associations in GC. A total of 163 genes were pinpointed among 3302 DEGs as significantly linked to TEC clusters, leading to the formulation of a risk signature comprising 8 genes. Furthermore, there was a notable correlation between the risk signature and the immune cell infiltration. Multivariate analysis findings indicated that the risk signature served as an independent prognostic factor for GC. Moreover, its efficacy in forecasting immune response was validated. CONCLUSION: TEC-based risk model is highly effective in predicting the survival outcomes of GC patients and can forecast the immune response. Targeting TECs may significantly inhibit tumor progression and enhance the efficacy of immunotherapy.

Identification and validation of a prognostic signature based on six immune-related genes for colorectal cancer.

BACKGROUND: Colorectal cancer (CRC) is a prevalent malignancy with high mortality and morbidity rates. Although the significant efficacy of immunotherapy is well established, it is only beneficial for a limited number of individuals with CRC. METHODS: Differentially expressed immune-related genes (DE-IRGs) were retrieved from The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO), and ImmPort databases. A prognostic signature comprising DE-IRGs was developed using univariate, LASSO, and multivariate Cox regression analyses. A nomogram integrating the independent prognostic factors was also developed. CIBERSORT was used to assess immune cell infiltration (ICI). Furthermore, wound-healing, colony formation, migration, and invasion assays were performed to study the involvement of ACTG1 in CRC. RESULTS: A signature including six DE-IRGs was developed. The overall survival (OS) rate was accurately estimated for TCGA and GSE38832 cohorts. The risk score (RS) of the signature was an independent factor for OS. Moreover, a nomogram encompassing age, RS, and pathological T stage accurately predicted the long-term OS probability of individuals with CRC. The high-risk group had an elevated proportion of patients treated with ICIs, including native B cells, relative to the low-risk group. Additionally, ACTG1 expression was upregulated, which supported the proliferation, migration, and invasion abilities of CRC cells. CONCLUSIONS: An immune-related prognostic signature was developed for predicting OS and for determining the immune status of individuals with CRC. The present study provides new insights into accurate immunotherapy for individuals with CRC. Moreover, ACTG1 may serve as a new immune biomarker.

Establishment of disulfidptosis-related LncRNA signature as biomarkers in colon adenocarcinoma.

PURPOSE: Metabolic reprogramming is a hallmark of cancer and plays a key role in precision oncology treatment. Long non-coding RNAs (lncRNAs) regulate cancer cell behavior, including metabolism. Disulfidptosis, a newly identified form of regulated cell death triggered by glucose starvation, has yet to be fully understood in colon adenocarcinoma (COAD). This study aimed to confirm the existence and role of disulfidptosis in COAD and identify disulfidptosis-related lncRNAs that may be targeted to induce disulfidptosis in COAD. METHODS: PI and F-actin staining were used to observe disulfidptosis in COAD cell lines. Disulfidptosis-related lncRNAs were identified based on the expression of disulfidptosis-associated genes in the TCGA-COAD database. A four-lncRNA signature for disulfidptosis was established. Subsequently, loss-of-function assays explored the roles of AC013652.1 and MCM3AP-AS1 in disulfidptosis. RESULTS: Disulfidptosis was observed in COAD cells under glucose starvation and could be reversed by agents that prevent disulfide stress, such as dithiothreitol (DTT) and tris-(2-carboxyethyl)-phosphine (TCEP). The prognostic value of disulfidptosis-associated genes in COAD patients was confirmed, with higher expression indicating longer survival. A disulfidptosis-related lncRNA signature comprising four lncRNAs was established based on the expression of these genes. Among these, AC013652.1 and MCM3AP-AS1 predicted worse prognoses. Furthermore, inhibiting AC013652.1 or MCM3AP-AS1 increased disulfidptosis-associated gene expression and cellular death, which could be reversed by DTT and TCEP. CONCLUSIONS: This study provides hitherto undocumented evidence of the existence of disulfidptosis and the prognostic value of disulfidptosis-associated genes in COAD. Importantly, we identified lncRNAs AC013652.1 and MCM3AP-AS1, which suppress disulfidptosis and may serve as potential therapeutic targets for COAD.

Integrated analysis of tumor-associated macrophages and M2 macrophages in CRC: unraveling molecular heterogeneity and developing a novel risk signature.

BACKGROUND: Emerging investigations have increasingly highlighted the critical role of tumor-associated macrophages (TAMs) and M2 macrophages in cancer development, progression, and metastasis, marking them as potential targets in various cancer types. The main objective of this research is to discover new biomarkers associated with TAM-M2 macrophages in colorectal cancer (CRC) and to dissect the molecular heterogeneity of CRC by combining single-cell RNA sequencing and bulk RNA-seq data. METHODS: By utilizing weighted gene co-expression network analysis (WGCNA), we acquired TAM-M2-associated genes by intersecting TAM marker genes obtained from scRNA-seq data with module genes of M2 macrophages derived from bulk RNA-seq data. We employed least absolute shrinkage and selection operator (LASSO) Cox analysis to select predictive biomarkers from these TAM-M2-related genes. Quantitative polymerase chain reaction (qPCR) was employed to validate the mRNA expression levels of the genes identified in the screening. This led to the development of the TAM-M2-related signature (TAMM2RS). We also conducted functional and immune landscape analyses of different risk groups. RESULTS: The combination of scRNA-seq and bulk RNA-seq analyses yielded 377 TAM-M2-related genes. DAPK1, NAGK, and TRAF1 emerged as key prognostic genes in CRC, which were identified through LASSO Cox analysis. Utilizing these genes, we constructed and validated the TAMM2RS, demonstrating its effectiveness in predicting survival in CRC patients. CONCLUSION: Our research offers a thorough investigation into the molecular mechanisms associated with TAM-M2 macrophages in CRC and unveils potential therapeutic targets, offering new insights for treatment strategies in colorectal cancer.

Glutathione metabolism-related gene signature predicts prognosis and treatment response in low-grade glioma.

Cancer cells can induce molecular changes that reshape cellular metabolism, creating specific vulnerabilities for targeted therapeutic interventions. Given the importance of reactive oxygen species (ROS) in tumor development and drug resistance, and the abundance of reduced glutathione (GSH) as the primary cellular antioxidant, we examined an integrated panel of 56 glutathione metabolism-related genes (GMRGs) across diverse cancer types. This analysis revealed a remarkable association between GMRGs and low-grade glioma (LGG) survival. Unsupervised clustering and a GMRGs-based risk score (GS) categorized LGG patients into two groups, linking elevated glutathione metabolism to poorer prognosis and treatment outcomes. Our GS model outperformed established clinical prognostic factors, acting as an independent prognostic factor. GS also exhibited correlations with pro-tumor M2 macrophage infiltration, upregulated immunosuppressive genes, and diminished responses to various cancer therapies. Experimental validation in glioma cell lines confirmed the critical role of glutathione metabolism in glioma cell proliferation and chemoresistance. Our findings highlight the presence of a unique metabolic susceptibility in LGG and introduce a novel GS system as a highly effective tool for predicting the prognosis of LGG.

The prognostic model and immune landscape based on cancer-associated fibroblast features for patients with locally advanced rectal cancer.

Background: This study aimed to construct a nomogram based on CAF features to predict the cancer-specific survival (CSS) rates of locally advanced rectal cancer (LARC) patients. Methods: The EPIC algorithm was employed to calculate the proportion of CAFs. based on the differentially expressed genes between the high and low CAF proportion subgroups, prognostic genes were identified via LASSO and Cox regression analyses. They were then used to construct a prognostic risk signature. Moreover, the GSE39582 and GGSE38832 datasets were used for external validation. Lastly, the level of immune infiltration was evaluated using ssGSEA, ESTIMATE, CIBERSORTx, and TIMER. Results: A higher level of CAF infiltration was associated with a worse prognosis. Additionally, the number of metastasized lymph nodes and distant metastases, as well as the level of immune infiltration were higher in the high CAF proportion subgroup. Five prognostic genes (SMOC2, TUBAL3, C2CD4A, MAP1B, BMP8A) were identified and subsequently incorporated into the prognostic risk signature to predict the 1-, 3-, and 5-year CSS rates in the training and validation sets. Differences in survival rates were also determined in the external validation cohort. Furthermore, independent prognostic factors, including TNM stage and risk score, were combined to established a nomogram. Notably, our results revealed that the proportions of macrophages and neutrophils and the levels of cytokines secreted by M2 macrophages were higher in the high-risk subgroup. Finally, the prognostic genes were significantly associated with the level of immune cell infiltration. Conclusion: Herein, a nomogram based on CAF features was developed to predict the CSS rate of LARC patients. The risk model was capable of reflecting differences in the level of immune cell infiltration.

Prognostic and immunological implications of glutathione metabolism genes in lung adenocarcinoma: A focus on the core gene SMS and its impact on M2 macrophage polarization.

Glutathione metabolism (GM) is a crucial part of various metabolic and pathophysiological processes. However, its role in lung adenocarcinoma (LUAD) has not been comprehensively studied. This study aimed to explore the potential relationship between GM genes, the prognosis, and the immune microenvironment of patients with LUAD. We constructed a risk signature model containing seven GM genes using Lasso combined Cox regression and validated it using six GEO datasets. Our analysis showed that it is an independent prognostic factor. Functional enrichment analysis revealed that the GM genes were significantly enriched in cell proliferation, cell cycle regulation, and metabolic pathways. Clinical and gene expression data of patients with LUAD were obtained from the TCGA database and patients were divided into high- and low-risk groups. The high-risk patient group had a poor prognosis, reduced immune cell infiltration, poor response to immunotherapy, high sensitivity to chemotherapy, and low sensitivity to targeted therapy. Subsequently, single-cell transcriptome analysis using the GSE143423 and GSE127465 datasets revealed that the core SMS gene was highly enriched in M2 Macrophages. Finally, nine GEO datasets and multiple fluorescence staining revealed a correlation between the SMS expression and M2 macrophage polarization. Our prognostic model in which the core SMS gene is closely related to M2 macrophage polarization is expected to become a novel target and strategy for tumor therapy.

Integrated Analysis Construct a Tumor-Associated Macrophage Novel Signature with Promising Implications in Predicting the Prognosis and Immunotherapeutic Response of Gastric Cancer Patients.

BACKGROUND: Gastric cancer (GC) remains one of the most prevalent malignant tumors worldwide. At present, tumor-associated macrophages (TAMs) are essential in the progression, metastasis, and drug resistance of tumors. Therefore, TAMs can be a crucial target for tumor treatment. AIMS: We intended to investigate the TAM characteristics in GC and develop a risk signature based on TAM to predict the prognosis of GC patients. METHODS: The single-cell RNA sequencing (scRNA-seq) and bulk RNA-seq data were acquired from a publicly available database. We utilized the Seurat pipeline to process the scRNA-seq data and determine TAM cell types using marker genes. Univariate Cox regression analysis was utilized to examine TAM-related prognostic genes, and then we employed Lasso-Cox regression analysis, and Multivariate Cox regression analysis established a novel risk profile to forecast the clinical value of the model with a new nomogram combining risk profiles and clinicopathological characteristics. RESULTS: The current study employed scRNA-seq data to identify five TAM clusters in GC, among which four were significantly associated with GC prognosis. Accordingly, we further developed a TAM-related risk signature utilizing nine genes. After evaluation, our model accurately predicted the prognosis of gastric cancer. Generally, GC patients with low TAMS scores exhibited a more favorable prognosis, greater benefits from immunotherapy, and higher levels of immune cell infiltration. CONCLUSIONS: The prognosis of GC can be effectively predicted by TAM-based risk signatures, and the signature may provide a new perspective for comprehensively guiding clinical diagnosis, prediction, and immunotherapy for gastric cancer.

Prognostic Risk Models Using Epithelial Cells Identify ß-Sitosterol as a Potential Therapeutic Target Against Esophageal Squamous Cell Carcinoma.

Background: Esophageal squamous cell carcinoma (ESCC) is an aggressive and fatal malignancy that leads to epithelial cancer. The association between epithelial cell heterogeneity, prognosis, and immune response in this cancer remains uncertain. This study aimed to investigate epithelial cell heterogeneity in ESCC and develop a predictive risk model using the identified cell types. Methods: Single-cell RNA sequencing (scRNA-seq) and differential ESCC gene data were accessed from the Gene Expression Omnibus. Functional enrichment analysis, inferCNV, cell development trajectories, and intercellular communication were analyzed following epithelial cell characterization. Differentially expressed ESCC (n = 773) and epithelial cell marker genes (n = 3407) were intersected to obtain core genes, and epithelial cell-related prognostic genes were identified. LASSO regression analysis was used to construct a prognostic model. The external dataset GSE53624 was used to further validate the stability of the model. Drug sensitivity predictions, and immune cell infiltration were analyzed. Molecular docking clarified the possible therapeutic role of ß-sitosterol in ESCC. Finally, wound healing assay, cell colony, and transwell assay were constructed to detect the effects of the core gene PDLIM2 on ESCC cell proliferation, invasion, and migration. Results: Eight cell clusters were identified, and epithelial cells were categorized into tumor and paratumor groups. The tumor group possessed more chromosomal variants than the paratumor group. Epithelial cells were associated with multiple cell types and significantly correlated with the Wnt, transforming growth factor, and epidermal growth factor signaling pathways. From 231 intersected genes, five core genes were screened for use in the risk model: CTSL, LAPTM4B, MYO10, NCF2, and PDLIM2. These genes may contribute to the cancerous transformation of normal esophageal epithelial cells and thereby act as biomarkers and potential therapeutic targets in patients with ESCC. ß-Sitosterol furthermore displayed excellent docking potential with these genes. Meanwhile, further experiments demonstrated that the gene PDLIM2 plays a major role in the progression of oesophageal squamous carcinoma. Conclusion: We successfully developed a risk model for the prognosis of ESCC based on epithelial cells that addresses the response of ESCC to immunotherapy and offers novel cancer treatment options.

Multi-omics and immunogenomics analysis revealed PFKFB3 as a targetable hallmark and mediates sunitinib resistance in papillary renal cell carcinoma: in silico study with laboratory verification.

Glycolysis-related metabolic reprogramming is a central hallmark of human cancers, especially in renal cell carcinoma. However, the regulatory function of glycolytic signature in papillary RCC has not been well elucidated. In the present study, the glycolysis-immune predictive signature was constructed and validated using WGCNA, glycolysis-immune clustering analysis. PPI network of DEGs was constructed and visualized. Functional enrichments and patients' overall survival were analyzed. QRT-PCR experiments were performed to detect hub genes' expression and distribution, siRNA technology was used to silence targeted genes; cell proliferation and migration assays were applied to evaluate the biological function. Glucose concentration, lactate secretion, and ATP production were measured. Glycolysis-Immune Related Prognostic Index (GIRPI) was constructed and combined analyzed with single-cell RNA-seq. High-GIRPI signature predicted significantly poorer outcomes and relevant clinical features of pRCC patients. Moreover, GIRPI also participated in several pathways, which affected tumor immune microenvironment and provided potential therapeutic strategy. As a key glycolysis regulator, PFKFB3 could promote renal cancer cell proliferation and migration in vitro. Blocking of PFKFB3 by selective inhibitor PFK-015 or glycolytic inhibitor 2-DG significantly restrained renal cancer cells' neoplastic potential. PFK-015 and sunitinib could synergistically inhibit pRCC cells proliferation. Glycolysis-Immune Risk Signature is closely associated with pRCC prognosis, progression, immune infiltration, and therapeutic response. PFKFB3 may serve as a pivotal glycolysis regulator and mediates Sunitinib resistance in pRCC patients.

Development of a prognostic model based on different disulfidptosis related genes typing for kidney renal clear cell carcinoma.

Background: Kidney renal clear cell carcinoma (KIRC) is a common and clinically significant subtype of kidney cancer. A potential therapeutic target in KIRC is disulfidptosis, a novel mode of cell death induced by disulfide stress. The aim of this study was to develop a prognostic model to explore the clinical significance of different disulfidptosis gene typings from KIRC. Methods: A comprehensive analysis of the chromosomal localization, expression patterns, mutational landscape, copy number variations, and prognostic significance of 10 disulfide death genes was conducted. Patients were categorized into distinct subtypes using the Non-negative Matrix Factorization (NMF) typing method based on disulfidptosis gene expression patterns. Weighted Gene Co-expression Network Analysis (WGCNA) was used on the KIRC dataset to identify differentially expressed genes between subtype clusters. A risk signature was created using LASSO-Cox regression and validated by survival analysis. An interaction between risk score and immune cell infiltration, tumor microenvironment characteristics and pathway enrichment analysis were investigated. Results: Initial findings highlight the differential expression of specific DRGs in KIRC, with genomic instability and somatic mutation analysis revealing key insights into their role in cancer progression. NMF clustering differentiates KIRC patients into subgroups with distinct survival outcomes and immune profiles, and hierarchical clustering identifies gene modules associated with key biological and clinical parameters, leading to the development of a risk stratification model (LRP8, RNASE2, CLIP4, HAS2, SLC22A11, and KCTD12) validated by survival analysis and predictive of immune infiltration and drug sensitivity. Pathway enrichment analysis further delineates the differential molecular pathways between high-risk and low-risk patients, offering potential targets for personalized treatment. Lastly, differential expression analysis of model genes between normal and KIRC cells provides insights into the molecular mechanisms underlying KIRC, highlighting potential biomarkers and therapeutic targets. Conclusion: This study contributes to the understanding of KIRC and provides a potential prognostic model using disulfidptosis gene for personalized management in KIRC patients. The risk signature shows clinical applicability and sheds light on the biological mechanisms associated with disulfide-induced cell death.

A novel risk signature based on liquid-liquid phase separation-related genes reveals prognostic and tumour microenvironmental features in clear cell renal cell carcinoma.

BACKGROUND: Clear cell renal cell carcinoma(ccRCC) is one of the most common malignancies. However, there are still many barriers to its underlying causes, early diagnostic techniques and therapeutic approaches. MATERIALS AND METHODS: The Cancer Genome Atlas (TCGA)- Kidney renal clear cell (KIRC) cohort differentially analysed liquid-liquid phase separation (LLPS)-related genes from the DrLLPS website. Univariate and multivariate Cox regression analyses and LASSO regression analyses were used to construct prognostic models. The E-MTAB-1980 cohort was used for external validation. Then, potential functions, immune infiltration analysis, and mutational landscapes were analysed for the high-risk and low-risk groups. Finally, quantitative real-time polymerase chain reaction (qRT-PCR) experiments as well as single-cell analyses validated the genes key to the model. RESULTS: We screened 174 LLPS-related genes in ccRCC and constructed a risk signature consisting of five genes (CLIC5, MXD3, NUF2, PABPC1L, PLK1). The high-risk group was found to be associated with worse prognosis in different subgroups. A nomogram constructed by combining age and tumour stage had a strong predictive power for the prognosis of ccRCC patients. In addition, there were differences in pathway enrichment, immune cell infiltration, and mutational landscapes between the two groups. The results of qRT-PCR in renal cancer cell lines and renal cancer tissues were consistent with the biosignature prediction. Three single-cell data of GSE159115, GSE139555, and GSE121636 were analysed for differences in the presence of these five genes in different cells. CONCLUSIONS: We developed a risk signature constructed based on the five LLPS-related genes and can have a high ability to predict the prognosis of ccRCC patients, further providing a strong support for clinical decision-making.

Immunogenic cell death related mRNAs associated signature to predict immunotherapeutic response in osteosarcoma.

Background: Immunogenic cell death (ICD) is related to cancer prognosis, which has a synergic effect in combination with chemotherapy or immunotherapy. Yet, the relationship between ICD and osteosarcoma remained unclear. Materials and methods: Three osteosarcoma datasets including therapeutically applicable research to generate effective treatments (TARGET), GSE126209 and GSE21257 datasets were included. A protein-protein interaction network was constructed based on ICD-related genes. We performed unsupervised consensus clustering to classify molecular subtypes (clusters). Survival analysis, Estimation of stromal and immune cells in malignant tumour tissues using expression data (ESTIMATE), Cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT), and differential analysis were employed to characterize the molecular differences between different clusters. Univariate Cox regression analysis was conducted to confirm prognostic genes. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) was used to demonstrate the aberrant expression of ICD-correlated signature genes in osteosarcoma. A series of cellular experiments, including cell counting kit-8 (CCK-8), transwell, and flow cytometry, were used to demonstrate the regulatory role of key genes in the ICD model on the malignant phenotype of osteosarcoma. Results: Three clusters (cluster1, 2, 3) were constructed and they showed distinct overall survival and immune infiltration. ICD-related genes were highly expressed in cluster1. Moreover, Cluster1 had the best prognosis, high immune score and high expression of human leukocyte antigen (HLA)-related genes. TLR4, LY96, IFNGR1, CD4, and CASP1 were identified as prognostic genes for establishing an ICD-related risk signature. According to the risk signature, two risk groups (high and low risks) showing differential prognosis and response to immunotherapy. The low risks group had a better prognosis but was not sensitive to immunotherapy. Molecular assays verified that prognostic genes were abnormally under-expressed in osteosarcoma. Cellular assays demonstrated that LY96, the most significantly down-regulated gene in osteosarcoma, inhibited the migration, invasion, and proliferation phenotypes of osteosarcoma cells and prolonged the cell cycle. Analysis of oxidative stress related pathway enrichment in tumor microenvironment was conducted by single-sample gene set enrichment analysis (ssGSEA). Conclusions: This study demonstrated the prognostic significance of ICD-correlated genes in osteosarcoma patients. The five-gene risk signature facilitate prognostic evaluation and prediction of osteosarcoma patients' response to immunotherapy. The risk signature also offered a possibility for the exploit of novel ICD-related treatment.