Search

Pancreatic cancer (PAC) is a highly fatal and aggressive type of cancer. Hypoxia is a common feature of PAC. The aim of this study was to develop a hypoxia status-related prognostic model for predicting the survival outcomes in PAC. The data sets of PAC from The Cancer Genome Atlas and the International Cancer Genome Consortium were used to construct and validate the signature. A 6 hypoxia status-related differential expression genes prognostic model for predicting the survival outcomes was established. The Kaplan–Meier analysis and Received operating characteristic curve indicated the good performance of the signature at predicting overall survival. Univariate and Multivariate Cox regression revealed that the signature was an independent prognostic factor in PAC. Weighted Gene Co-expression Network Analysis and immune infiltration analysis indicated that Immune-related pathways and immune cell infiltration was mostly enriched in the low-risk group, which presented a better prognosis. We also evaluated the predictive of the signature for immunotherapy and chemoradiotherapy. Risk gene LY6D may be a potential prognostic predictor of PAC. This model can be used as an independent prognostic factor for predicting clinical outcomes and a possible classifier for response to chemotherapy. [ABSTRACT FROM AUTHOR]

Background: Uveal melanoma (UVM) is the most common primary ocular malignancy, with a wide range of symptoms and outcomes. The programmed cell death (PCD) plays an important role in tumor development, diagnosis, and prognosis. There is still no research on the relationship between PCD-related genes and UVM. A novel PCD-associated prognostic model is urgently needed to improve treatment strategies. Objective: We aim to screen PCD-related prognostic signature and investigate its proliferation ability and apoptosis in UVM cells. Methods: The clinical information and RNA-seq data of the UVM patients were collected from the TCGA cohort. All the patients were classified using consensus clustering by the selected PCD-related genes. After univariate Cox regression and PPI network analysis, the prognostic PCD-related genes were then submitted to the LASSO regression analysis to build a prognostic model. The level of immune infiltration of 8-PCD signature in high- and low-risk patients was analyzed using xCell. The prediction on chemotherapy and immunotherapy response in UVM patients was assessed by GDSC and TIDE algorithm. CCK-8, western blot and Annexin V-FITC/PI staining were used to explore the roles of HMOX1 in UVM cells. Results: A total of 8-PCD signature was constructed and the risk score of the PCD signature was negatively correlated with the overall survival, indicating strong predictive ability and independent prognostic value. The risk score was positively correlated with CD8 Tcm, CD8 Tem and Th2 cells. Immune cells in high-risk group had poorer overall survival. The drug sensitivity demonstrated that cisplatin might impact the progression of UVM and better immunotherapy responsiveness in the high-risk group. Finally, Overespression HMOX1 (OE-HMOX1) decreased the cell viability and induced apoptosis in UVM cells. Recuse experiment results showed that ferrostatin-1 (fer-1) protected MP65 cells from apoptosis and necrosis caused by OE-HMOX1. Conclusion: The PCD signature may have a significant role in the tumor microenvironment, clinicopathological characteristics, prognosis and drug sensitivity. More importantly, HMOX1 depletion greatly induced tumor cell growth and inhibited cell apoptosis and fer-1 protected UVM cells from apoptosis and necrosis induced by OE-HMOX1. This work provides a foundation for effective therapeutic strategy in tumour treatment. [ABSTRACT FROM AUTHOR]

The progression of acute myeloid leukemia (AML) is influenced by the immune microenvironment in the bone marrow and dysregulated intracellular competing endogenous RNA (ceRNA) networks. Our study utilized data from UCSC Xena, The Cancer Genome Atlas Program, the Gene Expression Omnibus, and the Immunology Database and Analysis Portal. Using Cox regression analysis, we identified an immune-related prognostic signature. Genomic analysis of prognostic messenger RNA (mRNA) was conducted through Gene Set Cancer Analysis (GSCA), and a prognostic ceRNA network was constructed using the Encyclopedia of RNA Interactomes. Correlations between signature mRNAs and immune cell infiltration, checkpoints, and drug sensitivity were assessed using R software, gene expression profiling interactive analysis (GEPIA), and CellMiner, respectively. Adhering to the ceRNA hypothesis, we established a potential long noncoding RNA (lncRNA)/microRNA (miRNA)/mRNA regulatory axis. Our findings pinpointed 9 immune-related prognostic mRNAs (KIR2DL1, CSRP1, APOBEC3G, CKLF, PLXNC1, PNOC, ANGPT1, IL1R2, and IL3RA). GSCA analysis revealed the impact of copy number variations and methylation on AML. The ceRNA network comprised 14 prognostic differentially expressed lncRNAs (DE-lncRNAs), 6 prognostic DE-miRNAs, and 3 prognostic immune-related DE-mRNAs. Correlation analyses linked these mRNAs' expression to 22 immune cell types and 6 immune checkpoints, with potential sensitivity to 27 antitumor drugs. Finally, we identified a potential LINC00963/hsa-miR-431-5p/CSRP1 axis. This study offers innovative insights for AML diagnosis and treatment through a novel immune-related signature and ceRNA axis. Identified novel biomarkers, including 2 mRNAs (CKLF, PNOC), 1 miRNA (hsa-miR-323a-3p), and 10 lncRNAs (SNHG25, LINC01857, AL390728.6, AC127024.5, Z83843.1, AP002884.1, AC007038.1, AC112512, AC020659.1, AC005921.3) present promising candidates as potential targets for precision medicine, contributing to the ongoing advancements in the field. [ABSTRACT FROM AUTHOR]

This study intended to develop a new immunogenic cell death (ICD)-related prognostic signature for colorectal cancer (CRC) patients. The Non-Negative Matrix Factorization (NMF) algorithm was adopted to cluster tumor samples based on ICD gene expression to obtain ICD-related subtypes. Survival analysis and immune microenvironment analysis were conducted among different subtypes. Regression analysis was used to construct the model. Based on riskscore median, cancer patients were classified into high and low risk groups, and independent prognostic ability of the model was analyzed. The CIBERSORT algorithm was adopted to determine the immune infiltration level of both groups. We analyzed the differential genes between cluster 4 and cluster 1–3 and obtained 12 genes with the best prognostic features finally (NLGN1, SLC30A3, C3orf20, ADAD2, ATOH1, ATP6V1B1, KCNQ2, MUCL3, RGCC, CLEC17A, COL6A5, and INSL4). In addition, patients with lower risk had higher levels of infiltration of most immune cells, lower Tumor Immune Dysfunction and Exclusion (TIDE) level and higher immunophenscore (IPS) level than those with higher risk. This study constructed and validated the ICD feature signature predicting CRC prognosis and provide a reference criteria for guiding the prognosis and immunotherapy of CRC cancer patients. [ABSTRACT FROM AUTHOR]

CD8+ T lymphocytes are critical in the immune response against neoplasms, yet the prognostic relevance of CD8+ T cell-associated genes in hepatocellular carcinoma (HCC) is not fully understood. We sourced single-cell RNA-sequencing (scRNA-seq) and bulk RNA-seq data for HCC from the GSE98638 dataset and The Cancer Genome Atlas (TCGA) repository. We utilized Weighted Gene Correlation Network Analysis (WGCNA) to identify CD8+ T cell-related genes. A clinical prognostic model for risk stratification was then constructed via Cox-Lasso regression analysis. The Immunophenotypic Score (IPS) was utilized to evaluate the potential of immunotherapeutic interventions in the categorized cohorts. Validation of the expression of CD8+ T cell-associated risk genes was performed using quantitative reverse transcription PCR (qRT-PCR). Integrating scRNA-seq with RNA-seq data, we identified five CD8+ T cell-related signature genes: IKBKE, ATP1B3, MSC, ADA, and BATF. Notably, HCC patients in the high-risk group had markedly decreased overall survival. Elevated infiltration levels of CD8+ T cells, B cells, and macrophages were observed in the high-risk group. Moreover, there was a positive correlation between the risk score and immune checkpoints (ICPs), including PDCD1, CD274, and CTLA4. Patients within the high-risk group subject to PD1 and CTLA4 blockade exhibited higher IPS levels. Additionally, the expression of the five risk genes was upregulated in HCC cell lines and tissues compared to normal cells and tissues. Our findings establish a prognostic signature based on CD8+ T cells, offering a potent predictive model for clinical outcomes and responsiveness to immunotherapy in HCC patients. [ABSTRACT FROM AUTHOR]

The tumor microenvironment (TME) is crucial in tumor development, metastasis, and response to immunotherapy. DNA methylation can regulate the TME without altering the DNA sequence. However, research on the methylation-driven TME in clear-cell renal cell carcinoma (ccRCC) is still lacking. In this study, integrated DNA methylation and RNA-seq data were used to explore methylation-driven genes (MDGs). Immune scores were calculated using the ESTIMATE, which was employed to identify TME-related genes. A new signature connected with methylation-regulated TME using univariate, multivariate Cox regression and LASSO regression analyses was developed. This signature consists of four TME-MDGs, including AJAP1, HOXB9, MYH14, and SLC6A19, which exhibit high methylation and low expression in tumors. Validation was performed using qRT-PCR which confirmed their downregulation in ccRCC clinical samples. Additionally, the signature demonstrated stable predictive performance in different subtypes of ccRCC. Risk scores are positively correlated with TMN stages, immune cell infiltration, tumor mutation burden, and adverse outcomes of immunotherapy. Interestingly, the expression of four TME-MDGs are highly correlated with the sensitivity of first-line drugs in ccRCC treatment, especially pazopanib. Molecular docking indicates a high affinity binding between the proteins and pazopanib. In summary, our study elucidates the comprehensive role of methylation-driven TME in ccRCC, aiding in identifying patients sensitive to immunotherapy and targeted therapy, and providing new therapeutic targets for ccRCC treatment. [ABSTRACT FROM AUTHOR]

