A novel signature to guide osteosarcoma prognosis and immune microenvironment: Cuproptosis-related lncRNA.

Esophageal cancer (EC) ranks among the most prevalent malignancies globally and represents a significant and growing public health burden. This study aimed to construct a prognostic model leveraging anoikis-related genes (ARGs) to predict patient survival and elucidate the immunological microenvironment in EC. The findings are anticipated to enhance prognostic accuracy and inform therapeutic strategies, ultimately improving patient outcomes and treatment efficacy. A comprehensive analysis was conducted using 11 control samples and 159 EC samples obtained from The Cancer Genome Atlas (TCGA) database, alongside associated clinical features. A total of 794 ARGs were curated from GeneCards database. Functional enrichment analyses of EC-related differentially expressed ARGs were performed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). Prognostic differential ARGs associated with EC were identified through univariate Cox regression analysis, while LASSO regression was employed to minimize overfitting and construct a robust risk prognostic model. The EC cohort was stratified into training and testing groups for model development and verification. Model performance was evaluated through risk curves, survival curves, time-dependent receiver operating characteristic (ROC) curves, ROC curves for the riskscore and clinical features, and independent prognostic analysis. A nomogram with high predictive accuracy was also developed to estimate the prognosis of EC patients. To assess the impact of the risk prognosis model on the immune microenvironment of EC, analyses included tumor microenvironment analysis, single-sample gene set enrichment analysis (ssGSEA), immune cell infiltration correlation analysis, and differential analysis of immune checkpoint expression. Drug sensitivity profiling was conducted to identify potential therapeutic agents for EC. Finally, the expression of selected ARGs was validated at the mRNA level in EC cell lines using real-time quantitative PCR (RT-qPCR). The ARG-based risk prognostic model was constructed incorporating four high-risk ARGs (CDK1, IL17A, FOXC2, and OLFM3) and two low-risk ARGs (PIP5K1C and MAPK1). This model demonstrated strong predictive accuracy for the survival outcomes of EC patients. Immune correlation analyses revealed that the high-risk group exhibited significantly lower immunological scores compared to the low-risk group. Notably, immune cells such as macrophages and mast cells were markedly downregulated in the high-risk group. Additionally, key immunological functions, including APC co-inhibition, parainflammation, Type I IFN Response, and Type II IFN Response, were significantly suppressed in the high-risk group. Eight immune checkpoint-related genes (TNFRSF25, TNFRSF14, CD70, TNFSF15, TMIGD2, CD160, TNFSF18, and HHLA2) displayed distinct expression differences between high- and low-risk groups. The nomogram developed from this model demonstrated high efficacy in predicting EC patient prognosis. Furthermore, six potential therapeutic agents for EC were identified: BIRB.0796, Camptothecin, CHIR.99021, Methotrexate, PF.4708671, and Vorinostat. Finally, the mRNA expression levels of ARGs were validated using RT-qPCR in EC cell lines. Compared to normal esophageal epithelial cells (NE-2), CDK1 and MAPK1 were significantly upregulated in two EC cell lines (KYSE-30 and KYSE-180). This study provides valuable insights into the prognostic outcomes and immune microenvironment of EC through the analysis of ARGs. Furthermore, several potential therapeutic agents for EC were identified, offering promising avenues for treatment. These findings hold significant potential for enhancing the survival outcomes of EC patients and provide meaningful guidance for clinical decision-making in managing this malignancy.

Background: Colon adenocarcinoma (COAD), a malignant gastrointestinal tumor, has the characteristics of high mortality and poor prognosis. Even in the presence of oxygen, the Warburg effect, a major metabolic hallmark of almost all cancer cells, is characterized by increased glycolysis and lactate fermentation, which supports biosynthesis and provides energy to sustain tumor cell growth and proliferation. However, a thorough investigation into glycolysis- and lactate-related genes and their association with COAD prognosis, immune cell infiltration, and drug candidates is currently lacking.Methods: COAD patient data and glycolysis- and lactate-related genes were retrieved from The Cancer Genome Atlas (TCGA) and Gene Set Enrichment Analysis (GSEA) databases, respectively. After univariate Cox regression analysis, a nonnegative matrix factorization (NMF) algorithm was used to identify glycolysis- and lactate-related molecular subtypes. Least absolute shrinkage and selection operator (LASSO) Cox regression identified twelve glycolysis- and lactate-related genes (ADTRP, ALDOB, APOBEC1, ASCL2, CEACAM7, CLCA1, CTXN1, FLNA, NAT2, OLFM4, PTPRU, and SNCG) related to prognosis. The median risk score was employed to separate patients into high- and low-risk groups. The prognostic efficacy of the glycolysis- and lactate-related gene signature was assessed using Kaplan–Meier (KM) survival and receiver operating characteristic (ROC) curve analyses. The nomogram, calibration curves, decision curve analysis (DCA), and clinical impact curve (CIC) were employed to improve the clinical applicability of the prognostic signature. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed on differentially expressed genes (DEGs) from the high- and low-risk groups. Using CIBERSORT, ESTIMATE, and single-sample GSEA (ssGSEA) algorithms, the quantities and types of tumor-infiltrating immune cells were assessed. The tumor mutational burden (TMB) and cytolytic (CYT) activity scores were calculated between the high- and low-risk groups. Potential small-molecule agents were identified using the Connectivity Map (cMap) database and validated by molecular docking. To verify key core gene expression levels, quantitative real-time polymerase chain reaction (qRT–PCR) assays were conducted.Results: We identified four distinct molecular subtypes of COAD. Cluster 2 had the best prognosis, and clusters 1 and 3 had poor prognoses. High-risk COAD patients exhibited considerably poorer overall survival (OS) than low-risk COAD patients. The nomogram precisely predicted patient OS, with acceptable discrimination and excellent calibration. GO and KEGG pathway enrichment analyses of DEGs revealed enrichment mainly in the “glycosaminoglycan binding,” “extracellular matrix,” “pancreatic secretion,” and “focal adhesion” pathways. Patients in the low-risk group exhibited a larger infiltration of memory CD4+ T cells and dendritic cells and a better prognosis than those in the high-risk group. The chemotherapeutic agent sensitivity of patients categorized by risk score varied significantly. We predicted six potential small-molecule agents binding to the core target of the glycolysis- and lactate-related gene signature. ALDOB and APOBEC1 mRNA expression was increased in COAD tissues, whereas CLCA1 and OLFM4 mRNA expression was increased in normal tissues.Conclusion: In summary, we identified molecular subtypes of COAD and developed a glycolysis- and lactate-related gene signature with significant prognostic value, which benefits COAD patients by informing more precise and effective treatment decisions.

Esophageal cancer (EC) is characterized by a high degree of malignancy and poor prognosis. N6-methyladenosine (m6A), a prominent post-transcriptional modification of mRNA in mammalian cells, plays a pivotal role in regulating various cellular and biological processes. Similarly, cuproptosis has garnered attention for its potential implications in cancer biology. This study seeks to elucidate the impact of m6A- and cuproptosis-related long non-coding RNAs (m6aCRLncs) on the prognosis of patients with EC. The EC transcriptional data and corresponding clinical information were retrieved from The Cancer Genome Atlas (TCGA) database, comprising 11 normal samples and 159 EC samples. Data on 23 m6A regulators and 25 cuproptosis-related genes were sourced from the latest literature. The m6aCRLncs linked to EC were identified through co-expression analysis. Differentially expressed m6aCRLncs associated with EC prognosis were screened using the limma package in R and univariate Cox regression analysis. Subtype clustering was performed to classify EC patients, enabling the investigation of differences in clinical outcomes and immune microenvironment across patient clusters. A risk prognostic model was constructed using least absolute shrinkage and selection operator (LASSO) regression. Its robustness was evaluated through survival analysis, risk stratification curves, and receiver operating characteristic (ROC) curves. Additionally, the model's applicability across various clinical features and molecular subtypes of EC patients was assessed. To further explore the model's utility in predicting the immune microenvironment, single-sample gene set enrichment analysis (ssGSEA), immune cell infiltration analysis, and immune checkpoint differential expression analysis were conducted. Drug sensitivity analysis was performed to identify potential therapeutic agents for EC. Finally, the mRNA expression levels of m6aCRLncs in EC cell lines were validated using reverse transcription quantitative polymerase chain reaction (RT-qPCR). We developed a prognostic risk model based on five m6aCRLncs, namely ELF3-AS1, HNF1A-AS1, LINC00942, LINC01389, and MIR181A2HG, to predict survival outcomes and characterize the immune microenvironment in EC patients. Analysis of molecular subtypes and clinical features revealed significant differences in cluster distribution, disease stage, and N stage between high- and low-risk groups. Immune profiling further identified distinct immune cell populations and functional pathways associated with risk scores, including positive correlations with naive B cells, resting CD4+ T cells, and plasma cells, and negative correlations with macrophages M0 and M1. Additionally, we identified key immune checkpoint-related genes with significant differential expression between risk groups, including TNFRSF14, TNFSF15, TNFRSF18, LGALS9, CD44, HHLA2, and CD40. Furthermore, nine candidate drugs with potential therapeutic efficacy in EC were identified: Bleomycin, Cisplatin, Cyclopamine, PLX4720, Erlotinib, Gefitinib, RO.3306, XMD8.85, and WH.4.023. Finally, RT-qPCR validation of the mRNA expression levels of m6aCRLncs in EC cell lines demonstrated that ELF3-AS1 expression was significantly upregulated in the EC cell lines KYSE-30 and KYSE-180 compared to normal esophageal epithelial cells. This study elucidates the role of m6aCRLncs in shaping the prognostic outcomes and immune microenvironment of EC. Furthermore, it identifies potential therapeutic agents with efficacy against EC. These findings hold significant promise for enhancing the survival of EC patients and provide valuable insights to inform clinical decision-making in the management of this disease.