Urinary cancer is synonymous with clear cell renal cell carcinoma (ccRCC). Unfortunately, existing treatments for this illness are ineffective and unpromising. Finding novel ccRCC biomarkers is crucial to creating successful treatments. The Cancer Genome Atlas provided clear cell renal cell carcinoma transcriptome data. Functional enrichment analysis was performed on ccRCC and control samples' differentially expressed N6‐methyladenosine RNA methylation and ferroptosis‐related genes (DEMFRGs). Machine learning was used to find and model ccRCC patients' predicted genes. A nomogram was created for clear cell renal cell carcinoma patients. Prognostic genes were enriched. We examined patients' immune profiles by risk score. Our prognostic genes predicted ccRCC treatment drugs. We found 37 DEMFRGs by comparing 1913 differentially expressed ccRCC genes to 202 m6A RNA methylation FRGs. Functional enrichment analysis showed that hypoxia‐induced cell death and metabolism pathways were the most differentially expressed methylation functional regulating genes. Five prognostic genes were found by machine learning: TRIB3, CHAC1, NNMT, EGFR, and SLC1A4. An advanced renal cell carcinoma nomogram with age and risk score accurately predicted the outcome. These five prognostic genes were linked to various cancers. Immunological cell number and checkpoint expression differed between high‐ and low‐risk groups. The risk model successfully predicted immunotherapy outcome, showing high‐risk individuals had poor results. NIACIN, TAE‐684, ROCILETINIB, and others treat ccRCC. We found ccRCC prognostic genes that work. This discovery may lead to new ccRCC treatments. [ABSTRACT FROM AUTHOR]

Background: Disulfide, an essential compounds family, has diverse biological activity and can affect the dynamic balance between physiological and pathological states. A recently published study found that aberrant accumulation of disulfide had a lethal effect on cells. This mechanism of cell death, named disulfidptosis, differs from other known cell death mechanisms, including cuproptosis, apoptosis, necroptosis, and pyroptosis. The relationship between disulfidptosis and development of cancer, in particular endometrial carcinoma, remains unclear. Methods: To address this knowledge gap, we performed a preliminary analysis of samples from The Cancer Genome Atlas database. The samples were divided equally into a training group and a test group. A total of 2308 differentially expressed genes were extracted, and 11 were used to construct a prognostic model. Results: Based on the risk score calculated using the prognostic model, the samples were divided into a high‐risk group and a low‐risk group. Survival time, tumor mutation burden, and microsatellite instability scores differed significantly between the two groups. Furthermore, a between‐group difference in treatment effect was predicted. Comparison with other models in the literature indicated that this prognostic model had better predictive anility. Conclusion: The results of this study provide a general framework for understanding the relationship between disulfidptosis and endometrial cancer that could be used for clinical evaluation and selection of appropriate personalized treatment strategies. [ABSTRACT FROM AUTHOR]

Background: Patients with lung adenocarcinoma (LUAD) often develop a poor prognosis. Currently, researches on prognostic and immunotherapeutic capacity of aneuploidy-related genes in LUAD are limited. Methods: Genes related to aneuploidy were screened based on bulk RNA sequencing data from public databases using Spearman method. Next, univariate Cox and Lasso regression analyses were performed to establish an aneuploidy-related riskscore (ARS) model. Results derived from bioinformatics analysis were further validated using cellular experiments. In addition, typical LUAD cells were identified by subtype clustering, followed by SCENIC and intercellular communication analyses. Finally, ESTIMATE, ssGSEA and CIBERSORT algorithms were employed to analyze the potential relationship between ARS and tumor immune environment. Results: A five-gene ARS signature was developed. These genes were abnormally high-expressed in LUAD cell lines, and in particular the high expression of CKS1B promoted the proliferative, migratory and invasive phenotypes of LUAD cell lines. Low ARS group had longer overall survival time, higher degrees of inflammatory infiltration, and could benefit more from receiving immunotherapy. Patients in low ASR group responded more actively to traditional chemotherapy drugs (Erlotinib and Roscovitine). The scRNA-seq analysis annotated 17 cell subpopulations into seven cell clusters. Core transcription factors (TFs) such as CREB3L1 and CEBPD were enriched in high ARS cell group, while TFs such as BCLAF1 and UQCRB were enriched in low ARS cell group. CellChat analysis revealed that high ARS cell groups communicated with immune cells via SPP1 (ITGA4-ITGB1) and MK (MDK-NCl) signaling pathways. Conclusion: In this research, integrative analysis based on the ARS model provided a potential direction for improving the diagnosis and treatment of LUAD. [ABSTRACT FROM AUTHOR]

Background: Breast carcinoma (BRCA) is a life-threatening malignancy in women and shows a poor prognosis. Cuproptosis is a novel mode of cell death but its relationship with BRCA is unclear. This study attempted to develop a cuproptosis-relevant prognostic gene signature for BRCA. Methods: Cuproptosis-relevant subtypes of BRCA were obtained by consensus clustering. Differential expression analysis was implemented using the 'limma' package. Univariate Cox and multivariate Cox analyses were performed to determine a cuproptosis-relevant prognostic gene signature. The signature was constructed and validated in distinct datasets. Gene set variation analysis (GSVA) and gene set enrichment analysis (GSEA) were also conducted using the prognostic signature to uncover the underlying molecular mechanisms. ESTIMATE and CIBERSORT algorithms were applied to probe the linkage between the gene signature and tumor microenvironment (TME). Immunotherapy responsiveness was assessed using the Tumor Immune Dysfunction and Exclusion (TIDE) web tool. Real-time quantitative PCR (RT-qPCR) was performed to detect the expressions of cuproptosis-relevant prognostic genes in breast cancer cell lines. Results: Thirty-eight cuproptosis-associated differentially expressed genes (DEGs) in BRCA were mined by consensus clustering and differential expression analysis. Based on univariate Cox and multivariate Cox analyses, six cuproptosis-relevant prognostic genes, namely SAA1, KRT17, VAV3, IGHG1, TFF1, and CLEC3A, were mined to establish a corresponding signature. The signature was validated using external validation sets. GSVA and GSEA showed that multiple cell cycle-linked and immune-related pathways along with biological processes were associated with the signature. The results ESTIMATE and CIBERSORT analyses revealed significantly different TMEs between the two Cusig score subgroups. Finally, RT-qPCR analysis of cell lines further confirmed the expressional trends of SAA1, KRT17, IGHG1, and CLEC3A. Conclusion: Taken together, we constructed a signature for projecting the overall survival of BRCA patients and our findings authenticated the cuproptosis-relevant prognostic genes, which are expected to provide a basis for developing prognostic molecular biomarkers and an in-depth understanding of the relationship between cuproptosis and BRCA. [ABSTRACT FROM AUTHOR]

Objective: Patients with clear cell renal carcinoma (ccRCC) typically have poor prognosis. Glycosylation modification plays an important role in ccRCC. This study aimed to develop a novel signature for predicting ccRCC prognosis based on glycosylation-related genes (GRGs). Methods: Differentially expressed GRGs (DE_GRGs) were identified using The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus and used for constructing the risk model. The function of key genes was validated. Results: Twenty-two DE_GRGs were intersected between GSE53757 and TCGA. Patients with ccRCC were divided into two clusters based on the expression profile of these DE_GRGs. Significant differences in the infiltration of 10 immune cell types were observed between two subclusters. Subsequently, the prognostic signatures of nine DE_GRGs (CHST9, COLGALT1, FUT3, FUT6, HS3ST2, POMGNT2, ST8SIA4, UGT3A1, and UGT8) were established. Patients in the high-risk group showed a poorer prognosis relative to the low-risk group. According to univariate and multivariate Cox regression analyses, the risk score, stage, and grade could be independent prognostic factors. A nomogram incorporating information on gender, age, risk group, TNM stage, and clinical stage showed accurate prediction in the survival probability. Except for CHST9, HS3ST2, and ST8SIA4, expression patterns of the remaining 6 DE_GRGs in Caki-1 and 786-O cells were confirmed by quantitative real-time PCR. FUT6 overexpression resulted in the inhibition of proliferation, migration, and invasion of ccRCC cells. Conclusion: This study established a nine-DE_GRG-based prognostic signature, which independently predicted ccRCC prognosis. This finding emphasizes that GRGs are stratification factors for the precise prognosis of ccRCC.

ABSTRACT Background In the era of immunotherapy, there is a critical need for effective biomarkers to improve outcome prediction and guide treatment decisions for patients with lung squamous cell carcinoma (LUSC). We hypothesized that the immune contexture of LUSC may be influenced by tumor intrinsic events, such as autophagy. Aims We aimed to develop an autophagy‐related risk signature and assess its predictive value for immune phenotype. Methods and Results Expression profiles of autophagy‐related genes (ARGs) in LUSC samples were obtained from the TCGA and GEO databases. Survival analyses were conducted to identify survival‐related ARGs and construct a risk signature using the Random Forest algorithm. Four ARGs (CFLAR, RGS19, PINK1, and CTSD) with the most significant prognostic value were selected to construct the risk signature. Patients in the high‐risk group exhibited worse prognosis than those in the low‐risk group (p

BackgroundNatural killer (NK) cells play crucial roles in mediating anti-cancer activity in breast cancer (BRCA). However, the potential of NK cell-related molecules in predicting BRCA outcomes and guiding personalized therapy remains largely unexplored. This study focused on developing a prognostic and therapeutic prediction model for BRCA by incorporating NK cell-related genes.MethodsThe data analyzed primarily originated from the TCGA and GEO databases. The prognostic role of NK cells was evaluated, and marker genes of NK cells were identified via single-cell analysis. Module genes closely associated with immunotherapy resistance were identified by bulk transcriptome-based weighted correlation network analysis (WGCNA). Following taking intersection and LASSO regression, NK-related genes (NKRGs) relevant to BRCA prognosis were screened, and the NK-related prognostic signature was subsequently constructed. Analyses were further expanded to clinicopathological relevance, GSEA, tumor microenvironment (TME) analysis, immune function, immunotherapy responsiveness, and chemotherapeutics. Key NKRGs were screened by machine learning and validated by spatial transcriptomics (ST) and immunohistochemistry (IHC).ResultsTumor-infiltrating NK cells are a favorable prognostic factor in BRCA. By combining scRNA-seq and bulk transcriptomic analyses, we identified 7 NK-related prognostic NKRGs (CCL5, EFHD2, KLRB1, C1S, SOCS3, IRF1, and CCND2) and developed an NK-related risk scoring (NKRS) system. The prognostic reliability of NKRS was verified through survival and clinical relevance analyses across multiple cohorts. NKRS also demonstrated robust predictive power in various aspects, including TME landscape, immune functions, immunotherapy responses, and chemotherapeutic sensitivity. Additionally, KLRB1 and CCND2 emerged as key prognostic NKRGs identified through machine learning and external validation, with their expression correlation with NK cells confirmed in BRCA specimens by ST and IHC.ConclusionsWe developed a novel NK-related gene signature that has proven valuable for evaluating prognosis and treatment response in BRCA, expecting to advance precision medicine of BRCA.

Endoplasmic reticulum stress (ERS) and cuproptosis have remarkable effects on hepatocellular carcinoma (HCC) leading to a poor prognosis. The current study aimed to explore credible signature for predicting the prognosis of HCC based on ERS and cuproptosis-related lncRNAs. In our study, clinical and transcriptomic profiles of HCC patients were obtained from the Cancer Genome Atlas (TCGA) database. An ERS and cuproptosis-related lncRNA prognostic signature, including NRAV, SNHG3, LINC00839 and AC004687.1, was determined by correlation tests, Cox regression analysis, least absolute shrinkage, and selection operator (LASSO) methods. Survival and predictive value were evaluated using Kaplan–Meier and receiver operating characteristic (ROC) curves, while calibration and nomograms curves were developed. Besides the enrichment analyses for ERS and cuproptosis-related lncRNAs, mutational status and immune status were assessed with TMB and ESTIMATE. Additionally, consensus cluster analysis was employed to compare cancer subtype differences, while drug sensitivity and immunologic efficacy were evaluated for further exploration. qRT-PCR and CCK-8 were utilized to verify the alteration of the prognostic lncRNAs expression and proliferation in vitro. High-risk groups exhibited poorer prognosis. The signature exhibited robust predictive value as an independent prognostic indicator and showed significant correlation with clinicopathological features. In the enriched analysis, biological membrane pathways were enriched. Low-risk patients had lower TMB and higher immune status. A cluster analysis revealed that cluster 2 had the best clinical immunological efficacy and most active immune function. In brief, our constructed signature with ERS and cuproptosis-related lncRNAs was associated survival outcomes of HCC, and can be used to predict the clinical classification and curative effect.

Objective: Hepatocellular carcinoma (HCC) ranks as the third leading cause of cancer-related deaths, constituting 75%–85 % of all primary liver cancers. The objective of this study was to develop a necroptosis-related gene signature using single-cell and bulk RNA sequencing to predict HCC patient prognoses. Methods: A total of 25 key necroptosis regulators were identified from previous literature. We evaluated the necroptosis scores of different cell types using single-cell sequencing data from HCC and analyzed 168 necroptosis-related genes. The Cancer Genome Atlas Liver Hepatocellular Carcinoma (TCGA-LIHC) dataset served as the training set for establishing a novel necroptosis-related gene risk signature, employing univariate and multivariate Cox regression analyses. Additionally, the study examined the model's relevance in immunity and immunotherapy, and predicted chemosensitivity in HCC patients based on the gene signature. The key genes were validated by the biological experiments. Results: Compared to other cell types, hepatoma cells displayed the lowest necroptosis scores. A new six-gene necroptosis-related signature (S100A11, MAGEC2, MAGEA6, CTP2C9, SOX4, AKR1B10) was developed using the TCGA database and validated in the ICGC database. Patients in the high-risk category had poorer prognoses, with the risk score serving as an independent prognostic indicator beyond other clinical factors. These high-risk patients also exhibited greater immune infiltration but demonstrated a weaker anti-tumor response due to elevated expression of immune checkpoints. Pathways involving hypoxia, glycolysis, and P53, as well as the frequency of P53 somatic mutations, were notably heightened in the high-risk group. Additionally, the six genes in the model showed significantly different mRNA expression in hepatoma cell lines compared to normal hepatocytes, with the role of MAGEA6 in liver cancer being elucidated through critical experiments. Conclusions: This study successfully developed a six-gene necroptosis-related signature to predict prognoses in HCC patients. It further explored the roles of necroptosis in hepatoma cells and the tumor microenvironment.

Background: Hepatocellular carcinoma (HCC) is among the most prevalent digestive system malignancies and is associated with a poor prognosis. Necroptosis, a form of regulated death mediated by death receptors, exhibits characteristics of both necrosis and apoptosis. Long non-coding RNAs (lncRNAs) have been identified as crucial regulators in tumor necroptosis. This study aims to identify the necroptosis-related lncRNAs (np-lncRNA) in HCC and investigate their relationships with prognosis. Method: The RNA-sequencing data, along with clinicopathological and survival information of HCC patients were sourced from The Cancer Genome Atlas (TCGA) database. The np-lncRNAs were analyzed to assess their potential in predicting HCC prognosis. Prognostic signatures related to necroptosis were constructed using stepwise multivariate Cox regression analysis. The prognosis of patients was compared using Kaplan-Meier (KM) analysis. The accuracy of the prognostic signature was evaluated using Receiver operating characteristic (ROC) analysis and decision curve analysis (DCA). Quantitative real-time polymerase chain reaction(qPCR) was employed to validate the lncRNAs expression levels of lncRNAs among samples from an independent cohort. Results: The np-lncRNAs ZFPM2-AS1, AC099850.3, BACE1-AS, KDM4A-AS1 and MKLN1-AS were identified as potential prognostic biomarkers. The prognostic signature constructed from these np-lncRNAs achieved an Area Under the Curve (AUC) of 0.773. Based on the risk score derived from the signature, patients were divided into two groups, with the high-risk group exhibiting poorer overall survival. Gene Set Enrichment Analysis (GSEA) revealed significantly different between the low risk and high risk groups in tumor-related pathways (such as mTOR, MAPK and p53 signaling pathways) and immune-related functions (like T cell receptor signaling pathway and natural killer cell mediated cytotoxicity). The increased expression of np-lncRNAs was confirmed in another independent HCC cohort. Conclusions: This signature offers a dependable method for forecasting the prognosis of HCC patients. Our findings indicate a subset of np-lncRNA biomarkers that could be utilized for prognosis prediction and personalized treatment strategies of HCC patients.

Abstract The prognosis for patients with colorectal cancer (CRC) remains worse than expected due to metastasis, recurrence, and resistance to chemotherapy. Colorectal cancer stem cells (CRCSCs) play a vital role in tumor metastasis, recurrence, and chemotherapy resistance. However, there are currently no prognostic markers based on CRCSCs-related genes available for clinical use. In this study, single-cell transcriptome sequencing was employed to distinguish cancer stem cells (CSCs) in the CRC microenvironment and analyze their properties at the single-cell level. Subsequently, data from TCGA and GEO databases were utilized to develop a prognostic risk model for CRCSCs-related genes and validate its diagnostic performance. Additionally, functional enrichment, immune response, and chemotherapeutic drug sensitivity of the relevant genes in the risk model were investigated. Lastly, the key gene RPS17 in the risk model was identified as a potential prognostic marker and therapeutic target for further comprehensive studies. Our findings provide new insights into the prognostic treatment of CRC and offer novel perspectives for a systematic and comprehensive understanding of CRC development.

Abstract Accumulating evidence demonstrates that lncRNAs are involved in the regulation of the immune microenvironment and early tumor development. Immunogenic cell death occurs mainly through the release or increase of tumor-associated antigen and tumor-specific antigen, exposing “danger signals” to stimulate the body’s immune response. Given the recent development of immunotherapy in lung adenocarcinoma, we explored the role of tumor immunogenic cell death-related lncRNAs in lung adenocarcinoma for prognosis and immunotherapy benefit, which has never been uncovered yet. Based on the lung adenocarcinoma cohorts from the TCGA database and GEO database, the study developed the immunogenic cell death index signature by several machine learning algorithms and then validated the signature for prognosis and immunotherapy benefit of lung adenocarcinoma patients, which had a more stable performance compared with published signatures in predicting the prognosis, and demonstrated predictive value for benefiting from immunotherapy in multiple cohorts of multiple cancers, and also guided the utilization of chemotherapy drugs.

Abstract Background As a highly heterogeneous tumor, non-small cell lung cancer (NSCLC) is famous for its high incidence and mortality worldwide. Smoking can cause genetic changes, which leading to the occurrence and progress of NSCLC. Nevertheless, the function of smoking-related genes in NSCLC needs more research. Methods We downloaded transcriptome data and clinicopathological parameters from Gene Expression Omnibus (GEO) databases, and screened smoking-related genes. Lasso regression were applied to establish the 7-gene signature. The associations between the 7-gene signature and immune microenvironment analysis, survival analysis, drug sensitivity analysis and enriched molecular pathways were studied. Ultimately, cell function experiments were conducted to research the function of FCGBP in NSCLC. Results Through 7-gene signature, NSCLC samples were classified into high-risk group (HRG) and low-risk group (LRG). Significant difference in overall survival (OS) between HRG and LRG was found. Nomograms and ROC curves indicated that the 7-gene signature has a stable ability in predicting prognosis. Through the analysis of immune microenvironment, we found that LRG patients had better tumor immune activation. FCGBP showed the highest mutation frequency among the seven prognostic smoking related genes (LRRC31, HPGD, FCGBP, SPINK5, CYP24A1, S100P and FGG), and was notable down-regulated in NSCLC smokers compared with non-smoking NSCLC patients. The cell experiments confirmed that FCGBP knockdown promoting proliferation, migration, and invasion in NSCLC cells. Conclusion This smoking-related prognostic signature represents a promising tool for assessing prognosis and tumor microenvironment in smokers with NSCLC. The role of FCGBP in NSCLC was found by cell experiments, which can be served as diagnostic biomarker and immunotherapy target for NSCLC.

Abstract Purpose Cuproptosis, programmed cell death by intracellular copper-mediated lipoylated protein aggregation, is involved in various tumorigenesis and drug resistance abilities by mediating the tumor microenvironment. Previous studies have demonstrated that serum copper levels are higher in OC patients than in normal subjects. However, the exact relationship between cuproptosis and ovarian cancer progression remains to be further elucidated. Methods The Cancer Genome Atlas (TCGA) and gene expression omnibus (GEO) datasets were utilized to establish a cuproptosis-related prognostic signature in ovarian cancer. Subsequently, the bulk RNA-seq analysis and single-cell RNA-seq analysis were used to identify the relationship between signature with immune cell infiltration, chemotherapy, and cuproptosis-related scoring (CuRS) system. Finally, the potential biological functional roles of target genes in cuproptosis were validated in vitro. Results By using LASSO-Cox regression analysis to establish the cuproptosis-related prognostic model, our works demonstrated the accuracy and efficiency of our model in the TCGA (583 OC patients) and GEO (260 OC patients) OC cohorts, and the high-scoring groups showed worse survival outcomes. Notably, there were substantial differences between the high and low-risk groups in extensive respects, such as the activating transcription factors, cell pseudotime features, cell intercommunication patterns, immunocytes infiltration, chemotherapy response, and potential drug resistance. KIF26B was selected to construct a prognostic model from the identified 33 prognosis-related genes, and high expression of KIF26B predicted poorer prognosis in ovarian cancer. Ultimately, further in vitro experiments demonstrated that KIF26B participated in the proliferation and cisplatin resistance of OC cells. Knockdown of KIF26B increased the sensitivity of OC cells to elesclomol, a cuproptosis agonists. Conclusion This study constructed a new cuproptosis-related gene signature that has a good prognostic capacity in assessing the outcome of OC patients. This study enhances our understanding of cuproptosis associated with ovarian cancer aggressiveness, cross-talk with immunocytes, and serves as a novel chemotherapy strategy.

Abstract As a highly aggressive bone malignancy, osteosarcoma poses a significant therapeutic challenge, especially in the setting of metastasis or recurrence. This study aimed to investigate the potential of CD8-Tex cell-associated genes as prognostic biomarkers to reveal the immunogenomic profile of osteosarcoma and guide therapeutic decisions. mRNA expression data and clinical details of osteosarcoma patients were obtained from the TCGA database (TARGET-OS dataset). The GSE21257 dataset (from the GEO database) was used as an external validation set to provide additional information on osteosarcoma specimens. 84 samples from the TARGET-OS dataset were used as the training set, and 53 samples from the GSE21257 dataset served as the external validation cohort. Univariate Cox regression analysis was utilized to identify CD8 Tex cell genes associated with prognosis. The LASSO algorithm was performed for 1000 iterations to select the best subset to form the CD8 Tex cell gene signature (TRS). Final genes were identified using the multivariate Cox regression model of the LASSO algorithm. Risk scores were calculated to categorize patients into high- and low-risk groups, and clinical differences were explored by Kaplan–Meier survival analysis to assess model performance. Prediction maps were constructed to estimate 1-, 3-, and 5 year survival rates for osteosarcoma patients, including risk scores for CD8 Texcell gene markers and clinicopathologic factors. The ssGSEA algorithm was used to assess the differences in immune function between TRS-defined high- and low-risk groups. TME and immune cell infiltration were further assessed using the ESTIMATE and CIBERSORT algorithms. To explore the relationship between immune checkpoint gene expression levels and the two risk-defined groups. A CD8 Tex cell-associated gene signature was extracted from the TISCH database and prognostic markers including two genes were developed. The high-risk group showed lower survival, and model performance was validated by ROC curves and C-index. Predictive plots were constructed to demonstrate survival estimates, combining CD8 Tex cell gene markers and clinical factors. This study provides valuable insights into the molecular and immune characteristics of osteosarcoma and offers potential avenues for advances in therapeutic approaches.

Abstract Background The Warburg effect is a hallmark characteristic of colorectal cancer (CRC). Despite extensive research, the role of long non-coding RNAs (lncRNAs) in influencing the Warburg effect remains incompletely understood. Our study aims to identify lncRNAs that may modulate the Warburg effect by functioning as competing endogenous RNAs (ceRNAs). Methods Utilizing bioinformatics approaches, we extracted glycolysis-associated gene data from the Kyoto Encyclopedia of Genes and Genomes (KEGG) and identified 101 glycolysis-related lncRNAs in CRC. We employed Univariable Cox regression, Least Absolute Shrinkage and Selection Operator (LASSO) regression analysis, and Multivariable Cox regression to develop a prognostic model comprising four glycolysis-linked lncRNAs. We then constructed a prognostic nomogram integrating this lncRNA model with other relevant clinical parameters. Results The prognostic efficacy of our four-lncRNA signature and its associated nomogram was validated in both training and validation cohorts. Functional assays demonstrated significant glycolysis and hexokinase II (HK2) inhibition following the silencing of RUNDC3A − AS1, a key lncRNA in our prognostic signature, highlighting its regulatory importance in the Warburg effect. Conclusions Our research illuminates the critical role of glycolysis-centric lncRNAs in CRC. The developed prognostic model and nomogram underscore the pivotal prognostic and regulatory significance of the lncRNA RUNDC3A − AS1 in the Warburg effect in colorectal cancer.

Abstract Tumor-associated macrophages (TAMs) are a specific subset of macrophages that reside inside the tumor microenvironment. The dynamic interplay between TAMs and tumor cells plays a crucial role in the treatment response and prognosis of lung adenocarcinoma (LUAD). The study aimed to examine the association between TAMs and LUAD to advance the development of targeted strategies and immunotherapeutic approaches for treating this type of lung cancer. The study employed single-cell mRNA sequencing data to characterize the immune cell composition of LUAD and delineate distinct subpopulations of TAMs. The “BayesPrism” and “Seurat” R packages were employed to examine the association between these subgroups and immunotherapy and clinical features to identify novel immunotherapy biomarkers. Furthermore, a predictive signature was generated to forecast patient prognosis by examining the gene expression profile of immunotherapy-associated TAMs subsets and using 104 machine-learning techniques. A comprehensive investigation has shown the existence of a hitherto unidentified subgroup of TAMs known as RGS1 + TAMs, which has been found to have a strong correlation with the efficacy of immunotherapy and the occurrence of tumor metastasis in LUAD patients. CD83 was identified CD83 as a distinct biomarker for the expression of RGS1 + TAMs, showcasing its potential utility as an indicator for immunotherapeutic interventions. Furthermore, the prognostic capacity of the RTMscore signature, encompassing three specific mRNA (NR4A2, MMP14, and NPC2), demonstrated enhanced robustness when contrasted against the comprehensive collection of 104 features outlined in the published study. CD83 has potential as an immunotherapeutic biomarker. Meanwhile, The RTMscore signature established in the present study might be beneficial for survival prognostication.

Abstract Background The interplay between exosomes and the tumor microenvironment (TME) remains unclear. We investigated the influence of exosomes on the TME in hepatocellular carcinoma (HCC), focusing on their mRNA expression profile. Methods mRNA expression profiles of exosomes were obtained from exoRBase. RNA sequencing data from HCC patients’ tumors were acquired from The Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC). An exosome mRNA-related risk score model of prognostic value was established. The patients in the two databases were divided into high- and low-risk groups based on the median risk score value, and used to validate one another. Functional enrichment analysis was performed based on a differential gene prognosis model (DGPM). CIBERSORT was used to assess the abundance of immune cells in the TME. The correlation between the expression levels of immune checkpoint-related genes and DGPM was analyzed alongside the prediction value to drug sensitivity. Results A prognostic exosome mRNA-related 4-gene signature (DYNC1H1, PRKDC, CCDC88A, and ADAMTS5) was constructed and validated. A prognostic nomogram had prognostic ability for HCC. The genes for this model are involved in extracellular matrix, extracellular matrix (ECM)-receptor interaction, and the PI3K-Akt signaling pathway. Expression of genes here had a positive correlation with immune cell infiltration in the TME. Conclusions Our study results demonstrate that an exosome mRNA-related risk model can be established in HCC, highlighting the functional significance of the molecules in prognosis and risk stratification.