Autophagy has an important association with tumorigenesis, progression, and prognosis. However, the mechanism of autophagy-regulated genes on the risk prognosis of bladder cancer (BC) patients has not been fully elucidated yet. In this study, we created a prognostic model of BC risk based on autophagy-related genes, which further illustrates the value of genes associated with autophagy in the treatment of BC. We first downloaded human autophagy-associated genes and BC datasets from Human Autophagy Database and The Cancer Genome Atlas (TCGA) database, and finally obtained differential prognosis-associated genes for autophagy by univariate regression analysis and differential analysis of cancer versus normal tissues. Subsequently, we downloaded two datasets from Gene Expression Omnibus (GEO), GSE31684 and GSE15307, to expand the total number of samples. Based on these genes, we distinguished the molecular subtypes (C1, C2) and gene classes (A, B) of BC by consistent clustering analysis. Using the genes merged from TCGA and the two GEO datasets, we conducted least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression analysis to obtain risk genes and construct autophagy-related risk prediction models. The accuracy of this risk prediction model was assessed by receiver operating characteristic (ROC) and calibration curves, and then nomograms were constructed to predict the survival of bladder cancer patients at 1, 3, and 5years, respectively. According to the median value of the risk score, we divided BC samples into the high- and low-risk groups. Kaplan-Meier (K-M) survival analysis was performed to compare survival differences between subgroups. Then, we used single sample gene set enrichment analysis (ssGSEA) for immune cell infiltration abundance, immune checkpoint genes, immunotherapy response, gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis, and tumor mutation burden (TMB) analysis for different subgroups. We also applied quantitative real-time polymerase chain reaction (PCR) and immunohistochemistry (IHC) techniques to verify the expression of these six genes in the model. Finally, we chose the IMvigor210 dataset for external validation. Six risk genes associated with autophagy (SPOCD1, FKBP10, NAT8B, LDLR, STMN3, and ANXA2) were finally screened by LASSO regression algorithm and multivariate Cox regression analysis. ROC and calibration curves showed that the model established was accurate and reliable. Univariate and multivariate regression analyses were used to verify that the risk model was an independent predictor. K-M survival analysis indicated that patients in the high-risk group had significantly worse overall survival than those in the low-risk group. Analysis by algorithms such as correlation analysis, gene set variation analysis (GSVA), and ssGSEA showed that differences in immune microenvironment, enrichment of multiple biologically active pathways, TMB, immune checkpoint genes, and human leukocyte antigens (HLAs) were observed in the different risk groups. Then, we constructed nomograms that predicted the 1-, 3-, and 5-year survival rates of different BC patients. In addition, we screened nine sensitive chemotherapeutic drugs using the correlation between the obtained expression status of risk genes and drug sensitivity results. Finally, the external dataset IMvigor210 verified that the model is reliable and efficient. We established an autophagy-related risk prognostic model that is accurate and reliable, which lays the foundation for future personalized treatment of bladder cancer.

Head and neck squamous cell carcinoma (HNSCC) is a malignancy of epithelial origin and with poor prognosis. Exploring the biomarkers and prognostic models that can contribute to early tumor detection is meaningful. A comprehensive analysis was conducted according to the stage-related signature genes of HNSCC, and a prognostic model was developed to validate their ability to predict the prognosis. The transcriptome profiles and clinical information of HNSCC patients were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) respectively. mRNA expressions of differentially expressed genes (DEGs) were analyzed in stage I-II patients and stage III-IV patients from TCGA by R packages. A protein-protein interaction (PPI) network and core-gene network map were constructed, and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to examine pathway enrichment. Kaplan-Meier, least absolute shrinkage and selection operator (LASSO), and multivariate Cox regression were applied to establish a stage-associated signature model. A Spearman analysis was conducted to examine the correlations between the characteristic genes and immune cell infiltration. Kaplan-Meier analysis and a receiver operating characteristic (ROC) curve were used to test the effectiveness of the model. Univariate multivariate Cox regression analyses were used to assess whether the risk score was an independent prognostic indicator for HNSCC. In TCGA cohort, 5 genes (i.e., BRINP1, IL17A, ALB, FOXA2, and ZCCHC12) in the constructed prognostic risk model were associated with prognosis. Patients in the low-risk group had a better prognosis outcome than those in the high-risk group. The predictive power was good because all the area under the curve (AUC) of the risk score was higher than 0.6. Risk score [hazard ratio (HR) =1.985; P<0.001] was an independent risk factor for the prognosis of HNSCC. The results in the GEO cohort were consistent with those in the TCGA cohort. We constructed and verified a prognostic risk model of stage-related signature genes for HNSCC based on the GEO and TCGA data. Due to the good predictive accuracy of this model, the prognosis of and the tumor immune cell infiltration with patients can be estimated.

Liver hepatocellular carcinoma (LIHC) continues to pose a major global health concern and is characterized by elevated mortality rates and a lack of effective therapies. This study aimed to explore differential gene expression linked to cellular senescence and pyroptosis in LIHC and to develop a prognostic risk model for use in clinical settings. We acquired datasets from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). DESeq2 was used to identify differentially expressed genes associated with cell senescence and pyrodeath. The least absolute shrinkage and selection operator (LASSO) regression model was developed using cellular senescence- and pyroptosis-related differentially expressed genes (CSR&PRDEGs), and its predictive performance was evaluated with Kaplan-Meier survival analysis and time-dependent receiver operating characteristic (ROC) curves. We also performed various functional analyses of the genes. These findings were validated by real-time fluorescence quantitative polymerase chain reaction (PCR). Using bioinformatics analysis, we developed a prognostic risk framework incorporating six critical genes: ANXA2, APOA1, EZH2, IGF2BP3, SQSTM1, and TNFRSF11B.The model demonstrated a statistically significant difference in overall survival between the high-risk and low-risk groups (p < 0.05). Additionally, real-time fluorescence quantitative PCR confirmed that genes ANXA2, APOA1, EZH2, IGF2BP3, SQSTM1, and TNFRSF11B were significantly overexpressed in the peripheral blood of patients with LIHC in comparison to normal volunteers, thereby validating the prognostic risk model's accuracy. This study systematically elucidated the functions of genes associated with senescence and pyroptosis in LIHC cells. The constructed prognostic risk model serves to guide the development of personalized treatment plans, enhance patient management via risk stratification, facilitate the identification of high-risk patients, intensify monitoring or implement proactive interventions, thereby providing a novel perspective for the diagnosis and treatment of LIHC.