Abstract Background Dysregulation of iron metabolism has been shown to have significant implications for cancer development. We aimed to investigate the prognostic and immunological significance of iron metabolism-related genes (IMRGs) in nasopharyngeal carcinoma (NPC). Methods Multiple Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) datasets were analyzed to identify key IMRGs associated with prognosis. Additionally, the immunological significance of IMRGs was explored. Results A novel risk model was established using the LASSO regression algorithm, incorporating three genes (TFRC, SLC39A14, and ATP6V0D1).This model categorized patients into low and high-risk groups, and Kaplan–Meier analysis revealed significantly shorter progression-free survival for the high-risk group (P

Abstract Natural killer (NK) cells play essential roles in the tumor development, diagnosis, and prognosis of tumors. In this study, we aimed to establish a reliable signature based on marker genes in NK cells, thus providing a new perspective for assessing immunotherapy and the prognosis of patients with gastric cancer (GC). We analyzed a total of 1560 samples retrieved from the public database. We performed a comprehensive analysis of single-cell RNA-sequencing (scRNA-seq) data of gastric cancer and identified 377 marker genes for NK cells. By performing Cox regression analysis, we established a 12-gene NK cell-associated signature (NKCAS) for the Cancer Genome Atlas (TCGA) cohort, that assigned GC patients into a low-risk group (LRG) or a high-risk group (HRG). In the TCGA cohort, the areas under curve (AUC) value were 0.73, 0.81, and 0.80 at 1, 3, and 5 years. External validation of the predictive ability for the signature was then validated in the Gene Expression Omnibus (GEO) cohorts (GSE84437). The expression levels of signature genes were measured and validated in GC cell lines by real-time PCR. Moreover, NKCAS was identified as an independent prognostic factor by multivariate analysis. We combined this with a variety of clinicopathological characteristics (age, M stage, and tumor grade) to construct a nomogram to predict the survival outcomes of patients. Moreover, the LRG showed higher immune cell infiltration, especially CD8+ T cells and NK cells. The risk score was negatively associated with inflammatory activities. Importantly, analysis of the independent immunotherapy cohort showed that the LRG had a better prognosis and immunotherapy response when compared with the HRG. The identification of NK cell marker genes in this study suggests potential therapeutic targets. Additionally, the developed predictive signatures and nomograms may aid in the clinical management of GC.

Abstract Background GBM, also known as glioblastoma multiforme, is the most prevalent and lethal type of brain cancer. The cell proliferation, invasion, angiogenesis, and treatment of gliomas are significantly influenced by oxidative stress. Nevertheless, the connection between ORGs and GBM remains poorly comprehended. The objective of this research is to investigate the predictive significance of ORGs in GBM and their potential as targets for therapy. Methods We identified differentially expressed genes in glioma and ORGs from public databases. A risk model was established using LASSO regression and Cox analysis, and its performance was evaluated with ROC curves. We then performed consistent cluster analysis on the model, examining its correlation with immunity and drug response. Additionally, PCR, WB and IHC were employed to validate key genes within the prognostic model. Results 9 ORGs (H6PD, BMP2, SPP1, HADHA, SLC25A20, TXNIP, ACTA1, CCND1, EEF1A1) were selected via differential expression analysis, LASSO and Cox analysis, and incorporated into the risk model with high predictive accuracy. Enrichment analyses using GSVA and GSEA focused predominantly on malignancy-associated pathways. Subtype C of GBM had the best prognosis with the lowest risk score. Furthermore, the model exhibited a strong correlation with the infiltration of immune cells and had the capability to pinpoint potential targeted therapeutic medications for GBM. Ultimately, we selected HADHA for in vitro validation. The findings indicated that GBM exhibits a significant upregulation of HADHA. Knockdown of HADHA inhibited glioma cell proliferation and diminished their migration and invasion capacities and influenced the tumor growth in vivo. Conclusion The risk model, built upon 9 ORGs and the identification of GBM subtypes, suggests that ORGs have a broad application prospect in the clinical immunotherapy and targeted drug treatment of GBM. HADHA significantly influences the development of gliomas, both in vivo and in vitro.

Abstract Background Bladder cancer (BCa) stands out as a prevalent and highly lethal malignancy worldwide. Chemoresistance significantly contributes to cancer recurrence and progression. Traditional Tumor Node Metastasis (TNM) stage and molecular subtypes often fail to promptly identify treatment preferences based on sensitivity. Methods In this study, we developed a prognostic signature for BCa with uni-Cox + LASSO + multi-Cox survival analysis in multiple independent cohorts. Six machine learning algorithms were adopted to screen out the hub gene, RAC3. IHC staining was used to validate the expression of RAC3 in BCa tumor tissue. RT-qPCR and Western blot were performed to detect and quantify the mRNA and protein levels of RAC3. CCK8, colony formation, wound healing, and flow cytometry analysis of apoptosis were employed to determine cell proliferation, migration, and apoptosis. Molecular docking was used to find small target drugs, PIK-75. 3D cell viability assay was applied to evaluate the ATP viability of bladder cancer organoids before and after PIK-75 treated. Results The established clinical prognostic model, GIRS, comprises 13 genes associated with gemcitabine resistance and immunology. This model has demonstrated robust predictive capabilities for survival outcomes across various independent public cohorts. Additionally, the GIRS signature shows significant correlations with responses to both immunotherapy and chemotherapy. Leveraging machine learning algorithms, the hub gene, RAC3, was identified, and potential upstream transcription factors were screened through database analysis. IHC results showed that RAC3 was higher expressed in GEM-resistant BCa patients. Employing molecular docking, the small molecule drug PIK-75, as binding to RAC3, was identified. Experiments on cell lines, organoids and animals validated the biological effects of PIK-75 in bladder cancer. Conclusions The GIRS signature offers a valuable complement to the conventional anatomic TNM staging system and molecular subtype stratification in bladder cancer. The hub gene, RAC3, plays a crucial role in BCa and is significantly associated with resistance to gemcitabine. The small molecular drug, PIK-75 having the potential as a therapeutic agent in the context of gemcitabine-resistant and immune-related pathways.

Abstract Disulfidptosis, a new type of regulated cell death associated with the actin cytoskeleton, provides a new therapeutic tool for cancers. The direct relationship between disulfidptosis-related lncRNAs(DRLs) in liver hepatocellular carcinoma(HCC) remains unclear. We acquired transcriptomic data, corresponding clinical data, and tumor mutation data of HCC from the TCGA database. First of all, DRLs were determined through correlation analysis. Then, a prognostic model containing six DRLs was created by adopting univariate Cox regression, LASSO algorithm and multivariate Cox regression analysis. Based on the model, 424 HCC patients were divided into high- and low-risk groups. Next, we structured ROC curves and PCA through combining the model and clinical data. Enrichment analysis and immune infiltration analysis were adopted to further explore the relationship between the model and prognosis. In addition, we explored the relationship between the model and tumor mutation burden (TMB). There were significant differences between high- and low- risk groups, and patients in the high-risk group showed poor prognosis. Enrichment analysis suggested that metabolic progress was obviously different between the two groups. According to the analysis of immune infiltration, there were several differences in immune cells, function, and checkpoints. Patients with high-risk and high TMB demonstrated the least favorable prognosis. The two risk groups both manifested visiblly in chemotherapy drug sensitivity. To sum up, we set up a DRL-based signature and that may provide a predictable value for the prognosis and use of chemotherapy drugs for HCC patients.

Abstract Background Gastric cancer (GC) remains a malignant tumor with high morbidity and mortality, accounting for approximately 1,080,000 diagnosed cases and 770,000 deaths worldwide annually. Disulfidptosis, characterized by the stress-induced abnormal accumulation of disulfide, is a recently identified form of programmed cell death. Substantial studies have demonstrated the significant influence of immune clearance on tumor progression. Therefore, we aimed to explore the intrinsic correlations between disulfidptosis and immune-related genes (IRGs) in GC, as well as the potential value of disulfidptosis-related immune genes (DRIGs) as biomarkers. Methods This study incorporated the single-cell RNA sequencing (scRNA-seq) dataset GSE183904 and transcriptome RNA sequencing of GC from the TCGA database. Disulfidptosis-related genes (DRGs) and IRGs were derived from the representative literature on both cell disulfidptosis and immunity. The expression and distribution of DRGs were investigated at the single-cell level in different GC cell types. Pearson correlation analysis was used to identify the IRGs closely related to disulfidptosis. The prognostic signature of DRIGs was established using Cox and LASSO analyses. We then analyzed and evaluated the differences in long-term prognosis, Gene Set Enrichment Analysis (GSEA), immune infiltration, mutation profile, CD274 expression, and response to chemotherapeutic drugs between the two groups. A tissue array containing 63 paired GC specimens was used to verify the expression of 4 DRIGs and disulfidptosis regulator SLC7A11 through immunohistochemistry staining. Results The scRNA-seq analysis found that SLC7A11, SLC3A2, RPN1 and NCKAP1 were enriched in specific cell types and closely related to immune infiltration. Four DIRGs (GLA, HIF-1α, VPS35 and CDC37) were successfully identified to establish a signature to potently predict the survival time of GC patients. Patients with high risk scores generally experienced worse prognoses and exhibited greater resistant to classical chemotherapy drugs. Furthermore, the expression of GLA, HIF-1α, VPS35, CDC37 and SLC7A11 were elevated in GC tissues. A high expression of GLA, HIF-1α, VPS35 or CDC37 was associated with more advanced clinical stage of GC and increased SLC7A11 expression. Conclusion Current study first highlights the potential value of DRIGs as biomarkers in GC. We successfully constructed a robust model incorporating four DRIGs to accurately predict the survival time and clinicopathological characteristics of GC patients.

Abstract Background Numerous studies have shown that Schistosoma japonicum infection correlates with an increased risk of liver hepatocellular carcinoma (LIHC). However, data regarding the role of this infection in LIHC oncogenesis are scarce. This study aimed to investigate the potential mechanisms of hepatocarcinogenesis associated with Schistosoma japonicum infection. Methods By examining chronic liver disease as a mediator, we identified the genes contributing to Schistosoma japonicum infection and LIHC. We selected 15 key differentially expressed genes (DEGs) using weighted gene co-expression network analysis (WGCNA) and random survival forest models. Consensus clustering revealed two subgroups with distinct prognoses. Least Absolute Shrinkage and Selection Operator (LASSO) and Cox regression identified six prognostic DEGs, forming an Schistosoma japonicum infection-associated signature for strong prognosis prediction. This signature, which is an independent LIHC risk factor, was significantly correlated with clinical variables. Four DEGs, including BMI1, were selected based on their protein expression levels in cancerous and normal tissues. We confirmed BMI1's role in LIHC using Schistosoma japonicum-infected mouse models and molecular experiments. Results We identified a series of DEGs that mediate schistosomiasis, the parasitic disease caused by Schistosoma japonicum infection, and hepatocarcinogenesis, and constructed a suitable prognostic model. We analyzed the mechanisms by which these DEGs regulate disease and present the differences in prognosis between the different genotypes. Finally, we verified our findings using molecular biology experiments. Conclusion Bioinformatics and molecular biology analyses confirmed a relationship between schistosomiasis and liver hepatocellular cancer. Furthermore, we validated the role of a potential oncoprotein factor that may be associated with infection and carcinogenesis. These findings enhance our understanding of Schistosoma japonicum infection's role in LIHC carcinogenesis.

Abstract CD8+ exhausted T cells (CD8+ Tex) played a vital role in the progression and therapeutic response of cancer. However, few studies have fully clarified the characters of CD8+ Tex related genes in ovarian cancer (OC). The CD8+ Tex related prognostic signature (TRPS) was constructed with integrative machine learning procedure including 10 methods using TCGA, GSE14764, GSE26193, GSE26712, GSE63885 and GSE140082 dataset. Several immunotherapy benefits indicators, including Tumor Immune Dysfunction and Exclusion (TIDE) score, immunophenoscore (IPS), TMB score and tumor escape score, were used to explore performance of TRPS in predicting immunotherapy benefits of OC. The TRPS constructed by Enet (alpha = 0.3) method acted as an independent risk factor for OC and showed stable and powerful performance in predicting clinical outcome of patients. The C-index of the TRPS was higher than that of tumor grade, clinical stage, and many developed signatures. Low TRPS score indicated a higher level of CD8+ T cell, B cell, macrophage M1, and NK cells, representing a relative immunoactivated ecosystem in OC. OC patients with low risk score had a higher PD1&CTLA4 immunophenoscore, higher TMB score, lower TIDE score and lower tumor escape score, suggesting a better immunotherapy response. Moreover, higher TRPS score indicated a higher score of cancer-related hallmarks, including angiogenesis, EMT, hypoxia, glycolysis, and notch signaling. Vitro experiment showed that ARL6IP5 was downregulated in OC tissues and inhibited tumor cell proliferation. The current study constructed a novel TRPS for OC, which could serve as an indicator for predicting the prognosis, immune infiltration and immunotherapy benefits for OC patients.

Abstract 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.

Abstract Background Cholangiocarcinoma represents a malignant neoplasm originating from the hepatobiliary tree, with a subset of tumors developing inside the liver. Intrahepatic cholangiocarcinomas (ICC) commonly exhibit an asymptomatic presentation, rendering both diagnosis and treatment challenging. Cuproptosis, an emerging regulated cell death pathway induced by copper ions, has garnered attention recently. As cancer cells show altered copper metabolism and comparatively higher copper needs, cuproptosis may play a role in the development of ICC. However, studies investigating this possibility are currently lacking. Methods Single-cell and bulk RNA sequence data were analyzed, and correlations were established between the expression of cuproptosis-related molecules and ICC patient survival. Genes with predicting survival were used to create a CUPT score using Cox and LASSO regression and tumor mutation burden (TMB) analysis. The CIBERSORT software was employed to characterize immune cell infiltration within the tumors. Furthermore, immune infiltration prediction, biological function enrichment, and drug sensitivity analyses were conducted to explore the potential implications of the cuproptosis-related signature. The effects of silencing solute carrier family 39 member 4 gene (SLC39A4) expression using siRNA were investigated using assays measuring cell proliferation, colony formation, and cell migration. Key genes of cuproptosis were detected by western blotting. Results The developed CUPT score divided patients into high and low CUPT score groups. Those with a low score had significantly better prognosis and longer survival. In contrast, high CUPT scores were associated with worse clinical outcomes and significantly higher TMB. Comparisons of the two groups also indicated differences in the immune infiltrate present in the tumors. Finally, we were able to identify 95 drugs potentially affecting the cuproptosis pathway. Some of these might be effective in the treatment of ICC. The in vitro experiments revealed that suppressing the expression of SLC39A4 in ICC cell lines resulted in reduced cell proliferation, colony formation, and cell migration. It also led to an increase in cell death and the upregulation of key genes associated with cuproptosis, namely ferredoxin 1 (FDX1) and dihydrolipoyl transacetylase (DLAT). These findings strongly suggest that this cuproptosis-associated molecule may play a pivotal role in the development and metastasis of ICC. Conclusions Changes in the expression of a cuproptosis-related gene signature can predict the clinical prognosis of ICC with considerable accuracy. This supports the notion that cuproptosis influences the diversity and complexity of the immune microenvironment, mutational landscape, and biological behavior of ICC. Understanding this pathway better may hold promise for the development of innovative strategies in the management of this disease.

Background: N6-methylandenosine-related long non-coding RNAs (m6A-related lncRNAs) play a crucial role in the cancer progression and immunotherapeutic efficacy. The potential function of m6A-related lncRNAs signature in cervical cancer has not been systematically clarified. Methods: RNA-seq and the clinical data of cervical cancer were extracted from The Cancer Genome Atlas. All of the patients were randomly classified into training and testing cohorts. The m6A-related lncRNAs prognostic model was constructed by LASSO regression using data in the training cohort.The predictive value of the signature was validated in the whole cohort and testing cohort. Cervical cancer patients were divided into low- and high-risk subgroups by the median value of risk scores. Kaplan-Meier analysis, principal-component analysis (PCA), functional enrichment annotation, and nomogram were used for further evaluation. We also examined the immune response and potential drug sensitivity targeting this model. Results: Seventy-nine prognostic m6A-related lncRNAs were screened. The risk model comprising four m6A-related lncRNAs (AL139035.1, AC015922.2, AC073529.1, AC008124.1) was identified and verified as an independent prognostic predictor of cervical cancer. A nomogram based on age, tumor grade, clinical stage, TNM stage, and four m6A-related lncRNAs risk signatures was generated. It displayed good accuracy and reliability in predicting the overall survival of patients with CC. Based on our risk model, cervical cancer patients with potential immunotherapy benefits from the candidate drugs could be effectively screened. Conclusion: The four m6A-related lncRNAs signature may provide new targets and allow the prediction of immunotherapy response, which can assist developing individualized treatment for cervical cancer.

The prognosis for patients with colorectal cancer (CRC) remains worse than expected due to metastasis, recurrence, and resistance to chemotherapy. Colorectal cancer stem cells (CRCSCs) play a vital role in tumor metastasis, recurrence, and chemotherapy resistance. However, there are currently no prognostic markers based on CRCSCs-related genes available for clinical use. In this study, single-cell transcriptome sequencing was employed to distinguish cancer stem cells (CSCs) in the CRC microenvironment and analyze their properties at the single-cell level. Subsequently, data from TCGA and GEO databases were utilized to develop a prognostic risk model for CRCSCs-related genes and validate its diagnostic performance. Additionally, functional enrichment, immune response, and chemotherapeutic drug sensitivity of the relevant genes in the risk model were investigated. Lastly, the key gene RPS17 in the risk model was identified as a potential prognostic marker and therapeutic target for further comprehensive studies. Our findings provide new insights into the prognostic treatment of CRC and offer novel perspectives for a systematic and comprehensive understanding of CRC development. [ABSTRACT FROM AUTHOR]

Background: As a highly heterogeneous tumor, non-small cell lung cancer (NSCLC) is famous for its high incidence and mortality worldwide. Smoking can cause genetic changes, which leading to the occurrence and progress of NSCLC. Nevertheless, the function of smoking-related genes in NSCLC needs more research. Methods: We downloaded transcriptome data and clinicopathological parameters from Gene Expression Omnibus (GEO) databases, and screened smoking-related genes. Lasso regression were applied to establish the 7-gene signature. The associations between the 7-gene signature and immune microenvironment analysis, survival analysis, drug sensitivity analysis and enriched molecular pathways were studied. Ultimately, cell function experiments were conducted to research the function of FCGBP in NSCLC. Results: Through 7-gene signature, NSCLC samples were classified into high-risk group (HRG) and low-risk group (LRG). Significant difference in overall survival (OS) between HRG and LRG was found. Nomograms and ROC curves indicated that the 7-gene signature has a stable ability in predicting prognosis. Through the analysis of immune microenvironment, we found that LRG patients had better tumor immune activation. FCGBP showed the highest mutation frequency among the seven prognostic smoking related genes (LRRC31, HPGD, FCGBP, SPINK5, CYP24A1, S100P and FGG), and was notable down-regulated in NSCLC smokers compared with non-smoking NSCLC patients. The cell experiments confirmed that FCGBP knockdown promoting proliferation, migration, and invasion in NSCLC cells. Conclusion: This smoking-related prognostic signature represents a promising tool for assessing prognosis and tumor microenvironment in smokers with NSCLC. The role of FCGBP in NSCLC was found by cell experiments, which can be served as diagnostic biomarker and immunotherapy target for NSCLC. [ABSTRACT FROM AUTHOR]

Objectives: This study aimed to analyze the crosstalk network among T cells, epithelial cells, and fibroblasts in the tumor microenvironment of oral squamous cell carcinoma (OSCC) and to determine their prognostic values. Materials and Methods: Single‐cell subpopulation identification and communication analysis identified crosstalk markers. The least absolute shrinkage and selection operator Cox analysis identified key prognostic features by integrating the bulk transcriptome and clinical parameters. Functional analysis and immune infiltration were explored to determine possible mechanisms. Results: Interactions between epithelial cells and fibroblasts primarily involve MIF, MK, PTN, IGF, EGF, and PERIOSTIN, whereas T cells interact with epithelial cells and fibroblasts through MIF, CXCL, PAR, IFN, and EGF signals. We constructed a novel prognostic feature comprising 13 crosstalk genes: HBEGF, FGF7, GRN, ITGB5, CXCR6, ERBB2, AREG, F2RL2, NAMPT, KLK12, HMGB2, TUBA1B, and KLRD1. Patients were stratified based on the RiskScore. Functional analysis revealed that the high‐risk group was enriched in immunosuppressive pathways (p < 0.001). Immune checkpoints including PD‐1, PD‐L1, and CTLA4 were more highly expressed in the high‐risk group (p < 0.05). Conclusions: The crosstalk network among T cells, epithelial cells, and fibroblasts is complex and may have implications for prognosis and clinical treatments of OSCC patients. [ABSTRACT FROM AUTHOR]

Esophageal cancer (ESCA) has a bad prognosis. Long non-coding RNA (lncRNA) impacts on cell proliferation. However, the prognosis function of N6-methyladenosine (m6A)-associated lncRNAs (m6A-lncRNAs) in ESCA remains unknown. Univariate Cox analysis was applied to investigate prognosis related m6A-lncRNAs, based on which the samples were clustered. Wilcoxon rank and Chi-square tests were adopted to compare the clinical traits, survival, pathway activity and immune infiltration in different clusters where overall survival, clinical traits (N stage), tumor-invasive immune cells and pathway activity were found significantly different. Through least absolute shrinkage and selection operator and proportional hazard (Lasso-Cox) model, five m6A-lncRNAs were selected to construct the prognostic signature (m6A-lncSig) and risk score. To investigate the link between risk score and clinical traits or immunological microenvironments, Chi-square test and Spearman correlation analysis were utilized. Risk score was found connected with N stage, tumor stage, different clusters, macrophages M2, B cells naive and T cells CD4 memory resting. Risk score and tumor stage were found as independent prognostic variables. And the constructed nomogram model had high accuracy in predicting prognosis. The obtained m6A-lncSig could be taken as potential prognostic biomarker for ESCA patients. This study offers a theoretical foundation for clinical diagnosis and prognosis of ESCA. [ABSTRACT FROM AUTHOR]

Stage II colon cancer (CC) encompasses a heterogeneous group of patients with diverse survival experiences: 87% to 58% 5-year relative survival rates for stages IIA and IIC, respectively. While stage IIA patients are usually spared the adjuvant chemotherapy, some of them relapse and may benefit from it; thus, their timely identification is crucial. Current gene expression signatures did not specifically target this group nor did they find their place in clinical practice. Since processes at invasion front have also been linked to tumor progression, we hypothesize that aside from bulk tumor features, focusing on the invasion front may provide additional clues for this stratification. A retrospective matched case-control collection of 39 stage IIA microsatellite-stable (MSS) untreated CCs was analyzed to identify prognostic gene expression-based signatures. The endpoint was defined as relapse within 5 years vs. no relapse for at least 6 years. From the same tumors, three different classifiers (bulk tumor, invasion front, and constrained baseline on bulk tumor) were developed and their performance estimated. The baseline classifier, while the weakest, was validated in two independent data sets. The best performing signature was based on invasion front profiles [area under the receiver operating curve (AUC) = 0.931 (0.815--1.0)] and contained genes associated with KRAS pathway activation, apical junction complex, and heme metabolism. Its combination with bulk tumor classifier further improved the accuracy of the predictions. [ABSTRACT FROM AUTHOR]

Background: Colorectal cancer (CRC) presents a significant global health burden, characterized by a heterogeneous molecular landscape and various genetic and epigenetic alterations. Programmed cell death (PCD) plays a critical role in CRC, offering potential targets for therapy by regulating cell elimination processes that can suppress tumor growth or trigger cancer cell resistance. Understanding the complex interplay between PCD mechanisms and CRC pathogenesis is crucial. This study aims to construct a PCD‐related prognostic signature in CRC using machine learning integration, enhancing the precision of CRC prognosis prediction. Method: We retrieved expression data and clinical information from the Cancer Genome Atlas and Gene Expression Omnibus (GEO) datasets. Fifteen forms of PCD were identified, and corresponding gene sets were compiled. Machine learning algorithms, including Lasso, Ridge, Enet, StepCox, survivalSVM, CoxBoost, SuperPC, plsRcox, random survival forest (RSF), and gradient boosting machine, were integrated for model construction. The models were validated using six GEO datasets, and the programmed cell death score (PCDS) was established. Further, the model's effectiveness was compared with 109 transcriptome‐based CRC prognostic models. Result: Our integrated model successfully identified differentially expressed PCD‐related genes and stratified CRC samples into four subtypes with distinct prognostic implications. The optimal combination of machine learning models, RSF + Ridge, showed superior performance compared with traditional methods. The PCDS effectively stratified patients into high‐risk and low‐risk groups, with significant survival differences. Further analysis revealed the prognostic relevance of immune cell types and pathways associated with CRC subtypes. The model also identified hub genes and drug sensitivities relevant to CRC prognosis. Conclusion: The current study highlights the potential of integrating machine learning models to enhance the prediction of CRC prognosis. The developed prognostic signature, which is related to PCD, holds promise for personalized and effective therapeutic interventions in CRC. [ABSTRACT FROM AUTHOR]

Objective: Significant advancements have been made in hepatocellular carcinoma (HCC) therapeutics, such as immunotherapy for treating patients with HCC. However, there is a lack of reliable biomarkers for predicting the response of patients to therapy, which continues to be challenging. Cancer stem cells (CSCs) are involved in the oncogenesis, drug resistance, and invasion, as well as metastasis of HCC cells. Therefore, in this study, we aimed to create an mRNA expression-based stemness index (mRNAsi) model to predict the response of patients with HCC to immunotherapy. Methods: We retrieved gene expression and clinical data of patients with HCC from the GSE14520 dataset and the Cancer Genome Atlas (TCGA) database. Next, we used the "one-class logistic regression (OCLR)" algorithm to obtain the mRNAsi of patients with HCC. We performed "unsupervised consensus clustering" to classify patients with HCC based on the mRNAsi scores and stemness subtypes. The relationships between the mRNAsi model, clinicopathological features, and genetic profiles of patients were compared using various bioinformatic methods. We screened for differentially expressed genes to establish a stemness-based classifier for predicting the patient's prognosis. Next, we determined the effect of risk scores on the tumor immune microenvironment (TIME) and the response of patients to immune checkpoint blockade (ICB). Finally, we used qRT-PCR to investigate gene expression in patients with HCC. Results: We screened CSC-related genes using various bioinformatics tools in patients from the TCGA-LIHC cohort. We constructed a stemness classifier based on a nine-gene (PPARGC1A, FTCD, CFHR3, MAGEA6, CXCL8, CABYR, EPO, HMMR, and UCK2) signature for predicting the patient's prognosis and response to ICBs. Further, the model was validated in an independent GSE14520 dataset and performed well. Our model could predict the status of TIME, immunogenomic expressions, congenic pathway, and response to chemotherapy drugs. Furthermore, a significant increase in the proportion of infiltrating macrophages, Treg cells, and immune checkpoints was observed in patients in the high-risk group. In addition, tumor cells in patients with high mRNAsi scores could escape immune surveillance. Finally, we observed that the constructed model had a good expression in the clinical samples. The HCC tumor size and UCK2 genes expression were significantly alleviated and decreased, respectively, by treatments of anti-PD1 antibody. We also found knockdown UCK2 changed expressions of immune genes in HCC cell lines. Conclusion: The novel stemness-related model could predict the prognosis of patients and aid in creating personalized immuno- and targeted therapy for patients in HCC. [ABSTRACT FROM AUTHOR]

Purpose: Cuproptosis, programmed cell death by intracellular copper-mediated lipoylated protein aggregation, is involved in various tumorigenesis and drug resistance abilities by mediating the tumor microenvironment. Previous studies have demonstrated that serum copper levels are higher in OC patients than in normal subjects. However, the exact relationship between cuproptosis and ovarian cancer progression remains to be further elucidated. Methods: The Cancer Genome Atlas (TCGA) and gene expression omnibus (GEO) datasets were utilized to establish a cuproptosis-related prognostic signature in ovarian cancer. Subsequently, the bulk RNA-seq analysis and single-cell RNA-seq analysis were used to identify the relationship between signature with immune cell infiltration, chemotherapy, and cuproptosis-related scoring (CuRS) system. Finally, the potential biological functional roles of target genes in cuproptosis were validated in vitro. Results: By using LASSO-Cox regression analysis to establish the cuproptosis-related prognostic model, our works demonstrated the accuracy and efficiency of our model in the TCGA (583 OC patients) and GEO (260 OC patients) OC cohorts, and the high-scoring groups showed worse survival outcomes. Notably, there were substantial differences between the high and low-risk groups in extensive respects, such as the activating transcription factors, cell pseudotime features, cell intercommunication patterns, immunocytes infiltration, chemotherapy response, and potential drug resistance. KIF26B was selected to construct a prognostic model from the identified 33 prognosis-related genes, and high expression of KIF26B predicted poorer prognosis in ovarian cancer. Ultimately, further in vitro experiments demonstrated that KIF26B participated in the proliferation and cisplatin resistance of OC cells. Knockdown of KIF26B increased the sensitivity of OC cells to elesclomol, a cuproptosis agonists. Conclusion: This study constructed a new cuproptosis-related gene signature that has a good prognostic capacity in assessing the outcome of OC patients. This study enhances our understanding of cuproptosis associated with ovarian cancer aggressiveness, cross-talk with immunocytes, and serves as a novel chemotherapy strategy. [ABSTRACT FROM AUTHOR]

Although anoikis plays a role in cancer metastasis and aggressiveness, it has rarely been reported in diffuse large B cell lymphoma (DLBCL). We obtained RNA sequencing data and matched clinical data from the GEO database. An anoikis-related genes (ARGs)-based risk signature was developed in GSE10846 training cohort and validated in three other cohorts. Additionally, we predicted half-maximal inhibitory concentration (IC50) of drugs based on bioinformatics method and obtained the actual IC50 to some chemotherapy drugs via cytotoxicity assay. The high-risk group, as determined by our signature, was associated with worse prognosis and an immunosuppressive environment in DLBCL. Meanwhile, the nomogram based on eight variables had more accurate ability in forecasting the prognosis than the international prognostic index in DLBCL. The prediction of IC50 indicated that DLBCL patients in the high-risk group were more sensitive to doxorubicin, IPA-3, lenalidomide, gemcitabine, and CEP.701, while patients in the low-risk group were sensitive to cisplatin and dasatinib. Consistent with the prediction, cytotoxicity assay suggested the higher sensitivity to doxorubicin and gemcitabine and the lower sensitivity to dasatinib in the high-risk group in DLBCL. The ARG-based signature may provide a promising direction for prognosis prediction and treatment optimization for DLBCL patients. [ABSTRACT FROM AUTHOR]

Background: Anoikis is a cell death programmed to eliminate dysfunctional or damaged cells induced by detachment from the extracellular matrix. Utilizing an anoikis-based risk stratification is anticipated to understand melanoma's prognostic and immune landscapes comprehensively. Methods: Differential expression genes (DEGs) were analyzed between melanoma and normal skin tissues in The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression data sets. Next, least absolute shrinkage and selection operator, support vector machine–recursive feature elimination algorithm, and univariate and multivariate Cox analyses on the 308 DEGs were performed to build the prognostic signature in the TCGA–melanoma data set. Finally, the signature was validated in GSE65904 and GSE22155 data sets. NOTCH3, PIK3R2, and SOD2 were validated in our clinical samples by immunohistochemistry. Results: The prognostic model for melanoma patients was developed utilizing ten hub anoikis-related genes. The overall survival (OS) of patients in the high-risk subgroup, which was classified by the optimal cutoff value, was remarkably shorter in the TCGA–melanoma, GSE65904, and GSE22155 data sets. Low-risk patients exhibited low immune cell infiltration and high expression of immunophenoscores and immune checkpoints. They also demonstrated increased sensitivity to various drugs, including dasatinib and dabrafenib. NOTCH3, PIK3R2, and SOD2 were notably associated with OS by univariate Cox analysis in the GSE65904 data set. The clinical melanoma samples showed remarkably higher protein expressions of NOTCH3 (P = 0.003) and PIK3R2 (P = 0.009) than the para-melanoma samples, while the SOD2 protein expression remained unchanged. Conclusions: In this study, we successfully established a prognostic anoikis-connected signature using machine learning. This model may aid in evaluating patient prognosis, clinical characteristics, and immune treatment modalities for melanoma. [ABSTRACT FROM AUTHOR]

Lung adenocarcinoma (LUAD) is a major subtype of non‐small‐cell lung cancer and accompanies high mortality rates. While the role of bilirubin metabolism in cancer is recognized, its specific impact on LUAD and patient response to immunotherapy needs to be elucidated. This study aimed to develop a prognostic signature of bilirubin metabolism‐associated genes (BMAGs) to predict outcomes and efficacy of immunotherapy in LUAD. We analysed gene expression data from The Cancer Genome Atlas (TCGA) to identify survival‐related BMAGs and construct a prognostic model in LUAD. The prognostic efficacy of our model was corroborated by employing TCGA‐LUAD and five Gene Expression Omnibus datasets, effectively stratifying patients into risk‐defined cohorts with marked disparities in survival. The BMAG signature was indeed an independent prognostic determinant, outperforming established clinical parameters. The low‐risk group exhibited a more favourable response to immunotherapy, highlighted by increased immune checkpoint expression and immune cell infiltration. Further, somatic mutation profiling differentiated the molecular landscapes of the risk categories. Our screening further identified potential drug candidates preferentially targeting the high‐risk group. Our analysis of critical BMAGs showed the tumour‐suppressive role of FBP1, highlighting its suppression in LUAD and its inhibitory effects on tumour proliferation, migration and invasion, in addition to its involvement in cell cycle and apoptosis regulation. These findings introduce a potent BMAG‐based prognostic indicator and offer valuable insights for prognostication and tailored immunotherapy in LUAD. [ABSTRACT FROM AUTHOR]

As a highly aggressive bone malignancy, osteosarcoma poses a significant therapeutic challenge, especially in the setting of metastasis or recurrence. This study aimed to investigate the potential of CD8-Tex cell-associated genes as prognostic biomarkers to reveal the immunogenomic profile of osteosarcoma and guide therapeutic decisions. mRNA expression data and clinical details of osteosarcoma patients were obtained from the TCGA database (TARGET-OS dataset). The GSE21257 dataset (from the GEO database) was used as an external validation set to provide additional information on osteosarcoma specimens. 84 samples from the TARGET-OS dataset were used as the training set, and 53 samples from the GSE21257 dataset served as the external validation cohort. Univariate Cox regression analysis was utilized to identify CD8 Tex cell genes associated with prognosis. The LASSO algorithm was performed for 1000 iterations to select the best subset to form the CD8 Tex cell gene signature (TRS). Final genes were identified using the multivariate Cox regression model of the LASSO algorithm. Risk scores were calculated to categorize patients into high- and low-risk groups, and clinical differences were explored by Kaplan–Meier survival analysis to assess model performance. Prediction maps were constructed to estimate 1-, 3-, and 5 year survival rates for osteosarcoma patients, including risk scores for CD8 Texcell gene markers and clinicopathologic factors. The ssGSEA algorithm was used to assess the differences in immune function between TRS-defined high- and low-risk groups. TME and immune cell infiltration were further assessed using the ESTIMATE and CIBERSORT algorithms. To explore the relationship between immune checkpoint gene expression levels and the two risk-defined groups. A CD8 Tex cell-associated gene signature was extracted from the TISCH database and prognostic markers including two genes were developed. The high-risk group showed lower survival, and model performance was validated by ROC curves and C-index. Predictive plots were constructed to demonstrate survival estimates, combining CD8 Tex cell gene markers and clinical factors. This study provides valuable insights into the molecular and immune characteristics of osteosarcoma and offers potential avenues for advances in therapeutic approaches. [ABSTRACT FROM AUTHOR]