Objective To establish a prognostic risk model for uterine corpus endometrial carcinoma (UCEC) based on alternative splicing (AS) event data from The Cancer Genome Atlas (TCGA) and assess the accuracy of the model. Methods TCGA and SpliceSeq databases were used to acquire a summary of AS events and clinical data related to UCEC. Bioinformatic analysis was performed to identify differentially expressed AS events in UCEC. Least absolute shrinkage and selection operator (LASSO) regression and multivariate Cox regression analyses were used for constructing a prognostic risk model. Next, using the receiver operating characteristic (ROC) curve, Kaplan-Meier survival analysis, and independent prognostic analysis, we assessed the accuracy of the model. In addition, a splicing network was established based on the association between potential splicing factors and AS events. Results We downloaded clinical data and AS events of 527 UCEC cases from TCGA and SpliceSeq databases, respectively. We obtained 18,779 survival-associated AS events in UCEC using univariate Cox regression analysis and 487 AS events using LASSO regression analysis. Multivariate Cox regression analysis established a prognostic risk model for UCEC based on the percentage splicing value of 13 AS events. Independent prognostic effect on UCEC risk was then assessed using multivariate and univariate Cox regression analyses (P &lt; 0.001). The area under the curve was 0.827. The pathological stage and risk score were independent prognostic factors for UCEC. Herein, we established a regulatory network between alternative endometrial cancer-related splicing events and splicing factors. Conclusion We constructed a prognostic model of UCEC based on 13 AS events by analyzing datasets from TCGA and SpliceSeq databases with medium accuracy. The pathological stage and risk score were independent prognostic factors in the prognostic risk model.

Hepatocellular carcinoma (HCC) is a prevalent type of cancer with high incidence and mortality rates. It is the third most common cause of cancer-related deaths. CD8+ T cell exhaustion (TEX) is a progressive decline in T cell function due to sustained T cell receptor stimulation from continuous antigen exposure. Studies have shown that CD8+ TEX plays an important role in the anti-tumor immune process and is significantly correlated with patient prognosis. The aim of the research is to establish a reliable CD8+ TEX-based signature using single-cell RNA sequencing (scRNA-seq) and high-throughput RNA sequencing (RNA-seq), providing a new approach to evaluate HCC patient prognosis and immune microenvironment. The RNA-seq data of HCC patients were download from three different databases: The Cancer Genome Atlas (TCGA), the Gene Expression Omnibus (GEO), and the International Cancer Genome Consortium (ICGC). HCC's 10× scRNA data were acquired from GSE149614. Based on single-cell sequencing data, CD8+ TEX-related genes were identified using uniform manifold approximation and projection (UMAP) algorithm, singleR, and marker gene methods. Afterwards, we proceeded to construct CD8+ TEX signature using differential gene analysis, univariate Cox regression analysis, least absolute shrinkage and selection operator (LASSO) regression, and multivariate Cox regression analysis. We also validated the CD8+ TEX signature in GEO and ICGC external cohorts and investigated clinical characteristics, chemotherapy sensitivity, mutation landscape, functional analysis, and immune cell infiltration in different risk groups. The CD8+ TEX signature, consisting of 13 genes (HSPD1, UBB, DNAJB4, CALM1, LGALS3, BATF, COMMD3, IL7R, FDPS, DRAP1, RPS27L, PAPOLA, GPR171), was found to have a strong predictive effect on the prognosis of HCC. The Kaplan-Meier (KM) analysis showed that the overall survival (OS) rate of patients in the low-risk group was higher than that of patients in the high-risk group across different datasets and specific populations. The research findings suggested that the risk score was an independent predictor of HCC prognosis. The model based on clinical features and risk score has a strong predictive effect. We observed significant differences among various risk groups in terms of clinical characteristics, functional analysis, mutation landscape, chemotherapy sensitivity, and immune cell infiltration. We constructed a CD8+ TEX signature to predict the survival probability of patients with HCC. We also found that the model could predict the sensitivity of targeted drugs and immune cell infiltration, and the risk score was negatively correlated with CD8+ T cell infiltration. In summary, the CD8+ TEX signature of HCC was constructed for the prediction of prognosis and immune microenvironment by integrated analysis of bulk and scRNA-seq data.

Breast cancer (BRCA) is a malignant cancer that threatened the life of female with unsatisfactory prognosis. The aim of this study was to identify prognostic nuclear receptors (NRs) signature of BRCA. BRCA patient samples were collected from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database. Consensus clustering analysis, univariate Cox regression analysis and the least absolute shrinkage and selection operator (LASSO) Cox regression analysis were performed to evaluate, select NRs as prognostic factors and build Risk Score model. GSEAanalysis was explored to check signaling differences between High- and Low-Risk group. Nomogram model basing on age and Risk Score was established to predict the 1-, 3- and 5-year survival. Model performance was assessed by a time-dependent receiver operating characteristic (ROC) curve and calibration plot. CIBERSORT, ESTIMATE and TIMER algorithm were introduced to evaluate the immune landscape. NR3C1, NR4A3, THRA, RXRG, NR2F6, NR1D2 and RORB were optimized as a prognostic signature for BRCA. This seven-NR-based Risk Score could effectively predict overall survival status. The area under the curve (AUC) of 1-, 3- and 5-year overall survival are 0.702, 0.734 and 0.722 in TCGA training cohort, and 0.630, 0.721 and 0.823 in GEO validation cohort, respectively. Calibration plot demonstrated satisfactory agreement between predictive and observed outcomes. Nomogram model worked well on predicting survival probabilities. Multiple cancer-related pathways were highly enriched in High-Risk group. High- and Low-Risk groups showed significant differed immune cell infiltration. There exists an obvious connection between Risk Score and immune checkpoints LAG3, PD1 and TIM3. The seven-NR-based Risk Score represents a promising signature for estimating overall survival in patients with BRCA, and is correlated with the immune microenvironment.

Immune therapy has become first-line treatment option for patients with lung cancer, but some patients respond poorly to immune therapy, especially among patients with lung adenocarcinoma (LUAD). Novel tools are needed to screen potential responders to immune therapy in LUAD patients, to better predict the prognosis and guide clinical decision-making. Although many efforts have been made to predict the responsiveness of LUAD patients, the results were limited. During the era of immunotherapy, this study attempts to construct a novel prognostic model for LUAD by utilizing differentially expressed genes (DEGs) among patients with differential immune therapy responses. Transcriptome data of 598 patients with LUAD were downloaded from The Cancer Genome Atlas (TCGA) database, which included 539 tumor samples and 59 normal control samples, with a mean follow-up time of 29.69 months (63.1% of patients remained alive by the end of follow-up). Other data sources including three datasets from the Gene Expression Omnibus (GEO) database were analyzed, and the DEGs between immunotherapy responders and nonresponders were identified and screened. Univariate Cox regression analysis was applied with the TCGA cohort as the training set and GSE72094 cohort as the validation set, and least absolute shrinkage and selection operator (LASSO) Cox regression were applied in the prognostic-related genes which fulfilled the filter criteria to establish a prognostic formula, which was then tested with time-dependent receiver operating characteristic (ROC) analysis. Enriched pathways of the prognostic-related genes were analyzed with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses, and tumor immune microenvironment (TIME), tumor mutational burden, and drug sensitivity tests were completed with appropriate packages in R (The R Foundation of Statistical Computing). Finally, a nomogram incorporating the prognostic formula was established. A total of 1,636 DEGs were identified, 1,163 prognostic-related DEGs were extracted, and 34 DEGs were selected and incorporated into the immunotherapy responsiveness-related risk score (IRRS) formula. The IRRS formula had good performance in predicting the overall prognoses in patients with LUAD and had excellent performance in prognosis prediction in all LUAD subgroups. Moreover, the IRRS formula could predict anticancer drug sensitivity and immunotherapy responsiveness in patients with LUAD. Mechanistically, immune microenvironments varied profoundly between the two IRRS groups; the most significantly varied pathway between the high-IRRS and low-IRRS groups was ribonucleoprotein complex biogenesis, which correlated closely with the TP53 and TTN mutation burdens. In addition, we established a nomogram incorporating the IRRS, age, sex, clinical stage, T-stage, N-stage, and M-stage as predictors that could predict the prognoses of 1-year, 3-year, and 5-year survival in patients with LUAD, with an area under curve (AUC) of 0.718, 0.702, and 0.68, respectively. The model we established in the present study could predict the prognosis of LUAD patients, help to identify patients with good responses to anticancer drugs and immunotherapy, and serve as a valuable tool to guide clinical decision-making.

Pancreatic adenocarcinoma (PAAD) is one of the deadliest malignancies. Aging is described as the degeneration of physiological function, which is complexly correlated with cancer. It is significant to explore the influences of aging-related genes (ARGs) on PAAD. Based on The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) datasets, we used univariate Cox regression analysis and acquired eight differentially expressed ARGs with prognostic values. Two molecular subtypes were identified based on these ARGs to depict PAAD patients’ overall survival (OS) and immune microenvironments preliminarily. Cluster 1 had a poor OS as well as a worse immune microenvironment. Through least absolute shrinkage and selection operator (LASSO) regression analysis, we constructed a seven-ARG risk signature based on the TCGA dataset and verified it in Gene Expression Omnibus (GEO) and International Cancer Genome Consortium (ICGC) to predict the prognoses, immune microenvironments, signal pathways, tumor mutations, and drug sensitivity of PAAD patients. The high-risk group possessed an unfavorable OS compared with that of the low-risk group. We also verified the independence and clinical availability of the risk signature by Cox regression analyses and the establishment of a nomogram, respectively. The higher risk score was associated with several clinical factors such as higher grade and advanced tumor stage as well as lower immunoscore and cluster 1. The negative associations of risk scores with immune, stroma, and estimate scores proved the terrible immune microenvironment in the high-risk group. Relationships between risk score and immune checkpoint gene expression as well as signal pathways provided several therapeutic targets. PAAD patients in the low-risk group possessed lower tumor mutations as well as a higher susceptibility to axitinib and vorinostat. The high-risk group bore a higher TMB and cisplatin and dasatinib may be better options. We used immunohistochemistry and qPCR to confirm the expression of key ARGs with their influences on OS. In conclusion, we identified two ARG-mediated molecular subtypes and a novel seven-ARG risk signature to predict prognoses, immune microenvironments, signal pathways, tumor mutations, and drug sensitivity of PAAD patients.

Background: Colorectal cancer (CRC) is one gastrointestinal malignancy, accounting for 10% of cancer diagnoses and cancer-related deaths worldwide each year. Therefore, it is urgent to identify genes involved in CRC predicting the prognosis. Methods: CRC’s data were acquired from the Gene Expression Omnibus (GEO) database (GSE39582 and GSE41258 datasets) and The Cancer Genome Atlas (TCGA) database. The differentially expressed necroptosis-related genes (DENRGs) were sorted out between tumor and normal tissues. Univariate Cox regression analysis and least absolute shrinkage and selectionator operator (LASSO) analysis were applied to selected DENRGs concerning patients’ overall survival and to construct a prognostic biomarker. The effectiveness of this biomarker was assessed by the Kaplan–Meier curve and the receiver operating characteristic (ROC) analysis. The GSE39582 dataset was utilized as external validation for the prognostic signature. Moreover, using univariate and multivariate Cox regression analyses, independent prognostic factors were identified to construct a prognostic nomogram. Next, signaling pathways regulated by the signature were explored through the gene set enrichment analysis (GSEA). The single sample gene set enrichment analysis (ssGSEA) algorithm and tumor immune dysfunction and exclusion (TIDE) were used to explore immune correlation in the two groups, high-risk and low-risk ones. Finally, prognostic genes’ expression was examined in the GSE41258 dataset. Results: In total, 27 DENRGs were filtered, and a necroptosis-related prognostic signature based on 6 DENRGs was constructed, which may better understand the overall survival (OS) of CRC. The Kaplan–Meier curve manifested the effectiveness of the prognostic signature, and the ROC curve showed the same result. In addition, univariate and multivariate Cox regression analyses revealed that age, pathology T, and risk score were independent prognostic factors, and a nomogram was established. Furthermore, the prognostic signature was most significantly associated with the apoptosis pathway. Meanwhile, 24 immune cells represented significant differences between two groups, like the activated B cell. Furthermore, 32 immune checkpoints, TIDE scores, PD-L1 scores, and T-cell exclusion scores were significantly different between the two groups. Finally, a 6-gene prognostic signature represented different expression levels between tumor and normal samples significantly in the GSE41258 dataset. Conclusion: Our study established a signature including 6 genes and a prognostic nomogram that could significantly assess the prognosis of patients with CRC.

BackgroundPancreatic cancer (PC) is a highly fatal and aggressive disease with its incidence and mortality quite discouraging. An effective prediction model is urgently needed for the accurate assessment of patients’ prognosis to assist clinical decision-making.MethodsGene expression data and clinicopathological data of the samples were acquired from The Cancer Genome Atlas (TCGA), Genotype-Tissue Expression (GTEx), and Gene Expression Omnibus (GEO) databases. Differential expressed genes (DEGs) analysis, univariate Cox regression analysis, least absolute shrinkage and selection operator (LASSO) regression analysis, random forest screening and multivariate Cox regression analysis were applied to construct the risk signature. The effectiveness and independence of the model were validated by time-dependent receiver operating characteristic (ROC) curve, Kaplan–Meier (KM) survival analysis and survival point graph in training set, test set, TCGA entire set and GSE57495 set. The validity of the core gene was verified by immunohistochemistry and our own independent cohort. Meanwhile, functional enrichment analysis of DEGs between the high and low risk groups revealed the potential biological pathways. Finally, CMap database and drug sensitivity assay were utilized to identify potential small molecular drugs as the risk model-related treatments for PC patients.ResultsFour histone modification-related genes were identified to establish the risk signature, including CBX8, CENPT, DPY30 and PADI1. The predictive performance of risk signature was validated in training set, test set, TCGA entire set and GSE57495 set, with the areas under ROC curve (AUCs) for 3-year survival were 0.773, 0.729, 0.775 and 0.770 respectively. Furthermore, KM survival analysis, univariate and multivariate Cox regression analysis proved it as an independent prognostic factor. Mechanically, functional enrichment analysis showed that the poor prognosis of high-risk population was related to the metabolic disorders caused by inadequate insulin secretion, which was fueled by neuroendocrine aberration. Lastly, a cluster of small molecule drugs were identified with significant potentiality in treating PC patients.ConclusionsBased on a histone modification-related gene signature, our model can serve as a reliable prognosis assessment tool and help to optimize the treatment for PC patients. Meanwhile, a cluster of small molecule drugs were also identified with significant potentiality in treating PC patients.

BackgroundHepatocellular carcinoma (HCC) has one of the highest mortality rates worldwide. Abnormal glutamine metabolism (GM) has been reported to be involved in HCC progression. The current study sought to examine the predictive value of GM in HCC patient’s prognosis and therapy response.MethodsThe RNA-sequencing data and clinical information of HCC samples were obtained from The Cancer Genome Atlas (TCGA) database (N=377) and Gene Expression Omnibus (GEO) database (N=242). By analyzing a data set from TCGA, we showed that the GM landscape of HCC patients was developed based on the non-negative matrix factorization (NMF) algorithm. Univariate Cox regression and least absolute shrinkage and selection operator (LASSO)–penalized Cox regression analyses were used to construct a risk model. The accuracy of the model, which was based on the GM-related genes (GMRGs), was verified by Kaplan-Meier (K-M) and receiver operating characteristic (ROC) curves. We also verified the reliability of the model based on GEO data. Finally, the immune infiltration analysis, pathway enrichment analysis, and treatment response prediction results were compared to each other in the 2 risk groups.ResultsIn our study, the HCC samples were divided into 2 GM-related patterns; that is, C1 and C2. The multi-analysis revealed that the GM-related patterns were associated with the pathologic stage, T stages, N stages, histologic grade, and the tumor immune microenvironment (TIME). Next, the prognostic model containing 5 GMRGs (i.e., aldehyde dehydrogenase 5 family member A1, ASNSD1, carbamoyl-phosphate synthetase 1, GMPS, and PPAT) was constructed to calculate the risk score. The high-risk group of HCC patients had significantly worse overall survival (OS) than the low-risk group in both datasets (P<0.001). Multivariate Cox regression uncover the riskScores may serve as an independent prognostic marker for HCC patients [TCGA: hazard ratio (HR) =2.909 (1.940−4.362), P<0.001; GEO: HR =2.911 (1.753−5.848), P=0.043]. Finally, we found that the prognostic model was significantly correlated with the pathologic stage and TIME of the HCC patients in both databases. Moreover, the prognostic model may guide the immunotherapy, chemotherapy, and targeted drugs choice.ConclusionsIn summary, we developed a GM-related 5-gene risk-score model, which may be a useful tool for predicting prognosis and guiding the treatment of HCC patients.

Unfolded protein response signature unveils novel insights into breast cancer prognosis and tumor microenvironment